Discontinuity development patterns and the challenges for 3D discrete fracture network modeling on complicated exposed rock surfaces

Natural slopes usually display complicated exposed rock surfaces that are characterized by complex and substantial terrain undulation and ubiquitous undesirable phenomena such as vegetation cover and rockfalls.This study presents a systematic outcrop research of fracture pattern variations in a complicated rock slope,and the qualitative and quantitative study of the complex phenomena impact on threedimensional(3D)discrete fracture network(DFN)modeling.As the studies of the outcrop fracture pattern have been so far focused on local variations,thus,we put forward a statistical analysis of global variations.The entire outcrop is partitioned into several subzones,and the subzone-scale variability of fracture geometric properties is analyzed(including the orientation,the density,and the trace length).The results reveal significant variations in fracture characteristics(such as the concentrative degree,the average orientation,the density,and the trace length)among different subzones.Moreover,the density of fracture sets,which is approximately parallel to the slope surface,exhibits a notably higher value compared to other fracture sets across all subzones.To improve the accuracy of the DFN modeling,the effects of three common phenomena resulting from vegetation and rockfalls are qualitatively analyzed and the corresponding quantitative data processing solutions are proposed.Subsequently,the 3D fracture geometric parameters are determined for different areas of the high-steep rock slope in terms of the subzone dimensions.The results show significant variations in the same set of 3D fracture parameters across different regions with density differing by up to tenfold and mean trace length exhibiting differences of 3e4 times.The study results present precise geological structural information,improve modeling accuracy,and provide practical solutions for addressing complex outcrop issues.

Theoretical Modeling and Surface Roughness Prediction of Microtextured Surfaces in Ultrasonic Vibration-Assisted Milling

Textured surfaces with certain micro/nano structures have been proven to possess some advanced functions,such as reducing friction,improving wear and increasing wettability.Accurate prediction of micro/nano surface textures is of great significance for the design,fabrication and application of functional textured surfaces.In this paper,based on the kinematic analysis of cutter teeth,the discretization of ultrasonic machining process,transformation method of coordinate systems and the cubic spline data interpolation,an integrated theoretical model was established to characterize the distribution and geometric features of micro textures on the surfaces machined by different types of ultrasonic vibration-assisted milling(UVAM).Based on the theoretical model,the effect of key process parameters(vibration directions,vibration dimensions,cutting parameters and vibration parameters)on tool trajectories and microtextured surface morphology in UVAM is investigated.Besides,the effect of phase difference on the elliptical shape in 2D/3D ultrasonic elliptical vibration-assisted milling(UEVAM)was analyzed.Compared to conventional numerical models,the method of the cubic spline data interpolation is applied to the simulation of microtextured surface morphology in UVAM,which is more suitable for characterizing the morphological features of microtextured surfaces than traditional methods due to the presence of numerous micro textures.The prediction of surface roughness indicates that the magnitude of ultrasonic amplitude in z-direction should be strictly limited in 1D rotary UVAM,2D and 3D UEVAM due to the unfavorable effect of axial ultrasonic vibration on the surface quality.This study can provide theoretical guidance for the design and fabrication of microtextured surfaces in UVAM.

Relationships between Terrain Features and Forecasting Errors of Surface Wind Speeds in a Mesoscale Numerical Weather Prediction Model

Numerical weather prediction(NWP)models have always presented large forecasting errors of surface wind speeds over regions with complex terrain.In this study,surface wind forecasts from an operational NWP model,the SMS-WARR(Shanghai Meteorological Service-WRF ADAS Rapid Refresh System),are analyzed to quantitatively reveal the relationships between the forecasted surface wind speed errors and terrain features,with the intent of providing clues to better apply the NWP model to complex terrain regions.The terrain features are described by three parameters:the standard deviation of the model grid-scale orography,terrain height error of the model,and slope angle.The results show that the forecast bias has a unimodal distribution with a change in the standard deviation of orography.The minimum ME(the mean value of bias)is 1.2 m s^(-1) when the standard deviation is between 60 and 70 m.A positive correlation exists between bias and terrain height error,with the ME increasing by 10%−30%for every 200 m increase in terrain height error.The ME decreases by 65.6%when slope angle increases from(0.5°−1.5°)to larger than 3.5°for uphill winds but increases by 35.4%when the absolute value of slope angle increases from(0.5°−1.5°)to(2.5°−3.5°)for downhill winds.Several sensitivity experiments are carried out with a model output statistical(MOS)calibration model for surface wind speeds and ME(RMSE)has been reduced by 90%(30%)by introducing terrain parameters,demonstrating the value of this study.

