Capturing Finite Target Dynamics: Phase-Delayed Analytic Modeling of Multi-Layer Penetration Events

기술 요약: 위상 지연 해석 모델링을 통한 유한 표적 동역학 포착

문제 정의
Walker-Anderson 반무한 공간 관통 모델은 두꺼운 표적 내 강체 및 침식되는 로드(rod)의 관통 깊이를 계산하는 효율적인 도구로 오랫동안 사용되어 왔다. 그러나 이 모델을 유한한 표적으로 확장하려는 기존 시도들은 특히 얇은 표적에서 코-테일(nose-tail) 속도 프로파일을 정확하게 포착하는 데 어려움을 겪고 있다. 다층 구조의 얇은 벽 표적 스택에서 이러한 결함은 상세한 유체역학 시뮬레이션이나 실험적 예측과 비교했을 때 로드 침식 특성을 잘못 규정하는 결과를 초래한다. 구체적으로, 이전의 유한 영역 모델(Ravid 및 Walker의 모델 등)은 파동 전파 효과가 표적 두께를 가로지르는 데 필요한 시간을 고려하지 못하여, 표적 약화 및 투사체 회복 단계를 조기에 트리거하는 경성을 보인다.

방법론
저자들은 상세한 ALE3D 유체역학 시뮬레이션으로부터 얻은 통찰을 통합하여 유한 표적 내의 파동 전파를 올바르게 설명하는 업데이트된 Walker-Anderson 모델 수정안을 제안한다. 방법론은 세 가지 주요 기술적 개발을 포함한다:

1. 업데이트된 초기 조건: 본 논문은 반구형 헤드 투사체의 시작 조건을 수정한다. 전통적인 1차원 충격 점프 조건($P_{shock}$)에만 의존하는 대신, 저자들은 초기 충격을 위한 "정체 충격(stagnation shock)" 조건($P_{stag}$)을 도입한다. 이 조건은 표적이 더 빠르게 정체 흐름 형태를 채택한다고 가정하며, 원래의 충격 방정식에 비해 1/2의 계수를 도입한다. 이는 깨끗한(pristine) 충격 시 인터페이스 속도를 더 잘 초기화하며, 다중 표적 시나리오에서의 후속 충격에는 전통적인 $P_{shock}$ 조건이 유지된다.
2. 파동 전파 분석: 저자들은 ALE3D 시뮬레이션의 후처리 과정을 통해 유한 표적 내의 속도 및 소성 특성의 진화를 분석하였다. 저자들은 충격파 속도가 초기에는 과구동(over-driven)되지만, 소성 효과의 전파는 표적의 벌크 음속($c_0$)에 의해 가장 잘 근사된다는 것을 밝혀냈다. 이 분석을 통해 표적의 유한 효과(예: 배면 팽창 및 약화)가 투사체 인터페이스에 즉각적으로 영향을 미치는 것이 아니라, 파동 통과 시간에 의해 지연된다는 것을 확인하였다.
3. 위상 지연 결합(Phase-Delayed Coupling): 핵심 기여는 유한 표적 모델에 위상 지연을 구현하는 것이다. 저자들은 1차원 파동 방정식의 특성 형태를 활용하여 표적의 전면과 후면을 연결한다. 모델은 전방($t_+$) 및 후방($t_-$) 특성 시간을 계산하여, 인터페이스 항에 배면 특성(예: 배면 속도 $u_b$ 및 근접도 파라미터 $\tilde{\alpha}$)의 적용을 지연시킨다. 모델은 파동이 표적 두께를 통과할 충분한 시간이 흐를 때까지 무한 표적 가정에서 유한 표적 가정으로 전환되지 않도록 함으로써, 물리적으로 정당화되기 전까지 투사체 인터페이스가 배면 연화(back-face softening)의 영향을 받지 않도록 보장한다.

주요 기여

• 위상 지연 메커니즘: 벌크 파동 전파에 기반한 위상 지연의 도입은 해석 모델이 표적 약화의 시작 시점을 정확하게 맞출 수 있게 한다. 이는 이전의 유한 영역 모델에서 관찰되었던 투사체 코(nose)의 조기 재가속 현상을 방지한다.
• 추가 파라미터 없음: 업데이트된 모델은 새로운 튜닝 가능한 모델 파라미터를 요구하지 않으면서도, 얇은 표적에 대한 정확도를 개선하는 동시에 두꺼운 표적에 대해서도 양호한 거동을 유지한다. 이러한 수정 사항은 기존의 재료 물성(벌크 음속)과 충격 기하학에 의존한다.
• 개선된 초기 조건: 깨끗한 충격에 대한 정체 충격 초기화와 후속 충격에 대한 표준 충격 초기화의 구분은 더 정확한 초기 시간 속도 프로파일을 제공한다.
• 실험 검증: 본 모델은 2.90 cm 강철로 구성된 두꺼운 유한 표적에 대한 실험 데이터로 검증되었으며, 명시적인 투사체 회복 단계를 생략하고 자연스러운 운동량 균형에 의존하는 것이 강제적인 회복 트리거를 사용하는 모델보다 더 나은 일치도를 보임을 입증하였다.

결과
업데이트된 모델은 다양한 시나리오에 걸쳐 상세한 ALE3D 유체역학 시뮬레이션과 벤치마킹되었다:

• 얇은 표적: 0.5 cm 및 1 cm 두께의 표적에 대해, 새 모델은 기존의 Walker 배면 모델보다 성능이 현저히 우수하다. 새 모델은 코 속도 회복의 지연과 지속적인 표적 강도를 정확하게 포착하는 반면, 기존 모델은 조기 약화를 예측하였다.
• 다중 표적 충격: 6개의 4 cm 플레이트 스택이 포함된 시뮬레이션에서, 업데이트된 모델은 초기 반응과 지속적인 투사체 침식 모두에 대해 유체역학 시뮬레이션과 우수한 일치도를 보였다. 또한 특정 파단 변형률 기준에 대한 민감도 없이 표적 약화의 타이밍과 그에 따른 투사체 속도 프로파일을 성공적으로 포착하였다.
• 속도 범위: 모델은 파라미터 재조정 없이 1200 m/s에서 3000 m/s 사이의 넓은 충격 속도 범위에서 높은 일치도를 유지하였다.
• 실험 상관관계: 2.90 cm 강철 표적에 대한 실험적 출구 속도 데이터와의 비교를 통해, 특히 V50 영역 통과 및 회복 후 속도를 포착하는 데 있어 모델의 효용성을 확인하였다.

의의
본 논문은 이러한 수정 사항이 훨씬 더 많은 계산 자원과 전문가의 숙련도를 요구하는 고충실도 유체역학 시뮬레이션과 비교하여 "놀라운 정확도"를 제공한다고 주장한다. 상세한 시뮬레이션을 활용하여 반해석적(semi-analytic) 방법을 정교화함으로써, 저자들은 해석 모델의 속도와 단순성을 유지하면서도 주요 실패 모드인 얇은 표적 역학을 교정할 수 있는 도구를 제공한다. 이 연구는 기존의 Walker 무한 공간 모델 가정을 재해석할 수 있음을 시사한다: 즉, 모델은 투사체 인터페이스에 즉시 적용될 수 있지만(초기 시간의 즉각성 정당화), 표적 깊이 내의 실제적인 속도장 발달을 고려하기 위해 유한 배면 조건의 적용은 지연되어야 한다. 이 접근 방식은 엔지니어들이 전체 유체역학 코드를 사용하는 계산적 부담 없이도 고속 관통 이벤트에서의 결합된 충격체-표적 상호작용을 더 잘 이해할 수 있게 해준다.