Frost deformation and microstructure evolution of porous rock under uniform and unidirectional freeze-thaw conditions

The frost deterioration and deformation of porous rock are commonly investigated under uniform freeze-thaw(FT)conditions.However,the unidirectional FT condition,which is also prevalent in engineering practice,has received limited attention.Therefore,a comparative study on frost deformation and microstructure evolution of porous rock under both uniform and unidirectional FT conditions was performed.Firstly,frost deformation experiments of rock were conducted under cyclic uniform and unidirectional FT action,respectively.Results illustrate that frost deformation of saturated rock exhibits isotropic characteristics under uniform FT cycles,while it shows anisotropic characteristics under unidirectional FT condition with both the frost heaving strain and residual strain along FT direction much higher than those perpendicular to FT direction.Moreover,the peak value and residual value of cumulative frost strain vary as logarithmic functions with cycle number under both uniform and unidirectional FT conditions.Subsequently,the microstructure evolution of rock suffered cyclic uniform and unidirectional FT action were measured.Under uniform FT cycles,newly generated pores uniformly distribute in rock and pore structure of rock remains isotropic in micro scale,and thus the frost deformation shows isotropic characteristics in macro scale.Under unidirectional FT cycles,micro-cracks or pore belts generate with their orientation nearly perpendicular to the FT direction,and rock structure gradually becomes anisotropic in micro scale,resulting in the anisotropic characteristics of frost deformation in macro scale.

Machine Learning Techniques in Predicting Hot Deformation Behavior of Metallic Materials

In engineering practice,it is often necessary to determine functional relationships between dependent and independent variables.These relationships can be highly nonlinear,and classical regression approaches cannot always provide sufficiently reliable solutions.Nevertheless,Machine Learning(ML)techniques,which offer advanced regression tools to address complicated engineering issues,have been developed and widely explored.This study investigates the selected ML techniques to evaluate their suitability for application in the hot deformation behavior of metallic materials.The ML-based regression methods of Artificial Neural Networks(ANNs),Support Vector Machine(SVM),Decision Tree Regression(DTR),and Gaussian Process Regression(GPR)are applied to mathematically describe hot flow stress curve datasets acquired experimentally for a medium-carbon steel.Although the GPR method has not been used for such a regression task before,the results showed that its performance is the most favorable and practically unrivaled;neither the ANN method nor the other studied ML techniques provide such precise results of the solved regression analysis.

Microstructure evolution of 7050 aluminum forgings during surface cumulative plastic deformation

To elucidate the mechanisms of regulating the microstructure uniformity in 7050 aluminum forgings through surface cumulative plastic deformation(SCPD),the microstructure under different solution treatments was investigated using metallographic observation(OM),electron backscatter diffraction(EBSD),transmission electron microscopy(TEM),and X-ray diffraction(XRD).The findings demonstrate that the most uniform microstructure in the forgings is achieved with a solution treatment at 470℃for 30 min.The SCPD process generates a significant number of needle-shaped precipitates,resulting in a higher dislocation density and stored energy.Solution treatments alleviate the pinning effect of second-phase particles and facilitate static recrystallization(SRX)in forgings,leading to a reduction in grain size.Additionally,mechanical testing results demonstrate 7%−13%increase in tensile strength and more uniform elongation of the forgings in different directions.

Plastic deformation mechanism of γ-phase U–Mo alloy studied by molecular dynamics simulations

Uranium–molybdenum(U–Mo) alloys are critical for nuclear power generation and propulsion because of their superior thermal conductivity, irradiation stability, and anti-swelling properties. This study explores the plastic deformation mechanisms of γ-phase U–Mo alloys using molecular dynamics(MD) simulations. In the slip model, the generalized stacking fault energy(GSFE) and the modified Peierls–Nabarro(P–N) model are used to determine the competitive relationships among different slip systems. In the twinning model, the generalized plane fault energy(GPFE) is assessed to evaluate the competition between slip and twinning. The findings reveal that among the three slip systems, the {110}<111>slip system is preferentially activated, while in the {112}<111> system, twinning is favored over slip, as confirmed by MD tensile simulations conducted in various directions. Additionally, the impact of Mo content on deformation behavior is emphasized. Insights are provided for optimizing process conditions to avoid γ → α′′ transitions, thereby maintaining a higher proportion of γ-phase U–Mo alloys for practical applications.

