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.展开更多
Cables composed of rare-earth barium copper oxide(REBCO)tapes have been extensively used in various superconducting devices.In recent years,conductor on round core(CORC)cable has drawn the attention of researchers wit...Cables composed of rare-earth barium copper oxide(REBCO)tapes have been extensively used in various superconducting devices.In recent years,conductor on round core(CORC)cable has drawn the attention of researchers with its outstanding current-carrying capacity and mechanical properties.The REBCO tapes are wound spirally on the surface of CORC cable.Under extreme loadings,the REBCO tapes with layered composite structures are vulnerable,which can lead to degradation of critical current and even quenching of superconducting devices.In this paper,we simulate the deformation of CORC cable under external loads,and analyze the damage inside the tape with the cohesive zone model(CZM).Firstly,the fabrication and cabling of CORC are simulated,and the stresses and strains generated in the tape are extracted as the initial condition of the next step.Then,the tension and bending loads are applied to CORC cable,and the damage distribution inside the tape is presented.In addition,the effects of some parameters on the damage are discussed during the bending simulations.展开更多
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.展开更多
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.展开更多
Shale contains distributed directional bedding planes,which make the shale transverse isotropic.To model shale with consideration of the bedding planes,a cohesive finite element method(CFEM)is developed based on the r...Shale contains distributed directional bedding planes,which make the shale transverse isotropic.To model shale with consideration of the bedding planes,a cohesive finite element method(CFEM)is developed based on the randomized triangular mesh.The interface orientation generated from such mesh tends to be uniformly distributed with the element number increasing.To represent the bedding plane,the interfaces aligned with the bedding plane are assigned the cohesive law that characterizes the bedding plane while the other interfaces are assigned the cohesive law that characterizes the matrix.By this means,the anisotropy characteristics of the stiffness and the strength of shale are well represented.The simulation examples demonstrate that the bedding plane has a significant influence on the fracture trajectory,which is consistent with the observation in the experiment.It is suggested that this modeling method of shale is feasible.It provides an alternative approach to fracture simulation in shale.展开更多
In this paper,a coupled thermo-hydro-mechanical(THM)simulation in a faulted deformable porous medium is presented.This model involves solving the mass conservation,linear momentum balance,and energy balance equations ...In this paper,a coupled thermo-hydro-mechanical(THM)simulation in a faulted deformable porous medium is presented.This model involves solving the mass conservation,linear momentum balance,and energy balance equations which are derived from the Biot’s consolidation theory.Fluid pore pressure,solid displacement,and temperature are chosen as initial variables in these equations,and the finite element method in combination with the interface element is used for spatial discretization of continuous and discontinuities(fault)parts of the medium to solve the equations.The main purpose of this study is providing precise formulations,applicability,and ability of the triple-node zero-thickness interface element in THM modeling of faults.It should be noted that the system of equations is solved using a computer code written in Matlab program.In order to verify the developed method,simulations of index problems such as Mandel’s problem,and coupled modeling of a faulted porous medium and a faulted aquifer are presented.The modeling results obtained from the developed method show a very good agreement with those by other modeling methods,which indicates its accuracy.展开更多
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.展开更多
A new test method was proposed to evaluate the cohesive strength of composite laminates. Cohesive strength and the critical strain energy for Mode-II interlamiar fracture of E-glass/epoxy woven fabrication were determ...A new test method was proposed to evaluate the cohesive strength of composite laminates. Cohesive strength and the critical strain energy for Mode-II interlamiar fracture of E-glass/epoxy woven fabrication were determined from the single lap joint(SLJ) and end notch flexure(ENF) test, respectively. In order to verify their adequacy, a cohesive zone model simulation based on interface finite elements was performed. A closed form solution for determination of the penalty stiffness parameter was proposed. Modified form of Park-Paulino-Roesler traction-separation law was provided and conducted altogether with trapezoidal and bilinear mixed-mode damage models to simulate damage using Abaqus cohesive elements. It was observed that accurate damage prediction and numerical convergence were obtained using the proposed penalty stiffness. Comparison between three damage models reveals that good simulation of fracture process zone and delamination prediction were obtained using the modified PPR model as damage model. Cohesive zone length as a material property was determined. To ensure the sufficient dissipation of energy, it was recommended that at least 4 elements should span cohesive zone length.展开更多
An appropriate coupled cohesive law for predicting the mixed mode failure is established by combining normal separation and tangential separation of surfaces in the cohesive zone model (CZM) and the cohesive element...An appropriate coupled cohesive law for predicting the mixed mode failure is established by combining normal separation and tangential separation of surfaces in the cohesive zone model (CZM) and the cohesive element method. The Xu-Needleman exponential cohesive law with the fully shear failure mechanism is one of the most popular models. Based on the proposed consistently coupled rule/principle, the Xu-Needleman law with the fully shear failure mechanism is proved to be a non-consistently coupled cohesive law by analyzing the surface separation work. It is shown that the Xu-Needleman law is only valid in the mixed mode fracture when the normal separation work equals the tangential separation work. Based on the consistently coupled principle and the modification of the Xu-Needleman law, a consistently coupled cohesive (CCC) law is given. It is shown that the proposed CCC law has already overcome the non-consistency defect of the Xu-Needleman law with great promise in mixed mode analyses.展开更多
文摘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.
基金Project supported by the National Natural Science Foundation of China(Nos.U2241267,1217215511872195)。
文摘Cables composed of rare-earth barium copper oxide(REBCO)tapes have been extensively used in various superconducting devices.In recent years,conductor on round core(CORC)cable has drawn the attention of researchers with its outstanding current-carrying capacity and mechanical properties.The REBCO tapes are wound spirally on the surface of CORC cable.Under extreme loadings,the REBCO tapes with layered composite structures are vulnerable,which can lead to degradation of critical current and even quenching of superconducting devices.In this paper,we simulate the deformation of CORC cable under external loads,and analyze the damage inside the tape with the cohesive zone model(CZM).Firstly,the fabrication and cabling of CORC are simulated,and the stresses and strains generated in the tape are extracted as the initial condition of the next step.Then,the tension and bending loads are applied to CORC cable,and the damage distribution inside the tape is presented.In addition,the effects of some parameters on the damage are discussed during the bending simulations.
文摘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.
基金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.
基金supported by the National Natural Science Foundation of China(No.11772190)
文摘Shale contains distributed directional bedding planes,which make the shale transverse isotropic.To model shale with consideration of the bedding planes,a cohesive finite element method(CFEM)is developed based on the randomized triangular mesh.The interface orientation generated from such mesh tends to be uniformly distributed with the element number increasing.To represent the bedding plane,the interfaces aligned with the bedding plane are assigned the cohesive law that characterizes the bedding plane while the other interfaces are assigned the cohesive law that characterizes the matrix.By this means,the anisotropy characteristics of the stiffness and the strength of shale are well represented.The simulation examples demonstrate that the bedding plane has a significant influence on the fracture trajectory,which is consistent with the observation in the experiment.It is suggested that this modeling method of shale is feasible.It provides an alternative approach to fracture simulation in shale.
文摘In this paper,a coupled thermo-hydro-mechanical(THM)simulation in a faulted deformable porous medium is presented.This model involves solving the mass conservation,linear momentum balance,and energy balance equations which are derived from the Biot’s consolidation theory.Fluid pore pressure,solid displacement,and temperature are chosen as initial variables in these equations,and the finite element method in combination with the interface element is used for spatial discretization of continuous and discontinuities(fault)parts of the medium to solve the equations.The main purpose of this study is providing precise formulations,applicability,and ability of the triple-node zero-thickness interface element in THM modeling of faults.It should be noted that the system of equations is solved using a computer code written in Matlab program.In order to verify the developed method,simulations of index problems such as Mandel’s problem,and coupled modeling of a faulted porous medium and a faulted aquifer are presented.The modeling results obtained from the developed method show a very good agreement with those by other modeling methods,which indicates its accuracy.
文摘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.
文摘A new test method was proposed to evaluate the cohesive strength of composite laminates. Cohesive strength and the critical strain energy for Mode-II interlamiar fracture of E-glass/epoxy woven fabrication were determined from the single lap joint(SLJ) and end notch flexure(ENF) test, respectively. In order to verify their adequacy, a cohesive zone model simulation based on interface finite elements was performed. A closed form solution for determination of the penalty stiffness parameter was proposed. Modified form of Park-Paulino-Roesler traction-separation law was provided and conducted altogether with trapezoidal and bilinear mixed-mode damage models to simulate damage using Abaqus cohesive elements. It was observed that accurate damage prediction and numerical convergence were obtained using the proposed penalty stiffness. Comparison between three damage models reveals that good simulation of fracture process zone and delamination prediction were obtained using the modified PPR model as damage model. Cohesive zone length as a material property was determined. To ensure the sufficient dissipation of energy, it was recommended that at least 4 elements should span cohesive zone length.
基金Project supported by the National Natural Science Foundation of China(Nos.50878117 and 51038006)the China Scholarship Council Project(No.M.H.HE-2009621076)the Tsinghua University Initiative Scientific Research Program(No.20101081766)
文摘An appropriate coupled cohesive law for predicting the mixed mode failure is established by combining normal separation and tangential separation of surfaces in the cohesive zone model (CZM) and the cohesive element method. The Xu-Needleman exponential cohesive law with the fully shear failure mechanism is one of the most popular models. Based on the proposed consistently coupled rule/principle, the Xu-Needleman law with the fully shear failure mechanism is proved to be a non-consistently coupled cohesive law by analyzing the surface separation work. It is shown that the Xu-Needleman law is only valid in the mixed mode fracture when the normal separation work equals the tangential separation work. Based on the consistently coupled principle and the modification of the Xu-Needleman law, a consistently coupled cohesive (CCC) law is given. It is shown that the proposed CCC law has already overcome the non-consistency defect of the Xu-Needleman law with great promise in mixed mode analyses.