考虑多晶材料疲劳行为受晶粒拓扑结构的影响,建立了多晶材料Voronoi镶嵌的有限元微观结构模型,并在Voronoi边界插入内聚力单元模拟晶界特性。在此基础上,基于Abaqus双线性内聚力本构关系,建立了内聚力单元的疲劳累积损伤模型,定义了损...考虑多晶材料疲劳行为受晶粒拓扑结构的影响,建立了多晶材料Voronoi镶嵌的有限元微观结构模型,并在Voronoi边界插入内聚力单元模拟晶界特性。在此基础上,基于Abaqus双线性内聚力本构关系,建立了内聚力单元的疲劳累积损伤模型,定义了损伤变量,损伤判据以及损伤规律;通过Abaqus UMAT(User-defined material mechanical behavior)子程序编写了内聚力单元疲劳损伤的数值计算程序,实现了对微观结构下材料滚动接触疲劳裂纹萌生与扩展的行为模拟。以滚子与滚道的接触疲劳分析为例,将数值模拟结果与实测结果相比,表明了该数值仿真方法模拟多晶材料裂纹形成和扩展的有效性。展开更多
The effect of freeze-thaw(F-T)cycles on the mechanical behaviors and internal mechanism of rock mass is a critical research topic.In permafrost or seasonally frozen regions,F-T cycles have adverse effects on the mecha...The effect of freeze-thaw(F-T)cycles on the mechanical behaviors and internal mechanism of rock mass is a critical research topic.In permafrost or seasonally frozen regions,F-T cycles have adverse effects on the mechanical properties of rock mass,leading to many serious disasters in mining and geotechnical operations.In this paper,uniaxial compression tests are carried out on cyan sandstone after different F-T cycles.The failure modes and damage evolution of cyan sandstone under F-T cycles are studied.In addition,from the perspective of fracture and pore volume,the calculation equations of rock strain under frost heaving pressure and F-T cycles are established and verified with the corresponding laboratory tests.Subsequently,based on the classical damage theory,the F-T damage variables of cyan sandstone under different F-T cycles are calculated,and the meso-damage calculation model of cyan sandstone under F-T-loading coupling conditions is derived.Furthermore,through the discrete element numerical simulation software(PFC^(3D)),the microscopic damage evolution process of cyan sandstone under uniaxial compression after F-T cycles is studied,including the change of microcracks number,distribution of microcracks,and the acoustic emission(AE)count.The goal of this study is to investigate the damage evolution mechanism of rock from the mesoscopic and microscopic aspects,which has certain guiding value for accurately understanding the damage characteristics of rock in cold regions.展开更多
By considering the effect of hydraulic pressure filled in wing crack and the connected part of main crack on the stress intensity factor at wing crack tip, a new wing crack model exerted by hydraulic pressure and far ...By considering the effect of hydraulic pressure filled in wing crack and the connected part of main crack on the stress intensity factor at wing crack tip, a new wing crack model exerted by hydraulic pressure and far field stresses was proposed. By introducing the equivalent crack length lcq of wing crack, two terms make up the stress intensity factor K1 at wing crack tip： one is the component K（1） for a single isolated straight wing crack of length 2l subjected to hydraulic pressure in wing crack and far field stresses, and the other is the component K1^（2） due to the effective shear stress induced by the presence of the equivalent main crack. The FEM model of wing crack propagation subjected to hydraulic pressure and far field stresses was also established according to different side pressure coefficients and hydraulic pressures in crack. The result shows that a good agreement is found between theoretical model of wing crack proposed and finite element method （FEM）. In theory, an unstable crack propagation is shown if there is high hydraulic pressure and lateral tension. The wing crack model proposed can provide references for studying on hydraulic fracturing in rock masses.展开更多
The present research is focused on the numerical crack coalescence analysis of the micro-cracks and cracks produced during the cutting action of TBM disc cutters. The linear elastic fracture mechanics(LEFM) concepts a...The present research is focused on the numerical crack coalescence analysis of the micro-cracks and cracks produced during the cutting action of TBM disc cutters. The linear elastic fracture mechanics(LEFM) concepts and the maximum tangential stress criterion are used to investigate the micro crack propagation and its direction underneath the excavating discs. A higher order displacement discontinuity method with quadratic displacement discontinuity elements is used to estimate the stress intensity factors near the crack tips. Rock cutting mechanisms under single and double type discs are simulated by the proposed numerical method.The main purposes of the present modeling are to simulate the chip formation process of indented rocks by single and double discs.The effects of specific disc parameters(except speed) on the thrust force Ft, the rolling force Fr, and the specific energy ES are investigated. It has been shown that the specific energy(energy required to cut through a unit volume of rock) of the double disc is less than that of the single disc. Crack propagation in rocks under disc cutters is numerically modeled and the optimum ratio of disc spacing S to penetration depth Pd(i.e. S/Pd ratio) of about 10 is obtained, which is in good agreement with the theoretical and experimental results cited in the literature.展开更多
Based on the simplification of cutting process,a series of numerical simulations were conducted using a 2-D discrete element method to explore the effects of embedded cracks with different dip angles on the rock fragm...Based on the simplification of cutting process,a series of numerical simulations were conducted using a 2-D discrete element method to explore the effects of embedded cracks with different dip angles on the rock fragmentation process,cutting characteristics and breaking efficiency.The results show that the simulated results are in a good agreement with previous theoretical study.The main crack propagates to the top tip of embedded crack,except when the dip angle is 90°.Side cracks which are more fully developed in the rocks containing embedded cracks tend to propagate towards the free surface.According to the history of vertical cutting force,it is shown that the peak force is decreased by embedded cracks.The study on cutting efficiency was conducted by combining the quantity of crack and cutting energy.And the results show that breaking efficiency can be treated as a decreasing or a increasing function when the dip angle is less or larger than 30°,respectively.Breaking efficiency is higher than that in intact rock when the dip angle is larger than 45°.展开更多
文摘考虑多晶材料疲劳行为受晶粒拓扑结构的影响,建立了多晶材料Voronoi镶嵌的有限元微观结构模型,并在Voronoi边界插入内聚力单元模拟晶界特性。在此基础上,基于Abaqus双线性内聚力本构关系,建立了内聚力单元的疲劳累积损伤模型,定义了损伤变量,损伤判据以及损伤规律;通过Abaqus UMAT(User-defined material mechanical behavior)子程序编写了内聚力单元疲劳损伤的数值计算程序,实现了对微观结构下材料滚动接触疲劳裂纹萌生与扩展的行为模拟。以滚子与滚道的接触疲劳分析为例,将数值模拟结果与实测结果相比,表明了该数值仿真方法模拟多晶材料裂纹形成和扩展的有效性。
文摘页岩水力压裂技术是目前最常用的储层改造技术。页岩中存在具有方向性的不连续结构面,使其具有较为明显的层理效应,主要体现为页岩的横观各向同性。为了在水力压裂中反映层理效应,在考虑层理效应的黏结单元方法(cohesive finite element method,CFEM)基础上,考虑力学场和渗流场的相互作用,推导界面单元全流固耦合方程。最后,通过水力压裂试验结果与模拟计算结果的对比,展现在不同地应力条件下页岩层理性对水力裂纹扩展的影响作用。在应用该方法时可以根据给定层理方向,模型自动将与层理平行(或近于平行)的黏结面转化为层面界面,而且在整个模拟过程中,不需要引入额外断裂准则,也不需要网格修正,为页岩水压致裂问题的数值模拟提供了新途径。
基金Projects(52474167,52104109)supported by the National Natural Science Foundation of ChinaProject(2022JJ40602)supported by the Natural Science Foundation of Hunan Province,China。
文摘The effect of freeze-thaw(F-T)cycles on the mechanical behaviors and internal mechanism of rock mass is a critical research topic.In permafrost or seasonally frozen regions,F-T cycles have adverse effects on the mechanical properties of rock mass,leading to many serious disasters in mining and geotechnical operations.In this paper,uniaxial compression tests are carried out on cyan sandstone after different F-T cycles.The failure modes and damage evolution of cyan sandstone under F-T cycles are studied.In addition,from the perspective of fracture and pore volume,the calculation equations of rock strain under frost heaving pressure and F-T cycles are established and verified with the corresponding laboratory tests.Subsequently,based on the classical damage theory,the F-T damage variables of cyan sandstone under different F-T cycles are calculated,and the meso-damage calculation model of cyan sandstone under F-T-loading coupling conditions is derived.Furthermore,through the discrete element numerical simulation software(PFC^(3D)),the microscopic damage evolution process of cyan sandstone under uniaxial compression after F-T cycles is studied,including the change of microcracks number,distribution of microcracks,and the acoustic emission(AE)count.The goal of this study is to investigate the damage evolution mechanism of rock from the mesoscopic and microscopic aspects,which has certain guiding value for accurately understanding the damage characteristics of rock in cold regions.
基金Projects(10972238,51074071,50974059)supported by the National Natural Science Foundation of ChinaProject(10JJ3007)supported by the Natural Science Foundation of Hunan Province,China+1 种基金Project(11C0539)supported by Scientific Research Fund of Hunan Provincial Education Department,ChinaProject(200905)supported by Open Research Fund of Hunan Provincial Key of Safe Mining Techniques of Coal Mines,China
文摘By considering the effect of hydraulic pressure filled in wing crack and the connected part of main crack on the stress intensity factor at wing crack tip, a new wing crack model exerted by hydraulic pressure and far field stresses was proposed. By introducing the equivalent crack length lcq of wing crack, two terms make up the stress intensity factor K1 at wing crack tip： one is the component K（1） for a single isolated straight wing crack of length 2l subjected to hydraulic pressure in wing crack and far field stresses, and the other is the component K1^（2） due to the effective shear stress induced by the presence of the equivalent main crack. The FEM model of wing crack propagation subjected to hydraulic pressure and far field stresses was also established according to different side pressure coefficients and hydraulic pressures in crack. The result shows that a good agreement is found between theoretical model of wing crack proposed and finite element method （FEM）. In theory, an unstable crack propagation is shown if there is high hydraulic pressure and lateral tension. The wing crack model proposed can provide references for studying on hydraulic fracturing in rock masses.
文摘The present research is focused on the numerical crack coalescence analysis of the micro-cracks and cracks produced during the cutting action of TBM disc cutters. The linear elastic fracture mechanics(LEFM) concepts and the maximum tangential stress criterion are used to investigate the micro crack propagation and its direction underneath the excavating discs. A higher order displacement discontinuity method with quadratic displacement discontinuity elements is used to estimate the stress intensity factors near the crack tips. Rock cutting mechanisms under single and double type discs are simulated by the proposed numerical method.The main purposes of the present modeling are to simulate the chip formation process of indented rocks by single and double discs.The effects of specific disc parameters(except speed) on the thrust force Ft, the rolling force Fr, and the specific energy ES are investigated. It has been shown that the specific energy(energy required to cut through a unit volume of rock) of the double disc is less than that of the single disc. Crack propagation in rocks under disc cutters is numerically modeled and the optimum ratio of disc spacing S to penetration depth Pd(i.e. S/Pd ratio) of about 10 is obtained, which is in good agreement with the theoretical and experimental results cited in the literature.
基金Project(2013CB035401)supported by the National Basic Research Program of ChinaProject(51174228)supported by the National Natural Science Foundation of ChinaProject(71380100003)supported by Hunan Provincial Innovation Foundation for Postgraduate,China
文摘Based on the simplification of cutting process,a series of numerical simulations were conducted using a 2-D discrete element method to explore the effects of embedded cracks with different dip angles on the rock fragmentation process,cutting characteristics and breaking efficiency.The results show that the simulated results are in a good agreement with previous theoretical study.The main crack propagates to the top tip of embedded crack,except when the dip angle is 90°.Side cracks which are more fully developed in the rocks containing embedded cracks tend to propagate towards the free surface.According to the history of vertical cutting force,it is shown that the peak force is decreased by embedded cracks.The study on cutting efficiency was conducted by combining the quantity of crack and cutting energy.And the results show that breaking efficiency can be treated as a decreasing or a increasing function when the dip angle is less or larger than 30°,respectively.Breaking efficiency is higher than that in intact rock when the dip angle is larger than 45°.