An extended discontinuous deformation analysis for simulation of grouting reinforcement in a water-rich fractured rock tunnel

Grouting has been the most effective approach to mitigate water inrush disasters in underground engineering due to its ability to plug groundwater and enhance rock strength.Nevertheless,there is a lack of potent numerical tools for assessing the grouting effectiveness in water-rich fractured strata.In this study,the hydro-mechanical coupled discontinuous deformation analysis(HM-DDA)is inaugurally extended to simulate the grouting process in a water-rich discrete fracture network(DFN),including the slurry migration,fracture dilation,water plugging in a seepage field,and joint reinforcement after coagulation.To validate the capabilities of the developed method,several numerical examples are conducted incorporating the Newtonian fluid and Bingham slurry.The simulation results closely align with the analytical solutions.Additionally,a set of compression tests is conducted on the fresh and grouted rock specimens to verify the reinforcement method and calibrate the rational properties of reinforced joints.An engineering-scale model based on a real water inrush case of the Yonglian tunnel in a water-rich fractured zone has been established.The model demonstrates the effectiveness of grouting reinforcement in mitigating water inrush disaster.The results indicate that increased grouting pressure greatly affects the regulation of water outflow from the tunnel face and the prevention of rock detachment face after excavation.

Deformation mechanism and NPR anchor cable truss coupling support in tunnel through fault fracture zone

To address the issue of extensive deformation in the Tabaiyi Tunnel caused by the fault zone,nuclear magnetic resonance(NMR)technology was employed to analyze the physical and mechanical properties of waterabsorbing mudstone.This analysis aimed to understand the mechanism behind the significant deformations.Drawing from the principle of excavation stress compensation,a support scheme featuring NPR anchorcables and an asymmetric truss support system was devised.To validate the scheme,numerical analysis using a combination of the Discrete Element Method(DEM)-Finite Element Method(FEM)was conducted.Additionally,similar material model tests and engineering measurements were carried out.Field experiments were also performed to evaluate the NPR anchor-cable and truss support system,focusing on anchor cable forces,pressures between the truss and surrounding rock,pressures between the initial support and secondary lining,as well as the magnitude of settlement and convergence deformation in the surrounding rock.The results indicate that the waterinduced expansion of clay minerals,resulting from damage caused by fissure water,accelerated the softening of the mudstone's internal structure,leading to significant deformations in the Tabaiyi Tunnel under high tectonic stress.The original support design fell short as the length of the anchor rods was smaller than the expansion depth of the plastic zone.As a result,the initial support structure bore the entire load from the surrounding rock,and a non-coupled deformation contact was observed between the double-arch truss and the surrounding rock.The adoption of NPR asymmetric anchor-cable support effectively restrained the expansion and asymmetric distribution characteristics of the plastic zone.Considering the mechanical degradation caused by water absorption in mudstone,the rigid constraint provided by the truss proved crucial for controlling the stability of the surrounding rock.These research findings hold significant implications for managing large deformations in soft rock tunnels situated within fractured zones under high tectonic stress conditions.

Deformation monitoring of long-span railway bridges based on SBAS-InSAR technology

The deformation monitoring of long-span railway bridges is significant to ensure the safety of human life and property.The interferometric synthetic aperture radar(In SAR)technology has the advantage of high accuracy in bridge deformation monitoring.This study monitored the deformation of the Ganjiang Super Bridge based on the small baseline subsets(SBAS)In SAR technology and Sentinel-1A data.We analyzed the deformation results combined with bridge structure,temperature,and riverbed sediment scouring.The results are as follows:(1)The Ganjiang Super Bridge area is stable overall,with deformation rates ranging from-15.6 mm/yr to 10.7 mm/yr(2)The settlement of the Ganjiang Super Bridge deck gradually increases from the bridge tower toward the main span,which conforms to the typical deformation pattern of a cable-stayed bridge.(3)The sediment scouring from the riverbed cause the serious settlement on the bridge’s east side compared with that on the west side.(4)The bridge deformation negatively correlates with temperature,with a faster settlement at a higher temperature and a slow rebound trend at a lower temperature.The study findings can provide scientific data support for the health monitoring of long-span railway bridges.

