Discrete Element Modelling of Damage Evolution of Concrete Considering Meso-Structure of ITZ

The mechanical properties of interfacial transition zones(ITZs)have traditionally been simplified by reducing the stiffness of cement in previous simulation methods.A novel approach based on the discrete element method(DEM)has been developed for modeling concrete.This new approach efficiently simulates the meso-structure of ITZs,accurately capturing their heterogeneous properties.Validation against established uniaxial compression experiments confirms the precision of thismodel.The proposedmodel canmodel the process of damage evolution containing cracks initiation,propagation and penetration.Under increasing loads,cracks within ITZs progressively accumulate,culminating in macroscopic fractures that traverse themortarmatrix,forming the complex,serpentine path of cracks.This study reveals four distinct displacement patterns:tensile compliant,tensile opposite,mixed tensile-shear,and shear opposite patterns,each indicative of different stages in concrete’s damage evolution.The widening angle of these patterns delineates the progression of cracks,with the tensile compliant pattern signaling the initial crack appearance and the shear opposite pattern indicating the concrete model’s ultimate failure.

Damage evolution of rock-encased-backfill structure under stepwise cyclic triaxial loading

Rock-encased-backfill(RB)structures are common in underground mining,for example in the cut-andfill and stoping methods.To understand the effects of cyclic excavation and blasting activities on the damage of these RB structures,a series of triaxial stepwise-increasing-amplitude cyclic loading experiments was conducted with cylindrical RB specimens(rock on outside,backfill on inside)with different volume fractions of rock(VF=0.48,0.61,0.73,and 0.84),confining pressures(0,6,9,and 12 MPa),and cyclic loading rates(200,300,400,and 500 N/s).The damage evolution and meso-crack formation during the cyclic tests were analyzed with results from stress-strain hysteresis loops,acoustic emission events,and post-failure X-ray 3D fracture morphology.The results showed significant differences between cyclic and monotonic loadings of RB specimens,particularly with regard to the generation of shear microcracks,the development of stress memory and strain hardening,and the contact forces and associated friction that develops along the rock-backfill interface.One important finding is that as a function of the number of cycles,the elastic strain increases linearly and the dissipated energy increases exponentially.Also,compared with monotonic loading,the cyclic strain hardening characteristics are more sensitive to rising confining pressures during the initial compaction stage.Another finding is that compared with monotonic loading,more shear microcracks are generated during every reloading stage,but these microcracks tend to be dispersed and lessen the likelihood of large shear fracture formation.The transition from elastic to plastic behavior varies depending on the parameters of each test(confinement,volume fraction,and cyclic rate),and an interesting finding was that the transformation to plastic behavior is significantly lower under the conditions of 0.73 rock volume fraction,400 N/s cyclic loading rate,and 9 MPa confinement.All the findings have important practical implications on the ability of backfill to support underground excavations.

Effect of neutral polymeric bonding agent on tensile mechanical properties and damage evolution of NEPE propellant

Introducing Neutral Polymeric bonding agents(NPBA) into the Nitrate Ester Plasticized Polyether(NEPE)propellant could improve the adhesion between filler/matrix interface, thereby contributing to the development of new generations of the NEPE propellant with better mechanical properties. Therefore,understanding the effects of NPBA on the deformation and damage evolution of the NEPE propellant is fundamental to material design and applications. This paper studies the uniaxial tensile and stress relaxation responses of the NEPE propellant with different amounts of NPBA. The damage evolution in terms of interface debonding is further investigated using a cohesive-zone model(CZM). Experimental results show that the initial modulus and strength of the NEPE propellant increase with the increasing amount of NPBA while the elongation decreases. Meanwhile, the relaxation rate slows down and a higher long-term equilibrium modulus is reached. Experimental and numerical analyses indicate that interface debonding and crack propagation along filler-matrix interface are the dominant damage mechanism for the samples with a low amount of NPBA, while damage localization and crack advancement through the matrix are predominant for the ones with a high amount of NPBA. Finally, crosslinking density tests and simulation results also show that the effect of the bonding agent is interfacial rather than due to the overall crosslinking density change of the binder.

