• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.1
• /
• pp.17-26
• /
• 2014

This paper presents a dynamic crack propagation algorithm based on the Moving Least Squares(MLS) difference method. The derivative approximation for the MLS difference method is derived by Taylor expansion and moving least squares procedure. The method can analyze dynamic crack problems using only node model, which is completely free from the constraint of grid or mesh structure. The dynamic equilibrium equation is integrated by the Newmark method. When a crack propagates, the MLS difference method does not need the reconstruction of mode model at every time step, instead, partial revision of nodal arrangement near the new crack tip is carried out. A crack is modeled by the visibility criterion and dynamic energy release rate is evaluated to decide the onset of crack growth together with the corresponding growth angle. Mode I and mixed mode crack propagation problems are numerically simulated and the accuracy and stability of the proposed algorithm are successfully verified through the comparison with the analytical solutions and the Element-Free Galerkin method results.

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

• Computational Structural Engineering
• /
• v.10 no.3
• /
• pp.22-29
• /
• 1997

본 논문에서는 균열이 존재하는 구조부재에 충격이나 폭발하중이 가해진 경우 동적응력확대계수를 구하는 방법들은 논의하고 특히 코오스틱 실험법 및 수치적으로 코오스틱 곡선을 구하여 동적응력확대계수를 구하는 과정을 자세히 설명하였다. 폭발 및 충격에 의한 구조물의 파괴해석은 이와 같은 하중을 받는 압력용기, 빌딩, 초고속선, 해군 함정 등의 파괴강도설계 및 안전성 평가에 핵심기술로 대두되고 있으며 또한 우주항공산업, 고속전철, 암반역학 등의 여러 분야에서 중요한 의미를 갖는다. 따라서 앞으로도 균열진전 및 정지조건, 탄소성 동적파괴해석 및 재료의 충격거동 등에 대한 연구들이 계속되어져야 할 것으로 사료된다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.6
• /
• pp.583-590
• /
• 2016

This paper presents a dynamic crack propagation algorithm with Rayleigh damping effect based on the MLS(Moving Least Squares) Difference Method. Dynamic equilibrium equation and constitutive equation are derived by considering Rayliegh damping and governing equations are discretized by the MLS derivative approximation; the proportional damping, which has not been properly treated in the conventional strong formulations, was implemented in both the equilibrium equation and constitutive equation. Dynamic equilibrium equation including time relevant terms is integrated by the Central Difference Method and the discrete equations are simplified by lagging the velocity one step behind. A geometrical feature of crack is modeled by imposing the traction-free condition onto the nodes placed at crack surfaces and the effect of movement and addition of the nodes at every time step due to crack growth is appropriately reflected on the construction of total system. The robustness of the proposed numerical algorithm was proved by simulating single and multiple crack growth problems and the effect of proportional damping on the dynamic crack propagation analysis was effectively demonstrated.

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

• 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.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 Korea Concrete Institute
• /
• v.13 no.6
• /
• pp.637-644
• /
• 2001

A regularization method based on the Duvaut-Lions viscoplastic scheme for plastic-damage and continuum damage models, which provides mesh-independent and well-posed solutions in nonlinear earthquake analysis of concrete dams, is presented. A plastic-damage model regularized using the proposed rate-dependent viscosity method and its original rate-independent version are used for the earthquake damage analysis of a concrete dam to analyze the effect of the regualarization and mesh. The computational analysis shows that the regularized plastic-damage model gives well-posed solutions regardless mesh size and arrangement, while the rate-independent counterpart produces mesh-dependent ill-posed results.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.6 no.5
• /
• pp.45-51
• /
• 2002

The characteristics of nonlinear seismic behavior and failure mechanism of RC space frame in railroad viaducts have been studied by the numerical analysis in time domain. The structure concerned is modeled in 3 dimensional extent and the RC frame elements consisting of fibers are employed for the columns. The fibers are characterized as RC zone and PC one to distinguish the different energy release after cracking resulted from the bond characteristic between concrete and re-bar. Due to the deviation of the mass center and the stiffness center of the entire structure the complex behavior is shown under seismic actions. The excessive shear force is concentrated on the column beside flexible one relatively, which leads to the failure of bridge concerned.

• Tunnel and Underground Space
• /
• v.9 no.1
• /
• pp.64-71
• /
• 1999

A numerical analysis was performed to investigate the effects of notched blasthole in controlling the fracture plane. Analyzed were elastic and elasto-plastic response of rock, and fracture propagation under static and dynamic load conditions. Results showed that the region exceeding the tensile strength extended up to three times the radius of a normal blasthole in elastic analysis, while fifteen times in elasto-plastic analysis. It was shown that a crack was driven from the notch tip up to the distance of 23 times the hole radius in the case of a notched blasthole with a notch of 5 mm in depth and 30 mm in length.