Application of GNSS-PPP on Dynamic Deformation Monitoring of Offshore Platforms

The real-time dynamic deformation monitoring of offshore platforms under environmental excitation is crucial to their safe operation.Although Global Navigation Satellite System-Precise Point Positioning(GNSS-PPP)has been considered for this purpose,its monitoring accuracy is relatively low.Moreover,the influence of background noise on the dynamic monitoring accuracy of GNSS-PPP remains unclear.Hence,it is imperative to further validate the feasibility of GNSS-PPP for deformation monitoring of offshore platforms.To address these concerns,vibration table tests with different amplitudes and frequencies are conducted.The results demonstrate that GNSS-PPP can effectively monitor horizontal vibration displacement as low as±30 mm,which is consistent with GNSS-RTK.Furthermore,the spectral characteristic of background noise in GNSS-PPP is similar to that of GNSS-RTK(Real Time Kinematic).Building on this observation,an improved Complete Ensemble Empirical Mode Decomposition with Adaptive Noise(CEEMDAN)has been proposed to de-noise the data and enhance the dynamic monitoring accuracy of GNSS-PPP.Field monitoring application research is also undertaken,successfully extracting and analyzing the dynamic deformation of an offshore platform structure under environmental excitation using GNSS-PPP monitoring in conjunction with improved CEEMDAN de-noising.By comparing the de-noised dynamic deformation trajectories of the offshore platform during different periods,it is observed that the platform exhibits reversible alternating vibration responses under environmental excitation,with more pronounced displacement deformation in the direction of load action.The research results confirm the feasibility and potential of GNSS-PPP for dynamic deformation monitoring of offshore platforms.

Violent collisions can reveal hexadecapole deformation of nuclei

A new observable in heavy ion collision experiments was identified to be sensitive to the hexadecapole deformation of the colliding nuclei.Such deformation is difficult to measure in traditional nuclear electric transition measurements,as it is often overwhelmed by the nuclear quadrupole deformation.This opens the door to gain new insight into nuclear structure with experiments that were designed to study hot and dense nuclear matter.

Spatiotemporal deformation characteristics of Outang landslide and identification of triggering factors using data mining

Since the impoundment of Three Gorges Reservoir(TGR)in 2003,numerous slopes have experienced noticeable movement or destabilization owing to reservoir level changes and seasonal rainfall.One case is the Outang landslide,a large-scale and active landslide,on the south bank of the Yangtze River.The latest monitoring data and site investigations available are analyzed to establish spatial and temporal landslide deformation characteristics.Data mining technology,including the two-step clustering and Apriori algorithm,is then used to identify the dominant triggers of landslide movement.In the data mining process,the two-step clustering method clusters the candidate triggers and displacement rate into several groups,and the Apriori algorithm generates correlation criteria for the cause-and-effect.The analysis considers multiple locations of the landslide and incorporates two types of time scales:longterm deformation on a monthly basis and short-term deformation on a daily basis.This analysis shows that the deformations of the Outang landslide are driven by both rainfall and reservoir water while its deformation varies spatiotemporally mainly due to the difference in local responses to hydrological factors.The data mining results reveal different dominant triggering factors depending on the monitoring frequency:the monthly and bi-monthly cumulative rainfall control the monthly deformation,and the 10-d cumulative rainfall and the 5-d cumulative drop of water level in the reservoir dominate the daily deformation of the landslide.It is concluded that the spatiotemporal deformation pattern and data mining rules associated with precipitation and reservoir water level have the potential to be broadly implemented for improving landslide prevention and control in the dam reservoirs and other landslideprone areas.

Effect of deformation parameters on the austenite dynamic recrystallization behavior of a eutectoid pearlite rail steel

