Time-domain dynamic constitutive model suitable for mucky soil site seismic response

Soil nonlinear behavior displays noticeable effects on the site seismic response.This study proposes a new functional expression of the skeleton curve to replace the hyperbolic skeleton curve.By integrating shear modulus and combining the dynamic skeleton curve and the damping degradation coefficient,the constitutive equation of the logarithmic dynamic skeleton can be obtained,which considers the damping effect in a soil dynamics problem.Based on the finite difference method and the multi-transmitting boundary condition,a 1D site seismic response analysis program called Soilresp1D has been developed herein and used to analyze the time-domain seismic response in three types of sites.At the same time,this study also provides numerical simulation results based on the hyperbolic constitutive model and the equivalent linear method.The results verify the rationality of the new soil dynamic constitutive model.It can analyze the mucky soil site nonlinear seismic response,reflecting the deformation characteristics and damping effect of the silty soil.The hysteresis loop area is more extensive,and the residual strain is evident.

A statistical damage-based constitutive model for shearing of rock joints in brittle drop mode

Some rock joints exhibit significant brittleness,characterized by a sharp decrease in shear stress upon reaching the peak strength.However,existing models often fail to accurately represent this behavior and are encumbered by numerous parameters lacking clear mechanical significance.This study presents a new statistical damage constitutive model rooted in both damage mechanics and statistics,containing only three model parameters.The proposed model encompasses all stages of joint shearing,including the compaction stage,linear stage,plastic yielding stage,drop stage,strain softening stage,and residual strength stage.To derive the analytical expression of the constitutive model,three boundary conditions are introduced.Experimental data from both natural and artificial rock joints is utilized to validate the model,resulting in average absolute relative errors ranging from 3%to 8%.Moreover,a comparative analysis with established models illustrates that the proposed model captures stress drop and post-peak strain softening more effectively,with model parameters possessing clearer mechanical interpretations.Furthermore,parameter analysis is conducted to investigate the impacts of model parameters on the curves and unveil the relationship between these parameters and the mechanical properties of rock joints.Importantly,the proposed model is straightforward in form,and all model parameters can be obtained from direct shear tests,thus facilitating the utilization in numerical simulations.

Fatigue properties and damage constitutive model of salt rock based on CT scanning

To investigate the macroscopic fatigue properties and the mesoscopic pore evolution characteristics of salt rock under cyclic loading,fatigue tests under different upper-limit stresses were carried out on salt rock,and the mesoscopic pore structures of salt rock before and after fatigue tests and under different cycle numbers were measured using CT scanning instrument.Based on the test results,the effects of the cycle number and the upper-limit stress on the evolution of cracks,pore morphology,pore number,pore volume,pore size,plane porosity,and volume porosity of salt rock were analyzed.The failure path of salt rock specimens under cyclic loading was analyzed using the distribution law of plane porosity.The damage variable of salt rock under cyclic loading was defined on basis of the variation of volume porosity with cycle number.In order to describe the fatigue deformation behavior of salt rock under cyclic loading,the nonlinear Burgers damage constitutive model was further established.The results show that the model established can better reflect the whole development process of fatigue deformation of salt rock under cyclic loading.

Parameter calibration of the tensile-shear interactive damage constitutive model for sandstone failure

The tensile-shear interactive damage(TSID)model is a novel and powerful constitutive model for rock-like materials.This study proposes a methodology to calibrate the TSID model parameters to simulate sandstone.The basic parameters of sandstone are determined through a series of static and dynamic tests,including uniaxial compression,Brazilian disc,triaxial compression under varying confining pressures,hydrostatic compression,and dynamic compression and tensile tests with a split Hopkinson pressure bar.Based on the sandstone test results from this study and previous research,a step-by-step procedure for parameter calibration is outlined,which accounts for the categories of the strength surface,equation of state(EOS),strain rate effect,and damage.The calibrated parameters are verified through numerical tests that correspond to the experimental loading conditions.Consistency between numerical results and experimental data indicates the precision and reliability of the calibrated parameters.The methodology presented in this study is scientifically sound,straightforward,and essential for improving the TSID model.Furthermore,it has the potential to contribute to other rock constitutive models,particularly new user-defined models.

