Crack propagation in brittle material is not only crucial for structural safety evaluation,but also has a wideranging impact on material design,damage assessment,resource extraction,and scientific research.A thorough ...Crack propagation in brittle material is not only crucial for structural safety evaluation,but also has a wideranging impact on material design,damage assessment,resource extraction,and scientific research.A thorough investigation into the behavior of crack propagation contributes to a better understanding and control of the properties of brittle materials,thereby enhancing the reliability and safety of both materials and structures.As an implicit discrete elementmethod,the Discontinuous Deformation Analysis(DDA)has gained significant attention for its developments and applications in recent years.Among these developments,the particle DDA equipped with the bonded particle model is a powerful tool for predicting the whole process of material from continuity to failure.The primary objective of this research is to develop and utilize the particle DDAtomodel and understand the complex behavior of cracks in brittle materials under both static and dynamic loadings.The particle DDA is applied to several classical crack propagation problems,including the crack branching,compact tensile test,Kalthoff impact experiment,and tensile test of a rectangular plate with a hole.The evolutions of cracks under various stress or geometrical conditions are carefully investigated.The simulated results are compared with the experiments and other numerical results.It is found that the crack propagation patterns,including crack branching and the formation of secondary cracks,can be well reproduced.The results show that the particle DDA is a qualified method for crack propagation problems,providing valuable insights into the fracture mechanism of brittle materials.展开更多
Accurate dynamic modeling of landslides could help understand the movement mechanisms and guide disaster mitigation and prevention.Discontinuous deformation analysis(DDA)is an effective approach for investigating land...Accurate dynamic modeling of landslides could help understand the movement mechanisms and guide disaster mitigation and prevention.Discontinuous deformation analysis(DDA)is an effective approach for investigating landslides.However,DDA fails to accurately capture the degradation in shear strength of rock joints commonly observed in high-speed landslides.In this study,DDA is modified by incorporating simplified joint shear strength degradation.Based on the modified DDA,the kinematics of the Baige landslide that occurred along the Jinsha River in China on 10 October 2018 are reproduced.The violent starting velocity of the landslide is considered explicitly.Three cases with different violent starting velocities are investigated to show their effect on the landslide movement process.Subsequently,the landslide movement process and the final accumulation characteristics are analyzed from multiple perspectives.The results show that the violent starting velocity affects the landslide motion characteristics,which is found to be about 4 m/s in the Baige landslide.The movement process of the Baige landslide involves four stages:initiation,high-speed sliding,impact-climbing,low-speed motion and accumulation.The accumulation states of sliding masses in different zones are different,which essentially corresponds to reality.The research results suggest that the modified DDA is applicable to similar high-level rock landslides.展开更多
Since its introduction,discontinuous deformation analysis(DDA)has been widely used in different areas of rock mechanics.By dividing large blocks into subblocks and introducing artificial joints,DDA can be applied to r...Since its introduction,discontinuous deformation analysis(DDA)has been widely used in different areas of rock mechanics.By dividing large blocks into subblocks and introducing artificial joints,DDA can be applied to rock fracture simulation.However,parameter calibration,a fundamental issue in discontinuum methods,has not received enough attention in DDA.In this study,the parameter calibration of DDA for intact rock is carefully studied.To this end,a subblock DDA with Voronoi tessellation is presented first.Then,a modified contact constitutive law is introduced,in which the tensile and shear meso-strengths are modified to be independent of the bond lengths.This improvement can prevent the unjustified preferential failure of short edges.A method for imposing confining pressure is also introduced.Thereafter,sensitivity analysis is performed to investigate the influence of the calculated parameters and meso-parameters on the mechanical properties of modeled rock.Based on the sensitivity analysis,a unified calibration procedure is suggested for both cases with and without confining pressure.Finally,the calibration procedure is applied to two examples,including a biaxial compression test.The results show that the proposed Voronoi-based DDA can simulate rock fracture with and without confining pressure very well after careful parameter calibration.展开更多
The Niumiangou landslide was the largest landslide triggered by the 2008 Wenchuan earthquake,which was significantly affected by the amplification effect of seismic acceleration.The ringshear experiments indicated tha...The Niumiangou landslide was the largest landslide triggered by the 2008 Wenchuan earthquake,which was significantly affected by the amplification effect of seismic acceleration.The ringshear experiments indicated that the materials in the source area of the Niumiangou landslide were subjected to friction degradation under a big shear displacement,which may result in rapid movement of the landslide.In order to better understand the landslide movement and study the effect of the friction degradation on movement mechanisms,the dynamic process of Niumiangou landslide was simulated with a new numerical method,which combines the finite difference method(FDM)and the discontinuous deformation analysis(DDA).First,the FDM was used to study the initiation time,amplification effect and velocity of the landslide.Afterwards,these initiation velocities were applied to the blocks in the DDA model by corresponding coordination in the FDM model.A displacementdependent friction model of the sliding surface was incorporated into DDA code to further understand the kinetic behavior of the landslide.The results show that the displacement-dependent friction strongly decreases the friction coefficient of sliding surface under a big displacement,which can obviously promote the run-out and velocity of landslide.The model output well matches the topographic map formed by the landslide.This implies that the proposed model can be applied to the simulation of earthquake-induced landslides with amplification effect,and the friction degradation model is important to clarify the movement mechanism of high-speed and long-distance landslides.展开更多
The contact between blocks is treated by the open-close iteration in the conventional discontinuous deformation analysis (DDA),which needs to introduce spurious springs between two blocks in contact and to assume the ...The contact between blocks is treated by the open-close iteration in the conventional discontinuous deformation analysis (DDA),which needs to introduce spurious springs between two blocks in contact and to assume the normal stiffness and the tangential stiffness (the penalty factors). Unreasonable values of stiffness would result in numerical problems. To avoid the penalty factors and the open-close iteration,we reformulate the DDA as a mixed complementary problem (MiCP) and then choose the path Newton method (PNM) to solve the problem. Some examples including those originally designed by Shi are reanalyzed,which proves feasibility of the proposed procedure.展开更多
A discontinuous deformation and displacement(DDD) analysis method is proposed for modelling the rock failure process. This method combines the rock failure process analysis(RFPA) method(based on finite element method)...A discontinuous deformation and displacement(DDD) analysis method is proposed for modelling the rock failure process. This method combines the rock failure process analysis(RFPA) method(based on finite element method) and discontinuous deformation analysis(DDA) method. RFPA is used to simulate crack initiation, propagation and coalescence processes of rock during the small deformation state. The DDA method is used to simulate the movement of blocks created by the multiple cracks modelled by the RFPA. The newly developed DDD method is particularly suitable for modelling both crack propagation and block movement during the rock failure process because of the natural and convenient coupling of continuous and discontinuous deformation analyses. The proposed method has been used to simulate crack initiation, propagation and coalescence within a slope as well as the block movement during the landslide process. Numerical modelling results indicate that the proposed DDD method can automatically simulate crack propagation and block movement during the rock failure process without degrading accuracy.展开更多
This paper presents a modified method of discontinuous deformation analysis(DDA). In the presented method, open-close iteration may not be needed, small penetration is permitted among blocks, and springs are added bet...This paper presents a modified method of discontinuous deformation analysis(DDA). In the presented method, open-close iteration may not be needed, small penetration is permitted among blocks, and springs are added between contacting block pairs only when a penetration takes place. The three contact patterns(i.e. sliding, locking and opening) in original DDA method are not involved, and the recognition of these contact patterns and treatment of transformation among patterns are not required either,signi fi cantly saving the computing time. In a convex to concave contact, there are two candidate entrance edges which may cause uncertainty. In this case, we propose the angle bisector criterion to determine the entrance edge. The spring stiffness is much larger than Young's modulus in the original DDA, however we fi nd that the correct results can still be obtained when it is much smaller than Young's modulus. Finally,the penetrations by using penalty method and augmented Lagrangian method are compared. Penetration of the latter is 1/4 of the former. The range of spring stiffness for the latter is wider than the former,being 0.01-1 of the former. Both methods can lead to correct contact forces.展开更多
Rock blasting is a dynamic process accom panied with the propagations of shock waves and the dispersion of the explosion gas.This paper adopts the discontinuous deformation analysis(DDA)method to simulate the rock bla...Rock blasting is a dynamic process accom panied with the propagations of shock waves and the dispersion of the explosion gas.This paper adopts the discontinuous deformation analysis(DDA)method to simulate the rock blasting process.A dynamic parameter adjustment and the non-reflecting boundary condition are implemented in the DDA method.The sub-block DDA method to simulate fracture problems is used.The blasting process in jointed rock mass is simulated by application of the explosion gas pressure on the expanding borehole walls and induced connected fracture surfaces around the boreholes.The blast craters with different overburdens are derived.The whole process including the explosion gas dispersion,borehole expansion,rock mass failure and cast,and the formation of the final blasting piles in rock blasting are well reproduced numerically.Parametric study for different overburdens is carried out,and the results are analyzed and discussed.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.42372310).
文摘Crack propagation in brittle material is not only crucial for structural safety evaluation,but also has a wideranging impact on material design,damage assessment,resource extraction,and scientific research.A thorough investigation into the behavior of crack propagation contributes to a better understanding and control of the properties of brittle materials,thereby enhancing the reliability and safety of both materials and structures.As an implicit discrete elementmethod,the Discontinuous Deformation Analysis(DDA)has gained significant attention for its developments and applications in recent years.Among these developments,the particle DDA equipped with the bonded particle model is a powerful tool for predicting the whole process of material from continuity to failure.The primary objective of this research is to develop and utilize the particle DDAtomodel and understand the complex behavior of cracks in brittle materials under both static and dynamic loadings.The particle DDA is applied to several classical crack propagation problems,including the crack branching,compact tensile test,Kalthoff impact experiment,and tensile test of a rectangular plate with a hole.The evolutions of cracks under various stress or geometrical conditions are carefully investigated.The simulated results are compared with the experiments and other numerical results.It is found that the crack propagation patterns,including crack branching and the formation of secondary cracks,can be well reproduced.The results show that the particle DDA is a qualified method for crack propagation problems,providing valuable insights into the fracture mechanism of brittle materials.
基金supported by the National Natural Science Foundations of China(grant numbers U22A20601 and 52209142)the Opening fund of State Key Laboratory of Geohazard Prevention and Geoenvironment Protection(Chengdu University of Technology)(grant number SKLGP2022K018)+1 种基金the Science&Technology Department of Sichuan Province(grant number 2023NSFSC0284)the Science and Technology Major Project of Tibetan Autonomous Region of China(grant number XZ202201ZD0003G)。
文摘Accurate dynamic modeling of landslides could help understand the movement mechanisms and guide disaster mitigation and prevention.Discontinuous deformation analysis(DDA)is an effective approach for investigating landslides.However,DDA fails to accurately capture the degradation in shear strength of rock joints commonly observed in high-speed landslides.In this study,DDA is modified by incorporating simplified joint shear strength degradation.Based on the modified DDA,the kinematics of the Baige landslide that occurred along the Jinsha River in China on 10 October 2018 are reproduced.The violent starting velocity of the landslide is considered explicitly.Three cases with different violent starting velocities are investigated to show their effect on the landslide movement process.Subsequently,the landslide movement process and the final accumulation characteristics are analyzed from multiple perspectives.The results show that the violent starting velocity affects the landslide motion characteristics,which is found to be about 4 m/s in the Baige landslide.The movement process of the Baige landslide involves four stages:initiation,high-speed sliding,impact-climbing,low-speed motion and accumulation.The accumulation states of sliding masses in different zones are different,which essentially corresponds to reality.The research results suggest that the modified DDA is applicable to similar high-level rock landslides.
基金The authors would like to thank the National Natural Science Foundation of China(Grant Nos.51879184 and 52079091)for funding this work.
文摘Since its introduction,discontinuous deformation analysis(DDA)has been widely used in different areas of rock mechanics.By dividing large blocks into subblocks and introducing artificial joints,DDA can be applied to rock fracture simulation.However,parameter calibration,a fundamental issue in discontinuum methods,has not received enough attention in DDA.In this study,the parameter calibration of DDA for intact rock is carefully studied.To this end,a subblock DDA with Voronoi tessellation is presented first.Then,a modified contact constitutive law is introduced,in which the tensile and shear meso-strengths are modified to be independent of the bond lengths.This improvement can prevent the unjustified preferential failure of short edges.A method for imposing confining pressure is also introduced.Thereafter,sensitivity analysis is performed to investigate the influence of the calculated parameters and meso-parameters on the mechanical properties of modeled rock.Based on the sensitivity analysis,a unified calibration procedure is suggested for both cases with and without confining pressure.Finally,the calibration procedure is applied to two examples,including a biaxial compression test.The results show that the proposed Voronoi-based DDA can simulate rock fracture with and without confining pressure very well after careful parameter calibration.
基金supported by the National Natural Science Foundation of China(Nos.4167230041902290)Opening fund of State Key Laboratory of Geohazard Prevention and Geoenvironment Protection(Chengdu University of Technology)(No.SKLGP2017K015)。
文摘The Niumiangou landslide was the largest landslide triggered by the 2008 Wenchuan earthquake,which was significantly affected by the amplification effect of seismic acceleration.The ringshear experiments indicated that the materials in the source area of the Niumiangou landslide were subjected to friction degradation under a big shear displacement,which may result in rapid movement of the landslide.In order to better understand the landslide movement and study the effect of the friction degradation on movement mechanisms,the dynamic process of Niumiangou landslide was simulated with a new numerical method,which combines the finite difference method(FDM)and the discontinuous deformation analysis(DDA).First,the FDM was used to study the initiation time,amplification effect and velocity of the landslide.Afterwards,these initiation velocities were applied to the blocks in the DDA model by corresponding coordination in the FDM model.A displacementdependent friction model of the sliding surface was incorporated into DDA code to further understand the kinetic behavior of the landslide.The results show that the displacement-dependent friction strongly decreases the friction coefficient of sliding surface under a big displacement,which can obviously promote the run-out and velocity of landslide.The model output well matches the topographic map formed by the landslide.This implies that the proposed model can be applied to the simulation of earthquake-induced landslides with amplification effect,and the friction degradation model is important to clarify the movement mechanism of high-speed and long-distance landslides.
基金Supported by the National Natural Science Foundation of China (Grant No. 90715028)Research Project for Protection of High Cut Slopes in Three Gorges’ Reservoir (Grant No. 2008SXG01-2)
文摘The contact between blocks is treated by the open-close iteration in the conventional discontinuous deformation analysis (DDA),which needs to introduce spurious springs between two blocks in contact and to assume the normal stiffness and the tangential stiffness (the penalty factors). Unreasonable values of stiffness would result in numerical problems. To avoid the penalty factors and the open-close iteration,we reformulate the DDA as a mixed complementary problem (MiCP) and then choose the path Newton method (PNM) to solve the problem. Some examples including those originally designed by Shi are reanalyzed,which proves feasibility of the proposed procedure.
基金supported by the National Basic Research Program of China("973"Project)(Grant No.2014CB047100)the National Natural Science Foundation of China(Grant Nos.51421064,51474046 & 51174039)the Fundamental Research Funds for the Central Universities(Grant No.DUT14LK21)
文摘A discontinuous deformation and displacement(DDD) analysis method is proposed for modelling the rock failure process. This method combines the rock failure process analysis(RFPA) method(based on finite element method) and discontinuous deformation analysis(DDA) method. RFPA is used to simulate crack initiation, propagation and coalescence processes of rock during the small deformation state. The DDA method is used to simulate the movement of blocks created by the multiple cracks modelled by the RFPA. The newly developed DDD method is particularly suitable for modelling both crack propagation and block movement during the rock failure process because of the natural and convenient coupling of continuous and discontinuous deformation analyses. The proposed method has been used to simulate crack initiation, propagation and coalescence within a slope as well as the block movement during the landslide process. Numerical modelling results indicate that the proposed DDD method can automatically simulate crack propagation and block movement during the rock failure process without degrading accuracy.
基金supported by CRSRI Open Research Program (No.CKWV2014206/KY)the National Basic Research Program of China (No.2011CB710603)
文摘This paper presents a modified method of discontinuous deformation analysis(DDA). In the presented method, open-close iteration may not be needed, small penetration is permitted among blocks, and springs are added between contacting block pairs only when a penetration takes place. The three contact patterns(i.e. sliding, locking and opening) in original DDA method are not involved, and the recognition of these contact patterns and treatment of transformation among patterns are not required either,signi fi cantly saving the computing time. In a convex to concave contact, there are two candidate entrance edges which may cause uncertainty. In this case, we propose the angle bisector criterion to determine the entrance edge. The spring stiffness is much larger than Young's modulus in the original DDA, however we fi nd that the correct results can still be obtained when it is much smaller than Young's modulus. Finally,the penetrations by using penalty method and augmented Lagrangian method are compared. Penetration of the latter is 1/4 of the former. The range of spring stiffness for the latter is wider than the former,being 0.01-1 of the former. Both methods can lead to correct contact forces.
文摘Rock blasting is a dynamic process accom panied with the propagations of shock waves and the dispersion of the explosion gas.This paper adopts the discontinuous deformation analysis(DDA)method to simulate the rock blasting process.A dynamic parameter adjustment and the non-reflecting boundary condition are implemented in the DDA method.The sub-block DDA method to simulate fracture problems is used.The blasting process in jointed rock mass is simulated by application of the explosion gas pressure on the expanding borehole walls and induced connected fracture surfaces around the boreholes.The blast craters with different overburdens are derived.The whole process including the explosion gas dispersion,borehole expansion,rock mass failure and cast,and the formation of the final blasting piles in rock blasting are well reproduced numerically.Parametric study for different overburdens is carried out,and the results are analyzed and discussed.
