A beam--membrane biomechanical vocal fold model incorporating posturing and glottal conformation

이 논문은 고가의 고충실도 시뮬레이션을 대체할 수 있는 실용적인 대안으로서, 음성 생성 역학을 예측하고 음성 장애를 조사하기 위해 근육 주도 포스처링(posturing)과 성문 형태 변화를 통합한 계산 효율적인 빔-막 구조의 성대 생체역학 모델을 소개한다.

핵심

당신의 목소리가 복잡한 악기와 같다고 상상해 보세요. 다만 현이나 리드 대신, 목 안쪽에 있는 **성대(vocal folds)**라고 불리는 두 개의 살점이 있는 형태입니다. 당신이 말을 할 때, 공기가 이 두 플랩 사이의 틈을 통과하며 지나가면서 이들을 진동시켜 소리를 만들어냅니다.

이 논문은 성대가 어떻게 움직이고 진동하는지를 시뮬레이션하는 새롭고 영리한 컴퓨터 모델을 소개합니다. 저자들은 특정한 문제를 해결하고자 했습니다. 기존의 컴퓨터 모델들은 너무 단순하거나(만화 그림처럼), 혹은 너무 복잡하여(며칠이 걸리는 슈퍼컴퓨터 시뮬레이션처럼) 어느 쪽도 완벽하지 않았습니다. 그들의 목표는 "골디락스(Goldilocks)" 모델, 즉 빠르게 실행될 수 있을 만큼 빠르면서도 과학적으로 정확할 만큼 상세한 모델을 구축하는 것이었습니다.

다음은 이 연구 내용을 쉬운 비유를 사용하여 정리한 내용입니다.

1. 문제점: "너무 느린" 모델 vs "너무 단순한" 모델의 딜레마

목소리를 연구하는 것을 자동차 엔진이 어떻게 작동하는지 이해하려는 과정이라고 생각해보세요.

• "너무 단순한" 모델들은 어린아이의 장난감 자동차와 같습니다. 밀고 다니기는 쉽지만, 피스톤이나 연료가 어떻게 작동하는지는 알려주지 않습니다.
• "너무 복잡한" 모델들은 다이노미터 위에 놓인 실제 크기의 진짜 엔진과 같습니다. 믿을 수 없을 정도로 정확하지만, 시뮬레이션을 한 번 돌리려면 거대한 슈퍼컴퓨터가 필요하며, 소리의 아주 짧은 순간을 시뮬레이션하는 데에도 몇 주가 걸릴 수 있습니다. 이 때문에 수백 가지의 다른 시나리오(예: "이 근육을 더 조이면 어떻게 될까?")를 테스트하기가 어렵습니다.

저자들은 모델이 **고품질의 원격 제어 자동차(RC Car)**처럼 작동하기를 원했습니다. 즉, 현실적으로 움직이고 컨트롤에 반응하면서도, 한 오후 동안 수천 번은 테스트할 수 있을 만큼 빨라야 한다는 것입니다.

2. 해결책: "빔과 막(Beam and Membrane)" 샌드위치

이 모델을 만들기 위해, 저자들은 성대를 함께 작동하는 두 개의 뚜렷한 부분으로 이루어진 샌드위치처럼 취급했습니다.

• 빔 (척추/중심축): 그들은 깊은 층(근육과 인대)을 단단하고 휘어지는 '빔'으로 모델링했습니다. 이것은 유연한 자(ruler)와 같습니다. 자의 양 끝을 누르면 휘어지는 것처럼 말이죠. 이 부분은 모델에서 "포스처링(posturing, 자세 잡기)"을 담당합니다. 즉, 근육이 어떻게 성대를 늘리고 위치시키는지를 처리합니다.
• 막 (피부): 그들은 맨 위의 말랑말랑한 층(점막)을 얇고 신축성 있는 '막'으로 모델링했습니다. 이것은 풍선 껍질이나 드럼 헤드와 같습니다. 이 부분은 공기가 위로 흐를 때 물결치듯 출렁거립니다.

이 두 부분은 "스프링과 댐퍼(자동차의 쇼크 업소버 같은 것)"로 결합되어 있습니다. 이를 통해 단단한 빔은 휘어지면서도 부드러운 피부는 물결칠 수 있어, '점막파(mucosal wave)'라고 불리는 현실적인 파동 운동을 만들어냅니다.

