The analogy between the wave equation of liquid and the Navier equations of structural elasticity is examined in detail. By introducing appropriate parameters, the structural counterpart of the liquid sloshing model c...The analogy between the wave equation of liquid and the Navier equations of structural elasticity is examined in detail. By introducing appropriate parameters, the structural counterpart of the liquid sloshing model can be easily built. Therefore, the dynamic analysis of liquid sloshing can be reduced to that of structural elasticity, and the existing FEM structural analysis computer programs can be applied to liquid sloshing analysis without any modification. The present method also reveals the internal relationship between liquid sloshing and structural vibration. The effectiveness and reliability of the method is illustrated by the numerical example.展开更多
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 ...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.展开更多
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 con...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.展开更多
On the basis of the characteristics of broken rock zone, using the program of "discontinuous deformation analysis(DDA)", the changing law of influential factors of discontinuous rock mass in large broken roc...On the basis of the characteristics of broken rock zone, using the program of "discontinuous deformation analysis(DDA)", the changing law of influential factors of discontinuous rock mass in large broken rock zone was researched quantitatively for the first time. Based on the results of computation, the concept of "key part"of roadways and its stability criterion were brought forward, and it was pointed out that in inclined coal and rock seams the"key parts"of roadways are the upper side and the floor of surrounding rocks, especially the former.展开更多
Triaxial creep tests on CCG specimens were systematically performed using aself-made creep seepage experimental apparatus for determining the creep law of CCG.An improved triaxial creep model of CCG was established on...Triaxial creep tests on CCG specimens were systematically performed using aself-made creep seepage experimental apparatus for determining the creep law of CCG.An improved triaxial creep model of CCG was established on the basis of a Nishiharamodel and another visco-elasto-plastic model,parameters of which were fitted on test data.Furthermore,the creep model is validated according to the result of triaxial creep experiments,and the outcome shows that the proposed triaxial creep model can properly characterizethe properties of various creep deformation phases of CCG,especially the acceleratingcreep phase.At the same time,the instability conditions of CCG were presentedbased on the discussion of the improved model's stability in terms of stability theories ofdifferential equation solution.展开更多
We study a new model named the Green-Lindsay type therm-elastic model for nonhomogeneous media that consists of a system of dynamic thermoelasticity equations of displacement and dynamic heat conduction equation. We c...We study a new model named the Green-Lindsay type therm-elastic model for nonhomogeneous media that consists of a system of dynamic thermoelasticity equations of displacement and dynamic heat conduction equation. We construct the model based on the classical GL-model for homogeneous material. This system is coupled dynamic problem and the displacement field and heat field must be solved at the same time. By using Fadeo- Galerkin method, we proved that the problem we proposed exist unique weak solution under some regular assumption.展开更多
Deep rock mass possesses some unusual properties due to high earth stress,which further result in new problems that have not been well understood and explained up to date.In order to investigate the deformation mechan...Deep rock mass possesses some unusual properties due to high earth stress,which further result in new problems that have not been well understood and explained up to date.In order to investigate the deformation mechanism,the complete deformation process of deep rock mass,with a great emphasis on local shear deformation stage,was analyzed in detail.The quasi continuous shear deformation of the deep rock mass is described by a combination of smooth functions:the averaged distribution of the original deformation field,and the local discontinuities along the slip lines.Hence,an elasto-plastic model is established for the shear deformation process,in which the rotational displacement is taken into account as well as the translational component.Numerical analysis method was developed for case study.Deformation process of a tunnel under high earth stress was investigated for verification.展开更多
Without considering the influence of heat,existing fractal contact models are not applicable to analyze the contacts when the temperature changes.For this problem,the normal load model and the normal stiffness model o...Without considering the influence of heat,existing fractal contact models are not applicable to analyze the contacts when the temperature changes.For this problem,the normal load model and the normal stiffness model of thermal elasto-plastic contact of rough surfaces are developed respectively in this paper.The proposed model is based on the normal contact mechanics model of fractal theory of anisotropic and thermal elasto-plastic contact theory which can be used to characterize the rough surface thermodynamic properties.Then the validity of the model is verified.Finally,the influence of main parameters on the total normal load and the whole normal stiffness of thermal elasto-plastic contact at the interface is analyzed by contact simulation.The results show that the total normal load of thermal elasto-plastic contact increases with the increases of temperature.The whole normal stiffness of thermal elasto-plastic contact increases with increasing coefficient of linear expansion,scale factor,temperature difference or fractal dimension,but decreases with increasing fractal roughness.This model expands basic theory and applications of traditional models,and can be used to calculate and analyze the contacts when the temperature changes.展开更多
Numerous experimental studies reveal that the mechanical and deformational behaviors of sands are dependent on the combined effect of void ratio and stress. To predict this complex behavior of sands, a hypo-elastic mo...Numerous experimental studies reveal that the mechanical and deformational behaviors of sands are dependent on the combined effect of void ratio and stress. To predict this complex behavior of sands, a hypo-elastic model is developed based on the cross-anisotropic elasticity model, which involves four parameters: bulk module, tangent Young's module, volume deformation coefficient and Poisson ratio. A parameter defined as virtual peak deviatoric stress dependent on state parameter is introduced into hyperbolic stress strain relationship to determine tangent Young's module. In addition, an existing fitting equation for isotropic compression curves and an existing dilatancy equation, which can consider the effect of state of sands, are employed to determine bulk module and volume deformation coefficient. Thirteen model constants are involved in the proposed model, the values of which are fixed for a sand over a wide range of initial void ratios and initial confining pressures. Well known experimental data for drained and undrained triaxial compression tests of Toyoura sand are successfully modeled.展开更多
Existing seismic prediction methods struggle to effectively discriminate between fluids in tight gas reservoirs,such as those in the Sulige gas field in the Ordos Basin,where porosity and permeability are extremely lo...Existing seismic prediction methods struggle to effectively discriminate between fluids in tight gas reservoirs,such as those in the Sulige gas field in the Ordos Basin,where porosity and permeability are extremely low and the relationship between gas and water is complicated.In this paper,we have proposed a comprehensive seismic fluid identification method that combines ray-path elastic impedance(REI)inversion with fluid substitution for tight reservoirs.This approach is grounded in geophysical theory,forward modeling,and real data applications.We used geophysics experiments in tight gas reservoirs to determine that Brie's model is better suited to calculate the elastic parameters of mixed fluids than the conventional Wood’s model.This yielded a more reasonable and accurate fluid substitution model for tight gas reservoirs.We developed a forward model and carried out inversion of REI.which reduced the non-uniqueness problem that has plagued elastic impedance inversion in the angle domain.Our well logging forward model in the ray-path domain with different fluid saturations based on a fluid substitution model proved that REI identifies fluids more accurately when the ray parameters are large.The distribution of gas saturation can be distinguished from the crossplot of REI(p=0.10)and porosity.The inverted ray-path elastic impedance profile was further used to predict the porosity and gas saturation profile.Our new method achieved good results in the application of 2D seismic data in the western Sulige gas field.展开更多
Constitutive models aimed at predicting the mechanical response of lead-core bearing devices for passive seismic isolation are discussed in this paper. The study is focused on single-degree-of-freedom models which pro...Constitutive models aimed at predicting the mechanical response of lead-core bearing devices for passive seismic isolation are discussed in this paper. The study is focused on single-degree-of-freedom models which provide a relation between the shear displacement (shear strain) and the shear force (shear stress) in elastomeric and lad-core rubber bearings. Classical Bouc-Wen model along with a numerical procedure for identification of the model constants is described. Alternatively, a constitutive relation introducing a damage variable aimed at assessing the material degradation is also considered.展开更多
Analytical solutions are presented using method of separation of variables for the time periodic electroosmotic flow (EOF) of linear viscoelastic fluids in semicircular microchannel. The linear viscoelastic fluids use...Analytical solutions are presented using method of separation of variables for the time periodic electroosmotic flow (EOF) of linear viscoelastic fluids in semicircular microchannel. The linear viscoelastic fluids used here are described by the general Maxwell model. The solution involves analytically solving the linearized Poisson-Boltzmann (P -B) equation, together with the Cauchy momentum equation and the general Maxwell constitutive equation. By numerical computations, the influences of electric oscillating Reynolds number Re and Deborah number De on velocity amplitude are presented. For small Re, results show that the larger velocity amplitude is confined to the region near the charged wall when De is small. With the increase of the Deborah number De, the velocity far away the charged wall becomes larger for large Deborah number De. However, for larger Re, the oscillating characteristic of the velocity amplitude occurs and becomes significant with the increase of De, especially for larger Deborah number.展开更多
The Cellular Automaton(CA) modeling and simulation of solid dynamics is a long-standing difficult problem.In this paper we present a new two-dimensional CA model for solid dynamics.In this model the solid body is repr...The Cellular Automaton(CA) modeling and simulation of solid dynamics is a long-standing difficult problem.In this paper we present a new two-dimensional CA model for solid dynamics.In this model the solid body is represented by a set of white and black particles alternatively positioned in the x-and y-directions.The force acting on each particle is represented by the linear summation of relative displacements of the nearest-neighboring particles.The key technique in this new model is the construction of eight coefficient matrices.Theoretical and numerical analyses show that the present model can be mathematically described by a conservative system.So,it works for elastic material.In the continuum limit the CA model recovers the well-known Navier equation.The coefficient matrices are related to the shear module and Poisson ratio of the material body.Compared with previous CA model for solid body,this model realizes the natural coupling of deformations in the x-and y-directions.Consequently,the wave phenomena related to the Poisson ratio effects are successfully recovered.This work advances significantly the CA modeling and simulation in the field of computational solid dynamics.展开更多
文摘The analogy between the wave equation of liquid and the Navier equations of structural elasticity is examined in detail. By introducing appropriate parameters, the structural counterpart of the liquid sloshing model can be easily built. Therefore, the dynamic analysis of liquid sloshing can be reduced to that of structural elasticity, and the existing FEM structural analysis computer programs can be applied to liquid sloshing analysis without any modification. The present method also reveals the internal relationship between liquid sloshing and structural vibration. The effectiveness and reliability of the method is illustrated by the numerical example.
基金supported by the National Natural Science Foundation of China (No. 11002025, 40114066)the National Basic Research Program of China (973 Program) (No.2007CB209505)the RIPED Youth Innovation Foundation (No. 2010-A-26-01)
文摘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.
基金sponsored by the National Natural Science Foundation of China(under Grant Nos.41404090 and U1262208
文摘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.
文摘On the basis of the characteristics of broken rock zone, using the program of "discontinuous deformation analysis(DDA)", the changing law of influential factors of discontinuous rock mass in large broken rock zone was researched quantitatively for the first time. Based on the results of computation, the concept of "key part"of roadways and its stability criterion were brought forward, and it was pointed out that in inclined coal and rock seams the"key parts"of roadways are the upper side and the floor of surrounding rocks, especially the former.
基金Supported by the National Natural Science Foundation of China(50874124)the National Basic Research Program of China(973)(2005CB221502)+1 种基金the National Natural Science Foundation of China(50534080)the Natural Science Foundation Project of CQ CSTC(2008BA6028)
文摘Triaxial creep tests on CCG specimens were systematically performed using aself-made creep seepage experimental apparatus for determining the creep law of CCG.An improved triaxial creep model of CCG was established on the basis of a Nishiharamodel and another visco-elasto-plastic model,parameters of which were fitted on test data.Furthermore,the creep model is validated according to the result of triaxial creep experiments,and the outcome shows that the proposed triaxial creep model can properly characterizethe properties of various creep deformation phases of CCG,especially the acceleratingcreep phase.At the same time,the instability conditions of CCG were presentedbased on the discussion of the improved model's stability in terms of stability theories ofdifferential equation solution.
基金Foundation item: Supported by the National Natural Science Foundation of China(10771198)
文摘We study a new model named the Green-Lindsay type therm-elastic model for nonhomogeneous media that consists of a system of dynamic thermoelasticity equations of displacement and dynamic heat conduction equation. We construct the model based on the classical GL-model for homogeneous material. This system is coupled dynamic problem and the displacement field and heat field must be solved at the same time. By using Fadeo- Galerkin method, we proved that the problem we proposed exist unique weak solution under some regular assumption.
基金Project(50825403) supported by the National Science Fund for Distinguished Young ScholarsProject(2010CB732003) supported by the National Key Basic Research Program of ChinaProject(51021001) supported by the Science Fund for Creative Research Group of the National Natural Science Foundation of China
文摘Deep rock mass possesses some unusual properties due to high earth stress,which further result in new problems that have not been well understood and explained up to date.In order to investigate the deformation mechanism,the complete deformation process of deep rock mass,with a great emphasis on local shear deformation stage,was analyzed in detail.The quasi continuous shear deformation of the deep rock mass is described by a combination of smooth functions:the averaged distribution of the original deformation field,and the local discontinuities along the slip lines.Hence,an elasto-plastic model is established for the shear deformation process,in which the rotational displacement is taken into account as well as the translational component.Numerical analysis method was developed for case study.Deformation process of a tunnel under high earth stress was investigated for verification.
基金Project(52130501)supported by the National Natural Science Foundation of ChinaProject(LY20E050012)supported by the Natural Science Foundation of Zhejiang Province,ChinaProject(Y201942581)supported by the Scientific Research Project of Education Department of Zhejiang Province,China。
文摘Without considering the influence of heat,existing fractal contact models are not applicable to analyze the contacts when the temperature changes.For this problem,the normal load model and the normal stiffness model of thermal elasto-plastic contact of rough surfaces are developed respectively in this paper.The proposed model is based on the normal contact mechanics model of fractal theory of anisotropic and thermal elasto-plastic contact theory which can be used to characterize the rough surface thermodynamic properties.Then the validity of the model is verified.Finally,the influence of main parameters on the total normal load and the whole normal stiffness of thermal elasto-plastic contact at the interface is analyzed by contact simulation.The results show that the total normal load of thermal elasto-plastic contact increases with the increases of temperature.The whole normal stiffness of thermal elasto-plastic contact increases with increasing coefficient of linear expansion,scale factor,temperature difference or fractal dimension,but decreases with increasing fractal roughness.This model expands basic theory and applications of traditional models,and can be used to calculate and analyze the contacts when the temperature changes.
基金Project(2010BC732101)supported by the National Basic Research Program of China
文摘Numerous experimental studies reveal that the mechanical and deformational behaviors of sands are dependent on the combined effect of void ratio and stress. To predict this complex behavior of sands, a hypo-elastic model is developed based on the cross-anisotropic elasticity model, which involves four parameters: bulk module, tangent Young's module, volume deformation coefficient and Poisson ratio. A parameter defined as virtual peak deviatoric stress dependent on state parameter is introduced into hyperbolic stress strain relationship to determine tangent Young's module. In addition, an existing fitting equation for isotropic compression curves and an existing dilatancy equation, which can consider the effect of state of sands, are employed to determine bulk module and volume deformation coefficient. Thirteen model constants are involved in the proposed model, the values of which are fixed for a sand over a wide range of initial void ratios and initial confining pressures. Well known experimental data for drained and undrained triaxial compression tests of Toyoura sand are successfully modeled.
基金supported by the National Science and Technology Major Project(No.2016ZX05050 and 2017ZX05069)CNPC Major Technology Special Project(No.2016E-0503)
文摘Existing seismic prediction methods struggle to effectively discriminate between fluids in tight gas reservoirs,such as those in the Sulige gas field in the Ordos Basin,where porosity and permeability are extremely low and the relationship between gas and water is complicated.In this paper,we have proposed a comprehensive seismic fluid identification method that combines ray-path elastic impedance(REI)inversion with fluid substitution for tight reservoirs.This approach is grounded in geophysical theory,forward modeling,and real data applications.We used geophysics experiments in tight gas reservoirs to determine that Brie's model is better suited to calculate the elastic parameters of mixed fluids than the conventional Wood’s model.This yielded a more reasonable and accurate fluid substitution model for tight gas reservoirs.We developed a forward model and carried out inversion of REI.which reduced the non-uniqueness problem that has plagued elastic impedance inversion in the angle domain.Our well logging forward model in the ray-path domain with different fluid saturations based on a fluid substitution model proved that REI identifies fluids more accurately when the ray parameters are large.The distribution of gas saturation can be distinguished from the crossplot of REI(p=0.10)and porosity.The inverted ray-path elastic impedance profile was further used to predict the porosity and gas saturation profile.Our new method achieved good results in the application of 2D seismic data in the western Sulige gas field.
文摘Constitutive models aimed at predicting the mechanical response of lead-core bearing devices for passive seismic isolation are discussed in this paper. The study is focused on single-degree-of-freedom models which provide a relation between the shear displacement (shear strain) and the shear force (shear stress) in elastomeric and lad-core rubber bearings. Classical Bouc-Wen model along with a numerical procedure for identification of the model constants is described. Alternatively, a constitutive relation introducing a damage variable aimed at assessing the material degradation is also considered.
基金Supported by the National Natural Science Foundation of China under Grant Nos.11062005, 11202092Opening Fund of State Key Laboratory of Nonlinear Mechanics, the Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region+2 种基金the Natural Science Foundation of Inner Mongolia under Grant Nos.2010BS0107, 2012MS0107the Research Start up Fund for Excellent Talents at Inner Mongolia University under Grant No.Z20080211the Natural Science Key Fund of Inner Mongolia under Grant No.2009ZD01
文摘Analytical solutions are presented using method of separation of variables for the time periodic electroosmotic flow (EOF) of linear viscoelastic fluids in semicircular microchannel. The linear viscoelastic fluids used here are described by the general Maxwell model. The solution involves analytically solving the linearized Poisson-Boltzmann (P -B) equation, together with the Cauchy momentum equation and the general Maxwell constitutive equation. By numerical computations, the influences of electric oscillating Reynolds number Re and Deborah number De on velocity amplitude are presented. For small Re, results show that the larger velocity amplitude is confined to the region near the charged wall when De is small. With the increase of the Deborah number De, the velocity far away the charged wall becomes larger for large Deborah number De. However, for larger Re, the oscillating characteristic of the velocity amplitude occurs and becomes significant with the increase of De, especially for larger Deborah number.
基金Supported by the Science Foundations of China Academy of Engineering Physics under Grant Nos. 2012B0101014 and 2011A0201002National Natural Science Foundation of China under Grant Nos. 11075021,91130020,and 11202003Foundation of State Key Laboratory of Explosion Science and Technology
文摘The Cellular Automaton(CA) modeling and simulation of solid dynamics is a long-standing difficult problem.In this paper we present a new two-dimensional CA model for solid dynamics.In this model the solid body is represented by a set of white and black particles alternatively positioned in the x-and y-directions.The force acting on each particle is represented by the linear summation of relative displacements of the nearest-neighboring particles.The key technique in this new model is the construction of eight coefficient matrices.Theoretical and numerical analyses show that the present model can be mathematically described by a conservative system.So,it works for elastic material.In the continuum limit the CA model recovers the well-known Navier equation.The coefficient matrices are related to the shear module and Poisson ratio of the material body.Compared with previous CA model for solid body,this model realizes the natural coupling of deformations in the x-and y-directions.Consequently,the wave phenomena related to the Poisson ratio effects are successfully recovered.This work advances significantly the CA modeling and simulation in the field of computational solid dynamics.