Rainfall-triggered waste dump instability analysis based on surface 3D deformation in physical model test

Landslide is the second largest natural disaster after earthquake. It is of significance to study the evolution laws and failure mechanism of landslides based on its surface 3D deformation information. Based on the rainfall-triggered waste dump instability model test, we studied the failure mechanisms of the waste dump by integrating surface deformation and internal slope stress and proposed novel parameters for identifying landslide stability. We developed a noncontact measurement device, which can obtain millimeter-level 3D deformation data for surface scene in physical model test;Then we developed the similar materials and established a test model for a waste dump. Based on the failure characteristics of slope surface, internal stress of slope body and displacement contours during the whole process, we divided the slope instability process in model test into four stages: rainfall infiltration and surface erosion, shallow sliding, deep sliding, and overall instability. Based on the obtained surface deformation data, we calculated the volume change during slope instability process and compared it with the point displacement on slope surface. The results showed that the volume change can not only reflect the slow-ultra acceleration process of slope failure, but also fully reflect the above four stages and reduce the fluctuations caused by random factors. Finally, this paper proposed two stability identification parameters: the volume change rate above the slip surface and the relative velocity of volume change rate. According to the calculation of these two parameters in model test, they can be used for study the deformation and failure mechanism of slope stability.

A CFD Model to Evaluate Near-Surface Oil Spill from a Broken Loading Pipe in Shallow Coastal Waters

Oil spills continue to generate various issues and concerns regarding their effect and behavior in the marine environment,owing to the related potential for detrimental environmental,economic and social implications.It is essential to have a solid understanding of the ways in which oil interacts with the water and the coastal ecosystems that are located nearby.This study proposes a simplified model for predicting the plume-like transport behavior of heavy Bunker C fuel oil discharging downward from an acutely-angled broken pipeline located on the water surface.The results show that the spill overall profile is articulated in three major flow areas.The first,is the source field,i.e.,a region near the origin of the initial jet,followed by the intermediate or transport field,namely,the region where the jet oil flow transitions into an underwater oil plume flow and starts to move horizontally,and finally,the far-field,where the oil re-surface and spreads onto the shore at a significant distance from the spill site.The behavior of the oil in the intermediate field is investigated using a simplified injection-type oil spill model capable of mimicking the undersea trapping and lateral migration of an oil plume originating from a negatively buoyant jet spill.A rectangular domain with proper boundary conditions is used to implement the model.The Projection approach is used to discretize a modified version of the Navier-Stokes equations in two dimensions.A benchmark fluid flow issue is used to verify the model and the results indicate a reasonable relationship between specific gravity and depth as well as agreement with the aerial data and a vertical temperature profile plot.

Surface Defect Detection and Evaluation Method of Large Wind Turbine Blades Based on an Improved Deeplabv3+Deep Learning Model

The accumulation of defects on wind turbine blade surfaces can lead to irreversible damage,impacting the aero-dynamic performance of the blades.To address the challenge of detecting and quantifying surface defects on wind turbine blades,a blade surface defect detection and quantification method based on an improved Deeplabv3+deep learning model is proposed.Firstly,an improved method for wind turbine blade surface defect detection,utilizing Mobilenetv2 as the backbone feature extraction network,is proposed based on an original Deeplabv3+deep learning model to address the issue of limited robustness.Secondly,through integrating the concept of pre-trained weights from transfer learning and implementing a freeze training strategy,significant improvements have been made to enhance both the training speed and model training accuracy of this deep learning model.Finally,based on segmented blade surface defect images,a method for quantifying blade defects is proposed.This method combines image stitching algorithms to achieve overall quantification and risk assessment of the entire blade.Test results show that the improved Deeplabv3+deep learning model reduces training time by approximately 43.03%compared to the original model,while achieving mAP and MIoU values of 96.87%and 96.93%,respectively.Moreover,it demonstrates robustness in detecting different surface defects on blades across different back-grounds.The application of a blade surface defect quantification method enables the precise quantification of dif-ferent defects and facilitates the assessment of risk levels associated with defect measurements across the entire blade.This method enables non-contact,long-distance,high-precision detection and quantification of surface defects on the blades,providing a reference for assessing surface defects on wind turbine blades.

