Progressive fragmentation of granular assemblies within rockslides: Insights from discrete-continuous numerical modeling

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.

Impact Analysis of Microscopic Defect Types on the Macroscopic Crack Propagation in Sintered Silver Nanoparticles

Sintered silver nanoparticles(AgNPs)arewidely used in high-power electronics due to their exceptional properties.However,the material reliability is significantly affected by various microscopic defects.In this work,the three primary micro-defect types at potential stress concentrations in sintered AgNPs are identified,categorized,and quantified.Molecular dynamics(MD)simulations are employed to observe the failure evolution of different microscopic defects.The dominant mechanisms responsible for this evolution are dislocation nucleation and dislocation motion.At the same time,this paper clarifies the quantitative relationship between the tensile strain amount and the failure mechanism transitions of the three defect types by defining key strain points.The impact of defect types on the failure process is also discussed.Furthermore,traction-separation curves extracted from microscopic defect evolutions serve as a bridge to connect the macro-scale model.The validity of the crack propagation model is confirmed through tensile tests.Finally,we thoroughly analyze how micro-defect types influence macro-crack propagation and attempt to find supporting evidence from the MD model.Our findings provide a multi-perspective reference for the reliability analysis of sintered AgNPs.

Fractures interaction and propagation mechanism of multi-cluster fracturing on laminated shale oil reservoir

The continental shale reservoirs of Jurassic Lianggaoshan Formation in Sichuan Basin contain thin lamina,which is characterized by strong plasticity and developed longitudinal shell limestone interlayer.To improve the production efficiency of reservoirs by multi-cluster fracturing,it is necessary to consider the unbalanced propagation of hydraulic fractures and the penetration effect of fractures.This paper constructed a numerical model of multi-fracture propagation and penetration based on the finite element coupling cohesive zone method;considering the construction cluster spacing,pump rate,lamina strength and other parameters studied the influencing factors of multi-cluster fracture interaction propagation;combined with AE energy data and fracture mode reconstruction method,quantitatively characterized the comprehensive impact of the strength of thin interlayer rock interfaces on the initiation and propagation of fractures that penetrate layers,and accurately predicted the propagation pattern of hydraulic fractures through laminated shale oil reservoirs.Simulation results revealed that in the process of multi-cluster fracturing,the proportion of shear damage is low,and mainly occurs in bedding fractures activated by outer fractures.Reducing the cluster spacing enhances the fracture system's penetration ability,though it lowers the activation efficiency of lamina.The high plasticity of the limestone interlayer may impact the vertical propagation distance of the main fracture.Improving the interface strength is beneficial to the reconstruction of the fracture height,but the interface communication effect is limited.Reasonable selection of layers with moderate lamina strength for fracturing stimulation,increasing the pump rate during fracturing and setting the cluster spacing reasonably are beneficial to improve the effect of reservoir stimulation.

Improvements of cohesive zone model on artificial compliance and discontinuous force

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.

Finite Element Simulations of the Localized Failure and Fracture Propagation in Cohesive Materials with Friction

Strain localization frequently occurs in cohesive materials with friction(e.g.,composites,soils,rocks)and is widely recognized as a fundamental cause of progressive structural failure.Nonetheless,achieving high-fidelity simulation for this issue,particularly concerning strong discontinuities and tension-compression-shear behaviors within localized zones,remains significantly constrained.In response,this study introduces an integrated algorithmwithin the finite element framework,merging a coupled cohesive zone model(CZM)with the nonlinear augmented finite elementmethod(N-AFEM).The coupledCZMcomprehensively describes tension-compression and compressionshear failure behaviors in cohesive,frictional materials,while the N-AFEM allows nonlinear coupled intraelement discontinuities without necessitating extra nodes or nodal DoFs.Following CZM validation using existing experimental data,this integrated algorithm was utilized to analyze soil slope failure mechanisms involving a specific tensile strength and to assess the impact of mechanical parameters(e.g.,tensile strength,weighting factor,modulus)in soils.

Constitutive Behavior of the Interface between UHPC and Steel Plate without Shear Connector:From Experimental to Numerical Study

The application of ultra-high performance concrete(UHPC)as a covering layer for steel bridge decks has gained widespread popularity.By employing a connection without a shear connector between the steel plate and UHPC,namely,the sandblasted interface and the epoxy adhesive with sprinkled basalt aggregate interface,the installation cannot only be simplified but also the stress concentration resulting from the welded shear connectors can be eliminated.This study develops constitutive models for these two interfaces without shear connectors,based on the interfacial pull-off and push-out tests.For validation,three-point bending tests on the steel-UHPC composite plates are conducted.The results indicated that the proposed bilinear traction-separation model for the sandblasted interface and the trapezoidal traction-separation model for the epoxy adhesive with sprinkled basalt aggregate interface can generally calibrate the interfacial behavior.However,the utilization of the experimentally determined pure shear strength underestimates the load-carrying capacity of the composite plates in the case of three-point bending tests.By recalling the Mohr-Coulomb criterion,this underestimation is attributed to the enhancement of the interface shear strength by the presence of normal stress.