Damage constitutive model of lunar soil simulant geopolymer under impact loading

Lunar base construction is a crucial component of the lunar exploration program,and considering the dynamic characteristics of lunar soil is important for moon construction.Therefore,investigating the dynamic properties of lunar soil by establishing a constitutive relationship is critical for providing a theoretical basis for its damage evolution.In this paper,a split Hopkinson pressure bar(SHPB)device was used to perform three sets of impact tests under different pressures on a lunar soil simulant geopolymer(LSSG)with sodium silicate(Na_(2)SiO_(3))contents of 1%,3%,5%and 7%.The dynamic stressestrain curves,failure modes,and energy variation rules of LSSG under different pressures were obtained.The equation was modified based on the ZWT viscoelastic constitutive model and was combined with the damage variable.The damage element obeys the Weibull distribution and the constitutive equation that can describe the mechanical properties of LSSG under dynamic loading was obtained.The results demonstrate that the dynamic compressive strength of LSSG has a marked strain-rate strengthening effect.Na_(2)SiO_(3) has both strengthening and deterioration effects on the dynamic compressive strength of LSSG.As Na_(2)SiO_(3) grows,the dynamic compressive strength of LSSG first increases and then decreases.At a fixed air pressure,5%Na_(2)SiO_(3) had the largest dynamic compressive strength,the largest incident energy,the smallest absorbed energy,and the lightest damage.The ZWT equation was modified according to the stress response properties of LSSG and the range of the SHPB strain rate to obtain the constitutive equation of the LSSG,and the model’s correctness was confirmed.

A rate-dependent constitutive model for saturated frozen soil considering local breakage mechanism

A rate-dependent constitutive model for saturated frozen soil is vital in frozen soil mechanics,especially when simultaneously describing the nonlinearity,dilatancy and strain-softening characteristics.The distribution of the non-uniform strain rate of saturated frozen soil at the meso-scale due to the local icecementation breakage is described by a newly binary-medium-based homogenization equation.Based on the field-equation-based approach of the meso-mechanics theory,the interaction expression of the strain rate at macro-and meso-scale is derived,which can give the strain rate concentration tensor at different crushed degrees.With the thermodynamics and empirical assumption,a breakage ratio in the rate-dependent form is determined.This overcomes the limitations of the existing binary-medium-based models that are difficult to simulate rate-dependent mechanical response.Based on these assumptions,a newly binary-medium-based rate-dependent model is proposed considering both the ice bond breakage and material composition characteristics of saturated frozen soil.The proposed constitutive model has been validated by the test results on frozen soils with different temperatures and strain rates.

Damage constitutive model for strain-softening rock based on normal distribution and its parameter determination

Firstly, using the damage model for rock based on Lemaitre hypothesis about strain equivalence, a new technique for measuring strength of rock micro-cells by adopting the Mohr-Coulomb criterion was developed, and a statistical damage evolution equation was established based on the property that strength of micro-cells is consistent with normal distribution function, through discussing the characteristics of random distributions for strength of micro-cells, then a statistical damage constitutive model that can simulate the full process of rock strain softening under specific confining pressure was set up. Secondly, a new method to determine the model parameters which can be applied to the situations under different confining pressures was proposed, by deeply studying the relations between the model parameters and characteristic parameters of the full stress-strain curve under different confining pressures. Therefore, a unified statistical damage constitutive model for rock softening which can reflect the effect of different confining pressures was set up. This model makes the physical property of model parameters explicit, contains only conventional mechanical parameters, and leads its application more convenient. Finally, the rationality of this model and its parameters-determining method were identified via comparative analyses between theoretical and experimental curves.

