In recent years, there are growing demands of representing rock mechanics and rock engineering in a digital format that can be easily managed, manipulated, analyzed and shared. The objective of this paper is to give a...In recent years, there are growing demands of representing rock mechanics and rock engineering in a digital format that can be easily managed, manipulated, analyzed and shared. The objective of this paper is to give a comprehensive review of the status quo and future trends of digitization in rock mechanics and rock engineering. Research topics essential to the process of digitization are firstly discussed, including data acquisition, data standardization, geological modeling, visualization and digital-numerical integration. New techniques that will play an important role in digitization process but require further improvement are then briefly proposed. Finally, achievements of present methods and techniques for digitization in substantial rock mechanics and rock engineering are presented.展开更多
Floor heave is a common water-induced issue of loess tunnels.However,the mechanism of tunnel floor heave induced by repeated rainfall or irrigation remains poorly understood.This study analyzes moisture variation at t...Floor heave is a common water-induced issue of loess tunnels.However,the mechanism of tunnel floor heave induced by repeated rainfall or irrigation remains poorly understood.This study analyzes moisture variation at the tunnel base and invert deformation through a 15-year survey of an operational loess tunnel.An improved Green-Ampt model was developed to capture the effects of interannual repeated infiltration on water content,permeability,and matric suction at the tunnel base.The model was validated using soil column infiltration tests.The calculated wetting front depth was incorporated into the tunnel numerical simulation to analyze the development of floor heave and assess its several influencing factors.Results show the base average water content at the tunnel base increases exponentially with operational years.After 15 years,over 50%of the monitored sections reach saturation,leading to significant invert uplift.Based on the plastic limit,saturation water content,and liquid limit of loess,floor heave is categorized into four levels with deformation thresholds at 10,25,and 50 mm.The modified Green-Ampt model is suitable for repeated infiltration conditions and offers practical guidance.Neglecting hydraulic deterioration due to repeated infiltration may lead to significant prediction errors.Numerical simulations reveal that the inverted arch uplift height,width,and height-to-width ratio increase linearly with the rise in the number of infiltration events.The uplift deformation follows Gaussian distribution in both cross and longitudinal sections,forming aΛ-shaped floor heave.The influence of infiltration duration,water content,horizontal infiltration,asymmetric infiltration,burial depth,and the number of infiltration events enhanced successively on floor heave.This model facilitates the prediction of floor heave under dynamic hydrological conditions during long-term tunnel operation.展开更多
To evaluate the effects of possible ground explosion on a shallow-buried metro tunnel, this paper attempts to analyze the dynamic responses of the operating metro tunnel in soft soil, using a widely applied explicit d...To evaluate the effects of possible ground explosion on a shallow-buried metro tunnel, this paper attempts to analyze the dynamic responses of the operating metro tunnel in soft soil, using a widely applied explicit dynamic nonlinear finite element software ANSYS/LS-DYNA. The blast induced wave propagation in the soil and the tunnel, and the von Mises effective stress and acceleration of the tunnel lining were presented, and the safety of the tunnel lining was evaluated based on the failure criterion. Besides, the parametric study of the soil was also carried out. The numerical results indicate that the upper part of the tunnel lining cross-section with directions ranging from 0° to 22.5° and horizontal distances 0 to 7 m away from the explosive center are the vulnerable areas, and the metro tunnel might be safe when tunnel depth is more than 7 m and TNT charge on the ground is no more than 500 kg, and the selection of soil parameters should be paid more attentions to conduct a more precise analysis.展开更多
Although disintegrated dolomite,widely distributed across the globe,has conventionally been a focus of research in underground engineering,the issue of slope stability issues in disintegrated dolomite strata is gainin...Although disintegrated dolomite,widely distributed across the globe,has conventionally been a focus of research in underground engineering,the issue of slope stability issues in disintegrated dolomite strata is gaining increasing prominence.This is primarily due to their unique properties,including low strength and loose structure.Current methods for evaluating slope stability,such as basic quality(BQ)and slope stability probability classification(SSPC),do not adequately account for the poor integrity and structural fragmentation characteristic of disintegrated dolomite.To address this challenge,an analysis of the applicability of the limit equilibrium method(LEM),BQ,and SSPC methods was conducted on eight disintegrated dolomite slopes located in Baoshan,Southwest China.However,conflicting results were obtained.Therefore,this paper introduces a novel method,SMRDDS,to provide rapid and accurate assessment of disintegrated dolomite slope stability.This method incorporates parameters such as disintegrated grade,joint state,groundwater conditions,and excavation methods.The findings reveal that six slopes exhibit stability,while two are considered partially unstable.Notably,the proposed method demonstrates a closer match with the actual conditions and is more time-efficient compared with the BQ and SSPC methods.However,due to the limited research on disintegrated dolomite slopes,the results of the SMRDDS method tend to be conservative as a safety precaution.In conclusion,the SMRDDS method can quickly evaluate the current situation of disintegrated dolomite slopes in the field.This contributes significantly to disaster risk reduction for disintegrated dolomite slopes.展开更多
Freeze-sealing pipe roof method is applied in the Gongbei tunnel,which causes the ground surface uplift induced by frost heave.A frost heaving prediction approach based on the coefficient of cold expansion is proposed...Freeze-sealing pipe roof method is applied in the Gongbei tunnel,which causes the ground surface uplift induced by frost heave.A frost heaving prediction approach based on the coefficient of cold expansion is proposed to simulate the ground deformation of the Gongbei tunnel.The coefficient of cold expansion in the model and the frost heaving rate from the frost heave test under the hydration condition can achieve a good correspondence making the calculation result closer to the actual engineering.The ground surface uplift along the lateral and longitudinal direction are respectively analyzed and compared with the field measured data to validate the model.The results show that a good agreement between the frost heaving prediction model and the field measured data verifies the rationality and applicability of the proposed model.The maximum uplift of the Gongbei tunnel appears at the center of the model,gradually decreasing along with the lateral and longitudinal directions.The curve in the lateral direction presents a normal distribution due to the influence of the constraint of two sides,while the one along the lateral direction shapes like a parabola with the opening downward due to the temperature field distribution.The model provides a reference for frost heaving engineering calculation.展开更多
Disintegrated dolomite slope and tunnel disasters occur frequently due to poor water stability of disintegrated dolomite,primarily in a form of seepage failure.For engineering purposes,it is critical to determine the ...Disintegrated dolomite slope and tunnel disasters occur frequently due to poor water stability of disintegrated dolomite,primarily in a form of seepage failure.For engineering purposes,it is critical to determine the seepage properties of disintegrated dolomite within the strata.However,conventional experimental methods are time-consuming and expensive and may not be effective in investigating seepage characteristics due to the heterogeneity of disintegrated dolomite.In this study,pore network model(PNM)was established by the computerized tomography(CT)scanning technology to characterize the pores.Meanwhile,the seepage and coefficient of permeability under different inlet stress conditions based on the accurate pore model were realized by linking the commercial image processing software Avizo with the commercial multi-physics modeling package Comsol.The results show that the porosities of severely and completely disintegrated dolomites are 29.17% and 45.37%,respectively.The grade of pore development increases with disintegration grade,which facilitates seepage failure.Severely and completely disintegrated dolomites have the coefficients of permeability of 9.67×10^(-7) m/s and 1.61×10^(-6) m/s,respectively.Under conventional conditions,severely and completely disintegrated dolomites undergo seepage failure above a pressure difference of 6×10^(3) Pa and 5×10^(3) Pa,respectively.These results are consistent with both in situ water pressure tests in the borehole and laboratory tests with the constant-head method,demonstrating that CT scanning is an effective method for observing fractures and pores in disintegrated dolomite for seepage evaluation.展开更多
Partition of unity based numerical manifold method can solve continuous and discontinuous problems in a unified framework with a two-cover system,i.e.,the mathematical cover and physical cover.However,renewal of the t...Partition of unity based numerical manifold method can solve continuous and discontinuous problems in a unified framework with a two-cover system,i.e.,the mathematical cover and physical cover.However,renewal of the topology of the two-cover system poses a challenge for multiple crack propagation problems and there are few references.In this study,a robust and efficient strategy is proposed to update the cover system of the numerical manifold method in simulation of multiple crack propagation problems.The proposed algorithm updates the cover system with a bottom-up process:1)identification of fractured manifold elements according to the previous and latest crack tip position;and 2)local topological update of the manifold elements,physical patches,block boundary loops,and non-persistent joint loops according to the scenario classification of the propagating crack.The proposed crack tracking strategy and classification of the renewal cases promote a robust and efficient cover renewal algorithm for multiple crack propagation analysis.Three crack propagation examples show that the proposed algorithm performs well in updating the cover system.This cover renewal methodology can be extended for numerical manifold method with polygonal mathematical covers.展开更多
The widespread utilisation of tunnel boring machines(TBMs)in underground construction engineering requires a detailed investigation of the cutter-rock interaction.In this paper,we conduct a series of largescale standi...The widespread utilisation of tunnel boring machines(TBMs)in underground construction engineering requires a detailed investigation of the cutter-rock interaction.In this paper,we conduct a series of largescale standing rotary cutting tests on granite in conjunction with high-fidelity numerical simulations based on a particle-type discrete element method(DEM)to explore the effects of key cutting parameters on the TBM cutter performance and the distribution of cutter-rock contact stresses.The assessment results of cutter performance obtained from the cutting tests and numerical simulations reveal similar dependencies on the key cutting parameters.More specifically,the normal and rolling forces exhibit a positive correlation with penetration but are slightly influenced by the cutting radius.In contrast,the side force decreases as the cutting radius increases.Additionally,the side force shows a positive relationship with the penetration for smaller cutting radii but tends to become negative as the cutting radius increases.The cutter's relative effectiveness in rock breaking is significantly impacted by the penetration but shows little dependency on the cutting radius.Consequently,an optimal penetration is identified,leading to a low boreability index and specific energy.A combined Hertz-Weibull function is developed to fit the cutter-rock contact stress distribution obtained in DEM simulations,whereby an improved CSM(Colorado School of Mines)model is proposed by replacing the original monotonic cutting force distribution with this combined Hertz-Weibull model.The proposed model outperforms the original CSM model as demonstrated by a comparison of the estimated cutting forces with those from the tests/simulations.The findings from this work that advance our understanding of TBM cutter performance have important implications for improving the efficiency and reliability of TBM tunnelling in granite.展开更多
Due to the large number of finite element mesh generated,it is difficult to use full-scale model to simulate largesection underground engineering,especially considering the coupling effect.A regional model is attempte...Due to the large number of finite element mesh generated,it is difficult to use full-scale model to simulate largesection underground engineering,especially considering the coupling effect.A regional model is attempted to achieve this simulation.A variable boundary condition method for hybrid regional model is proposed to realize the numerical simulation of large-section tunnel construction.Accordingly,the balance of initial ground stress under asymmetric boundary conditions achieves by applying boundary conditions step by step with secondary development ofDynaflowscripts,which is the key issue of variable boundary conditionmethod implementation.In this paper,Gongbei tunnel based on hybrid regional model involvingmulti-field coupling is simulated.Meanwhile,the variable boundary condition method for regional model is verified against model initialization and the ground deformation due to tunnel excavation is predicted via the proposed hybrid regional model.Compared with the monitoring data of actual engineering,the results indicated that the hybrid regional model has a good prediction effect.展开更多
Under the conditions of underground coal gasification(UCG),the pore structure and mechanical properties of overlying rocks are crucial for the stability of surrounding strata.Therefore,it is of great significance to i...Under the conditions of underground coal gasification(UCG),the pore structure and mechanical properties of overlying rocks are crucial for the stability of surrounding strata.Therefore,it is of great significance to investigate the thermal damage of surrounding rocks.In this study,the variations in apparent characteristics,mass,P-wave velocity,porosity and strength of sandstone after high temperature from room temperature to 800◦C.The results indicate that the physical and mechanical properties of sandstone exhibit different degrees of variation after being subjected to high temperatures.Macropores complexity peaks at 400℃ then weakens due to crack network simplification,while mesopores regularity linearly increases with cementation loss.Porosity growth rate shows three-stage escalation through thresholds at 200◦C and 600℃,rising from 3.81%(25℃)to 18.70%(800℃).The mechanism of sandstone damage caused by high temperatures has been explored based on microscopic thin-section analysis,and the findings of the research can offer insights for the evaluation of rock damage in underground coal gasification.展开更多
Particle crushing commonly occurs in granular materials and affects their structures and mechanical properties.Unlike idealized particles in experimental single particle crushing tests with two loading points,natural ...Particle crushing commonly occurs in granular materials and affects their structures and mechanical properties.Unlike idealized particles in experimental single particle crushing tests with two loading points,natural particles are crushed under multicontact loading.To date,the criteria and patterns of par-ticle crushing under multicontact conditions are not fully understood.By using the three-dimensional discrete element method,this report explores the effect of multicontact loading on the crushing criterion of a single particle,the crushing pattern,and the relationship between the particle crushing strength and loading distribution.The particles are modelled as aggregates of glued Voronoi polyhedra.The numer-ical results indicate that the logarithm of the mean principal stress has a good linear correlation with the coordination number.For a specific coordination number,the number of child particles presents a significant normal distribution.For a specific number of child particles,the volumes of child particles can be statistically described as normal or gamma distribution.Three typical models are proposed to qual-itatively analyse the relationship between the loading distribution and crushing strength.The relevant conclusions can be helpful in engineering practice and in further studies on crushable granular materials via the discrete element method.展开更多
In discrete element method(DEM)simulations,multi-sphere(MS)clumped and convex particles are two main particle models that are used to study the mechanical behaviours of granular materials.Of interest is the evaluation...In discrete element method(DEM)simulations,multi-sphere(MS)clumped and convex particles are two main particle models that are used to study the mechanical behaviours of granular materials.Of interest is the evaluation of the effect of multiple contacts between clumped particles or single contacts between convex particles on the mechanical behaviours of granular materials.In this context,a series of drained triaxial compression tests were conducted on convex true(CT)ellipsoids and MS ellipsoids with aspect ratios(ARs)ranging from 1.0-2.0.The microscale results indicate that at a given AR,the critical friction angleφ_(c)changes with the particle type,whereas the peak friction angleφ_(p)is nearly independent of the particle type.The anisotropic analysis provides underlying mechanisms of the shear strength evolution from two perspectives.First,the anisotropies of granular materials are essential to shear strength as the deviatoric(q)-to-effective mean(p′)stress ratio can be expressed as the sum of the anisotropies,i.e.,q/p'≈0.4a_(c)+0.4a_(n)+0.6a_(t),where ac,an and at are the normal contact anisotropy,normal contact force anisotropy and tangential contact force anisotropy,respectively.For all samples,a_(c)and a_(n)underpin the shear strength and are influenced by the particle type.The similarφ_(p)displayed by the CT and MS ellipsoids does not translate to similar a_(n)and a_(c)but similar a_(c)+a_(n)for the two particle types.In addition,owing to their larger a_(c)+a_(n),the CT ellipsoids have a higherφ_(c)than the MS ellipsoids.Second,there is a satisfactory linear relationship between q/p'and ac within strong and non-sliding(sn)contacts a_(c)^(sn)(i.e.,q/p′=ka_(c)^(sn)),where k is the fitting parameter.Accordingly,in the peak state,the subtle difference in shear strength is attributed to the greater acsn in the CT ellipsoids than in the MS ellipsoids that is counteracted by the smaller k.However,in the critical state,the greater difference in a_(c)^(sn)between the CT and MS ellipsoids is partially offset by the smaller difference in k,causing a higherφ_(c)in the CT ellipsoids than in the MS ellipsoids.展开更多
Multi-sphere clumps are commonly used to simulate non-spherical particles in discrete element method simulations.It is of interest whether the degree of local non-convexity λ affects the mechanical behaviour of granu...Multi-sphere clumps are commonly used to simulate non-spherical particles in discrete element method simulations.It is of interest whether the degree of local non-convexity λ affects the mechanical behaviour of granular materials with the same non-convexity η.A series of discrete-element-method biaxial shear tests are conducted on rough particle packings with rη=0.075 and different λ values(ranging from 0.134 to 0.770).The microscale results show that the contact type changes with an increase in λ.However,the critical strength is independent of λ.The evaluation of the contributions of different contact types to the critical shear strength and a detailed analysis of the anisotropies help clarify the microscopic mechanisms that result in the independence of the critical shear strength from λ.展开更多
The soil-rock mixture(SRM) is highly heterogeneous. Before carrying out numerical analysis,a structure model should be generated. A reliable way to obtain such structure is by generating random aggregate structure bas...The soil-rock mixture(SRM) is highly heterogeneous. Before carrying out numerical analysis,a structure model should be generated. A reliable way to obtain such structure is by generating random aggregate structure based on random sequential addition(RSA). The classical RSA is neither efficient nor robust since valid positions to place new inclusions are formulated by trial, which involves repetitive overlapping tests. In this paper, the algorithm of Entrance block between block A and B(EAB)is synergized with background mesh to redesign RSA so that permissible positions to place new inclusions can be predicted,resulting in dramatic improvement in efficiency and robustness.展开更多
Discrete element modeling was used to investigate the effect of particle size distribution on the small strain shear stiffness of granular soils and explore the fundamental mechanism controlling this small strain shea...Discrete element modeling was used to investigate the effect of particle size distribution on the small strain shear stiffness of granular soils and explore the fundamental mechanism controlling this small strain shear stiffness at the particle level. The results indicate that the mean particle size has a negligible effect on the small strain shear modulus. The observed increase of the shear modulus with increasing particle size is caused by a scale effect. It is suggested that the ratio of sample size to the mean particle size should be larger than 11.5 to avoid this possible scale effect. At the same confining pressure and void ratio, the small strain shear modulus decreases as the coefficient of uniformity of the soil increases. The Poisson＇s ratio decreases with decreasing void ratio and increasing confining pressure instead of being constant as is commonly assumed. Microscopic analyses indicate that the small strain shear stiffness and Poisson＇s ratio depend uniquely on the soil＇s coordination number.展开更多
The hydrogeologic conditions in Tianjin,China,are characterized by a system of multiple alternating aquifers and aquitards and a high groundwater level.Generally,the groundwater system to a depth of 70 m includes one ...The hydrogeologic conditions in Tianjin,China,are characterized by a system of multiple alternating aquifers and aquitards and a high groundwater level.Generally,the groundwater system to a depth of 70 m includes one phreatic aquifer and three semiconfined aquifers(AqI–AqIII).The distribution of aquifers varies greatly.Accurate estimates of aquifer parameters from pumping tests are essential in geotechnical engineering practice.In this study,the applicability of different models for estimating hydraulic parameters in a multiaquifer system is analyzed.To avoid errors in the graphic-analytical process and to enhance the efficiency and accuracy,a hybrid algorithm,called GALMA(a genetic algorithm(GA)combined with the Levenberg–Marquardt(LM)algorithm),is used with the Neuman and Witherspoon model and ratio method.The estimated values in the second and third semiconfined aquifers(AqII and AqIII)from various solutions are compared with the measured data.Further validations based on the prediction of the drawdown values are performed via a three dimensional numerical simulation.展开更多
The critical solid fraction(ϕJ),which marks the transition between the solid and liquid phases in the jamming diagram,is influenced by several factors.In this study,the dependency ofϕJ on strain rate and boundary cond...The critical solid fraction(ϕJ),which marks the transition between the solid and liquid phases in the jamming diagram,is influenced by several factors.In this study,the dependency ofϕJ on strain rate and boundary conditions is examined through discrete element method simulations considering a frictionless polydisperse granular system.Different approaches are used to determineϕJ.The observed boundary effect is due to the nonuniform solid fraction distribution induced by the clustering of particles close to rigid-wall boundaries at high compression rates.The solid fraction distribution within the sample in the rigid-wall simulations approaches that in the periodic-boundary simulations as the compression rate decreases.With increasing compression rate,the major force transmission network contains fewer mechanically stable particles and a less stable force transmission network.This causes jamming of the granular assembly at a lower solid fraction.These force transmission networks,however,are fragile and disintegrate quickly upon relaxation.展开更多
Significant movement of in-situ retaining walls is usually assumed to begin with bulk excavation.However,an increasing number of case studies show that lowering the pore water pressures inside a diaphragm wall-type ba...Significant movement of in-situ retaining walls is usually assumed to begin with bulk excavation.However,an increasing number of case studies show that lowering the pore water pressures inside a diaphragm wall-type basement enclosure prior to bulk excavation can cause wall movements in the order of some centimeters.This paper describes the results of a laboratory-scale experiment carried out to explore mechanisms of in situ retaining wall movement associated with dewatering inside the enclosure prior to bulk excavation.Dewatering reduces the pore water pressures inside the enclosure more than outside,resulting in the wall moving as an unpropped cantilever supported only by the soil.Lateral effective stresses in the shallow soil behind the wall are reduced,while lateral effective stresses in front of the wall increase.Although the associated lateral movement was small in the laboratory experiment,the movement could be proportionately larger in the field with a less stiff soil and a potentially greater dewatered depth.The implementation of a staged dewatering system,coupled with the potential for phased excavation and propping strategies,can effectively mitigate dewatering-induced wall and soil movements.This approach allows for enhanced stiffness of the wall support system,which can be dynamically adjusted based on real-time displacement monitoring data when necessary.展开更多
Underground debris flows,arising from the complex interplay of anthropogenic activities and rainfall-induced hydromechanical processes,present significant geotechnical hazards that remain poorly understood due to thei...Underground debris flows,arising from the complex interplay of anthropogenic activities and rainfall-induced hydromechanical processes,present significant geotechnical hazards that remain poorly understood due to their hidden nature and dynamic multiphase triggers.Focusing on underground debris flow in a mining area in Southwest China,this study advances an integrated framework combining air-ground transient electromagnetic method(AGTEM)and computational fluid dynamics coupled with the discrete element method(CFD-DEM),revealing the migration mechanism in which microscale multiphase hydraulic erosion drives the macroscopic initiation of underground debris flow.Key findings include:(1)The identification of three transport phases(rapid erosion,slow erosion,and stabilization)provides actionable thresholds for monitoring and mitigation.(2)The coupled feedback between hydraulic conductivity anisotropy and the formation of preferential flow is the primary driver of large-scale debris transportation.(3)Linking mining-induced seismic energy to vibration-induced liquefaction via DEM simulations offers a physics-based explanation for flow mobilization triggers.The integrated geophysical-numerical framework offers new capabilities for predicting initiation thresholds and developing physics-based mitigation strategies in mining-affected terrains.展开更多
文摘In recent years, there are growing demands of representing rock mechanics and rock engineering in a digital format that can be easily managed, manipulated, analyzed and shared. The objective of this paper is to give a comprehensive review of the status quo and future trends of digitization in rock mechanics and rock engineering. Research topics essential to the process of digitization are firstly discussed, including data acquisition, data standardization, geological modeling, visualization and digital-numerical integration. New techniques that will play an important role in digitization process but require further improvement are then briefly proposed. Finally, achievements of present methods and techniques for digitization in substantial rock mechanics and rock engineering are presented.
基金financially supported by the National Natural Science Foundation of China(grant numbers:42077265,41927806)the Science and Technology Project of Gansu Provincial Department of Transportation(grant number:2023-12)the Gansu Science and Technology Department of Key Projects(grant number:21YF5FA002)。
文摘Floor heave is a common water-induced issue of loess tunnels.However,the mechanism of tunnel floor heave induced by repeated rainfall or irrigation remains poorly understood.This study analyzes moisture variation at the tunnel base and invert deformation through a 15-year survey of an operational loess tunnel.An improved Green-Ampt model was developed to capture the effects of interannual repeated infiltration on water content,permeability,and matric suction at the tunnel base.The model was validated using soil column infiltration tests.The calculated wetting front depth was incorporated into the tunnel numerical simulation to analyze the development of floor heave and assess its several influencing factors.Results show the base average water content at the tunnel base increases exponentially with operational years.After 15 years,over 50%of the monitored sections reach saturation,leading to significant invert uplift.Based on the plastic limit,saturation water content,and liquid limit of loess,floor heave is categorized into four levels with deformation thresholds at 10,25,and 50 mm.The modified Green-Ampt model is suitable for repeated infiltration conditions and offers practical guidance.Neglecting hydraulic deterioration due to repeated infiltration may lead to significant prediction errors.Numerical simulations reveal that the inverted arch uplift height,width,and height-to-width ratio increase linearly with the rise in the number of infiltration events.The uplift deformation follows Gaussian distribution in both cross and longitudinal sections,forming aΛ-shaped floor heave.The influence of infiltration duration,water content,horizontal infiltration,asymmetric infiltration,burial depth,and the number of infiltration events enhanced successively on floor heave.This model facilitates the prediction of floor heave under dynamic hydrological conditions during long-term tunnel operation.
基金Supported by the National Natural Science Foundation of China (40874074, 50950110347)the National High Technology Research and Development Program (863 Program) of China (2006AA11ZAA8)Shanghai Science and Technology Development Funds (07ZR14117)
文摘To evaluate the effects of possible ground explosion on a shallow-buried metro tunnel, this paper attempts to analyze the dynamic responses of the operating metro tunnel in soft soil, using a widely applied explicit dynamic nonlinear finite element software ANSYS/LS-DYNA. The blast induced wave propagation in the soil and the tunnel, and the von Mises effective stress and acceleration of the tunnel lining were presented, and the safety of the tunnel lining was evaluated based on the failure criterion. Besides, the parametric study of the soil was also carried out. The numerical results indicate that the upper part of the tunnel lining cross-section with directions ranging from 0° to 22.5° and horizontal distances 0 to 7 m away from the explosive center are the vulnerable areas, and the metro tunnel might be safe when tunnel depth is more than 7 m and TNT charge on the ground is no more than 500 kg, and the selection of soil parameters should be paid more attentions to conduct a more precise analysis.
基金supported by the National Natural Science Foundation of China(Grant No.42162026)the Applied Basic Research Foundation of Yunnan Province(Grant No.202201AT070083).
文摘Although disintegrated dolomite,widely distributed across the globe,has conventionally been a focus of research in underground engineering,the issue of slope stability issues in disintegrated dolomite strata is gaining increasing prominence.This is primarily due to their unique properties,including low strength and loose structure.Current methods for evaluating slope stability,such as basic quality(BQ)and slope stability probability classification(SSPC),do not adequately account for the poor integrity and structural fragmentation characteristic of disintegrated dolomite.To address this challenge,an analysis of the applicability of the limit equilibrium method(LEM),BQ,and SSPC methods was conducted on eight disintegrated dolomite slopes located in Baoshan,Southwest China.However,conflicting results were obtained.Therefore,this paper introduces a novel method,SMRDDS,to provide rapid and accurate assessment of disintegrated dolomite slope stability.This method incorporates parameters such as disintegrated grade,joint state,groundwater conditions,and excavation methods.The findings reveal that six slopes exhibit stability,while two are considered partially unstable.Notably,the proposed method demonstrates a closer match with the actual conditions and is more time-efficient compared with the BQ and SSPC methods.However,due to the limited research on disintegrated dolomite slopes,the results of the SMRDDS method tend to be conservative as a safety precaution.In conclusion,the SMRDDS method can quickly evaluate the current situation of disintegrated dolomite slopes in the field.This contributes significantly to disaster risk reduction for disintegrated dolomite slopes.
基金supported by the financial support from National Natural Science Foundation of China(No.51478340)Natural Science Foundation of Jiangsu Province(No.BK20200707)+4 种基金The Natural Science Foundation of the Jiangsu Higher Education Institutions of China(No.20KJB560029)China Postdoctoral Science Foundation(No.2020M671670)Key Laboratory of Soft Soils and Geoenvironmental Engineering(Zhejiang University)Ministry of Education(No.2020P04)the support above is gratefully acknowledged.
文摘Freeze-sealing pipe roof method is applied in the Gongbei tunnel,which causes the ground surface uplift induced by frost heave.A frost heaving prediction approach based on the coefficient of cold expansion is proposed to simulate the ground deformation of the Gongbei tunnel.The coefficient of cold expansion in the model and the frost heaving rate from the frost heave test under the hydration condition can achieve a good correspondence making the calculation result closer to the actual engineering.The ground surface uplift along the lateral and longitudinal direction are respectively analyzed and compared with the field measured data to validate the model.The results show that a good agreement between the frost heaving prediction model and the field measured data verifies the rationality and applicability of the proposed model.The maximum uplift of the Gongbei tunnel appears at the center of the model,gradually decreasing along with the lateral and longitudinal directions.The curve in the lateral direction presents a normal distribution due to the influence of the constraint of two sides,while the one along the lateral direction shapes like a parabola with the opening downward due to the temperature field distribution.The model provides a reference for frost heaving engineering calculation.
基金supported by the National Natural Science Foundation of China(Grant No.42162026)the Basic Research Program in Yunnan Province,China(Grant No.202401AT070328)the Young Talents Project of“Xingdian Talent Support Program"in Yunnan Province,China(Grant No.YNWR-QNBJ-2020-019).
文摘Disintegrated dolomite slope and tunnel disasters occur frequently due to poor water stability of disintegrated dolomite,primarily in a form of seepage failure.For engineering purposes,it is critical to determine the seepage properties of disintegrated dolomite within the strata.However,conventional experimental methods are time-consuming and expensive and may not be effective in investigating seepage characteristics due to the heterogeneity of disintegrated dolomite.In this study,pore network model(PNM)was established by the computerized tomography(CT)scanning technology to characterize the pores.Meanwhile,the seepage and coefficient of permeability under different inlet stress conditions based on the accurate pore model were realized by linking the commercial image processing software Avizo with the commercial multi-physics modeling package Comsol.The results show that the porosities of severely and completely disintegrated dolomites are 29.17% and 45.37%,respectively.The grade of pore development increases with disintegration grade,which facilitates seepage failure.Severely and completely disintegrated dolomites have the coefficients of permeability of 9.67×10^(-7) m/s and 1.61×10^(-6) m/s,respectively.Under conventional conditions,severely and completely disintegrated dolomites undergo seepage failure above a pressure difference of 6×10^(3) Pa and 5×10^(3) Pa,respectively.These results are consistent with both in situ water pressure tests in the borehole and laboratory tests with the constant-head method,demonstrating that CT scanning is an effective method for observing fractures and pores in disintegrated dolomite for seepage evaluation.
基金Project(51321065,51479191,11672360)supported by the National Natural Science Foundation of China。
文摘Partition of unity based numerical manifold method can solve continuous and discontinuous problems in a unified framework with a two-cover system,i.e.,the mathematical cover and physical cover.However,renewal of the topology of the two-cover system poses a challenge for multiple crack propagation problems and there are few references.In this study,a robust and efficient strategy is proposed to update the cover system of the numerical manifold method in simulation of multiple crack propagation problems.The proposed algorithm updates the cover system with a bottom-up process:1)identification of fractured manifold elements according to the previous and latest crack tip position;and 2)local topological update of the manifold elements,physical patches,block boundary loops,and non-persistent joint loops according to the scenario classification of the propagating crack.The proposed crack tracking strategy and classification of the renewal cases promote a robust and efficient cover renewal algorithm for multiple crack propagation analysis.Three crack propagation examples show that the proposed algorithm performs well in updating the cover system.This cover renewal methodology can be extended for numerical manifold method with polygonal mathematical covers.
基金supported by the National Natural Science Foundation of China(Grant Nos.52278407 and 52378407)the China Postdoctoral Science Foundation(Grant No.2023M732670)the support by the Postdoctoral Fellowship Program of China Postdoctoral Science Foundation.
文摘The widespread utilisation of tunnel boring machines(TBMs)in underground construction engineering requires a detailed investigation of the cutter-rock interaction.In this paper,we conduct a series of largescale standing rotary cutting tests on granite in conjunction with high-fidelity numerical simulations based on a particle-type discrete element method(DEM)to explore the effects of key cutting parameters on the TBM cutter performance and the distribution of cutter-rock contact stresses.The assessment results of cutter performance obtained from the cutting tests and numerical simulations reveal similar dependencies on the key cutting parameters.More specifically,the normal and rolling forces exhibit a positive correlation with penetration but are slightly influenced by the cutting radius.In contrast,the side force decreases as the cutting radius increases.Additionally,the side force shows a positive relationship with the penetration for smaller cutting radii but tends to become negative as the cutting radius increases.The cutter's relative effectiveness in rock breaking is significantly impacted by the penetration but shows little dependency on the cutting radius.Consequently,an optimal penetration is identified,leading to a low boreability index and specific energy.A combined Hertz-Weibull function is developed to fit the cutter-rock contact stress distribution obtained in DEM simulations,whereby an improved CSM(Colorado School of Mines)model is proposed by replacing the original monotonic cutting force distribution with this combined Hertz-Weibull model.The proposed model outperforms the original CSM model as demonstrated by a comparison of the estimated cutting forces with those from the tests/simulations.The findings from this work that advance our understanding of TBM cutter performance have important implications for improving the efficiency and reliability of TBM tunnelling in granite.
基金supported by the financial support from National Natural Sci-ence Foundation of China(No.51478340)Natural Science Foundation of Jiangsu Province(No.BK20200707)+2 种基金Natural Science Foundation of the Jiangsu Higher Education Institutions of China(No.20KJB560029)China Postdoctoral Science Foundation(No.2020M671670)Key Laboratory of Soft Soils and Geoenvironmental Engineering(Zhejiang University),Ministry of Education(No.2020P04).
文摘Due to the large number of finite element mesh generated,it is difficult to use full-scale model to simulate largesection underground engineering,especially considering the coupling effect.A regional model is attempted to achieve this simulation.A variable boundary condition method for hybrid regional model is proposed to realize the numerical simulation of large-section tunnel construction.Accordingly,the balance of initial ground stress under asymmetric boundary conditions achieves by applying boundary conditions step by step with secondary development ofDynaflowscripts,which is the key issue of variable boundary conditionmethod implementation.In this paper,Gongbei tunnel based on hybrid regional model involvingmulti-field coupling is simulated.Meanwhile,the variable boundary condition method for regional model is verified against model initialization and the ground deformation due to tunnel excavation is predicted via the proposed hybrid regional model.Compared with the monitoring data of actual engineering,the results indicated that the hybrid regional model has a good prediction effect.
基金the National Natural Science Foundation of China(No.42307244 and 42230704)Xuzhou Research Program Youth Science and Technology Project(KC23044)Open Fund of State Key Laboratory of Geological Disaster Prevention and Geological Environment Protection(SKLGP2024K022).
文摘Under the conditions of underground coal gasification(UCG),the pore structure and mechanical properties of overlying rocks are crucial for the stability of surrounding strata.Therefore,it is of great significance to investigate the thermal damage of surrounding rocks.In this study,the variations in apparent characteristics,mass,P-wave velocity,porosity and strength of sandstone after high temperature from room temperature to 800◦C.The results indicate that the physical and mechanical properties of sandstone exhibit different degrees of variation after being subjected to high temperatures.Macropores complexity peaks at 400℃ then weakens due to crack network simplification,while mesopores regularity linearly increases with cementation loss.Porosity growth rate shows three-stage escalation through thresholds at 200◦C and 600℃,rising from 3.81%(25℃)to 18.70%(800℃).The mechanism of sandstone damage caused by high temperatures has been explored based on microscopic thin-section analysis,and the findings of the research can offer insights for the evaluation of rock damage in underground coal gasification.
基金supported by the National Natural Science Foundation of China(grant Nos.51809292 and 51978531).
文摘Particle crushing commonly occurs in granular materials and affects their structures and mechanical properties.Unlike idealized particles in experimental single particle crushing tests with two loading points,natural particles are crushed under multicontact loading.To date,the criteria and patterns of par-ticle crushing under multicontact conditions are not fully understood.By using the three-dimensional discrete element method,this report explores the effect of multicontact loading on the crushing criterion of a single particle,the crushing pattern,and the relationship between the particle crushing strength and loading distribution.The particles are modelled as aggregates of glued Voronoi polyhedra.The numer-ical results indicate that the logarithm of the mean principal stress has a good linear correlation with the coordination number.For a specific coordination number,the number of child particles presents a significant normal distribution.For a specific number of child particles,the volumes of child particles can be statistically described as normal or gamma distribution.Three typical models are proposed to qual-itatively analyse the relationship between the loading distribution and crushing strength.The relevant conclusions can be helpful in engineering practice and in further studies on crushable granular materials via the discrete element method.
基金The study was financially supported by the National Natural Science Foundation of China,(Nos.51809292,51478481 and 51508141)Postdoctoral Fund of Central South University,China(No.205455)Beijing Municipal Science and Technology Project:Research and Application of Design and Construction Technology of Railway Engineering Traveling the Rift Valley,China(No.Z181100003918005).
文摘In discrete element method(DEM)simulations,multi-sphere(MS)clumped and convex particles are two main particle models that are used to study the mechanical behaviours of granular materials.Of interest is the evaluation of the effect of multiple contacts between clumped particles or single contacts between convex particles on the mechanical behaviours of granular materials.In this context,a series of drained triaxial compression tests were conducted on convex true(CT)ellipsoids and MS ellipsoids with aspect ratios(ARs)ranging from 1.0-2.0.The microscale results indicate that at a given AR,the critical friction angleφ_(c)changes with the particle type,whereas the peak friction angleφ_(p)is nearly independent of the particle type.The anisotropic analysis provides underlying mechanisms of the shear strength evolution from two perspectives.First,the anisotropies of granular materials are essential to shear strength as the deviatoric(q)-to-effective mean(p′)stress ratio can be expressed as the sum of the anisotropies,i.e.,q/p'≈0.4a_(c)+0.4a_(n)+0.6a_(t),where ac,an and at are the normal contact anisotropy,normal contact force anisotropy and tangential contact force anisotropy,respectively.For all samples,a_(c)and a_(n)underpin the shear strength and are influenced by the particle type.The similarφ_(p)displayed by the CT and MS ellipsoids does not translate to similar a_(n)and a_(c)but similar a_(c)+a_(n)for the two particle types.In addition,owing to their larger a_(c)+a_(n),the CT ellipsoids have a higherφ_(c)than the MS ellipsoids.Second,there is a satisfactory linear relationship between q/p'and ac within strong and non-sliding(sn)contacts a_(c)^(sn)(i.e.,q/p′=ka_(c)^(sn)),where k is the fitting parameter.Accordingly,in the peak state,the subtle difference in shear strength is attributed to the greater acsn in the CT ellipsoids than in the MS ellipsoids that is counteracted by the smaller k.However,in the critical state,the greater difference in a_(c)^(sn)between the CT and MS ellipsoids is partially offset by the smaller difference in k,causing a higherφ_(c)in the CT ellipsoids than in the MS ellipsoids.
基金financially supported by the National Natural Science Foundation of China,China(No.51809292,51478481 and 51508141)Postdoctoral Fund of Central South University,China(No.205455)Beijing Municipal Science and Technology Project:Research and Application of Design and Construction Technology of Railway Engineering Traveling the Rift Valley,China(No.Z181100003918005).
文摘Multi-sphere clumps are commonly used to simulate non-spherical particles in discrete element method simulations.It is of interest whether the degree of local non-convexity λ affects the mechanical behaviour of granular materials with the same non-convexity η.A series of discrete-element-method biaxial shear tests are conducted on rough particle packings with rη=0.075 and different λ values(ranging from 0.134 to 0.770).The microscale results show that the contact type changes with an increase in λ.However,the critical strength is independent of λ.The evaluation of the contributions of different contact types to the critical shear strength and a detailed analysis of the anisotropies help clarify the microscopic mechanisms that result in the independence of the critical shear strength from λ.
基金supported by the National Basic Research Program of China(973 Program)(Grant No.2014CB047100)the National Natural Science Foundation of China(Grant Nos.11572009,51538001 and 51609240)
文摘The soil-rock mixture(SRM) is highly heterogeneous. Before carrying out numerical analysis,a structure model should be generated. A reliable way to obtain such structure is by generating random aggregate structure based on random sequential addition(RSA). The classical RSA is neither efficient nor robust since valid positions to place new inclusions are formulated by trial, which involves repetitive overlapping tests. In this paper, the algorithm of Entrance block between block A and B(EAB)is synergized with background mesh to redesign RSA so that permissible positions to place new inclusions can be predicted,resulting in dramatic improvement in efficiency and robustness.
基金The work presented in this paper was supported by the National Natural Science Foundation of China （Grant Nos. 51308408, 41272291,51238009） and the Fundamental Research Funds for the Central Universities, and the Open Foundation of State Key Labo- ratory of Hydrology-Water Resources and Hydraulic Engineering （Grant No. 2014492311 ）.
文摘Discrete element modeling was used to investigate the effect of particle size distribution on the small strain shear stiffness of granular soils and explore the fundamental mechanism controlling this small strain shear stiffness at the particle level. The results indicate that the mean particle size has a negligible effect on the small strain shear modulus. The observed increase of the shear modulus with increasing particle size is caused by a scale effect. It is suggested that the ratio of sample size to the mean particle size should be larger than 11.5 to avoid this possible scale effect. At the same confining pressure and void ratio, the small strain shear modulus decreases as the coefficient of uniformity of the soil increases. The Poisson＇s ratio decreases with decreasing void ratio and increasing confining pressure instead of being constant as is commonly assumed. Microscopic analyses indicate that the small strain shear stiffness and Poisson＇s ratio depend uniquely on the soil＇s coordination number.
基金This research was supported by the National Natural Science Foundation of China(Nos.51708405,51378345)the Project of Tianjin Science and Technology Plan(No.16YDLJSF00040).
文摘The hydrogeologic conditions in Tianjin,China,are characterized by a system of multiple alternating aquifers and aquitards and a high groundwater level.Generally,the groundwater system to a depth of 70 m includes one phreatic aquifer and three semiconfined aquifers(AqI–AqIII).The distribution of aquifers varies greatly.Accurate estimates of aquifer parameters from pumping tests are essential in geotechnical engineering practice.In this study,the applicability of different models for estimating hydraulic parameters in a multiaquifer system is analyzed.To avoid errors in the graphic-analytical process and to enhance the efficiency and accuracy,a hybrid algorithm,called GALMA(a genetic algorithm(GA)combined with the Levenberg–Marquardt(LM)algorithm),is used with the Neuman and Witherspoon model and ratio method.The estimated values in the second and third semiconfined aquifers(AqII and AqIII)from various solutions are compared with the measured data.Further validations based on the prediction of the drawdown values are performed via a three dimensional numerical simulation.
基金This work was supported by the National Natural Science Foun-dation of China(grant numbers 41672262,41877227)the National Key Research and Development Program of China(grantnumber 2017YFC0806004).
文摘The critical solid fraction(ϕJ),which marks the transition between the solid and liquid phases in the jamming diagram,is influenced by several factors.In this study,the dependency ofϕJ on strain rate and boundary conditions is examined through discrete element method simulations considering a frictionless polydisperse granular system.Different approaches are used to determineϕJ.The observed boundary effect is due to the nonuniform solid fraction distribution induced by the clustering of particles close to rigid-wall boundaries at high compression rates.The solid fraction distribution within the sample in the rigid-wall simulations approaches that in the periodic-boundary simulations as the compression rate decreases.With increasing compression rate,the major force transmission network contains fewer mechanically stable particles and a less stable force transmission network.This causes jamming of the granular assembly at a lower solid fraction.These force transmission networks,however,are fragile and disintegrate quickly upon relaxation.
基金supported by the National Natural Science Foundation of China(Grant Nos.52478342 and 52238009)the Science and Technology Innovation Program of Hunan Province(Grant No.2022RC1172)the Natural Science Foundation of Jiangxi Province(Grant No.20223BBG71018),Finally,we deeply appreciate the warm and efficient work of editors and reviewers.
文摘Significant movement of in-situ retaining walls is usually assumed to begin with bulk excavation.However,an increasing number of case studies show that lowering the pore water pressures inside a diaphragm wall-type basement enclosure prior to bulk excavation can cause wall movements in the order of some centimeters.This paper describes the results of a laboratory-scale experiment carried out to explore mechanisms of in situ retaining wall movement associated with dewatering inside the enclosure prior to bulk excavation.Dewatering reduces the pore water pressures inside the enclosure more than outside,resulting in the wall moving as an unpropped cantilever supported only by the soil.Lateral effective stresses in the shallow soil behind the wall are reduced,while lateral effective stresses in front of the wall increase.Although the associated lateral movement was small in the laboratory experiment,the movement could be proportionately larger in the field with a less stiff soil and a potentially greater dewatered depth.The implementation of a staged dewatering system,coupled with the potential for phased excavation and propping strategies,can effectively mitigate dewatering-induced wall and soil movements.This approach allows for enhanced stiffness of the wall support system,which can be dynamically adjusted based on real-time displacement monitoring data when necessary.
基金support from the National Natural Science Foundation of China(Grant Nos.42377170 and 42407212)the Postdoctoral Fellowship Program of China Postdoctoral Science Foundation(Grant No.GZB20230606)+1 种基金the Postdoctoral Research Foundation of China(Grant No.2024M752679)the Sichuan Natural Science Foundation(Grant No.2025ZNSFSC1205).
文摘Underground debris flows,arising from the complex interplay of anthropogenic activities and rainfall-induced hydromechanical processes,present significant geotechnical hazards that remain poorly understood due to their hidden nature and dynamic multiphase triggers.Focusing on underground debris flow in a mining area in Southwest China,this study advances an integrated framework combining air-ground transient electromagnetic method(AGTEM)and computational fluid dynamics coupled with the discrete element method(CFD-DEM),revealing the migration mechanism in which microscale multiphase hydraulic erosion drives the macroscopic initiation of underground debris flow.Key findings include:(1)The identification of three transport phases(rapid erosion,slow erosion,and stabilization)provides actionable thresholds for monitoring and mitigation.(2)The coupled feedback between hydraulic conductivity anisotropy and the formation of preferential flow is the primary driver of large-scale debris transportation.(3)Linking mining-induced seismic energy to vibration-induced liquefaction via DEM simulations offers a physics-based explanation for flow mobilization triggers.The integrated geophysical-numerical framework offers new capabilities for predicting initiation thresholds and developing physics-based mitigation strategies in mining-affected terrains.
