In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hy...In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.展开更多
A two-dimensional axisymmetric finite element model based on an improved cohesive element method was developed to simulate interfacial debonding, sliding friction, and residual thermal stresses in SiC composites durin...A two-dimensional axisymmetric finite element model based on an improved cohesive element method was developed to simulate interfacial debonding, sliding friction, and residual thermal stresses in SiC composites during single-fiber push-out tests to extract the interfacial bond strength and frictional stress. The numerical load–displacement curves agree well with experimental curves,indicating that this cohesive element method can be used for calculating the interfacial properties of SiC composites.The simulation results show that cracks are most likely to occur at the ends of the experimental sample, where the maximum shear stress is observed and that the interfacial shear strength and constant sliding friction stress decrease with an increase in temperature. Moreover, the load required to cause complete interfacial failure increases with the increase in critical shear strength, and the composite materials with higher fiber volume fractions have higher bearing capacities. In addition, the initial failure load increases with an increase in interphase thickness.展开更多
冰荷载是影响海上风力机安全运行的重要决定性因素,严重时会致使海上风力机结构发生冰激振动破坏和冰激疲劳失效。该研究基于粘聚单元(cohesive element method,CEM)-有限元(finite element method,FEM)耦合方法,通过非线性分布弹簧考虑...冰荷载是影响海上风力机安全运行的重要决定性因素,严重时会致使海上风力机结构发生冰激振动破坏和冰激疲劳失效。该研究基于粘聚单元(cohesive element method,CEM)-有限元(finite element method,FEM)耦合方法,通过非线性分布弹簧考虑桩-土相互作用,建立海冰-寒区单桩海上风力机结构在风-冰联合作用下的整体耦合冰激振动非线性有限元模型。进而,基于非线性数值仿真工具LS-DYNA,分别模拟冰与直立结构和带有抗冰锥的基础结构相互作用过程,并与现有的挤压和弯曲冰力模型进行对比,验证该研究模拟动冰荷载的准确性,讨论两种冰破坏模式下动冰荷载的变化规律。最后,为解决海上风力机发生的强烈冰激振动问题,分别采用振动控制方法和施加抗冰锥的方式,开展风、冰联合作用下海上风力机动冰力和动力响应研究,对比分析以上两种减振方式的减振机理和减振效果差异。结果表明,虽然抗冰锥可明显降低冰-海上风力机相互作用的动冰荷载幅值,然而采用振动控制策略的海上风力机减振效果明显优于抗冰锥。因此,在海上风力机的冰激结构损伤研究和抗冰设计中必须分别考虑以上两种减振方式对结构的影响。展开更多
Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient applicat...Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.展开更多
Cohesive element is developed from the Dugdal-Barenblatt model in the field of fracture mechanics. The mechanical constitutive relation of cohesive element can be artificially assumed depending on the specific applica...Cohesive element is developed from the Dugdal-Barenblatt model in the field of fracture mechanics. The mechanical constitutive relation of cohesive element can be artificially assumed depending on the specific applications. It has been successfully applied in the study of crystal plasticity/brittle fracture process and decohesion between delaminations. In this paper, tensile experiments of large steel plate with different length of pre-existing cracks are conducted. Based on commercial software ABAQUS, cohesive element is adopted to simulate the tensile tests, and appropriate parameter values are obtained by fitting displacement-force curves. Using these parameters, a numerical method is presented by applying cohesive element to thermo-elastic-plastic finite element method (TEP-FEM) to simulate plate rigid restraint cracking (PRRC) tests. By changing constitutive relation of cohesive element, dimensions of the model and welding conditions, the influence of welding restraint intensity and welding conditions on the crack propagation are discussed, respectively. Three types of welding cold cracking are simulated. Significant influence of welding cold cracking on resistant stress in welding line is captured by this numerical method.展开更多
In order to investigate the fatigue behavior of asphalt concrete, a new numerical approach based on a bi-linear cohesive zone model (CZM) is developed. Integrated with the CZM, a fatigue damage evolution model is es...In order to investigate the fatigue behavior of asphalt concrete, a new numerical approach based on a bi-linear cohesive zone model (CZM) is developed. Integrated with the CZM, a fatigue damage evolution model is established to indicate the gradual degradation of cohesive properties of asphalt concrete under cyclic loading. Then the model is implemented in the finite element software ABAQUS through a user-defined subroutine. Based on the proposed model, an indirect tensile fatigue test is finally simulated. The fatigue lives obtained through numerical analysis show good agreement with laboratory results. Fatigue damage accumulates in a nonlinear manner during the cyclic loading process and damage initiation phase is the major part of fatigue failure. As the stress ratio increases, the time of the steady damage growth stage decreases significantly. It is found that the proposed fatigue damage evolution model can serve as an accurate and efficient tool for the prediction of fatigue damage of asphalt concrete.展开更多
The gradual increase in shipping and drilling activities in the Arctic regions has resulted in the increased importance of studying the structural safety of polar ships in various ice conditions.Rafted ice refers to a...The gradual increase in shipping and drilling activities in the Arctic regions has resulted in the increased importance of studying the structural safety of polar ships in various ice conditions.Rafted ice refers to a type of accumulated and overlapped sea ice;it is driven by external forces,such as wind and waves,and may exert high loads on ships and threaten their structural safety.Therefore,the properties of rafted ice and the construction of numerical models should be studied before exploring the interaction and collision between ships and rafted ice.Based on the nonlinear finite-element method,this paper introduces the cohesive element model for the simulation of rafted ice.The interaction between ships and rafted ice is studied,and the ice force of the hull is obtained.Numerical simulation results are compared with model test findings,and the effectiveness of the cohesive element method in the construction of the model of rafted ice materials is verified.On this basis,a multilayer rafted ice model is constructed,and its interaction with the ship is studied.The research unveils that rafted ice parts impede crack generation and slow down crack propagation to a certain extent.展开更多
Indirect fracturing in the roof of broken soft coal seams has been demonstrated to be a feasible technology.In this work,the No.5 coal seam in the Hancheng block was taken as the research object.Based on the findings ...Indirect fracturing in the roof of broken soft coal seams has been demonstrated to be a feasible technology.In this work,the No.5 coal seam in the Hancheng block was taken as the research object.Based on the findings of true triaxial hydraulic fracturing experiments and field pilot under this technology and the cohesive element method,a 3D numerical model of indirect fracturing in the roof of broken soft coal seams was established,the fracture morphology propagation and evolution law under different conditions was investigated,and analysis of main controlling factors of fracture parameters was conducted with the combination weight method,which was based on grey incidence,analytic hierarchy process and entropy weight method.The results show that“士”-shaped fractures,T-shaped fractures,cross fractures,H-shaped fractures,and“干”-shaped fractures dominated by horizontal fractures were formed.Different parameter combinations can form different fracture morphologies.When the coal seam permeability is lower and the minimum horizontal principal stress difference between layers and fracturing fluid injection rate are both larger,it tends to form“士”-shaped fractures.When the coal seam permeability and minimum horizontal principal stress between layers and perforation position are moderate,cross fractures are easily generated.Different fracture parameters have different main controlling factors.Engineering factors of perforation location,fracturing fluid injection rate and viscosity are the dominant factors of hydraulic fracture shape parameters.This study can provide a reference for the design of indirect fracturing in the roof of broken soft coal seams.展开更多
The shear mechanical behavior is regarded as an essential factor affecting the stability of the surrounding rocks in underground engineering.The shear strength and failure mechanisms of layered rock are significantly ...The shear mechanical behavior is regarded as an essential factor affecting the stability of the surrounding rocks in underground engineering.The shear strength and failure mechanisms of layered rock are significantly affected by the foliation angles.Direct shear tests were conducted on cubic slate samples with foliation angles of 0°,30°,45°,60°,and 90°.The effect of foliation angles on failure patterns,acoustic emission(AE)characteristics,and shear strength parameters was analyzed.Based on AE characteristics,the slate failure process could be divided into four stages:quiet period,step-like increasing period,dramatic increasing period,and remission period.A new empirical expression of cohesion for layered rock was proposed,which was compared with linear and sinusoidal cohesion expressions based on the results made by this paper and previous experiments.The comparative analysis demonstrated that the new expression has better prediction ability than other expressions.The proposed empirical equation was used for direct shear simulations with the combined finite-discrete element method(FDEM),and it was found to align well with the experimental results.Considering both computational efficiency and accuracy,it was recommended to use a shear rate of 0.01 m/s for FDEM to carry out direct shear simulations.To balance the relationship between the number of elements and the simulation results in the direct shear simulations,the recommended element size is 1 mm.展开更多
文摘In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.
基金supported by the National Natural Science Foundation of China(No.11405124)Science Challenge Project(No.TZ2018004)+1 种基金Natural Science Basic Research Plan in Shaanxi Province of China(No.2015JQ1030)the Shaanxi Province Postdoctoral Science Foundation(2014)
文摘A two-dimensional axisymmetric finite element model based on an improved cohesive element method was developed to simulate interfacial debonding, sliding friction, and residual thermal stresses in SiC composites during single-fiber push-out tests to extract the interfacial bond strength and frictional stress. The numerical load–displacement curves agree well with experimental curves,indicating that this cohesive element method can be used for calculating the interfacial properties of SiC composites.The simulation results show that cracks are most likely to occur at the ends of the experimental sample, where the maximum shear stress is observed and that the interfacial shear strength and constant sliding friction stress decrease with an increase in temperature. Moreover, the load required to cause complete interfacial failure increases with the increase in critical shear strength, and the composite materials with higher fiber volume fractions have higher bearing capacities. In addition, the initial failure load increases with an increase in interphase thickness.
文摘冰荷载是影响海上风力机安全运行的重要决定性因素,严重时会致使海上风力机结构发生冰激振动破坏和冰激疲劳失效。该研究基于粘聚单元(cohesive element method,CEM)-有限元(finite element method,FEM)耦合方法,通过非线性分布弹簧考虑桩-土相互作用,建立海冰-寒区单桩海上风力机结构在风-冰联合作用下的整体耦合冰激振动非线性有限元模型。进而,基于非线性数值仿真工具LS-DYNA,分别模拟冰与直立结构和带有抗冰锥的基础结构相互作用过程,并与现有的挤压和弯曲冰力模型进行对比,验证该研究模拟动冰荷载的准确性,讨论两种冰破坏模式下动冰荷载的变化规律。最后,为解决海上风力机发生的强烈冰激振动问题,分别采用振动控制方法和施加抗冰锥的方式,开展风、冰联合作用下海上风力机动冰力和动力响应研究,对比分析以上两种减振方式的减振机理和减振效果差异。结果表明,虽然抗冰锥可明显降低冰-海上风力机相互作用的动冰荷载幅值,然而采用振动控制策略的海上风力机减振效果明显优于抗冰锥。因此,在海上风力机的冰激结构损伤研究和抗冰设计中必须分别考虑以上两种减振方式对结构的影响。
文摘Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.
文摘Cohesive element is developed from the Dugdal-Barenblatt model in the field of fracture mechanics. The mechanical constitutive relation of cohesive element can be artificially assumed depending on the specific applications. It has been successfully applied in the study of crystal plasticity/brittle fracture process and decohesion between delaminations. In this paper, tensile experiments of large steel plate with different length of pre-existing cracks are conducted. Based on commercial software ABAQUS, cohesive element is adopted to simulate the tensile tests, and appropriate parameter values are obtained by fitting displacement-force curves. Using these parameters, a numerical method is presented by applying cohesive element to thermo-elastic-plastic finite element method (TEP-FEM) to simulate plate rigid restraint cracking (PRRC) tests. By changing constitutive relation of cohesive element, dimensions of the model and welding conditions, the influence of welding restraint intensity and welding conditions on the crack propagation are discussed, respectively. Three types of welding cold cracking are simulated. Significant influence of welding cold cracking on resistant stress in welding line is captured by this numerical method.
基金The Open Research Fund of Key Laboratory of Highway Engineering of Sichuan Province of Southw est Jiaotong University (No.LHTE002201102)
文摘In order to investigate the fatigue behavior of asphalt concrete, a new numerical approach based on a bi-linear cohesive zone model (CZM) is developed. Integrated with the CZM, a fatigue damage evolution model is established to indicate the gradual degradation of cohesive properties of asphalt concrete under cyclic loading. Then the model is implemented in the finite element software ABAQUS through a user-defined subroutine. Based on the proposed model, an indirect tensile fatigue test is finally simulated. The fatigue lives obtained through numerical analysis show good agreement with laboratory results. Fatigue damage accumulates in a nonlinear manner during the cyclic loading process and damage initiation phase is the major part of fatigue failure. As the stress ratio increases, the time of the steady damage growth stage decreases significantly. It is found that the proposed fatigue damage evolution model can serve as an accurate and efficient tool for the prediction of fatigue damage of asphalt concrete.
基金The National Natural Science Foundation of China(Grant Nos.52192693,52192690,51979051,51979056 and U20A20327)the National Key Research and Development Program of China(Grant No.2021YFC2803400)。
文摘The gradual increase in shipping and drilling activities in the Arctic regions has resulted in the increased importance of studying the structural safety of polar ships in various ice conditions.Rafted ice refers to a type of accumulated and overlapped sea ice;it is driven by external forces,such as wind and waves,and may exert high loads on ships and threaten their structural safety.Therefore,the properties of rafted ice and the construction of numerical models should be studied before exploring the interaction and collision between ships and rafted ice.Based on the nonlinear finite-element method,this paper introduces the cohesive element model for the simulation of rafted ice.The interaction between ships and rafted ice is studied,and the ice force of the hull is obtained.Numerical simulation results are compared with model test findings,and the effectiveness of the cohesive element method in the construction of the model of rafted ice materials is verified.On this basis,a multilayer rafted ice model is constructed,and its interaction with the ship is studied.The research unveils that rafted ice parts impede crack generation and slow down crack propagation to a certain extent.
基金National Natural Science Foundation of China(11672333).
文摘Indirect fracturing in the roof of broken soft coal seams has been demonstrated to be a feasible technology.In this work,the No.5 coal seam in the Hancheng block was taken as the research object.Based on the findings of true triaxial hydraulic fracturing experiments and field pilot under this technology and the cohesive element method,a 3D numerical model of indirect fracturing in the roof of broken soft coal seams was established,the fracture morphology propagation and evolution law under different conditions was investigated,and analysis of main controlling factors of fracture parameters was conducted with the combination weight method,which was based on grey incidence,analytic hierarchy process and entropy weight method.The results show that“士”-shaped fractures,T-shaped fractures,cross fractures,H-shaped fractures,and“干”-shaped fractures dominated by horizontal fractures were formed.Different parameter combinations can form different fracture morphologies.When the coal seam permeability is lower and the minimum horizontal principal stress difference between layers and fracturing fluid injection rate are both larger,it tends to form“士”-shaped fractures.When the coal seam permeability and minimum horizontal principal stress between layers and perforation position are moderate,cross fractures are easily generated.Different fracture parameters have different main controlling factors.Engineering factors of perforation location,fracturing fluid injection rate and viscosity are the dominant factors of hydraulic fracture shape parameters.This study can provide a reference for the design of indirect fracturing in the roof of broken soft coal seams.
基金support from the Natural Science Foundation of China(Grant Nos.41941018,U21A20153,42177140).
文摘The shear mechanical behavior is regarded as an essential factor affecting the stability of the surrounding rocks in underground engineering.The shear strength and failure mechanisms of layered rock are significantly affected by the foliation angles.Direct shear tests were conducted on cubic slate samples with foliation angles of 0°,30°,45°,60°,and 90°.The effect of foliation angles on failure patterns,acoustic emission(AE)characteristics,and shear strength parameters was analyzed.Based on AE characteristics,the slate failure process could be divided into four stages:quiet period,step-like increasing period,dramatic increasing period,and remission period.A new empirical expression of cohesion for layered rock was proposed,which was compared with linear and sinusoidal cohesion expressions based on the results made by this paper and previous experiments.The comparative analysis demonstrated that the new expression has better prediction ability than other expressions.The proposed empirical equation was used for direct shear simulations with the combined finite-discrete element method(FDEM),and it was found to align well with the experimental results.Considering both computational efficiency and accuracy,it was recommended to use a shear rate of 0.01 m/s for FDEM to carry out direct shear simulations.To balance the relationship between the number of elements and the simulation results in the direct shear simulations,the recommended element size is 1 mm.