Creep characteristics and constitutive model for frozen mixed soils

The uniform settle caused by the permafrost creep is still the primary problem of engineering construction in cold region.To investigate the creep characteristics of frozen soils mixed with silty clay and coarse-grained sand,several triaxial creep tests of frozen mixed soils under different conditions(temperature,confining pressure,coarse-grained particle content) were performed,and the effects of the temperature,confining pressure and coarse particle content on the creep characteristics of frozen mixed soils were also analyzed.The results of the experiments indicated that when the confining pressure was low,the specimen exhibited an attenuation creep under a low-stress level(0.4-0.7) and a non-attenuation creep under a high-stress level(0.7-0.9).In contrast,when the confining pressure was high,the specimen had both initial and stable creep stages,but no the accelerated creep stage.The higher the content of coarse grains,the shorter the stable creep stage,and the easier to enter the accelerated creep stage for the specimen.Further,the long-term strength of frozen soils decreased with an increase in the content of coarse grains.Finally,a newly improved Nishihara model was proposed to consider both the hardening effect and damage effect by introducing both the hardening and damage variables,which can model the entire creep process of frozen soils modeled relatively easily.It was found that with the increasing content of coarse particle,both the strengthening and damaged effects in the creep process are reduced。

An elasto-plastic and viscoplastic damage constitutive model for dilatancy and fracturing behavior of soft rock squeezing deformation

Soft rock squeezing deformation mainly consists of pre-peak damage-dilatancy and post-peak fracture-bulking at the excavation unloading instant,and creep-dilatancy caused by time-dependent damage and fracturing.Based on the classic elastoplastic and Perzyna over-stress viscoplastic theories,as well as triaxial unloading confining pressure test and triaxial unloading creep test results,an elastoplastic and viscoplastic damage constitutive model is established for the short-and long-term dilatancy and fracturing behavior of soft rock squeezing deformation.Firstly,the criteria for each deformation and failure stage are expressed as a linear function of confining pressure.Secondly,the total damage evolution equation considering time-dependent damage is proposed,including the initial damage produced at the excavation instant,in which the damage variable increases exponentially with the lateral strain,and creep damage.Thirdly,a transient five-stages elasto-plastic constitutive equation for the short-term deformation after excavation that comprised of elasticity,pre-peak damage-dilatancy,post-peak brittle-drop,linear strain-softening,and residual perfectly-plastic regimes is developed based on incremental elasto-plastic theory and the nonassociated flow rule.Fourthly,regarding the timedependent properties of soft rock,based on the Perzyna viscoplastic over-stress theory,a viscoplastic damage model is set up to capture creep damage and dilatancy behavior.Viscoplastic strain is produced when the stress exceeds the initial static yield surface fs;the distance between the static yield surface fs and the dynamic yield surface fd determines the viscoplastic strain rate.Finally,the established constitutive model is numerically implemented and field applied to the-848 m belt conveyer haulage roadway of Huainan Panyidong Coal Mine.Laboratory test results and in-situ monitoring results validate the rationality of the established constitutive model.The presented model takes both the transient and time-dependent damage and fracturing into consideration.

Identification of structure and parameters of rheological constitutive model for rocks using differential evolution algorithm

To determine structure and parameters of a rheological constitutive model for rocks,a new method based on differential evolution(DE) algorithm combined with FLAC3D(a numerical code for geotechnical engineering) was proposed for identification of the global optimum coupled of model structure and its parameters.At first,stochastic coupled mode was initialized,the difference in displacement between the numerical value and in-situ measurements was regarded as fitness value to evaluate quality of the coupled mode.Then the coupled-mode was updated continually using DE rule until the optimal parameters were found.Thus,coupled-mode was identified adaptively during back analysis process.The results of applications to Jinping tunnels in China show that the method is feasible and efficient for identifying the coupled-mode of constitutive structure and its parameters.The method overcomes the limitation of the traditional method and improves significantly precision and speed of displacement back analysis process.

