This study investigated the failure mechanism associated with the rock mass structure and the dynamic fragmentation process of blocky rocks of the 2018 Daanshan rockslide that occurred on 11 August,2018.It was found t...This study investigated the failure mechanism associated with the rock mass structure and the dynamic fragmentation process of blocky rocks of the 2018 Daanshan rockslide that occurred on 11 August,2018.It was found that the initially collapsed rock of this rockslide was partitioned along the unconformity and strata interfaces.We analyzed how the unique rock mass structure,coupled with the road cut and the antecedent rainfall,jointly resulted in its failure.Based on the rock types and geological structures,the initial stratified configuration of a discrete element model was setup to reveal the influences of the local structure.The numerical model was divided into three parts.Part 1 is the basalt of the Nandaling Formation,the normal and shear stiffnesses of the basalt particles are set as 80 MPa and 40 MPa.Parts 2 and 3 are the sandstones interbedded with mudstone and sandstone of the Shihezi Formation,and the normal and shear stiffnesses of these parts were set as 6 MPa and 10 MPa,respectively.The dynamic process of the rockslide,particularly the rock fragmentation process,was numerically analyzed using a 3D discrete element method.The numerical results were compared with real-time videos and field investigations.The results show that the rock fragmentation and the final deposition range match well with the real disaster phenomenon,and the calculation accuracy of the rockslide reaches 82.41%.Moreover,a parameter sensitivity analysis was conducted,and classical uniform models under different bonding forces were established;the stratified model can better restore the true state of the fragmentation,movement,and deposition processes of rockslides.Therefore,for complicated rocks with significant differences in lithology,clarifying the rock mass stratigraphy is essential for an accurate reconstruction of the dynamic process of rockslides.展开更多
A catastrophic landslide occurred at Xinmo village in Maoxian County, Sichuan Province,China, on June 24, 2017. A 2.87×106 m3 rock mass collapsed and entrained the surface soil layer along the landslide path. Eig...A catastrophic landslide occurred at Xinmo village in Maoxian County, Sichuan Province,China, on June 24, 2017. A 2.87×106 m3 rock mass collapsed and entrained the surface soil layer along the landslide path. Eighty-three people were killed or went missing and more than 103 houses were destroyed. In this paper, the geological conditions of the landslide are analyzed via field investigation and high-resolution imagery. The dynamic process and runout characteristics of the landslide are numerically analyzed using a depth-integrated continuum method and Mac Cormack-TVD finite difference algorithm.Computational results show that the evaluated area of the danger zone matchs well with the results of field investigation. It is worth noting that soil sprayed by the high-speed blast needs to be taken into account for such kind of large high-locality landslide. The maximum velocity is about 55 m/s, which is consistent with most cases. In addition, the potential danger zone of an unstable block is evaluated. The potential risk area evaluated by the efficient depthintegrated continuum method could play a significant role in disaster prevention and secondary hazard avoidance during rescue operations.展开更多
Many rock avalanches were triggered by the Wenchuan earthquake on May 12,2008 in southwest China.Protection galleries covered with a single soil layer are usually used to protect against rockfall.Since one-layer prote...Many rock avalanches were triggered by the Wenchuan earthquake on May 12,2008 in southwest China.Protection galleries covered with a single soil layer are usually used to protect against rockfall.Since one-layer protection galleries do not have sufficient buffer capacity,a two-layered absorbing system has been designed.This study aims to find whether an expanded poly-styrol(EPS) cushion,which is used in the soil-covered protection galleries for shock absorption,could be positioned under dynamic loadings.The dynamic impacts of the two-layered absorbing system under the conditions ofrock avalanches are numerically simulated through a 2D discrete element method.By selecting reasonable parameters,a series of numerical experiments were conducted to find the best combination for the twolayered absorbing system.The values of the EPS layer area as a percentage of the total area were set as 0%(S1),22%(S2),and 70%(S3).22% of the area of the EPS layer was found to be a reasonable value,and experiments were conducted to find the best position of the EPS layer in the two-layered absorbing system.The numerical results yield useful conclusions regarding the interaction between the impacting avalanches and the two-layered absorbing system.The soil layer can absorb the shock energy effectively and S2(0.4-m thick EPS cushion covered with soillayer) is the most efficient combination,which can reduce the impact force,compared with the other combinations.展开更多
In this paper, a generalized limit equilibrium method of solving the active earth pressure problem behind a retaining wall is proposed.Differing from other limit equilibrium methods, an arbitrary slip surface shape wi...In this paper, a generalized limit equilibrium method of solving the active earth pressure problem behind a retaining wall is proposed.Differing from other limit equilibrium methods, an arbitrary slip surface shape without any assumptions of pre-defined shapes is needed in the current framework, which is verified to find the most probable failure slip surface. Based on the current computational framework, numerical comparisons with experiment, discrete element method and other methods are carried out. In addition, the influences of the inclination of the wall, the soil cohesion, the angle of the internal friction of the soil, the slope inclination of the backfill soil on the critical pressure coefficient of the soil, the point of application of the resultant earth pressure and the shape of the slip surface are also carefully investigated. The results demonstrate that limit equilibrium solution from predefined slip plane assumption, including Coulomb solution, is a special case of current computational framework. It is well illustrated that the current method is feasible to evaluate the characteristics of earth pressure problem.展开更多
A catastrophic landslide occurred at Hongao dumpsite in Guangming New District of Shenzhen, South China, on December 20, 2015. An estimated total volume of 2.73×106 m3 of construction spoils was mobilized during ...A catastrophic landslide occurred at Hongao dumpsite in Guangming New District of Shenzhen, South China, on December 20, 2015. An estimated total volume of 2.73×106 m3 of construction spoils was mobilized during this event. The landslide traveled a long distance on a low-relief terrain. The affected area was approximately 1100 m in length and 630 m in width. This landslide made 33 buildings destroyed, 73 people died and 4 people lost. Due to the special dumping history and other factors, soil in this landfill is of high initial water content. To identify the major factors that attribute to the long runout character, a two-phase flow model of Iverson and George was used to simulate the dynamics of this landslide. The influence of initial hydraulic permeability, initial dilatancy, and earth pressure coefficient was examined through numerical simulations. We found that pore pressure has the most significant effect on the dynamic characteristics of Shenzhen landslides. Average pore pressure ratio ofthe whole basal surface was used to evaluate the degree of liquefaction for the sliding material. The evolution and influence factors of this ratio were analyzed based on the computational results. An exponential function was proposed to fit the evolution curve of the average pore pressure ratio, which can be used as a reasonable and simplified evaluation of the pore pressure. This fitting function can be utilized to improve the single-phase flow model.展开更多
The great diversity and complexity of geological hazards in terms of flowing materials,environment,triggering mechanisms and physical processes during the flow bring great difficulties to the numerical parameter selec...The great diversity and complexity of geological hazards in terms of flowing materials,environment,triggering mechanisms and physical processes during the flow bring great difficulties to the numerical parameter selection for the discrete element method.In order to identity the significance of individual parameters on the landslides dynamic process and provide valuable contribution to the runout analysis of similar landslide,the dynamic process and associated microscopic mechanism of the Turnoff Creek rock avalanche in Canada are simulated.The present numerical results are compared with the field survey data and the results of depth-integrated continuum method.The final deposit range matches well with the field survey data.It is illustrated that the discrete element method is robust and feasible to capture the dynamic characteristics of large rock avalanche over a complex terrain.Besides,a new method to assess the landslide hazard level based on the discrete element method is proposed.According to the parameter sensitivity analysis,it is demonstrated that the basal friction coefficient and bond strength are essential to the final deposit while rolling coefficient and restitution coefficient have little effects on it.展开更多
基金funded by the Strategic Priority Research Program of CAS(Grant No.XDA23090303)the NSFC(Grant No.42022054)+1 种基金Sichuan Science and Technology Program(Grant No.2022YFS0543)the Science Foundation for Distinguished Young Scholars of Sichuan Province(Grant No.2020JDJQ0044)State Key Laboratory of Geohazard Prevention and Geoenvironment Protection Independent Research Project(SKLGP2019Z013).
文摘This study investigated the failure mechanism associated with the rock mass structure and the dynamic fragmentation process of blocky rocks of the 2018 Daanshan rockslide that occurred on 11 August,2018.It was found that the initially collapsed rock of this rockslide was partitioned along the unconformity and strata interfaces.We analyzed how the unique rock mass structure,coupled with the road cut and the antecedent rainfall,jointly resulted in its failure.Based on the rock types and geological structures,the initial stratified configuration of a discrete element model was setup to reveal the influences of the local structure.The numerical model was divided into three parts.Part 1 is the basalt of the Nandaling Formation,the normal and shear stiffnesses of the basalt particles are set as 80 MPa and 40 MPa.Parts 2 and 3 are the sandstones interbedded with mudstone and sandstone of the Shihezi Formation,and the normal and shear stiffnesses of these parts were set as 6 MPa and 10 MPa,respectively.The dynamic process of the rockslide,particularly the rock fragmentation process,was numerically analyzed using a 3D discrete element method.The numerical results were compared with real-time videos and field investigations.The results show that the rock fragmentation and the final deposition range match well with the real disaster phenomenon,and the calculation accuracy of the rockslide reaches 82.41%.Moreover,a parameter sensitivity analysis was conducted,and classical uniform models under different bonding forces were established;the stratified model can better restore the true state of the fragmentation,movement,and deposition processes of rockslides.Therefore,for complicated rocks with significant differences in lithology,clarifying the rock mass stratigraphy is essential for an accurate reconstruction of the dynamic process of rockslides.
基金Financial support from National Nature Science Foundation of China (Grant No. 41572303, 41520104002)Chinese Academy of Sciences “Light of West China” Program and Youth Innovation Promotion Association
文摘A catastrophic landslide occurred at Xinmo village in Maoxian County, Sichuan Province,China, on June 24, 2017. A 2.87×106 m3 rock mass collapsed and entrained the surface soil layer along the landslide path. Eighty-three people were killed or went missing and more than 103 houses were destroyed. In this paper, the geological conditions of the landslide are analyzed via field investigation and high-resolution imagery. The dynamic process and runout characteristics of the landslide are numerically analyzed using a depth-integrated continuum method and Mac Cormack-TVD finite difference algorithm.Computational results show that the evaluated area of the danger zone matchs well with the results of field investigation. It is worth noting that soil sprayed by the high-speed blast needs to be taken into account for such kind of large high-locality landslide. The maximum velocity is about 55 m/s, which is consistent with most cases. In addition, the potential danger zone of an unstable block is evaluated. The potential risk area evaluated by the efficient depthintegrated continuum method could play a significant role in disaster prevention and secondary hazard avoidance during rescue operations.
基金financial support from the Project of National Science Foundation of China(Grant No.41272346)the National Outstanding Youth Funds(Grant No.41225011)+2 种基金financial support from the Science & Technology Research Plan of China Railway Eryuan Engineering Group CO.LTD (Grant No.13164196(13-15))the Project of National Science Foundation of China(Grant Nos. 41472293,91430105)"hundred talents" program of CAS
文摘Many rock avalanches were triggered by the Wenchuan earthquake on May 12,2008 in southwest China.Protection galleries covered with a single soil layer are usually used to protect against rockfall.Since one-layer protection galleries do not have sufficient buffer capacity,a two-layered absorbing system has been designed.This study aims to find whether an expanded poly-styrol(EPS) cushion,which is used in the soil-covered protection galleries for shock absorption,could be positioned under dynamic loadings.The dynamic impacts of the two-layered absorbing system under the conditions ofrock avalanches are numerically simulated through a 2D discrete element method.By selecting reasonable parameters,a series of numerical experiments were conducted to find the best combination for the twolayered absorbing system.The values of the EPS layer area as a percentage of the total area were set as 0%(S1),22%(S2),and 70%(S3).22% of the area of the EPS layer was found to be a reasonable value,and experiments were conducted to find the best position of the EPS layer in the two-layered absorbing system.The numerical results yield useful conclusions regarding the interaction between the impacting avalanches and the two-layered absorbing system.The soil layer can absorb the shock energy effectively and S2(0.4-m thick EPS cushion covered with soillayer) is the most efficient combination,which can reduce the impact force,compared with the other combinations.
基金Financial support from the Key Research Program of Chinese Academy of Sciences (Grant No. KZZD-EW-05-01)the NSFC (Grant Nos. 41101008, 41272346)the Youth Talent Team Program of the Institute of Mountain Hazards and Environment, Chinese Academy of Sciences, Opening Fund of State Key Laboratory of Geohazard Prevention and Geoenvironment Protection (Chengdu University of Technology) (Grant No. SKLGP2011K010)
文摘In this paper, a generalized limit equilibrium method of solving the active earth pressure problem behind a retaining wall is proposed.Differing from other limit equilibrium methods, an arbitrary slip surface shape without any assumptions of pre-defined shapes is needed in the current framework, which is verified to find the most probable failure slip surface. Based on the current computational framework, numerical comparisons with experiment, discrete element method and other methods are carried out. In addition, the influences of the inclination of the wall, the soil cohesion, the angle of the internal friction of the soil, the slope inclination of the backfill soil on the critical pressure coefficient of the soil, the point of application of the resultant earth pressure and the shape of the slip surface are also carefully investigated. The results demonstrate that limit equilibrium solution from predefined slip plane assumption, including Coulomb solution, is a special case of current computational framework. It is well illustrated that the current method is feasible to evaluate the characteristics of earth pressure problem.
基金supported by the National Key R&D Program of China(Grant Nos.2017YFC1502502,2017YFC1502506)National Nature Science Foundation of China(Grant Nos.41672318,51679229,41372331)+1 种基金135 Strategic Program of the Institute of Mountain Hazards and Environment,CAS(Grant No.SDS-135-1701)supported by Youth Innovation Promotion Association of the Chinese Academy of Sciences(2018405)
文摘A catastrophic landslide occurred at Hongao dumpsite in Guangming New District of Shenzhen, South China, on December 20, 2015. An estimated total volume of 2.73×106 m3 of construction spoils was mobilized during this event. The landslide traveled a long distance on a low-relief terrain. The affected area was approximately 1100 m in length and 630 m in width. This landslide made 33 buildings destroyed, 73 people died and 4 people lost. Due to the special dumping history and other factors, soil in this landfill is of high initial water content. To identify the major factors that attribute to the long runout character, a two-phase flow model of Iverson and George was used to simulate the dynamics of this landslide. The influence of initial hydraulic permeability, initial dilatancy, and earth pressure coefficient was examined through numerical simulations. We found that pore pressure has the most significant effect on the dynamic characteristics of Shenzhen landslides. Average pore pressure ratio ofthe whole basal surface was used to evaluate the degree of liquefaction for the sliding material. The evolution and influence factors of this ratio were analyzed based on the computational results. An exponential function was proposed to fit the evolution curve of the average pore pressure ratio, which can be used as a reasonable and simplified evaluation of the pore pressure. This fitting function can be utilized to improve the single-phase flow model.
基金Financial support from the National Natural Science Foundation of China(Grant No.41520104002,41572303)the Strategic Priority Research Program of CAS(Grant No.XDA23090303)the National Key Research and Development Program of China(Project No.2017YFC1501000)。
文摘The great diversity and complexity of geological hazards in terms of flowing materials,environment,triggering mechanisms and physical processes during the flow bring great difficulties to the numerical parameter selection for the discrete element method.In order to identity the significance of individual parameters on the landslides dynamic process and provide valuable contribution to the runout analysis of similar landslide,the dynamic process and associated microscopic mechanism of the Turnoff Creek rock avalanche in Canada are simulated.The present numerical results are compared with the field survey data and the results of depth-integrated continuum method.The final deposit range matches well with the field survey data.It is illustrated that the discrete element method is robust and feasible to capture the dynamic characteristics of large rock avalanche over a complex terrain.Besides,a new method to assess the landslide hazard level based on the discrete element method is proposed.According to the parameter sensitivity analysis,it is demonstrated that the basal friction coefficient and bond strength are essential to the final deposit while rolling coefficient and restitution coefficient have little effects on it.