A Fourier series solution for the free in-plane vibration of a rectangular plate with a straight through crack

거대한 평평한 금속판, 예를 들어 비행기의 외피나 배의 바닥과 같은 것을 상상해 보십시오. 이 판이 진동할 때, 그것은 단순히 트램펄린처럼 위아래로 까닥거리는 것에 그치지 않습니다. 또한 옆으로 흔들리거나 길이 방향으로 늘어나기도 합니다. 이러한 '면내(in-plane)' 흔들림은 소리와 에너지의 숨겨진 고속도로와 같습니다. 만약 기계에서 들려오는 이상한 웅웅거림을 들어본 적이 있거나, 특정 속도에서 자동차가 덜덜거리는 것을 느낀 적이 있다면, 여러분은 아마도 이 보이지 않는 파동을 마주했던 것입니다. 엔지니어들은 이 판이 정확히 어떻게 흔들릴지 예측해야 하며, 그렇지 않으면 구조물을 파괴할 정도로 격렬하게 흔들리는 '공진' 지점에 우연히 도달하게 될 수도 있습니다.

이제 똑같은 금속판에 금(crack)이 가 있다고 상상해 보십시오. 금은 재료의 직물에 생긴 작고 들쭉날쭉한 찢어진 틈과 같습니다. 이는 판을 더 약하게 만듭니다. 그리고 파동이 통과하는 방식도 엉망이고 예측 불가능하게 만듭니다. 과학자들은 수십 년 동안 금이 간 판이 위아래로 흔들리는 방식(눈에 보이는 쉬운 종류의 진동)을 밝혀내기 위해 노력해 왔지만, 옆으로 흔들리는 진동은 훨씬 더 풀기 어려운 숙제였습니다. 그것은 마치 매듭지어진 두 밧줄을 푸는 것과 같습니다. 금이 두 방향의 움직임을 동시에 어지럽히기 때문에 수학적으로 매우 복잡해지기 때문입니다. 이러한 주파수를 계산하는 좋은 방법이 없었기에, 엔지니어들은 대략적인 추측이나 매우 느린 컴퓨터 시뮬레이션에 의존해야 했으며, 이는 설계를 위해 이상적인 방법이 아닙니다.

이 논문은 이 까다로운 숙제를 해결하기 위해 등장했습니다. 저자인 Xutao Sun은 금이 간 직사각형 판이 면내에서 어떻게 진동하는지 정확히 예측할 수 있는 영리하고 새로운 수학적 레시피인 '푸리에 급수 해법(Fourier series solution)'을 개발했습니다. 이 판을 거대한, 복잡한 드럼 가죽이라고 생각하고, 이를 관리 가능한 작은 직사각형 조각들로 나누어 본다고 상상해 보십시오. 금은 조각들이 서로 어떻게 움직여야 할지 합의하지 못하는 벽 역할을 합니다. 이를 해결하기 위해 저자는 '영역 분할(domain decomposition)'이라는 기법을 사용하는데, 이는 엉망이 된 피자를 네 개 또는 여섯 개의 깔끔한 조각으로 잘라 각 조각을 따로 분석하는 것과 같습니다.

진정한 마법은 저자가 각 조각의 움직임을 설명하는 방식에서 일어납니다. 단순한 파동을 사용하는 대신, 이 해법은 '이중 푸리에 코사인 급수(double Fourier cosine series)'를 사용합니다. 울퉁불퉁한 도로를 묘사하기 위해 매끄럽고 완만한 언덕들을 겹겹이 쌓아 올린다고 상상해 보십시오. 저자는 이러한 수학적 언덕들을 쌓아서 진동의 형태를 맞춥니다. 하지만 핵심은, 금이 움직임의 갑작스러운 도약이나 '불연속성'(도로 위의 계단과 같은 것)을 만들어낸다는 점입니다. 따라서 일반적인 언덕만으로는 충분하지 않습니다. 그래서 저자는 특수한 '보조 함수(auxiliary functions)'를 추가합니다. 이것들은 조각들이 만나는 지점과 금이 있는 곳의 가장자리를 부드럽게 만들어주는 수학적 경사로나 다리 역할을 합니다. 이를 통해 판이 금을 가로질러 움직이려 할 때 수학적 모델이 무너지지 않도록 보장합니다.

이 논문은 다양한 유형의 균열(가장자리에서 시작되는 측면 균열과 내부를 떠다니는 내부 균열)이 있는 판을 시뮬레이션함으로써 이 방법을 테스트합니다. 결과는 인상적입니다. 이 방법은 빠르게 수렴합니다. 즉, 저자가 쌓아 올리는 '언덕'을 더 많이 추가할수록 답이 매우 빠르게 정밀한 숫자로 안정된다는 뜻입니다. 균열이 없는 판에 대한 기존 데이터 및 균열이 있는 판에 대한 새로운 컴퓨터 시뮬레이션과 비교했을 때, 결과는 거의 완벽하게 일치했습니다.

가장 흥러운 발견 중 하나는 금이 모든 진동에 똑같이 영향을 미치지는 않는다는 점입니다. 이는 기타 줄과 같습니다. 줄의 움직임이 거의 없는 지점에 손가락을 대면 음정이 변하지 않습니다. 마찬가지로, 금이 특정 진동 모드의 '데드 존(dead zone)'에 위치하면 주파수는 거의 변하지 않습니다. 그러나 금이 판을 강하게 늘리거나 전단(shear)하는 지점에 위치하면, 주파수는 급격히 떨어질 수 있습니다. 판의 가장자리가 어떻게 고정되어 있는지에 따라 때로는 거의 40% 또는 50%까지 떨어지기도 합니다. 예를 들어, 고정된 가장자리(꽉 잡힌 상태) 근이나 자유로운 가장자리(느슨하게 흔들리는 상태) 근처에 금이 있는 판은, 모든 면이 단순히 지지된 판보다 손상에 훨씬 더 민감합니다.

또한 이 논문은 판이 만드는 실제 움직임 패턴인 '모드 형상(mode shapes)'을 그려냅니다. 금이 길어질수록 이러한 패턴은 왜곡됩니다. 온전한 판의 매끄럽고 대칭적인 파동이 깨지고, 움직임 벡터(금의 끝부분 주변에서 금속이 어느 방향으로 움직이는지를 보여주는 화살표)는 튀거나 방향을 바꿉니다. 저자는 균열의 길이와 위치가 다른 판들에 대한 방대한 양의 새로운 주파수 수치 표와 시각적 지도를 제공하며, 이전에는 상세한 데이터가 존재하지 않았던 과학적 기록의 공백을 채웠습니다.

요약하자면, 이 연구는 금이 간 금속판이 어떻게 흔들리는지 계산할 수 있는 강력하고 정확하며 빠른 방법을 제공합니다. 금이 일반적으로 진동 주파수를 낮추기는 하지만, 그 하락 폭은 금의 위치, 길이, 그리고 판이 어떻게 고정되어 있는지에 따라 크게 달라진다는 것을 확인해 줍니다. 이 새로운 수학적 도구는 엔지니어들에게 구조물의 건강 상태를 점검하고 위험한 진동을 피하도록 설계할 수 있는 신뢰할 수 있는 방법을 제공하여, 우리가 의지하는 비행기, 배, 기계들이 약간의 흉터가 생기더라도 안전하게 유지될 수 있도록 합니다.