Anisotropic time-dependent behaviors of shale under direct shearing and associated empirical creep models

Understanding the anisotropic creep behaviors of shale under direct shearing is a challenging issue.In this context,we conducted shear-creep and steady-creep tests on shale with five bedding orientations (i.e.0°,30°,45°,60°,and 90°),under multiple levels of direct shearing for the first time.The results show that the anisotropic creep of shale exhibits a significant stress-dependent behavior.Under a low shear stress,the creep compliance of shale increases linearly with the logarithm of time at all bedding orientations,and the increase depends on the bedding orientation and creep time.Under high shear stress conditions,the creep compliance of shale is minimal when the bedding orientation is 0°,and the steady-creep rate of shale increases significantly with increasing bedding orientations of 30°,45°,60°,and 90°.The stress-strain values corresponding to the inception of the accelerated creep stage show an increasing and then decreasing trend with the bedding orientation.A semilogarithmic model that could reflect the stress dependence of the steady-creep rate while considering the hardening and damage process is proposed.The model minimizes the deviation of the calculated steady-state creep rate from the observed value and reveals the behavior of the bedding orientation's influence on the steady-creep rate.The applicability of the five classical empirical creep models is quantitatively evaluated.It shows that the logarithmic model can well explain the experimental creep strain and creep rate,and it can accurately predict long-term shear creep deformation.Based on an improved logarithmic model,the variations in creep parameters with shear stress and bedding orientations are discussed.With abovementioned findings,a mathematical method for constructing an anisotropic shear creep model of shale is proposed,which can characterize the nonlinear dependence of the anisotropic shear creep behavior of shale on the bedding orientation.

An elasto-plastic constitutive model for soft rock considering mobilization of strength

A new elasto-plastic constitutive model is presented in the framework of plasticity theory. The strength characteristics of a diatomaceous soft rock is investigated. The friction angle and cohesion of soft rock are mobilized as a function of plastic strain. A hyperbolic hardening function for the mobilized friction and a mixed parabolic and exponential equation for the mobilized cohesion are proposed. In view of the unified strength theory and the mobilizations of strength components, a yield function is given. A plastic potential function is determined by using the non-associated plastic flow rule. An elasto-plastic constitutive model is developed and verified. The results indicate that the proposed model can predict the behavior of soft rock accurately. The advantages of the proposed constitutive model are analyzed. The evidences support that the proposed constitutive model is a mixed hardening/softening model. A hump hardening/softening function for mobilized friction is extended to a more generalized condition.

A bounding surface visco-plasticity model considering generalized spacing ratio of soils

The non-unique critical state of soils with time-dependent behaviors is a significant issue in geotechnical engineering problems.However,previous bounding surface plasticity models cannot predict accurately the non-unique critical state of soils,because the distance between the compression line and critical state line charged by strain-rate effect is basically neglected.To fill this gap,a generalized spacing ratio of soils is defined in the elasto-viscoplastic framework,and a bounding surface visco-plasticity model is formulated and verified,which can consider the generalized spacing ratio.Specifically,the generalized spacing ratio of soils reflects the distance between the compression line and the critical state line of soils with time-dependent behaviors.Then,the generalized spacing ratio is introduced into an improved anisotropic bounding surface.A new expression of the visco-plastic multiplier is derived by solving the consistency equation of an anisotropic bounding surface.In the expression,a strain rate index is proposed to account for the strain-rate effect on visco-plastic strain increment,and a visco-plastic hardening modulus is derived to predict the visco-plastic response of soils in overconsolidation conditions.The model is then verified through constant strain rate tests and creep tests.Notably,it can capture the non-unique critical states of soils with time-dependent behaviors due to the generalized spacing ratio and the creep rupture of soils due to the visco-plastic multiplier that considers the stress ratio and visco-plastic strain rate.