Comparative Analysis of Statistical Thickness Models for the Determination of the External Specific Surface and the Surface of the Micropores of Materials: The Case of a Clay Concrete Stabilized Using Sugar Cane Molasses

In this work, four empirical models of statistical thickness, namely the models of Harkins and Jura, Hasley, Carbon Black and Jaroniec, were compared in order to determine the textural properties (external surface and surface of micropores) of a clay concrete without molasses and clay concretes stabilized with 8%, 12% and 16% molasses. The results obtained show that Hasley’s model can be used to obtain the external surfaces. However, it does not allow the surface of the micropores to be obtained, and is not suitable for the case of simple clay concrete (without molasses) and for clay concretes stabilized with molasses. The Carbon Black, Jaroniec and Harkins and Jura models can be used for clay concrete and stabilized clay concrete. However, the Carbon Black model is the most relevant for clay concrete and the Harkins and Jura model is for molasses-stabilized clay concrete. These last two models augur well for future research.

Seasonal constraint of dynamic water temperature on riverine dissolved inorganic nitrogen transport in land surface modeling

水体温度变化对河流可溶性无机氮(DIN)输送有着强烈控制作用.然而,在全球尺度上河流DIN输送量对水温度变化的响应尚不清楚.因此,本文基于陆面过程模式,耦合河流水温估算和DIN传输方案,设定有,无动态水温情景,对比研究陆面模拟中水温变化对河流DIN通量变化的影响.结果表明:在考虑水温动态变化后,在30°N和30°S之间, DIN通量年振幅减小5%–25%.在中国东部地区,水温动态变化使河流DIN通量在夏季减少1%–3%,在冬季增加1%–5%,对DIN通量具有明显的季节性约束作用,表明动态水温的表达在河流DIN输送模拟中的重要性.

Elucidating Dominant Factors Affecting Land Surface Hydrological Simulations of the Community Land Model over China

In order to compare the impacts of the choice of land surface model(LSM)parameterization schemes,meteorological forcing,and land surface parameters on land surface hydrological simulations,and explore to what extent the quality can be improved,a series of experiments with different LSMs,forcing datasets,and parameter datasets concerning soil texture and land cover were conducted.Six simulations are run for the Chinese mainland on 0.1°×0.1°grids from 1979 to 2008,and the simulated monthly soil moisture(SM),evapotranspiration(ET),and snow depth(SD)are then compared and assessed against observations.The results show that the meteorological forcing is the most important factor governing output.Beyond that,SM seems to be also very sensitive to soil texture information;SD is also very sensitive to snow parameterization scheme in the LSM.The Community Land Model version 4.5(CLM4.5),driven by newly developed observation-based regional meteorological forcing and land surface parameters(referred to as CMFD_CLM4.5_NEW),significantly improved the simulations in most cases over the Chinese mainland and its eight basins.It increased the correlation coefficient values from 0.46 to 0.54 for the SM modeling and from 0.54 to 0.67 for the SD simulations,and it decreased the root-mean-square error(RMSE)from 0.093 to 0.085 for the SM simulation and reduced the normalized RMSE from 1.277 to 0.201 for the SD simulations.This study indicates that the offline LSM simulation using a refined LSM driven by newly developed observation-based regional meteorological forcing and land surface parameters can better model reginal land surface hydrological processes.

Analysis of the interaction between bolt-reinforced rock and surface support in tunnels based on convergence-confinement method

