The multi-scale modeling combined with the cohesive zone model(CZM)and the molecular dynamics(MD)method were preformed to simulate the crack propagation in NiTi shape memory alloys(SMAs).The metallographic microscope ...The multi-scale modeling combined with the cohesive zone model(CZM)and the molecular dynamics(MD)method were preformed to simulate the crack propagation in NiTi shape memory alloys(SMAs).The metallographic microscope and image processing technology were employed to achieve a quantitative grain size distribution of NiTi alloys so as to provide experimental data for molecular dynamics modeling at the atomic scale.Considering the size effect of molecular dynamics model on material properties,a reasonable modeling size was provided by taking into account three characteristic dimensions from the perspective of macro,meso,and micro scales according to the Buckinghamπtheorem.Then,the corresponding MD simulation on deformation and fracture behavior was investigated to derive a parameterized traction-separation(T-S)law,and then it was embedded into cohesive elements of finite element software.Thus,the crack propagation behavior in NiTi alloys was reproduced by the finite element method(FEM).The experimental results show that the predicted initiation fracture toughness is in good agreement with experimental data.In addition,it is found that the dynamics initiation fracture toughness increases with decreasing grain size and increasing loading velocity.展开更多
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.展开更多
The cohesive zone model(CZM)has been used widely and successfully in fracture propagation,but some basic problems are still to be solved.In this paper,artificial compliance and discontinuous force in CZM are investiga...The cohesive zone model(CZM)has been used widely and successfully in fracture propagation,but some basic problems are still to be solved.In this paper,artificial compliance and discontinuous force in CZM are investigated.First,theories about the cohesive element(local coordinate system,stiffness matrix,and internal nodal force)are presented.The local coordinate system is defined to obtain local separation;the stiffness matrix for an eight-node cohesive element is derived from the calculation of strain energy;internal nodal force between the cohesive element and bulk element is obtained from the principle of virtual work.Second,the reason for artificial compliance is explained by the effective stiffnesses of zero-thickness and finite-thickness cohesive elements.Based on the effective stiffness,artificial compliance can be completely removed by adjusting the stiffness of the finite-thickness cohesive element.This conclusion is verified from 1D and 3D simulations.Third,three damage evolution methods(monotonically increasing effective separation,damage factor,and both effective separation and damage factor)are analyzed.Under constant unloading and reloading conditions,the monotonically increasing damage factor method without discontinuous force and healing effect is a better choice than the other two methods.The proposed improvements are coded in LS-DYNA user-defined material,and a drop weight tear test verifies the improvements.展开更多
Amid urbanization and the continuous expansion of transportation networks,the necessity for tunnel construction and maintenance has become paramount.Addressing this need requires the investigation of efficient,economi...Amid urbanization and the continuous expansion of transportation networks,the necessity for tunnel construction and maintenance has become paramount.Addressing this need requires the investigation of efficient,economical,and robust tunnel reinforcement techniques.This paper explores fiber reinforced polymer(FRP)and steel fiber reinforced concrete(SFRC)technologies,which have emerged as viable solutions for enhancing tunnel structures.FRP is celebrated for its lightweight and high-strength attributes,effectively augmenting load-bearing capacity and seismic resistance,while SFRC’s notable crack resistance and longevity potentially enhance the performance of tunnel segments.Nonetheless,current research predominantly focuses on experimental analysis,lacking comprehensive theoretical models.To bridge this gap,the cohesive zone model(CZM),which utilizes cohesive elements to characterize the potential fracture surfaces of concrete/SFRC,the rebar-concrete interface,and the FRP-concrete interface,was employed.A modeling approach was subsequently proposed to construct a tunnel segment model reinforced with either SFRC or FRP.Moreover,the corresponding mixed-mode constitutive models,considering interfacial friction,were integrated into the proposed model.Experimental validation and numerical simulations corroborated the accuracy of the proposed model.Additionally,this study examined the reinforcement design of tunnel segments.Through a numerical evaluation,the effectiveness of innovative reinforcement schemes,such as substituting concrete with SFRC and externally bonding FRP sheets,was assessed utilizing a case study from the Fuzhou Metro Shield Tunnel Construction Project.展开更多
A new cyclic cohesive zone fatigue damage model is proposed to address the fatigue problem spanning highand low cycle stages.The new damage model is integrated with the damage extrapolation technique to improvecalcula...A new cyclic cohesive zone fatigue damage model is proposed to address the fatigue problem spanning highand low cycle stages.The new damage model is integrated with the damage extrapolation technique to improvecalculation efficiency.The model’s effectiveness in regulating the low-cycle fatigue evolution rate,overall fatiguedamage evolution rate,and stress level at the fatigue turning point is assessed through the comparison of the S-Ncurves.The fatigue damage model’s high precision is proved based on the minor deviation of stress at the turningpoint of the S-N curve from the actual scenario.Finally,the fatigue damage evolution is simulated consideringthe effects of pre-load pressure and welding residual stress.It is observed that laser welding induces a significantresidual tensile stress,accelerating fatigue damage evolution,while compressive loading impedes fatigue damageprogression.展开更多
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.展开更多
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.展开更多
This paper concentrates on simulating fracture in thin walled single-lap joints connected by resistance spot-welding(RSW)process which were subjected to tensile loading.For this purpose,three sets of lap-joints with d...This paper concentrates on simulating fracture in thin walled single-lap joints connected by resistance spot-welding(RSW)process which were subjected to tensile loading.For this purpose,three sets of lap-joints with different spot configurations were tested to achieve the joints’tensile behavior.To simulate the joints tensile behavior,firstly a 2D axisymmetric finite element(FE)model was used to calculate residual stresses induced during the welding process.Then the results were transferred to 3D models as pre-stress.In this step,cohesive zone model(CZM)technique was used to simulate fracture in the models under tensile load.Cohesive zone parameters were extracted using coach-peel and shear lap specimens.The results were employed to simulate deformation and failure in single lap spot weld samples.It has been shown that considering the residual stresses in simulating deformation and fracture load enables quite accurate predictions.展开更多
The softening-melting characteristics of ferrous burden play a crucial role in the thickness and position of the cohesive zone.The influence of the basicity and experimental atmosphere on the softening-melting behavio...The softening-melting characteristics of ferrous burden play a crucial role in the thickness and position of the cohesive zone.The influence of the basicity and experimental atmosphere on the softening-melting behavior of primary slag under slag-coke interaction was investigated using in situ visualization method.The mechanism was analyzed using the FactSage software,X-ray diffraction,and electron probe microanalysis.The softening and melting temperatures of the samples increased with increasing basicity under different atmospheres.The difference between softening and melting temperatures is smaller in a H_(2) atmosphere than in a CO atmosphere;in H_(2) atmosphere,the range of softening zone in the cohesive zone was significantly thinner.The formed low-melting-point FeO-bearing phases decrease when H_(2) was used as the reducing agent.In addition,according to FactSage calculations,the high content of metallic iron reduced in the H_(2) atmosphere raised the softening temperature of the primary slag.It also narrowed and moved downward the cohesive zone due to an increase in softening and melting temperatures.Meanwhile,the increase in basicity promoted the decrease in liquid ratio and improved the permeability of cohesive zone.展开更多
Based on the principle of discrete element method （DEM）, a 2D slot model of a COREX melter gasifier was established to analyze the influence of cohesive zone shape on solid flow, including mass distribution, velocit...Based on the principle of discrete element method （DEM）, a 2D slot model of a COREX melter gasifier was established to analyze the influence of cohesive zone shape on solid flow, including mass distribution, velocity distribution, normal force distribution and porosity distribution at a microscopic level. The results show that the co- hesive zone shape almost does not affect the particle movement in the upper shaft and deadman shape. The particles in the lower central bottom experience large normal force to support the particles above them, while particles around the raceway and in the fast flow zone exhibit weak force network. The porosity distribution was also examined under three kinds of cohesive zones. Like the velocity distribution, the whole packed bed can be divided into four main re- gions. With the increase of cohesive zone position, the low porosity region located in the root of cohesive zone increa- ses. And the porosity distribution becomes asymmetric in the case of biased cohesive zone.展开更多
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.展开更多
The key parameters of the adhesive layer of a reinforcing patch are of great significance and affect the ability to suppress crack propagation in an Al–Li alloy patch-reinforced structure.This paper proposes a method...The key parameters of the adhesive layer of a reinforcing patch are of great significance and affect the ability to suppress crack propagation in an Al–Li alloy patch-reinforced structure.This paper proposes a method to determine the key parameters of the adhesive layer of adhesively bonded joints in the Al–Li alloy patch-reinforced structure.A zero-thickness cohesive zone model(CZM)was selected to simulate the adhesive layer’s fracture process,and an orthogonal simulation was designed to compare against the test results.A three-dimensional progressive damage model of an Al–Li alloy patch-reinforced structure with single-lap adhesively bonded joints was developed.The simulation’s results closely agree with the test results,demonstrating that this method of determining the key parameters is likely accurate.The results also verify the correctness of the cohesive strength and fracture energy,the two key parameters of the cohesive zone model.The model can accurately predict the strength and fracture process of adhesively bonded joints,and can be used in research to suppress crack propagation in Al–Li alloy patch-reinforced structures.展开更多
Investigation of paper cutting process is vital for the design of cutting tools,but the fracture mechanism of paper cutting is still unclear.Here,we focus on the cutting process of paper,including the key parameters o...Investigation of paper cutting process is vital for the design of cutting tools,but the fracture mechanism of paper cutting is still unclear.Here,we focus on the cutting process of paper,including the key parameters of cohesive zone model(CZM)for the orthotropic paper,to simulate the shear fracture process.Firstly,the material constants of the orthotropic paper are determined by longitudinal and transverse tensile test.Secondly,based on the tensile stressstrain curves,combined with damage theory and numerical simulations,the key parameters of the CZM for the orthotropic paper are obtained.Finally,a model III fracture is simulated to verify the accuracy of the model.Results show that the load-displacement curves obtained by the simulation is consistent with the test results.展开更多
Cohesive zone model was used to simulate two-dimensional plane strain crack propagation at the grain level model including grain boundary zones. Simulated results show that the original crack-tip may not be separated ...Cohesive zone model was used to simulate two-dimensional plane strain crack propagation at the grain level model including grain boundary zones. Simulated results show that the original crack-tip may not be separated firstly in an elastic-viscoplastic polycrystals. The grain interior＇s material properties （e.g. strain rate sensitivity） characterize the competitions between plastic and cohesive energy dissipation mechanisms. The higher the strain rate sensitivity is, the larger amount of the external work is transformed into plastic dissipation energy than into cohesive energy, which delays the cohesive zone rupturing. With the strain rate sensitivity decreased, the material property tends to approach the elastic-plastic responses. In this case, the plastic dissipation energy decreases and the cohesive dissipation energy increases which accelerates the cohesive zones debonding. Increasing the cohesive strength or the critical separation displacement will reduce the stress triaxiality at grain interiors and grain boundaries. Enhancing the cohesive zones ductility can improve the matrix materials resistance to void damage.展开更多
An improved interface cohesive zone model is developed for the simulation of interface contact, under mixed-mode loading. A new debonding initiation criterion and propagation of debonding law, taking into account the ...An improved interface cohesive zone model is developed for the simulation of interface contact, under mixed-mode loading. A new debonding initiation criterion and propagation of debonding law, taking into account the pressure stress influence on contact shear strength, is proposed. The model is implemented in a finite-element program using subroutine VUINTER of ABAQUS Explicit. An edge-notch four-point bending process and laminated vibration damping steel sheet punch forming test are simulated with the improved model in ABAQUS Explicit. The numerical predictions agree satisfactorily with the corresponding experimental results.展开更多
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.展开更多
The changes in the softening and melting behaviors of ferrous burden in the cohesive zone and the characteristics of the slag–iron–coke interface in a blast furnace were investigated by simulating an actual blast fu...The changes in the softening and melting behaviors of ferrous burden in the cohesive zone and the characteristics of the slag–iron–coke interface in a blast furnace were investigated by simulating an actual blast furnace under hydrogen-rich conditions.According to the variation in the transient shrinkage of the burden under different atmospheres,the shrinkage start temperature of the sinter was higher than that of the pellets.The negative shrinkage rate of the pellets was greater than that of the sinter.Additionally,the softening start temperature in the blast furnace decreased under hydrogen-rich conditions,giving the blast furnace a broader range of softening zones.The softening start temperatures of the pellets and sinter decreased from 1102 to 949℃ and 1152 to 1080℃,respectively.The hydrogen-rich traditional blast furnace conditions narrowed the melting zone temperature range and shifted it toward the high-temperature zone,significantly improving the burden layer permeability.However,under the hydrogen-rich oxygen blast furnace conditions,there were a decrease in the melting start temperature,a shift of the melting zone location to the low-temperature zone,and an increase in the burden layer permeability and pressure difference.A comparison of the slag–iron–coke interface characteristics under different atmospheric conditions showed that the carbon content in metallic iron decreased under hydrogen-rich traditional blast furnace conditions compared with traditional blast furnace conditions.Contrastingly,under hydrogen-rich oxygen blast furnace conditions,the carbon content in metallic iron increased compared with oxygen blast furnace conditions.These findings provide guidance for the development of low-carbon ironmaking processes in blast furnaces.展开更多
To reduce greenhouse gas emissions in ironmaking,the steel industry is advancing innovative low-carbon blast furnace(BF)technologies.A critical challenge for implementing such innovations lies in optimizing permeabili...To reduce greenhouse gas emissions in ironmaking,the steel industry is advancing innovative low-carbon blast furnace(BF)technologies.A critical challenge for implementing such innovations lies in optimizing permeability within the BF's cohesive zone(CZ),which directly impacts operational stability and efficiency.This study employs a coupled computational fluid dynamics-discrete element method(CFD-DEM)to calibrate Young's modulus by respectively fitting the relationship between Young's modulus and temperature,as well as pressure drop,based on a reported lab-scale softening and smelting experimental data of ore-coke heterogeneous alternating layer packed beds resembling BFs,and develops a softening particle model(SPM).The SPM establishes a temperature-dependent relationship between mechanical properties of softened ore particles and CZ conditions in industrial-scale BFs.Simulations of particle shrinkage behavior and pressure drop trends using the CFD-DEM-SPM framework demonstrate strong correlation with experimental data,validating its accuracy for predictive analysis.Furthermore,this study investigates how layer arrangement configurations,size ratios between ore and coke particles,and coke blending proportions influence CZ characteristics.Key findings identify an optimal batch weight configuration to enhance permeability within the CZ while maintaining operational stability.Additionally,results indicate that increasing the relative particle size of ore compared to coke or enhancing the proportion of blended coke in burden mixes improves CZ permeability,offering actionable strategies for reducing carbon intensity in BF operations.These insights provide critical guidance toward developing low-carbon BF processes compatible with global climate targets.展开更多
Cohesive zone plays a vital role in the stable operation of a blast furnace(BF),yet the complex phase change process of iron ore particles in this zone is still not well understood.In this study,a novel one-dimensiona...Cohesive zone plays a vital role in the stable operation of a blast furnace(BF),yet the complex phase change process of iron ore particles in this zone is still not well understood.In this study,a novel one-dimensional(1D)unsteady phase change model was developed to elucidate the heat transfer and melting mechanisms of iron ore particles.After model validation,the effects of several key operating parameters(e.g.,particle diameter,gas velocity,initial temperature)on the phase change behavior of iron ore particles were analyzed,and the joint effect of multiple parameters was discussed.The results show that larger-sized iron ore particles possess lower specific surface areas,which in turn reduces their convective heat absorption capacity.Consequently,the distance from the solid-liquid phase interface to the particle surface increases,thereby slowing down the movement of the phase interface and pro-longing the melting duration of the particles.Increasing the gas velocity and the initial temperature does not have a significant impact on reducing the duration of the complete melting process.Under the specified conditions,it is observed that increasing the gas velocity by 3-fold and 9-fold results in a reduction of the melting duration by 2.4%and 8.3%,respectively.Elevating the initial temperature of iron ore particles results in a decrease in the core-to-surface temperature difference,a slower heating rate,and a shorter duration to achieve melting.Among the factors affecting the melting process,the particle diameter is found to be the most significant in terms of the liquid phase precipitation,mushy zone thickness,and core-to-surface temperature difference of iron ore particles.展开更多
A coupled thermal-hydro-mechanical cohesive phase-field model for hydraulic fracturing in deep coal seams is presented.Heat exchange between the cold fluid and the hot rock is considered,and the thermal contribution t...A coupled thermal-hydro-mechanical cohesive phase-field model for hydraulic fracturing in deep coal seams is presented.Heat exchange between the cold fluid and the hot rock is considered,and the thermal contribution terms between the cold fluid and the hot rock are derived.Heat transfer obeys Fourier's law,and porosity is used to relate the thermodynamic parameters of the fracture and matrix domains.The net pressure difference between the fracture and the matrix is neglected,and thus the fluid flow is modeled by the unified fluid-governing equations.The evolution equations of porosity and Biot's coefficient during hydraulic fracturing are derived from their definitions.The effect of coal cleats is considered and modeled by Voronoi polygons,and this approach is shown to have high accuracy.The accuracy of the proposed model is verified by two sets of fracturing experiments in multilayer coal seams.Subsequently,the differences in fracture morphology,fluid pressure response,and fluid pressure distribution between direct fracturing of coal seams and indirect fracturing of shale interlayers are explored,and the effects of the cluster number and cluster spacing on fracture morphology for multi-cluster fracturing are also examined.The numerical results show that the proposed model is expected to be a powerful tool for the fracturing design and optimization of deep coalbed methane.展开更多
基金Funded by the National Natural Science Foundation of China Academy of Engineering Physics and Jointly Setup"NSAF"Joint Fund(No.U1430119)。
文摘The multi-scale modeling combined with the cohesive zone model(CZM)and the molecular dynamics(MD)method were preformed to simulate the crack propagation in NiTi shape memory alloys(SMAs).The metallographic microscope and image processing technology were employed to achieve a quantitative grain size distribution of NiTi alloys so as to provide experimental data for molecular dynamics modeling at the atomic scale.Considering the size effect of molecular dynamics model on material properties,a reasonable modeling size was provided by taking into account three characteristic dimensions from the perspective of macro,meso,and micro scales according to the Buckinghamπtheorem.Then,the corresponding MD simulation on deformation and fracture behavior was investigated to derive a parameterized traction-separation(T-S)law,and then it was embedded into cohesive elements of finite element software.Thus,the crack propagation behavior in NiTi alloys was reproduced by the finite element method(FEM).The experimental results show that the predicted initiation fracture toughness is in good agreement with experimental data.In addition,it is found that the dynamics initiation fracture toughness increases with decreasing grain size and increasing loading velocity.
文摘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.
文摘The cohesive zone model(CZM)has been used widely and successfully in fracture propagation,but some basic problems are still to be solved.In this paper,artificial compliance and discontinuous force in CZM are investigated.First,theories about the cohesive element(local coordinate system,stiffness matrix,and internal nodal force)are presented.The local coordinate system is defined to obtain local separation;the stiffness matrix for an eight-node cohesive element is derived from the calculation of strain energy;internal nodal force between the cohesive element and bulk element is obtained from the principle of virtual work.Second,the reason for artificial compliance is explained by the effective stiffnesses of zero-thickness and finite-thickness cohesive elements.Based on the effective stiffness,artificial compliance can be completely removed by adjusting the stiffness of the finite-thickness cohesive element.This conclusion is verified from 1D and 3D simulations.Third,three damage evolution methods(monotonically increasing effective separation,damage factor,and both effective separation and damage factor)are analyzed.Under constant unloading and reloading conditions,the monotonically increasing damage factor method without discontinuous force and healing effect is a better choice than the other two methods.The proposed improvements are coded in LS-DYNA user-defined material,and a drop weight tear test verifies the improvements.
基金funded by the Scientific research startup Foundation of Fujian University of Technology(GY-Z21067 and GY-Z21026).
文摘Amid urbanization and the continuous expansion of transportation networks,the necessity for tunnel construction and maintenance has become paramount.Addressing this need requires the investigation of efficient,economical,and robust tunnel reinforcement techniques.This paper explores fiber reinforced polymer(FRP)and steel fiber reinforced concrete(SFRC)technologies,which have emerged as viable solutions for enhancing tunnel structures.FRP is celebrated for its lightweight and high-strength attributes,effectively augmenting load-bearing capacity and seismic resistance,while SFRC’s notable crack resistance and longevity potentially enhance the performance of tunnel segments.Nonetheless,current research predominantly focuses on experimental analysis,lacking comprehensive theoretical models.To bridge this gap,the cohesive zone model(CZM),which utilizes cohesive elements to characterize the potential fracture surfaces of concrete/SFRC,the rebar-concrete interface,and the FRP-concrete interface,was employed.A modeling approach was subsequently proposed to construct a tunnel segment model reinforced with either SFRC or FRP.Moreover,the corresponding mixed-mode constitutive models,considering interfacial friction,were integrated into the proposed model.Experimental validation and numerical simulations corroborated the accuracy of the proposed model.Additionally,this study examined the reinforcement design of tunnel segments.Through a numerical evaluation,the effectiveness of innovative reinforcement schemes,such as substituting concrete with SFRC and externally bonding FRP sheets,was assessed utilizing a case study from the Fuzhou Metro Shield Tunnel Construction Project.
基金supported by the National Natural Science Foundation of China(Grant Nos.11932006,12202314,12172121,and 12002118).
文摘A new cyclic cohesive zone fatigue damage model is proposed to address the fatigue problem spanning highand low cycle stages.The new damage model is integrated with the damage extrapolation technique to improvecalculation efficiency.The model’s effectiveness in regulating the low-cycle fatigue evolution rate,overall fatiguedamage evolution rate,and stress level at the fatigue turning point is assessed through the comparison of the S-Ncurves.The fatigue damage model’s high precision is proved based on the minor deviation of stress at the turningpoint of the S-N curve from the actual scenario.Finally,the fatigue damage evolution is simulated consideringthe effects of pre-load pressure and welding residual stress.It is observed that laser welding induces a significantresidual tensile stress,accelerating fatigue damage evolution,while compressive loading impedes fatigue damageprogression.
基金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.
文摘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.
文摘This paper concentrates on simulating fracture in thin walled single-lap joints connected by resistance spot-welding(RSW)process which were subjected to tensile loading.For this purpose,three sets of lap-joints with different spot configurations were tested to achieve the joints’tensile behavior.To simulate the joints tensile behavior,firstly a 2D axisymmetric finite element(FE)model was used to calculate residual stresses induced during the welding process.Then the results were transferred to 3D models as pre-stress.In this step,cohesive zone model(CZM)technique was used to simulate fracture in the models under tensile load.Cohesive zone parameters were extracted using coach-peel and shear lap specimens.The results were employed to simulate deformation and failure in single lap spot weld samples.It has been shown that considering the residual stresses in simulating deformation and fracture load enables quite accurate predictions.
基金The authors gratefully acknowledge the financial support of the National Natural Science Foundation of China(Nos.U1960205 and 51804024).
文摘The softening-melting characteristics of ferrous burden play a crucial role in the thickness and position of the cohesive zone.The influence of the basicity and experimental atmosphere on the softening-melting behavior of primary slag under slag-coke interaction was investigated using in situ visualization method.The mechanism was analyzed using the FactSage software,X-ray diffraction,and electron probe microanalysis.The softening and melting temperatures of the samples increased with increasing basicity under different atmospheres.The difference between softening and melting temperatures is smaller in a H_(2) atmosphere than in a CO atmosphere;in H_(2) atmosphere,the range of softening zone in the cohesive zone was significantly thinner.The formed low-melting-point FeO-bearing phases decrease when H_(2) was used as the reducing agent.In addition,according to FactSage calculations,the high content of metallic iron reduced in the H_(2) atmosphere raised the softening temperature of the primary slag.It also narrowed and moved downward the cohesive zone due to an increase in softening and melting temperatures.Meanwhile,the increase in basicity promoted the decrease in liquid ratio and improved the permeability of cohesive zone.
基金Item Sponsored by the Fundamental Research Funds for the Central Universities of China(N090402021)
文摘Based on the principle of discrete element method （DEM）, a 2D slot model of a COREX melter gasifier was established to analyze the influence of cohesive zone shape on solid flow, including mass distribution, velocity distribution, normal force distribution and porosity distribution at a microscopic level. The results show that the co- hesive zone shape almost does not affect the particle movement in the upper shaft and deadman shape. The particles in the lower central bottom experience large normal force to support the particles above them, while particles around the raceway and in the fast flow zone exhibit weak force network. The porosity distribution was also examined under three kinds of cohesive zones. Like the velocity distribution, the whole packed bed can be divided into four main re- gions. With the increase of cohesive zone position, the low porosity region located in the root of cohesive zone increa- ses. And the porosity distribution becomes asymmetric in the case of biased cohesive zone.
文摘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(51575535)supported by the National Natural Science Foundation of ChinaProject(2015CX002)supported by the Innovation-driven Plan in Central South University,China+2 种基金Project(zzyjkt2013-09B)supported by the Fund of the State Key Laboratory of High Performance Manufacturing,ChinaProject(2017zzts638)supported by the Fundamental Research Funds for the Central Universities,ChinaProject(2016RS2015)supported by the Scientific and Technological Leading Talent Projects of Hunan Province,China
文摘The key parameters of the adhesive layer of a reinforcing patch are of great significance and affect the ability to suppress crack propagation in an Al–Li alloy patch-reinforced structure.This paper proposes a method to determine the key parameters of the adhesive layer of adhesively bonded joints in the Al–Li alloy patch-reinforced structure.A zero-thickness cohesive zone model(CZM)was selected to simulate the adhesive layer’s fracture process,and an orthogonal simulation was designed to compare against the test results.A three-dimensional progressive damage model of an Al–Li alloy patch-reinforced structure with single-lap adhesively bonded joints was developed.The simulation’s results closely agree with the test results,demonstrating that this method of determining the key parameters is likely accurate.The results also verify the correctness of the cohesive strength and fracture energy,the two key parameters of the cohesive zone model.The model can accurately predict the strength and fracture process of adhesively bonded joints,and can be used in research to suppress crack propagation in Al–Li alloy patch-reinforced structures.
基金supported by the National Natural Science Foundation of China(No.11702147)。
文摘Investigation of paper cutting process is vital for the design of cutting tools,but the fracture mechanism of paper cutting is still unclear.Here,we focus on the cutting process of paper,including the key parameters of cohesive zone model(CZM)for the orthotropic paper,to simulate the shear fracture process.Firstly,the material constants of the orthotropic paper are determined by longitudinal and transverse tensile test.Secondly,based on the tensile stressstrain curves,combined with damage theory and numerical simulations,the key parameters of the CZM for the orthotropic paper are obtained.Finally,a model III fracture is simulated to verify the accuracy of the model.Results show that the load-displacement curves obtained by the simulation is consistent with the test results.
基金Project supported by the National Natural Science Foundation of China (No.50371042) the Post Doctoral Science Foundation of China (No.20040350031)
文摘Cohesive zone model was used to simulate two-dimensional plane strain crack propagation at the grain level model including grain boundary zones. Simulated results show that the original crack-tip may not be separated firstly in an elastic-viscoplastic polycrystals. The grain interior＇s material properties （e.g. strain rate sensitivity） characterize the competitions between plastic and cohesive energy dissipation mechanisms. The higher the strain rate sensitivity is, the larger amount of the external work is transformed into plastic dissipation energy than into cohesive energy, which delays the cohesive zone rupturing. With the strain rate sensitivity decreased, the material property tends to approach the elastic-plastic responses. In this case, the plastic dissipation energy decreases and the cohesive dissipation energy increases which accelerates the cohesive zones debonding. Increasing the cohesive strength or the critical separation displacement will reduce the stress triaxiality at grain interiors and grain boundaries. Enhancing the cohesive zones ductility can improve the matrix materials resistance to void damage.
基金the National Natural Science Foundation of China(No.50775142).
文摘An improved interface cohesive zone model is developed for the simulation of interface contact, under mixed-mode loading. A new debonding initiation criterion and propagation of debonding law, taking into account the pressure stress influence on contact shear strength, is proposed. The model is implemented in a finite-element program using subroutine VUINTER of ABAQUS Explicit. An edge-notch four-point bending process and laminated vibration damping steel sheet punch forming test are simulated with the improved model in ABAQUS Explicit. The numerical predictions agree satisfactorily with the corresponding experimental results.
基金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.
基金support of the National Natural Science Foundation of China(Nos.U1960205 and 51804024)China Baowu Low Carbon Metallurgy Innovation Foundation(BWLCF202101 and BWLCF202104)China Minmetals Science and Technology Special Plan Foundation(2020ZXA01).
文摘The changes in the softening and melting behaviors of ferrous burden in the cohesive zone and the characteristics of the slag–iron–coke interface in a blast furnace were investigated by simulating an actual blast furnace under hydrogen-rich conditions.According to the variation in the transient shrinkage of the burden under different atmospheres,the shrinkage start temperature of the sinter was higher than that of the pellets.The negative shrinkage rate of the pellets was greater than that of the sinter.Additionally,the softening start temperature in the blast furnace decreased under hydrogen-rich conditions,giving the blast furnace a broader range of softening zones.The softening start temperatures of the pellets and sinter decreased from 1102 to 949℃ and 1152 to 1080℃,respectively.The hydrogen-rich traditional blast furnace conditions narrowed the melting zone temperature range and shifted it toward the high-temperature zone,significantly improving the burden layer permeability.However,under the hydrogen-rich oxygen blast furnace conditions,there were a decrease in the melting start temperature,a shift of the melting zone location to the low-temperature zone,and an increase in the burden layer permeability and pressure difference.A comparison of the slag–iron–coke interface characteristics under different atmospheric conditions showed that the carbon content in metallic iron decreased under hydrogen-rich traditional blast furnace conditions compared with traditional blast furnace conditions.Contrastingly,under hydrogen-rich oxygen blast furnace conditions,the carbon content in metallic iron increased compared with oxygen blast furnace conditions.These findings provide guidance for the development of low-carbon ironmaking processes in blast furnaces.
基金The authors acknowledge the financial support from the National Natural Science Foundation of China(grant Nos.52274328 and 52074079)China Baowu Low Carbon Metallurgy Innovation Foundation(grant No.BWLCF202403)for this work.
文摘To reduce greenhouse gas emissions in ironmaking,the steel industry is advancing innovative low-carbon blast furnace(BF)technologies.A critical challenge for implementing such innovations lies in optimizing permeability within the BF's cohesive zone(CZ),which directly impacts operational stability and efficiency.This study employs a coupled computational fluid dynamics-discrete element method(CFD-DEM)to calibrate Young's modulus by respectively fitting the relationship between Young's modulus and temperature,as well as pressure drop,based on a reported lab-scale softening and smelting experimental data of ore-coke heterogeneous alternating layer packed beds resembling BFs,and develops a softening particle model(SPM).The SPM establishes a temperature-dependent relationship between mechanical properties of softened ore particles and CZ conditions in industrial-scale BFs.Simulations of particle shrinkage behavior and pressure drop trends using the CFD-DEM-SPM framework demonstrate strong correlation with experimental data,validating its accuracy for predictive analysis.Furthermore,this study investigates how layer arrangement configurations,size ratios between ore and coke particles,and coke blending proportions influence CZ characteristics.Key findings identify an optimal batch weight configuration to enhance permeability within the CZ while maintaining operational stability.Additionally,results indicate that increasing the relative particle size of ore compared to coke or enhancing the proportion of blended coke in burden mixes improves CZ permeability,offering actionable strategies for reducing carbon intensity in BF operations.These insights provide critical guidance toward developing low-carbon BF processes compatible with global climate targets.
基金the National Natural Science Foundation of China project(grant No.52264042,51904122)the China Postdoctoral Science Foundation Funding(grant No.2021M690975)+2 种基金the Jiangxi Provincial Natural Science Foundation(grant No.20212BDH81001,20212BAB214023,20223AAG01009,20214BBG74005)Opening Research Projects of State Key Laboratory of Advanced Metallurgy(grant No.K22-03&04)Research Project from Xiangtan Iron and Steel Corporation of Hunan Hualing and Nanchang 100-100 Plan Foundation for financial support to this work.
文摘Cohesive zone plays a vital role in the stable operation of a blast furnace(BF),yet the complex phase change process of iron ore particles in this zone is still not well understood.In this study,a novel one-dimensional(1D)unsteady phase change model was developed to elucidate the heat transfer and melting mechanisms of iron ore particles.After model validation,the effects of several key operating parameters(e.g.,particle diameter,gas velocity,initial temperature)on the phase change behavior of iron ore particles were analyzed,and the joint effect of multiple parameters was discussed.The results show that larger-sized iron ore particles possess lower specific surface areas,which in turn reduces their convective heat absorption capacity.Consequently,the distance from the solid-liquid phase interface to the particle surface increases,thereby slowing down the movement of the phase interface and pro-longing the melting duration of the particles.Increasing the gas velocity and the initial temperature does not have a significant impact on reducing the duration of the complete melting process.Under the specified conditions,it is observed that increasing the gas velocity by 3-fold and 9-fold results in a reduction of the melting duration by 2.4%and 8.3%,respectively.Elevating the initial temperature of iron ore particles results in a decrease in the core-to-surface temperature difference,a slower heating rate,and a shorter duration to achieve melting.Among the factors affecting the melting process,the particle diameter is found to be the most significant in terms of the liquid phase precipitation,mushy zone thickness,and core-to-surface temperature difference of iron ore particles.
基金Project supported by the National Natural Science Foundation of China(No.42202314)。
文摘A coupled thermal-hydro-mechanical cohesive phase-field model for hydraulic fracturing in deep coal seams is presented.Heat exchange between the cold fluid and the hot rock is considered,and the thermal contribution terms between the cold fluid and the hot rock are derived.Heat transfer obeys Fourier's law,and porosity is used to relate the thermodynamic parameters of the fracture and matrix domains.The net pressure difference between the fracture and the matrix is neglected,and thus the fluid flow is modeled by the unified fluid-governing equations.The evolution equations of porosity and Biot's coefficient during hydraulic fracturing are derived from their definitions.The effect of coal cleats is considered and modeled by Voronoi polygons,and this approach is shown to have high accuracy.The accuracy of the proposed model is verified by two sets of fracturing experiments in multilayer coal seams.Subsequently,the differences in fracture morphology,fluid pressure response,and fluid pressure distribution between direct fracturing of coal seams and indirect fracturing of shale interlayers are explored,and the effects of the cluster number and cluster spacing on fracture morphology for multi-cluster fracturing are also examined.The numerical results show that the proposed model is expected to be a powerful tool for the fracturing design and optimization of deep coalbed methane.
