In recent years,there has been an increasing need for climate information across diverse sectors of society.This demand has arisen from the necessity to adapt to and mitigate the impacts of climate variability and cha...In recent years,there has been an increasing need for climate information across diverse sectors of society.This demand has arisen from the necessity to adapt to and mitigate the impacts of climate variability and change.Likewise,this period has seen a significant increase in our understanding of the physical processes and mechanisms that drive precipitation and its variability across different regions of Africa.By leveraging a large volume of climate model outputs,numerous studies have investigated the model representation of African precipitation as well as underlying physical processes.These studies have assessed whether the physical processes are well depicted and whether the models are fit for informing mitigation and adaptation strategies.This paper provides a review of the progress in precipitation simulation overAfrica in state-of-the-science climate models and discusses the major issues and challenges that remain.展开更多
The dissolution of MgO-refractory into the slag had an obvious influence on the steel-slag reaction and the slag property,especially for high-aluminum steels.The dissolution behavior of MgO-refractory was investigated...The dissolution of MgO-refractory into the slag had an obvious influence on the steel-slag reaction and the slag property,especially for high-aluminum steels.The dissolution behavior of MgO-refractory was investigated under various conditions,including the temperature,the initial steel composition,and the initial slag composition.A steel-slag-refractory kinetic model for high-aluminum steel was developed,which incorporated the process of MgO-refractory dissolution.The dependence of the MgO mass transfer coefficient k_(MgO)^(r)on temperature T during MgO-refractory dissolution process was established,as described by ln k_(MgO)^(r)=63,754/T+24.38524.It was indicated that the MgO dissolution rate was significantly influenced by the temperature.A higher temperature increased the dissolution rate of MgO.The initial steel composition had a slight impact on the MgO dissolution rate.Additionally,the initial slag composition strongly impacted the MgO saturation concentration and the dissolution rate.A lower initial Al_(2)O_(3)/SiO_(2)ratio increased the MgO dissolution rate.The steel-slag-refractory kinetic model accurately predicted the dissolution of MgO-refractory and the influence of dissolved MgO on the viscosity and composition change during steel-slag-refractory reactions.It was suggested that a higher temperature can hardly reduce the viscosity due to the dissolution of the MgO-refractory.展开更多
In materials science and engineering design,high-fidelity and high-efficiency numerical simulation has become a driving force for innovation and practical implementation.To address longstanding bottlenecks in the deve...In materials science and engineering design,high-fidelity and high-efficiency numerical simulation has become a driving force for innovation and practical implementation.To address longstanding bottlenecks in the development of conventional material constitutive models—such as lengthy modeling cycles and difficulties in numerical implementation—this study proposes an intelligent modeling and code generation approach powered by large languagemodels.A structured knowledge base integrating constitutive theory,numerical algorithms,and UMAT(User Material)interface specifications is constructed,and a retrieval-augmented generation strategy is employed to establish an end-to-end workflow spanning experimental data parsing,constitutive model formulation,and automatic UMAT subroutine generation.Experimental results show that the method achieves high accuracy for both a classical Johnson–Cookmodel and a physics-informed neural network(PINN)model,with key parameter identification errors below 5%.Moreover,the automatically generated UMAT subroutines yield finite element simulation results in Abaqus that are highly consistent with theoretical predictions(coefficient of determination R2>0.98)while maintaining good numerical stability.This framework is currently focused on the automatic construction of rate-dependent elastoplastic material models,and its core method also provides a clear path for extending to other constitutive categories such as hyperelasticity and viscoelasticity.This work provides an effective technical route for the rapid development and reliable numerical implementation of material constitutive models,significantly advancing the intelligence level of computational mechanics research and improving engineering application efficiency.展开更多
(Quasi-)closed-form results for the statistical properties of unmanned aerial vehicle(UAV)airto-ground channels are derived for the first time using a novel spatial-vector-based method from a threedimensional(3-D)arbi...(Quasi-)closed-form results for the statistical properties of unmanned aerial vehicle(UAV)airto-ground channels are derived for the first time using a novel spatial-vector-based method from a threedimensional(3-D)arbitrary-elevation one-cylinder model.The derived results include a closed-form expression for the space-time correlation function and some quasi-closed-form ones for the space-Doppler power spectrum density,the level crossing rate,and the average fading duration,which are shown to be the generalizations of those previously obtained from the two-dimensional(2-D)one-ring model and the 3-D low-elevation one-cylinder model for terrestrial mobile-to-mobile channels.The close agreements between the theoretical results and the simulations as well as the measurements validate the utility of the derived channel statistics.Based on the derived expressions,the impacts of some parameters on the channel characteristics are investigated in an effective,efficient,and explicable way,which leads to a general guideline on the manual parameter estimation from the measurement description.展开更多
The India-Asia collision resulted in the formation of Qinghai-Tibet Plateau.Lower crustal flow model was proposed to explain the mechanism of Cenozoic tectonic deformation of Qinghai-Tibet Plateau.In this study,we pro...The India-Asia collision resulted in the formation of Qinghai-Tibet Plateau.Lower crustal flow model was proposed to explain the mechanism of Cenozoic tectonic deformation of Qinghai-Tibet Plateau.In this study,we propose a new approach by combining centrifugal analog modeling with numerical simulation to simulate the tectonic uplift history of the plateau based on the lower crustal flow model,and to investigate the material migration characteristics and the influence of crustal motion velocity and ductile layer viscosity on the plateau tectonic geomorphology.The models reproduce steep-sided flat-topped geomorphic features and clockwise rotation of the material at eastern Himalayan Syntaxis,verifying the rationality of the models.The results show that the greater the crustal motion velocity and the greater the ductile layer viscosity,the steeper the terrain change;and conversely,the smaller the crustal motion velocity and the smaller the ductile layer viscosity,the gentler the terrain change.This study further indicates that the weak lower crust plays an important role in the formation of geomorphic features and material migration characteristics of Qinghai-Tibet Plateau,and provides a new insight for the study of the uplift mechanism of the Tibetan Plateau.展开更多
In rock engineering,natural cracks in rock masses subjected to external loads tend to initiate and propagate,leading to potential safety hazards.To investigate the effect of cracking behavior on the mechanical propert...In rock engineering,natural cracks in rock masses subjected to external loads tend to initiate and propagate,leading to potential safety hazards.To investigate the effect of cracking behavior on the mechanical properties of rocks,the cracking processes of pre-cracked rocks have been extensively studied using numerical modeling methods.The peridynamics(PD)exhibits advantages over other numerical methods due to the absence of the requirements for remeshing and external crack growth criterion.However,for modeling pre-cracked rock cracking processes under impact,current PD implementations lack generally applicable rock constitutive models and impact contact models,which leads to difficulties in determining rock material parameters and efficiently calculating impact loads.This paper proposes a non-ordinary state-based peridynamics(NOSBPD)modeling method integrating the Drucker-Prager(DP)plasticity model and an efficient contact model to address the above problems.In the proposed method,the Drucker-Prager plasticity model is integrated into the NOSBPD,thereby equipping NOSBPD with the capability to accurately characterize the nonlinear stress-strain relationship inherent in rocks.An efficient contact model between particles and meshes is designed to calculate the impact loads,which is essentially a coupling method of PD with the finite element method(FEM).The effectiveness of the proposed NOSBPD modeling method is verified by comparison with other numerical methods and experiments.Experimental results indicate that the proposed method can effectively and accurately predict the 3D cracking processes of pre-cracked cracks under impact loading,and the maximum principal stress is the key driver behind wing crack formation in pre-cracked rocks.展开更多
In-space cable-driven manipulators exhibit several advantages,such as a large range of motion,high dexterity,and lightweight structure.However,kinematic and dynamic analysis play an essential role in designing a cable...In-space cable-driven manipulators exhibit several advantages,such as a large range of motion,high dexterity,and lightweight structure.However,kinematic and dynamic analysis play an essential role in designing a cable-driven manipulator.In this paper,the kinematic analysis of a type of cable-driven manipulator is performed,and a motion planning scheme is conducted to actuate this manipulator.Moreover,a flexible multi-body dynamic model of a cable-driven manipulator considering the frictional contact between the cables and pulleys is established.To describe properties such as flexibility,vibration,and variable length of the cable,this paper utilizes reducedorder beam elements of the Absolute Nodal Coordinates Formulation(ANCF)in Arbitrary Lagrangian Eulerian(ALE)framework.Additionally,a virtual element is introduced to model the contact segment in the cable-pulley system.A tension decay factor is employed to account for the friction in the contact segment.To validate the proposed method,a semi-analytical model based on D'Alembert's principle is established.Cross-verification is performed to validate the accuracy of both models.The model is further applied to simulate the rotation of the cable-driven manipulator with different structural parameters and frictional factors.The results from the analyses provide valuable guidance for the design and motion control of the in-space cable-driven manipulator.Finally,a prototype of a single module is manufactured and tested.Ground experiments are carried out to verify the kinematic and dynamic models.展开更多
Predicting rock blasting outcomes in mining has been crucial since its inception.Blasting remains the most energy-and cost-efficient method for rock breaking and is often the only practical option.However,the mechanis...Predicting rock blasting outcomes in mining has been crucial since its inception.Blasting remains the most energy-and cost-efficient method for rock breaking and is often the only practical option.However,the mechanism is complex,influenced by various rock properties,explosives,and blast design parameters,making their effects difficult to quantify.Traditional stress-based models struggle with many parameters,such as stress and Poisson's ratio,which are challenging to measure in the field.Empirical models,though simpler,often oversimplify blast conditions.Both types of models are limited to simulating a few blastholes and cannot handle full-scale blasts involving hundreds of blastholes.However,modeling full-scale blasts with all blast design parameters is most required for modern mining applications.This paper presents a novel strain-based modeling approach for blasting and geomechanical applications,utilizing measurable variables such as particle velocity,strain,and displacement.By bypassing complex constitutive relations,strain-based models capture critical blasting trends and simulate full-scale blasts with full-blast design parameters with minimal calibration.The framework encompasses field strain measurements,model construction based on measurable variables,and laboratoryderived strain-failure criteria,each offering potential for future enhancement.Additionally,a standardized field test for site characterization is recommended.The approach is demonstrated through the Multiple Blasthole Fragmentation model,which simulates rock fragmentation and fragment strain during blasting,highlighting the practicality and effectiveness of strain-based modeling for multiple blasthole blasts.Moreover,this approach extends beyond blasting,with potential applications in highwall stability monitoring and other geomechanical applications.Strain-based modeling provides a simplified yet effective solution,avoiding the complexities of rock constitutive relations and field stress measurements while enabling full-blast design simulations for large-scale field blasts.展开更多
Glassy polymers are widely used in biomedical applications in a solvent environment,yet their long-term performance is governed by the competing effects of physical aging and solvent-induced plasticization.Here,we dev...Glassy polymers are widely used in biomedical applications in a solvent environment,yet their long-term performance is governed by the competing effects of physical aging and solvent-induced plasticization.Here,we develop a constitutive model that explicitly couples the solvent concentration,structural relaxation,and mechanical response.This framework is built on a multiplicative decomposition of deformation and an Eyring-type flow rule,with structural evolution described by an effective temperature.A generalized shift factor is introduced to quantify how the solvent concentration and effective temperature jointly affect the relaxation time,thereby integrating physical aging and plasticization.The model is subsequently applied to methacrylate(MA)-based copolymer networks immersed in phosphate-buffered saline for up to nine months.Simulations accurately capture key experimental features,including the strong softening of highly swellable networks,the partial recovery due to aging,and the mitigating role of hydrophobic crosslinking in reducing solvent uptake.While the current single-mode description cannot reproduce the full relaxation spectrum,it establishes an efficient framework for predicting the long-term mechanical performance under coupled environmental and mechanical loading.This study provides a constitutive description of solvent-swollen glassy polymers,offering mechanistic insight into the interplay between plasticization and aging.Beyond biomedical MA networks,this framework establishes a foundation for predicting the long-term performance of polymer glasses under coupled aqueous environmental and mechanical loading.展开更多
Huperzine A(HupA) is a highly selective, reversible acetylcholinesterase(AChE) inhibitor that exhibits neuroprotective effects and is clinically used to manage benign memory decline.However, the specific relationship ...Huperzine A(HupA) is a highly selective, reversible acetylcholinesterase(AChE) inhibitor that exhibits neuroprotective effects and is clinically used to manage benign memory decline.However, the specific relationship between the pharmacokinetic(PK) profile of HupA and cerebral acetylcholine(ACh) dynamics remains poorly characterized. Here, we characterize the PK-pharmacodynamic(PD) properties of HupA in rats under both physiological and pathological conditions. Following a single intramuscular injection, HupA exhibits a short halflife but rapid brain penetration, while multiple dosing significantly enhances its brain exposure. In a middle cerebral artery occlusion(MCAO) rat model, HupA demonstrates increased brain distribution. Furthermore, HupA elevates ACh concentrations across multiple brain regions, concurrently modulating several monoamine neurotransmitters. Using a minimal physiologically based pharmacokinetic-pharmacodynamic(mPBPK-PD) modeling approach,cerebral ACh dynamics were accurately predicted based on the pharmacokinetics of HupA in systemic circulation. The developed mPBPK-PD model exhibits robust predictive performance and holds potential for guiding the optimization of clinical dosing regimens and improving the therapeutic efficacy of HupA.展开更多
Activation pruning reduces neural network complexity by eliminating low-importance neuron activations,yet identifying the critical pruning threshold—beyond which accuracy rapidly deteriorates—remains computationally...Activation pruning reduces neural network complexity by eliminating low-importance neuron activations,yet identifying the critical pruning threshold—beyond which accuracy rapidly deteriorates—remains computationally expensive and typically requires exhaustive search.We introduce a thermodynamics-inspired framework that treats activation distributions as energy-filtered physical systems and employs the free energy of activations as a principled evaluation metric.Phase-transition-like phenomena in the free-energy profile—such as extrema,inflection points,and curvature changes—yield reliable estimates of the critical pruning threshold,providing a theoretically grounded means of predicting sharp accuracy degradation.To further enhance efficiency,we propose a renormalized free energy technique that approximates full-evaluation free energy using only the activation distribution of the unpruned network.This eliminates repeated forward passes,dramatically reducing computational overhead and achieving speedups of up to 550×for MLPs.Extensive experiments across diverse vision architectures(MLP,CNN,ResNet,MobileNet,Vision Transformer)and text models(LSTM,BERT,ELECTRA,T5,GPT-2)on multiple datasets validate the generality,robustness,and computational efficiency of our approach.Overall,this work establishes a theoretically grounded and practically effective framework for activation pruning,bridging the gap between analytical understanding and efficient deployment of sparse neural networks.展开更多
Separation bubbles forming on airfoils significantly influence aerodynamic behavior,particularly at low Reynolds numbers,making their accurate prediction a critical challenge in transition modelling.This study investi...Separation bubbles forming on airfoils significantly influence aerodynamic behavior,particularly at low Reynolds numbers,making their accurate prediction a critical challenge in transition modelling.This study investigates numerical modeling of a separation bubble and the effects of airfoil thickness and camber variation on the formation of the bubble dynamics at low Reynolds numbers.The numerical results were compared with the experimental results obtained from surface pressure distribution measurements,oil flow visualisation,and surface shear measurements to analyse the detailed flow behavior.The combination of pressure and flow visualisation techniques provided complementary insights,enabling a detailed characterisation of bubble formation.The results reveal that both the thickness and camber of the airfoil significantly influence the location,length,and stability of the bubble.At low Reynolds number flows(Re=0.5×10^(5)),particularly for highly cambered profiles,closer to the leading edge,separation and long bubbles were observed.As the Reynolds number increased,the separation point shifted to the leading edge,and reattachment became more likely.In numerical studies,transition models can accurately model the bubble initiation point;however,they often fail to model the bubble reattachment points accurately.This is due to the inadequacy of models that use empirical expressions for turbulence modelling,particularly in low Reynolds number flows,in their viscous modelling.In this study,it was concluded that transition onset terms,which specifically affect bubble formation,should be modified for more accurate modeling.展开更多
Interaction between the converter and the grid may lead to harmonic oscillations.The impedance-based method is an effective way to deal with the stability issue.In this study,the impedance-based method is used to inve...Interaction between the converter and the grid may lead to harmonic oscillations.The impedance-based method is an effective way to deal with the stability issue.In this study,the impedance-based method is used to investigate the small-signal stability of a cascaded 12-pulse line-commutated converter-based high-voltage direct current(LCC-HVDC)transmission system.In the modeling part,the impedance models of the single rectifier and inverter are established respectively with consideration to the effect of frequency coupling,which has improved the accuracy of the models.Based on the models,the AC impedance models of the cascaded LCC-HVDC transmission system are established both on the rectifier and inverter side.In the stability analysis part,the stability of the system is analyzed under different working conditions.The simulation results reveal that the established impedance model can properly represent the stability of this system.The findings of this study can provide a theoretical reference for the stability design and oscillation suppression strategy of LCC-HVDC transmission systems and LCC interconnected systems.展开更多
Kinetic impact is the most practical planetary-defense technique,with momentum-transfer efficiency central to deflection design.We present a Monte Carlo photometric framework that couples ejecta sampling,dynamical evo...Kinetic impact is the most practical planetary-defense technique,with momentum-transfer efficiency central to deflection design.We present a Monte Carlo photometric framework that couples ejecta sampling,dynamical evolution,and image synthesis to compare directly with HST,LICIACube,ground-based and Lucy observations of the DART impact.Decomposing ejecta into(1)a highvelocity(~1600 m/s)plume exhibiting Na/K resonance,(2)a low-velocity(~1 m/s)conical component shaped by binary gravity and solar radiation pressure,and(3)meter-scale boulders,we quantify each component’s mass and momentum.Fitting photometric decay curves and morphological evolution yields size-velocity distributions and,via scaling laws,estimates of Dimorphos’bulk density,cratering parameters,and cohesive strength that agree with dynamical constraints.Photometric ejecta modeling therefore provides a robust route to constrain momentum enhancement and target properties,improving predictive capability for kinetic-deflection missions.展开更多
Xylogenesis,the process through which wood cells are formed,results in the long-term storage of carbon in woody biomass,making it a key component of the global carbon cycle.Understanding how environmental drivers infl...Xylogenesis,the process through which wood cells are formed,results in the long-term storage of carbon in woody biomass,making it a key component of the global carbon cycle.Understanding how environmental drivers influence xylogenesis during the growing season is therefore of great interest.However,studying shortterm drivers of wood production using xylogenetic data is complicated by the usual sampling scheme and the influence of eccentric growth,i.e.,heterogeneous growth around the stem.In this study,we improve xylogenesis research by introducing a statistical approach that explicitly considers seasonal phenology,short-term growth rates,and growth eccentricity.To this end,we developed Bayesian models of xylogenesis and compared them with a conventional method based on the use of Gompertz functions.Our results show that eccentricity generated high temporal autocorrelation between successive samples,and that explicitly taking it into account improved both the representativeness of phenology and intra-ring variability.We observed consistent short-term patterns in the model residuals,suggesting the influence of an unaccounted-for environmental variable on cell production.The proposed models offer several advantages over traditional methods,including robust confidence intervals around predictions,consistency with phenology,and reduced sensitivity to extreme observations at the end of the growing season,often linked to eccentric growth.These models also provide a benchmark for mechanistic testing of short-term drivers of wood formation.展开更多
In the converter steelmaking process,the flow dynamics is closely related to the refractory lining structure of the bath,such as hearth height-to-diameter(H/D)ratio and lining erosion at different campaign stages.The ...In the converter steelmaking process,the flow dynamics is closely related to the refractory lining structure of the bath,such as hearth height-to-diameter(H/D)ratio and lining erosion at different campaign stages.The step of pre-processing in computational fluid dynamics(CFD)simulation is time-consuming for different lining structures,and usually takes around a week per case using the traditional direct modeling method.A parametric modeling tool has been developed to quickly generate various converter structures with quality structured grids within seconds,based on Python and OpenFOAM software.CFD simulations were established and validated using hydraulic modeling to investigate the flow dynamics and lining erosion characteristics in a 100 t top–bottom combined blowing converter under different H/D ratios and campaign stages(initial,middle,and late).The results show that the average molten bath velocity is positively correlated with bath depth.An increase in bath depth extends the path length for kinetic energy transfer of combined blowing gas streams.Excessively large bath depth or diameter will deteriorate the flow pattern and result into corresponding dead zones.Furnace wall and bottom erosion intensifies at higher H/D ratios but decreases in the late campaign stages.The H/D ratio of 1.67 is recommended in the initial design stage considering the flow characteristics.In the late campaign stage,increasing the bottom-blowing flow rate and carrying out furnace maintenance operations are recommended to maintain metallurgical efficiency and lining safety.展开更多
Mount Kandil is situated in the eastern sector of the EAHP(Eastern Anatolian High Plateau),to the south of the Lesser Caucasus.The mountain lies at the westernmost end of the Aras Mountains,which extends approximately...Mount Kandil is situated in the eastern sector of the EAHP(Eastern Anatolian High Plateau),to the south of the Lesser Caucasus.The mountain lies at the westernmost end of the Aras Mountains,which extends approximately 80 km along a NW-SE axis.With a summit reaching~3214 m(a.s.l.),Mount Kandil is a stratovolcano that,like many other peaks within the EAHP and the Lesser Caucasus,experienced significant environmental changes during Late Pleistocene.Among these,glacial processes stand out as the most profound,having distinctly shaped the mountains geomorphic landscape.This study presents,for the first time,a comprehensive analysis of the glacial morphology of Mount Kandil based on multiple datasets.Field-based morphological observations indicate that an area of approximately 32.62 km^(2)has been sculpted by glacial activity.Within six glaciated regions on Mount Kandil,25 cirques and 6 glacial valleys have been identified.In addition,moraines in various locations exhibit characteristic morphologies.Furthermore,valley glaciers are inferred to have descended to altitudes as low as~2000 m.The paleoequilibrium line(p ELA)was estimated to use AABR method within GIS,yielding a mean pELA of~2730 m.Ice thickness modelling indicates that the thickness of glaciers in the Kandil Mountain valleys reaches up to~350 m.Due to its orographic extension,Mount Kandil is exposed to humid northwest winds and receives substantial frontal precipitation(~686 mm annually).The compiled geomorphic,cartographic and morphometric parameters suggest that the glaciation dynamics of Mount Kandil—situated between the Southeastern Taurus and the Lesser Caucasus—closely resemble those observed in the Lesser Caucasus.This indicates that glaciation was primarily governed by northern atmospheric systems with additional influences from southerly or westerly winds.The integrated data also underscores the role of multiple atmospheric systems in controlling the glaciation regime around the Lesser Caucasus.Additionally,findings on regional pELA question the common belief that pELA increases eastward in EAHP.展开更多
In deep coal mining,surrounding rock is subjected to both high in-situ stress and intense mining disturbances,leading to significant time-dependent behavior.Accurately capturing this behavior is essential for predicti...In deep coal mining,surrounding rock is subjected to both high in-situ stress and intense mining disturbances,leading to significant time-dependent behavior.Accurately capturing this behavior is essential for predicting long-term roadway stability,necessitating the development of a reliable constitutive creep model and numerical simulation approach.In this study,creep experiments were conducted on pre-damaged rock with varying initial damage levels to investigate the time-dependent mechanical properties.Based on the experimental results,an accelerated-creep criterion was proposed,and an elastic-viscoplastic creep damage model(EVPCD)was established that simultaneously considers the effects of time-dependent damage and instantaneous damage caused by stress disturbances on rock creep behavior.Subsequently,the effectiveness of the proposed creep model was verified using experimental data,and the secondary development of the EVPCD model was completed based on the FLAC3D platform.Following this,a long-term stability analysis method of deep surrounding rock that accounts for excavation-and mining-induced disturbances was proposed.Using the main roadway of Xutuan Coal Mine as a case study,numerical simulations were carried out to investigate the time-dependent deformation and failure characteristics of the surrounding rock following excavation and mining disturbance.Combined with on-site monitoring of the surrounding rock damage areas,the results indicate that the EVPCD outperforms the CVISC and Nishihara models in predicting the time-dependent behavior of deep surrounding rock.展开更多
Slopes are likely to fail in areas with frequent rainfall and earthquakes.The deformation characteristics of unsaturated slopes subjected to post-rainfall earthquakes are investigated using centrifuge model tests and ...Slopes are likely to fail in areas with frequent rainfall and earthquakes.The deformation characteristics of unsaturated slopes subjected to post-rainfall earthquakes are investigated using centrifuge model tests and finite element analyses.Three tests of the slope deformation under earthquake and post-rainfall earthquakes are first studied using image analysis techniques.Then,based on an elastoplastic constitutive model,numerical simulations are carried out using the finite element method and compared with the centrifuge test results.Finally,a parametric study is performed to clarify the effects of antecedent rainfall on earthquake-induced slope deformation.The results show that slope deformation caused by post-rainfall earthquakes differs from that caused by earthquakes without antecedent rainfall.The seepage flow and soil strength of the slope are affected by previous rainfall conditions,such as intensity and duration,which directly influence the slope deformation caused by the subsequent earthquake.Soil displacement and strain become greater and the slip surface is more noticeable during the post-rainfall earthquake of higher intensity.In addition,the time interval between the rainfall and the earthquake has a considerable impact on the detailed characteristics of the slope deformation,and the significant deformation occurs at the time of lowest soil strength when seepage flow reaches the lower part of the slope.Moreover,the repeated intermittent rainfall greatly affects the subsequent earthquake-induced slope deformation,the main characteristics of which are closely related to the changes in saturation and strength of the slope.However,with the prolonged time gap between each round of rainfall,the earthquake-induced slope deformation becomes insignificant.展开更多
基金the World Climate Research Programme(WCRP),Climate Variability and Predictability(CLIVAR),and Global Energy and Water Exchanges(GEWEX)for facilitating the coordination of African monsoon researchsupport from the Center for Earth System Modeling,Analysis,and Data at the Pennsylvania State Universitythe support of the Office of Science of the U.S.Department of Energy Biological and Environmental Research as part of the Regional&Global Model Analysis(RGMA)program area。
文摘In recent years,there has been an increasing need for climate information across diverse sectors of society.This demand has arisen from the necessity to adapt to and mitigate the impacts of climate variability and change.Likewise,this period has seen a significant increase in our understanding of the physical processes and mechanisms that drive precipitation and its variability across different regions of Africa.By leveraging a large volume of climate model outputs,numerous studies have investigated the model representation of African precipitation as well as underlying physical processes.These studies have assessed whether the physical processes are well depicted and whether the models are fit for informing mitigation and adaptation strategies.This paper provides a review of the progress in precipitation simulation overAfrica in state-of-the-science climate models and discusses the major issues and challenges that remain.
基金support from the National Key R&D Program of China(Grant No.2023YFB3709901)the National Natural Science Foundation of China(Grant No.U22A20171)+1 种基金China Baowu Low Carbon Metallurgy Innovation Foundation(Grant No.BWLCF202315)the High Steel Center(HSC)at North China University of Technology and University of Science and Technology Beijing,China.
文摘The dissolution of MgO-refractory into the slag had an obvious influence on the steel-slag reaction and the slag property,especially for high-aluminum steels.The dissolution behavior of MgO-refractory was investigated under various conditions,including the temperature,the initial steel composition,and the initial slag composition.A steel-slag-refractory kinetic model for high-aluminum steel was developed,which incorporated the process of MgO-refractory dissolution.The dependence of the MgO mass transfer coefficient k_(MgO)^(r)on temperature T during MgO-refractory dissolution process was established,as described by ln k_(MgO)^(r)=63,754/T+24.38524.It was indicated that the MgO dissolution rate was significantly influenced by the temperature.A higher temperature increased the dissolution rate of MgO.The initial steel composition had a slight impact on the MgO dissolution rate.Additionally,the initial slag composition strongly impacted the MgO saturation concentration and the dissolution rate.A lower initial Al_(2)O_(3)/SiO_(2)ratio increased the MgO dissolution rate.The steel-slag-refractory kinetic model accurately predicted the dissolution of MgO-refractory and the influence of dissolved MgO on the viscosity and composition change during steel-slag-refractory reactions.It was suggested that a higher temperature can hardly reduce the viscosity due to the dissolution of the MgO-refractory.
基金funded by the National Natural Science Foundation of China,grant number 52405341Foundation of National Key Laboratory of Computational Physics,grant number 6142A05QN24012+1 种基金Chongqing Science and Technology Committee,grant number CSTB2023NSCQ-MSX0363The Science and Technology Research Program of Chongqing Municipal Education Commission,grant number KJQN202301117.
文摘In materials science and engineering design,high-fidelity and high-efficiency numerical simulation has become a driving force for innovation and practical implementation.To address longstanding bottlenecks in the development of conventional material constitutive models—such as lengthy modeling cycles and difficulties in numerical implementation—this study proposes an intelligent modeling and code generation approach powered by large languagemodels.A structured knowledge base integrating constitutive theory,numerical algorithms,and UMAT(User Material)interface specifications is constructed,and a retrieval-augmented generation strategy is employed to establish an end-to-end workflow spanning experimental data parsing,constitutive model formulation,and automatic UMAT subroutine generation.Experimental results show that the method achieves high accuracy for both a classical Johnson–Cookmodel and a physics-informed neural network(PINN)model,with key parameter identification errors below 5%.Moreover,the automatically generated UMAT subroutines yield finite element simulation results in Abaqus that are highly consistent with theoretical predictions(coefficient of determination R2>0.98)while maintaining good numerical stability.This framework is currently focused on the automatic construction of rate-dependent elastoplastic material models,and its core method also provides a clear path for extending to other constitutive categories such as hyperelasticity and viscoelasticity.This work provides an effective technical route for the rapid development and reliable numerical implementation of material constitutive models,significantly advancing the intelligence level of computational mechanics research and improving engineering application efficiency.
基金supported in part by the National Key Research and Development Program of China(2021YFB2900501)in part by the Shaanxi Science and Technology Innovation Team(2023-CX-TD-03)+3 种基金in part by the Science and Technology Program of Shaanxi Province(2021GXLH-Z-038)in part by the Natural Science Foundation of Hunan Province(2023JJ40607 and 2023JJ50045)in part by the Scientific Research Foundation of Hunan Provincial Education Department(23B0713 and 24B0603)in part by the National Natural Science Foundation of China(62401371,62101275,and 62372070).
文摘(Quasi-)closed-form results for the statistical properties of unmanned aerial vehicle(UAV)airto-ground channels are derived for the first time using a novel spatial-vector-based method from a threedimensional(3-D)arbitrary-elevation one-cylinder model.The derived results include a closed-form expression for the space-time correlation function and some quasi-closed-form ones for the space-Doppler power spectrum density,the level crossing rate,and the average fading duration,which are shown to be the generalizations of those previously obtained from the two-dimensional(2-D)one-ring model and the 3-D low-elevation one-cylinder model for terrestrial mobile-to-mobile channels.The close agreements between the theoretical results and the simulations as well as the measurements validate the utility of the derived channel statistics.Based on the derived expressions,the impacts of some parameters on the channel characteristics are investigated in an effective,efficient,and explicable way,which leads to a general guideline on the manual parameter estimation from the measurement description.
基金supported by Excellent Research Group Project for Multiphase Evolution in Hyper-Gravity of the National Natural Science Foundation of China(No.52588202)。
文摘The India-Asia collision resulted in the formation of Qinghai-Tibet Plateau.Lower crustal flow model was proposed to explain the mechanism of Cenozoic tectonic deformation of Qinghai-Tibet Plateau.In this study,we propose a new approach by combining centrifugal analog modeling with numerical simulation to simulate the tectonic uplift history of the plateau based on the lower crustal flow model,and to investigate the material migration characteristics and the influence of crustal motion velocity and ductile layer viscosity on the plateau tectonic geomorphology.The models reproduce steep-sided flat-topped geomorphic features and clockwise rotation of the material at eastern Himalayan Syntaxis,verifying the rationality of the models.The results show that the greater the crustal motion velocity and the greater the ductile layer viscosity,the steeper the terrain change;and conversely,the smaller the crustal motion velocity and the smaller the ductile layer viscosity,the gentler the terrain change.This study further indicates that the weak lower crust plays an important role in the formation of geomorphic features and material migration characteristics of Qinghai-Tibet Plateau,and provides a new insight for the study of the uplift mechanism of the Tibetan Plateau.
基金support from the National Natural Science Foundation of China(Grant Nos.42277161 and 42230709).
文摘In rock engineering,natural cracks in rock masses subjected to external loads tend to initiate and propagate,leading to potential safety hazards.To investigate the effect of cracking behavior on the mechanical properties of rocks,the cracking processes of pre-cracked rocks have been extensively studied using numerical modeling methods.The peridynamics(PD)exhibits advantages over other numerical methods due to the absence of the requirements for remeshing and external crack growth criterion.However,for modeling pre-cracked rock cracking processes under impact,current PD implementations lack generally applicable rock constitutive models and impact contact models,which leads to difficulties in determining rock material parameters and efficiently calculating impact loads.This paper proposes a non-ordinary state-based peridynamics(NOSBPD)modeling method integrating the Drucker-Prager(DP)plasticity model and an efficient contact model to address the above problems.In the proposed method,the Drucker-Prager plasticity model is integrated into the NOSBPD,thereby equipping NOSBPD with the capability to accurately characterize the nonlinear stress-strain relationship inherent in rocks.An efficient contact model between particles and meshes is designed to calculate the impact loads,which is essentially a coupling method of PD with the finite element method(FEM).The effectiveness of the proposed NOSBPD modeling method is verified by comparison with other numerical methods and experiments.Experimental results indicate that the proposed method can effectively and accurately predict the 3D cracking processes of pre-cracked cracks under impact loading,and the maximum principal stress is the key driver behind wing crack formation in pre-cracked rocks.
基金co-supported by the National Natural Science Foundation of China(Nos.12102034 and 12125201)the Open Fund of State Key Laboratory of Robotics and Systems(HIT),China。
文摘In-space cable-driven manipulators exhibit several advantages,such as a large range of motion,high dexterity,and lightweight structure.However,kinematic and dynamic analysis play an essential role in designing a cable-driven manipulator.In this paper,the kinematic analysis of a type of cable-driven manipulator is performed,and a motion planning scheme is conducted to actuate this manipulator.Moreover,a flexible multi-body dynamic model of a cable-driven manipulator considering the frictional contact between the cables and pulleys is established.To describe properties such as flexibility,vibration,and variable length of the cable,this paper utilizes reducedorder beam elements of the Absolute Nodal Coordinates Formulation(ANCF)in Arbitrary Lagrangian Eulerian(ALE)framework.Additionally,a virtual element is introduced to model the contact segment in the cable-pulley system.A tension decay factor is employed to account for the friction in the contact segment.To validate the proposed method,a semi-analytical model based on D'Alembert's principle is established.Cross-verification is performed to validate the accuracy of both models.The model is further applied to simulate the rotation of the cable-driven manipulator with different structural parameters and frictional factors.The results from the analyses provide valuable guidance for the design and motion control of the in-space cable-driven manipulator.Finally,a prototype of a single module is manufactured and tested.Ground experiments are carried out to verify the kinematic and dynamic models.
文摘Predicting rock blasting outcomes in mining has been crucial since its inception.Blasting remains the most energy-and cost-efficient method for rock breaking and is often the only practical option.However,the mechanism is complex,influenced by various rock properties,explosives,and blast design parameters,making their effects difficult to quantify.Traditional stress-based models struggle with many parameters,such as stress and Poisson's ratio,which are challenging to measure in the field.Empirical models,though simpler,often oversimplify blast conditions.Both types of models are limited to simulating a few blastholes and cannot handle full-scale blasts involving hundreds of blastholes.However,modeling full-scale blasts with all blast design parameters is most required for modern mining applications.This paper presents a novel strain-based modeling approach for blasting and geomechanical applications,utilizing measurable variables such as particle velocity,strain,and displacement.By bypassing complex constitutive relations,strain-based models capture critical blasting trends and simulate full-scale blasts with full-blast design parameters with minimal calibration.The framework encompasses field strain measurements,model construction based on measurable variables,and laboratoryderived strain-failure criteria,each offering potential for future enhancement.Additionally,a standardized field test for site characterization is recommended.The approach is demonstrated through the Multiple Blasthole Fragmentation model,which simulates rock fragmentation and fragment strain during blasting,highlighting the practicality and effectiveness of strain-based modeling for multiple blasthole blasts.Moreover,this approach extends beyond blasting,with potential applications in highwall stability monitoring and other geomechanical applications.Strain-based modeling provides a simplified yet effective solution,avoiding the complexities of rock constitutive relations and field stress measurements while enabling full-blast design simulations for large-scale field blasts.
基金the funding support from the Smart Medicine and Engineering Interdisciplinary Innovation Project of Ningbo University(No.ZHYG003)。
文摘Glassy polymers are widely used in biomedical applications in a solvent environment,yet their long-term performance is governed by the competing effects of physical aging and solvent-induced plasticization.Here,we develop a constitutive model that explicitly couples the solvent concentration,structural relaxation,and mechanical response.This framework is built on a multiplicative decomposition of deformation and an Eyring-type flow rule,with structural evolution described by an effective temperature.A generalized shift factor is introduced to quantify how the solvent concentration and effective temperature jointly affect the relaxation time,thereby integrating physical aging and plasticization.The model is subsequently applied to methacrylate(MA)-based copolymer networks immersed in phosphate-buffered saline for up to nine months.Simulations accurately capture key experimental features,including the strong softening of highly swellable networks,the partial recovery due to aging,and the mitigating role of hydrophobic crosslinking in reducing solvent uptake.While the current single-mode description cannot reproduce the full relaxation spectrum,it establishes an efficient framework for predicting the long-term mechanical performance under coupled environmental and mechanical loading.This study provides a constitutive description of solvent-swollen glassy polymers,offering mechanistic insight into the interplay between plasticization and aging.Beyond biomedical MA networks,this framework establishes a foundation for predicting the long-term performance of polymer glasses under coupled aqueous environmental and mechanical loading.
基金supported by the National Key Research and Development Program of China (No. 2024YFA1308200)the National Natural Science Foundation of China (Nos. 82274009 and81973556)。
文摘Huperzine A(HupA) is a highly selective, reversible acetylcholinesterase(AChE) inhibitor that exhibits neuroprotective effects and is clinically used to manage benign memory decline.However, the specific relationship between the pharmacokinetic(PK) profile of HupA and cerebral acetylcholine(ACh) dynamics remains poorly characterized. Here, we characterize the PK-pharmacodynamic(PD) properties of HupA in rats under both physiological and pathological conditions. Following a single intramuscular injection, HupA exhibits a short halflife but rapid brain penetration, while multiple dosing significantly enhances its brain exposure. In a middle cerebral artery occlusion(MCAO) rat model, HupA demonstrates increased brain distribution. Furthermore, HupA elevates ACh concentrations across multiple brain regions, concurrently modulating several monoamine neurotransmitters. Using a minimal physiologically based pharmacokinetic-pharmacodynamic(mPBPK-PD) modeling approach,cerebral ACh dynamics were accurately predicted based on the pharmacokinetics of HupA in systemic circulation. The developed mPBPK-PD model exhibits robust predictive performance and holds potential for guiding the optimization of clinical dosing regimens and improving the therapeutic efficacy of HupA.
基金output of a research project implemented as part of the Basic Research Program at HSE University。
文摘Activation pruning reduces neural network complexity by eliminating low-importance neuron activations,yet identifying the critical pruning threshold—beyond which accuracy rapidly deteriorates—remains computationally expensive and typically requires exhaustive search.We introduce a thermodynamics-inspired framework that treats activation distributions as energy-filtered physical systems and employs the free energy of activations as a principled evaluation metric.Phase-transition-like phenomena in the free-energy profile—such as extrema,inflection points,and curvature changes—yield reliable estimates of the critical pruning threshold,providing a theoretically grounded means of predicting sharp accuracy degradation.To further enhance efficiency,we propose a renormalized free energy technique that approximates full-evaluation free energy using only the activation distribution of the unpruned network.This eliminates repeated forward passes,dramatically reducing computational overhead and achieving speedups of up to 550×for MLPs.Extensive experiments across diverse vision architectures(MLP,CNN,ResNet,MobileNet,Vision Transformer)and text models(LSTM,BERT,ELECTRA,T5,GPT-2)on multiple datasets validate the generality,robustness,and computational efficiency of our approach.Overall,this work establishes a theoretically grounded and practically effective framework for activation pruning,bridging the gap between analytical understanding and efficient deployment of sparse neural networks.
基金the Scientific and Technological Research Council of Turkey(TÜB˙ITAK)for support under project number:122M826to the Scientific Research Projects Unit of Erciyes University under contract No.:FYL-2023-13162 and FYL-2024-13701.
文摘Separation bubbles forming on airfoils significantly influence aerodynamic behavior,particularly at low Reynolds numbers,making their accurate prediction a critical challenge in transition modelling.This study investigates numerical modeling of a separation bubble and the effects of airfoil thickness and camber variation on the formation of the bubble dynamics at low Reynolds numbers.The numerical results were compared with the experimental results obtained from surface pressure distribution measurements,oil flow visualisation,and surface shear measurements to analyse the detailed flow behavior.The combination of pressure and flow visualisation techniques provided complementary insights,enabling a detailed characterisation of bubble formation.The results reveal that both the thickness and camber of the airfoil significantly influence the location,length,and stability of the bubble.At low Reynolds number flows(Re=0.5×10^(5)),particularly for highly cambered profiles,closer to the leading edge,separation and long bubbles were observed.As the Reynolds number increased,the separation point shifted to the leading edge,and reattachment became more likely.In numerical studies,transition models can accurately model the bubble initiation point;however,they often fail to model the bubble reattachment points accurately.This is due to the inadequacy of models that use empirical expressions for turbulence modelling,particularly in low Reynolds number flows,in their viscous modelling.In this study,it was concluded that transition onset terms,which specifically affect bubble formation,should be modified for more accurate modeling.
基金supported in part by the National Natural Science Foundation of China under 52125704 and 51937001.
文摘Interaction between the converter and the grid may lead to harmonic oscillations.The impedance-based method is an effective way to deal with the stability issue.In this study,the impedance-based method is used to investigate the small-signal stability of a cascaded 12-pulse line-commutated converter-based high-voltage direct current(LCC-HVDC)transmission system.In the modeling part,the impedance models of the single rectifier and inverter are established respectively with consideration to the effect of frequency coupling,which has improved the accuracy of the models.Based on the models,the AC impedance models of the cascaded LCC-HVDC transmission system are established both on the rectifier and inverter side.In the stability analysis part,the stability of the system is analyzed under different working conditions.The simulation results reveal that the established impedance model can properly represent the stability of this system.The findings of this study can provide a theoretical reference for the stability design and oscillation suppression strategy of LCC-HVDC transmission systems and LCC interconnected systems.
基金supported by the National Natural Science Foundation of China(Grant No.12272018)the National Key Basic Research Project(2022JCJQZD20600).
文摘Kinetic impact is the most practical planetary-defense technique,with momentum-transfer efficiency central to deflection design.We present a Monte Carlo photometric framework that couples ejecta sampling,dynamical evolution,and image synthesis to compare directly with HST,LICIACube,ground-based and Lucy observations of the DART impact.Decomposing ejecta into(1)a highvelocity(~1600 m/s)plume exhibiting Na/K resonance,(2)a low-velocity(~1 m/s)conical component shaped by binary gravity and solar radiation pressure,and(3)meter-scale boulders,we quantify each component’s mass and momentum.Fitting photometric decay curves and morphological evolution yields size-velocity distributions and,via scaling laws,estimates of Dimorphos’bulk density,cratering parameters,and cohesive strength that agree with dynamical constraints.Photometric ejecta modeling therefore provides a robust route to constrain momentum enhancement and target properties,improving predictive capability for kinetic-deflection missions.
基金supported by the Discovery Grants program of the Natural Sciences and Engineering Research Council of Canada(No.RGPIN-2021-03553)the Canadian Research Chair in dendroecology and dendroclimatology(CRC-2021-00368)+3 种基金the Ministère des Ressources Naturelles et des Forèts(MRNF,Contract no.142332177-D)the Natural Sciences and Engineering Research Council of Canada(Alliance Grant No.ALLRP 557148-20,obtained in partnership with the MRNF and Resolute Forest Products)the Fonds de recherche du Qu ebec–Nature et technologies(Partnership Research Program on the Contribution of the Forestry Sector to Climate Change MitigationGrant No.2022-0FC-309064)。
文摘Xylogenesis,the process through which wood cells are formed,results in the long-term storage of carbon in woody biomass,making it a key component of the global carbon cycle.Understanding how environmental drivers influence xylogenesis during the growing season is therefore of great interest.However,studying shortterm drivers of wood production using xylogenetic data is complicated by the usual sampling scheme and the influence of eccentric growth,i.e.,heterogeneous growth around the stem.In this study,we improve xylogenesis research by introducing a statistical approach that explicitly considers seasonal phenology,short-term growth rates,and growth eccentricity.To this end,we developed Bayesian models of xylogenesis and compared them with a conventional method based on the use of Gompertz functions.Our results show that eccentricity generated high temporal autocorrelation between successive samples,and that explicitly taking it into account improved both the representativeness of phenology and intra-ring variability.We observed consistent short-term patterns in the model residuals,suggesting the influence of an unaccounted-for environmental variable on cell production.The proposed models offer several advantages over traditional methods,including robust confidence intervals around predictions,consistency with phenology,and reduced sensitivity to extreme observations at the end of the growing season,often linked to eccentric growth.These models also provide a benchmark for mechanistic testing of short-term drivers of wood formation.
基金fundings of National Natural Science Foundation of China(No.52374322)Baowu Group are highly appreciated.
文摘In the converter steelmaking process,the flow dynamics is closely related to the refractory lining structure of the bath,such as hearth height-to-diameter(H/D)ratio and lining erosion at different campaign stages.The step of pre-processing in computational fluid dynamics(CFD)simulation is time-consuming for different lining structures,and usually takes around a week per case using the traditional direct modeling method.A parametric modeling tool has been developed to quickly generate various converter structures with quality structured grids within seconds,based on Python and OpenFOAM software.CFD simulations were established and validated using hydraulic modeling to investigate the flow dynamics and lining erosion characteristics in a 100 t top–bottom combined blowing converter under different H/D ratios and campaign stages(initial,middle,and late).The results show that the average molten bath velocity is positively correlated with bath depth.An increase in bath depth extends the path length for kinetic energy transfer of combined blowing gas streams.Excessively large bath depth or diameter will deteriorate the flow pattern and result into corresponding dead zones.Furnace wall and bottom erosion intensifies at higher H/D ratios but decreases in the late campaign stages.The H/D ratio of 1.67 is recommended in the initial design stage considering the flow characteristics.In the late campaign stage,increasing the bottom-blowing flow rate and carrying out furnace maintenance operations are recommended to maintain metallurgical efficiency and lining safety.
基金supported by Van Yüzüncü Yıl University,Scientific Research Projects Coordination Unit(Project No:SDK-2025-11935)Van Yüzüncü Yıl University,Scientific Research Projects Coordination Unit for supporting the study。
文摘Mount Kandil is situated in the eastern sector of the EAHP(Eastern Anatolian High Plateau),to the south of the Lesser Caucasus.The mountain lies at the westernmost end of the Aras Mountains,which extends approximately 80 km along a NW-SE axis.With a summit reaching~3214 m(a.s.l.),Mount Kandil is a stratovolcano that,like many other peaks within the EAHP and the Lesser Caucasus,experienced significant environmental changes during Late Pleistocene.Among these,glacial processes stand out as the most profound,having distinctly shaped the mountains geomorphic landscape.This study presents,for the first time,a comprehensive analysis of the glacial morphology of Mount Kandil based on multiple datasets.Field-based morphological observations indicate that an area of approximately 32.62 km^(2)has been sculpted by glacial activity.Within six glaciated regions on Mount Kandil,25 cirques and 6 glacial valleys have been identified.In addition,moraines in various locations exhibit characteristic morphologies.Furthermore,valley glaciers are inferred to have descended to altitudes as low as~2000 m.The paleoequilibrium line(p ELA)was estimated to use AABR method within GIS,yielding a mean pELA of~2730 m.Ice thickness modelling indicates that the thickness of glaciers in the Kandil Mountain valleys reaches up to~350 m.Due to its orographic extension,Mount Kandil is exposed to humid northwest winds and receives substantial frontal precipitation(~686 mm annually).The compiled geomorphic,cartographic and morphometric parameters suggest that the glaciation dynamics of Mount Kandil—situated between the Southeastern Taurus and the Lesser Caucasus—closely resemble those observed in the Lesser Caucasus.This indicates that glaciation was primarily governed by northern atmospheric systems with additional influences from southerly or westerly winds.The integrated data also underscores the role of multiple atmospheric systems in controlling the glaciation regime around the Lesser Caucasus.Additionally,findings on regional pELA question the common belief that pELA increases eastward in EAHP.
基金funded by the National Natural Science Foundation of China(Nos.52004098,U24B2041,and 52274079)the Key Research and Development Program of Henan Province(No.251111320400)+1 种基金the Key Research Project Plan for Higher Education Institutions in Henan Province(Nos.24A570006 and 25A570002)the Scientific and Technological Research Project in Henan Province(No.242102320061).
文摘In deep coal mining,surrounding rock is subjected to both high in-situ stress and intense mining disturbances,leading to significant time-dependent behavior.Accurately capturing this behavior is essential for predicting long-term roadway stability,necessitating the development of a reliable constitutive creep model and numerical simulation approach.In this study,creep experiments were conducted on pre-damaged rock with varying initial damage levels to investigate the time-dependent mechanical properties.Based on the experimental results,an accelerated-creep criterion was proposed,and an elastic-viscoplastic creep damage model(EVPCD)was established that simultaneously considers the effects of time-dependent damage and instantaneous damage caused by stress disturbances on rock creep behavior.Subsequently,the effectiveness of the proposed creep model was verified using experimental data,and the secondary development of the EVPCD model was completed based on the FLAC3D platform.Following this,a long-term stability analysis method of deep surrounding rock that accounts for excavation-and mining-induced disturbances was proposed.Using the main roadway of Xutuan Coal Mine as a case study,numerical simulations were carried out to investigate the time-dependent deformation and failure characteristics of the surrounding rock following excavation and mining disturbance.Combined with on-site monitoring of the surrounding rock damage areas,the results indicate that the EVPCD outperforms the CVISC and Nishihara models in predicting the time-dependent behavior of deep surrounding rock.
基金supported by the China Postdoctoral Science Foundation(CPSF)(Grant No.2024M762769)the Natural Science Basic Research Program of Shaanxi(Grant No.2024JC-YBQN-0333)the Postdoctoral Fellowship Program of CPSF(Grant No.GZC20232230).
文摘Slopes are likely to fail in areas with frequent rainfall and earthquakes.The deformation characteristics of unsaturated slopes subjected to post-rainfall earthquakes are investigated using centrifuge model tests and finite element analyses.Three tests of the slope deformation under earthquake and post-rainfall earthquakes are first studied using image analysis techniques.Then,based on an elastoplastic constitutive model,numerical simulations are carried out using the finite element method and compared with the centrifuge test results.Finally,a parametric study is performed to clarify the effects of antecedent rainfall on earthquake-induced slope deformation.The results show that slope deformation caused by post-rainfall earthquakes differs from that caused by earthquakes without antecedent rainfall.The seepage flow and soil strength of the slope are affected by previous rainfall conditions,such as intensity and duration,which directly influence the slope deformation caused by the subsequent earthquake.Soil displacement and strain become greater and the slip surface is more noticeable during the post-rainfall earthquake of higher intensity.In addition,the time interval between the rainfall and the earthquake has a considerable impact on the detailed characteristics of the slope deformation,and the significant deformation occurs at the time of lowest soil strength when seepage flow reaches the lower part of the slope.Moreover,the repeated intermittent rainfall greatly affects the subsequent earthquake-induced slope deformation,the main characteristics of which are closely related to the changes in saturation and strength of the slope.However,with the prolonged time gap between each round of rainfall,the earthquake-induced slope deformation becomes insignificant.
