A viscoelastic metamaterial beam for integrated vibration isolation and energy harvesting

Locally resonant metamaterials have low-frequency band gaps and the capability of converging vibratory energy in the band gaps at resonant cells.It has been demonstrated by several researchers that the dissipatioin of vibratory energy within the band gap can be improved by using viscoelastic materials.This paper designs an integrated viscoelastic metamaterial for energy harvesting and vibration isolation.The viscoelastic metamaterial is achieved by a viscoelastic beam periodically arrayed with spatial ball-pendulum nonlinear energy harvesters.The nonlinear resonator with an energy harvesting function is achieved by placing a free-rolling magnetic ball in a spherical cavity with an additional induction coil.The dynamic equations of viscoelastic metamaterials under transverse excitation are established,and the energy harvesting and vibration isolation characteristics within the dispersion relation of viscoelastic metamaterials are analyzed.The results show that the vibrations of the main body of the viscoelastic metamaterial beam are significantly suppressed in the frequency range of the local resonance band gap.At the same time,the elastic waves are limited in the nonlinear resonator with an energy harvesting function,which improves the energy output.Finally,an experimental platform of viscoelastic metamaterial vibration is established for validation purposes.

A robust viscoelastic surfactant tolerating 20% HCl up to 150℃ for oil well stimulation

Viscoelastic surfactants(VES)are often used as viscous diverters in acidizing stimulation to prolong the acid consumption time and maximize zonal coverage of the acid for improving well productivity.However,the ceiling temperature of commercial VES cannot exceed 120℃in practical use because of the poor thermal stability and fragile molecular structure,hindering their implementation in hightemperature oil reservoirs,i.e.,≥150℃.Here we synthesized a novel C22-tailed diamine,N-erucaminopropyl-N,N-dimethylamine(EDPA),and examined comparatively its rheological behavior,assemblies morphology and molecular stability in 20 wt%HCl with a commercial VES,erucyl dimethyl amidopropyl betaine(EDAB).The feasibility of EDPA for acidizing stimulation was assessed by acid etching of carbonate rock with its HCl solution at 150℃.Rheological results showed that the 2.5 wt%EDPA—20 wt%HCl solution maintains stable viscosity of 90 m Pa s at 150℃for 60 min,while that of 2.0 wt%EDAB HCl solution is just 1 m Pa s under identical conditions.1H NMR spectra and cryo-TEM observations revealed that the chemical structure and self-assembled architectures of EDPA remained intact in such context,but the EDAB suffered from degradation due to the hydrolysis of the amide group,accounting for the poor heat-resistance and acid-tolerance.The reaction rate of 2.5 wt%EDPA HCl solution with carbonate rock was one order of magnitude lower than that of 20 wt%HCl solution at 150℃,underpinning the potential of EDPA to be used in the high-temperature reservoirs acidizing.This work improved the thermal tolerance of VES in highly concentrated HCl solution,paving a feasible way for the acidization of high-temperature reservoir environments(~150℃).

The effect of reactive plasticizer on viscoelastic and mechanical properties of solid rocket propellants based on different types of HTPB resin

Conventional plasticizers deteriorate mechanical and viscoelastic properties of the propellants due to their migration upon aging and long-term storage,which affects reliability and safety properties during exploitation.To address this issue,conventional plasticizer,dioctyl adipate(DOA),is replaced by reactive one,castor oil(CO).In addition,three different types of HTPB were used to obtain propellants with designed viscoelastic and mechanical properties.The CO increased propellants viscosity,without a significant impact on the propellant processability,regardless to the type of prepolymer.Conversely,mechanical properties were different depending on the type of resin,which were further analyzed by gel permeation chromatography(GPC).Addition of CO formed a denser polymer network and shifted T_(g) to higher values,compared to the compositions with DOA.The tensile strength of CO-containing propellants was lower at +20℃ and +50℃ compared to the reference compositions,while the strain at maximum load and strain at break were significantly increased with pronounced plastic deformation,especially for samples at -30℃.The inclusion of CO in the propellants composition gives more room for adjusting a wide range of mechanical properties.

Development of a Three-Dimensional Multiscale Octree SBFEM for Viscoelastic Problems of Heterogeneous Materials

