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

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

Support design method for deep soft-rock tunnels in non-hydrostatic high in-situ stress field

Due to the long-term plate tectonic movements in southwestern China,the in-situ stress field in deep formations is complex.When passing through deep soft-rock mass under non-hydrostatic high in-situ stress field,tunnels will suffer serious asymmetric deformation.There is no available support design method for tunnels under such a situation in existing studies to clarify the support time and support stiffness.This study first analyzed the mechanical behavior of tunnels in non-hydrostatic in-situ stress field and derived the theoretical equations of the ground squeezing curve(GSC)and ground loosening curve(GLC).Then,based on the convergence confinement theory,the support design method of deep soft-rock tunnels under non-hydrostatic high in-situ stress field was established considering both squeezing and loosening pressures.In addition,this method can provide the clear support time and support stiffness of the second layer of initial support.The proposed design method was applied to the Wanhe tunnel of the China-Laos railway in China.Monitoring data indicated that the optimal support scheme had a good effect on controlling the tunnel deformation in non-hydrostatic high in-situ stress field.Field applications showed that the secondary lining could be constructed properly.

Effect of Modulus Heterogeneity on the Equilibrium Shape and Stress Field ofαPrecipitate in Ti-6Al-4V

For media with inclusions(e.g.,precipitates,voids,reinforcements,and others),the difference in lattice parameter and the elastic modulus between the matrix and inclusions cause stress concentration at the interfaces.These stress fields depend on the inclusions’size,shape,and distribution and will respond instantly to the evolving microstructure.This study develops a phase-field model concerningmodulus heterogeneity.The effect of modulus heterogeneity on the growth process and equilibrium state of theαplate in Ti-6Al-4V during precipitation is evaluated.Theαprecipitate exhibits strong anisotropy in shape upon cooling due to the interplay of the elastic strain and interfacial energy.The calculated orientation of the habit plane using the homogeneous modulus ofαphase shows the smallest deviation fromthat of the habit plane observed in the experiment,compared to the case where the homogeneous modulus ofβphase is adopted.In addition,the equilibrium volume ofαphase within the systemusing homogeneousβmodulus exhibits the largest dependency on the applied stresses.The stress fields across theα/βinterface are further calculated under the assumption of modulus heterogeneity and compared to those using homogeneous modulus of eitherαorβphase.This study provides an essential theoretical basis for developing mechanics models concerning systems with heterogeneous structures.

The 2023 Turkey earthquake doublet: Earthquake relocation, seismic tomography, and stress field inversion

On February 6,2023,two earthquakes with magnitudes of M_(W) 7.8 and M_(W) 7.5 struck southeastern Turkey,causing significant casualties and economic losses.These seismic events occurred along the East Anatolian Fault Zone,a convergent boundary between the Arabian Plate and the Anatolian Subplate.In this study,we analyze the M_(W) 7.8 and M_(W) 7.5 earthquakes by comparing their aftershock relocations,tomographic images,and stress field inversions.The earthquakes were localized in the upper crust and exhibited steep dip angles.Furthermore,the aftershocks occurred either close to the boundaries of low and high P-wave velocity anomaly zones or within the low P-wave velocity anomaly zones.The East Anatolia Fault,associated with the M_(W) 7.8 earthquake,and the SürgüFault,related to the M_(W) 7.5 earthquake,predominantly experienced shear stress.However,their western sections experienced a combination of strike-slip and tensile stresses in addition to shear stress.The ruptures of the M_(W) 7.8 and M_(W) 7.5 earthquakes appear to have bridged a seismic gap that had seen sparse seismicity over the past 200 years prior to the 2023 Turkey earthquake sequence.

Stress corrosion cracking behavior of buried oil and gas pipeline steel under the coexistence of magnetic field and sulfate-reducing bacteria

Magnetic field and microorganisms are important factors influencing the stress corrosion cracking(SCC)of buried oil and gas pipelines. Once SCC occurs in buried pipelines, it will cause serious hazards to the soil environment. The SCC behavior of X80 pipeline steel under the magnetic field and sulfate-reducing bacteria(SRB) environment was investigated by immersion tests, electrochemical tests, and slow strain rate tensile(SSRT) tests. The results showed that the corrosion and SCC sensitivity of X80 steel decreased with increasing the magnetic field strength in the sterile environment. The SCC sensitivity was higher in the biotic environment inoculated with SRB, but it also decreased with increasing magnetic field strength, which was due to the magnetic field reduces microbial activity and promotes the formation of dense film layer. This work provided theoretical guidance on the prevention of SCC in pipeline steel under magnetic field and SRB coexistence.

