동작 캡처만으로 지면 반력 추정하기

초록

본 논문은 마커 기반 3차원 동작 캡처 데이터만을 이용해 전신의 무게중심(CoM)을 계산하고, 이를 기반으로 뉴턴의 제2법칙을 적용해 지면 반력(GRF)을 추정·분해하는 방법을 제시한다. 16개 신체 구간 모델링, 최소 토크 변화 원칙을 활용한 이중지지 단계 해석, 그리고 실험적으로 수집한 10 m 보행 데이터와 기존 상용 소프트웨어(Visual3D)와의 비교를 통해 CoM 추정 정확도와 GRF 재구성의 타당성을 검증하였다.

상세 분석

이 연구는 임상 현장에서 비용·공간·전문 인력이 제한적인 상황에서도 보행 역학을 정량화할 수 있는 ‘힘판‑프리(force‑plate‑free)’ 프레임워크를 구축한다는 점에서 의미가 크다. 먼저 저자들은 Dumas 등(2011)의 16구간 전신 모델을 채택하고, 각 구간의 질량 비율과 무게중심 위치를 최신 인체인체계계수(anthropometric tables)로 보정하였다. 손 구간은 마커가 부족한 경우 Winter et al.(1998)의 비례법을 적용해 추정함으로써, 실제 실험에서 흔히 발생하는 상지 마커 결손 문제를 해결한다.

CoM은 각 구간의 무게중심 좌표를 질량 가중 평균하여 실시간(프레임 단위)으로 계산되며, 이때 사용된 좌표계는 국제생체역학학회(ISB) 표준을 변형해 전후‑측‑수직 축을 정의한다. 이렇게 얻어진 CoM 위치와 2차 미분을 통해 가속도를 구하고, 중력 가속도를 보정한 뒤 뉴턴의 제2법칙 ( \sum \mathbf{R}=m\mathbf{a} ) 을 적용한다. 단일 지지 단계에서는 한쪽 다리의 GRF가 0이므로 직접적으로 구할 수 있지만, 이중 지지 단계에서는 두 다리의 GRF가 동시에 작용하므로 방정식이 하나 부족하다. 이를 보완하기 위해 저자들은 ‘최소 토크 변화(minimum torque‑change)’ 원칙을 최소화 문제로 정식화하고, 라그랑주 승수를 이용해 닫힌 형태의 해를 도출하였다. 이 접근법은 기존 연구에서 제시된 역학적 최적화와 유사하지만, 구현이 간단하고 실시간 적용이 가능하다는 장점이 있다.

실험에서는 4명의 건강한 피험자를 대상으로 10 m 보행 테스트를 수행했으며, 마커 데이터는 200 Hz, 힘판 데이터는 2000 Hz로 수집하였다. CoM 추정 결과는 Visual3D와 비교했을 때 전후축(X)에서 3.9–6.2 cm, 측축(Y)에서 0.31–0.52 cm, 수직축(Z)에서 4.3–8.9 cm의 평균 오차를 보였다. 특히 Z축에서 일정한 상향 편향이 있었는데, 이는 사용된 인체 모델(Dumas)과 Visual3D가 적용한 Dempster·Hanavan 모델 간의 질량 분포 차이에서 기인한다. 편향을 7.2 cm 보정하면 잔여 오차는 0.01–2.93 cm 수준으로 감소한다.

GRF 재구성 측면에서는, 추정된 CoM 가속도와 최소 토크 변화 해를 이용해 수직·전후·측방 성분을 각각 계산했으며, 보행 주기 동안의 GRF 파형이 실제 힘판 측정값과 유사한 형태를 보였다. 특히 발 착지·이륙 시점(heel‑strike, toe‑off) 검출은 sacrum 기준의 전후 변위 극값을 이용한 Zeni 방법과 일치했으며, 이는 추정된 GRF를 이용한 이벤트 검출이 충분히 신뢰할 수 있음을 시사한다.

한계점으로는(1) 피험자 수가 적고 연령·성별 다양성이 제한적이며, (2) 실험 환경이 깔끔한 직선 보행에 국한돼 있어 급격한 방향 전환·계단·불규칙 지면 등 복합 움직임에 대한 일반화가 아직 검증되지 않았다. 또한 CoM 추정 시 상지 마커 결손을 보완하기 위해 사용한 비례법은 개인별 손의 형태 차이를 반영하지 못하므로, 고정밀 역학 분석이 요구되는 경우 오차가 누적될 가능성이 있다. 마지막으로 최소 토크 변화 원칙이 실제 근육·관절 제어 전략을 완전히 대변하지는 않으며, 다른 최적화 기준(예: 최소 에너지, 최소 힘)과의 비교 연구가 필요하다.

전반적으로 이 논문은 마커 기반 동작 캡처만으로도 충분히 정확한 CoM과 GRF를 추정할 수 있음을 실증적으로 보여주며, 임상·재활 현장에서 비용 효율적인 보행 분석 도구로 활용될 가능성을 제시한다.