Understandings of the effect of hot deformation parameters close to the practical production line on grain refinement are crucial for enhancing both the strength and toughness of future rail steels.In this work,the austenite dynamic recrystallization(DRX)behaviors of a eutectoid pearlite rail steel were studied using a thermo-mechanical simulator with hot deformation parameters frequently employed in rail production lines.The single-pass hot deformation results reveal that the prior austenite grain sizes(PAGSs)for samples with different deformation reductions decrease initially with an increase in deformation temperature.However,once the deformation temperature is beyond a certain threshold,the PAGSs start to increase.It can be attributed to the rise in DRX volume fraction and the increase of DRX grain with deformation temperature,respectively.Three-pass hot deformation results show that the accumulated strain generated in the first and second deformation passes can increase the extent of DRX.In the case of complete DRX,PAGS is predominantly determined by the deformation temperature of the final pass.It suggests a strategic approach during industrial production where part of the deformation reduction in low temperature range can be shifted to the medium temperature range to release rolling mill loads.

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.

Deformation Characteristics of Hydrate-Bearing Sediments

The safe and efficient development of natural gas hydrate requires a deep understanding of the deformation behaviors of reservoirs.In this study,a series of triaxial shearing tests are carried out to investigate the deformation properties of hydrate-bearing sediments.Variations of volumetric and lateral strains versus hydrate saturation are analyzed comprehensively.Results indicate that the sediments with high hydrate saturation show dilative behaviors,which lead to strain-softening characteristics during shearing.The volumetric strain curves have a tendency to transform gradually from dilatation to compression with the increase in effective confining pressure.An easy prediction model is proposed to describe the relationship between volumetric and axial strains.The model coefficientβis the key dominating factor for the shape of volumetric strain curves and can be determined by the hydrate saturation and stress state.Moreover,a modified model is established for the calculation of lateral strain.The corresponding determination method is provided for the easy estimation of model coefficients for medium sand sediments containing hydrate.This study provides a theoretical and experimental reference for deformation estimation in natural gas hydrate development.

Enhancing constitutive description and workability characterization of Mg alloy during hot deformation using machine learning-based Arrhenius-type model

Hot deformation is a commonly employed processing technique to enhance the ductility and workability of Mg alloy.However,the hot deformation of Mg alloy is highly sensitive to factors such as temperature,strain rate,and strain,leading to complex flow behavior and an exceptionally narrow processing window for Mg alloy.To overcome the shortcomings of the conventional Arrhenius-type(AT)model,this study developed machine learning-based Arrhenius-type(ML-AT)models by combining the genetic algorithm(GA),particle swarm optimization(PSO),and artificial neural network(ANN).Results indicated that when describing the flow behavior of the AQ80 alloy,the PSO-ANN-AT model demonstrates the most prominent prediction accuracy and generalization ability among all ML-AT and AT models.Moreover,an activation energy-processing(AEP)map was established using the reconstructed flow stress and activation energy fields based on the PSO-ANN-AT model.Experimental validations revealed that this AEP map exhibits superior predictive capability for microstructure evolution compared to the one established by the traditional interpolation methods,ultimately contributing to the precise determination of the optimum processing window.These findings provide fresh insights into the accurate constitutive description and workability characterization of Mg alloy during hot deformation.

Origin of nucleation and growth of extension twins in grains unsuitably oriented for twinning during deformation of Mg-1%Al

A large number of anomalous extension twins,with low or even negative twinning Schmid factors,were found to nucleate and grow in a strongly textured Mg-1Al alloy during tensile deformation along the extruded direction.The deformation mechanisms responsible for this behaviour were investigated through in-situ electron back-scattered diffraction,grain reference orientation deviation,and slip trace-modified lattice rotation.It was found that anomalous extension twins nucleated mainly at the onset of plastic deformation at or near grain boundary triple junctions.They were associated with the severe strain incompatibility between neighbour grains as a result from the differentbasal slip-induced lattice rotations.Moreover,the anomalous twins were able to grow with the applied strain due to the continuous activation ofbasal slip in different neighbour grains,which enhanced the strain incompatibility.These results reveal the complexity of the deformation mechanisms in Mg alloys at the local level when deformed along hard orientations.

