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.

Mechanical behaviors of backfill-rock composites: Physical shear test and back-analysis

The shear behavior of backfill-rock composites is crucial for mine safety and the management of surface subsidence.For exposing the shear failure mechanism of backfill-rock composites,we conducted shear tests on backfill-rock composites under three constant normal loads,compared with the unfilled rock.To investigate the macro-and meso-failure characteristics of the samples in the shear tests,the cracking behavior of samples was recorded by a high-speed camera and acoustic emission monitoring.In parallel with the experimental test,the numerical models of backfill-rock composites and unfilled rock were established using the discrete element method to analyze the continuous-discontinuous shearing process.Based on the damage mechanics and statistics,a novel shear constitutive model was proposed to describe mechanical behavior.The results show that backfill-rock composites had a special bimodal phenomenon of shearing load-deformation curve,i.e.the first shearing peak corresponded to rock break and the second shearing peak induced by the broken of aeolian sand-cement/fly ash paste backfill.Moreover,the shearing characteristic curves of the backfill-rock composites could be roughly divided into four stages,i.e.the shear failure of the specimens experienced:stage I:stress concentration;stage II:crack propagation;stage III:crack coalescence;stage IV:shearing friction.The numerical simulation shows that the existence of aeolian sand-cement/fly ash paste backfill inevitably altered the coalescence type and failure mode of the specimens and had a strengthening effect on the shear strength of backfillrock composites.Based on damage mechanics and statistics,a shear constitutive model was proposed to describe the shear fracture characteristics of specimens,especially the bimodal phenomenon.Finally,the micro-and meso-mechanisms of shear failure were discussed by combining the micro-test and numerical results.The research can advance the better understanding of the shear behavior of backfill-rock composites and contribute to the safety of mining engineering.

Inverse design of mechanical metamaterial achieving a prescribedconstitutive curve

Besides exhibiting excellent capabilities such as energy absorption,phase-transforming metamaterials offer a vast design space for achieving nonlinear constitutive relations.This is facilitated by switching between different patterns under deformation.However,the related inverse design problem is quite challenging,due to the lack of appropriate mathematical formulation and the convergence issue in the post-buckling analysis of intermediate designs.In this work,periodic unit cells are explicitly described by the moving morphable voids method and effectively analyzed by eliminating the degrees of freedom in void regions.Furthermore,by exploring the Pareto frontiers between error and cost,an inverse design formulation is proposed for unit cells.This formulation aims to achieve a prescribed constitutive curve and is validated through numerical examples and experimental results.The design approach presented here can be extended to the inverse design of other types of mechanical metamaterials with prescribed nonlinear effective properties.

Flow characteristics and hot workability of a typical low-alloy high-strength steel during multi-pass deformation

Heavy components of low-alloy high-strength(LAHS) steels are generally formed by multi-pass forging. It is necessary to explore the flow characteristics and hot workability of LAHS steels during the multi-pass forging process, which is beneficial to the formulation of actual processing parameters. In the study, the multi-pass hot compression experiments of a typical LAHS steel are carried out at a wide range of deformation temperatures and strain rates. It is found that the work hardening rate of the experimental material depends on deformation parameters and deformation passes, which is ascribed to the impacts of static and dynamic softening behaviors. A new model is established to describe the flow characteristics at various deformation passes. Compared to the classical Arrhenius model and modified Zerilli and Armstrong model, the newly proposed model shows higher prediction accuracy with a confidence level of 0.98565. Furthermore, the connection between power dissipation efficiency(PDE) and deformation parameters is revealed by analyzing the microstructures. The PDE cannot be utilized to reflect the efficiency of energy dissipation for microstructure evolution during the entire deformation process, but only to assess the efficiency of energy dissipation for microstructure evolution in a specific deformation parameter state.As a result, an integrated processing map is proposed to better study the hot workability of the LAHS steel, which considers the effects of instability factor(IF), PDE, and distribution and size of grains. The optimized processing parameters for the multi-pass deformation process are the deformation parameters of 1223–1318 K and 0.01–0.08 s^(-1). Complete dynamic recrystallization occurs within the optimized processing parameters with an average grain size of 18.36–42.3 μm. This study will guide the optimization of the forging process of heavy components.

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.