Effect of intermittent joint distribution on the mechanical and acoustic behavior of rock masses

The mechanical characteristics and acoustic behavior of rock masses are greatly influenced by stochastic joints.In this study,numerical models of rock masses incorporating intermittent joints with different numbers and dip angles were produced using the finite element method(FEM)with the intrinsic cohesive zone model(ICZM).Then,the uniaxial compressive and wave propagation simulations were performed.The results indicate that the joint number and dip angle can affect the mechanical and acoustic properties of the models.The uniaxial compressive strength(UCS)and wave velocity of rock masses decrease monotonically as the joint number increases.However,the wave velocity grows monotonically as the joint dip angle increases.When the joint dip angle is 45°–60°,the UCS of the rock mass is lower than that of other dip angles.The wave velocity parallel to the joints is greater than that perpendicular to the joints.When the dip angle of joints remains unchanged,the UCS and wave velocity are positively related.When the joint dip angle increases,the variation amplitude of the UCS regarding the wave velocity increases.To reveal the effect of the joint distribution on the velocity,a theoretical model was also proposed.According to the theoretical wave velocity,the change in wave velocity of models with various joint numbers and dip angles was consistent with the simulation results.Furthermore,a theoretical indicator(i.e.fabric tensor)was adopted to analyze the variation of the wave velocity and UCS.

Investigation of FRP and SFRC Technologies for Efficient Tunnel Reinforcement Using the Cohesive Zone Model

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.

Research on Fatigue Damage Behavior of Main Beam Sub-Structure of Composite Wind Turbine Blade

Given the difficulty in accurately evaluating the fatigue performance of large composite wind turbine blades(referred to as blades),this paper takes the main beam structure of the blade with a rectangular cross-sectionas the simulation object and establishes a composite laminate rectangular beam structure that simultaneouslyincludes the flange,web,and adhesive layer,referred to as the blade main beam sub-structure specimen,throughthe definition of blade sub-structures.This paper examines the progressive damage evolution law of the compositelaminate rectangular beam utilizing an improved 3D Hashin failure criterion,cohesive zone model,B-K failurecriterion,and computer simulation technology.Under static loading,the layup angle of the anti-shear web hasa close relationship with the static load-carrying capacity of the composite laminate rectangular beam;under fatigueloading,the fatigue damage will first occur in the lower flange adhesive area of the whole composite laminaterectangular beam and ultimately result in the fracture failure of the entire structure.These results provide a theoreticalreference and foundation for evaluating and predicting the fatigue performance of the blade main beamstructure and even the full-size blade.

COHESIVE ZONE FINITE ELEMENT-BASED MODELING OF HYDRAULIC FRACTURES

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.

Failure simulation in resistance spot-welded lap-joints using cohesive zone modeling

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.

Mesoscale Modeling of Hooked-End Steel Fiber Reinforced Concrete under Uniaxial Compression Using Cohesive Elements

Based on the cohesive zone model, the 2D mesostructures were developed for numerical studies of multi-phase hooked-end steel fiber reinforced concrete under uniaxial compression. The zero-thickness cohesive interface elements were inserted within the mortar, on interfaces of mortar and aggregates and interfaces of mortar and fibers to simulate the failure process of fiber reinforced concrete. The results showed that the numerical results matched well the experimental results in both failure modes and stress-strain behavior. Hooked-end steel fiber reinforced concrete exhibited ductile failure and maintained integrity during a whole failure process. Compared with normal concrete, HES fiber reinforced concrete was greater stiffness and compressive strength;the descending branch of the stress-strain curve was significantly flatter;the residual stress was higher.

Numerical simulation of the mechanical behavior of superconducting tape in conductor on round core cable using the cohesive zone model

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.

Numerical study of fatigue damage of asphalt concrete using cohesive zone model

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.

INTERFACE DAMAGE ANALYSIS OF FIBER REINFORCED COMPOSITES WITH DUCTILE MATRIX

