• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.3
• /
• pp.309-316
• /
• 2015

A state-based peridynamic model is able to describe a general constitutive model from the standard continuum theory. The response of a material at a point is dependent on the deformation of all bonds connected to the point within the nonlocal horizon region. Therefore, the state-based peridynamic model permits both the volume and shear changes of the material which is promising to reproduce the complicated dynamic brittle fracture phenomena, such as crack branching, secondary cracks, cascade cracks, crack coalescence, etc. In this paper, the two-dimensional state-based peridynamic model for a linear elastic plane stress solid is employed. The damage model incorporates the energy release rate and the peridynamic energy potential. For brittle glass materials, the impact of the crack-parallel compressive stress waves on the crack branching pattern is investigated. The peridynamic solution for this problem captures the main features, observed experimentally, of dynamic crack propagation and branching. Cascade cracks under strong tensile loading and secondary cracks are also well reproduced with the state-based peridynamic simulations.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.2
• /
• pp.87-94
• /
• 2014

In solid mechanics, the peridynamics theory has provided a suitable framework for material failure and damage propagation simulation. Peridynamics is computationally expensive since it is required to solve enormous nonlocal interactions based upon integro-differential equations. Thus, multiscale coupling methods with other local models are of interest for efficient and accurate implementations of peridynamics. In this study, peridynamic models are restricted to regions where discontinuities or stress concentrations are present. In the domains characterized by smooth displacements, classical local models can be employed. We introduce a recently developed blending scheme to concurrently couple bond-based peridynamic models and the Navier equation of classical elasticity. We demonstrate numerically that the proposed blended model is suitable for point loads and static fracture, suggesting an alternative framework for cases where peridynamic models are too expensive, while classical local models are not accurate enough.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.4
• /
• pp.257-263
• /
• 2019

Multilayered structures include lamination by relatively thick plies and thin interlayers. For efficient peridynamic analysis of dynamic fracturing multilayered structures, the interlayer is modeled using ghost peridynamic particles while the ply is formulated via real peridynamics. With the nonlocal ghost interlayer, one may keep the discretization resolution low for the ply. In this study, the characteristics of dynamic wave propagation through the nonlocal ghost interlayer in peridynamic analysis are investigated. It is observed that the interlayer not only binds adjacent plies, but also significantly influences energy transfer between plies, and thereby their deformation and motion. In addition, near a surface or boundary, peridynamic particles do not have a full nonlocal neighborhoods. This causes the effective material properties near the surface to be different from those in the bulk. Surface correction based on neighborhood volumes is employed. The impact of surface correction on wave propagation in multilayered structures is investigated.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.2
• /
• pp.105-114
• /
• 2018

In this paper, the computer simulations are carried out by using the peridynamic theory model with various conditions including quasi-static loads, dynamic loads and crack propagation, branching crack pattern and isotropic materials, orthotropic materials. Three examples, a plate with a hole under quasi-static loading, a plate with a pre-existing crack under dynamic loading and a lamina with a pre-existing crack under quasi-static loading are analyzed by computational simulations. In order to simulate the quasi-static load, an adaptive dynamic relaxation technique is used. In the orthotropic material analysis, a homogenization method is used considering the strain energy density ratio between the classical continuum mechanics and the peridynamic. As a result, crack propagation and branching cracks are observed successfully and the direction and initiation of the crack are also captured within the peridynamic modeling. In case of applying peridynamic used homogenization method to a relatively complicated orthotropic material, it is also verified by comparing with experimental results.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.3
• /
• pp.151-157
• /
• 2021

This paper is devoted to a convergence study of the nonlocal integral operator in peridynamics. The implicit formulation can be an efficient approach to obtain the static/quasi-static solution of crack propagation problems. Implicit methods require constly large-matrix operations. Therefore, convergence is important for improving computational efficiency. When the radial influence function is utilized in the nonlocal integral equation, the fractional Laplacian integral equation is obtained. It has been mathematically proved that the condition number of the system matrix is affected by the order of the radial influence function and nonlocal horizon size. We formulate the static crack problem with peridynamics and utilize Newton-Raphson methods with a preconditioned conjugate gradient scheme to solve this nonlinear stationary system. The convergence behavior and the computational time for solving the implicit algebraic system have been studied with respect to the order of the radial influence function and nonlocal horizon size.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.3
• /
• pp.301-307
• /
• 2015