To investigate the interaction of the bolt-reinforced rock and the surface support,an analytical model of the convergence-confinement type is proposed,considering the sequential installation of the fully grouted rockbolts and the surface support.The rock mass is assumed to be elastic-brittle-plastic material,obeying the linear Mohr-Coulomb criterion or the non-linear Hoek-Brown criterion.According to the strain states of the tunnel wall at bolt and surface support installation and the relative magnitude between the bolt length and the plastic depth during the whole process,six cases are categorized upon solving the problem.Each case is divided into three stages due to the different effects of the active rockbolts and the passive surface support.The fictitious pressure is introduced to quantify the threedimensional(3D)effect of the tunnel face,and thus,the actual physical location along the tunnel axis of the analytical section can be considered.By using the bolt-rock strain compatibility and the rocksurface support displacement compatibility conditions,the solutions of longitudinal tunnel displacement and the reaction pressure of surface support along the tunnel axis are obtained.The proposed analytical solutions are validated by a series of 3D numerical simulations.Extensive parametric studies are conducted to examine the effect of the typical parameters of rockbolts and surface support on the tunnel displacement and the reaction pressure of the surface support under different rock conditions.The results show that the rockbolts are more effective in controlling the tunnel displacement than the surface support,which should be installed as soon as possible with a suitable length.For tunnels excavated in weak rocks or with restricted displacement control requirements,the surface support should also be installed or closed timely with a certain stiffness.The proposed method provides a convenient alternative approach for the optimization of rockbolts and surface support at the preliminary stage of tunnel design.

Analytical Modeling and Mechanism Analysis of Time-Varying Excitation for Surface Defects in Rolling Element Bearings

Surface defects,including dents,spalls,and cracks,for rolling element bearings are the most common faults in rotating machinery.The accurate model for the time-varying excitation is the basis for the vibration mechanism analysis and fault feature extraction.However,in conventional investigations,this issue is not well and fully addressed from the perspective of theoretical analysis and physical derivation.In this study,an improved analytical model for time-varying displacement excitations(TVDEs)caused by surface defects is theoretically formulated.First and foremost,the physical mechanism for the effect of defect sizes on the physical process of rolling element-defect interaction is revealed.According to the physical interaction mechanism between the rolling element and different types of defects,the relationship between time-varying displacement pulse and defect sizes is further analytically derived.With the obtained time-varying displacement pulse,the dynamic model for the deep groove bearings considering the internal excitation caused by the surface defect is established.The nonlinear vibration responses and fault features induced by surface defects are analyzed using the proposed TVDE model.The results suggest that the presence of surface defects may result in the occurrence of the dual-impulse phenomenon,which can serve as indexes for surface-defect fault diagnosis.

Variation in the surface heat flux on the north and south slopes of Mount Qomolangma

The distinctive conditions present on the north and south slopes of Mount Qomolangma,along with the intricate variations in the underlying surfaces,result in notable variations in the surface energy flux patterns of the two slopes.In this paper,data from TESEBS(Topographical Enhanced Surface Energy Balance System),remote sensing data from eight cloud-free scenarios,and observational data from nine stations are utilized to examine the fluctuations in the surface heat flux on both slopes.The inclusion of MCD43A3 satellite data enhances the surface albedo,contributing to more accurate simulation outcomes.The model results are validated using observational data.The RMSEs of the net radiation,ground heat,sensible heat,and latent heat flux are 40.73,17.09,33.26,and 30.91 W m^(−2),respectively.The net radiation flux is greater on the south slope and exhibits a rapid decline from summer to autumn.Due to the influence of the monsoon,on the north slope,the maximum sensible heat flux occurs in the pre-monsoon period in summer and the maximum latent heat flux occurs during the monsoon.The south slope experiences the highest latent heat flux in summer.The dominant flux on the north slope is sensible heat,while it is latent heat on the south slope.The seasonal variations in the ground heat flux are more pronounced on the south slope than on the north slope.Except in summer,the ground heat flux on the north slope surpasses that on the south slope.

BCC_AVIM陆面模式不同土壤垂直离散化方案对土壤水热输送的数值模拟

