• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.4
• /
• pp.279-285
• /
• 2012

In this paper a computational method for a homogenized peridynamics description of unidirectional fiber-reinforced composites is presented. For these materials, dynamic brittle fracture and damage are simulated with the proposed peridynamic model. Compared with observations from dynamic experiments by Coker et al.(2001), the peridynamic computational model can reproduce various characteristics of dynamic fracture and supersonic or intersonic crack growth in asymmetrically loaded unidirectional fiber-reinforced composite plates. Also we analyze the same model in the symmetric loading condition and figure out that the asymmetric loading leads to a much higher propagation speed. Consistent results have been reported in the experiments.

• 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.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.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.

• Computers and Concrete
• /
• v.33 no.5
• /
• pp.605-619
• /
• 2024

The increase of reclaimed asphalt pavement (RAP) content in recycled asphalt concrete (RAC) is accompanied by the degradation of low-temperature cracking resistance, which has become an obstacle to the development of RAC. This paper aims to reveal the meso-scale mechanisms of the low-temperature fracture behavior of RAC and provide a theoretical basis for the economical recycling of RAP. For this purpose, micromechanical heterogeneous peridynamic model of RAC was established and validated by comparing three-point bending (TPB) test results against corresponding numerical simulation results of RAC with 50% RAP content. Furthermore, the models with different aggregate shapes (i.e., average aggregates circularity (${\bar{C_r}}=1.00$, 0.75, and 0.50) and RAP content (i.e., 0%, 15%, 30%, 50%, 75%, and 100%) were constructed to investigate the effect of aggregate shape and RAP content on the low-temperature cracking resistance. The results show that peridynamic models can accurately simulate the low-temperature fracture behavior of RAC, with only 2.9% and 13.9% differences from the TPB test in flexural strength and failure strain, respectively. On the meso-scale, the damage in the RAC is mainly controlled by horizontal tensile stress and the stress concentration appears in the interface transition zone (ITZ). Aggregate shape has a significant effect on the low-temperature fracture resistance, i.e., higher aggregate circularity leads to better low-temperature performance. The large number of microcracks generated during the damage evolution process for the peridynamic model with circular aggregates contributes to slowing down the fracture, whereas the severe stress concentration at the corners leads to the fracture of the aggregates with low circularity under lower stress levels. The effect of RAP content below 30% or above 50% is not significant, but a substantial reduction (16.9% in flexural strength and 16.4% in failure strain) is observed between the RAP content of 30% and 50%. This reduction is mainly attributed to the fact that the damage in the ITZ region transfers significantly to the aggregates, especially the RAP aggregates, when the RAP content ranges from 30% to 50%.

• Structural Engineering and Mechanics
• /
• v.75 no.3
• /
• pp.339-356
• /
• 2020

In this article, a developed bond-based peridynamic model for functionally graded materials (FGMs) is proposed to simulate the dynamic fracture behaviors in FGMs. In the developed bond-based peridynamic model for FGMs, bonds are categorized into three different types, including transverse directionally peridynamic bond, gradient directionally peridynamic bond and arbitrary directionally peridynamic bond, according to the geometrical relationship between directions of peridynamic bonds and gradient bonds in FGMs. The peridynamic micromodulus in the gradient directionally and arbitrary directionally peridynamic bonds can be determined using the weighted projection method. Firstly, the standard bond-based peridynamic simulations of crack propagation and branching in the homogeneous PMMA plate are performed for validations, and the results are in good agreement with the previous experimental observations and the previous phase-field numerical results. Then, the numerical study of crack initiation, propagation and branching in FGMs are conducted using the developed bond-based peridynamic model, and the influence of gradient direction on the dynamic fracture behaviors, such as crack patterns and crack tip propagation speed, in FGMs is systematically studied. Finally, numerical results reveal that crack branching in FGMs under dynamic loading conditions is easier to occur as the gradient angle decreases, which is measured by the gradient direction and direction of the initial crack.

• Geomechanics and Engineering
• /
• v.15 no.5
• /
• pp.1047-1059
• /
• 2018

This paper focuses on the cracking and fragmentation process in rock materials containing a pair of non-parallel flaws, which are through the specimen thickness, under vertical compression. Several numerical experiments are conducted with varying flaw arrangements that affect the initiation and tensile wing cracks, shear crack growth, and crack coalescing behaviors. To obtain realistic numerical results, a parallelized peridynamics formulation coupled with a finite element method, which is able to capture arbitrarily occurring cracks, is employed. From previous studies, crack initiation and propagation of tensile wing cracks, horsetail cracks, and anti-wing cracks are well understood along with the coalescence between two parallel flaws. In this study, the coalescence behaviors, their fragmentation sequences, and the role of an x-shaped shear band in rock material containing two non-parallel flaws are discussed in detail on the basis of simulation results strongly correlated with previous experimental results. Firstly, crack initiation and propagation of tensile wing cracks and shear cracks between non-parallel flaws are investigated in time-history and then sequential coalescing behavior is analyzed. Secondly, under the effect of varying inclination angles of two non-parallel flaws and overlapping ratios between a pair of non-parallel flaws, the cracking patterns including crack coalescence, fragmentation, and x-shaped shear band are investigated. These numerical results, which are in good agreement with reported physical test results, are expected to provide insightful information of the fracture mechanism of rock with non-parallel flaws.

• 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.