• Computational Structural Engineering
• /
• v.26 no.4
• /
• pp.19-24
• /
• 2013

• Computational Structural Engineering
• /
• v.26 no.4
• /
• pp.25-30
• /
• 2013

• Computational Structural Engineering
• /
• v.31 no.1
• /
• pp.23-31
• /
• 2018

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.2
• /
• pp.209-218
• /
• 2017

Sea ice is the main factor affecting the safety of the Arctic engineering. However, traditional numerical methods derived from classical continuum mechanics have difficulties in resolving discontinuous problems like ice damage. In this paper, a non-local, meshfree numerical method called "peridynamics", which is based on integral form, was applied to simulate the interaction between level ice and a cylindrical, vertical, rigid structure at different velocities. Ice in the simulation was freshwater ice and simplified as elastic-brittle material with a linear elastic constitutive model and critical equivalent strain criterion for material failure in state-based peridynamics. The ice forces obtained from peridynamic simulation are in the same order as experimental data. Numerical visualization shows advantages of applying peridynamics on ice damage. To study the repetitive nature of ice force, damage zone lengths of crushing failure were computed and conclude that damage zone lengths are 0.15-0.2 times as ice thickness.

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

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

• Coupled systems mechanics
• /
• v.4 no.1
• /
• pp.1-18
• /
• 2015

The state-based peridynamics is considered a nonlocal method in which the equations of motion utilize integral form as opposed to the partial differential equations in the classical continuum mechanics. As a result, the enforcement of boundary conditions in solid mechanics analyses cannot follow the standard way as in a classical continuum theory. In this paper, a new approach for the boundary condition enforcement in the state-based peridynamic formulation is presented. The new method is first formulated based on a convex kernel approximation to restore the Kronecker-delta property on the boundary in 1-D case. The convex kernel approximation is further localized near the boundary to meet the condition that recovers the correct boundary particle forces. The new formulation is extended to the two-dimensional problem and is shown to reserve the conservation of linear momentum and angular momentum. Three numerical benchmarks are provided to demonstrate the effectiveness and accuracy of the proposed approach.

• Interaction and multiscale mechanics
• /
• v.7 no.1
• /
• pp.525-542
• /
• 2014

The state-based peridynamics is considered a nonlocal method in which the equations of motion utilize integral form as opposed to the partial differential equations in the classical continuum mechanics. As a result, the enforcement of boundary conditions in solid mechanics analyses cannot follow the standard way as in a classical continuum theory. In this paper, a new approach for the boundary condition enforcement in the state-based peridynamic formulation is presented. The new method is first formulated based on a convex kernel approximation to restore the Kronecker-delta property on the boundary in 1-D case. The convex kernel approximation is further localized near the boundary to meet the condition that recovers the correct boundary particle forces. The new formulation is extended to the two-dimensional problem and is shown to reserve the conservation of linear momentum and angular momentum. Three numerical benchmarks are provided to demonstrate the effectiveness and accuracy of the proposed approach.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.6
• /
• pp.637-643
• /
• 2011

The bond-based peridynamic model is able to capture many of the essential characteristics of dynamic brittle fracture observed in experiments: crack branching, crack-path instability, asymmetries of crack paths, successive branching, secondary cracking at right angles from existing crack surfaces, etc. In this paper we investigate the influence of the stress waves on the crack branching angle and the velocity profile. We observe that crack branching in peridynamics evolves as the phenomenology proposed by the experimental evidence: when a crack reaches a critical stage(macroscopically identified by its stress intensity factor) it splits into two or more branches, each propagating with the same speed as the parent crack, but with a much reduced process zone.