陆面模式中的土壤如何分层至关重要,一般情况下认为在受大气—地面界面影响较大的土壤表层应该细分,土壤越往深层,土壤层厚度的分层可以相应加大。已有研究表明:在不同天气和气候积分条件下,土壤厚度在根带以下,其对模拟结果有不同的影响。这表示根据不同的研究要求,应当改变土壤的分层方式。然而什么是陆面模式中的最优土壤分层方式仍不确定。因此本文主要利用BCC_AVIM陆面模式探究不同土壤垂直离散化方案对于土壤水热属性,表面辐射通量与感热、潜热通量模拟的敏感性,从而达到提升模式模拟效果的目的。BCC_AVIM陆面模式中原土壤分层方案为10层,将土壤各层的节点深度,土壤层厚度以及土壤各层的界面深度进行插值,由原来的10层土壤层插值到20层,在本文中称为方案一;参考模式CLM5.0中的土壤垂直离散化方案并改进到BCC_AVIM陆面模式中,原土壤层次也由10层增加到20层,为本文的方案二。将改进后的方案一和方案二与原方案结果对比分析发现:(1)方案一和方案二对土壤温度的模拟与实测数据更吻合,对于各层土壤温度的数值大小与变化趋势的模拟效果有所提升,其中方案一对浅层土壤温度的模拟效果更好。(2)在土壤湿度的模拟上,三种方案对浅层土壤湿度的模拟效果较好,对深层土壤湿度的模拟效果相对较差,其中方案一对土壤各层湿度曲线变化趋势与数值大小的模拟更加贴近实测数据。(3)方案一对土壤各层是否发生冻结或消融的时间判定更合理,更加接近实测数据。整体上,方案一的模拟效果最佳。由此,得出结论:方案一相较于原方案对土壤温、湿度的模拟效果有所改善,这表明在同样的土壤深度下,更密集更细致的土壤分层有利于提升模式对土壤水热输送的模拟能力。同时,方案一的模拟效果总体优于方案二,这也表明,在相同的土壤层次下,浅中层拥有更密集的土壤分层,对于提升模式对土壤水热输送的模拟能力有一定的积极影响。

Soil-water characteristic surface model of soil-rock mixture

The relationship between the water content or saturation of unsaturated soils and its matrix suction is commonly described by the soilwater characteristic curve(SWCC).Currently,study on the SWCC model is focused on fine-grained soils like clay and silty soils,but the SWCC model for grinding soil-rock mixture(SRM)is less studied.Considering that the SRM is in a certain compaction state in the actual project,this study established a surface model with three variables of coupling compaction degree-substrate suction-moisture content based on the Cavalcante-Zornberg soil-water characteristic curve model.Then,the influence of each fitting parameter on the curve was analyzed.For the common SRM,the soil-water characteristic test was conducted.Moreover,the experimental measurements exhibit remarkable consistency with the mode surface.The analysis shows that the surface model intuitively describes the soil-water characteristics of grinding SRM,which can provide the SWCC of soils with bimodal pore characteristics under specific compaction degrees.Furthermore,it can reflect the influence of compaction degrees on the SWCC of rock-soil mass and has a certain prediction effect.The SWCC of SRM with various soil-rock ratios have a double-step shape.With the increase in compaction degree,the curves as a whole tend toward decreasing mass moisture content.The curve changes are mainly concentrated in the large pore section.

Contact Angle Prediction Model for Underwater Oleophobic Surfaces Based on Multifractal Theory

Traditional microstructure scale parameters have difficulty describing the structure and distribution of a roughmaterial’s surface morphology comprehensively and quantitatively. This study constructs hydrophilic and underwateroleophobic surfaces based on polyvinylidene fluoride (PVDF) using a chemical modification method, and the fractaldimension and multifractal spectrum are used to quantitatively characterize the microscopic morphology. A new contactangle prediction model for underwater oleophobic surfaces is established. The results show that the fractal dimension ofthe PVDF surface first increases and then decreases with the reaction time. The uniformity characterized by the multifractalspectrum was generally consistent with scanning electron microscope observations. The contact angle of water droplets onthe PVDF surface is negatively correlated with the fractal dimension, and oil droplets in water are positively correlated.When the fractal dimension is 2.0975, the new contact angle prediction model has higher prediction accuracy. Themaximum and minimum relative deviations of the contact angle between the theoretical and measured data are 18.20%and 0.72%, respectively. For water ring transportation, the larger the fractal dimension and spectral width of the materialsurface, the smaller the absolute value of the spectral difference, the stronger the hydrophilic and oleophobic properties, andthe better the water ring transportation stability.

An accurate analytical surface potential model of heterojunction tunnel FET