Crack propagation and damage evolution of metallic cylindrical shells under internal explosive loading

This paper investigates the three-dimensional crack propagation and damage evolution process of metallic column shells under internal explosive loading.The calibration of four typical failure parameters for 40CrMnSiB steel was conducted through experiments and subsequently applied to simulations.The numerical simulation results employing the four failure criteria were compared with the differences and similarities observed in freeze-recovery tests and ultra-high-speed tests.This analysis addressed the critical issue of determining failure criteria for the fracture of a metal shell under internal explosive loads.Building upon this foundation,the damage parameter D_(c),linked to the cumulative crack density,was defined based on the evolution characteristics of a substantial number of cracks.The relationship between the damage parameter and crack velocity over time was established,and the influence of the internal central pressure on the damage parameter and crack velocity was investigated.Variations in the fracture modes were found under different failure criteria,with the principal strain failure criterion proving to be the most effective for simulating 3D crack propagation in a pure shear fracture mode.Through statistical analysis of the shell penetration fracture radius data,it was determined that the fracture radius remained essentially constant during the crack evolution process and could be considered a constant.The propagation velocity of axial cracks ranged between 5300 m/s and 12600 m/s,surpassing the Rayleigh wave velocity of the shell material and decreasing linearly with time.The increase in shell damage exhibited an initial rapid phase,followed by deceleration,demonstrating accelerated damage during the propagation stage of the blast wave and decelerated damage after the arrival of the rarefaction wave.This study provides an effective approach for investigating crack propagation and damage evolution.The derived crack propagation and damage evolution law serves as a valuable reference for the development of crack velocity theory and the construction of shell damage evolution modes.

Micro-macro evolution of mechanical behaviors of thermally damaged rock:A state-of-the-art review

The influence of thermal damage on macroscopic and microscopic characteristics of different rocks has received much attention in the field of rock engineering.When the rocks are subjected to thermal treatment,the change of macroscopic characteristics and evolution of micro-structure would be induced,ultimately resulting in different degrees of thermal damage in rocks.To better understand the thermal damage mechanism of different rocks and its effect on the rock performance,this study reviews a large number of test results of rock specimens experiencing heating and cooling treatment in the laboratory.Firstly,the variations of macroscopic behaviors,including physical parameters,mechanical parameters,thermal conductivity and permeability,are examined.The variations of mechanical parameters with thermal treatment variables(i.e.temperature or the number of thermal cycles)are divided into four types.Secondly,several measuring methods for microstructure,such as polarizing microscopy,fluorescent method,scanning electron microscopy(SEM),X-ray computerized tomography(CT),acoustic emission(AE)and ultrasonic technique,are introduced.Furthermore,the effect of thermal damage on the mechanical parameters of rocks in response to different thermal treatments,involving temperature magnitude,cooling method and thermal cycle,are discussed.Finally,the limitations and prospects for the research of rock thermal damage are proposed.

A vector sum analysis method for stability evolution of expansive soil slope considering shear zone damage softening

Slope stability analysis is a classical mechanical problem in geotechnical engineering and engineering geology.It is of great significance to study the stability evolution of expansive soil slopes for engineering construction in expansive soil areas.Most of the existing studies evaluate the slope stability by analyzing the limit equilibrium state of the slope,and the analysis method for the stability evolution considering the damage softening of the shear zone is lacking.In this study,the large deformation shear mechanical behavior of expansive soil was investigated by ring shear test.The damage softening characteristic of expansive soil in the shear zone was analyzed,and a shear damage model reflecting the damage softening behavior of expansive soil was derived based on the damage theory.Finally,by skillfully combining the vector sum method and the shear damage model,an analysis method for the stability evolution of the expansive soil slope considering the shear zone damage softening was proposed.The results show that the shear zone subjected to large displacement shear deformation exhibits an obvious damage softening phenomenon.The damage variable equation based on the logistic function can be well used to describe the shear damage characteristics of expansive soil,and the proposed shear damage model is in good agreement with the ring shear test results.The vector sum method considering the damage softening behavior of the shear zone can be well applied to analyze the stability evolution characteristics of the expansive soil slope.The stability factor of the expansive soil slope decreases with the increase of shear displacement,showing an obvious progressive failure behavior.