The multiscale method provides an effective approach for the numerical analysis of heterogeneous viscoelastic materials by reducing the degree of freedoms(DOFs).A basic framework of the Multiscale Scaled Boundary Finite Element Method(MsSBFEM)was presented in our previous works,but those works only addressed two-dimensional problems.In order to solve more realistic problems,a three-dimensional MsSBFEM is further developed in this article.In the proposed method,the octree SBFEM is used to deal with the three-dimensional calculation for numerical base functions to bridge small and large scales,the three-dimensional image-based analysis can be conveniently conducted in small-scale and coarse nodes can be flexibly adjusted to improve the computational accuracy.Besides,the Temporally Piecewise Adaptive Algorithm(TPAA)is used to maintain the computational accuracy of multiscale analysis by adaptive calculation in time domain.The results of numerical examples show that the proposed method can significantly reduce the DOFs for three-dimensional viscoelastic analysis with good accuracy.For instance,the DOFs can be reduced by 9021 times compared with Direct Numerical Simulation(DNS)with an average error of 1.87%in the third example,and it is very effective in dealing with three-dimensional complex microstructures directly based on images without any geometric modelling process.

Intraocular pressure before and after capsulorhexis using two viscoelastic substances and two surgical approaches in enucleated porcine eyes

AIM:To investigate the influence of ophthalmic viscoelastic devices(OVDs)and different surgical approaches on the intraocular pressure(IOP)before and after creation of the curvilinear circular capsulorhexis(CCC)as a measure for anterior chamber stability during this maneuver.METHODS:Prospective experimental WetLab study carried out on enucleated porcine eyes.IOP was measured before and after CCC with the iCare Rebound tonometer(iCare ic200;iCare Finland Oy,Vantaa,Finland).The OVDs used were a cohesive one[Z-Hyalin,Carl Zeiss Meditec AG,Germany;hyaluronic acid(HA)]and a dispersive[Z-Celcoat,Carl Zeiss Meditec AG,Germany;hydroxy propylmethylcellulosis(HPMC)].The CCC was created using Utrata forceps or 23 g microforceps in different combinations with the OVDs.RESULTS:Using the Utrata forceps the IOP dropped from 63.65±6.44 to 11.25±3.63 mm Hg during the CCC.The use of different OVDs made no difference.Using the 23 g microforceps the IOP dropped from 65.35±8.15 to 36.55±6.09 mm Hg.The difference between IOP drop using either Utrata forceps or 23 g microforceps was highly significant regardless of the OVD used.CONCLUSION:Using the sideport for the creation of the capsulorhexis leads to a lesser drop in IOP during this maneuver compared to the main incision in enucleated porcine eyes.The use of different OVD has no significant influence on IOP drop.

A flexible multiscale algorithm based on an improved smoothed particle hydrodynamics method for complex viscoelastic flows

Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the nonlinear dynamic behaviors of viscoelastic fluids.However,traditional grid-based multiscale methods are confined to simple viscoelastic flows with short relaxation time,and there is a lack of uniform multiscale scheme available for coupling different solvers in the simulations of viscoelastic fluids.In this paper,a universal multiscale method coupling an improved smoothed particle hydrodynamics(SPH)and multiscale universal interface(MUI)library is presented for viscoelastic flows.The proposed multiscale method builds on an improved SPH method and leverages the MUI library to facilitate the exchange of information among different solvers in the overlapping domain.We test the capability and flexibility of the presented multiscale method to deal with complex viscoelastic flows by solving different multiscale problems of viscoelastic flows.In the first example,the simulation of a viscoelastic Poiseuille flow is carried out by two coupled improved SPH methods with different spatial resolutions.The effects of exchanging different physical quantities on the numerical results in both the upper and lower domains are also investigated as well as the absolute errors in the overlapping domain.In the second example,the complex Wannier flow with different Weissenberg numbers is further simulated by two improved SPH methods and coupling the improved SPH method and the dissipative particle dynamics(DPD)method.The numerical results show that the physical quantities for viscoelastic flows obtained by the presented multiscale method are in consistence with those obtained by a single solver in the overlapping domain.Moreover,transferring different physical quantities has an important effect on the numerical results.

Spatially-graded 3D-printed viscoelastic truss metamaterials for impact trajectory control and energy absorption