A cohesive zone model is employed to simulate the fiber/matrix interface damage of composites with ductile matrix. The study is carried out to investigate the dependence of the interface damage and the composite tensile strength on the micro parameters of the composite. These parameters contain fiber packing pattern, fiber volume fraction, and the modulus ratio of the fiber to the matrix. The investigation reveals that though the high fiber vo lume fraction, the high fiber′s modulus and the square fiber packing can supply strong reinforcement to the composite, the interface damage is susceptible in these cases. The tensile strength of the composite is dominated by the interface strength when the interface debonding occurs.

Analysis on the Forecast Characteristics of K Index and Low-altitude Jet Stream in the Rainstorm Process

By using the durative rainstorm data in South China during May-early June in 2010,the forecast characteristics of K index and low level jet were analyzed.The results found that K2 had the good indication,advancement and relativity on the intensity and falling zone forecast of regional rainstorm in future 24 h,and the positive relative coefficient reached 0.987.The low level jet also had the same advancement and indication significance on the intensity and influence scope of regional rainstorm in 24 h in the future,and the relative coefficient reached above 0.8.K2 and the low level jet were selected as the main factors,and the basic conceptual model of rainstorm falling zone was established.The model has passed the computer program and realized the business automation.K2 provided the important basis for the forecast of rainstorm intensity and falling zone.

Numerical Simulation of Particle/Matrix Interface Failure in Composite Propellant

Interface debonding between particle and matrix in composite propellant influences its macroscopic mechanical properties greatly. For this, the laws of interface cohesive damage and failure were analyzed. Then, its microscopic computational model was established. The interface mechanical response was modeled by the bilinear cohesive zone model. The effects of interface properties and particle sizes on the macroscopic mechanical behavior were investigated. Numerical simulation of debonding damage evolution of composite propellant under finite deformation was carried out. The debonding damage nucleation, propagation mechanism and non-uniform distribution of microscopic stress-strain fields were discussed. The results show that the finite element simulation method based on microstructure model can effectively predict the trend of macroscopic mechanical behavior and particle/matrix debonding evolution process. It can be used for damage simulation and failure assessment of composite propellants.

Mechanical characteristic variation of ballastless track in highspeed railway:effect of train–track interaction and environment loads

Due to the fact that ballastless tracks in highspeed railways are not only subjected to repeated train–track dynamic interaction loads,but also suffer from complex environmental loads,the fundamental understanding of mechanical performance of ballastless tracks under sophisticated service conditions is an increasingly demanding and challenging issue in high-speed railway networks.This work aims to reveal the effect of train–track interaction and environment loads on the mechanical characteristic variation of ballastless tracks in high-speed railways,particularly focusing on the typical interface damage evolution between track layers.To this end,a finite element model of a double-block ballastless track involving the cohesive zone model for the track interface is first established to analyze the mechanical properties of the track interface under the loading–unloading processes of the negative temperature gradient load(TGL)followed by the same cycle of the positive TGL.Subsequently,the effect of wheel–rail longitudinal interactions on the nonlinear dynamic characteristics of the track interface is investigated by using a vehicle-slab track vertical-longitudinal coupled dynamics model.Finally,the influence of dynamic water pressure induced by vehicle dynamic load on the mechanical characteristics and damage evolution of the track interface is elucidated using a fluid–solid coupling method.Results show that the loading history of the positive and negative TGLs has a great impact on the nonlinear development and distribution of the track interface stress and damage;the interface damage could be induced by the wheel–rail longitudinal vibrations at a high vehicle running speed owing to the dynamic amplification effect caused by short wave irregularities;the vehicle dynamic load could produce considerable water pressure that presents nonlinear spatial–temporal characteristics at the track interface,which would lead to the interface failure under a certain condition due to the coupled dynamic effect of vehicle load and water pressure.

Delamination analysis of woven fabrication laminates using cohesive zone model

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.

Analysis of pre-crack kinking and delamination in terms of SLB test of honeycomb sandwich structure

The flatwise tension(FWT) and single leg bending(SLB) tests were used to investigate the fracture behavior of honeycomb sandwich specimens.In the FWT test,only the interlaminar delamination was observed.The test results show that the interfacial peel strength is higher than the interlaminar peel strength.In terms of SLB experiment,a new fracture mode was found,namely IKP(initiation of interlaminar delamination,kinking into facesheet and propagation of interlaminar delamination).The virtual crack closure technique was applied to separate the values of Mode I and Mode II components of the strain energy release rate in SLB experiment.The finite element analysis result shows that the Mode I strain energy release rate is higher than the Mode II strain energy release rate.To simulate the failure of SLB test of honeycomb sandwich specimens,a new computational model based on the Tsai-Hill failure criterion and the cohesive zone model is proposed.In comparison with experimental results,it can be concluded that the computational model can validly simulate the IKP of a honeycomb sandwich structures with reasonable accuracy.