Analytical model for predicting time-dependent lateral deformation of geosynthetics-reinforced soil walls with modular block facing

To date,few models are available in the literature to consider the creep behavior of geosynthetics when predicting the lateral deformation(d)of geosynthetics-reinforced soil(GRS)retaining walls.In this study,a general hyperbolic creep model was first introduced to describe the long-term deformation of geosynthetics,which is a function of elapsed time and two empirical parameters a and b.The conventional creep tests with three different tensile loads(Pr)were conducted on two uniaxial geogrids to determine their creep behavior,as well as the a-Pr and b-Pr relationships.The test results show that increasing Pr accelerates the development of creep deformation for both geogrids.Meanwhile,a and b respectively show exponential and negatively linear relationships with Pr,which were confirmed by abundant experimental data available in other studies.Based on the above creep model and relationships,an accurate and reliable analytical model was then proposed for predicting the time-dependent d of GRS walls with modular block facing,which was further validated using a relevant numerical investigation from the previous literature.Performance evaluation and comparison of the proposed model with six available prediction models were performed.Then a parametric study was carried out to evaluate the effects of wall height,vertical spacing of geogrids,unit weight and internal friction angle of backfills,and factor of safety against pullout on d at the end of construction and 5 years afterwards.The findings show that the creep effect not only promotes d but also raises the elevation of the maximum d along the wall height.Finally,the limitations and application prospects of the proposed model were discussed and analyzed.

Liquefaction susceptibility and deformation characteristics of saturated coral sandy soils subjected to cyclic loadings-a critical review

Coral sandy soils widely exist in coral island reefs and seashores in tropical and subtropical regions.Due to the unique marine depositional environment of coral sandy soils,the engineering characteristics and responses of these soils subjected to monotonic and cyclic loadings have been a subject of intense interest among the geotechnical and earthquake engineering communities.This paper critically reviews the progress of experimental investigations on the undrained behavior of coral sandy soils under monotonic and cyclic loadings over the last three decades.The focus of coverage includes the contractive-dilative behavior,the pattern of excess pore-water pressure(EPWP)generation and the liquefaction mechanism and liquefaction resistance,the small-strain shear modulus and strain-dependent shear modulus and damping,the cyclic softening feature,and the anisotropic characteristics of undrained responses of saturated coral sandy soils.In particular,the advances made in the past decades are reviewed from the following aspects:(1)the characterization of factors that impact the mechanism and patterns of EPWP build-up;(2)the identification of liquefaction triggering in terms of the apparent viscosity and the average flow coefficient;(3)the establishment of the invariable form of strain-based,stress-based,or energy-based EPWP ratio formulas and the unique relationship between the new proxy of liquefaction resistance and the number of cycles required to reach liquefaction;(4)the establishment of the invariable form of the predictive formulas of small strain modulus and strain-dependent shear modulus;and(5)the investigation on the effects of stress-induced anisotropy on liquefaction susceptibility and dynamic deformation characteristics.Insights gained through the critical review of these advances in the past decades offer a perspective for future research to further resolve the fundamental issues concerning the liquefaction mechanism and responses of coral sandy sites subjected to cyclic loadings associated with seismic events in marine environments.

Coseismic deformation and fault slip distribution of the 2023 M_(W)7.8 and M_(W)7.6 earthquakes in Türkiye