A novel indirect optical method for rock stress measurement using microdeformation field analysis

Stress measurement plays a crucial role in geomechanics and rock engineering,especially for the design and construction of large-scale rock projects.This paper presents a novel method,based on the traditional stress relief approach,for indirectly measuring rock stress using optical techniques.The proposed method allows for the acquisition of full-field strain evolution on the borehole’s inner wall before and after disturbance,facilitating the determination of three-dimensional(3D)stress information at multiple points within a single borehole.The study focuses on presenting the method’s theoretical framework,laboratory validation results,and equipment design conception.The theoretical framework comprises three key components:the optical imaging method of the borehole wall,the digital image correlation(DIC)method,and the stress calculation procedure.Laboratory validation tests investigate strain field distribution on the borehole wall under varying stress conditions,with stress results derived from DIC strain data.Remarkably,the optical method demonstrates better measurement accuracy during the unloading stage compared to conventional strain gauge methods.At relatively high stress levels,the optical method demonstrates a relative error of less than 7%and an absolute error within 0.5 MPa.Furthermore,a comparative analysis between the optical method and the conventional contact resistance strain gauge method highlights the optical method’s enhanced accuracy and stability,particularly during the unloading stage.The proposed optical stress measurement device represents a pioneering effort in the application of DIC technology to rock engineering,highlighting its potential to advance stress measurement techniques in the field.

Impact of well placement and flow rate on production efficiency and stress field in the fractured geothermal reservoirs

Geothermal energy has gained wide attention as a renewable alternative for mitigating greenhouse gas emissions.The advancements in enhanced geothermal system technology have enabled the exploitation of previously inaccessible geothermal resources.However,the extraction of geothermal energy from deep reservoirs poses many challenges due to high‐temperature and high‐geostress conditions.These factors can significantly impact the surrounding rock and its fracture formation.A comprehensive understanding of the thermal–hydraulic–mechanical(THM)coupling effect is crucial to the safe and efficient exploitation of geothermal resources.This study presented a THM coupling numerical model for the geothermal reservoir of the Yangbajing geothermal system.This proposed model investigated the geothermal exploitation performance and the stress distribution within the reservoir under various combinations of geothermal wells and mass flow rates.The geothermal system performance was evaluated by the criteria of outlet temperature and geothermal productivity.The results indicate that the longer distance between wells can increase the outlet temperature of production wells and improve extraction efficiency in the short term.In contrast,the shorter distance between wells can reduce the heat exchange area and thus mitigate the impact on the reservoir stress.A larger mass flow rate is conducive to the production capacity enhancement of the geothermal system and,in turn causes a wider range of stress disturbance.These findings provide valuable insights into the optimization of geothermal energy extraction while considering reservoir safety and long‐term sustainability.This study deepens the understanding of the THM coupling effects in geothermal systems and provides an efficient and environmentally friendly strategy for a geothermal energy system.

Simulation Analysis of Stress Field of Walnut Shell Composite Powder in Laser Additive Manufacturing Forming

A calculation model of stress field in laser additive manufacturing of walnut shell composite powder(walnut shell/Co-PES powder)was established.The DFLUX subroutine was used to implement the moveable application of a double ellipsoid heat source by considering the mechanical properties varying with temperature.The stress field was simulated by the sequential coupling method,and the experimental results were in good accordance with the simulation results.In addition,the distribution and variation of stress and strain field were obtained in the process of laser additive manufacturing of walnut shell composite powder.The displacement of laser additive manufacturing walnut shell composite parts gradually decreased with increasing preheating temperature,decreasing laser power and increasing scanning speed.During the cooling process,the displacement of laser additive manufacturing of walnut shell composite parts gradually increased with the increasing preheating temperature,decreasing scanning speed and increasing laser power.

Effect of hydraulic fracturing induced stress field on weak surface activation during unconventional reservoir development

Unconventional reservoirs usually contain many weak surfaces such as faults,laminae and natural fractures,and effective activation and utilization of these weak surfaces in reservoirs can significantly improve the extraction effect.In hydraulic fracturing,when the artificial fracture approaches the natural fracture,the natural fracture would be influenced by both the original in-situ stress field and the hydraulic fracturing-induced stress field.In this paper,the hydraulic fracturing-induced stress field is calculated based on the relative position of hydraulic fracture and natural fracture,the original in-situ stress,the net pressure inside the hydraulic fracture and the pore pressure of the formation.Furthermore,the stability model of the natural fracture is established by combining the Mohr-Coulomb rupture criterion,and extensive parametric studies are conducted to explore the impact of each parameter on the stability of the natural fracture.The validity of the proposed model is verified by comparing with the reservoir characteristics and fracturing process of the X-well 150e155 formation in the Songliao Basin.It is found that the stress field induced by the hydraulic fracture inhibits the activation of the natural fracture after the artificial fracture crossed the natural fracture.Therefore,for similar reservoirs as X-well 150e155,it is suggested to connect natural fractures with hydraulic fractures first and then activate natural fractures which can effectively utilize the natural fractures and form a complex fracture network.

