Peridynamic modeling of the crack velocity dependence via an incubation time fracture criterion

본 연구는 Homalite-100 에 대한 Ravi-Chandar 와 Knauss 의 실험을 모델링하기 위해 잠복 시간 파괴 기준을 적용한 페리다이나믹 접근법을 활용하여, 일정한 균열 속도에서 모드-I 응력확대계수의 변화와 더 높은 속도에서의 미세 분지 발생이 동적 파괴에서 균열 속도 의존성의 본질에 대한 새로운 통찰을 제공함을 규명하였다.

핵심

상상해 보세요. Homalite-100 시트와 같은 취성 플라스틱 조각 위를 균열이 질주하는 모습을 말입니다. 물리학의 옛날에는 과학자들이 균열의 이동 속도를 알면 그 균열을 밀어내는 정확한 '응력' (또는 압력) 을 계산할 수 있다고 생각했습니다. 마치 자동차가 시속 60 마일로 달리고 있다면 엔진이 정확히 100 마력을 생산하고 있어야 한다고 생각하는 것과 같았습니다. 간단하죠?

하지만 1980 년대의 실험들은 이것이 사실이 아님을 보여주었습니다. 때로는 균열이 정확히 같은 속도로 이동하고 있었지만, 그것을 밀어내는 압력은 극적으로 달랐습니다. 마치 두 대의 자동차가 모두 시속 60 마일로 달리고 있었지만, 한 대는 작은 엔진을 가지고 있었고 다른 한 대는 로켓 부스터를 장착하고 있었던 것과 같았습니다. 과학자들은 당황했습니다: 왜 같은 속도가 서로 다른 '밀어내는 힘'을 가질까요?

이 논문은 저자들이 새로운 종류의 컴퓨터 시뮬레이션을 사용하여 이 미스터리를 해결하는 탐정 이야기입니다.

탐정 도구: 퍼리다이나믹스 (Peridynamics)

균열에 대한 대부분의 컴퓨터 모델은 도미노 사슬과 같습니다. 한 도미노가 쓰러지면 다음 도미노를 밀어냅니다. 하지만 도미노 하나가 사라지면 (균열이 생기면) 사슬이 끊어지고 수학이 멈춰 버립니다.

저자들은 퍼리다이나믹스라는 방법을 사용했습니다. 이를 사슬이 아니라 벌떼로 생각하세요. 각 벌은 중간에 간격이 있더라도 일정 거리 내의 다른 모든 벌과 대화할 수 있습니다. 한 벌이 날아가 버리면 (균열이 형성되면) 다른 벌들은 그 벌과 대화를 멈추지만, 나머지 벌떼는 완벽하게 움직임을 계속합니다. 이를 통해 컴퓨터는 혼란에 빠지지 않고 파괴와 균열을 처리할 수 있습니다.

비밀 재료: '부화 시간'

이 논문에서의 진정한 돌파구는 균열이 실제로 언제 깨져야 하는지를 결정하는 방식에 있습니다.

옛날 방식에서는 압력이 충분히 높아지면 재료가 즉시 깨졌습니다. 하지만 저자들은 **부화 시간 기준 (Incubation Time Criterion)**이라는 규칙을 사용했습니다.

마른 나뭇가지를 꺾으려 한다고 상상해 보세요. 당기자마자 바로 꺾이지는 않습니다. 당기고, 섬유가 늘어나고 약해지는 동안 찰나의 순간을 유지한 뒤, 그제야 꺾입니다. 그 찰나의 순간이 바로 '부화 시간'입니다.

저자들은 컴퓨터 벌떼에게 압력의 마지막 몇 마이크로초를 기억하도록 프로그래밍했습니다. 재료는 그 짧은 '부화' 기간 동안의 평균 압력이 충분히 높을 때만 깨집니다. 이는 재료가 깨지기로 '결정'하는 데는 아주 조금의 시간이 필요하다는 사실을 반영한 것입니다.

그들이 발견한 것

그들은 실제 실험과 마찬가지로 플라스틱 판을 잡아당기는 시뮬레이션을 실행했습니다. 그들이 발견한 것은 다음과 같습니다:

1. 속도 대 압력 퍼즐: 실제 실험과 마찬가지로, 그들의 컴퓨터는 동일한 균열 속도에 대해 압력 (응력 강도 계수) 이 단일 숫자가 아님을 보여주었습니다. 그것은 일정한 범위를 가졌습니다. 때로는 낮았고 때로는 높았습니다.
2. '미세 분기' 효과: 균열이 천천히 이동할 때는 직진했습니다. 하지만 속도가 빨라져서 (초당 400 미터 이상) 떨리기 시작했습니다. 그것은 나뭇가지가 작은 가지로 갈라지듯 미세한 측면 균열을 틔우기 시작했습니다.
   • 비유: 스프린트를 하는 러너를 상상해 보세요. 일정한 조깅을 할 때는 직선으로 달립니다. 하지만 최고 속도로 스프린트를 할 때는 균형을 유지하기 위해 약간 비틀거리고 지그재그로 움직이기 시작합니다.
   • 결과: 이러한 작은 '비틀거림' (미세 분기) 이 압력 수치를 극적으로 오르내리게 만들었습니다. 이것이 주어진 속도에 대해 압력이 유일하지 않았던 이유를 설명해 주었습니다. 균열은 질주하면서 물리적으로 그 모양을 약간 변화시키고 있었던 것입니다.

결론

이 논문은 동일한 균열 속도에 대해 서로 다른 압력 값을 보는 이유는 균열이 매끄럽고 완벽한 선이 아니기 때문이라고 결론 내립니다. 그것은 요동치는 혼란스러운 살아 있는 존재입니다.

• 낮은 속도에서: 균열은 안정적이며 압력도 상대적으로 안정적입니다.
• 높은 속도에서: 균열은 '미세 분기' (작은 측면 균열을 틔움) 를 시작합니다. 이 혼란은 압력을 요동치게 만들어 실험에서 관찰된 산포를 발생시킵니다.

이 '벌떼' (퍼리다이나믹스) 와 '기다림 기간' (부화 시간) 을 결합함으로써 저자들은 수십 년간 실제 실험에서 보여준 균열 속도와 압력 사이의 messy(불규칙한) 그리고 유일하지 않은 관계를 성공적으로 재현했습니다. 그들은 데이터의 '노이즈'가 실수가 아니라, 빠르게 이동하는 균열이 어떻게 행동하는지에 대한 실제 물리적 특징임을 증명했습니다.

기술 요약

기술 요약: 잠복 시간 파괴 기준을 통한 균열 속도 의존성에 대한 퍼리동역학 모델링

문제 제기
본 연구는 동적 파괴 역학의 핵심적인 미해결 문제인 순간 모드 I 응력 확대 계수 (SIF, $K_I$) 와 균열 전파 속도 ($\dot{a}$) 간의 의존성을 다룬다. 고전적인 선형 탄성 파괴 역학 (LEFM) 과 Freund 및 Rosakis 기반의 동적 기준들은 종종 $K_I$와 $\dot{a}$ 사이에 고유한 관계가 존재함을 시사한다. 그러나 Ravi-Chandar 와 Knauss 가 수행한 취성 Homalite-100 고분자 판에 대한 실험 데이터는 단일 균열 속도가 SIF 값의 범위에 대응되는 비고유적 의존성을 보여준다. furthermore, 더 높은 속도에서는 미세 분기 (micro-branching) 가 발생하여 SIF 측정치에 상당한 산란을 초래한다. 과제는 경험적으로 조정된 속도 의존성 재료 물성에 의존하지 않고 이러한 비고유성과 관련된 변동을 수치적으로 재현하는 데 있다.

방법론
저자들은 Ravi-Chandar 와 Knauss 의 실험 조건을 재현하여 Homalite-100 에서의 동적 균열 전파를 시뮬레이션하기 위해 비국소적 퍼리동역학 (peridynamic) 프레임워크를 활용한다. 방법론적 핵심 혁신은 퍼리동역학 공식화에 잠복 시간 파괴 기준 (ITFC) 을 통합한 것이다.