High-speed penetration of ogive-nose projectiles into thick concrete targets:Tests and a projectile nose evolution model

The majority of the projectiles used in the hypersonic penetration study are solid flat-nosed cylindrical projectiles with a diameter of less than 20 mm.This study aims to fill the gap in the experimental and analytical study of the evolution of the nose shape of larger hollow projectiles under hypersonic penetration.In the hypersonic penetration test,eight ogive-nose AerMet100 steel projectiles with a diameter of 40 mm were launched to hit concrete targets with impact velocities that ranged from 1351 to 1877 m/s.Severe erosion of the projectiles was observed during high-speed penetration of heterogeneous targets,and apparent localized mushrooming occurred in the front nose of recovered projectiles.By examining the damage to projectiles,a linear relationship was found between the relative length reduction rate and the initial kinetic energy of projectiles in different penetration tests.Furthermore,microscopic analysis revealed the forming mechanism of the localized mushrooming phenomenon for eroding penetration,i.e.,material spall erosion abrasion mechanism,material flow and redistribution abrasion mechanism and localized radial upsetting deformation mechanism.Finally,a model of highspeed penetration that included erosion was established on the basis of a model of the evolution of the projectile nose that considers radial upsetting;the model was validated by test data from the literature and the present study.Depending upon the impact velocity,v0,the projectile nose may behave as undistorted,radially distorted or hemispherical.Due to the effects of abrasion of the projectile and enhancement of radial upsetting on the duration and amplitude of the secondary rising segment in the pulse shape of projectile deceleration,the predicted DOP had an upper limit.

A methodology for damage evaluation of underground tunnels subjected to static loading using numerical modeling

We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensities.To account for these variations,we utilized a Monte Carlo Simulation(MCS)technique coupled with the finite difference code FLAC^(3D),to conduct two thousand seven hundred numerical simulations of a horseshoe tunnel located within a rock mass with different geological strength index system(GSIs)and subjected to different states of static loading.To quantify the severity of damage within the rock mass,we selected one stress-based(brittle shear ratio(BSR))and one strain-based failure criterion(plastic damage index(PDI)).Based on these criteria,we then developed fragility curves.Additionally,we used mathematical approximation techniques to produce vulnerability functions that relate the probabilities of various damage states to loading intensities for different quality classes of blocky rock mass.The results indicated that the fragility curves we obtained could accurately depict the evolution of the inner and outer shell damage around the tunnel.Therefore,we have provided engineers with a tool that can predict levels of damages associated with different failure mechanisms based on variations in rock mass quality and in situ stress state.Our method is a numerically developed,multi-variate approach that can aid engineers in making informed decisions about the robustness of underground tunnels.

Modeling time-dependent mechanical behavior of hard rock considering excavation-induced damage and complex 3D stress states

To investigate the long-term stability of deep rocks,a three-dimensional(3D)time-dependent model that accounts for excavation-induced damage and complex stress state is developed.This model comprises three main components:a 3D viscoplastic isotropic constitutive relation that considers excavation damage and complex stress state,a quantitative relationship between critical irreversible deformation and complex stress state,and evolution characteristics of strength parameters.The proposed model is implemented in a self-developed numerical code,i.e.CASRock.The reliability of the model is validated through experiments.It is indicated that the time-dependent fracturing potential index(xTFPI)at a given time during the attenuation creep stage shows a negative correlation with the extent of excavationinduced damage.The time-dependent fracturing process of rock demonstrates a distinct interval effect of the intermediate principal stress,thereby highlighting the 3D stress-dependent characteristic of the model.Finally,the influence of excavation-induced damage and intermediate principal stress on the time-dependent fracturing characteristics of the surrounding rocks around the tunnel is discussed.