This study demonstrates that two-and three-dimensional spatially graded,truss-based polymeric-material metamaterials can be designed for beneficial impact mitigation and energy absorption capabilities.Through a combination of numerical and experimental techniques,we highlight the broad property space of periodic viscoelastic trusses,realized using 3D printing via selective laser sintering.Extending beyond periodic designs,we investigate the impact response of spatially variant viscoelastic lattices in both two and three dimensions.Our result reveal that introducing spatial variations in lattice topology allows for redirecting of the impact trajectory,opening new opportunities for engineering and tailoring lightweight materials with target impact functionality.This is achieved through the combined selection of base material and metamaterial design.

Tensile Shock Physics in Compressible Thermoviscoelastic Solid Medium

This paper addresses tensile shock physics in thermoviscoelastic (TVE) solids without memory. The mathematical model is derived using conservation and balance laws (CBL) of classical continuum mechanics (CCM), incorporating the contravariant second Piola-Kirchhoff stress tensor, the covariant Green’s strain tensor, and its rates up to order n. This mathematical model permits the study of finite deformation and finite strain compressible deformation physics with an ordered rate dissipation mechanism. Constitutive theories are derived using conjugate pairs in entropy inequality and the representation theorem. The resulting mathematical model is both thermodynamically and mathematically consistent and has closure. The solution of the initial value problems (IVPs) describing evolutions is obtained using a variationally consistent space-time coupled finite element method, derived using space-time residual functional in which the local approximations are in hpk higher-order scalar product spaces. This permits accurate description problem physics over the discretization and also permits precise a posteriori computation of the space-time residual functional, an accurate measure of the accuracy of the computed solution. Model problem studies are presented to demonstrate tensile shock formation, propagation, reflection, and interaction. A unique feature of this research is that tensile shocks can only exist in solid matter, as their existence requires a medium to be elastic (presence of strain), which is only possible in a solid medium. In tensile shock physics, a decrease in the density of the medium caused by tensile waves leads to shock formation ahead of the wave. In contrast, in compressive shocks, an increase in density and the corresponding compressive waves result in the formation of compression shocks behind of the wave. Although these are two similar phenomena, they are inherently different in nature. To our knowledge, this work has not been reported in the published literature.

Hermite Finite Element Method for Vibration Problem of Euler-Bernoulli Beam on Viscoelastic Pasternak Foundation

Viscoelastic foundation plays a very important role in civil engineering. It can effectively disperse the structural load into the foundation soil and avoid the damage caused by the concentrated load. The model of Euler-Bernoulli beam on viscoelastic Pasternak foundation can be used to analyze the deformation and response of buildings under complex geological conditions. In this paper, we use Hermite finite element method to get the numerical approximation scheme for the vibration equation of viscoelastic Pasternak foundation beam. Convergence and error estimation are rigourously established. We prove that the fully discrete scheme has convergence order O(τ2+h4), where τis time step size and his space step size. Finally, we give four numerical examples to verify the validity of theoretical analysis.

BISQ model based on a Kelvin-Voigt viscoelastic frame in a partially saturated porous medium

Taking into account three important porous media mechanisms during wave propagation （the Biot-flow, squirt-flow, and solid-skeleton viscoelastic mechanisms）, we introduce water saturation into the dynamic governing equations of wave propagation by analyzing the effective medium theory and then providing a viscoelastic Biot/squirt （BISQ） model which can analyze the wave propagation problems in a partially viscous pore fluid saturated porous media. In this model, the effects of pore fluid distribution patterns on the effective bulk modulus at different frequencies are considered. Then we derive the wave dynamic equations in the time-space domain. The phase velocity and the attenuation coefficient equations of the viscoelatic BISQ model in the frequency-wavenumber domain are deduced through a set of plane harmonic solution assumptions. Finally, by means of numerical simulations, we investigate the effects of water saturation, permeability, and frequency on compressional wave velocity and attenuation. Based on tight sandstone and carbonate experimental observed data, the compressional wave velocities of partially saturated reservoir rocks are calculated. The compressional wave velocity in carbonate reservoirs is more sensitive to gas saturation than in sandstone reservoirs.

Modeling and analysis of azimuthal AVO responses from a viscoelastic anisotropic reflector

