유한요소해석

유한요소법(FEM)의 개념적 기반은 1940년대 초로 거슬러 올라간다. 1941년 브리티시컬럼비아 대학의 러시아계 캐나다인 구조 공학자 알렉산더 흐렌니코프(Alexander Hrennikoff)는 연속 막·판 구조를 이산 격자 프레임워크로 근사하는 논문을 발표하였다. 1943년 수학자 리하르트 쿠란트(Richard Courant)는 각주 막대의 비틀림 문제를 풀기 위해 삼각형 및 직사각형 하위 영역(사실상 초기의 '유한요소')을 도입하고 변분법을 적용하여 편미분방정식의 근사해를 구하였다. 그러나 FEM을 실용적인 계산 도구로 발전시킨 것은 공학 분야였다. 1956년 보잉사의 M.J. 터너(Turner), R.W. 클러프(Clough), H.C. 마틴(Martin), L.J. 톱(Topp)은 이산 삼각 요소를 이용한 직접 강성법(direct stiffness method)을 제시하는 획기적인 논문을 발표하였다. 레이 W. 클러프는 1960년 ASCE 전자 계산 학회에서 발표한 'The Finite Element Method in Plane Stress Analysis'라는 논문에서 공식적으로 '유한요소법(finite element method)'이라는 용어를 만들었다. 1960~70년대에 걸쳐 스완지 대학의 올기에르드 지엔키에비치(Olgierd Zienkiewicz)와 슈투트가르트 대학의 존 H. 아지리스(John H. Argyris)가 수학적·전산적 프레임워크를 발전시켜 FEM을 토목, 기계, 항공우주 공학 전반에 광범위하게 적용할 수 있게 하였다. 1980~90년대에는 컴퓨팅 성능의 향상과 ANSYS, Abaqus, Strand7 등 상용 FEA 소프트웨어의 개발로 생물학자와 고생물학자를 포함한 훨씬 넓은 연구 커뮤니티에서 이 방법을 활용하게 되었다.

FEA가 공학에서 생물학으로 전이된 과정은 정형외과 및 치과 생체역학 분야에서 시작되었으며, 이후 기능적 하중 조건에서의 뼈, 임플란트, 치아의 역학적 거동 연구에 FE 모델이 사용되었다. 1990년대 후반이 되자 이 방법은 척추동물 기능형태학 연구 문제에 적용되기 시작하였다. FEA의 고생물학 적용에서 획기적인 연구는 2001년 에밀리 J. 레이필드(Emily J. Rayfield)와 동료들이 학술지 Nature에 발표한 것으로, 쥬라기의 대형 수각류 공룡 알로사우루스 프라길리스(Allosaurus fragilis)의 두개골을 분석하였다. 이 연구는 당시 생물체에 대해 구축된 것 중 가장 기하학적으로 복잡한 FE 모델을 생성하여, 알로사우루스의 두개골이 근육 구동 교합력의 추정치를 훨씬 초과하는 힘을 견딜 수 있음을 보여줌으로써 두개골이 구조적으로 '과설계(overbuilt)'되어 있었음을 시사하였다. 2004년 레이필드는 티라노사우루스 렉스(Tyrannosaurus rex) 두개골에 대한 분석을 후속 발표하여 교합 및 인열(tearing) 하중을 시뮬레이션함으로써 '천공-당김(puncture-pull)' 섭식 가설을 검증하였다. 결과는 티라노사우루스의 두개골이 두 종류의 하중 모두에 동등하게 적응되어 있음을 보여주었고, 압축 및 전단 응력이 (알로사우루스처럼 전두-두정부가 아닌) 비골에 집중됨을 밝혀 티라노사우루스과 비골의 견고함에 대한 생체역학적 근거를 제공하였다. 또한 상악-관골 봉합선이 인장 충격 흡수 기능을 수행하여 국소 인장 응력을 줄이지만 전체적으로는 두개골을 약간 약화시키는 것으로 나타났다. 레이필드의 2007년 Annual Review of Earth and Planetary Sciences 리뷰 논문은 고생물학에서의 FEA에 대한 포괄적 프레임워크를 제시하여 이후 연구들의 기반이 되는 참고문헌이 되었다.

전형적인 고생물학적 FEA 작업 흐름은 여러 단계로 이루어진다. 먼저 화석 구조의 3차원 기하학 정보를 CT 스캔(마이크로 CT 또는 의료용 CT)으로 획득하며, 사진측량법이나 표면 레이저 스캐닝도 사용된다. 획득한 체적 데이터를 Avizo, Mimics, 또는 3D Slicer 같은 오픈소스 대안을 사용하여 세그멘테이션하여 주변 매트릭스로부터 관심 골격 구조를 분리한다. 세그멘테이션된 표면 메시는 구멍, 비다양체(non-manifold) 모서리, 겹치는 면, 지나치게 불균일한 요소 크기 등의 결함을 제거하기 위해 정리 및 보수 과정을 거친다. 정리된 표면 메시는 일반적으로 4절점 사면체 요소(tet4)로 구성된 체적 메시로 변환되며, 특정 응용에서는 고차 요소나 셸, 빔, 판 등 다른 기하학적 형태의 요소도 사용된다.

