An Elastoplastic Fracture Model Based on Bond-Based Peridynamics

Fracture in ductile materials often occurs in conjunction with plastic deformation.However,in the bond-based peridynamic(BB-PD)theory,the classic mechanical stress is not defined inherently.This makes it difficult to describe plasticity directly using the classical plastic theory.To address the above issue,a unified bond-based peridynamics model was proposed as an effective tool to solve elastoplastic fracture problems.Compared to the existing models,the proposed model directly describes the elastoplastic theory at the bond level without the need for additional calculation means.The results obtained in the context of this model are shown to be consistent with FEM results in regard to force-displacement curves,displacement fields,stress fields,and plastic deformation regions.The model exhibits good capability of capturing crack propagation in ductile material failure problems.

The Boundary Element Method for Ordinary State-Based Peridynamics

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.

ABAQUS and ANSYS Implementations of the Peridynamics-Based Finite Element Method(PeriFEM)for Brittle Fractures

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.

Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures

Molecular dynamics(MD)has served as a powerful tool for designing materials with reduced reliance on laboratory testing.However,the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach.In this work,we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets.Firstly,we develop a novel coarse-graining method,to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets.Inspired by the weighted essentially non-oscillatory(WENO)scheme,the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil,then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities.Then,based on the coarse-grained MD data,a two-phase optimizationbased learning approach is proposed to infer the optimal peridynamics model with damage criterion.In the first phase,we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties.Then,in the second phase,the material damage criterion is learnt as a smoothed step function from the data with fractures.As a result,a peridynamics surrogate is obtained.As a continuum model,our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training,and hence allows for substantial reductions in computational cost compared with MD.We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene.Our tests show that the proposed data-driven model is robust and generalizable,in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.

Study on Crack Propagation Parameters of Tunnel Lining Structure Based on Peridynamics

The numerical simulation results utilizing the Peridynamics(PD)method reveal that the initial crack and crack propagation of the tunnel concrete lining structure agree with the experimental data compared to the Japanese prototype lining test.The load structure model takes into account the cracking process and distribution of the lining segment under the influence of local bias pressure and lining thickness.In addition,the influence of preset cracks and lining section formon the crack propagation of the concrete lining model is studied.This study evaluates the stability and sustainability of tunnel structure by the Peridynamics method,which provides a reference for the analysis of the causes of lining cracks,and also lays a foundation for the prevention,reinforcement and repair of tunnel lining cracks.

The Coupled Thermo-Chemo-Mechanical Peridynamics for ZrB_(2) Ceramics Ablation Behavior

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.

Numerical Simulations of the Ice Load of a Ship Navigating in Level Ice Using Peridynamics

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.

A Dual-Support Smoothed Particle Hydrodynamics for Weakly Compressible Fluid Inspired By the Dual-Horizon Peridynamics

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.

Numerical Simulation of Dynamic Interaction Between Ice and Wide Vertical Structure Based on Peridynamics

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.

A Possible Reason About Origin of Singularity and Anomalous Dispersion in Peridynamics

In the benchmark problems of peridynamics,there are some eccentric results,for example,singularity of uniaxial tension and anomalous dispersion of wave.The reasons to give rise to these results are investigated.We calculated local tension and wave of an infinite rod after adding a divergence of local stress in the peridynamic motion equation.The acquired results verify that the singularity in the peridynamic solution of local tension problem and anomalous dispersion of peridynamic wave are all eliminated.Therefore,the anomalous features of some peridynamic solutions likely stem from the lack of local stress characterizing contact interactions.

The Multi-Horizon Peridynamics

We present a novel refinement approach in peridynamics(PD).The proposed approach takes advantage of the PD flexibility in choosing the shape of the horizon by introducing multiple domains(with no intersections)to the nodes of the refinement zone.We will show that no ghost forces are needed when changing the horizon sizes in both subdomains.The approach is applied to both bond-based and state-based peridynamics and verified for a simple wave propagation refinement problem illustrating the efficiency of the method.

Coupled Digital Image Correlation and Peridynamics for Full-Field Deformation Measurement and Local Damage Prediction

Digital image correlation(DIC)measurement technique and peridynamics(PD)method have been applied in specific fields extensively owing to their respective advantages in obtaining full-field deformation and local failure of loaded materials and structures.This study provides a simple way to couple DIC measurements with PD simulations,which can circumvent the difficulties of DIC in dealing with discontinuous deformations.Taking the failure analysis of a compact tension specimen of aluminum alloy and a static three-point bending concrete beam as examples,the DIC experimental system firstly measures the full-field displacements,and then the PD simulation is applied on potential damage regions determined according to the correlation coefficients,to track the micro-crack evolution and macro-crack propagation.As results,the coupled DIC and PD approach can effectively measure the full-field displacement and the localized damage accumulation and crack propagation.

Coupling of Peridynamics and Numerical Substructure Method for Modeling Structures with Local Discontinuities

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.

Modelling of Contact Damage in Brittle Materials Based on Peridynamics

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.

A coupled thermo-mechanical peridynamic model for fracture behavior of granite subjected to heating and water-cooling processes