A damage constitutive model of rock-like materials containing a single crack under the action of chemical corrosion and uniaxial compression

To describe the deformation and strength characteristics of the corroded rock-like specimens containing a single crack under uniaxial compression,a damage constitutive model combining hydro-chemical damage with coupling damage of micro-flaws and macro-cracks is proposed.Firstly,based on phenomenological theory,the damage variable of the rock-like specimens subjected to water environment erosion and chemical corrosion is obtained.Secondly,a coupled damage variable for cracked rock-like specimens is derived based on the Lemaitre strain equivalence hypothesis,which combines the Weibull statistical damage model for micro-flaws and the fracture mechanics model for a macro single crack.Then,considering the residual strength characteristics of the rock-like materials,the damage variable is modified by introducing the correction coefficient,and the damage constitutive model of the corroded rock-like specimens with a single crack under uniaxial compression is established.The model is verified by comparing the experimental stress−strain curves,and the results are in good agreement with those provided in the literature.Finally,the correction coefficient of the damage variable proposed in this paper is discussed.The damage constitutive model developed in this paper provides an effective method to describe the stress−strain relationship and residual strength characteristics of the corroded rock-like specimens with a single crack under uniaxial compression.

Recent developments of generalized plasticity models for saturated and unsaturated soils

Soil undergoes both elastic and plastic deformations under different loading conditions. A relatively accurate constitutive model of soil behaviors should be capable of predicting the elastic and plastic deformations properly. Among a large number of elastoplastic constitutive models developed over the last several decades, constitutive models based on generalized plasticity have been successfully utilized in modeling the mechanical behavior of various soils. This paper attempts to present a review of the most recent developments of generalized plasticity models for geotechnical problems. After a brief review of generalized plasticity theories and constitutive models, limitations of the original Pastor-Zienkiewicz model in practical application are summarized. Afterwards, recent achievements in the generalized plasticity models for both saturated and unsaturated soils and their applicability are analyzed, and a general approach for modification of generalized plasticity models is highlighted.

An elasto-plastic constitutive model incorporating strain softening and dilatancy for interface thin-layer element and its verification

The mechanical behaviors of the interface between coarse-grained soil and concrete were investigated by simple shear tests under condition of mixed soil slurry (bentonite mixed with cement grout).For comparison,the interfaces both without slurry and with bentonite slurry were analyzed.The experimental results show that different slurries exert much influence on the strength and deformation of soil/structure interface.Under mixed soil slurry,strain softening and shear dilatation are observed,while shear dilatation appears under the small normal stress of the interface without slurry,and shear contraction is significant under the condition of the bentonite slurry.The thickness of the interface was determined by analyzing the disturbed height of the sample with both simple shear test and particle flow code (PFC).An elasto-plastic constitutive model incorporating strain softening and dilatancy for thin layer element of interface was formulated in the framework of generalized potential theory.The relation curves of shear stress and shear strain,as well as the relation curves of normal strain and shear strain,were fitted by a piecewise function composed by hyperbolic functions and resembling normal functions.The entire model parameters can be identified by tests.The new model is verified by comparing the measured data of indoor cut-off wall model tests with the predictions from finite element method (FEM).The FEM results indicate that the stress of wall calculated by using Goodman element is too large,and the maximum deviation between the test data and prediction is about 45%.While the prediction from the proposed model is close to the measured data,and the error is generally less than 10%.

Principle of Interaction between Plastic Volumetric and Shear Strains and the Constitutive Model for Geotechnical Materials