Thermo-hydro-mechanical (THM) coupled simulation of the land subsidence due to aquifer thermal energy storage (ATES) system in soft soils

Aquifer thermal energy storage(ATES)system has received attention for heating or cooling buildings.However,it is well known that land subsidence becomes a major environmental concern for ATES projects.Yet,the effect of temperature on land subsidence has received practically no attention in the past.This paper presents a thermo-hydro-mechanical(THM)coupled numerical study on an ATES system in Shanghai,China.Four water wells were installed for seasonal heating and cooling of an agriculture greenhouse.The target aquifer at a depth of 74e104.5 m consisted of alternating layers of sand and silty sand and was covered with clay.Groundwater level,temperature,and land subsidence data from 2015 to 2017 were collected using field monitoring instruments.Constrained by data,we constructed a field scale three-dimensional(3D)model using TOUGH(Transport of Unsaturated Groundwater and Heat)and FLAC3D(Fast Lagrangian Analysis of Continua)equipped with a thermo-elastoplastic constitutive model.The effectiveness of the numerical model was validated by field data.The model was used to reproduce groundwater flow,heat transfer,and mechanical responses in porous media over three years and capture the thermo-and pressure-induced land subsidence.The results show that the maximum thermoinduced land subsidence accounts for about 60%of the total subsidence.The thermo-induced subsidence is slightly greater in winter than that in summer,and more pronounced near the cold well area than the hot well area.This study provides some valuable guidelines for controlling land subsidence caused by ATES systems installed in soft soils.

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.

Development and Application of a Power Law Constitutive Model for Eddy Current Dampers

Eddy current dampers (ECDs) have emerged as highly desirable solutions for vibration control due to theirexceptional damping performance and durability. However, the existing constitutive models present challenges tothe widespread implementation of ECD technology, and there is limited availability of finite element analysis (FEA)software capable of accurately modeling the behavior of ECDs. This study addresses these issues by developing anewconstitutivemodel that is both easily understandable and user-friendly for FEAsoftware. By utilizing numericalresults obtained from electromagnetic FEA, a novel power law constitutive model is proposed to capture thenonlinear behavior of ECDs. The effectiveness of the power law constitutive model is validated throughmechanicalproperty tests and numerical seismic analysis. Furthermore, a detailed description of the application process ofthe power law constitutive model in ANSYS FEA software is provided. To facilitate the preliminary design ofECDs, an analytical derivation of energy dissipation and parameter optimization for ECDs under harmonicmotionis performed. The results demonstrate that the power law constitutive model serves as a viable alternative forconducting dynamic analysis using FEA and optimizing parameters for ECDs.

A macro-mesoscopic constitutive model for porous and cracked rock under true triaxial conditions

The complex mechanical and damage mechanisms of rocks are intricately tied to their diverse mineral compositions and the formation of pores and cracks under external loads.Numerous rock tests reveal a complex interplay between the closure of porous defects and the propagation of induced cracks,presenting challenges in accurately representing their mechanical properties,especially under true triaxial stress conditions.This paper proposes a conceptualization of rock at the mesoscopic level as a two-phase composite,consisting of a bonded medium matrix and frictional medium inclusions.The bonded medium is characterized as a mesoscopic elastic material,encompassing various minerals surrounding porous defects.Its mechanical properties are determined using the mixed multi-inclusion method.Transformation of the bonded medium into the frictional medium occurs through crack extension,with its elastoplastic properties defined by the DruckerePrager yield criterion,accounting for hardening,softening,and extension.MorieTanaka and Eshelby’s equivalent inclusion methods are applied to the bonded and frictional media,respectively.The macroscopic mechanical properties of the rock are derived from these mesoscopic media.Consequently,a True Triaxial Macro-Mesoscopic(TTMM)constitutive model is developed.This model effectively captures the competitive effect and accurately describes the stress-deformation characteristics of granite.Utilizing the TTMM model,the strains resulting from porous defect closure and induced crack extension are differentiated,enabling quantitative determination of the associated damage evolution.

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.

Phase-field modeling for anisotropic ductile damage of magnesium alloys at finite deformations