메시 요소에는 재료 물성값이 할당되어야 한다. 대부분의 고생물학 연구에서 뼈는 균질, 등방성, 선형 탄성 재료로 모델링되며, 영률(E)은 약 17,000~20,000 MPa, 푸아송비(ν)는 약 0.3으로 설정하는 것이 일반적이다. 이 값들은 현생 주룡류(archosaur) 뼈에 대한 실험 데이터에서 도출된다. 일부 연구에서는 비균질 또는 이방성 물성을 할당하며, 최근에는 치아 에나멜(E ≈ 50,000~80,000 MPa)과 상아질(E ≈ 18,000 MPa)을 구분하여 모델링하는 효과를 탐구하는 연구도 이루어지고 있다.

이어서 경계 조건을 정의한다. 근력은 모델 표면의 근육 부착 영역에 적용되며(건조 두개골법, 근육 재구성, 생리학적 단면적 계산 등으로 추정), 제약 조건(고정 변위점)은 측두하악관절이나 치아-먹이 접촉점 등 생물학적으로 적절한 위치에 배치된다. 하중 시나리오는 다양한 치아 위치에서의 편측 또는 양측 교합 같은 행동을 시뮬레이션한다. 모델은 상용 소프트웨어(ANSYS, Abaqus, Strand7) 또는 점차 증가하는 오픈소스 도구(Fossils, FEBio)의 선형 정적 솔버를 사용하여 풀린다. 출력은 폰 미제스 응력, 주응력 및 변형률, 변형, 변형 에너지, 제약 노드에서의 반력(교합력) 분포를 포함한다.

고생물학적 FE 모델에서 일반적으로 추출하는 여러 역학 지표가 있다. 폰 미제스 응력은 모든 응력 성분의 효과를 결합한 스칼라 값으로, 뼈와 같은 연성 재료의 항복을 예측하는 데 유용하여 가장 널리 보고된다. 주응력(최대 및 최소)은 인장 대 압축 하중의 방향과 크기를 보여준다. 변형 에너지는 구조를 변형시키는 데 투입된 총 일로, 입력 근력 중 출력 교합력으로 전환되지 않고 변형에 소모되는 에너지가 적을수록 구조적으로 효율적인 것으로 평가된다. 역학적 효율(mechanical efficiency)은 교합 반력 대 총 입력 근력의 비율로 정의되며, 분류군 간 섭식 성능을 직접 비교할 수 있게 한다.

고생물학적 FEA에서 핵심적인 도전은 모델 검증이다. 멸종 생물의 역학적 거동을 실험적으로 검증하는 것은 불가능하기 때문이다. 현생 유사 생물을 사용한 검증 연구가 수행되어 왔다. 스트레이트 등(Strait et al., 2005)은 교합 중 마카크(Macaca) 두개골에서의 생체 내(in vivo) 변형률 게이지 데이터와 FE 예측치를 비교하여, 절대적인 변형률 크기는 때때로 부정확했지만 근력과 재료 물성이 정확하게 반영되었을 때 변형률 분포의 일반적인 패턴과 방향은 모델에 의해 잘 예측됨을 발견하였다. 미시시피악어(Alligator mississippiensis) 하악골(Porro et al., 2011), 돼지 두개골(Bright and Gröning, 2011), 도마뱀 두개골(Dutel et al., 2021) 등에서도 유사한 검증 연구가 수행되었다. 이러한 연구들은 응력과 변형률의 절대 크기는 신중하게 해석해야 하지만, 동일한 프로토콜로 구축된 모델 간의 응력 분포 패턴에 대한 상대적 비교는 견고하고 과학적으로 의미 있음을 종합적으로 입증하였다.

고생물학적 FE 모델은 수많은 단순화와 가정 위에 구축되며, 그 영향을 검토하는 상당한 양의 민감도 분석 문헌이 축적되어 있다. 모델 출력에 영향을 미치는 주요 매개변수는 메시 밀도(조잡한 메시는 부정확한 응력 추정을 초래할 수 있음), 재료 물성 할당(에나멜, 상아질, 봉합선 포함 여부가 국소 응력 패턴과 교합력 추정에 영향), 근력 크기 및 방향(근육 감싸기 vs. 직선 힘 벡터의 효과 포함), 경계 조건 배치(제약 노드의 수와 위치), 두개골 봉합선, 연골두개, 치주인대, 내부 해부학적 구조 등의 포함 또는 배제이다. 허브스트 등(Herbst et al., 2021)은 치아를 모두 뼈로 단순화한 모델과 실제 에나멜 층을 모델링한 모델 사이에서 전체 응력 분포 패턴에는 눈에 띄는 차이가 없었지만, 절대 교합력 값은 최대 약 28%까지 영향을 받았음을 보여주었다. 레이필드(Rayfield, 2019)는 두개골 봉합선의 부적절한 모델링이나 생략이 응력, 변형률, 변형에서 부정확한 결과를 초래할 수 있음을 지적하였다.