Based on the accurate and efficient thermal injection method, we develop a fully analytical surface potential model for the heterojunction tunnel field-effect transistor(H-TFET). This model accounts for both the effects of source depletion and inversion charge, which are the key factors influencing the charge, capacitance and current in H-TFET. The accuracy of the model is validated against TCAD simulation and is greatly improved in comparison with the conventional model based on Maxwell–Boltzmann approximation. Furthermore, the dependences of the surface potential and electric field on biases are well predicted and thoroughly analyzed.

Effect of Fibre Size on Mechanical Properties and Surface Roughness of PLA Composites by Using Fused Deposition Modelling (FDM)

Naturalfibre as a reinforcing agent has been widely used in many industrial applications.Nevertheless,several factors need to be considered,such as the size and weight percentage of thefibre used in binding.Using fused deposition modelling(FDM),this factor was investigated by varying the size of naturalfibre as the responding variable with afixed weight percentage of kenaffibre.The process of modifying the naturalfibre in terms of size might increase the dispersion of kenaffibre in the polymer matrix and increase the adhesion bonding between thefibre and matrix of composites,subsequently improving the interfacial bonding between these two phases.In this paper,the effect offibre size was evaluated by performing the mechanical test,Scanning Electron Micrograph(SEM)to observe the morphology of the composites,and also by surface analysis.The surface roughness was visualised using a 3D profilometer and thefigure was illustrated as colour shading in the image.The composite withfibre size≤100μm displayed better tensile andflexural strength,compared to other sizes.In conclusion,by reducing the size of thefibre,the composites could develop high strength performance for industrial applications.

Evaluation of Nonbreaking Wave-Induced Mixing Parameterization Schemes Based on a One-Dimensional Ocean Model

Surface waves have a considerable effect on vertical mixing in the upper ocean.In the past two decades,the vertical mixing induced through nonbreaking surface waves has been used in ocean and climate models to improve the simulation of the upper ocean.Thus far,several nonbreaking wave-induced mixing parameterization schemes have been proposed;however,no quantitative comparison has been performed among them.In this paper,a one-dimensional ocean model was used to compare the performances of five schemes,including those of Qiao et al.(Q),Hu and Wang(HW),Huang and Qiao(HQ),Pleskachevsky et al.(P),and Ghantous and Babanin(GB).Similar to previous studies,all of these schemes can decrease the simulated sea surface temperature(SST),increase the subsurface temperature,and deepen the mixed layer,thereby alleviating the common thermal deviation problem of the ocean model for upper ocean simulation.Among these schemes,the HQ scheme exhibited the weakest wave-induced mixing effect,and the HW scheme exhibited the strongest effect;the other three schemes exhibited roughly the same effect.In particular,the Q and P schemes exhibited nearly the same effect.In the simulation based on observations from the Ocean Weather Station Papa,the HQ scheme exhibited the best performance,followed by the Q scheme.In the experiment with the HQ scheme,the root-mean-square deviation of the simulated SST from the observations was 0.43℃,and the mixed layer depth(MLD)was 2.0 m.As a contrast,the deviations of the SST and MLD reached 1.25℃ and 8.4 m,respectively,in the experiment without wave-induced mixing.

Effects of carrier density and interactions on pairing symmetry in a t_(2g) model

By utilizing the fluctuation exchange approximation method,we perform a study on the superconducting pairing symmetry in a t_(2g) three-orbital model on the square lattice.Although the tight-binding parameters of the model are based on Sr_(2)RuO_(4),we have systematically studied the evolution of superconducting pairing symmetry with the carrier density and interactions,making our findings relevant to a broader range of material systems.Under a moderate Hund’s coupling,we find that spin fluctuations dominate the superconducting pairing,leading to a prevalent spin-singlet pairing with a d_(x^(2)-y^(2))-wave symmetry for the carrier density within the range of n=1.5-4 per site.By reducing the Hund’s coupling,the charge fluctuations are enhanced and play a crucial role in determining the pairing symmetry,leading to a transition of the pairing symmetry from the spin-singlet d_(x^(2)-y^(2))-wave to the spin-triplet p-wave.Furthermore,we find that the superconducting pairings are orbital dependent.As the carrier density changes from n=4 to n=1.5,the active orbitals for superconducting pairing shift from the quasi-two-dimensional orbital dxy to the quasi-one-dimensional orbitals d_(xz) and d_(yz).

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.