Thermal damage and thermal fracture of rocks are two important indicators in geothermal mining projects.This paper investigates the effects of heating and water-cooling on granite specimens at various temperatures.The laboratory uniaxial compression experiments were also conducted.Then,a coupled thermo-mechanical ordinary state-based peridynamic(OSB-PD)model and corresponding numerical scheme were developed to simulate the damage of rocks after the heating and cooling processes,and the change of crack evolution process was predicted.The results demonstrate that elevated heating temperatures exacerbate the thermal damage to the specimens,resulting in a decrease in peak strength and an increase in ductility of granite.The escalating occurrence of thermal-induced cracks significantly affects the crack evolution process during the loading phase.The numerical results accurately reproduce the damage and fracture characteristics of the granite under different final heating temperatures(FHTs),which are consistent with the test results in terms of strength,crack evolution process,and failure mode.

Sub-Homogeneous Peridynamic Model for Fracture and Failure Analysis of Roadway Surrounding Rock

The surrounding rock of roadways exhibits intricate characteristics of discontinuity and heterogeneity.To address these complexities,this study employs non-local Peridynamics(PD)theory and reconstructs the kernel function to represent accurately the spatial decline of long-range force.Additionally,modifications to the traditional bondbased PD model are made.By considering the micro-structure of coal-rock materials within a uniform discrete model,heterogeneity characterized by bond random pre-breaking is introduced.This approach facilitates the proposal of a novel model capable of handling the random distribution characteristics of material heterogeneity,rendering the PD model suitable for analyzing the deformation and failure of heterogeneous layered coal-rock mass structures.The established numerical model and simulation method,termed the sub-homogeneous PD model,not only incorporates the support effect but also captures accurately the random heterogeneous micro-structure of roadway surrounding rock.The simulation results obtained using this model show good agreement with field measurements from the Fucun coal mine,effectively validating the model’s capability in accurately reproducing the deformation and failure mode of surrounding rock under bolt-supported(anchor cable).The proposed subhomogeneous PD model presents a valuable and effective simulation tool for studying the deformation and failure of roadway surrounding rock in coal mines,offering new insights and potential advancements.

A fatigue damage-cumulative model in peridynamics

While the present structural integrity evaluation method is based on the philosophy of assumed similitude, Fatigue and Damage Tolerance(F&DT) evaluations for next generation of air-vehicles require high-fidelity physical models within cyberspace. To serve the needs of F&DT evaluation in digital twin paradigm, a fatigue damage-cumulative model within peridynamic framework is proposed in this paper. Based on the concept of fatigue element block and damage accumulation law in form of Coffin-Manson relationship, the proposed model applies to both fatigue crack initiation and fatigue crack growth;fatigue crack growth rates under constant-amplitude and simple variable-amplitude block loading cases can be well predicted for three common structural materials without inputs of Paris law parameters. Additionally, the proposed model can also be easily extended to a probabilistic version;for verification, multiple-site-damage problems are simulated and the statistic nature of fatigue process in experiments can be well captured. In the end, main features of the proposed model are summarized, and distinctions from the other models are discussed. There may be a potential for the peridynamic damage-cumulative model proposed in this work to numerically predict fatigue problems in digital twin paradigm for future generations of aerospace vehicles.

A Coupled Thermomechanical Crack Propagation Behavior of Brittle Materials by Peridynamic Differential Operator

This study proposes a comprehensive,coupled thermomechanical model that replaces local spatial derivatives in classical differential thermomechanical equations with nonlocal integral forms derived from the peridynamic differential operator(PDDO),eliminating the need for calibration procedures.The model employs a multi-rate explicit time integration scheme to handle varying time scales in multi-physics systems.Through simulations conducted on granite and ceramic materials,this model demonstrates its effectiveness.It successfully simulates thermal damage behavior in granite arising from incompatible mineral expansion and accurately calculates thermal crack propagation in ceramic slabs during quenching.To account for material heterogeneity,the model utilizes the Shuffle algorithm andWeibull distribution,yielding results that align with numerical simulations and experimental observations.This coupled thermomechanical model shows great promise for analyzing intricate thermomechanical phenomena in brittle materials.

Finite Difference-Peridynamic Differential Operator for Solving Transient Heat Conduction Problems

Transient heat conduction problems widely exist in engineering.In previous work on the peridynamic differential operator(PDDO)method for solving such problems,both time and spatial derivatives were discretized using the PDDO method,resulting in increased complexity and programming difficulty.In this work,the forward difference formula,the backward difference formula,and the centered difference formula are used to discretize the time derivative,while the PDDO method is used to discretize the spatial derivative.Three new schemes for solving transient heat conduction equations have been developed,namely,the forward-in-time and PDDO in space(FT-PDDO)scheme,the backward-in-time and PDDO in space(BT-PDDO)scheme,and the central-in-time and PDDO in space(CT-PDDO)scheme.The stability and convergence of these schemes are analyzed using the Fourier method and Taylor’s theorem.Results show that the FT-PDDO scheme is conditionally stable,whereas the BT-PDDO and CT-PDDO schemes are unconditionally stable.The stability conditions for the FT-PDDO scheme are less stringent than those of the explicit finite element method and explicit finite difference method.The convergence rate in space for these three methods is two.These constructed schemes are applied to solve one-dimensional and two-dimensional transient heat conduction problems.The accuracy and validity of the schemes are verified by comparison with analytical solutions.

Euler’s First-Order Explicit Method–Peridynamic Differential Operator for Solving Population Balance Equations of the Crystallization Process

Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.