We propose a method for mOdeling azimuthal AVO responses from a fractured i reflector. The method calculates the integrated reflected wavetrains, and the wavetrains contain elastodynamic information including the contrast in impedance and anelasticity i across interfaces, the intemal anisotropic propagation, the dispersion and attenuation along i the wave path, and tuning and interference. The results suggest that for large angles of incidence, the velocity dispersion and attenuation increase the amplitudes of PP waves from the top and decrease those from the bottom. For azimuthal responses at specific angles of incidence, the reflected wavetrains of PP waves tend to have longer duration with increasing azimuth. In contrast, model-converted PSV and PSH reflections show stable azimuthal features and are less affected by the reflector thickness. The amplitudes of PSV reflections increase with increasing azimuth; moreover, the waves have no reflection energy at 0° and 90° azimuth and maximum amplitude at 45° azimuth.

A generalized viscoelastic Biot/squirt model for clay-bearing sandstones in a wide range of permeabilities

To simultaneously take into account the Biot-flow mechanism, the squirt-flow mechanism, and the frame-viscoelasticity mechanism, a generalized viscoelastic BISQ （Biot/squirt） model is developed for wave propagation in clay-bearing sandstones based on Dvorkin＇s elastic BISQ model. The present model is extended to a wide range of permeabilities （k 〉 0.05 mD） by introducing a dimensionless correction factor for viscoelastic parameters, defined as a function of the permeability and the clay content. We describe the frame＇s stress-strain relationship of the clay-bearing sandstones by the differential constitutive equations of generalized viscoelasticity and then derive the viscoelastic-wave dynamic equations. With the assumption of a plane-wave solution, we finally yield the phase velocities and the attenuation coefficients by solving the dynamic wave equations in the frequency and wave number domain. The comparison of numerical results and experimental data shows that the generalized viscoelastic BISQ model is applicable for modeling the wave propagation in most of the sandstones mainly bearing kaolinite clay.

Effect of viscoelasticity on nonlinear vibration of single microbubble

Based on the Church-Hoff model, the nonlinear oscillations of a single encapsulated microbubble with a finite thickness shell are theoretically studied. The effects of viscoelasticity on radial oscillations and the fundamental and harmonic components are researched. The peaks of radial oscillations and magnitudes of power spectra of the fundamental and harmonic components all increase gradually with the shear modulus of shell varying from 0 to 10 MPa by an interval of 0. 1 MPa at the same shear viscosity, while they decrease as the shear viscosity increases from 0 to 1 Pa · s by an interval of 0. 01 Pa · s at the same shear modulus. The fluctuation ranges of subharmonic and ultraharmonic signals are much larger than both the fundamental and second harmonic components. It means that the effect of viscoelasticity on the subharmonic and ultraharmonic signals is greater than that on the fundamental and second harmonic components. So adjusting the viscoelasticity of the shell is a potential method to obtain a perfect microbubble contrast agent used for the subharmonic and ultraharmonic imaging. Four points with significant fundamental and harmonic components are chosen as an example： a shear viscosity of 0. 39 Pa · s with shear modulus of 3.9, 6. 6, and 8.6 MPa, respectively; a shear modulus of 6.6 MPa with a shear viscosity of 0.42 Pa · s.

Reflective cracking viscoelastic response of asphalt concrete under dynamic vehicle loading

In order to investigate the mechanical response of reflective cracking in asphalt concrete pavement under dynamic vehicle loading, a finite element model is established in ABAQUS. The viscoelastic behavior is described by a prony series which is calculated through nonlinear fitting to the creep test data obtained in the laboratory. Based on the viscoelastic theory, the time-temperature equivalence principle, fracture mechanics and the dynamic finite element method, both the Jintegral and the mix-mode stress intensity factor are utilized as fracture evaluation parameters, and a half-sine dynamic loading is used to simulate the vehicle loading. Finally, the mechanical response of the pavement reflective cracking is analyzed under different vehicle speeds, different environmental conditions and various damping factors. The results indicate that increasing either the vehicle speed or the structure damping factor decreases the maximum values of fracture parameters, while the structure temperature has little effect on the fracture parameters. Due to the fact that the vehicle speed can be enhanced by improving the road traffic conditions, and the pavement damping factor can become greater by modifying the components of materials, the development of reflective cracking can be delayed and the asphalt pavement service life can be effectively extended through both of these ways.

Elastic and Viscoelastic Poisson’s Ratios: The Theoretical Mechanics Perspective

