Rock residual strength,as an important input parameter,plays an indispensable role in proposing the reasonable and scientific scheme about stope design,underground tunnel excavation and stability evaluation of deep ch...Rock residual strength,as an important input parameter,plays an indispensable role in proposing the reasonable and scientific scheme about stope design,underground tunnel excavation and stability evaluation of deep chambers.Therefore,previous residual strength models of rocks established were reviewed.And corresponding related problems were stated.Subsequently,starting from the effects of bedding and whole life-cycle evolution process,series of triaxial mechanical tests of deep bedded sandstone with five bedding angles were conducted under different confining pressures.Then,six residual strength models considering the effects of bedding and whole life-cycle evolution process were established and evaluated.Finally,a cohesion loss model for determining residual strength of deep bedded sandstone was verified.The results showed that the effects of bedding and whole life-cycle evolution process had both significant influences on the evolution characteristic of residual strength of deep bedded sandstone.Additionally,residual strength parameters:residual cohesion and residual internal friction angle of deep bedded sandstone were not constant,which both significantly changed with increasing bedding angle.Besides,the cohesion loss model was the most suitable for determining and estimating the residual strength of bedded rocks,which could provide more accurate theoretical guidance for the stability control of deep chambers.展开更多
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.展开更多
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.展开更多
Rock fragmentation plays a critical role in rock avalanches,yet conventional approaches such as classical granular flow models or the bonded particle model have limitations in accurately characterizing the progressive...Rock fragmentation plays a critical role in rock avalanches,yet conventional approaches such as classical granular flow models or the bonded particle model have limitations in accurately characterizing the progressive disintegration and kinematics of multi-deformable rock blocks during rockslides.The present study proposes a discrete-continuous numerical model,based on a cohesive zone model,to explicitly incorporate the progressive fragmentation and intricate interparticle interactions inherent in rockslides.Breakable rock granular assemblies are released along an inclined plane and flow onto a horizontal plane.The numerical scenarios are established to incorporate variations in slope angle,initial height,friction coefficient,and particle number.The evolutions of fragmentation,kinematic,runout and depositional characteristics are quantitatively analyzed and compared with experimental and field data.A positive linear relationship between the equivalent friction coefficient and the apparent friction coefficient is identified.In general,the granular mass predominantly exhibits characteristics of a dense granular flow,with the Savage number exhibiting a decreasing trend as the volume of mass increases.The process of particle breakage gradually occurs in a bottom-up manner,leading to a significant increase in the angular velocities of the rock blocks with increasing depth.The simulation results reproduce the field observations of inverse grading and source stratigraphy preservation in the deposit.We propose a disintegration index that incorporates factors such as drop height,rock mass volume,and rock strength.Our findings demonstrate a consistent linear relationship between this index and the fragmentation degree in all tested scenarios.展开更多
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.展开更多
In reinforced concrete structures,corrosion of the rebar produces 2–6 times more corrosion product than the original material,creating pressure on the surrounding concrete,leading to cracking.The study of corrosion a...In reinforced concrete structures,corrosion of the rebar produces 2–6 times more corrosion product than the original material,creating pressure on the surrounding concrete,leading to cracking.The study of corrosion and cracking in reinforced concrete structures is therefore of great importance for enhancing the durability of concrete.Unlike many previous studies,we used ribbed rebar similar to that used commercially and considered the mechanical behavior of the interface transition zone(ITZ)between the aggregate and mortar to simulate the processes of corrosion and cracking of reinforced concrete structures.We explored the failure mode of the interface layer under uniform corrosion and the influence of different factors on the corrosion expansion cracking and the shedding mode of a concrete cover.This was achieved by establishing a three-phase meso-scale model of concrete based on secondary development of ABAQUS,simulating the mechanical behavior of the ITZ using a cohesive element,and establishing a rust expansion cracking model for single and multiple rebars.The results showed that:(1)Under uniform rust expansion,concrete cracks are distributed in a cross pattern with a slightly shorter lower limb.(2)When the corrosion rate is low,the ITZ is not damaged.With an increase in the corrosion rate,the proportion of elements with tensile damage in the ITZ first increases and then decreases.(3)In the case of a single rebar,the larger the cover thickness,the higher the corrosion rate corresponding to ITZ failure,and the arrangement of the rebar has little influence on the ITZ failure mode.(4)In the case of multiple rebars,the concrete cover cracks when the rebar spacing is small,and wedge-shaped spalling occurs when the spacing is large.展开更多
This study investigates the bond performance at the interfacial region shared by Ultra-High Performance Concrete(UHPC)and steel tubes through push-out tests.This study examines how changes in steel fiber volumetric ra...This study investigates the bond performance at the interfacial region shared by Ultra-High Performance Concrete(UHPC)and steel tubes through push-out tests.This study examines how changes in steel fiber volumetric ratio and thickness of steel tube influence the bond strength characteristics.The results show that as the enhancement of the steel tube wall thickness,the ultimate bond strength at the interface improves significantly,whereas the initial bond strength exhibits only slight variations.The influence of steel fiber volumetric ratio presents a nonlinear trend,with initial bond strength decreasing at low fiber content and increasing significantly as fiber content rises.Additionally,finite element(FE)simulations were applied to replicate the experimental conditions,and the outcomes showed strong correlation with the experimental data,confirming the exactitude of the FE model in predicting the bond behavior at the UHPC-Steel interface.These findings provide valuable insights for optimizing the design of UHPC-Filled steel tubes in high-performance structure.展开更多
A comprehensive numerical investigation into mixed⁃mode delamination is presented in this study.It aims to assess the impact of thermal and friction effects through mixed⁃mode flexure crack propagation testing.Finite ...A comprehensive numerical investigation into mixed⁃mode delamination is presented in this study.It aims to assess the impact of thermal and friction effects through mixed⁃mode flexure crack propagation testing.Finite element analysis was employed to model the delamination process,incorporating a contact cohesive zone model.This model couples the traction⁃separation law,the contact law,and the Coulomb friction law simultaneously.The thermomechanical analysis in this study is performed using a sequentially coupled approach,implemented with the finite element software ABAQUS.The findings underscore the importance of this study.展开更多
The development of geological lamination in shale reservoirs influences fracture propagation during hydraulic stimulation,and the fracture generation mechanism as it propagates through the laminated interface is close...The development of geological lamination in shale reservoirs influences fracture propagation during hydraulic stimulation,and the fracture generation mechanism as it propagates through the laminated interface is closely related to fracturing effects.In this paper,the laminated shale was selected to conduct three-point bending experiments using digital image correlation(DIC)and acoustic emission(AE)techniques,which revealed that the propagation path of cross-layer fractures exhibits dislocation features.The cohesive fracture mechanism of cross-layer fractures is investigated from the viewpoint of the fracture process zone(FPZ),which displays the characteristics of intermittence and dislocation during fracture development.A computational criterion for predicting the dislocation of cross-layer fracture at the interface is proposed,which shows that the maximum dislocation range does not exceed 72%of the FPZ length.Considering the mechanical differences between adjacent layers of laminated shale,the cohesive zone model of cross-layer fracture is discussed,from which the constitutive relationship and fracture energy during FPZ development are characterized,and the discontinuous nature of the constitutive relationship is found.This study improves the understanding of the geometry and cohesive fracture mechanism of the cross-layer fracture and provides valuable insights for field fracturing in shale reservoirs.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
The objective of this paper is to develop a methodology for calibration of a discrete element grain-based model(GBM)to replicate the hydro-mechanical properties of a brittle rock measured in the laboratory,and to appl...The objective of this paper is to develop a methodology for calibration of a discrete element grain-based model(GBM)to replicate the hydro-mechanical properties of a brittle rock measured in the laboratory,and to apply the calibrated model to simulating the formation of excavation damage zone(EDZ)around underground excavations.Firstly,a new cohesive crack model is implemented into the universal distinct element code(UDEC)to control the fracturing behaviour of materials under various loading modes.Next,a methodology for calibration of the components of the UDEC-Voronoi model is discussed.The role of connectivity of induced microcracks on increasing the permeability of laboratory-scale samples is investigated.The calibrated samples are used to investigate the influence of pore fluid pressure on weakening the drained strength of the laboratory-scale rock.The validity of the Terzaghi’s effective stress law for the drained peak strength of low-porosity rock is tested by performing a series of biaxial compression test simulations.Finally,the evolution of damage and pore pressure around two unsupported circular tunnels in crystalline granitic rock is studied.展开更多
Models to describe the damage and fracture behaviors of the interface between the fuel foil and cladding in UMo/Zr monolithic fuel plates were established and numerically implemented.The effects of the interfacial coh...Models to describe the damage and fracture behaviors of the interface between the fuel foil and cladding in UMo/Zr monolithic fuel plates were established and numerically implemented.The effects of the interfacial cohesive strength and cohesive energy on the irradiationinduced thermal-mechanical behaviors of fuel plates were investigated.The results indicated that for heterogeneously irradiated fuel plates:(1)interfacial damage and failure were predicted to be initiated near the fuel foil corner with higher fission densities,accompanied by the formation of a large gap after interface failure,which was consistent with some experimental observations;high tensile stresses in the fuel foil occurred near the edges of the failed interface,attributed to through-thickness cracking of the fuel foil,as found in some post-irradiation examinations;(2)the cohesive strength and cohesive energy of the interface both influenced the in-pile evolution behaviors of fuel plates;a lower cohesive strength or cohesive energy resulted in faster interfacial damage;(3)after interface fracture,the thickness of the whole plate increased to a greater degree(by~20%)than that of the samples without interfacial damage,which was attributed to the locally enhanced Mises stresses and the nearby creep deformations around the cracked interface.This study provided a theoretical basis for assessing failure in fuel elements.展开更多
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.展开更多
基金Projects(2024YFC3013801,2022YFC3004602)supported by the National Key R&D Program of ChinaProjects(U23B2093,52034009)supported by the National Natural Science Foundation of China。
文摘Rock residual strength,as an important input parameter,plays an indispensable role in proposing the reasonable and scientific scheme about stope design,underground tunnel excavation and stability evaluation of deep chambers.Therefore,previous residual strength models of rocks established were reviewed.And corresponding related problems were stated.Subsequently,starting from the effects of bedding and whole life-cycle evolution process,series of triaxial mechanical tests of deep bedded sandstone with five bedding angles were conducted under different confining pressures.Then,six residual strength models considering the effects of bedding and whole life-cycle evolution process were established and evaluated.Finally,a cohesion loss model for determining residual strength of deep bedded sandstone was verified.The results showed that the effects of bedding and whole life-cycle evolution process had both significant influences on the evolution characteristic of residual strength of deep bedded sandstone.Additionally,residual strength parameters:residual cohesion and residual internal friction angle of deep bedded sandstone were not constant,which both significantly changed with increasing bedding angle.Besides,the cohesion loss model was the most suitable for determining and estimating the residual strength of bedded rocks,which could provide more accurate theoretical guidance for the stability control of deep chambers.
基金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.
基金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.
基金support from the National Key R&D plan(Grant No.2022YFC3004303)the National Natural Science Foundation of China(Grant No.42107161)+3 种基金the State Key Laboratory of Hydroscience and Hydraulic Engineering(Grant No.2021-KY-04)the Open Research Fund Program of State Key Laboratory of Hydroscience and Engineering(sklhse-2023-C-01)the Open Research Fund Program of Key Laboratory of the Hydrosphere of the Ministry of Water Resources(mklhs-2023-04)the China Three Gorges Corporation(XLD/2117).
文摘Rock fragmentation plays a critical role in rock avalanches,yet conventional approaches such as classical granular flow models or the bonded particle model have limitations in accurately characterizing the progressive disintegration and kinematics of multi-deformable rock blocks during rockslides.The present study proposes a discrete-continuous numerical model,based on a cohesive zone model,to explicitly incorporate the progressive fragmentation and intricate interparticle interactions inherent in rockslides.Breakable rock granular assemblies are released along an inclined plane and flow onto a horizontal plane.The numerical scenarios are established to incorporate variations in slope angle,initial height,friction coefficient,and particle number.The evolutions of fragmentation,kinematic,runout and depositional characteristics are quantitatively analyzed and compared with experimental and field data.A positive linear relationship between the equivalent friction coefficient and the apparent friction coefficient is identified.In general,the granular mass predominantly exhibits characteristics of a dense granular flow,with the Savage number exhibiting a decreasing trend as the volume of mass increases.The process of particle breakage gradually occurs in a bottom-up manner,leading to a significant increase in the angular velocities of the rock blocks with increasing depth.The simulation results reproduce the field observations of inverse grading and source stratigraphy preservation in the deposit.We propose a disintegration index that incorporates factors such as drop height,rock mass volume,and rock strength.Our findings demonstrate a consistent linear relationship between this index and the fragmentation degree in all tested scenarios.
文摘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.
基金the National Natural Science Foundation of China(Nos.U1934213 and 51878572)。
文摘In reinforced concrete structures,corrosion of the rebar produces 2–6 times more corrosion product than the original material,creating pressure on the surrounding concrete,leading to cracking.The study of corrosion and cracking in reinforced concrete structures is therefore of great importance for enhancing the durability of concrete.Unlike many previous studies,we used ribbed rebar similar to that used commercially and considered the mechanical behavior of the interface transition zone(ITZ)between the aggregate and mortar to simulate the processes of corrosion and cracking of reinforced concrete structures.We explored the failure mode of the interface layer under uniform corrosion and the influence of different factors on the corrosion expansion cracking and the shedding mode of a concrete cover.This was achieved by establishing a three-phase meso-scale model of concrete based on secondary development of ABAQUS,simulating the mechanical behavior of the ITZ using a cohesive element,and establishing a rust expansion cracking model for single and multiple rebars.The results showed that:(1)Under uniform rust expansion,concrete cracks are distributed in a cross pattern with a slightly shorter lower limb.(2)When the corrosion rate is low,the ITZ is not damaged.With an increase in the corrosion rate,the proportion of elements with tensile damage in the ITZ first increases and then decreases.(3)In the case of a single rebar,the larger the cover thickness,the higher the corrosion rate corresponding to ITZ failure,and the arrangement of the rebar has little influence on the ITZ failure mode.(4)In the case of multiple rebars,the concrete cover cracks when the rebar spacing is small,and wedge-shaped spalling occurs when the spacing is large.
基金supported by grants from the Natural Science Foundation of Fujian Province(2021J011062)Minjiang Scholars Funding(GY-633Z21067).
文摘This study investigates the bond performance at the interfacial region shared by Ultra-High Performance Concrete(UHPC)and steel tubes through push-out tests.This study examines how changes in steel fiber volumetric ratio and thickness of steel tube influence the bond strength characteristics.The results show that as the enhancement of the steel tube wall thickness,the ultimate bond strength at the interface improves significantly,whereas the initial bond strength exhibits only slight variations.The influence of steel fiber volumetric ratio presents a nonlinear trend,with initial bond strength decreasing at low fiber content and increasing significantly as fiber content rises.Additionally,finite element(FE)simulations were applied to replicate the experimental conditions,and the outcomes showed strong correlation with the experimental data,confirming the exactitude of the FE model in predicting the bond behavior at the UHPC-Steel interface.These findings provide valuable insights for optimizing the design of UHPC-Filled steel tubes in high-performance structure.
文摘A comprehensive numerical investigation into mixed⁃mode delamination is presented in this study.It aims to assess the impact of thermal and friction effects through mixed⁃mode flexure crack propagation testing.Finite element analysis was employed to model the delamination process,incorporating a contact cohesive zone model.This model couples the traction⁃separation law,the contact law,and the Coulomb friction law simultaneously.The thermomechanical analysis in this study is performed using a sequentially coupled approach,implemented with the finite element software ABAQUS.The findings underscore the importance of this study.
基金financiallysupported by the Excellent Young Fund of Sinopec Petroleum Exploration and Production Research Institute(Grant No.YK2024009)the National Natural Science Foundation of China(Grant Nos.U23B6004 and 51925405).
文摘The development of geological lamination in shale reservoirs influences fracture propagation during hydraulic stimulation,and the fracture generation mechanism as it propagates through the laminated interface is closely related to fracturing effects.In this paper,the laminated shale was selected to conduct three-point bending experiments using digital image correlation(DIC)and acoustic emission(AE)techniques,which revealed that the propagation path of cross-layer fractures exhibits dislocation features.The cohesive fracture mechanism of cross-layer fractures is investigated from the viewpoint of the fracture process zone(FPZ),which displays the characteristics of intermittence and dislocation during fracture development.A computational criterion for predicting the dislocation of cross-layer fracture at the interface is proposed,which shows that the maximum dislocation range does not exceed 72%of the FPZ length.Considering the mechanical differences between adjacent layers of laminated shale,the cohesive zone model of cross-layer fracture is discussed,from which the constitutive relationship and fracture energy during FPZ development are characterized,and the discontinuous nature of the constitutive relationship is found.This study improves the understanding of the geometry and cohesive fracture mechanism of the cross-layer fracture and provides valuable insights for field fracturing in shale reservoirs.
文摘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.
文摘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.
基金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(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.
文摘The objective of this paper is to develop a methodology for calibration of a discrete element grain-based model(GBM)to replicate the hydro-mechanical properties of a brittle rock measured in the laboratory,and to apply the calibrated model to simulating the formation of excavation damage zone(EDZ)around underground excavations.Firstly,a new cohesive crack model is implemented into the universal distinct element code(UDEC)to control the fracturing behaviour of materials under various loading modes.Next,a methodology for calibration of the components of the UDEC-Voronoi model is discussed.The role of connectivity of induced microcracks on increasing the permeability of laboratory-scale samples is investigated.The calibrated samples are used to investigate the influence of pore fluid pressure on weakening the drained strength of the laboratory-scale rock.The validity of the Terzaghi’s effective stress law for the drained peak strength of low-porosity rock is tested by performing a series of biaxial compression test simulations.Finally,the evolution of damage and pore pressure around two unsupported circular tunnels in crystalline granitic rock is studied.
基金supported by the National Natural Science Foundation of China(Nos.12132005,11772095,12102094)the Foundation from Science and Technology on Reactor System Design Technology Laboratory。
文摘Models to describe the damage and fracture behaviors of the interface between the fuel foil and cladding in UMo/Zr monolithic fuel plates were established and numerically implemented.The effects of the interfacial cohesive strength and cohesive energy on the irradiationinduced thermal-mechanical behaviors of fuel plates were investigated.The results indicated that for heterogeneously irradiated fuel plates:(1)interfacial damage and failure were predicted to be initiated near the fuel foil corner with higher fission densities,accompanied by the formation of a large gap after interface failure,which was consistent with some experimental observations;high tensile stresses in the fuel foil occurred near the edges of the failed interface,attributed to through-thickness cracking of the fuel foil,as found in some post-irradiation examinations;(2)the cohesive strength and cohesive energy of the interface both influenced the in-pile evolution behaviors of fuel plates;a lower cohesive strength or cohesive energy resulted in faster interfacial damage;(3)after interface fracture,the thickness of the whole plate increased to a greater degree(by~20%)than that of the samples without interfacial damage,which was attributed to the locally enhanced Mises stresses and the nearby creep deformations around the cracked interface.This study provided a theoretical basis for assessing failure in fuel elements.
基金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.