• 퍼리동역학 프레임워크: 시뮬레이션은 오픈소스 소프트웨어 Peridigm 에 구현된 일반 상태 기반 (ordinary state-based) 퍼리동역학 모델을 사용한다. 재료는 구형 호라이즌 ($\delta$) 내에서 상호작용하는 점들로 이산화되어, 공간 미분 연산자 없이 균열과 같은 불연속성을 자연스럽게 처리할 수 있게 한다.
• 잠복 시간 파괴 기준 (ITFC): 고전적 기준이 순간적인 파괴를 가정하는 것과 달리, ITFC 는 파괴가 미세균열과 같은 준비 과정에 유한한 시간 ($\tau$) 을 필요로 한다고 주장한다. 이 기준은 잠복 기간 $\tau$와 특징 길이 $d$에 걸쳐 시간 평균 응력을 평가한다.
  • 특징 길이 $d$는 임계 정적 SIF ($K_{Ic}$) 와의 일관성을 보장하기 위해 준정적 Irwin 기준에서 유도된다.
  • 퍼리동역학 구현은 원격 응력 (RS) 파괴 기준을 적용한다. 결합 파괴는 잠복 시간 $\tau$에 걸쳐 수직 결합 응력을 적분하여 결정된다 (식 11).
  • 계산 복잡성을 관리하고 전체 거시 분기보다는 미세 분기에 초점을 맞추기 위해, 기준은 주 균열 축에 평행 ($0^\circ$) 하고 수직 ($90^\circ$) 인 평면에서만 평가된다.
• 시뮬레이션 설정: 모델은 50 mm 프리노치를 가진 300 mm $\times$ 500 mm 시편을 시뮬레이션한다. 동적 하중은 25 $\mu$s 상승 시간을 가진 사다리꼴 체적 힘 펄스를 통해 인가된다. 동적 SIF 는 동적 J 적분을 위한 이동 직사각형 윤곽선을 사용하여 계산되는데, 이는 적분 경로가 주 균열 선단을 포함하되 발달 중인 미세 분기를 제외하도록 보장한다.

주요 결과
수치 시뮬레이션은 $K_I$-$\dot{a}$ 관계와 관련하여 실험적으로 관찰된 현상을 성공적으로 재현한다:

1. 비고유적 $K_I$-$\dot{a}$ 의존성: 결과는 균열 전파 속도가 거의 일정할 때 순간 SIF 가 상당한 변동을 보임을 확인한다. 결과적으로, 모델은 실험에서 관찰된 비고유성을 반영하여 주어진 속도에 대한 SIF 값의 범위를 예측한다.
2. 속도 안정화 및 과도 현상: 실험에서는 균열 속도가 거의 즉시 안정화되는 것으로 나타나지만, 시뮬레이션은 정상 상태에 도달하기 전에 약 30–40 $\mu$s 동안 지속되는 초기 가속 단계를 보인다. 이 과도 단계 동안 SIF 이력은 뚜렷한 변동을 나타낸다.
3. 미세 분기 및 산란: 더 높은 하중률 (속도 > 400 m/s 로 이어짐) 에서 시뮬레이션은 미세 분기의 시작을 보인다. 이 현상은 SIF 진폭의 현저한 증가와 상관관계가 있다. 고속에서의 SIF 값에 대한 수치적 산란 (0.53–1.25 MPa$\sqrt{m}$) 은 실험 범위 (0.45–1.5 MPa$\sqrt{m}$) 와 합리적으로 일치한다.
4. 메커니즘 검증: 본 연구는 관찰된 SIF 산란이 측정 오차의 인공물이 아니라, 균열 선단의 역학, 특히 가속 단계와 미세 분기의 시작에 의해 물리적으로 정당화됨을 보여준다.

의의 및 주장
본 논문은 ITFC 의 퍼리동역학 구현이 광범위한 하중률에 걸쳐 동적 파괴 인성을 평가하고 균열 전파를 모델링하는 가치 있는 접근법을 제공한다고 주장한다. 주요 의의는 임의의 속도 의존성 재료 매개변수를 도입하지 않고도 취성 파괴의 속도 민감성과 $K_I$-$\dot{a}$ 관계의 비고유성을 포착할 수 있는 모델의 능력에 있다.

저자들은 모델이 실험적 경향을 질적으로 포착하지만, 속도 안정화의 타이밍과 SIF 변동의 정확한 크기 측면에서 정량적 불일치가 남아 있음을 겸손하게 지적한다. 그들은 비고유적 의존성의 성공적인 예측이 응력의 시간적 이력을 고려하는 잠복 시간의 포함과 동적 파괴와 관련된 특이장 및 불연속성을 자연스럽게 처리하는 퍼리동역학의 비국소적 성질에 기인한다고 본다. 연구는 균열 가속 단계, SIF 변동, 그리고 미세 분기 시작 사이의 상호작용이 실험적으로 관찰된 $K_I$-$\dot{a}$ 거동 뒤에 있는 지배적 메커니즘이라고 결론 내린다.