A bond-based peridynamic model has been shown to be capable of analyzing many of dynamic brittle fracture phenomena. However, there have been issued limitations on handling constitutive models of various materials. Especially, it assumes bonds act independently of each other, so that Poisson's ratio for 3D model is fixed as 1/4 as well as taking only account the bond stretching results in a volume change not a shear change. In this paper a state-based peridynamic model of dynamic brittle fracture is presented. The state-based peridynamic model is a generalized peridynamic model that is able to directly use a constitutive model from the standard theory. It permits the response of a material at a point to depend collectively on the deformation of all bonds connected to the point. Thus, the volume and shear changes of the material can be reproduced by the state-based peridynamic theory. For a linearly elastic solid, a plane stress model is introduced and the damage model suitable for the state-based peridynamic model is discussed. Through a convergence study under decreasing the peridynamic nonlocal region($\delta$-convergence), the dynamic fracture model is verified. It is also shown that the state-based peridynamic model is reliable for modeling dynamic crack propagatoin.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.5
• /
• pp.389-396
• /
• 2017

Peridynamics has been widely used in the dynamic fracture analysis of brittle materials. Recently, various crack patterns(compact region, floret, Hertz-type crack, etc.) of multilayered glass structures in experiments(Bless et al. 2010) were implemented with a bond-based peridynamic simulation(Bobaru et al.. 2012). The actual glass layers are bound with thin elastic interlayer material while the interlayer is missing from the peridynamic model used in the previous numerical study. In this study, the peridynamic interlayer modeling for the multilayered structures is proposed. It requires enormous computational time and memory to explicitly model very thin interlayer materials. Instead of explicit modeling, fictitious peridynamic particles are introduced for modeling interlayer materials. The computational efficiency and accuracy of the proposed peridynamic interlayer model are verified through numerical tests. Furthermore, preventing penetration scheme based on short-range interaction force is employed for the multilayered structure under compression and verified through parametric tests.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.6
• /
• pp.373-380
• /
• 2018

We present the peridynamic dynamic fracture analysis to solve impact fracturing of multilayered glass impacted by a high-velocity object. In the most practical multilayered glass structures, main layers are glued by thin elastic masking films. Thus, it is difficult and expensive to construct the numerical model for such a multilayered structure. In this paper, we employ efficient numerical modeling of multilayered structures with a nonlocal ghost interlayer model in which ghost particles are distributed between main layers and they are interacting with each other in peridynamic way. We also consider a simple nonlocal contact condition in peridynamic frameworks to solve impact and penetration of the high-velocity impactor to the multilayered structure. Finally we can confirm the fracture capabilities of the method using a multilayered glass model in which 7 glass layers and a single elastic backing layer are affixed by polyvinyl butyral films.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.33 no.3
• /
• pp.171-178
• /
• 2020

In this study, we develop MPI-OpenMP hybrid parallelization for multibody peridynamic simulations. Peridynamics is suitable for analyzing complicated dynamic fractures and various discontinuities. However, compared with a conventional finite element method, nonlocal interactions in peridynamics cost more time and memory. In multibody peridynamic analysis, the costs increase due to the additional interactions that occur when computing the nonlocal contact and ghost interlayer models between adjacent bodies. The costs become excessive when further refinement and smaller time steps are required in cases of high-velocity impact fracturing or similar instances. Thus, high computational efficiency and performance can be achieved by parallelization and optimization of multibody peridynamic simulations. The analytical code is developed using an Intel Fortran MPI compiler and OpenMP in NURION of the KISTI HPC center and parallelized through MPI-OpenMP hybrid parallelization. Further parallelization is conducted by hybridizing with OpenMP threads in each MPI process. We also try to minimize communication operations by model-based decomposition of MPI processes. The numerical results for the impact fracturing of multiple bodies show that the computing performance improves significantly with MPI-OpenMP hybrid parallelization.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.1
• /
• pp.37-44
• /
• 2019

A bond-based peridynamic model has been reported dynamic fracture characteristic of brittle materials through a simple constitutive model. In the model, each bond is assumed to be a simple spring operating independently. As a result, this simple bond interaction modeling restricts the material behavior having a fixed Poisson's ratio of 1/4 and not being capable of expressing shear deformation. We consider a state-based peridynamics as a generalized peridynamic model. Constitutive models in the state-based peridynamics are corresponding to those in continuum theory. In state-based peridynamics, thus, the response of a material particle depends collectively on deformation of all bonds connected to other particles. So, a state-based peridynamic theory can represent the volume and shear changes of the material. In this paper, the perfect plasticity is considered to express plastic deformation of material by the state-based peridynamic constitutive model with perfect plastic flow rule. The elastic-plastic behavior of the material is verified through the stress-strain curves of the flat plate example. Furthermore, we simulate the high-speed impact on 3D granite model with a nonlocal contact modeling. It is observed that the damage patterns obtained by peridynamics are similar to experimental observations.