3. "근육 원격 제로 컨트롤"

이 모델의 가장 멋진 기능 중 하나는 근육을 다루는 방식입니다. 현실 세계에서는 뇌가 목 안의 미세한 근육들에게 수축하라는 명령을 내려 성대의 모양을 변화시킵니다.

• 저자들은 **"포스처링 모델(Posturing Model)"**을 만들었는데, 이는 마치 원격 제어기와 같습니다.
• 버튼을 누르면(근육을 활성화하면), 모델은 "자(beam)"가 어떻게 휘고 늘어나는지 계산합니다.
• 이렇게 휘어진 모양은 깔때기 모양(앞은 좁고 뒤는 넓은 형태)이나 활 모양(미소 짓듯 굽은 형태)과 같은 특정 형태를 만들어냅니다.
• 그런 다음 모델은 이 형태들을 가져와서 "소리" 시뮬레이션을 실행합니다.

4. 발견한 점 (결과)

저자들은 자신들의 모델이 실제 인간의 목소리 생성을 얼마나 잘 모사하는지 확인하기 위해 모델을 실행했습니다. 그들은 자신들의 "원격 제어 자동차" 결과를 실제 실험 결과 및 "슈퍼컴퓨터" 모델들과 비교했습니다.

• 성공적임: 이 모델은 복잡한 음성 행동을 성공적으로 재현했습니다. 예를 들어, 모델에게 특정 근육을 활성화하라고 "명령"했을 때, 모델은 의사들이 실제 환자들에게서 보는 특이한 모양들(예: 모래시계형 틈이나 활 모양의 굽음)을 자연스럽게 만들어냈습니다.
• "하부 가장자리 선행(Inferior Edge Lead)": 실제 상황에서는 성대의 아래쪽 가장자리가 진동하는 동안 위쪽 가장단보다 약간 앞서서 움직이는 경우가 많습니다. 이전의 단순한 모델들은 이를 인위적으로 설정해주어야 했습니다. 하지만 이 새로운 모델에서는 빔과 막이 연결된 방식 덕분에 이 현상이 자연스럽게 일어납니다. 이는 마치 실제 깃발이 펄럭이는 것과 같습니다. 바람 때문에 아래쪽 부분이 먼저 펄럭이도록 따로 프로그래밍할 필요 없이, 물리 법칙에 의해 자연스럽게 발생하는 것과 같습니다.
• 속도: 가장 큰 승리는 속도입니다. 고정밀 모델이 소리의 아주 짧은 순간을 시뮬레이션하는 데 **1,200시간(50일!)**이 걸릴 수 있는 반면, 이 새로운 모델은 표준 노트북을 사용하여 1분 미만에 동일한 작업을 수행할 수 있습니다.

5. 왜 중요한가 (논문에 따르면)

이 논문은 이 도구가 음성 장애를 이해하는 데 있어 획기적인 돌파구라고 주장합니다.

• 모델이 매우 빠르기 때문에, 연구자들은 이제 수천 번의 "만약 ~한다면(what-if)" 시나리오를 실행할 수 있습니다. 그들은 서로 다른 근육 활성화 패턴이 어떻게 비효율적인 목소리나 조직 손상(예: 성대를 너무 세게 부딪히는 경우)으로 이어지는지 테스트할 수 있습니다.
• 이는 특정 목소리 문제가 왜 발생하는지를 설명하는 데 도움을 줍니다. 예를 들어, 성대의 뒷부분이 열려 있는 상태(후방 틈)가 되면 성대의 충돌 방식이 어떻게 변하고, 이것이 잠재적으로 어떻게 부상으로 이어지는지를 보여주었습니다.

요 요약

요약하자면, 저자들은 성대의 빠르고, 스마트하며, 물리적으로 현실적인 컴퓨터 시뮬레이션을 구축했습니다. 그들은 성대를 휘어지는 빔과 그 위를 덮고 있는 물결치는 피부로 취급하였고, 이를 가상의 근육으로 제어하도록 만들었습니다. 이 모델은 슈퍼컴퓨터 없이도 목소리 생성의 복잡한 춤을 포착해내며, 우리의 목소리가 어떻게 작동하고 왜 때때로 문제가 생기는지를 연구할 수 있는 효율적인 새로운 방법을 제시합니다.

기술 요약