고생물학적 FEA 연구의 대다수는 작은 변형, 선형 탄성 재료(훅의 법칙), 정적 평형을 가정하는 선형 정적 해석을 사용한다. 이 접근법은 계산 비용이 적고 골격 요소에 관한 많은 생체역학적 질문에 충분하다. 그러나 마르세-노게(Marcé-Nogué, 2022)가 리뷰한 바와 같이, 비선형 FEA는 고생물학에서 충분히 활용되지 못하고 있는 중요한 역량을 제공한다. 재료 비선형성은 초탄성(hyperelastic) 연조직(인대, 연골, 치주인대)과 뼈의 소성 변형을 모델링할 수 있게 한다. 기하학적 비선형성(대변위)은 대형 공룡의 긴 사지골이나 익룡 날개뼈 같은 가늘고 긴 구조의 좌굴(buckling) 해석을 가능케 한다. 접촉 비선형성은 손목뼈, 청소골, 척추 등에서 골면 간의 분리와 미끄러짐이 가능한 관절 다체 시스템의 시뮬레이션을 허용한다. 이러한 비선형 접근법은 뉴턴-랩슨(Newton-Raphson)법 등의 반복 솔버와 상당히 큰 계산 자원을 필요로 하지만, 속성학적으로 변형된 화석의 역변형(retrodeformation)이나 극한 하중 하에서의 골격 요소 구조 파괴 등 선형 모델로는 다룰 수 없는 문제를 해결할 수 있다.

FEA는 비교 기능형태학의 필수 도구로 자리 잡아, 골격 형태의 진화적 최적화에 관한 가설을 검증할 수 있게 해 준다. 하나의 분기군 내 여러 분류군에 대해 FE 모델을 구축하고 표준화된 하중 조건을 적용함으로써, 생체역학적 성능 차이를 정량화하고 이를 식이, 섭식 전략, 이동 양식 등의 생태적 변수와 연관시킬 수 있다. 수각류 공룡(티라노사우루스과의 극단적 골식증(osteophagy)으로 이어지는 대진화적 경향 포함), 식육목 포유류, 악어형류, 분추류 양서류 등 다양한 그룹에서 두개골 생체역학 연구가 수행되었다. 2005년에서 2020년 사이에 750건 이상의 논문이 FEA를 생물학적·진화적 질문에 적용하였다(Tseng, 2021). 그러나 첸(Tseng, 2021)은 중요한 방법론적 주의점을 입증하였다: 3D FE 연구에서 분류군 표본 크기가 작은 경우(약 11~12개 분류군 미만) 생체역학적 성능과 섭식 생태 간에 인위적으로 부풀려진 상관관계가 나타나 위양성(false positive) 결과를 생성할 수 있다는 것이다. 이 발견은 대표적인 분류군 표본 추출과 계통 구조에 대한 고려를 포함한 신중한 연구 설계, 또는 대안적으로 이론적 형태공간(theoretical morphospace) 접근법의 채택을 권장하게 되었다.

고생물학적 FEA는 전통적으로 Strand7, Abaqus, ANSYS, Hypermesh 등의 상용 공학 소프트웨어에 의존해 왔다. 최근에는 생물학적 응용에 특화된 오픈소스 대안들이 개발되고 있다. Fossils 소프트웨어(Chatar et al., 2023)는 뼈에 대한 근육 구동 생체역학적 하중을 시뮬레이션하는 무료의 효과적인 도구를 제공하며, 상용 패키지와 동등하거나 그 이상의 성능을 보인다. BFEX 애드온(Díaz de León Muñoz et al., 2025)은 Fossils를 오픈소스 3D 소프트웨어인 블렌더(Blender)와 통합하여 모델 준비 작업 흐름을 간소화한다. FEBio(Maas et al., 2012)는 생물학적 유한요소해석을 위한 전문 기능을 제공한다. 이러한 발전들은 FEA에 대한 접근을 민주화하여, 고가의 소프트웨어 라이선스 없이도 엄밀한 생체역학적 분석을 수행할 수 있게 해주고 있다.

최근의 방법론적 진보로는 자동화된 CT 세그멘테이션을 위한 딥러닝의 통합(Zhang et al., 2025), 정적 스냅샷이 아닌 동적 섭식 시퀀스를 시뮬레이션하기 위한 FEA와 다체동역학 해석(MDA)의 결합, FEA 결과를 기하학적 형태계측학 및 계통비교방법론과 연결하여 전체 분기군에 걸친 형태-기능 진화를 연구하는 방법이 있다. 이 분야는 비선형 재료 모델, 다체 접촉 시뮬레이션, 현생 근연종의 계통학적 브라케팅(phylogenetic bracketing)을 통한 연조직 물성의 더 나은 특성화를 통해 모델의 현실성을 개선하려는 지속적인 노력과 함께 빠르게 성장하고 있다.