On February 6,2023,a devastating earthquake with a moment magnitude of M_(W)7.8 struck the town of Pazarcik in south-central Türkiye,followed by another powerful earthquake with a moment magnitude of M_(W)7.6 that struck the nearby city of Elbistan 9 h later.To study the characteristics of surface deformation caused by this event and the influence of fault rupture,this study calculated the static coseismic deformation of 56 stations and dynamic displacement waveforms of 15 stations using data from the Turkish national fixed global navigation satellite system(GNSS)network.A maximum static coseismic displacement of 0.38 m for the M_(W)7.8 Kahramanmaras earthquake was observed at station ANTE,36 km from the epicenter,and a maximum dynamic coseismic displacement of 4.4 m for the M_(W)7.6 Elbistan earthquake was observed at station EKZ1,5 km from the epicenter.The rupture-slip distributions of the two earthquakes were inverted using GNSS coseismic deformation as a constraint.The results showed that the Kahramanmaras earthquake rupture segment was distinct and exposed on the ground,resulting in significant rupture slip along the Amanos and Pazarcik fault segments of the East Anatolian Fault.The maximum slip in the Pazarcik fault segment was 10.7 m,and rupture occurred at depths of 0–15 km.In the Cardak fault region,the Elbistan earthquake caused significant ruptures at depths of 0–12 km,with the largest amount of slip reaching 11.6 m.The Coulomb stress change caused by the Kahramanmaras earthquake rupture along the Cardak fault segment was approximately 2 bars,and the area of increased Coulomb stress corresponded to the subsequent rupture region of the M_(W)7.6 earthquake.Thus,it is likely that the M_(W)7.8 earthquake triggered or promoted the M_(W)7.6 earthquake.Based on the cumulative stress impact of the M_(W)7.8 and M_(W)7.6 events,the southwestern segment of the East Anatolian Fault,specifically the Amanos fault segment,experienced a Coulomb rupture stress change exceeding 2 bars,warranting further attention to assess its future seismic hazard risk.

Application of modified discontinuous deformation analysis to dynamic modelling of the Baige landslide in the Jinsha River

Accurate dynamic modeling of landslides could help understand the movement mechanisms and guide disaster mitigation and prevention.Discontinuous deformation analysis(DDA)is an effective approach for investigating landslides.However,DDA fails to accurately capture the degradation in shear strength of rock joints commonly observed in high-speed landslides.In this study,DDA is modified by incorporating simplified joint shear strength degradation.Based on the modified DDA,the kinematics of the Baige landslide that occurred along the Jinsha River in China on 10 October 2018 are reproduced.The violent starting velocity of the landslide is considered explicitly.Three cases with different violent starting velocities are investigated to show their effect on the landslide movement process.Subsequently,the landslide movement process and the final accumulation characteristics are analyzed from multiple perspectives.The results show that the violent starting velocity affects the landslide motion characteristics,which is found to be about 4 m/s in the Baige landslide.The movement process of the Baige landslide involves four stages:initiation,high-speed sliding,impact-climbing,low-speed motion and accumulation.The accumulation states of sliding masses in different zones are different,which essentially corresponds to reality.The research results suggest that the modified DDA is applicable to similar high-level rock landslides.

Anisotropic strength,deformation and failure of gneiss granite under high stress and temperature coupled true triaxial compression

The anisotropic mechanical behavior of rocks under high-stress and high-temperature coupled conditions is crucial for analyzing the stability of surrounding rocks in deep underground engineering.This paper is devoted to studying the anisotropic strength,deformation and failure behavior of gneiss granite from the deep boreholes of a railway tunnel that suffers from high tectonic stress and ground temperature in the eastern tectonic knot in the Tibet Plateau.High-temperature true triaxial compression tests are performed on the samples using a self-developed testing device with five different loading directions and three temperature values that are representative of the geological conditions of the deep underground tunnels in the region.Effect of temperature and loading direction on the strength,elastic modulus,Poisson’s ratio,and failure mode are analyzed.The method for quantitative identification of anisotropic failure is also proposed.The anisotropic mechanical behaviors of the gneiss granite are very sensitive to the changes in loading direction and temperature under true triaxial compression,and the high temperature seems to weaken the inherent anisotropy and stress-induced deformation anisotropy.The strength and deformation show obvious thermal degradation at 200℃due to the weakening of friction between failure surfaces and the transition of the failure pattern in rock grains.In the range of 25℃ 200℃,the failure is mainly governed by the loading direction due to the inherent anisotropy.This study is helpful to the in-depth understanding of the thermal-mechanical behavior of anisotropic rocks in deep underground projects.