Here is proposed the principle of interaction between plastic volumetric and shear strains, revealing the main origin of generating the complexity and variety of deformations for geotechnical materials. Here are also explained the manners of the interaction between plastic volumetric and shear strains and the conditions of generating shear dilatancy. It is demonstrated that dependency of the stress path exists and is a combination of effects of this interaction. According to this principle, it is theoretically proved that the space critical state line exists, and is unique and independent of the stress history. Based on this principle, the constitutive models that are able completely and accurately to characterize the basic behavior features for geotechnical materials have been constructed within the framework of thermodynamics. What is determined is a general expression of the constitutive relation as well as the inequality of the dissipative potential increment for obeying the second law of thermodynamics.

THEORETICAL MODEL OF EFFECTIVE STRESS COEFFICIENT FOR ROCK/SOIL-LIKE POROUS MATERIALS

Physical mechanisms and influencing factors on the effective stress coefficient for rock/soil-like porous materials are investigated, based on which equivalent connectivity index is proposed. The equivalent connectivity index, relying on the meso-scale structure of porous material and the property of liquid, denotes the connectivity of pores in Representative Element Area （REA）. If the conductivity of the porous material is anisotropic, the equivalent connectivity index is a second order tensor. Based on the basic theories of continuous mechanics and tensor analysis, relationship between area porosity and volumetric porosity of porous materials is deduced. Then a generalized expression, describing the relation between effective stress coefficient tensor and equivalent connectivity tensor of pores, is proposed, and the expression can be applied to isotropic media and also to anisotropic materials. Furthermore, evolution of porosity and equivalent connectivity index of the pore are studied in the strain space, and the method to determine the corresponding functions in expressions above is proposed using genetic algorithm and genetic programming. Two applications show that the results obtained by the method in this paper perfectly agree with the test data. This paper provides an important theoretical support to the coupled hydro-mechanical research.

Structural bonding-breakage constitutive model for natural unsaturated clayey soils

The natural clayey soils are usually structural and unsaturated,which makes their mechanical properties quite different from the remolded saturated soils.A structural constitutive model is proposed to simulate the bonding-breakage micro-mechanism.In this model,the unsaturated soil element is divided into a cementation element and a friction element according to the binary medium theory,and the stress-strain coordination for these two elements is obtained. The cementation element is regarded as elastic,whereas the friction element is regarded as elastoplastic which can be described with the Gallipoli＇s model.The theoretical formulation is verified with the comparative experiments of isotropic compressions on the saturated and unsaturated structural soils.Parametric analyses of the effects of damage variables on the model predictions are further carried out,which show that breakage deformation of natural clayey soils increases with the rising amount of initial defects.

Joints in unsaturated rocks:Thermo-hydro-mechanical formulation and constitutive behaviour

A formulation for the coupled analysis of thermo-hydro-mechanical （THM） problems in joints is first presented. The work involves the establishment of equilibrium and mass and energy balance equations. Balance equations were formulated taking into account two phases： water and air. The joint element developed was implemented in a general purpose finite element computer code for THM analysis of porous media （Code_Bright）. The program was then used to study a number of cases ranging from laboratory tests to large scale in situ tests. A numerical simulation of coupled hydraulic shear tests of rough granite joints is first presented. The tests as well as the model show the coupling between permeability and the deformation of thejoints. The experimental investigation was focused on the effects of suction on the mechanical behaviour of rock joints. Laboratory tests were performed in a direct shear cell equipped with suction control. Suction was imposed using a vapour forced convection circuit connected to the cell and controlled by an air pump. Artificial joints of Lilla claystone were prepared.Joint roughness of varying intensity was created by carving the surfaces in contact in such a manner that rock ridges of different tip angles were formed. These angles ranged from 0° （smooth joint） to 45° （very rough joint profile）. The geometric profiles of the two surfaces in contact were initially positioned in a ＂matching＂ situation. Several tests were performed for different values of suctions （200, 100, and 20 MPa） and for different values of vertical stresses （30, 60, and 150 kPa）. A constitutive model including the effects of suction and joint roughness is proposed to simulate the unsaturated behaviour of rock joints. The new constitutive law was incorporated in the code and experimental results were numerically simulated.