A whole process damage constitutive model for layered sandstone under uniaxial compression based on Logistic function

Bedding structural planes significantly influence the mechanical properties and stability of engineering rock masses.This study conducts uniaxial compression tests on layered sandstone with various bedding angles(0°,15°,30°,45°,60°,75°and 90°)to explore the impact of bedding angle on the deformational mechanical response,failure mode,and damage evolution processes of rocks.It develops a damage model based on the Logistic equation derived from the modulus’s degradation considering the combined effect of the sandstone bedding dip angle and load.This model is employed to study the damage accumulation state and its evolution within the layered rock mass.This research also introduces a piecewise constitutive model that considers the initial compaction characteristics to simulate the whole deformation process of layered sandstone under uniaxial compression.The results revealed that as the bedding angle increases from 0°to 90°,the uniaxial compressive strength and elastic modulus of layered sandstone significantly decrease,slightly increase,and then decline again.The corresponding failure modes transition from splitting tensile failure to slipping shear failure and back to splitting tensile failure.As indicated by the modulus’s degradation,the damage characteristics can be categorized into four stages:initial no damage,damage initiation,damage acceleration,and damage deceleration termination.The theoretical damage model based on the Logistic equation effectively simulates and predicts the entire damage evolution process.Moreover,the theoretical constitutive model curves closely align with the actual stress−strain curves of layered sandstone under uniaxial compression.The introduced constitutive model is concise,with fewer parameters,a straightforward parameter determination process,and a clear physical interpretation.This study offers valuable insights into the theory of layered rock mechanics and holds implications for ensuring the safety of rock engineering.

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.

Novel damage constitutive models and new quantitative identification method for stress thresholds of rocks under uniaxial compression

Four key stress thresholds exist in the compression process of rocks,i.e.,crack closure stress(σ_(cc)),crack initiation stress(σ_(ci)),crack damage stress(σ_(cd))and compressive strength(σ_(c)).The quantitative identifications of the first three stress thresholds are of great significance for characterizing the microcrack growth and damage evolution of rocks under compression.In this paper,a new method based on damage constitutive model is proposed to quantitatively measure the stress thresholds of rocks.Firstly,two different damage constitutive models were constructed based on acoustic emission(AE)counts and Weibull distribution function considering the compaction stages of the rock and the bearing capacity of the damage element.Then,the accumulative AE counts method(ACLM),AE count rate method(CRM)and constitutive model method(CMM)were introduced to determine the stress thresholds of rocks.Finally,the stress thresholds of 9 different rocks were identified by ACLM,CRM,and CMM.The results show that the theoretical stress−strain curves obtained from the two damage constitutive models are in good agreement with that of the experimental data,and the differences between the two damage constitutive models mainly come from the evolutionary differences of the damage variables.The results of the stress thresholds identified by the CMM are in good agreement with those identified by the AE methods,i.e.,ACLM and CRM.Therefore,the proposed CMM can be used to determine the stress thresholds of rocks.

An extended micromechanical-based plastic damage model for understanding water effects on quasi-brittle rocks

Water effects on the mechanical properties of rocks have been extensively investigated through experiments and numerical models.However,few studies have established a comprehensive link between the microscopic mechanisms of water-related micro-crack and the constitutive behaviors of rocks.In this work,we shall propose an extended micromechanical-based plastic damage model for understanding weakening effect induced by the presence of water between micro-crack’s surfaces on quasi-brittle rocks,based on the Mori-Tanaka homogenization and irreversible thermodynamics framework.Regarding the physical mechanism,water strengthens micro-crack propagation,which induces damage evolution during the pre-and post-stage,and weakens the elastic effective properties of rock matrix.After proposing a special calibration procedure for the determination of model parameters based on the laboratory compression tests,the proposed micromechanical-based model is verified by comparing the model predictions to the experimental results.The model effectively captures the mechanical behaviors of quasibrittle rocks subjected to the weakening effects of water.