기술 요약: 포스처링(Posturing) 및 성문 형태(Glottal Conformation)를 통합한 빔-막 구조의 생체역학적 성대 모델

문제 제기
성대(VFs)의 포스처(posture)는 후두 근육의 활성화에 의해 결정되며, 이는 음성 생성 역학의 주요 결정 요인이다. 비정상적인 성대 구성은 빈번하게 비효율적인 발성과 음성 장애와 연관된다. 임상적 관찰을 통해 다양한 성문 폐쇄 패턴(예: 후방 개방, 보우형(bowed) 형태, 모래시계 형상)이 확인되었으나, 이를 지배하는 생체역학적 메커니즘은 완전히 이해되지 않은 상태로 남아 있다. 내재 근육의 효과를 포함하는 기존의 고충실도 유한요소 모델은 계산 비용이 매우 높아 대규모 파라미터 조사에 활용하기에 제한적이다. 반대로, 축약 차수 모델(reduced-order models)은 근육 활성화를 기계적 특성에 연결할 때 경험적인 규칙에 의존하거나 단순화된 기하학적 구조(예: 직사각형 또는 삼각형)를 가정함으로써, 생리학적으로 현실적인 성대 모양과 복잡한 폐쇄 패턴을 예측하는 데 실패한다. 따라서 생체역학적 해석 가능성을 유지하면서 성문 형태가 발성 역학에 미치는 영향을 포착할 수 있는 계산 효율적인 프레임워크가 필요하다.

방법론
저자들은 성대의 체부(body)와 피복층(cover layers)을 각각 복합 빔(composite beam)과 결합된 막(coupled membrane)으로 취급하는 계산 비용이 저렴한 성대 모델을 제안한다. 이 프레임워크는 두 가지 주요 구성 요소를 통합한다:

1. 포스처링 모델(Posturing Model): 축약된 근육 제어 프레임워크를 기반으로 하는 이 구성 요소는 내재 후두 근육(갑상피열근[TA], 윤상갑상근[CT], 측윤상피열근[LCA], 피열간근[IA], 후윤상피열근[PCA])의 정규화된 활성화 수준을 사전 발성 구성 파라미터로 매핑한다. 구체적으로, 이 모델은 윤상갑상 관절과 피열 연골의 회전 및 병진 운동을 시뮬레이션하여 명목 성대 변형률($\bar{\varepsilon}$)과 성문 반각($\theta_g$)을 계산한다. 이 파라미터들은 내부 굽힘 모멘트를 도입하여 성문 형태에 영향을 주는 입력값 역할을 한다.
2. 빔-막 동역학 모델(Beam–Membrane Dynamic Model): 각 성대는 세 가지 해부학적 층(점막(막), 성대 인대, TA 근육(빔))으로 구성된 직육면체로 모델링된다.
   • 빔 구성 요소는 1차원 오일러-베르누이(Euler-Bernoulli) 공식을 사용하여 굽힘 모멘트를 전달할 수 있는 성대 인대와 TA 근육을 나타낸다.
   • 막 구성 요소는 점막층을 2차원 표면으로 나타낸다.
   • 이 구성 요소들은 점탄성 상호작용을 모델링하기 위해 분포된 스프링-댐퍼 요소에 의해 기계적으로 결합된다.
   • 시스템은 공기역학적 하중(점성 보정이 포함된 이상 베르누이 흐름으로 모델링됨)과 충돌 압력(페널티 기반 접촉 공식으로 모델링됨)의 영향을 받는다.
   • 지배 방정식은 Matlab을 이용한 유한차분 이산화 방식을 사용하여 해결된다. 모델의 출력값에는 성문 면적 파형, 유량, 방사 음압, 그리고 기본 주파수($f_0$), 음압 레벨(SPL), 폐쇄 지수(CQ), 충돌 압력과 같은 생체역학적 지표가 포함된다.

주요 기여

