A 2D discrete model with a phase-field regularisation applied to dynamic crack propagation

유리나 경화강과 같은 취성 재료가 깨질 때, 그것은 단순히 두 조각으로 툭 끊어지는 것이 아닙니다. 대신, 균열이 고체 내부를 가로질러 빠르게 진행하며, 종종 복잡한 파쇄망으로 갈라져 물체를 여러 조각으로 산산조각 냅니다. 이 균열이 정확히 어디로 이동할지, 그리고 얼마나 빨리 움직일지를 예측하는 것은 공학에서 가장 어려운 과제 중 중 하나입니다. 만약 이 파괴 과정을 시뮬레이션하는 데 사용되는 모델에 결함이 있다면, 그것은 다리가 안전한 방향으로 무너질 것이라고 제시할 수도 있지만 실제로는 치명적으로 붕괴할 수도 있습니다. 수십 년 동안 과학자들은 재료의 원자 구조를 모방하기 위해, 재료를 작은 점들이 스프링으로 연결된 격자로 취급하는 컴퓨터 모델을 구축하려고 노력해 왔습니다. 이러한 모델들은 강력하지만, 한 가지 기묘한 결함을 가지고 있습니다. 바로 균열이 가는 경로가 재료 자체의 물리 법칙이 아니라, 컴퓨터 격자가 어떻게 회전되어 있는지에 따라 달라진다는 점입니다. 이는 마치 실제 유리 조각에 생긴 균열이 엔지니어가 컴퓨터 화면을 돌렸다는 이유만으로 갑자기 방향을 바꾸는 것과 같습니다.

이 문제를 해결하기 위해, 스페인 네브리하 대학교의 연구원 마티아스 브라운(Matías Braun)은 동적 파괴를 시뮬레이션하는 새로운 방법을 개발했습니다. 그의 접근 방식은 전통적인 점 격자와 '위상장(phase-field)' 기법을 결합한 것입니다. 위상장 기법은 균열을 점에서 점으로 건너뛰는 날카로운 선이 아니라, 재료를 통해 퍼져 나가는 부드럽고 모호한 손상 영역으로 취급합니다. 이 두 가지 아이디어를 혼합함으로써, 브라운은 컴퓨터 격자의 방향과 상관없이 균열이 어떤 방향으로든 이동할 수 있는 모델을 만들었습니다. 그는 이 새로운 방법을 세 가지 뚜로 다른 시나리오, 즉 고속 충격을 받은 강판, 잡아당겨서 끊어질 때까지 노치가 있는 판, 그리고 내부에서 외부로 폭발하는 두꺼운 금속 실린더에 대해 테스트했습니다. 모든 경우에서, 새 모델은 실제 실험 결과와 일치하는 결과를 나타냈으며, 컴퓨터 격자로 인한 인위적인 편향을 제거하여 재료가 실제로 어떻게 파괴되는지에 대한 더 명확하고 신뢰할 수 있는 관점을 제공했습니다.

브라운이 해결한 문제의 핵심은 컴퓨터가 원자의 보이지 않는 세계를 처리하는 방식에 있습니다. 그의 초기 연구에서는 재료를 점과 점을 잇는 선으로 연결된 삼각형 격자로 표현하는 이산 모델을 사용했습니다. 재료에 응력이 가해지면 점들 사이의 연결이 끊어지며 균열을 시뮬레이션했습니다. 그러나 점들이 고정된 패턴으로 배열되어 있었기 때문에, 균열은 격자의 선을 따라가도록 강요되었습니다. 만약 격자가 회전하면, 재료 자체는 변하지 않았음에도 불구하고 균열의 경로가 바뀌었습니다. 이는 실제 균열이 컴퓨터 메쉬의 임의적인 각도가 아니라 재료의 특성에 의해 결정된 최소 저항 경로를 따라야 한다는 점에서 매우 중요한 문제입니다.

브라운의 해결책은 위상장 정규화(phase-field regularisation)를 도입하는 것이었습니다. 그는 특정 점들 사이의 연결을 끊도록 강제하는 대신, 각 점에 재료가 얼마나 손상되었는지를 나타내는 값을 할당했습니다. 이 값은 재료가 완벽하게 온전한 상태를 의미하는 0에서부터 완전히 파괴된 상태를 의미하는 1까지의 범위를 갖습니다. 균열은 이 값이 점진적으로 증가하는 영역이 되어, 톱니 모양의 선 대신 부드러운 전이를 만들어냅니다. 이를 통해 균열은 격자의 경직된 구조를 무시하고 어떤 방향으로든 휘어지거나 갈라질 수 있습니다. 결정적으로, 브라운은 이 시스템을 재료의 두 가지 근본적인 특성인 인장 강도(균열을 시작하는 데 필요한 힘)와 파괴 에너지(균열을 계속 성장시키는 데 필요한 에너지)를 사용하여 보정했습니다. 이 시뮬레이션을 실제 물리적 상수들에 고정함으로써, 모델의 결과가 수학이 아닌 재료에 의존하도록 보장했습니다.