Tensile Fractures and in situ Stress Measurement Data Constraints on Cretaceous-Present Tectonic Stress Field Evolution of the Tanlu Fault Zone in Shandong Province,North China Craton

Tectonic stress fields are the key drivers of tectonic events and the evolution of regional structures.The tectonic stress field evolution of the Tanlu fault zone in Shandong Province,located in the east of the North China Craton(NCC),may have preserved records of the NCC’s tectonic history.Borehole television survey and hydraulic fracturing were conducted to analyze the paleo and present tectonic stress fields.Three groups of tensile fractures were identified via borehole television,their azimuths being NNW-SSE,NW-SE and NE-SW,representing multiple stages of tectonic events.Hydraulic fracturing data indicates that the study region is experiencing NEE-SWW-oriented compression and nearly-N-Soriented extension,in accordance with strike-slip and compression.Since the Cretaceous,the orientation of the extensional stress has evolved counterclockwise and sequentially from nearly-NW-SE-oriented to NE-SW-oriented and even nearly N-S-oriented,the stress state having transitioned from strike-slip-extension to strike-slip-compression,in association with the rotating and oblique subduction of the Pacific Plate beneath the NCC,with the participation of the Indian Plate.

Development of a time-dependent energy model to calculate the mining-induced stress over gates and pillars

Generally, longwall mining-induced stress results from the stress relaxation due to destressed zone that occurs above the mined panel. Knowledge of induced stress is very important for accurate design of adjacent gateroads and intervening pillars which helps to raise the safety and productivity of longwall mining operations. This study presents a novel time-dependent analytical model for determination of the longwall mining-induced stress and investigates the coefficient of stress concentration over adjacent gates and pillars. The model is developed based on the strain energy balance in longwall mining incorporated to a rheological constitutive model of caved materials with time-varying parameters. The study site is the Tabas coal mine of Iran. In the proposed model, height of destressed zone above the mined panel, total longwall mining-induced stress, abutment angle, induced vertical stress, and coefficient of stress concentration over neighboring gates and intervening pillars are calculated. To evaluate the effect of proposed model parameters on the coefficient of stress concentration due to longwall mining, sensitivity analysis is performed based on the field data and experimental constants. Also, the results of the proposed model are compared with those of existing models. The comparative results confirm a good agreement between the proposed model and the in situ measurements. According to the obtained results, it is concluded that the proposed model can be successfully used to calculate the longwall mining-induced stress. Therefore, the optimum design of gate supports and pillar dimensions would be attainable which helps to increase the mining efficiency.

Evaluation of stress voltage on off-state time-dependent breakdown for GaN MIS-HEMT with SiNx gate dielectric

Stress voltages on time-dependent breakdown characteristics of GaN MIS-HEMTs during negative gate bias stress (with VGS < 0, VD = VS = 0) and off-state stress (VG < VTh, VDS > 0, VS = 0) are investigated. For negative bias stress, the breakdown time distribution (β) decreases with the increasing negative gate voltage, while β is larger for higher drain voltage at off-state stress. Two humps in the time-dependent gate leakage occurred under both breakdown conditions, which can be ascribed to the dielectric breakdown triggered earlier and followed by the GaN layer breakdown. Combining the electric distribution from simulation and long-term monitoring of electric parameter, the peak electric fields under the gate edges at source and drain sides are confirmed as the main formation locations for per-location paths during negative gate voltage stress and off-state stress, respectively.

3D digital-image correlation insight into generalized relaxation behavior of sandstone under stress and pore pressure coupling