Finite element model updating for structural damage detection using transmissibility data

This paper presents a new finite element model updating method for estimating structural parameters and detecting structural damage location and severity based on the structural responses(output-only data).The method uses the sensitivity relation of transmissibility data through a least-squares algorithm and appropriate normalization of the extracted equations.The proposed transmissibility-based sensitivity equation produces a more significant number of equations than the sensitivity equations based on the frequency response function(FRF),which can estimate the structural parameters with higher accuracy.The abilities of the proposed method are assessed by using numerical data of a two-story two-bay frame model and a plate structure model.In evaluating different damage cases,the number,location,and stiffness reduction of the damaged elements and the severity of the simulated damage have been accurately identified.The reliability and stability of the presented method against measurement and modeling errors are examined using error-contaminated data.The parameter estimation results prove the method’s capabilities as an accurate model updating algorithm.

Extension of Flow Behaviour and Damage Models for Cast Iron Alloys with Strain Rate Effect

Cast iron alloys with low production cost and quite good mechanical properties are widely used in the automotive industry.To study the mechanical behavior of a typical ductile cast iron(GJS-450)with nodular graphite,uni-axial quasi-static and dynamic tensile tests at strain rates of 10^(-4),1,10,100,and 250 s^(-1)were carried out.In order to investigate the influence of stress state on the deformation and fracture parameters,specimens with various geometries were used in the experiments.Stress strain curves and fracture strains of the GJS-450 alloy in the strain rate range of 10^(-4)to 250 s^(-1)were obtained.A strain rate-dependent plastic flow model was proposed to describe the mechanical behavior in the corresponding strain-rate range.The available damage model was extended to take the strain rate into account and calibrated based on the analysis of local fracture strains.Simulations with the proposed plastic flow model and the damage model were conducted to observe the deformation and fracture process.The results show that the strain rate has obviously nonlinear effects on the yield stress and fracture strain of GJS-450 alloys.The predictions with the proposed plastic flow and damage models at various strain rates agree well with the experimental results,which illustrates that the rate-dependent plastic flow and damage models can be used to describe the mechanical behavior of cast iron alloys at elevated strain rates.The proposed plastic flow and damage models can be used to describe the deformation and fracture analysis of materials with similar properties.

Constitutive model of viscoelastic dynamic damage for the material of gas obturator in modular-charge howitzer

In order to investigate the mechanical response behavior of the gas obturator of the breech mechanism,made of polychloroprene rubber(PCR), uniaxial compression experiments were carried out by using a universal testing machine and a split Hopkinson pressure bar(SHPB), obtaining stress-strain responses at different temperatures and strain rates. The results revealed that, in comparison to other polymers, the gas obturator material exhibited inconspicuous strain softening and hardening effects;meanwhile, the mechanical response was more affected by the strain rate than by temperature. Subsequently, a succinct viscoelastic damage constitutive model was developed based on the ZWT model, including ten undetermined parameters, formulated with incorporating three parallel components to capture the viscoelastic response at high strain rate and further enhanced by integrating a three-parameter Weibull function to describe the damage. Compared to the ZWT model, the modified model could effectively describe the mechanical response behavior of the gas obturator material at high strain rates. This research laid a theoretical foundation for further investigation into the influence of chamber sealing issues on artillery firing.

Mechanical behavior and damage constitutive model of sandstone under hydro-mechanical (H-M) coupling