이 새로운 방법이 효과적임을 증명하기 위해, 브라운은 세 가지 엄격한 테스트를 수행했습니다. 첫 번째는 미리 파인 노치가 있는 강판을 투사체가 고속으로 타격하는 고전적인 실험인 칼토프-윙클러(Kalthoff–Winkler) 테스트였습니다. 이 시나리오에서 균열은 노치로부터 특정 각도로 튀어나옵니다. 기존 모델들은 격자가 어떻게 회전되었는지에 따라 잘못된 각도를 예측하는 경우가 많았습니다. 브라운의 시뮬레이션에서 균열은 약 70도의 각도로 일관되게 이동했으며, 이는 실험 관찰 결과와 완벽하게 일치했습니다. 이 결과는 격자가 판과 정렬되어 있든, 15도 또는 30도 회전되어 있든 동일하게 유지되었습니다. 균열 속도 또한 재료 내 음속의 약 70%로 안정적으로 유지되었으며, 격자의 방향과는 독립적이었습니다.

두 번째 테스트는 인장 하에서의 단일 가장자리 노치 판(single-edge-notched plate)을 이용한 것으로, 균열이 일정한 힘으로 잡아당길 때 어떻게 천천히 성장하는지를 연구하는 설정입니다. 여기서 목표는 균열이 노치로부터 반대쪽 가장까지 직선으로 이동하는지 확인하는 것이었습니다. 전통적인 격자 모델은 격자의 정렬에 따라 균열이 지그재그로 움직이거나 경로를 벗어나도록 강제하는 경우가 많았습니다. 그러나 브라운의 모델은 격자의 방향에 관계없이 모든 방향에서 직선 균열을 생성했습니다. 판을 파괴하는 데 필요한 힘 또한 격자가 회전되더라도 2% 미만의 오차만을 보이며 일관되게 나타났습니다. 이는 모델이 컴퓨터 메쉬의 인위적인 영향 없이 재료의 경로와 강도를 모두 정확하게 예측할 수 있음을 입증했습니다.

마지막으로 가장 복잡한 테스트는 내부 폭발에 노출된 두꺼운 실린더의 파편화였습니다. 이 시나리오에서는 노치에 의해 유도되는 단일 균열이 존재하는 것이 아니라, 많은 균열이 자발적으로 핵을 형성하여 바깥쪽으로 질주하며 실린더를 수십 조각으로 산산조각 냅니다. 이는 격자의 대칭성이 파편의 패턴을 결정할 수 있는 혼돈스러운 사건입니다. 브라운이 실린더의 축에 맞춰 격자를 중앙에 배치했을 때, 균열은 대칭적인 패턴을 형성했습니다. 그러나 격자를 약간 이동시켜 더 이상 완벽하게 중앙에 위치하지 않게 했을 때, 균열은 불규칙하게 분포되었고 파편의 개수도 변했습니다. 이는 모델이 단순히 격자의 대칭성에 반응하는 것이 아니라 물리적 설정에 민감하게 반응함을 확인시켜 주었습니다. 시뮬레이션은 12개에서 19개의 파편을 예측했는데, 이는 다른 고급 방법들의 범위 안에 들어가는 수치이지만 일부보다 낮은 수치였는데, 이는 아마도 모델의 부드러운 손상 영역이 너무 가까이 있는 균열들의 형성을 방지했기 때문일 것입니다.

이 모든 테스트 과정에서 연구진은 시스템의 에너지를 면밀히 추적했습니다. 그들은 재료가 파괴되면서 소실된 에너지와 새로운 균열 표면을 만드는 데 필요한 에너지를 비교했습니다. 어떤 경우에는 손상 모델에 의해 계산된 에너지와 끊어진 결합에 의해 실제로 방출된 에너지 사이에 작은 불일치가 발생했습니다. 이는 하이브리드 접근 방식에서 나타나는 알려진 한계로, 균열을 일으키는 구동력과 끊어진 결합의 에너지가 서로 다르게 계산되기 때문입니다. 그러나 기계적 에너지 균형은 0.5% 이내의 오차로 정확하게 유지되었으며, 이는 시뮬레이션이 물리적으로 타당함을 확인해 주었습니다. 결과는 이 모델이 완벽하지는 않을지라도, 역동적 파괴를 이해하는 데 있어 중요한 진전임을 시사합니다.

이 연구의 함의는 단순히 더 나은 컴퓨터 그래픽을 만드는 것에 그치지 않습니다. 격자 방향에 대한 의존성을 제거함으로써, 엔지니어들은 충돌 테스트, 구조적 실패 또는 재료 테스트에 대한 시뮬레이션이 재료의 실제 물리 법칙을 반영하고 있다고 신뢰할 수 있습니다. 이 모델은 원자의 이산적인 특성과 고체 물체의 연속적인 거동 사이의 간극을 성공적으로 메웠습니다. 이는 적절한 수학적 도구를 사용한다면, 컴퓨터의 격자가 결과를 좌우하게 하지 않고도 파괴라는 혼돈스럽고 격렬한 과정을 시뮬레이션하는 것이 가능하다는 것을 보여줍니다. 이 연구는 이 접근 방식이 극한 조건에서 재료가 어떻게 파괴되는지를 예측하는 견고한 방법을 제공하며, 더 안전하고 탄력적인 구조물을 설계할 수 있는 명확한 길을 제시한다고 결론짓습니다.