In the framework of finite element meshes,a novel continuous/discontinuous deformation analysis(CDDA)method is proposed in this paper for modeling of crack problems.In the present CDDA,simple polynomial interpolations...In the framework of finite element meshes,a novel continuous/discontinuous deformation analysis(CDDA)method is proposed in this paper for modeling of crack problems.In the present CDDA,simple polynomial interpolations are defined at the deformable block elements,and a link element is employed to connect the adjacent block elements.The CDDA is particularly suitable for modeling the fracture propagation because the switch from continuous deformation analysis to discontinuous deformation analysis is natural and convenient without additional procedures.The SIFs(stress intensity factors)for various types of cracks,such as kinked cracks or curved cracks,can be easily computed in the CDDA by using the virtual crack extension technique(VCET).Both the formulation and implementation of the VCET in CDDA are simple and straightforward.Numerical examples indicate that the present CDDA can obtain high accuracy in SIF results with simple polynomial interpolations and insensitive to mesh sizes,and can automatically simulate the crack propagation without degrading accuracy.展开更多
A method of continuous-discontinuous cellular automaton for modeling the growth and coalescence of multiple cracks in brittle material is presented. The method uses the level set to track arbitrary discontinuities, an...A method of continuous-discontinuous cellular automaton for modeling the growth and coalescence of multiple cracks in brittle material is presented. The method uses the level set to track arbitrary discontinuities, and calculation grids are independent of the discontinuities and no remeshing are required with the crack growing. Based on Grif- fith fracture theory and Mohr-Coulumb criterion, a mixed fracture criterion for multiple cracks growth in brittle mate- rial is proposed. The method treats the junction and coales- cence of multiple cracks, and junction criterion and coales- cence criterion for brittle material are given, too. Besides, in order to overcome the tracking error in the level set ap- proximation for crack junction and coalescence, a dichotomy searching algorithm is proposed. Introduced the above the- ories into continuous-discontinuous cellular automaton, the present method can be applied to solving multiple crack growth in brittle material, and only cell stiffness is needed and no assembled global stiffness is needed. Some numerical examples are given to shown that the present method is efficient and accurate for crack junction, coalescence and percolation problems.展开更多
A plane mechanical model of rock breaking process by double disc cutter at the center of the cutterhead is established based on contact mechanics to analyze the stress evolution in the rock broken by cutters with diff...A plane mechanical model of rock breaking process by double disc cutter at the center of the cutterhead is established based on contact mechanics to analyze the stress evolution in the rock broken by cutters with different spacings. A continuous-discontinuous coupling numerical method based on zero-thickness cohesive elements is developed to simulate rock breaking using double cutters. The process, mechanism,and characteristics of rock breaking are comprehensively analyzed from five aspects: peak force, breaking form, breaking efficiency, crack mode, and breaking degree. The results show that under the penetrating action of cutters, dense cores are formed due to shear failure under respective cutters. The tensile cracks propagate in the rock, and then rock chips form with increasing penetration depth. When the cutter spacing is increased from 10 to 80 mm, the peak force gradually increases, the rock breaking range increases first and then decreases, the specific energy decreases first and then rises, and the breaking coefficient of intermediate rock decreases from 0.955 to 0.788. The area of rock breaking is positively correlated with the length of the tensile crack. Furthermore, the length of the tensile crack accounts for 14.4%–33.6% of the total crack length.展开更多
Thermal cracking of rocks can significantly affect the durability of underground structures in engineering practices such as geothermal energy extraction,storage of nuclear waste and tunnelling in freezeethaw cycle in...Thermal cracking of rocks can significantly affect the durability of underground structures in engineering practices such as geothermal energy extraction,storage of nuclear waste and tunnelling in freezeethaw cycle induced areas.It is a scenario of strong coupled thermomechanical process involving discontinuity behaviours of rocks.In this context,a numerical model was proposed to investigate the thermal cracking of rocks,in a framework of the continuous-discontinuous element method(CDEM)for efficiently capturing the initiation and propagation of multiple cracks.A simplex integration strategy was adopted to account for the influences of temperature-dependent material properties.Several benchmark tests were considered and the obtained results were compared with analytical solutions and numerical results from the literature.The results show that the fracture degree of the cases when considering temperature-dependent material parameters had 10%differences approximately compared with the cases with constant parameters.展开更多
The numerical simulation of internal features,such as inclusions and voids,is important to analyze their impact on the performance of composite materials.However,the complex geometries of internal features and the ind...The numerical simulation of internal features,such as inclusions and voids,is important to analyze their impact on the performance of composite materials.However,the complex geometries of internal features and the induced continuous-discontinuous(C-D)deformation fields are challenges to their numerical simulation.In this study,a 3D approach using a simple mesh to simulate irregular internal geometries is developed for the first time.With the help of a developed voxel crack model,image models that are efficient when recording complex geometries are directly imported into the simulation.Surface reconstructions,which are usually labor-intensive,are excluded from this approach.Moreover,using image models as the geometric input,image processing techniques are applied to detect material interfaces and develop contact pairs.Then,the C-D deformations of the complex internal features are directly calculated based on the numerical manifold method.The accuracy and convergence of the developed3D approach are examined based on multiple benchmarks.Successful 3D C-D simulation of sandstones with naturally formed complex microfeatures demonstrates the capability of the developed approach.展开更多
A continuous-discontinuous cellular automaton method is developed for rock initiation and propagation simulations, in which the level set method, discontinuous enrichment shape functions and discontinuous cellular aut...A continuous-discontinuous cellular automaton method is developed for rock initiation and propagation simulations, in which the level set method, discontinuous enrichment shape functions and discontinuous cellular automaton are combined. No renmshing is needed for crack growth analysis, and all calculations are restricted to cells without an assembled global stiffness matrix. The frictional contact theory is employed to construct the contact model of normal pressure and tangential shear on crack surfaces. A discontinuous cellular automaton updating rule suitable for frictional contact of rock is proposed simultaneously with Newton＇s iteration method for nonlinear iteration. Besides, a comprehensive fracturing criterion for brittle rock under compression-shear loading is developed. The accuracy and effectivenesss of the proposed method is proved by numerical simulation.展开更多
基金The authors gratefully acknowledge the support of Nature Science Foundation of China(Grant No.41130751)National Basic Research Program of China(Grant No.2011CB013800)New Century Excellent Talents Project in China(NCET-12-0415).
文摘In the framework of finite element meshes,a novel continuous/discontinuous deformation analysis(CDDA)method is proposed in this paper for modeling of crack problems.In the present CDDA,simple polynomial interpolations are defined at the deformable block elements,and a link element is employed to connect the adjacent block elements.The CDDA is particularly suitable for modeling the fracture propagation because the switch from continuous deformation analysis to discontinuous deformation analysis is natural and convenient without additional procedures.The SIFs(stress intensity factors)for various types of cracks,such as kinked cracks or curved cracks,can be easily computed in the CDDA by using the virtual crack extension technique(VCET).Both the formulation and implementation of the VCET in CDDA are simple and straightforward.Numerical examples indicate that the present CDDA can obtain high accuracy in SIF results with simple polynomial interpolations and insensitive to mesh sizes,and can automatically simulate the crack propagation without degrading accuracy.
基金supported by the National Key Basic Research Program of China(2013CB036405)the National Natural Science Foundation of China(11002154,41272349,and 41372315)the CAS/SAFEA International Partnership Program for Creative Research Teams(KZCX2-YW-T12)
文摘A method of continuous-discontinuous cellular automaton for modeling the growth and coalescence of multiple cracks in brittle material is presented. The method uses the level set to track arbitrary discontinuities, and calculation grids are independent of the discontinuities and no remeshing are required with the crack growing. Based on Grif- fith fracture theory and Mohr-Coulumb criterion, a mixed fracture criterion for multiple cracks growth in brittle mate- rial is proposed. The method treats the junction and coales- cence of multiple cracks, and junction criterion and coales- cence criterion for brittle material are given, too. Besides, in order to overcome the tracking error in the level set ap- proximation for crack junction and coalescence, a dichotomy searching algorithm is proposed. Introduced the above the- ories into continuous-discontinuous cellular automaton, the present method can be applied to solving multiple crack growth in brittle material, and only cell stiffness is needed and no assembled global stiffness is needed. Some numerical examples are given to shown that the present method is efficient and accurate for crack junction, coalescence and percolation problems.
基金funded by the National Key Research and Development Program of China (No. 2021YFB3401501)the Fundamental Research Funds for the Central Universities (No2022JCCXLJ01)。
文摘A plane mechanical model of rock breaking process by double disc cutter at the center of the cutterhead is established based on contact mechanics to analyze the stress evolution in the rock broken by cutters with different spacings. A continuous-discontinuous coupling numerical method based on zero-thickness cohesive elements is developed to simulate rock breaking using double cutters. The process, mechanism,and characteristics of rock breaking are comprehensively analyzed from five aspects: peak force, breaking form, breaking efficiency, crack mode, and breaking degree. The results show that under the penetrating action of cutters, dense cores are formed due to shear failure under respective cutters. The tensile cracks propagate in the rock, and then rock chips form with increasing penetration depth. When the cutter spacing is increased from 10 to 80 mm, the peak force gradually increases, the rock breaking range increases first and then decreases, the specific energy decreases first and then rises, and the breaking coefficient of intermediate rock decreases from 0.955 to 0.788. The area of rock breaking is positively correlated with the length of the tensile crack. Furthermore, the length of the tensile crack accounts for 14.4%–33.6% of the total crack length.
基金the financial support from the Natural Science Foundation of Hebei Province(Grant No.E2020050012)the National Natural Science Foundation of China(NSFC)(Grant No.52178324)the National Key Research and Development Project of China,the Ministry of Science and Technology of China(Grant No.2018YFC1505504).
文摘Thermal cracking of rocks can significantly affect the durability of underground structures in engineering practices such as geothermal energy extraction,storage of nuclear waste and tunnelling in freezeethaw cycle induced areas.It is a scenario of strong coupled thermomechanical process involving discontinuity behaviours of rocks.In this context,a numerical model was proposed to investigate the thermal cracking of rocks,in a framework of the continuous-discontinuous element method(CDEM)for efficiently capturing the initiation and propagation of multiple cracks.A simplex integration strategy was adopted to account for the influences of temperature-dependent material properties.Several benchmark tests were considered and the obtained results were compared with analytical solutions and numerical results from the literature.The results show that the fracture degree of the cases when considering temperature-dependent material parameters had 10%differences approximately compared with the cases with constant parameters.
基金supported by the National Natural Science Foundation of China(Grant Nos.41807277,42172306,and U1965204)the Natural Science Foundation of Hebei Province(Grant No.D2019202440)。
文摘The numerical simulation of internal features,such as inclusions and voids,is important to analyze their impact on the performance of composite materials.However,the complex geometries of internal features and the induced continuous-discontinuous(C-D)deformation fields are challenges to their numerical simulation.In this study,a 3D approach using a simple mesh to simulate irregular internal geometries is developed for the first time.With the help of a developed voxel crack model,image models that are efficient when recording complex geometries are directly imported into the simulation.Surface reconstructions,which are usually labor-intensive,are excluded from this approach.Moreover,using image models as the geometric input,image processing techniques are applied to detect material interfaces and develop contact pairs.Then,the C-D deformations of the complex internal features are directly calculated based on the numerical manifold method.The accuracy and convergence of the developed3D approach are examined based on multiple benchmarks.Successful 3D C-D simulation of sandstones with naturally formed complex microfeatures demonstrates the capability of the developed approach.
基金supported by the National Key Technologies R&D Program of China(No.2013BAB02B01)the National Natural Science Foundation of China(Nos.41272349,41172284 and 51322906)
文摘A continuous-discontinuous cellular automaton method is developed for rock initiation and propagation simulations, in which the level set method, discontinuous enrichment shape functions and discontinuous cellular automaton are combined. No renmshing is needed for crack growth analysis, and all calculations are restricted to cells without an assembled global stiffness matrix. The frictional contact theory is employed to construct the contact model of normal pressure and tangential shear on crack surfaces. A discontinuous cellular automaton updating rule suitable for frictional contact of rock is proposed simultaneously with Newton＇s iteration method for nonlinear iteration. Besides, a comprehensive fracturing criterion for brittle rock under compression-shear loading is developed. The accuracy and effectivenesss of the proposed method is proved by numerical simulation.