The occurrence of geological hazards and the instability of geotechnical engineering structures are closely related to the time-dependent behavior of rock.However,the idealization boundary condition for constant stress in creep or constant strain in relaxation is not usually attained in natural geological systems.Therefore,generalized relaxation tests that explore the simultaneous changes of stress and strain with time under different stress levels with constant pore-water pressure are conducted in this study.The results show that in area Ⅰ,area Ⅱ,and area Ⅲ,the stress and strain both change synchronously with time and show similar evolutionary laws as the strain-time curve for creep or the stress-time curve for relaxation.When the applied stress level surpasses the δ_(ci) or δ_(cd) threshold,the variations in stress and strain and their respective rates of change exhibit a significant increase.The radial deformation and its rate of change exhibit greater sensitivity in response to stress levels.The apparent strain deforms homogeneously at the primary stage,and subsequently,gradually localizes due to the microcrack development at the secondary stage.Ultimately,interconnection of the microcracks causes the formation of a shear-localization zone at the tertiary stage.The strain-time responses inside and outside the localization zone are characterized by local strain accumulation and inelastic unloading during the secondary and tertiary stages,respectively.The width of the shear-localization zone is found to range from 4.43 mm to 7.08 mm and increased with a longer time-to-failure.Scanning electron microscopy(SEM)reveals a dominant coalescence of intergranular cracks on the fracture surface,and the degree of physiochemical deterioration caused by water-rock interaction is more severe under a longer lifetime.The brittle sandstone’s time-dependent deformation is essentially controlled by microcrack development during generalized relaxation,and its expectancy-life is determined by its initial microstructural state and the rheological path.

Mechanism of principal stress rotation and deformation failure behavior induced by excavation in roadways

The failure modes of rock after roadway excavation are diverse and complex.A comprehensive investigation of the internal stress field and the rotation behavior of the stress axis in roadways is essential for elucidating the mechanism of roadway failure.This study aimed to examine the spatial relationship between roadways and stress fields.The law of stress axis rotation under three-dimensional(3D)stress has been extensively studied.A stress model of roadways in the spatial stress field was established,and the far-field stress state at different spatial positions of the roadways was analyzed.A mechanical model of roadways under a 3D stress state was established using far-field stress solutions as boundary conditions.The distribution of principal stressesσ1,σ2 andσ3 around the roadways and the variation of the stress principal axis were solved.It was found that the stability boundary of the stress principal axis exhibits hysteresis when compared with that of the principal stress magnitudes.A numerical analysis model for spatial roadways was established to validate the distribution of principal stress and the mechanism of principal axis rotation.Research has demonstrated that the stress axis undergoes varying degrees of spatial rotation in different orientations and radial depths.Based on the distribution of principal stress and the rotation law of the stress principal axis,the entire evolution mechanism of the two stress adjustments to form the final failure form after roadway excavation has been revealed.The on-site detection results also corroborate the findings presented in this paper.The results provide a basis for the analysis of the failure mechanism under a 3D stress state.

Calculation of Mass Concrete Temperature and Creep Stress under the Influence of Local Air Heat Transfer

Temperature-induced cracking during the construction of mass concrete is a significant concern.Numerical simulations of concrete temperature have primarily assumed that the concrete is placed in an open environment.The problem of heat transfer between the air and concrete has been simplified to the concrete’s heat dissipation boundary.However,in the case of tubular concrete structures,where air inlet and outlet are relatively limited,the internal air temperature does not dissipate promptly to the external environment as it rises.To accurately simulate the temperature and creep stress in tubular concrete structures with enclosed air spaces during construction,we establish an air–concrete coupled heat transfer model according to the principles of conjugate heat transfer,and the accuracy of the model is verified through experiments.Furthermore,we conduct a case study to analyze the impact of airflow within the ship lock corridor on concrete temperature and creep stress.The results demonstrate that enhancing airflow within the corridor can significantly reduce the maximum concrete temperature.Compared with cases in which airflow within the corridor is neglected,the maximum concrete temperature and maximum tensile stress can be reduced by 12.5℃ and 0.7 MPa,respectively,under a wind speed of 4 m/s.The results of the traditional calculation method are relatively close to those obtained at a wind speed of 1 m/s.However,the temperature reduction process in the traditional method is faster,and the method yields greater tensile stress values for the corridor location.

Acetaminophen overdose-induced acute liver injury can be alleviated by static magnetic field