A recent review publication presented an extensive and comprehensive assessment of the phenomenological relations of Poisson’s ratios (PRs) to the behavior and responses of contemporary materials under specific loading conditions. The present review and analysis paper is intended as a theoretical mechanics complement covering mathematical and physical modeling of a single original elastic and of six time and process (i.e. path and stress) dependent viscoelastic PR definitions as well as a seventh special path independent one. The implications and consequences of such models on material characterization are analyzed and summarized. Indeed, PRs based on experimentally obtained 2-D strains under distinct creep and/or relaxation processes exhibit radically different time responses for identical material specimen. These results confirm the PR’s implicit path dependence in addition to their separate intrinsic time reliance. Such non-uniqueness of viscoelastic PRs renders them unsuitable as universal material descriptors. Analytical formulations and experimental measurements also examine the physical impossibility of instantaneously achieving time independent loads or strains or their rates thus making certain PR definitions based on constant state variables, while mathematically valid, physically unrealistic and unachievable. A newly developed theoretical/experimental protocol for the determination of the time when loading patterns reach stead-state conditions based on strain accelerations demonstrates the capability to measure this time from experimental data. Due to the process dependent PRs, i.e. stress and stress history paths, the non-existence of a unique viscoelastic PR and of a universal elastic-viscoelastic correspondence principle or analogy (EVCP) in terms of PRs is demonstrated. Additionally and independently, the required double convolution integral construction of linear viscoelastic constitutive relations with the inclusion of PRs is cumbersome analytically and computationally needlessly highly CPU intensive. Furthermore, there is no theoretical fundamental hint as to what loading path is required to produce a unique universal viscoelastic PR definition necessary for formulating a PR based constitutive relation or an EVCP protocol. The analysis associated with an additional Class VII viscoelastic PR establishes it as a universal representation which is loading path and strain independent while still remaining time dependent. This Class PR can be the one used if it is desired to express constitutive relations in terms of PRs, subject to the caveat applying to all PR Classes regarding the CPU intensiveness in the time space due to triple product and double convolution integral constitutive relations. However, the use PRs is unnecessary since any set of material behavior can be uniquely and completely defined in terms of only moduli and/or compliances. The mathematical model of instantaneous initial loading paths, based on Heavi-side functions, is examined in detail and shown to lead to infinite velocities and accelerations. Additionally, even if non-instantaneous gradual loading functions are employed the resulting PRs are still load and load history dependent. Consequently, they represent specialized PR responses applicable and limited to those particular load and history combinations. Although the analyses contained herein are generalized to non-homogeneous linear viscoelastic materials, the main focus is on PR time and process dependence. The non-homogeneous material results and conclusions presented herein apply equally to homogeneous viscoelasticity and per se do not influence the results or conclusions of the analytical development regarding viscoelastic PRs. In short, these PR analyses apply to all linear viscoelastic material characterization.

On wave dispersion of rotating viscoelastic nanobeam based on general nonlocal elasticity in thermal environment

The present research focuses on the analysis of wave propagation on a rotating viscoelastic nanobeam supported on the viscoelastic foundation which is subject to thermal gradient effects.A comprehensive and accurate model of a viscoelastic nanobeam is constructed by using a novel nonclassical mechanical model.Based on the general nonlocal theory(GNT),Kelvin-Voigt model,and Timoshenko beam theory,the motion equations for the nanobeam are obtained.Through the GNT,material hardening and softening behaviors are simultaneously taken into account during wave propagation.An analytical solution is utilized to generate the results for torsional(TO),longitudinal(LA),and transverse(TA)types of wave dispersion.Moreover,the effects of nonlocal parameters,Kelvin-Voigt damping,foundation damping,Winkler-Pasternak coefficients,rotating speed,and thermal gradient are illustrated and discussed in detail.

Non-harmonic resonance of viscoelastic structures subjected to time-dependent exponentially decreasing transverse distributed loads

In this paper,the non-harmonic resonance of Bernoulli viscoelastic beams,Kirchhoff viscoelastic plates,Timoshenko viscoelastic beams,and Mindlin viscoelastic plates subjected to time-dependent exponentially decreasing transverse distributed load is investigated for the first time.The constitutive equations are expressed utilizing Boltzmann integral law with a constant bulk modulus.The displacement vector is approximated by employing the separation of variables method.The Laplace transformation is used to transfer equations from the time domain to the Laplace domain and vice versa.The novel point of the proposed method is to express,prove and calculate the critical time in which the displacement will be several times the displacement at time zero.In addition,this new method calculates the maximum deflection at the critical time,explicitly and exactly,without any need to follow the time-displacement curve with a low computational cost.Additionally,the proposed method introduces the critical range of time so that the responses are greater than the responses at time zero.

