• 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.27 no.4
• /
• pp.297-303
• /
• 2014

Using the bond-based peridynamics and the parallel computation with binary decomposition, an adjoint shape design sensitivity analysis(DSA) method is developed for the dynamic crack propagation problems. The peridynamics includes the successive branching of cracks and employs the explicit scheme of time integration. The adjoint variable method is generally not suitable for path-dependent problems but employed since the path of response analysis is readily available. The accuracy of analytical design sensitivity is verified by comparing it with the finite difference one. The finite difference method is susceptible to the amount of design perturbations and could result in inaccurate design sensitivity for highly nonlinear peridynamics problems with respect to the design. It turns out that $C^1$-continuous volume fraction is necessary for the accurate evaluation of shape design sensitivity in peridynamic discretization.

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

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2011.04a
• /
• pp.561-564
• /
• 2011

다양한 공학/산업적 측면에서 동적 취성 파괴 현상은 매우 중요하다. 취성 균열은 다른 균열 전파에 비해 그 전파 속도가 매우 빠르고 전파 범위가 넓기 때문에 대규모의 파괴 현상을 일으킨다. 동적 전파 중인 취성 균열 거동을 모델화하기 위해 오랜 기간 동안 많은 연구가 진행되었지만, 여전히 많은 부분들이 해석되지 못한 채 남아있다. 특히 균열 생성 및 전파를 위해 인위적인 조건들을 도입해야 하는 것은 기존 방법론들이 가지는 공통적인 문제점이다. 본 연구는 peridynamics를 동적 분기 균열 문제 해석에 도입한다. Peridynamics는 전통적인 연속체 이론에 기반한 수치해석 모델화 기법으로 균열과 같은 비연속성이 있는 문제의 모델화에 강점이 있으며, 인위적인 조건 없이 매우 간단한 방법으로 파괴 현상을 해석할 수 있다. 본 연구에서는 peridynamics 모델이 실험적으로 관측된 분기균열 형상과 균열 전파 속도를 매우 잘 예측해 낼 수 있음을 보인다. 또한 균열팁 주변에 높은 응력이 발생할 때 나타나는 연쇄 분기 현상도 해석할 수 있다. 이와 같은 연구를 통해 응력파가 균열 전파 속도를 변화시키고 전파 방향에도 영향을 주는 것을 알 수 있었다. 수치해석 결과도 또한 실험 결과들과 잘 부합함을 확인하였다.

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

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

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.4
• /
• pp.425-431
• /
• 2015

Based on a bond-based peridynamics theory for dynamic crack propagation problems, this paper presents a design sensitivity analysis and optimization method. Peridynamics has a peculiar advantage over the existing continuum theory in the mathematical modelling of problems where discontinuities arise. For the design optimization of the crack propagation problems, a non-shape design sensitivity is derived using the adjoint variable method. The obtained adjoint sensitivity of displacement and strain energy turns out to be very accurate and efficient compared to the finite different sensitivity. The obtained design sensitivities are futher utilized to optimally control the position of bifurcation point in the design optimization of crack propagation in a plate under tension. A numerical experiment demonstrates that the optimal distribution of material density could delay the position of bifurcation.