With climate change,high-altitude areas have been frequently observed with rising temperature and humidity levels,causing an increased likelihood of collapse of ice-rich slopes and threatening downstream human settlem...With climate change,high-altitude areas have been frequently observed with rising temperature and humidity levels,causing an increased likelihood of collapse of ice-rich slopes and threatening downstream human settlements and infrastructural assets.For example,two giant glaciers collapsed in 2016 in the Aru Range,Xizang,China,killing nine herders.Thus,developing numerical methodologies for stability analysis and reproducing the collapse and subsequent movement of landslide debris is imperative for proactively managing disaster risks.This study focuses on the two collapse events within the Aru Range,to numerically analyze the pre-collapse stability of the slopes and their movement processes after collapse.Compared with previous research,this study considers the impact of various environmental factors on the temperature and stability of the two Aru glaciers,especially the heat flux caused by subglacial seepage and geothermal activity.In addition to proving similar stability between the two slopes before the collapse and simulating the positions of headwalls after collapse,this study demonstrates the need of selecting the slope region for simulation,and clarifies the influence of subglacial water flow on the positions of headwalls.Finally,this study reproduces the transport distance of the sliding body and simulates the tsunami caused by the Aru glacial debris rushing into Aru Co Lake.An effective friction coefficient of 0.10-0.11 between the glacier debris and the terrain is proposed.This provides a reference for stability analyses and collapse consequence predictions of ice-rich slopes,aiding in developing strategies for hazard mitigation.展开更多
Exploring dynamic mechanical responses and failure behaviors of hot dry rock(HDR)is significant for geothermal exploitation and stability assessment.In this study,via the split Hopkinson pressure bar(SHPB)system,a ser...Exploring dynamic mechanical responses and failure behaviors of hot dry rock(HDR)is significant for geothermal exploitation and stability assessment.In this study,via the split Hopkinson pressure bar(SHPB)system,a series of dynamic compression tests were conducted on granite treated by cyclic thermal shocks at different temperatures.We analyzed the effects of cyclic thermal shock on the thermal-related physical and dynamic mechanical behaviors of granite.Specifically,the P-wave velocity,dynamic strength,and elastic modulus of the tested granite decrease with increasing temperature and cycle number,while porosity and peak strain increase.The degradation law of dynamic mechanical properties could be described by a cubic polynomial.Cyclic thermal shock promotes shear cracks propagation,causing dynamic failure mode of granite to transition from splitting to tensile-shear composite failure,accompanied by surface spalling and debris splashing.Moreover,the thermal shock damage evolution and coupled failure mechanism of tested granite are discussed.The evolution of thermal shock damage with thermal shock cycle numbers shows an obvious S-shaped surface,featured by an exponential correlation with dynamic mechanical parameters.In addition,with increasing thermal shock temperature and cycles,granite mineral species barely change,but the length and width of thermal cracks increase significantly.The non-uniform expansion of minerals,thermal shock-induced cracking,and water-rock interaction are primary factors for deteriorating dynamic mechanical properties of granite under cyclic thermal shock.展开更多
This paper deals with the singular chemotaxis-Navier-Stokes system with indirect signal consumption n_(t)+u·▽v=△n-Х▽·(n/v▽u);v_(t)+u·▽v=△v-uw;w_(t)+u·▽w=△w-w+n;u_(t)+(u·▽)u=△u-▽P+...This paper deals with the singular chemotaxis-Navier-Stokes system with indirect signal consumption n_(t)+u·▽v=△n-Х▽·(n/v▽u);v_(t)+u·▽v=△v-uw;w_(t)+u·▽w=△w-w+n;u_(t)+(u·▽)u=△u-▽P+n▽Ф;▽·u=0,x∈Ω,t>0 in a bounded and smooth domainΩ⊂ℝ2 with no-flux/no-flux/no-flux/no-slip boundary conditions,whereΦ∈W2,∞(Ω).A recent literature[Dai F,Liu B.J Differential Equations,2023,369:115–155]has proved that for all reasonably regular initial data,the associated initial-boundary value problem possesses a global classical solution,but qualitative information on the behavior of solution has never been touched so far.In stark contrast to the positive effect of indirect signal consumption mechanism on the global solvability of system,the analysis of asymptotic behavior of solution to the system with indirect signal consumption is essentially complicated than that with direct signal consumption because the favorable coupled structure between cells and signal is broken down by the indirect signal consumption mechanism.The present study shows that the global classical solution exponentially stabilizes toward the corresponding spatially homogeneous equilibria under a smallness condition on the initial cell mass.In comparison to the previously known result concerning the uniform convergence of solution to the system with direct signal consumption,our result inter alia provides a more in-depth understanding on the asymptotic behavior of solution.展开更多
In order to improve the accuracy of the photogrammetric joint roughness coefficient(JRC)value,the present study proposed a novel method combining an autonomous shooting parameter selection algorithm with a composite e...In order to improve the accuracy of the photogrammetric joint roughness coefficient(JRC)value,the present study proposed a novel method combining an autonomous shooting parameter selection algorithm with a composite error model.Firstly,according to the depth map-based photogrammetric theory,the estimation of JRC from a three-dimensional(3D)digital surface model of rock discontinuities was presented.Secondly,an automatic shooting parameter selection algorithm was novelly proposed to establish the 3D model dataset of rock discontinuities with varying shooting parameters and target sizes.Meanwhile,the photogrammetric tests were performed with custom-built equipment capable of adjusting baseline lengths,and a total of 36 sets of JRC data was gathered via a combination of laboratory and field tests.Then,by combining the theory of point cloud coordinate computation error with the equation of JRC calculation,a composite error model controlled by the shooting parameters was proposed.This newly proposed model was validated via the 3D model dataset,demonstrating the capability to correct initially obtained JRC values solely based on shooting parameters.Furthermore,the implementation of this correction can significantly reduce errors in JRC values obtained via photographic measurement.Subsequently,our proposed error model was integrated into the shooting parameter selection algorithm,thus improving the rationality and convenience of selecting suitable shooting parameter combinations when dealing with target rock masses with different sizes.Moreover,the optimal combination of three shooting parameters was offered.JRC values resulting from various combinations of shooting parameters were verified by comparing them with 3D laser scan data.Finally,the application scope and limitations of the newly proposed approach were further addressed.展开更多
Rock masses in high-elevation or cold regions are vulnerable to the combined effects of freeze-thaw(F-T)weathering and dynamic mixed-mode loading,posing a serious threaten to the safety and stability of geotechnical e...Rock masses in high-elevation or cold regions are vulnerable to the combined effects of freeze-thaw(F-T)weathering and dynamic mixed-mode loading,posing a serious threaten to the safety and stability of geotechnical engineering.In this study,a series of dynamic fracture tests were conducted on notched semi-circular bend(NSCB)granite specimens subjected to different mixed-mode loading and F-T cycles using a split Hopkinson pressure bar(SHPB)test system.The effects of F-T treatment and dynamic mixed-mode loading on the fracture properties of granite,including effective fracture toughness,progressive fracture process,and macroscopic morphology of fracture surface,were comprehensively revealed.The experimental results suggest that the dynamic effective fracture toughness of NSCB specimens is dependent on the loading rate,particularly when the mode I loading is dominant.Additionally,the fracture toughness decreases as the number of F-T cycles increases,with an inflection point at 30 F-T cycles.All granite specimens subjected to mixed-mode loading exhibit a curved fracture path,with faster crack propagation speed and more fine cracks in specimens exposed to higher F-T cycles.Macroscopic morphology of fracture surface obtained using three-dimensional(3D)scanning indicates that the fractal dimension of the fracture surface increases with increasing F-T cycles,and the increment is more pronounced for specimens subjected to a higher mode II loading component.Moreover,this study compared the fracture resistance of F-T treated granite subjected to dynamic mixed loading using the maximum tangential stress(MTS)criterion and the generalized maximum tangential stress-based semi-analytical(SA-GMTS)criterion.Compared with the MTS criterion,the SA-GMTS criterion shows a more reasonable consistency with the experimental results,with a root mean square error within±7%.展开更多
The surrounding rock of underground space is always affected by external dynamic disturbance from the side position,such as blasting vibration from a stope at the same level or seismic waves from adjacent strata.A ser...The surrounding rock of underground space is always affected by external dynamic disturbance from the side position,such as blasting vibration from a stope at the same level or seismic waves from adjacent strata.A series of laboratory tests,numerical simulations and theoretical analyses were carried out in this study to disclose the sliding mechanism of roof rock blocks under lateral disturbance.Firstly,the experiments on trapezoidal key block under various clamping loads and disturbance were conducted,followed by numerical simulations using the fast Lagrangian analysis of continua(FLAC3D).Then,based on the conventional wave propagation model and the classical shear-slip constitutive model,a theoretical model was proposed to capture the relative displacement between blocks and the sliding displacement of the key block.The results indicate that the sliding displacement of the key block increased linearly with the disturbance energy and decreased exponentially with the clamping load when the key block was disturbed to slide(without instability).Meanwhile,when the key block was disturbed to fall,two types of instability process may appear as immediate type or delayed type.In addition,the propagation of stress waves in the block system exhibited obvious low-velocity and lowfrequency characteristics,resulting in the friction reduction effect appearing at the contact interface,which is the essential reason for the sliding of rock blocks.The results can be applied to practical underground engineering and provide valuable guidance for the early detection and prevention of rockfalling disasters.展开更多
Rock engineering is highly susceptible to cyclic loads resulting from earthquakes,quarrying or rockbursts.Acquiring the fatigue properties and failure mechanism of rocks is pivotal for long-term stability assessment o...Rock engineering is highly susceptible to cyclic loads resulting from earthquakes,quarrying or rockbursts.Acquiring the fatigue properties and failure mechanism of rocks is pivotal for long-term stability assessment of rock engineering structures.So far,significant progress has been gained on the mechanical characteristics of rocks subjected to cyclic loading.For providing a global insight of typical results and main features of rocks under cyclic loading conditions,this study comprehensively reviews the state-ofthe-art of deformation and failure mechanism and fatigue constitutive relationship of rocks subjected to cyclic loading in the past 60 years.Firstly,cyclic tests on rocks are classified into different types based on loading paths,loading parameters,loading types and environment conditions.Secondly,representative results are summarized and highlighted in terms of the fatigue response of rocks,including the deformation degradation,energy dissipation,damage evolution and failure characteristics;both laboratory testing and numerical results are presented,and various measurement techniques such as X-ray microcomputed tomography(micro-CT)and digital image correlation(DIC)are considered.Thirdly,the influences of cyclic loads on the mechanical characteristics of rocks are discussed,including the cyclic stress,frequency,amplitude and waveform.Subsequently,constitutive relationships for rocks subjected to cyclic loading are outlined,in which typical fatigue constitutive models are compared and analyzed,regarding the elastoplastic model,the internal variable model,the energy-based damage model and the discrete element-based model.Finally,some ambiguous questions and prospective research are interpreted and discussed.展开更多
Objective：To compare radiofrequency ablation （RFA） or microwave ablation （MWA） and transcatheter arterial chemoembolization （TACE） with RFA or MWA monotherapy in hepatocellular carcinoma （HCC）.Methods：A pros...Objective：To compare radiofrequency ablation （RFA） or microwave ablation （MWA） and transcatheter arterial chemoembolization （TACE） with RFA or MWA monotherapy in hepatocellular carcinoma （HCC）.Methods：A prospective,randomized,controlled trial was conducted on 94 patients with HCC ≤7 cm at a single tertiary referral center from June 2008 to June 2010 at the Department of Hepatobiliary Surgery,the Second Affiliated Hospital of Southeast University.The patients were randomly assigned into the TACERFA or TACE-MWA （combined treatment group） and the RFA-alone or MWA-alone groups （control group）.The primary end point was overall survival.The secondary end point was recurrence-free survival,and the tertiary end point was adverse effects.Results：Until the time of censor,17 patients in the TACE-RFA or TACE-MWA group had died.The median follow-up time of the patients who were still alive for the TACE-RFA or TACE-MWA group was 47.5±11.3 months （range,29 to 62 months）.The 1-,3-and 5-year overall survival for the TACE-RFA or TACE-MWA group was 93.6％,68.1％ and 61.7％,respectively.Twenty-five patients in the RFA or MWA group had died.The median follow-up time of the patients who were still alive for the RFA or MWA group was 47.0±12.9 months （range,28 to 62 months）.The 1-,3-and 5-year overall survival for the RFA or MWA group was 85.1％,59.6％ and 44.7％,respectively.The patients in the TACE-RFA or TACE-MWA group had better overall survival than the RFA or MWA group [hazard ratio （HR）,0.526; 95％ confidence interval （95％ CO,0.334-0.823; P=0.002],and showed better recurrence-free survival than the RFA or MWA group （HR,0.582; 95％ CI,0.368-0.895; P=0.008）.Conclusions：RFA or MWA combined with TACE in the treatment of HCC ≤7 cm was superior to RFA or MWA alone in improving survival by reducing arterial and portal blood flow due to TACE with iodized oil before RFA.展开更多
Precisely understanding the dynamic mechanical properties and failure modes of rocks subjected to true triaxial stress state(σ1>σ2>σ3,whereσ1,σ2,andσ3 are the major principal stress,intermediate principal ...Precisely understanding the dynamic mechanical properties and failure modes of rocks subjected to true triaxial stress state(σ1>σ2>σ3,whereσ1,σ2,andσ3 are the major principal stress,intermediate principal stress,and minor principal stress,respectively)is essential to the safety of underground engineering.However,in the laboratory,it is difficult to maintain the constant true triaxial stress state of rocks during the dynamic testing process.Herein,a numerical servo triaxial Hopkinson bar(NSTHB)was developed to study the dynamic responses of rocks confronted with a true triaxial stress state,in which lateral stresses can maintain constant.The results indicate that the dynamic strength and elastic modulus of rocks increase with the rise of intermediate principal stressσ2,while the dynamic elastic modulus is independent of the dynamic strain rate.Simulated acoustic emission distributions indicate that the intermediate principal stressσ2 dramatically affects dynamic failure modes of triaxial confined rocks.Asσ2 increases,the failure pattern switches from a single diagonal shear zone into two parallel shear zones with a small slant.Moreover,a recent triaxial Hopkinson bar experimental system using three bar pairs is also numerically established,and the measuring discrepancies are identified between the two numerical bar systems.The proposed NSTHB system provides a controllable tool for studying the dynamic triaxial behavior of rocks.展开更多
Microseismic(MS)event locations are vital aspect of MS monitoring technology used to delineate the damage zone inside the surrounding rock mass.However,complex geological conditions can impose significantly adverse ef...Microseismic(MS)event locations are vital aspect of MS monitoring technology used to delineate the damage zone inside the surrounding rock mass.However,complex geological conditions can impose significantly adverse effects on the final location results.To achieve a high-accuracy location in a complex cavern-containing structure,this study develops an MS location method using the fast marching method(FMM)with a second-order difference approach(FMM2).Based on the established velocity model with three-dimensional(3D)discrete grids,the realization of the MS location can be achieved by searching the minimum residual between the theoretical and actual first arrival times.Moreover,based on the calculation results of FMM2,the propagation paths from the MS sources to MS sensors can be obtained using the linear interpolation approach and the Runge–Kutta method.These methods were validated through a series of numerical experiments.In addition,our proposed method was applied to locate the recorded blasting and MS events that occurred during the excavation period of the underground caverns at the Houziyan hydropower station.The location results of the blasting activities show that our method can effectively reduce the location error compared with the results based on the uniform velocity model.Furthermore,the obtained MS location was verified through the occurrence of shotcrete fractures and spalling,and the monitoring results of the in-situ multipoint extensometer.Our proposed method can offer a more accurate rock fracture location and facilitate the delineation of damage zones inside the surrounding rock mass.展开更多
Accurate prediction of compressive strength of rocks relies on the rate-dependent behaviors of rocks, and correlation among the geometrical, physical, and mechanical properties of rocks. However, these properties may ...Accurate prediction of compressive strength of rocks relies on the rate-dependent behaviors of rocks, and correlation among the geometrical, physical, and mechanical properties of rocks. However, these properties may not be easy to control in laboratory experiments, particularly in dynamic compression experiments. By training three machine learning models based on the support vector machine(SVM), backpropagation neural network(BPNN), and random forest(RF) algorithms, we isolated different input parameters, such as static compressive strength, P-wave velocity, specimen dimension, grain size, bulk density, and strain rate, to identify their importance in the strength prediction. Our results demonstrated that the RF algorithm shows a better performance than the other two algorithms. The strain rate is a key input parameter influencing the performance of these models, while the others(e.g. static compressive strength and P-wave velocity) are less important as their roles can be compensated by alternative parameters. The results also revealed that the effect of specimen dimension on the rock strength can be overshadowed at high strain rates, while the effect on the dynamic increase factor(i.e. the ratio of dynamic to static compressive strength) becomes significant. The dynamic increase factors for different specimen dimensions bifurcate when the strain rate reaches a relatively high value, a clue to improve our understanding of the transitional behaviors of rocks from low to high strain rates.展开更多
The volume of influence of excavation at the right bank slope of Dagangshan Hydropower Station, southwest China, is essentially determined from microseismic monitoring, numerical modeling and conventional measurements...The volume of influence of excavation at the right bank slope of Dagangshan Hydropower Station, southwest China, is essentially determined from microseismic monitoring, numerical modeling and conventional measurements as well as in situ observations. Microseismic monitoring is a new application technique for investigating microcrackings in rock slopes. A micro- seismic monitoring network has been systematically used to monitor rock masses unloading relaxation due to continuous exca- vation of rock slope and stress redistribution caused by dam impoundment later on, and to identify and delineate the potential slippage regions since May, 2010. An important database of seismic source locations is available. The analysis of microseismic events showed a particular tempo-spatial distribution. Seismic events predominantly occurred around the upstream slope of 1180 m elevation, especially focusing on the hanging wall of fault XL316-1. Such phenomenon was interpreted by numerical modeling using RFPA-SRM code （realistic failure process analysis-strength reduction method）. By comparing microseismic activity and results of numerical simulation with in site observation and conventional measurements results, a strong correlation can he obtained between seismic source locations and excavation-induced stress distribution in the working areas. The volume of influence of the rock slope is thus determined. Engineering practices show microseismic monitoring can accurately diagnose magnitude, intensity and associated tempo-spatial characteristics of tectonic activities such as faults and unloading zones. The integrated technique combining seismic monitoring with numerical modeling, as well as in site observation and conventional surveying, leads to a better understanding of the internal effect and relationship between microseismic activity and stress field in the right bank slope from different perspectives.展开更多
Evaluating the fracture resistance of rocks is essential for predicting and preventing catastrophic failure of cracked structures in rock engineering.This investigation developed a brittle fracture model to predict te...Evaluating the fracture resistance of rocks is essential for predicting and preventing catastrophic failure of cracked structures in rock engineering.This investigation developed a brittle fracture model to predict tensile mode(mode I)failure loads of cracked rocks.The basic principle of the model is to estimate the reference crack corresponding to the fracture process zone(FPZ)based on the maximum normal strain(MNSN)ahead of the crack tip,and then use the effective crack to calculate the fracture toughness.We emphasize that the non-singular stress/strain terms should be considered in the description of the MNSN.In this way,the FPZ,non-singular terms and the biaxial stress state at the crack tip are simul-taneously considered.The principle of the model is explicit and easy to apply.To verify the proposed model,laboratory experiments were performed on a rock material using six groups of specimens.The model predicted the specimen geometry dependence of the measured fracture toughness well.More-over,the potential of the model in analyzing the size effect of apparent fracture toughness was discussed and validated through experimental data reported in the literature.The model was demonstrated su-perior to some commonly used fracture models and is an excellent tool for the safety assessment of cracked rock structures.展开更多
Accurately characterizing the mechanical responses and cracking mechanism of three-dimensional confined fractured rocks under coupled static-dynamic loading is of paramount importance for underground engineering const...Accurately characterizing the mechanical responses and cracking mechanism of three-dimensional confined fractured rocks under coupled static-dynamic loading is of paramount importance for underground engineering construction.Using a modified split Hopkinson pressure bar(SHPB)system,five groups of single-flawed specimens with the axial prestress ratio from 0 to 0.8 are tested at the strain rates in the range of 65-205 s-1under a fixed radial prestress.Our results indicate that both the dynamic strength and total strength show significant positive linear correlations with the strain rate,and the dynamic strength shows more strain rate sensitivity under higher axial prestress.The dynamic strength and corresponding failure strain decrease with increasing axial prestress,while the total strength is barely affected by the axial prestress.The dynamic elastic modulus initially increases before the axial prestress ratio reaches 0.6 and then decreases.The failure pattern of tested specimens changes from single diagonal failure to an“X”shaped conjugated failure as axial prestress increases.Furthermore,the progressive cracking processes of confined single-flawed specimens under different axial prestresses are numerically visualized by the discrete element method(DEM).Based on the displacement trend lines on both sides of cracking surface,five crack types are identified and classified in our simulation.The displacement field distributions of the DEM models reveal that the macroscopic single diagonal failure under lower axial prestress is mainly controlled by mixed tensile-shear cracks,while the“X”shaped conjugated failure under higher axial prestress is shear dominated.展开更多
Rock mass is a fractured porous medium usually subjected to complex geostress and fluid pressure simultaneously.Moreover,the properties of rock mass change in time and space due to mining-induced fractures.Therefore,i...Rock mass is a fractured porous medium usually subjected to complex geostress and fluid pressure simultaneously.Moreover,the properties of rock mass change in time and space due to mining-induced fractures.Therefore,it is always challenging to accurately measure rock mass properties.In this study,a three-dimensional(3D)microseismic(MS)data-driven damage model for jointed rock mass under hydro-mechanical coupling conditions is proposed.It is a 3D finite element model that takes seepage,damage and stress field effects into account jointly.Multiple factors(i.e.joints,water and microseismicity)are used to optimize the rock mass mechanical parameters at different scales.The model is applied in Shirengou iron mine to study the damage evolution of rock mass and assess the crown pillar stability during the transition from open-pit to underground mining.It is found that the damage pattern is mostly controlled by the structure,water and rock mass parameters.The damage pattern is evidently different from the two-dimensional result and is more consistent with the field observations.This difference is caused by the MS-derived damage acting on the rock mass.MS data are responsible for gradually correcting the damage zone,changing the direction in which it expands,and promoting it to evolve close to reality.For the crown pillar,the proposed model yields a more trustworthy safety factor.In order to guarantee the stability of the pillar,it is suggested to take waterproof and reinforcement measures in areas with a high degree of damage.展开更多
In this study,a novel testing method is proposed to characterize the dynamic shear property and failure mechanism of rocks by introducing the short beam compression(SBC)specimen into the split Hopkinson pressure bar(S...In this study,a novel testing method is proposed to characterize the dynamic shear property and failure mechanism of rocks by introducing the short beam compression(SBC)specimen into the split Hopkinson pressure bar(SHPB)system.Firstly,the stress distribution of SBC specimen is comprehensively analyzed by finite element method(FEM),and the results show that the optimal notch separation ratio of SBC specimen is C/H?0.2 to achieve successful dynamic simple-shear tests.Then,dynamic shear tests are conducted on sandstone using the SBC-SHPB method.Via careful pulse shaping technique,the dynamic force balance is guaranteed for SBC specimens,and the testing results show that the dynamic shear strength of sandstone is significantly rate-dependent.Combining the results of dynamic compression and tension tests,the failure envelopes of sandstone under different loading rates are obtained in the principle stress plane.It is found that the failure envelope of sandstone constantly expands outwards with increasing loading rate.Moreover,the energy partition of SBC specimen is quantified by virtue of high-speed digital image correlation(DIC)technique.The results show that the kinetic energy portion is non-negligible,and the shear fracture energy increases with increasing loading rate.In addition,the microscopic shear cracking mechanism of SBC specimen is analyzed by the thin section observation:the intra-granular(TG)fracture of minerals dominates the dynamic shear failure of sandstone,and the portion of TG fracture increases with increasing loading rate.This study provides a convenient and reliable method to investigate the dynamic shear property and failure mechanism of rocks.展开更多
基金the financial support from the National Natural Science Foundation of China(Grant Nos.52039007 and 42477189)the Sichuan Science and Technology Program(Grant No.2024YFHZ0341).
文摘With climate change,high-altitude areas have been frequently observed with rising temperature and humidity levels,causing an increased likelihood of collapse of ice-rich slopes and threatening downstream human settlements and infrastructural assets.For example,two giant glaciers collapsed in 2016 in the Aru Range,Xizang,China,killing nine herders.Thus,developing numerical methodologies for stability analysis and reproducing the collapse and subsequent movement of landslide debris is imperative for proactively managing disaster risks.This study focuses on the two collapse events within the Aru Range,to numerically analyze the pre-collapse stability of the slopes and their movement processes after collapse.Compared with previous research,this study considers the impact of various environmental factors on the temperature and stability of the two Aru glaciers,especially the heat flux caused by subglacial seepage and geothermal activity.In addition to proving similar stability between the two slopes before the collapse and simulating the positions of headwalls after collapse,this study demonstrates the need of selecting the slope region for simulation,and clarifies the influence of subglacial water flow on the positions of headwalls.Finally,this study reproduces the transport distance of the sliding body and simulates the tsunami caused by the Aru glacial debris rushing into Aru Co Lake.An effective friction coefficient of 0.10-0.11 between the glacier debris and the terrain is proposed.This provides a reference for stability analyses and collapse consequence predictions of ice-rich slopes,aiding in developing strategies for hazard mitigation.
基金The authors are grateful for the financial support from the National Natural Science Foundation of China(Grant Nos.52225904 and 52039007)the Natural Science Foundation of Sichuan Province(Grant No.2023NSFSC0377)supported by the New Cornerstone Science Foundation through the XPLORER PRIZE.
文摘Exploring dynamic mechanical responses and failure behaviors of hot dry rock(HDR)is significant for geothermal exploitation and stability assessment.In this study,via the split Hopkinson pressure bar(SHPB)system,a series of dynamic compression tests were conducted on granite treated by cyclic thermal shocks at different temperatures.We analyzed the effects of cyclic thermal shock on the thermal-related physical and dynamic mechanical behaviors of granite.Specifically,the P-wave velocity,dynamic strength,and elastic modulus of the tested granite decrease with increasing temperature and cycle number,while porosity and peak strain increase.The degradation law of dynamic mechanical properties could be described by a cubic polynomial.Cyclic thermal shock promotes shear cracks propagation,causing dynamic failure mode of granite to transition from splitting to tensile-shear composite failure,accompanied by surface spalling and debris splashing.Moreover,the thermal shock damage evolution and coupled failure mechanism of tested granite are discussed.The evolution of thermal shock damage with thermal shock cycle numbers shows an obvious S-shaped surface,featured by an exponential correlation with dynamic mechanical parameters.In addition,with increasing thermal shock temperature and cycles,granite mineral species barely change,but the length and width of thermal cracks increase significantly.The non-uniform expansion of minerals,thermal shock-induced cracking,and water-rock interaction are primary factors for deteriorating dynamic mechanical properties of granite under cyclic thermal shock.
文摘This paper deals with the singular chemotaxis-Navier-Stokes system with indirect signal consumption n_(t)+u·▽v=△n-Х▽·(n/v▽u);v_(t)+u·▽v=△v-uw;w_(t)+u·▽w=△w-w+n;u_(t)+(u·▽)u=△u-▽P+n▽Ф;▽·u=0,x∈Ω,t>0 in a bounded and smooth domainΩ⊂ℝ2 with no-flux/no-flux/no-flux/no-slip boundary conditions,whereΦ∈W2,∞(Ω).A recent literature[Dai F,Liu B.J Differential Equations,2023,369:115–155]has proved that for all reasonably regular initial data,the associated initial-boundary value problem possesses a global classical solution,but qualitative information on the behavior of solution has never been touched so far.In stark contrast to the positive effect of indirect signal consumption mechanism on the global solvability of system,the analysis of asymptotic behavior of solution to the system with indirect signal consumption is essentially complicated than that with direct signal consumption because the favorable coupled structure between cells and signal is broken down by the indirect signal consumption mechanism.The present study shows that the global classical solution exponentially stabilizes toward the corresponding spatially homogeneous equilibria under a smallness condition on the initial cell mass.In comparison to the previously known result concerning the uniform convergence of solution to the system with direct signal consumption,our result inter alia provides a more in-depth understanding on the asymptotic behavior of solution.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52225904 and 52039007)the Fundamental Research Funds for the Central Universities,CHD(Grant No.300102212207).
文摘In order to improve the accuracy of the photogrammetric joint roughness coefficient(JRC)value,the present study proposed a novel method combining an autonomous shooting parameter selection algorithm with a composite error model.Firstly,according to the depth map-based photogrammetric theory,the estimation of JRC from a three-dimensional(3D)digital surface model of rock discontinuities was presented.Secondly,an automatic shooting parameter selection algorithm was novelly proposed to establish the 3D model dataset of rock discontinuities with varying shooting parameters and target sizes.Meanwhile,the photogrammetric tests were performed with custom-built equipment capable of adjusting baseline lengths,and a total of 36 sets of JRC data was gathered via a combination of laboratory and field tests.Then,by combining the theory of point cloud coordinate computation error with the equation of JRC calculation,a composite error model controlled by the shooting parameters was proposed.This newly proposed model was validated via the 3D model dataset,demonstrating the capability to correct initially obtained JRC values solely based on shooting parameters.Furthermore,the implementation of this correction can significantly reduce errors in JRC values obtained via photographic measurement.Subsequently,our proposed error model was integrated into the shooting parameter selection algorithm,thus improving the rationality and convenience of selecting suitable shooting parameter combinations when dealing with target rock masses with different sizes.Moreover,the optimal combination of three shooting parameters was offered.JRC values resulting from various combinations of shooting parameters were verified by comparing them with 3D laser scan data.Finally,the application scope and limitations of the newly proposed approach were further addressed.
基金support from the National Natural Science Foundation of China(Grant Nos.52225904 and 52039007)the Natural Science Foundation of Sichuan Province(Grant No.2023NSFSC0377).
文摘Rock masses in high-elevation or cold regions are vulnerable to the combined effects of freeze-thaw(F-T)weathering and dynamic mixed-mode loading,posing a serious threaten to the safety and stability of geotechnical engineering.In this study,a series of dynamic fracture tests were conducted on notched semi-circular bend(NSCB)granite specimens subjected to different mixed-mode loading and F-T cycles using a split Hopkinson pressure bar(SHPB)test system.The effects of F-T treatment and dynamic mixed-mode loading on the fracture properties of granite,including effective fracture toughness,progressive fracture process,and macroscopic morphology of fracture surface,were comprehensively revealed.The experimental results suggest that the dynamic effective fracture toughness of NSCB specimens is dependent on the loading rate,particularly when the mode I loading is dominant.Additionally,the fracture toughness decreases as the number of F-T cycles increases,with an inflection point at 30 F-T cycles.All granite specimens subjected to mixed-mode loading exhibit a curved fracture path,with faster crack propagation speed and more fine cracks in specimens exposed to higher F-T cycles.Macroscopic morphology of fracture surface obtained using three-dimensional(3D)scanning indicates that the fractal dimension of the fracture surface increases with increasing F-T cycles,and the increment is more pronounced for specimens subjected to a higher mode II loading component.Moreover,this study compared the fracture resistance of F-T treated granite subjected to dynamic mixed loading using the maximum tangential stress(MTS)criterion and the generalized maximum tangential stress-based semi-analytical(SA-GMTS)criterion.Compared with the MTS criterion,the SA-GMTS criterion shows a more reasonable consistency with the experimental results,with a root mean square error within±7%.
基金This work was financially supported by National Key Research and Development Program of China(Grant No.2022YFC2903903)National Natural Science Foundation of China(Grant No.52304132)Yunnan Major Scientific and Technological Projects(Grant No.202202AG050014).These support is gratefully acknowledged.
文摘The surrounding rock of underground space is always affected by external dynamic disturbance from the side position,such as blasting vibration from a stope at the same level or seismic waves from adjacent strata.A series of laboratory tests,numerical simulations and theoretical analyses were carried out in this study to disclose the sliding mechanism of roof rock blocks under lateral disturbance.Firstly,the experiments on trapezoidal key block under various clamping loads and disturbance were conducted,followed by numerical simulations using the fast Lagrangian analysis of continua(FLAC3D).Then,based on the conventional wave propagation model and the classical shear-slip constitutive model,a theoretical model was proposed to capture the relative displacement between blocks and the sliding displacement of the key block.The results indicate that the sliding displacement of the key block increased linearly with the disturbance energy and decreased exponentially with the clamping load when the key block was disturbed to slide(without instability).Meanwhile,when the key block was disturbed to fall,two types of instability process may appear as immediate type or delayed type.In addition,the propagation of stress waves in the block system exhibited obvious low-velocity and lowfrequency characteristics,resulting in the friction reduction effect appearing at the contact interface,which is the essential reason for the sliding of rock blocks.The results can be applied to practical underground engineering and provide valuable guidance for the early detection and prevention of rockfalling disasters.
基金the financial support from the National Natural Science Foundation of China(Grant Nos.52039007 and 52009086)the Sichuan Province Youth Science and Technology Innovation Team,China(Grant No.2020JDTD0001).
文摘Rock engineering is highly susceptible to cyclic loads resulting from earthquakes,quarrying or rockbursts.Acquiring the fatigue properties and failure mechanism of rocks is pivotal for long-term stability assessment of rock engineering structures.So far,significant progress has been gained on the mechanical characteristics of rocks subjected to cyclic loading.For providing a global insight of typical results and main features of rocks under cyclic loading conditions,this study comprehensively reviews the state-ofthe-art of deformation and failure mechanism and fatigue constitutive relationship of rocks subjected to cyclic loading in the past 60 years.Firstly,cyclic tests on rocks are classified into different types based on loading paths,loading parameters,loading types and environment conditions.Secondly,representative results are summarized and highlighted in terms of the fatigue response of rocks,including the deformation degradation,energy dissipation,damage evolution and failure characteristics;both laboratory testing and numerical results are presented,and various measurement techniques such as X-ray microcomputed tomography(micro-CT)and digital image correlation(DIC)are considered.Thirdly,the influences of cyclic loads on the mechanical characteristics of rocks are discussed,including the cyclic stress,frequency,amplitude and waveform.Subsequently,constitutive relationships for rocks subjected to cyclic loading are outlined,in which typical fatigue constitutive models are compared and analyzed,regarding the elastoplastic model,the internal variable model,the energy-based damage model and the discrete element-based model.Finally,some ambiguous questions and prospective research are interpreted and discussed.
文摘Objective：To compare radiofrequency ablation （RFA） or microwave ablation （MWA） and transcatheter arterial chemoembolization （TACE） with RFA or MWA monotherapy in hepatocellular carcinoma （HCC）.Methods：A prospective,randomized,controlled trial was conducted on 94 patients with HCC ≤7 cm at a single tertiary referral center from June 2008 to June 2010 at the Department of Hepatobiliary Surgery,the Second Affiliated Hospital of Southeast University.The patients were randomly assigned into the TACERFA or TACE-MWA （combined treatment group） and the RFA-alone or MWA-alone groups （control group）.The primary end point was overall survival.The secondary end point was recurrence-free survival,and the tertiary end point was adverse effects.Results：Until the time of censor,17 patients in the TACE-RFA or TACE-MWA group had died.The median follow-up time of the patients who were still alive for the TACE-RFA or TACE-MWA group was 47.5±11.3 months （range,29 to 62 months）.The 1-,3-and 5-year overall survival for the TACE-RFA or TACE-MWA group was 93.6％,68.1％ and 61.7％,respectively.Twenty-five patients in the RFA or MWA group had died.The median follow-up time of the patients who were still alive for the RFA or MWA group was 47.0±12.9 months （range,28 to 62 months）.The 1-,3-and 5-year overall survival for the RFA or MWA group was 85.1％,59.6％ and 44.7％,respectively.The patients in the TACE-RFA or TACE-MWA group had better overall survival than the RFA or MWA group [hazard ratio （HR）,0.526; 95％ confidence interval （95％ CO,0.334-0.823; P=0.002],and showed better recurrence-free survival than the RFA or MWA group （HR,0.582; 95％ CI,0.368-0.895; P=0.008）.Conclusions：RFA or MWA combined with TACE in the treatment of HCC ≤7 cm was superior to RFA or MWA alone in improving survival by reducing arterial and portal blood flow due to TACE with iodized oil before RFA.
基金the financial support from the National Natural Science Foundation of China(Nos.52039007 and 52009086)the Sichuan Province Youth Science and Technology Innovation Team(No.2020JDTD0001)。
文摘Precisely understanding the dynamic mechanical properties and failure modes of rocks subjected to true triaxial stress state(σ1>σ2>σ3,whereσ1,σ2,andσ3 are the major principal stress,intermediate principal stress,and minor principal stress,respectively)is essential to the safety of underground engineering.However,in the laboratory,it is difficult to maintain the constant true triaxial stress state of rocks during the dynamic testing process.Herein,a numerical servo triaxial Hopkinson bar(NSTHB)was developed to study the dynamic responses of rocks confronted with a true triaxial stress state,in which lateral stresses can maintain constant.The results indicate that the dynamic strength and elastic modulus of rocks increase with the rise of intermediate principal stressσ2,while the dynamic elastic modulus is independent of the dynamic strain rate.Simulated acoustic emission distributions indicate that the intermediate principal stressσ2 dramatically affects dynamic failure modes of triaxial confined rocks.Asσ2 increases,the failure pattern switches from a single diagonal shear zone into two parallel shear zones with a small slant.Moreover,a recent triaxial Hopkinson bar experimental system using three bar pairs is also numerically established,and the measuring discrepancies are identified between the two numerical bar systems.The proposed NSTHB system provides a controllable tool for studying the dynamic triaxial behavior of rocks.
基金the Key Program of National Natural Science Foundation of China(52039007)for providing financial support.
文摘Microseismic(MS)event locations are vital aspect of MS monitoring technology used to delineate the damage zone inside the surrounding rock mass.However,complex geological conditions can impose significantly adverse effects on the final location results.To achieve a high-accuracy location in a complex cavern-containing structure,this study develops an MS location method using the fast marching method(FMM)with a second-order difference approach(FMM2).Based on the established velocity model with three-dimensional(3D)discrete grids,the realization of the MS location can be achieved by searching the minimum residual between the theoretical and actual first arrival times.Moreover,based on the calculation results of FMM2,the propagation paths from the MS sources to MS sensors can be obtained using the linear interpolation approach and the Runge–Kutta method.These methods were validated through a series of numerical experiments.In addition,our proposed method was applied to locate the recorded blasting and MS events that occurred during the excavation period of the underground caverns at the Houziyan hydropower station.The location results of the blasting activities show that our method can effectively reduce the location error compared with the results based on the uniform velocity model.Furthermore,the obtained MS location was verified through the occurrence of shotcrete fractures and spalling,and the monitoring results of the in-situ multipoint extensometer.Our proposed method can offer a more accurate rock fracture location and facilitate the delineation of damage zones inside the surrounding rock mass.
基金supported by National Research Foundation,Singapore under its Virtual Singapore R&D Programme (Award No.NRF2019VSG-GMS-001)。
文摘Accurate prediction of compressive strength of rocks relies on the rate-dependent behaviors of rocks, and correlation among the geometrical, physical, and mechanical properties of rocks. However, these properties may not be easy to control in laboratory experiments, particularly in dynamic compression experiments. By training three machine learning models based on the support vector machine(SVM), backpropagation neural network(BPNN), and random forest(RF) algorithms, we isolated different input parameters, such as static compressive strength, P-wave velocity, specimen dimension, grain size, bulk density, and strain rate, to identify their importance in the strength prediction. Our results demonstrated that the RF algorithm shows a better performance than the other two algorithms. The strain rate is a key input parameter influencing the performance of these models, while the others(e.g. static compressive strength and P-wave velocity) are less important as their roles can be compensated by alternative parameters. The results also revealed that the effect of specimen dimension on the rock strength can be overshadowed at high strain rates, while the effect on the dynamic increase factor(i.e. the ratio of dynamic to static compressive strength) becomes significant. The dynamic increase factors for different specimen dimensions bifurcate when the strain rate reaches a relatively high value, a clue to improve our understanding of the transitional behaviors of rocks from low to high strain rates.
基金supported by the National Natural Science Foundation of China (Nos. 50820125405, 50909013 and 50804006)the National Basic Research Program (973) of China (No. 2007CB209404)
文摘The volume of influence of excavation at the right bank slope of Dagangshan Hydropower Station, southwest China, is essentially determined from microseismic monitoring, numerical modeling and conventional measurements as well as in situ observations. Microseismic monitoring is a new application technique for investigating microcrackings in rock slopes. A micro- seismic monitoring network has been systematically used to monitor rock masses unloading relaxation due to continuous exca- vation of rock slope and stress redistribution caused by dam impoundment later on, and to identify and delineate the potential slippage regions since May, 2010. An important database of seismic source locations is available. The analysis of microseismic events showed a particular tempo-spatial distribution. Seismic events predominantly occurred around the upstream slope of 1180 m elevation, especially focusing on the hanging wall of fault XL316-1. Such phenomenon was interpreted by numerical modeling using RFPA-SRM code （realistic failure process analysis-strength reduction method）. By comparing microseismic activity and results of numerical simulation with in site observation and conventional measurements results, a strong correlation can he obtained between seismic source locations and excavation-induced stress distribution in the working areas. The volume of influence of the rock slope is thus determined. Engineering practices show microseismic monitoring can accurately diagnose magnitude, intensity and associated tempo-spatial characteristics of tectonic activities such as faults and unloading zones. The integrated technique combining seismic monitoring with numerical modeling, as well as in site observation and conventional surveying, leads to a better understanding of the internal effect and relationship between microseismic activity and stress field in the right bank slope from different perspectives.
基金he authors thank the financial support fromthe Key Program of National Natural Science Foundation of China(GrantNo.52039007)the Youth Science and Technology Innovation Research Team Fund of Sichuan Province(Grant No.2020JDTD0001).
文摘Evaluating the fracture resistance of rocks is essential for predicting and preventing catastrophic failure of cracked structures in rock engineering.This investigation developed a brittle fracture model to predict tensile mode(mode I)failure loads of cracked rocks.The basic principle of the model is to estimate the reference crack corresponding to the fracture process zone(FPZ)based on the maximum normal strain(MNSN)ahead of the crack tip,and then use the effective crack to calculate the fracture toughness.We emphasize that the non-singular stress/strain terms should be considered in the description of the MNSN.In this way,the FPZ,non-singular terms and the biaxial stress state at the crack tip are simul-taneously considered.The principle of the model is explicit and easy to apply.To verify the proposed model,laboratory experiments were performed on a rock material using six groups of specimens.The model predicted the specimen geometry dependence of the measured fracture toughness well.More-over,the potential of the model in analyzing the size effect of apparent fracture toughness was discussed and validated through experimental data reported in the literature.The model was demonstrated su-perior to some commonly used fracture models and is an excellent tool for the safety assessment of cracked rock structures.
基金financial support from the National Natural Science Foundation of China(Grant Nos.52039007and 52009086)the Youth Science and Technology Innovation Team of Sichuan Province,China(Grant No.2020JDTD0001)。
文摘Accurately characterizing the mechanical responses and cracking mechanism of three-dimensional confined fractured rocks under coupled static-dynamic loading is of paramount importance for underground engineering construction.Using a modified split Hopkinson pressure bar(SHPB)system,five groups of single-flawed specimens with the axial prestress ratio from 0 to 0.8 are tested at the strain rates in the range of 65-205 s-1under a fixed radial prestress.Our results indicate that both the dynamic strength and total strength show significant positive linear correlations with the strain rate,and the dynamic strength shows more strain rate sensitivity under higher axial prestress.The dynamic strength and corresponding failure strain decrease with increasing axial prestress,while the total strength is barely affected by the axial prestress.The dynamic elastic modulus initially increases before the axial prestress ratio reaches 0.6 and then decreases.The failure pattern of tested specimens changes from single diagonal failure to an“X”shaped conjugated failure as axial prestress increases.Furthermore,the progressive cracking processes of confined single-flawed specimens under different axial prestresses are numerically visualized by the discrete element method(DEM).Based on the displacement trend lines on both sides of cracking surface,five crack types are identified and classified in our simulation.The displacement field distributions of the DEM models reveal that the macroscopic single diagonal failure under lower axial prestress is mainly controlled by mixed tensile-shear cracks,while the“X”shaped conjugated failure under higher axial prestress is shear dominated.
基金We acknowledge the combined support from the National Natural Science Foundation of China(Grant Nos.52039007 and 42102325)Tiandi Science and Technology Co.,Ltd.(Grant No.2022-2-TD-MS012).
文摘Rock mass is a fractured porous medium usually subjected to complex geostress and fluid pressure simultaneously.Moreover,the properties of rock mass change in time and space due to mining-induced fractures.Therefore,it is always challenging to accurately measure rock mass properties.In this study,a three-dimensional(3D)microseismic(MS)data-driven damage model for jointed rock mass under hydro-mechanical coupling conditions is proposed.It is a 3D finite element model that takes seepage,damage and stress field effects into account jointly.Multiple factors(i.e.joints,water and microseismicity)are used to optimize the rock mass mechanical parameters at different scales.The model is applied in Shirengou iron mine to study the damage evolution of rock mass and assess the crown pillar stability during the transition from open-pit to underground mining.It is found that the damage pattern is mostly controlled by the structure,water and rock mass parameters.The damage pattern is evidently different from the two-dimensional result and is more consistent with the field observations.This difference is caused by the MS-derived damage acting on the rock mass.MS data are responsible for gradually correcting the damage zone,changing the direction in which it expands,and promoting it to evolve close to reality.For the crown pillar,the proposed model yields a more trustworthy safety factor.In order to guarantee the stability of the pillar,it is suggested to take waterproof and reinforcement measures in areas with a high degree of damage.
基金The authors thank the financial support from the National Natural Science Foundation of China(Grant.Nos.52039007 and 52225904)the Youth Science and Technology Innovation Research Team Fund of Sichuan Province(Grant.No.2020JDTD0001).
文摘In this study,a novel testing method is proposed to characterize the dynamic shear property and failure mechanism of rocks by introducing the short beam compression(SBC)specimen into the split Hopkinson pressure bar(SHPB)system.Firstly,the stress distribution of SBC specimen is comprehensively analyzed by finite element method(FEM),and the results show that the optimal notch separation ratio of SBC specimen is C/H?0.2 to achieve successful dynamic simple-shear tests.Then,dynamic shear tests are conducted on sandstone using the SBC-SHPB method.Via careful pulse shaping technique,the dynamic force balance is guaranteed for SBC specimens,and the testing results show that the dynamic shear strength of sandstone is significantly rate-dependent.Combining the results of dynamic compression and tension tests,the failure envelopes of sandstone under different loading rates are obtained in the principle stress plane.It is found that the failure envelope of sandstone constantly expands outwards with increasing loading rate.Moreover,the energy partition of SBC specimen is quantified by virtue of high-speed digital image correlation(DIC)technique.The results show that the kinetic energy portion is non-negligible,and the shear fracture energy increases with increasing loading rate.In addition,the microscopic shear cracking mechanism of SBC specimen is analyzed by the thin section observation:the intra-granular(TG)fracture of minerals dominates the dynamic shear failure of sandstone,and the portion of TG fracture increases with increasing loading rate.This study provides a convenient and reliable method to investigate the dynamic shear property and failure mechanism of rocks.
