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.

상태 기반 페리다이나믹 모델은 일반적인 재료 구성 모델을 구현할 수 있고 비국부 영역 내에서 연결된 모든 결합의 변형을 통해 각 절점의 재료 응답이 결정되기 때문에 체적 및 전단 변형을 모두 표현할 수 있다. 따라서 상태 기반 모델은 복잡한 동적 취성 파괴 현상(분기균열, 2차 균열, 계단균열, 균열 유착 등)을 해석하는데 유용하다. 본 논문에서는 평면응력 탄성체에 대해 2차원 상태 기반 페리다이나믹 모델을 적용하고 에너지해방율과 페리다이나믹 에너지 포텐셜로부터 손상 모델을 구성하였다. 페리다이나믹 파괴 해석 모델을 통해 취성 유리 재료에 대해 균열 면에 평행한 압축 응력파가 균열 분기 패턴에 미치는 영향에 대해 조사하였다. 실험을 통해 관찰된 동적 균열 진전 및 분기 패턴에 대한 주요 특성들이 페리다이나믹 해석을 통해 확인되었다. 또한 강한 인장 하중 하의 계단균열과 이차균열 등이 상태 기반 페리다이나믹 시뮬레이션을 통해 잘 모사가 되는 것을 확인할 수 있었다.

본 논문은 MLS(Moving Least Squares) 차분법을 바탕으로 동적균열전파 해석을 수행하기 위한 알고리즘을 제시한다. MLS 차분법은 절점만으로 이루어진 수치모델을 사용하며, 이동최소제곱법을 이용하여 전개한 Taylor 다항식을 기초로 미분근사식을 유도하기 때문에, 요소망의 제약에서 완벽하게 벗어난 절점해석이 가능하다. 시간항을 포함하는 동적 평형방정식은 Newmark 방법으로 시간적분 하였다. 동적하중을 받는 균열이 전파할 때, 매 시간단계마다 절점모델을 재구성하지 않고 균열선단 주변에서 국부적인 수정을 통해 해석이 가능하다. 동적균열을 묘사하기 위해 가시한계법(visibility criterion)을 적용하였고, 동적 에너지해방률을 산정하여 균열의 진전유무와 그에 상응하는 진전방향을 결정하였다. 모드Ⅰ 상태와 혼합모드 상태에서 균열이 진전하는 현상을 모사하였고, 이론해와 Element-Free Galerkin법으로 계산한 결과와의 비교를 통해 개발된 알고리즘의 정확성과 안정성을 검증하였다.

이 논문에서는 단일방향 및 크로스-플라이 섬유강화 복합재에서의 동적 균열 전파를 모사하기 위해 특별히 고안된 스펙트랄 방법의 정식화와 수치적 구현 방법에 대해 제시한다. 이 방법은 균열면에서 작용하는 힘과 변위 사이의 스펙트랄 관계식에 기초하고 있고, 재료는 횡등방성 고체로 가정된다. 본 논문에서 제안된 방법에 의해 섬유보강 복합재에서 전파하는 균열문제가 대해 검토 해석되며, 실험 및 문헌에 나와 있는 결과와 비교된다. 이 방법은 FRP 보강 철근 콘크리트 구조물에서의 균열해석문제에 직접 적용이 가능하다.