• 계산 효율성: 제안된 프레임워크는 고충실도 유체-구조 상호작용(FSI) 모델에 비해 계산 비용을 대폭 절감한다. 일반적인 1초 시뮬레이션은 표준 노트북에서 1분 미만이 소요되는 반면, 유사한 고충실도 시뮬레이션은 수천 프로세서 시간을 요구할 수 있다.
• 생체역학적 해석 가능성: 경험적인 기하학적 제약을 부과하는 럼프드 매스(lumped-mass) 모델과 달리, 이 연속체 기반 모델은 근본적인 기계적 원리와 조직 재료 특성으로부터 강성 및 관성 항을 도출한다. 이 모델은 근육 활성화가 성대 모양을 변화시키는 핵심 메커니즘인 굽힘 모멘트의 전달을 자연스럽게 포착한다.
• 동적 성문 형태: 모델은 경험적인 기하학적 규칙 없이도 근육 활성화 패턴으로부터 복잡하고 임상적으로 관찰되는 정적 및 동적 성문 형태(예: 전후방 개방, 내측 돌출, 수렴/발산 프로파일)를 성공적으로 재현한다.
• 검증: 본 프레임워크는 고충실도 계산 연구 및 실험적 관찰 결과와 비교 검증되었으며, 정적 구성 및 발성 측정치에서 질적인 일관성을 입증하였다.

결과
수치 시뮬레이션은 다양한 근육 활성화 시나리오에 대한 모델의 예측 능력을 보여준다:

• 정적 구성: 모델은 PCA 활성화에 의한 후방 성문 개방, 낮은 TA 활성화에 의한 전방 개방, 증가된 TA 활성화에 의한 내측 돌출(오목한 형태) 등 알려진 임상 패턴을 재현한다. 또한 CT 활성화에 의한 직선화 효과를 포착한다.
• 모달 발성(Modal Phonation): 지속적인 모달 발성 사례에서 모델은 비대칭 파형을 가진 주기적인 성문 개방 및 폐쇄를 생성한다. 모델은 속도 지수(SQ) 1.52와 폐쇄 지수(CQ) 0.49를 나타내며, 이는 건강한 모달 발성과 일치한다. 시뮬레이션은 교대하는 수렴 및 발산 성문 형태와 유체-구조 상호작용에 의해 구동되는 하단에서 상단 가장자리로의 접촉 압력 전파를 보여준다.
• 파라미터 연구 (CT 및 TA 활성화): 활성화 맵은 $f_0$가 주로 CT 활성화(증가)와 TA 활성화(감소)에 의해 지배됨을 보여주며, 이는 고충실도 모델의 경향과 질적으로 일치한다. 모델은 근육 활성화에 따른 SPL 및 유량의 비선형적 의존성을 포착한다.
• 후방 성문 개방 (PCA 활성화): PCA 활성도가 높아지면 후방 성문 간격이 유도된다. 모델은 이것이 평균 성문 기하 구조를 변화시킬 뿐만 아니라, 공기역학적 및 접촉 압력 분포에 상당한 전후방 비대칭을 유도함을 밝혀낸다. 높은 PCA 활성화 수준에서 모델은 고차 진동 모드의 출현과 최대 충돌 압력의 비단조적 변화를 예측하며, 이는 기하 구조와 충돌 역학 사이의 비선형적 상호작용을 강조한다.

의의 및 주장
본 논문은 제안된 프레임워크가 발성 생체역학을 조사하기 위한 실용적이고 계산 가능한 도구를 제공한다고 주장한다. 고충실도 유한요소 모델과 단순화된 축약 차수 모델 사이의 간극을 메움으로써, 효율성과 생리학적 현실성 사이의 균형을 맞춘다. 저자들은 모델이 고충실도 유한요소 모델 및 임상 연구에서 보고된 질적 경향을 재현함으로써 그 예측 능력을 뒷받침한다고 기술한다.

본 연구의 의의는 비정상적인 근육 활성화 및 비효효율적인 성대 기능(예: 근긴장성 발성장애, 성대 과기능)과 관련된 음성 장애의 메커니즘을 조사하기 위한 대규모 파라미터 조사를 용이하게 할 수 있다는 점에 있다. 이 프레임워크는 조직 외상과 관련된 공간 정보를 보존하면서 근육 활성화 패턴, 조직 특성 및 성문 구성이 발성 결과에 어떻게 영향을 미치는지에 대한 체계적인 검토를 가능하게 한다. 저자들은 성도와의 양방향 음향 결합 부재 및 일부 다른 연구에서 관찰된 낮은 TA 활성화 단계에서의 초기 $f_0$ 증가 현상 부재와 같은 한계점을 언급하며, 이는 향후 개선이 필요한 영역임을 시사한다. 그러나 이 모델이 발성을 지배하는 필수적인 생체역학적 메커니즘을 포착하고 있다는 점을 강조한다.