An approach to viscoelastic c haracterization of dispersive media by inversion of a general wave propagation model

In the characterization of elastic properties of tissue using dynarmic optical coherence elasto-graphy,shear/surface waves are propagated and tracked in order to estimate speed and Y oung's modulus.However,for dispersive tssues,the displacement pulse is highly damped and distorted during propagation,diminishing the ffectiveness of peak tracking approaches,and leading to biased cstimates of wave speed.Further,plane wave propagation is sometimes assumed,which contributes to estimation erors.Therefore,we invert a wave propagation model that incorpo-rates propagation,decay,and distortion of pulses in a dispersive media in order to accurately estimate its elastic and viscous components.The model uses a general first-order approximation of dispersion,avoiding the use of any particular rheological model of tisue.Experiments are conducted in elastic and viscoelastic tissue mimicking phantoms by producing a Gaussian push using acoustic radiation force excitation and measuring the wave propagation using a Fourier domain optical coherence tomography system.Results confirmed the ffectiveness of the inversion method in est imat ing viscoelastic parameters in both the viscoelastic and elastic phantoms when compared to mechanical measurements.Finally,the viscoelastic characterization of a fresh porcine comea was conducted.Preliminary results validate this approach when compared to other methods.

Computational-Analysis of the Non-Isothermal Dynamics of the Gravity-Driven Flow of Viscoelastic-Fluid-Based Nanofluids Down an Inclined Plane

The paper explores the gravity-driven flow of the thin film of a viscoelastic-fluid-based nanofluids(VFBN)along an inclined plane under non-isothermal conditions and subjected to convective cooling at the free-surface.The Newton’s law of cooling is used to model the convective heat-exchange with the ambient at the free-surface.The Giesekus viscoelastic constitutive model,with appropriate modifications to account for non-isothermal effects,is employed to describe the polymeric effects.The unsteady and coupled non-linear partial differential equations(PDEs)describing the model problem are obtained and solved via efficient semi-implicit numerical schemes based on finite difference methods(FDM)implemented in Matlab.The response of the VFBN velocity,temperature,thermal-conductivity and polymeric-stresses to variations in the volume-fraction of embedded nanoparticles is investigated.It is shown that these quantities all increase as the nanoparticle volume-fraction becomes higher.

On Estimating Magnitude of a Maximum Sequent Earthquake by Viscoelastic Coulomb Stress Change and a Discussion of the Relationship between the M_S7.3 Earthquakes in Yutian 2008 and 2014

On the basis of the previous studies of the layered crustal model in the Yutian area,combined with the field GPS continuous observation data,we roughly estimate the viscous coefficient of each layer. With the viscoelastic horizontal layer model,we calculate the viscoelastic co-seismic Coulomb stress change caused by the Yutian M_S7. 3 earthquakes 2008 and 2014 respectively. Based on the Coulomb stress change,using the calculation method of "direct "aftershock frequency,we come up with the theoretical earthquake frequency directly related to the mainshock and the co-seismic Coulomb stress change in the study area. Then we put forward a method,based on the comparison of theoretical and actual earthquake frequency or the comparison between theoretical and practical earthquake frequency-distance decay curve fitting residuals,to estimate the magnitude of a maximum sequent earthquake,directly related to the mainshock co-seismic Coulomb stress change. Results calculated by different methods show that the maximum follow-up earthquake magnitude caused by the coseismic Coulomb stress change lies from M_S7. 2 to M_S7. 5 following Yutian M_S7. 3 earthquake in 2008; but that of the 2014 Yutian M_S7. 3 earthquake is M_S6. 3. The former is very close to the Yutian M_S7. 3 earthquake in 2014.Because of the same magnitude,relatively close spatial distance,short time interval,the same region of the external force,the strong correlation between two seismic tectonic and a clear stress interaction,we thus consider that the two Yutian M_S7. 3 earthquakes in 2008 and 2014 constitute a pair of generalized double shock type earthquake. This is consistent with the sequence type characteristic of past "double shock"earthquakes in the region. In this paper,the influence of the magnitude lower limit and the b-value in the relationship of G-R on the results is discussed. As a result,when the viscoelastic coseismic Coulomb stress variation is determined,the lower limit of magnitude has little effect on the maximum sequent earthquake magnitude estimation,but b-value of G-R has a greater impact on the results.