Acetaminophen(APAP),the most frequently used mild analgesic and antipyretic drug worldwide,is implicated in causing 46%of all acute liver failures in the USA and between 40%and 70%in Europe.The predominant pharmacological intervention approved for mitigating such overdose is the antioxidant N-acetylcysteine(NAC);however,its efficacy is limited in cases of advanced liver injury or when administered at a late stage.In the current study,we discovered that treatment with a moderate intensity static magnetic field(SMF)notably reduced the mortality rate in mice subjected to high-dose APAP from 40%to 0%,proving effective at both the initial liver injury stage and the subsequent recovery stage.During the early phase of liver injury,SMF markedly reduced APAPinduced oxidative stress,free radicals,and liver damage,resulting in a reduction in multiple oxidative stress markers and an increase in the antioxidant glutathione(GSH).During the later stage of liver recovery,application of vertically downward SMF increased DNA synthesis and hepatocyte proliferation.Moreover,the combination of NAC and SMF significantly mitigated liver damage induced by high-dose APAP and increased liver recovery,even 24 h post overdose,when the effectiveness of NAC alone substantially declines.Overall,this study provides a noninvasive non-pharmaceutical tool that offers dual benefits in the injury and repair stages following APAP overdose.Of note,this tool can work as an alternative to or in combination with NAC to prevent or minimize liver damage induced by APAP,and potentially other toxic overdoses.

Experimental study on time-dependent stress and strain of in-plane shear(ModeⅡ) fracture process of rock

Shear-box test with strain measurement was used to study time-dependent stress and strain of in-plane shear(Mode Ⅱ) fracture process of rock and to reveal the mechanism of Mode Ⅱ fracture.Numerical results show that the maximum shear stress τmax at the crack tip is much larger than the maximum tensile stress σ1 and the ratio of τmax/σ1 is about 5,which favors Mode Ⅱ fracture of rock.Test results indicate that the strain-time curve comprises three stages:the linear deformation stage,the micro-cracking stage and the macroscopic crack propagation.The strain in the direction of the original notch plane is negative,due to restraining effect of compressive loading applied to the original notch plane.Both σ1 and τmax are increased as the load increases,but the slope of τmax is larger than that of σ1 and the value of τmax is always larger than that of σ1.Therefore,τmax reaches its limited value at peak load before σ1 and results in Mode Ⅱ fracture of rock.Shear-box(i.e.compression-shear) test becomes a potential standard method for achieving the true Mode Ⅱ fracture and determining Mode Ⅱ fracture toughness of rock.

Photodetachment dynamics of H^- ion in a harmonic potential plus a time-dependent oscillating electric field

The photodetachment dynamics of H^- ion in a harmonic potential plus an oscillating electric field is studied using the time-dependent closed orbit theory. An analytical formula for calculating the photodetachment cross section of this system is put forward. It is found that the photodetachment cross section of this system is nearly unaffected for the weak oscillating electric field strength, but oscillates complicatedly when the oscillating electric field strength turns strong. In addition, the frequency of the harmonic potential and the oscillating electric field （the frequency of the harmonic potential and the frequency of the oscillating electric field are the same in the paper, unless otherwise stated.） can also affect the photodetachment dynamics of this system. With the increase of the frequency in the harmonic potential and the oscillating electric field, the number of the closed orbits for the detached electrons increased, which makes the oscillatory structure in the photodetachment cross section much more complex. Our study presents an intuitive understanding of the photodetachment dynamics driven by a harmonic potential plus an oscillating electric field from a space and time dependent viewpoint. This study is very useful in guiding the future experimental research for the photodetachment dynamics in the electric field both changing with space and time.

Efficient solver for time-dependent Schrodinger equation with interaction between atoms and strong laser field

We present a parallel numerical method of simulating the interaction of atoms with a strong laser field by solving the time-depending Schr?dinger equation(TDSE) in spherical coordinates. This method is realized by combining constructing block diagonal matrices through using the real space product formula(RSPF) with splitting out diagonal sub-matrices for short iterative Lanczos(SIL) propagator. The numerical implementation of the solver guarantees efficient parallel computing for the simulation of real physical problems such as high harmonic generation(HHG) in these interaction systems.

Residual elastic stress strain field and geometrically necessary dislocation density distribution around nano-indentation in TA15 titanium alloy

Nanoindentation and high resolution electron backscatter diffraction（EBSD） were combined to examine the elastic modulus and hardness of α and β phases,anisotropy in residual elastic stress strain fields and distributions of geometrically necessary dislocation（GND） density around the indentations within TA15 titanium alloy.The nano-indention tests were conducted on α and β phases,respectively.The residual stress strain fields surrounding the indentation were calculated through crosscorrelation method from recorded patterns.The GND density distribution around the indentation was calculated based on the strain gradient theories to reveal the micro-mechanism of plastic deformation.The results indicate that the elastic modulus and hardness for α p hase are 129.05 GPas and 6.44 GPa,while for β phase,their values are 109.80 GPa and 4.29 GPa,respectively.The residual Mises stress distribution around the indentation is relatively heterogeneous and significantly influenced by neighboring soft β phase.The region with low residual stress around the indentation is accompanied with markedly high a type and prismatic-GND density.