Underground engineering often passes through water-rich fractured rock masses, which are prone to fracture and instability under the long-term coupling of in-situ stress field and pore water(P-W) pressure, ultimately threatening the stability of underground structures. In order to explore the mechanical properties of rocks under H-M coupling, the corresponding damage constitutive(D-C) model has become the focus of attention. Considering the inadequacy of the current research on rock strength parameters,energy evolution characteristics and D-C model under H-M coupling, the mechanical properties of typical sandstone samples are discussed based on laboratory tests. The results show that the variation of characteristic stresses of sandstone under H-M coupling conforms to the normalized attenuation equation and Mohr-Coulomb(M-C) criterion. The P-W pressure mechanism of sandstone exhibits a dynamic change from softening effect to H-M fracturing effect. The closure stress is mainly provided by cohesive strength, while the initiation stress, damage stress, and peak stress are jointly dominated by cohesive strength and friction strength. In addition, residual stress is attributed to the friction strength formed by the bite of the fracture surface. Subsequently, the energy evolution characteristics of sandstone under H-M coupling were studied, and it was found that P-W pressure weakened the energy storage capacity and energy dissipation capacity of sandstone, and H-M fracturing was an important factor in reducing its energy storage efficiency. Finally, combined with energy dissipation theory and statistical damage theory, two types of D-C models considering P-W pressure are proposed accordingly, and the model parameters can be determined by four methods. The application results indicate that the proposed and modified D-C models have high reliability, and can characterize the mechanical behavior of sandstone under H-M coupling, overcome the inconvenience of existing D-C models due to excessive mechanical parameters,and can be applied to the full-range stress–strain process. The results are conducive to revealing the deformation and damage mechanisms of rocks under H-M coupling, and can provide theoretical guidance for related engineering problems.

A creep model for ultra-deep salt rock considering thermal-mechanical damage under triaxial stress conditions

To investigate the specific creep behavior of ultra-deep buried salt during oil and gas exploitation,a set of triaxial creep experiments was conducted at elevated temperatures with constant axial pressure and unloading confining pressure conditions.Experimental results show that the salt sample deforms more significantly with the increase of applied temperature and deviatoric loading.The accelerated creep phase is not occurring until the applied temperature reaches 130℃,and higher temperature is beneficial to the occurrence of accelerated creep.To describe the specific creep behavior,a novel three-dimensional(3D)creep constitutive model is developed that incorporates the thermal and mechanical variables into mechanical elements.Subsequently,the standard particle swarm optimization(SPSO)method is adopted to fit the experimental data,and the sensibility of key model parameters is analyzed to further illustrate the model function.As a result,the model can accurately predict the creep behavior of salt under the coupled thermo-mechanical effect in deep-buried condition.Based on the research results,the creep mechanical behavior of wellbore shrinkage is predicted in deep drilling projects crossing salt layer,which has practical implications for deep rock mechanics problems.

Regional Climate Damage Quantification and Its Impacts on Future Emission Pathways Using the RICE Model

This study quantified the regional damages resulting from temperature and sea level changes using the Regional Integrated of Climate and Economy(RICE)model,as well as the effects of enabling and disabling the climate impact module on future emission pathways.Results highlight varied damages depending on regional economic development and locations.Specifically,China and Africa could suffer the most serious comprehensive damages caused by temperature change and sea level rise,followed by India,other developing Asian countries(OthAsia),and other high-income countries(OHI).The comprehensive damage fractions for China and Africa are projected to be 15.1%and 12.5%of gross domestic product(GDP)in 2195,with corresponding cumulative damages of 124.0 trillion and 87.3 trillion United States dollars(USD)from 2005 to 2195,respectively.Meanwhile,the comprehensive damage fractions in Japan,Eurasia,and Russia are smaller and projected to be lower than 5.6%of GDP in 2195,with cumulative damages of 6.8 trillion,4.2 trillion,and 3.3 trillion USD,respectively.Additionally,coastal regions like Africa,the European Union(EU),and OHI show comparable damages for sea level rise and temperature change.In China,however,sea level-induced damages are projected to exceed those from temperature changes.Moreover,this study indicates that switching the damage modules on or off affects the regional and global emission trajectories,but the magnitude is relatively small.By 2195,global emissions under the experiments with all of the damage modules switched off,only the sea level damage module switched on,and only the temperature damage module switched on,were 3.5%,2.3%and 1.2%higher than those with all of the damage modules switched on,respectively.

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.