The damage anisotropy of an extruded ZK60 Mg alloy is characterized using tensile tests and scanning electronic microscopy.The accumulation of anisotropic deformations leads to the great differences of the dimple evolution and strains at fracture along different loading directions.To introduce the anisotropic deformation information into the damage constitutive relationship,a thermodynamically consistent phase-field model of ductile damage fully coupled with elastoplastic finite deformations is developed in this study.Using the user-defined constitutive relationship and displacement-temperature coupling element,the finite element simulations are conducted.The results show that:(1)ZK60 Mg alloys presents clear R-value difference in 0°,45°,and 90°tests of intact specimens.The 45°test possesses the greatest R-value(1.50)and the greatest strain at fracture,however,the R-value for 0°is less than 1,indicating the thinning is preferential.(2)The higher ultimate stress leads to a larger average dimension of the dimples,whereas the higher density correlates with a larger elongation ratio at the fracture.The disappearance of the stress-bearing area indicates that the phase-field assumption on stress degradation is completely compatible with the dimple analysis on fractography.(3)The simulation results of the stress-strain relationships and damage paths correlate well with the experimental ductile damage of magnesium alloys at 200◦C.Slight errors are basically attributed to the modeling parameters and finite element iteration algorithm.The proposed model presents fine applicability and reliability for the predictions of plastic deformations,ductile damage,and fracture of anisotropic Mg alloys.

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.

Data-driven modeling on anisotropic mechanical behavior of brain tissue with internal pressure

Brain tissue is one of the softest parts of the human body,composed of white matter and grey matter.The mechanical behavior of the brain tissue plays an essential role in regulating brain morphology and brain function.Besides,traumatic brain injury(TBI)and various brain diseases are also greatly influenced by the brain's mechanical properties.Whether white matter or grey matter,brain tissue contains multiscale structures composed of neurons,glial cells,fibers,blood vessels,etc.,each with different mechanical properties.As such,brain tissue exhibits complex mechanical behavior,usually with strong nonlinearity,heterogeneity,and directional dependence.Building a constitutive law for multiscale brain tissue using traditional function-based approaches can be very challenging.Instead,this paper proposes a data-driven approach to establish the desired mechanical model of brain tissue.We focus on blood vessels with internal pressure embedded in a white or grey matter matrix material to demonstrate our approach.The matrix is described by an isotropic or anisotropic nonlinear elastic model.A representative unit cell(RUC)with blood vessels is built,which is used to generate the stress-strain data under different internal blood pressure and various proportional displacement loading paths.The generated stress-strain data is then used to train a mechanical law using artificial neural networks to predict the macroscopic mechanical response of brain tissue under different internal pressures.Finally,the trained material model is implemented into finite element software to predict the mechanical behavior of a whole brain under intracranial pressure and distributed body forces.Compared with a direct numerical simulation that employs a reference material model,our proposed approach greatly reduces the computational cost and improves modeling efficiency.The predictions made by our trained model demonstrate sufficient accuracy.Specifically,we find that the level of internal blood pressure can greatly influence stress distribution and determine the possible related damage behaviors.

Experimental study of the damage characteristics of rocks containing non-penetrating cracks under cyclic loading

The damage evolution process of non-penetrating cracks often causes some unexpected engineering disasters.Gypsum specimens containing non-penetrating crack(s)are used to study the damage evolution and characteristics under cyclic loading.The results show that under cyclic loading,the relationship between the number of non-penetrating crack(s)and the characteristic parameters(cyclic number,peak stress,peak strain,failure stress,and failure strain)of the pre-cracked specimens can be represented by a decreasing linear function.The damage evolution equation is fitted by calibrating the accumulative plastic strain for each cycle,and the damage constitutive equation is proposed by the concept of effective stress.Additionally,non-penetrating cracks are more likely to cause uneven stress distribution,damage accumulation,and local failure of specimen.The local failure can change the stress distribution and relieve the inhibition of non-penetrating crack extension and eventually cause a dramatic destruction of the specimen.Therefore,the evolution process caused by non-penetrating cracks can be regarded as one of the important reasons for inducing rockburst.These results are expected to improve the understanding of the process of spalling formation and rockburst and can be used to analyze the stability of rocks or rock structures.

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.

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.

Time-dependent behavior and permeability evolution of limestone under hydro-mechanical coupling

