This study demonstrates a homogenization approach via a modified state-based peridynamic(PD)method to predict the effective elastic properties of composite materials with periodic microstructure.The procedure of model...This study demonstrates a homogenization approach via a modified state-based peridynamic(PD)method to predict the effective elastic properties of composite materials with periodic microstructure.The procedure of modeling the PD unit cell(UC)of continuous fiber-reinforced composite is presented.Periodic boundary conditions are derived and implemented through the Lagrange multiplier method.A matrix-dominated approach for modeling the interphase properties between dissimilar materials is proposed.The periodicity and continuity assumptions are employed to determine the stress and strain fields,as well as the effective elastic properties.The PD-UCs of square and hexagonal packs as well as the 0/90 laminate microstructure are modeled and compared with the analytical,numerical and experimental results from the literature.Good agreement of predicted effective properties can be observed.Unlike other PD homogenization approaches,the effective material properties can be directly and individually obtained from simple loading conditions.展开更多
The peridynamics(PD),as a promising nonlocal continuum mechanics theory,shines in solving discontinuous problems.Up to now,various numerical methods,such as the peridynamic mesh-free particlemethod(PD-MPM),peridynamic...The peridynamics(PD),as a promising nonlocal continuum mechanics theory,shines in solving discontinuous problems.Up to now,various numerical methods,such as the peridynamic mesh-free particlemethod(PD-MPM),peridynamic finite element method(PD-FEM),and peridynamic boundary element method(PD-BEM),have been proposed.PD-BEM,in particular,outperforms other methods by eliminating spurious boundary softening,efficiently handling infinite problems,and ensuring high computational accuracy.However,the existing PD-BEM is constructed exclusively for bond-based peridynamics(BBPD)with fixed Poisson’s ratio,limiting its applicability to crack propagation problems and scenarios involving infinite or semi-infinite problems.In this paper,we address these limitations by introducing the boundary element method(BEM)for ordinary state-based peridynamics(OSPD-BEM).Additionally,we present a crack propagationmodel embeddedwithin the framework ofOSPD-BEM to simulate crack propagations.To validate the effectiveness of OSPD-BEM,we conduct four numerical examples:deformation under uniaxial loading,crack initiation in a double-notched specimen,wedge-splitting test,and threepoint bending test.The results demonstrate the accuracy and efficiency of OSPD-BEM,highlighting its capability to successfully eliminate spurious boundary softening phenomena under varying Poisson’s ratios.Moreover,OSPDBEMsignificantly reduces computational time and exhibits greater consistencywith experimental results compared to PD-MPM.展开更多
In rock engineering,natural cracks in rock masses subjected to external loads tend to initiate and propagate,leading to potential safety hazards.To investigate the effect of cracking behavior on the mechanical propert...In rock engineering,natural cracks in rock masses subjected to external loads tend to initiate and propagate,leading to potential safety hazards.To investigate the effect of cracking behavior on the mechanical properties of rocks,the cracking processes of pre-cracked rocks have been extensively studied using numerical modeling methods.The peridynamics(PD)exhibits advantages over other numerical methods due to the absence of the requirements for remeshing and external crack growth criterion.However,for modeling pre-cracked rock cracking processes under impact,current PD implementations lack generally applicable rock constitutive models and impact contact models,which leads to difficulties in determining rock material parameters and efficiently calculating impact loads.This paper proposes a non-ordinary state-based peridynamics(NOSBPD)modeling method integrating the Drucker-Prager(DP)plasticity model and an efficient contact model to address the above problems.In the proposed method,the Drucker-Prager plasticity model is integrated into the NOSBPD,thereby equipping NOSBPD with the capability to accurately characterize the nonlinear stress-strain relationship inherent in rocks.An efficient contact model between particles and meshes is designed to calculate the impact loads,which is essentially a coupling method of PD with the finite element method(FEM).The effectiveness of the proposed NOSBPD modeling method is verified by comparison with other numerical methods and experiments.Experimental results indicate that the proposed method can effectively and accurately predict the 3D cracking processes of pre-cracked cracks under impact loading,and the maximum principal stress is the key driver behind wing crack formation in pre-cracked rocks.展开更多
The non-ordinary state-based peridynamic(NOSB PD)model has the capability of incorporating existing constitutive relationships in the classical continuum mechanics.In the present work,we first develop an NOSB PD model...The non-ordinary state-based peridynamic(NOSB PD)model has the capability of incorporating existing constitutive relationships in the classical continuum mechanics.In the present work,we first develop an NOSB PD model corresponding to the Johnson–Holmquist II(JH-2)constitutive damage model,which can describe the severe damage of concrete under intense impact compression.Besides,the numerical oscillation problem of the NOSB PD caused by zero-energy mode is analyzed and hence a bond-associated non-ordinary state-based peridynamic(BA-NOSB PD)model is adopted to remove the oscillation.Then,the elastic deformation of a three-dimensional bar is analyzed to verify the capability of BA-NOSB PD in eliminating the numerical oscillation.Furthermore,concrete spalling caused by the interaction of incident compression wave and reflected tension wave is simulated.The dynamic tensile fracture process of concrete multiple spalling is accurately reproduced for several examples according to the spalling number and spalling thickness analysis,illustrating the approach can well simulate and analyze the concrete spalling discontinuities.展开更多
Strain hardening and strain rate play an important role in dynamic deformation and failure problems such as high-velocity impact cases.In this paper,a non-ordinary state-based peridynamic model for failure and damage ...Strain hardening and strain rate play an important role in dynamic deformation and failure problems such as high-velocity impact cases.In this paper,a non-ordinary state-based peridynamic model for failure and damage of concrete materials subjected to impacting condition is proposed,taking the advantages of both damage model and nonlocal peridynamic method.The Holmquist-Johnson-Cook(HJC)model describing the mechanical character and damage of concrete materials under large strain,high strain rate and high hydrostatic pressure was reformulated in the framework of non-ordinary statebased peridynamic theory,and the corresponding numerical approach was developed.The proposed model and numerical approach were validated through simulating typical impacting examples and comparing the results with available experimental observations and results in literature.展开更多
In this study,a new state-based peridynamic formulation is developed for functionally graded Euler-Bernoulli beams.The equation of motion is developed by using Lagrange’s equation and Taylor series.Both axial and tra...In this study,a new state-based peridynamic formulation is developed for functionally graded Euler-Bernoulli beams.The equation of motion is developed by using Lagrange’s equation and Taylor series.Both axial and transverse displacements are taken into account as degrees of freedom.Four different boundary conditions are considered including pinned support-roller support,pinned support-pinned support,clamped-clamped and clamped-free.Peridynamic results are compared against finite element analysis results for transverse and axial deformations and a very good agreement is observed for all different types of boundary conditions.展开更多
In this study,a recently developed peridynamic lattice model called the“Statebased Peridynamic Lattice Model”(SPLM)is improved and demonstrated.In the SPLM,rather than as a continuum,solids are simulated using a clo...In this study,a recently developed peridynamic lattice model called the“Statebased Peridynamic Lattice Model”(SPLM)is improved and demonstrated.In the SPLM,rather than as a continuum,solids are simulated using a close-packed lattice of peridynamically interacting particles.The new SPLM approach advances the SPLM model by improving the damage and plasticity models.Elasticity,plasticity and damage are coupled in this approach.A robust method for damage initiation is developed.A new damage model called the“two-spring damage model”allows damage to localize to a single lattice particle,thus allowing highly localized damage(cracks)to emerge in a realistic manner.A plasticity model that includes hardening,softening,and damage due to plasticity is proposed and demonstrated.Peridynamic boundary effects are modeled efficiently and reasonably.The improved SPLM method is then employed to simulate three common concrete laboratory tests:uniaxial tension,uniaxial compression,and the Brazilian split cylinder test.The SPLM results are then compared with results from the earlier SPLM model,with simplified classical predictions,and with laboratory results.By solving the same benchmark problems using various lattice rotations and lattice spacings,the approach is demonstrated to be sufficiently objective to be a useful engineering tool to predict the essentially random behavior of concrete laboratory specimens.The improved SPLM demonstrates significant improvements over the previously published version and is found to simulate concrete structures accurately and efficiently using far less computational effort than comparable computational simulation methods.展开更多
We establish the a priori convergence rate for finite element approximations of a class of nonlocal nonlinear fracture models.We consider state-based peridynamic models where the force at a material point is due to bo...We establish the a priori convergence rate for finite element approximations of a class of nonlocal nonlinear fracture models.We consider state-based peridynamic models where the force at a material point is due to both the strain between two points and the change in volume inside the domain of the nonlocal interaction.The pairwise interactions between points are mediated by a bond potential of multi-well type while multi-point interactions are associated with the volume change mediated by a hydrostatic strain potential.The hydrostatic potential can either be a quadratic function,delivering a linear force–strain relation,or a multi-well type that can be associated with the material degradation and cavitation.We first show the well-posedness of the peridynamic formulation and that peridynamic evolutions exist in the Sobolev space H2.We show that the finite element approximations converge to the H2 solutions uniformly as measured in the mean square norm.For linear continuous fi nite elements,the convergence rate is shown to be Ct Δt+Csh2/ε2,where𝜖is the size of the horizon,his the mesh size,and Δt is the size of the time step.The constants Ct and Cs are independent of Δt and h and may depend on ε through the norm of the exact solution.We demonstrate the stability of the semi-discrete approximation.The stability of the fully discrete approximation is shown for the linearized peridynamic force.We present numerical simulations with the dynamic crack propagation that support the theoretical convergence rate.展开更多
We solve the local uniaxial tension of an infinite rod in the framework of non-ordinary state-based peridynamics.The singular solutions of stress and displacement are acquired.When the influencing range of the window ...We solve the local uniaxial tension of an infinite rod in the framework of non-ordinary state-based peridynamics.The singular solutions of stress and displacement are acquired.When the influencing range of the window function approaches zero,these two solutions will return to the solutions of the classical elasticity.The analysis shows that the singularities of the solutions stem from such a feature of the window function that must be represented by a rapidly decreasing function in physics.Contrary to the classical elasticity,the stress solution of peridynamics is smoother than the displacement solution.In addition,a criterion used to select the window function is proposed in this paper.展开更多
Sea ice exhibits complex mechanical properties,and no unified constitutive model currently exists.This study establishes an elastoplastic sea ice constitutive model based on non-ordinary state-based Peridynamics(PD)an...Sea ice exhibits complex mechanical properties,and no unified constitutive model currently exists.This study establishes an elastoplastic sea ice constitutive model based on non-ordinary state-based Peridynamics(PD)and the TsaiWu yield criterion,applying force state calculations to sea ice collisions while mitigating zero energy modes.A Fortran program implements the elastic-plastic constitutive equation of PD to simulate spherical ice-steel plate collisions.The program's accuracy in simulating sea ice collisions is validated through comparison with finite element results.Using the established model,this study simulates collisions between vertical structures and layer ice,analyzing the effects of impact velocity,vertical structure size,and critical elongation on sea ice load.The findings demonstrate positive correlations between collision force and impact velocity,vertical structure size,and critical elongation.At high velocities,impact significantly affects collision force,primarily following a quadratic function,while vertical structure effects exhibit a linear relationship.展开更多
Peridynamics(PD)is an effective method for simulating the spontaneous initiation and propagation of tensile cracks in materials.However,it faces great challenges in simulating compression-shear cracking of geomaterial...Peridynamics(PD)is an effective method for simulating the spontaneous initiation and propagation of tensile cracks in materials.However,it faces great challenges in simulating compression-shear cracking of geomaterials due to the lack of efficient contact-friction models.This paper introduces an original contact-friction model that leverages twin mesh and potential function principles within PD to model rock cracking under tensile and compressive stresses.The contact detection algorithm,based on space segmentation axis-aligned bounding box(AABB)tree data structure,is used to address the significant challenge of highly efficient contact detection in compression and shear problems.In this method,the twin mesh and potential function are utilized to quantify contact detection and contact degree,as well as friction behavior.This is in contrast to the distance and circular contact area model,which lacks physical significance in the classical PD method.As demonstrated by the tests on specimens containing cracks,the proposed model can capture 8 types of secondary fractures,reduce the contact detection error by about 29%e56%,and increase the contact retrieval efficiency by over 1600 times compared to the classic PD models.This significantly enhances the capability of PD to simulate the initiation,expansion,and coalescence of intricate compression-shear cracks.展开更多
As a typical brittle material,glass is widely used in construction,transportation,shipbuilding,aviation,aerospace and other industries.The unsafe factors of glass mainly come from its rupture.Thus,establishing a set o...As a typical brittle material,glass is widely used in construction,transportation,shipbuilding,aviation,aerospace and other industries.The unsafe factors of glass mainly come from its rupture.Thus,establishing a set of prediction models for the cracks growth of glass under dynamic load is necessary.This paper presents a contact damage model for glass based on the ordinary state-based peridynamic theory by introducing a contact force function.The Hertz contact(nonembedded contact)problem is simulated,and the elastic contact force is determined by adjusting the penalty factor.The proposed model verifies the feasibility of penalty-based method to simulate the contact problem of glass.The failure process of glass specimen under impact is simulated,where two loading methods,the drop ball test and the split Hopkinson pressure bar are considered.Numerical results agree well with the experimental observations,thereby verifying the effectiveness of the proposed model.展开更多
In this study,a numerical method was developed based on peridynamics to determine the ice loads for a ship navigating in level ice.Convergence analysis of threedimensional ice specimen with tensile and compression loa...In this study,a numerical method was developed based on peridynamics to determine the ice loads for a ship navigating in level ice.Convergence analysis of threedimensional ice specimen with tensile and compression loading are carried out first.The effects of ice thickness,sailing speed,and ice properties on the mean ice loads were also investigated.It is observed that the ice fragments resulting from the icebreaking process will interact with one another as well as with the water and ship hull.The ice fragments may rotate,collide,or slide along the ship hull,and these ice fragments will eventually drift away from the ship.The key characteristics of the icebreaking process can be obtained using the peridynamic model such as the dynamic generation of cracks in the ice sheet,propagation and accumulation of ice fragments,as well as collision,rotation,and sliding of the ice fragments along the ship hull.The simulation results obtained for the ice loads and icebreaking process were validated against those determined from the Lindqvist empirical formula and there is good agreement between the results.展开更多
In this study,we propose the first unified implementation strategy for peridynamics in commercial finite element method(FEM)software packages based on their application programming interface using the peridynamics-bas...In this study,we propose the first unified implementation strategy for peridynamics in commercial finite element method(FEM)software packages based on their application programming interface using the peridynamics-based finite element method(PeriFEM).Using ANSYS and ABAQUS as examples,we present the numerical results and implementation details of PeriFEM in commercial FEM software.PeriFEM is a reformulation of the traditional FEM for solving peridynamic equations numerically.It is considered that the non-local features of peridynamics yet possesses the same computational framework as the traditional FEM.Therefore,this implementation benefits from the consistent computational frameworks of both PeriFEM and the traditional FEM.An implicit algorithm is used for both ANSYS and ABAQUS;however,different convergence criteria are adopted owing to their unique features.In ANSYS,APDL enables users to conveniently obtain broken-bond information from UPFs;thus,the convergence criterion is chosen as no new broken bond.In ABAQUS,obtaining broken-bond information is not convenient for users;thus,the default convergence criterion is used in ABAQUS.The codes integrated into ANSYS and ABAQUS are both verified through benchmark examples,and the computational convergence and costs are compared.The results show that,for some specific examples,ABAQUS is more efficient,whereas the convergence criterion adopted in ANSYS is more robust.Finally,3D examples are presented to demonstrate the ability of the proposed approach to deal with complex engineering problems.展开更多
In the ice-covered oceanic region,the collision between sea ice and offshore structures will occur,causing the crushing failure of ice and the vibration of structures.The vibration can result in fatigue damage of stru...In the ice-covered oceanic region,the collision between sea ice and offshore structures will occur,causing the crushing failure of ice and the vibration of structures.The vibration can result in fatigue damage of structure and even endanger the crews’health.It is no doubt that this ice-structure interaction has been noted with great interest by the academic community for a long time and numerous studies have been done through theoretical analysis,experimental statistics and numerical simulation.In this paper,the bond-based Peridynamics method is applied to simulate the interaction between sea ice and wide vertical structures,where sea ice is modeled as elastic-plastic material,with a certain yield condition and failure criterion.Oscillation equation of single-degree-of-freedom is considered to investigate the vibration features of the structure during the interaction process.The damage of ice,ice forces and vibration responses of structure in the duration are obtained through numerical simulation.A parametric investigation is undertaken to identify the key parameters,such as ice thickness,the diameter of structure and relative velocity that trigger the ice crushing,ice forces and vibration responses of the structure.Results indicate that all three parameters have a positive correlation with the overall level of ice force and vibration displacement.Besides,a velocity coefficient is proposed to predict the vibration displacement based on its relation with ice speed.展开更多
Peridynamics(PD)is a widely used theory to simulate discontinuities,but its application in real-world structural problems is somewhat limited due to the relatively low-efficiency.The numerical substructure method(NSM)...Peridynamics(PD)is a widely used theory to simulate discontinuities,but its application in real-world structural problems is somewhat limited due to the relatively low-efficiency.The numerical substructure method(NSM)presented by the authors and co-workers provides an efficient approach for modeling structures with local nonlinearities,which is usually restricted in problems of continuum mechanics.In this paper,an approach is presented to couple the PD theory with the NSM for modeling structures with local discontinuities,taking advantage of the powerful capability of the PD for discontinuities simulation and high computational efficiency of the NSM.The structure is simulated using liner elastic finite element(FE)model while the local cracking regions are isolated and simulated using a PD substructure model.A force corrector calculated from the PD model is applied on the FE model to consider the effect of discontinuities.The PD is integrated in the substructure model using interface elements with embedded PD nodes.The equations of motions of both the NSM system and the PD substructure are solved using the central difference method.Three examples of two-dimensional(2D)concrete cantilever beams under the concentrated force are investigated to verify the proposed coupling approach.展开更多
In this paper,a peridynamics-based finite element method(PeriFEM)is proposed for the quasi-static fracture analysis,which is of the consistent computational framework with the classical finite element method(FEM).Firs...In this paper,a peridynamics-based finite element method(PeriFEM)is proposed for the quasi-static fracture analysis,which is of the consistent computational framework with the classical finite element method(FEM).First,the integral domain of peridynamics is reconstructed,and a new type of element called peridynamic element(PE)is defined.Although PEs are generated by the continuous elements(CEs)of classical FEM,they do not affect each other.Then,spatial discretization is performed based on PEs and CEs,and the linear equations about nodal displacement are established according to the principle of minimum potential energy.Besides,cracks are characterized as degradation of the mechanical properties of PEs.Finally,the validity of the proposed method is demonstrated through numerical examples.展开更多
A dual-support smoothed particle hydrodynamics(DS-SPH)that allows variable smoothing lengths while satisfying the conservations of linear momentum,angular momentum and energy is developed.The present DS-SPH is inspire...A dual-support smoothed particle hydrodynamics(DS-SPH)that allows variable smoothing lengths while satisfying the conservations of linear momentum,angular momentum and energy is developed.The present DS-SPH is inspired by the dual-support,a concept introduced from dual-horizon peridynamics from the authors and applied here to SPH so that the unbalanced interactions between the particles with different smoothing lengths can be correctly considered and computed.Conventionally,the SPH formulation employs either the influence domain or the support domain.The concept of dual-support identifies that the influence domain and the support domain involves the duality and should be simultaneously in the SPH formulation when variable smoothing lengths are used.The DS-SPH formulation can be implemented into conventional SPH codes with minimal changes and also without compromising the computational efficiency.A number of numerical examples involving weakly compressible.fluid are presented to demonstrate the capability of the method.展开更多
The ablation of ultra-high-temperature ceramics(UTHCs)is a complex physicochemical process including mechanical behavior,temperature effect,and chemical reactions.In order to realize the structural optimization and fu...The ablation of ultra-high-temperature ceramics(UTHCs)is a complex physicochemical process including mechanical behavior,temperature effect,and chemical reactions.In order to realize the structural optimization and functional design of ultra-high temperature ceramics,a coupled thermo-chemo-mechanical bond-based peridynamics(PD)model is proposed based on the ZrB_(2) ceramics oxidation kinetics model and coupled thermomechanical bond-based peridynamics.Compared with the traditional coupled thermo-mechanical model,the proposedmodel considers the influenceof chemical reactionprocessonthe ablation resistanceof ceramicmaterials.In order to verify the reliability of the proposed model,the thermo-mechanical coupling model,damage model and oxidation kinetic model are established respectively to investigate the applicability of the proposedmodel proposed in dealing with thermo-mechanical coupling,crack propagation,and chemical reaction,and the results show that the model is reliable.Finally,the coupled thermo-mechanical model and coupled thermo-chemo-mechanical model are used to simulate the crack propagation process of the plate under the thermal shock load,and the results show that the oxide layer plays a good role in preventing heat transfer and protecting the internal materials.Based on the PD fully coupled thermo-mechanical model,this paper innovatively introduces the oxidation kinetic model to analyze the influence of parameter changes caused by oxide layer growth and chemical growth strain on the thermal protection ability of ceramics.The proposed model provides an effective simulation technology for the structural design of UTHCs.展开更多
基金This work is supported by the National Natural Science Foundation of China under Grant Nos.1190219711972234 and is sponsored by Shanghai Sailing Program under Contract No.19YF1421700.
文摘This study demonstrates a homogenization approach via a modified state-based peridynamic(PD)method to predict the effective elastic properties of composite materials with periodic microstructure.The procedure of modeling the PD unit cell(UC)of continuous fiber-reinforced composite is presented.Periodic boundary conditions are derived and implemented through the Lagrange multiplier method.A matrix-dominated approach for modeling the interphase properties between dissimilar materials is proposed.The periodicity and continuity assumptions are employed to determine the stress and strain fields,as well as the effective elastic properties.The PD-UCs of square and hexagonal packs as well as the 0/90 laminate microstructure are modeled and compared with the analytical,numerical and experimental results from the literature.Good agreement of predicted effective properties can be observed.Unlike other PD homogenization approaches,the effective material properties can be directly and individually obtained from simple loading conditions.
基金supported by the National Key R&D Program of China(2020YFA0710500).
文摘The peridynamics(PD),as a promising nonlocal continuum mechanics theory,shines in solving discontinuous problems.Up to now,various numerical methods,such as the peridynamic mesh-free particlemethod(PD-MPM),peridynamic finite element method(PD-FEM),and peridynamic boundary element method(PD-BEM),have been proposed.PD-BEM,in particular,outperforms other methods by eliminating spurious boundary softening,efficiently handling infinite problems,and ensuring high computational accuracy.However,the existing PD-BEM is constructed exclusively for bond-based peridynamics(BBPD)with fixed Poisson’s ratio,limiting its applicability to crack propagation problems and scenarios involving infinite or semi-infinite problems.In this paper,we address these limitations by introducing the boundary element method(BEM)for ordinary state-based peridynamics(OSPD-BEM).Additionally,we present a crack propagationmodel embeddedwithin the framework ofOSPD-BEM to simulate crack propagations.To validate the effectiveness of OSPD-BEM,we conduct four numerical examples:deformation under uniaxial loading,crack initiation in a double-notched specimen,wedge-splitting test,and threepoint bending test.The results demonstrate the accuracy and efficiency of OSPD-BEM,highlighting its capability to successfully eliminate spurious boundary softening phenomena under varying Poisson’s ratios.Moreover,OSPDBEMsignificantly reduces computational time and exhibits greater consistencywith experimental results compared to PD-MPM.
基金support from the National Natural Science Foundation of China(Grant Nos.42277161 and 42230709).
文摘In rock engineering,natural cracks in rock masses subjected to external loads tend to initiate and propagate,leading to potential safety hazards.To investigate the effect of cracking behavior on the mechanical properties of rocks,the cracking processes of pre-cracked rocks have been extensively studied using numerical modeling methods.The peridynamics(PD)exhibits advantages over other numerical methods due to the absence of the requirements for remeshing and external crack growth criterion.However,for modeling pre-cracked rock cracking processes under impact,current PD implementations lack generally applicable rock constitutive models and impact contact models,which leads to difficulties in determining rock material parameters and efficiently calculating impact loads.This paper proposes a non-ordinary state-based peridynamics(NOSBPD)modeling method integrating the Drucker-Prager(DP)plasticity model and an efficient contact model to address the above problems.In the proposed method,the Drucker-Prager plasticity model is integrated into the NOSBPD,thereby equipping NOSBPD with the capability to accurately characterize the nonlinear stress-strain relationship inherent in rocks.An efficient contact model between particles and meshes is designed to calculate the impact loads,which is essentially a coupling method of PD with the finite element method(FEM).The effectiveness of the proposed NOSBPD modeling method is verified by comparison with other numerical methods and experiments.Experimental results indicate that the proposed method can effectively and accurately predict the 3D cracking processes of pre-cracked cracks under impact loading,and the maximum principal stress is the key driver behind wing crack formation in pre-cracked rocks.
基金supported by the Fundamental Research Funds for the Central Universities(Grant B200202231)the National Natural Science Foundation of China(Grants 11932006,11672101,U1934206,and 12002118)+1 种基金the National Key Research&Development Plan of China(Grants 2018 YFC0406703 and 2017YFC1502603)the China Postdoctoral Science Foundation(Grant 2019M651667).
文摘The non-ordinary state-based peridynamic(NOSB PD)model has the capability of incorporating existing constitutive relationships in the classical continuum mechanics.In the present work,we first develop an NOSB PD model corresponding to the Johnson–Holmquist II(JH-2)constitutive damage model,which can describe the severe damage of concrete under intense impact compression.Besides,the numerical oscillation problem of the NOSB PD caused by zero-energy mode is analyzed and hence a bond-associated non-ordinary state-based peridynamic(BA-NOSB PD)model is adopted to remove the oscillation.Then,the elastic deformation of a three-dimensional bar is analyzed to verify the capability of BA-NOSB PD in eliminating the numerical oscillation.Furthermore,concrete spalling caused by the interaction of incident compression wave and reflected tension wave is simulated.The dynamic tensile fracture process of concrete multiple spalling is accurately reproduced for several examples according to the spalling number and spalling thickness analysis,illustrating the approach can well simulate and analyze the concrete spalling discontinuities.
基金The authors acknowledge the support from the National Key R&D Program of China(No.2018YFC0406703)National Natural Science Foundation of China(No.51679077)the Fundamental Research Funds for the Central Universities in China(No.2015B18314,2017B13014).
文摘Strain hardening and strain rate play an important role in dynamic deformation and failure problems such as high-velocity impact cases.In this paper,a non-ordinary state-based peridynamic model for failure and damage of concrete materials subjected to impacting condition is proposed,taking the advantages of both damage model and nonlocal peridynamic method.The Holmquist-Johnson-Cook(HJC)model describing the mechanical character and damage of concrete materials under large strain,high strain rate and high hydrostatic pressure was reformulated in the framework of non-ordinary statebased peridynamic theory,and the corresponding numerical approach was developed.The proposed model and numerical approach were validated through simulating typical impacting examples and comparing the results with available experimental observations and results in literature.
文摘In this study,a new state-based peridynamic formulation is developed for functionally graded Euler-Bernoulli beams.The equation of motion is developed by using Lagrange’s equation and Taylor series.Both axial and transverse displacements are taken into account as degrees of freedom.Four different boundary conditions are considered including pinned support-roller support,pinned support-pinned support,clamped-clamped and clamped-free.Peridynamic results are compared against finite element analysis results for transverse and axial deformations and a very good agreement is observed for all different types of boundary conditions.
文摘In this study,a recently developed peridynamic lattice model called the“Statebased Peridynamic Lattice Model”(SPLM)is improved and demonstrated.In the SPLM,rather than as a continuum,solids are simulated using a close-packed lattice of peridynamically interacting particles.The new SPLM approach advances the SPLM model by improving the damage and plasticity models.Elasticity,plasticity and damage are coupled in this approach.A robust method for damage initiation is developed.A new damage model called the“two-spring damage model”allows damage to localize to a single lattice particle,thus allowing highly localized damage(cracks)to emerge in a realistic manner.A plasticity model that includes hardening,softening,and damage due to plasticity is proposed and demonstrated.Peridynamic boundary effects are modeled efficiently and reasonably.The improved SPLM method is then employed to simulate three common concrete laboratory tests:uniaxial tension,uniaxial compression,and the Brazilian split cylinder test.The SPLM results are then compared with results from the earlier SPLM model,with simplified classical predictions,and with laboratory results.By solving the same benchmark problems using various lattice rotations and lattice spacings,the approach is demonstrated to be sufficiently objective to be a useful engineering tool to predict the essentially random behavior of concrete laboratory specimens.The improved SPLM demonstrates significant improvements over the previously published version and is found to simulate concrete structures accurately and efficiently using far less computational effort than comparable computational simulation methods.
文摘We establish the a priori convergence rate for finite element approximations of a class of nonlocal nonlinear fracture models.We consider state-based peridynamic models where the force at a material point is due to both the strain between two points and the change in volume inside the domain of the nonlocal interaction.The pairwise interactions between points are mediated by a bond potential of multi-well type while multi-point interactions are associated with the volume change mediated by a hydrostatic strain potential.The hydrostatic potential can either be a quadratic function,delivering a linear force–strain relation,or a multi-well type that can be associated with the material degradation and cavitation.We first show the well-posedness of the peridynamic formulation and that peridynamic evolutions exist in the Sobolev space H2.We show that the finite element approximations converge to the H2 solutions uniformly as measured in the mean square norm.For linear continuous fi nite elements,the convergence rate is shown to be Ct Δt+Csh2/ε2,where𝜖is the size of the horizon,his the mesh size,and Δt is the size of the time step.The constants Ct and Cs are independent of Δt and h and may depend on ε through the norm of the exact solution.We demonstrate the stability of the semi-discrete approximation.The stability of the fully discrete approximation is shown for the linearized peridynamic force.We present numerical simulations with the dynamic crack propagation that support the theoretical convergence rate.
基金the National Natural Science Foundation of China (11672129)the Research Fund of State Key Laboratory of Mechanics and Control of Mechanical Structures (Nanjing University of Aeronautics and Astronautics, MCMS-I-0218G01)
文摘We solve the local uniaxial tension of an infinite rod in the framework of non-ordinary state-based peridynamics.The singular solutions of stress and displacement are acquired.When the influencing range of the window function approaches zero,these two solutions will return to the solutions of the classical elasticity.The analysis shows that the singularities of the solutions stem from such a feature of the window function that must be represented by a rapidly decreasing function in physics.Contrary to the classical elasticity,the stress solution of peridynamics is smoother than the displacement solution.In addition,a criterion used to select the window function is proposed in this paper.
基金financially supported by the Aeronautical Science Foundation of China(Grant No.023M031077001)Shandong Provincial Natural Science Foundation(Grant No.ZR2022QE092)the Open Fund Project of Key Laboratory of Ocean Observation Technology,MNR(Grant No.2022klootA03)。
文摘Sea ice exhibits complex mechanical properties,and no unified constitutive model currently exists.This study establishes an elastoplastic sea ice constitutive model based on non-ordinary state-based Peridynamics(PD)and the TsaiWu yield criterion,applying force state calculations to sea ice collisions while mitigating zero energy modes.A Fortran program implements the elastic-plastic constitutive equation of PD to simulate spherical ice-steel plate collisions.The program's accuracy in simulating sea ice collisions is validated through comparison with finite element results.Using the established model,this study simulates collisions between vertical structures and layer ice,analyzing the effects of impact velocity,vertical structure size,and critical elongation on sea ice load.The findings demonstrate positive correlations between collision force and impact velocity,vertical structure size,and critical elongation.At high velocities,impact significantly affects collision force,primarily following a quadratic function,while vertical structure effects exhibit a linear relationship.
基金supported by the National Natural Science Foundation of China(Grant No.52278333)the China Scholarship Council(CSC)and the Science and Technology Department of Liaoning Province(Grant No.2024JH2/102500069).
文摘Peridynamics(PD)is an effective method for simulating the spontaneous initiation and propagation of tensile cracks in materials.However,it faces great challenges in simulating compression-shear cracking of geomaterials due to the lack of efficient contact-friction models.This paper introduces an original contact-friction model that leverages twin mesh and potential function principles within PD to model rock cracking under tensile and compressive stresses.The contact detection algorithm,based on space segmentation axis-aligned bounding box(AABB)tree data structure,is used to address the significant challenge of highly efficient contact detection in compression and shear problems.In this method,the twin mesh and potential function are utilized to quantify contact detection and contact degree,as well as friction behavior.This is in contrast to the distance and circular contact area model,which lacks physical significance in the classical PD method.As demonstrated by the tests on specimens containing cracks,the proposed model can capture 8 types of secondary fractures,reduce the contact detection error by about 29%e56%,and increase the contact retrieval efficiency by over 1600 times compared to the classic PD models.This significantly enhances the capability of PD to simulate the initiation,expansion,and coalescence of intricate compression-shear cracks.
基金This study was funded by National Natural Science Foundation of China(Nos.11932006,U1934206)Recipient:Qing Zhang.And National Natural Science Foundation of China(No.12002118)+1 种基金Recipient:Xin Gu.And Natural Science Foundation of the Jiangsu Higher Education Institutions of China(No.20KJB580015)Recipient:Runpu Li。
文摘As a typical brittle material,glass is widely used in construction,transportation,shipbuilding,aviation,aerospace and other industries.The unsafe factors of glass mainly come from its rupture.Thus,establishing a set of prediction models for the cracks growth of glass under dynamic load is necessary.This paper presents a contact damage model for glass based on the ordinary state-based peridynamic theory by introducing a contact force function.The Hertz contact(nonembedded contact)problem is simulated,and the elastic contact force is determined by adjusting the penalty factor.The proposed model verifies the feasibility of penalty-based method to simulate the contact problem of glass.The failure process of glass specimen under impact is simulated,where two loading methods,the drop ball test and the split Hopkinson pressure bar are considered.Numerical results agree well with the experimental observations,thereby verifying the effectiveness of the proposed model.
基金This project has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement(Grant No.2017YFE0111400)the National Key R&D Program Strategic International Science and Technology Innovation Cooperation Key Specialities(Grant No.2016YFE0202700)+3 种基金the National Natural Science Foundation of China(Grant Nos.51579054 and 51639004)the Ministry of Industry and Information Technology’s High-tech Ship Research Project(Grant No.2017-614)Mr.Renwei Liu is supported by a two-year visiting student fellowship in University of California,Berkeley from Chinese Scholar Council(Grant No.201706680104)this support is gratefully acknowledged.The authors also graciously acknowledge Professor Shaofan Li of University of California,Berkeley and Fei Han of Dalian University of Technology for their guidance and fruitful discussion regarding this work.
文摘In this study,a numerical method was developed based on peridynamics to determine the ice loads for a ship navigating in level ice.Convergence analysis of threedimensional ice specimen with tensile and compression loading are carried out first.The effects of ice thickness,sailing speed,and ice properties on the mean ice loads were also investigated.It is observed that the ice fragments resulting from the icebreaking process will interact with one another as well as with the water and ship hull.The ice fragments may rotate,collide,or slide along the ship hull,and these ice fragments will eventually drift away from the ship.The key characteristics of the icebreaking process can be obtained using the peridynamic model such as the dynamic generation of cracks in the ice sheet,propagation and accumulation of ice fragments,as well as collision,rotation,and sliding of the ice fragments along the ship hull.The simulation results obtained for the ice loads and icebreaking process were validated against those determined from the Lindqvist empirical formula and there is good agreement between the results.
基金the financial support received from the National Natural Science Foundation of China(12272082,11872016)the National Key Laboratory of Shock Wave and Detonation Physics(JCKYS2021212003).
文摘In this study,we propose the first unified implementation strategy for peridynamics in commercial finite element method(FEM)software packages based on their application programming interface using the peridynamics-based finite element method(PeriFEM).Using ANSYS and ABAQUS as examples,we present the numerical results and implementation details of PeriFEM in commercial FEM software.PeriFEM is a reformulation of the traditional FEM for solving peridynamic equations numerically.It is considered that the non-local features of peridynamics yet possesses the same computational framework as the traditional FEM.Therefore,this implementation benefits from the consistent computational frameworks of both PeriFEM and the traditional FEM.An implicit algorithm is used for both ANSYS and ABAQUS;however,different convergence criteria are adopted owing to their unique features.In ANSYS,APDL enables users to conveniently obtain broken-bond information from UPFs;thus,the convergence criterion is chosen as no new broken bond.In ABAQUS,obtaining broken-bond information is not convenient for users;thus,the default convergence criterion is used in ABAQUS.The codes integrated into ANSYS and ABAQUS are both verified through benchmark examples,and the computational convergence and costs are compared.The results show that,for some specific examples,ABAQUS is more efficient,whereas the convergence criterion adopted in ANSYS is more robust.Finally,3D examples are presented to demonstrate the ability of the proposed approach to deal with complex engineering problems.
基金This work is supported financially by the National Key R&D Program of China[2018YFC1406000,2016YFE0202700]Supported by the National Natural Science Foundation of China(NSFC)[Grant Nos.51809061,51639004]+1 种基金Supported by the Natural Science Foundation of Heilongjiang Province of China[LC2018021]Supported by the Fundamental Research Funds for the Central Universities[HEUCFM180111].
文摘In the ice-covered oceanic region,the collision between sea ice and offshore structures will occur,causing the crushing failure of ice and the vibration of structures.The vibration can result in fatigue damage of structure and even endanger the crews’health.It is no doubt that this ice-structure interaction has been noted with great interest by the academic community for a long time and numerous studies have been done through theoretical analysis,experimental statistics and numerical simulation.In this paper,the bond-based Peridynamics method is applied to simulate the interaction between sea ice and wide vertical structures,where sea ice is modeled as elastic-plastic material,with a certain yield condition and failure criterion.Oscillation equation of single-degree-of-freedom is considered to investigate the vibration features of the structure during the interaction process.The damage of ice,ice forces and vibration responses of structure in the duration are obtained through numerical simulation.A parametric investigation is undertaken to identify the key parameters,such as ice thickness,the diameter of structure and relative velocity that trigger the ice crushing,ice forces and vibration responses of the structure.Results indicate that all three parameters have a positive correlation with the overall level of ice force and vibration displacement.Besides,a velocity coefficient is proposed to predict the vibration displacement based on its relation with ice speed.
基金Financial support by the National Key Research and Development program of China under Grant No.2016YFC0701106the National Natural Science Foundation of China under grants No.51578473the program of China Scholarship Council(CSC,No.201606060083)are gratefully acknowledged.
文摘Peridynamics(PD)is a widely used theory to simulate discontinuities,but its application in real-world structural problems is somewhat limited due to the relatively low-efficiency.The numerical substructure method(NSM)presented by the authors and co-workers provides an efficient approach for modeling structures with local nonlinearities,which is usually restricted in problems of continuum mechanics.In this paper,an approach is presented to couple the PD theory with the NSM for modeling structures with local discontinuities,taking advantage of the powerful capability of the PD for discontinuities simulation and high computational efficiency of the NSM.The structure is simulated using liner elastic finite element(FE)model while the local cracking regions are isolated and simulated using a PD substructure model.A force corrector calculated from the PD model is applied on the FE model to consider the effect of discontinuities.The PD is integrated in the substructure model using interface elements with embedded PD nodes.The equations of motions of both the NSM system and the PD substructure are solved using the central difference method.Three examples of two-dimensional(2D)concrete cantilever beams under the concentrated force are investigated to verify the proposed coupling approach.
基金The authors gratefully acknowledge the financial support from the National Natural Science Foundation of China(11872016)Fundamental Research Funds for the Central Universities(DUT20RC(5)005,DUT20LAB203)the Key Research and Development Project of Liaoning Province(2020JH2/10500003).
文摘In this paper,a peridynamics-based finite element method(PeriFEM)is proposed for the quasi-static fracture analysis,which is of the consistent computational framework with the classical finite element method(FEM).First,the integral domain of peridynamics is reconstructed,and a new type of element called peridynamic element(PE)is defined.Although PEs are generated by the continuous elements(CEs)of classical FEM,they do not affect each other.Then,spatial discretization is performed based on PEs and CEs,and the linear equations about nodal displacement are established according to the principle of minimum potential energy.Besides,cracks are characterized as degradation of the mechanical properties of PEs.Finally,the validity of the proposed method is demonstrated through numerical examples.
基金The authors acknowledge the supports from the ERC-CoG(Computational Modeling and Design of Lithium-ion Batteries(COMBAT)),RISE-BESTOFRAC and National Science Foundation of China(51474157).
文摘A dual-support smoothed particle hydrodynamics(DS-SPH)that allows variable smoothing lengths while satisfying the conservations of linear momentum,angular momentum and energy is developed.The present DS-SPH is inspired by the dual-support,a concept introduced from dual-horizon peridynamics from the authors and applied here to SPH so that the unbalanced interactions between the particles with different smoothing lengths can be correctly considered and computed.Conventionally,the SPH formulation employs either the influence domain or the support domain.The concept of dual-support identifies that the influence domain and the support domain involves the duality and should be simultaneously in the SPH formulation when variable smoothing lengths are used.The DS-SPH formulation can be implemented into conventional SPH codes with minimal changes and also without compromising the computational efficiency.A number of numerical examples involving weakly compressible.fluid are presented to demonstrate the capability of the method.
基金support from the National Natural Science Foun-dation of China(11972267).
文摘The ablation of ultra-high-temperature ceramics(UTHCs)is a complex physicochemical process including mechanical behavior,temperature effect,and chemical reactions.In order to realize the structural optimization and functional design of ultra-high temperature ceramics,a coupled thermo-chemo-mechanical bond-based peridynamics(PD)model is proposed based on the ZrB_(2) ceramics oxidation kinetics model and coupled thermomechanical bond-based peridynamics.Compared with the traditional coupled thermo-mechanical model,the proposedmodel considers the influenceof chemical reactionprocessonthe ablation resistanceof ceramicmaterials.In order to verify the reliability of the proposed model,the thermo-mechanical coupling model,damage model and oxidation kinetic model are established respectively to investigate the applicability of the proposedmodel proposed in dealing with thermo-mechanical coupling,crack propagation,and chemical reaction,and the results show that the model is reliable.Finally,the coupled thermo-mechanical model and coupled thermo-chemo-mechanical model are used to simulate the crack propagation process of the plate under the thermal shock load,and the results show that the oxide layer plays a good role in preventing heat transfer and protecting the internal materials.Based on the PD fully coupled thermo-mechanical model,this paper innovatively introduces the oxidation kinetic model to analyze the influence of parameter changes caused by oxide layer growth and chemical growth strain on the thermal protection ability of ceramics.The proposed model provides an effective simulation technology for the structural design of UTHCs.