Excessively high pore water pressure presents unpredictable risks to the safety of rock tunnels in mountainous regions that are predominantly composed of limestone. Investigating the creep characteristics and permeability evolution of limestone under varying hydrated conditions is crucial for a better understanding of the delayed deformation mechanisms of limestone rock tunnels. To this end, this paper initially conducts a series of multi-stage triaxial creep tests on limestone samples under varying pore water pressures. The experiment examines how pore water pressure affects limestone’s creep strain, strain rate, long-term strength, lifespan, and permeability, all within the context of hydraulicmechanical(HM) coupling. To better describe the creep behavior associated with pore water pressure, this paper proposes a new nonlinear fractional creep constitutive model. This constitutive model depicts the initial, steady-state, and accelerated phases of limestone’s creep behavior. Finally, the proposed model is applied to the numerical realization of deformation in limestone tunnel, validating the effectiveness of the proposed constitutive model in predicting tunnel’s creep deformation. This research enhances our understanding of limestone’s creep characteristics and permeability evolution under HM coupling, laying a foundation for assessing the longterm stability of mountain tunnels.

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.

Effect of Gd Addition on Hot Deformation Behavior and Microstructure Evolution of 7075 Aluminum Alloy

In order to clarify the effect of rare earth Gd on the microstructure evolution and deformation behavior of 7075 aluminum alloy during hot compression,uniaxial compression tests of Al-Zn-Mg-Cu-0.5%Gd were conducted at strain rates of 0.001,0.01,0.1,and 1 s^(-1)with the temperatures ranging from 350 to 450℃.The microstructural evolution during deformation was characterized using optical microscopy and electron backscatter diffraction(EBSD)techniques.The experimental results indicate that the addition of the rare earth element Gd significantly increases the peak flow stress and thermal activation energy of the alloy.Due to the pinning effect of rare earth phases,dislocation movement is hindered,leading to an increased level of work hardening in the alloy.However,the dynamic recrystallization of the alloy is complicated.At a high Z(Zener-Hollomon parameter)values,recrystallization occurs in the form of DDRX(Discontinuous Dynamic Recrystallization),making it easier to nucleate at grain boundaries.As the Z value decreases gradually,the recrystallization mechanism transitions from discontinuous dynamic recrystallization(DDRX)to continuous dynamic recrystallization(CDRX).At a low Z values with the strain rate of 0.001 s^(-1),the inhibitory effect of rare earths weakens,resulting in a comparable recrystallization ratio between Al-Zn-Mg-Cu-Gd alloy and 7075 aluminum alloy.Moreover,the average grain size of the aluminum alloy with Gd addition is only half that of 7075 aluminum.The addition of Gd provides Orowan and substructure strengthening for the alloy,which greatly improves the work-hardening of the alloy compared with 7075 aluminum alloy and improves the strength of the alloy.

Experiment and constitutive modelling of creep deformation in the frozen silt-concrete interface

To ensure the long-term safety and stability of bridge pile foundations in permafrost regions,it is necessary to investigate the rheological effects on the pile tip and pile side bearing capacities.The creep characteristics of the pile-frozen soil interface are critical for determining the long-term stability of permafrost pile foundations.This study utilized a self-developed large stress-controlled shear apparatus to investigate the shear creep characteristics of the frozen silt-concrete interface,and examined the influence of freezing temperatures(−1,−2,and−5°C),contact surface roughness(0,0.60,0.75,and 1.15 mm),normal stress(50,100,and 150 kPa),and shear stress on the creep characteristics of the contact surface.By incorporating the contact surface’s creep behavior and development trends,we established a creep constitutive model for the frozen silt-concrete interface based on the Nishihara model,introducing nonlinear elements and a damage factor.The results revealed significant creep effects on the frozen silt-concrete interface under constant load,with creep displacement at approximately 2-15 times the instantaneous displacement and a failure creep displacement ranging from 6 to 8 mm.Under different experimental conditions,the creep characteristics of the frozen silt-concrete interface varied.A larger roughness,lower freezing temperatures,and higher normal stresses resulted in a longer sample attenuation creep time,a lower steady-state creep rate,higher long-term creep strength,and stronger creep stability.Building upon the Nishihara model,we considered the influence of shear stress and time on the viscoelastic viscosity coefficient and introduced a damage factor to the viscoplasticity.The improved model effectively described the entire creep process of the frozen silt-concrete interface.The results provide theoretical support for the interaction between pile and soil in permafrost regions.