기술 요약: 다중 접촉 이동에서의 진자 구조의 출현

문제 제기

선형 역진자 모델 (LIPM) 은 다리형 이동 제어에서 보편적인 가정으로, 종종 제어기의 목표에서 자연스럽게 도출되기보다는 모델링 제약으로 부과됩니다. 가변 높이나 각운동량에 대한 확장들이 존재하지만, 진자 형태는 일반적으로 '수동적으로' 추가됩니다. 본 논문은 더 좁은 질문을 다룹니다: 어떤 특정 목적 함수가 진자 힘 패턴을 최적 해로 자연스럽게 출현시키며, 이 그림이 어디서 무너지는가?

저자들은 각운동량 변화율 ($\|\dot{H}_G\|$) 을 관성력 및 수직 가속도에 대한 표준 항과 함께 패널티로 부과하는 작은 질량 중심 최적 제어 문제 (OCP) 를 분석합니다. 목표는 각운동량 변화율 최소화가 시스템을 자연스럽게 진자 다양체 ($\dot{H}_G = 0$) 로 유도하는지 확인하고, 다양한 지지 자세 (풀랭크 vs $N=2$) 및 작업 요구 사항 하에서 이 행동의 한계를 규명하는 것입니다.

방법론

본 논문은 다음으로 구성된 비용 함수를 최소화하는 Net-Wrench OCP(정의 1) 를 수립합니다:

1. 각운동량 패널티: $\alpha \|\dot{H}_G\|^2$
2. 작업 추적: $\lambda \|\dot{H}_G - \dot{H}_{task}\|^2$ (비균형 작업의 경우)
3. 수직 가속도 패널티: $\beta (n^\top \ddot{c})^2$
4. 힘 정규화: $\gamma \sum \|f_i\|^2$

분석은 점질량 4 족동물과 Unitree Go1 로봇에 대한 MuJoCo 시뮬레이션 (개루프 QP 및 폐루프 토크 수준 MPC 모두) 을 통한 수치 실험으로 검증된 세 가지 이론적 단계를 거칩니다:

1. 풀랭크 지지 자세 분석: 모멘트 야코비안의 특이값 분해 (SVD) 를 사용하여 각운동량 가중치 $\alpha \to \infty$일 때의 점근적 행동을 유도합니다.
2. $N=2$ 지지 자세 분석: 모멘트 연산자의 랭크가 1 인 랭크 결손 경우 (예: 대각선 트로트) 를 조사하여 불가피한 잔차와 마찰 원뿔 제약을 분석합니다.
3. 작업 주도 분석: 작업 추적 항 ($\lambda > 0$) 을 추가하는 것이 최적해를 진자 다양체에서 어떻게 벗어나게 하는지 검토합니다.

주요 기여 및 이론적 결과

1. 풀랭크 지지 자세에서의 진자 구조 출현

• 정리 1: 풀랭크 지지 자세에서 최적 해는 $\dot{H}_G \to 0$인 진자 힘 패턴을 향해 $O(1/\alpha)$ 속도로 점근적으로 수렴합니다.
• 기하학적 상수: 비례 상수는 모멘트 야코비안의 SVD 에 의해 결정되며, 이는 직접적으로 발간격 기하학에 해당합니다.
• 고전 모델과의 연결: 이 진자 다양체에서 제로 모멘트 점 (ZMP) 은 자연스럽게 가상 피벗 포인트와 같아지며, 발산 운동 성분 (DCM) 궤도 방정식은 가정이 아닌 결과로서 도출됩니다.

2. $N=2$ 지지 자세 (대각선 트로트) 의 근본적 한계

• 정리 5 (기하학적 바닥): $N=2$ 지지 자세에서 모멘트 야코비안은 랭크 1 입니다. 힘 재분배와 무관하게 두 발을 연결하는 선 ($\hat{D}$) 방향을 따라 잔류 각운동량 변화율이 남습니다. 이는 $\|\dot{H}_G\|$에 대한 0 이 아닌 하한 (바닥) 을 생성합니다.
• 정리 6 (마찰 원뿔 굴곡): 두 번째 장애물은 마찰 원뿔에서 발생합니다. 수직 모멘트 성분을 상쇄하는 데 필요한 힘이 마찰 원뿔 바깥에 있으면, 임계 수평 가속도 ($a_x^*$) 에서 "비연속적인 굴곡"이 발생합니다. 이 지점을 넘어서면 달성 가능한 $\|\dot{H}_G\|$는 기하학적 바닥을 엄격히 초과합니다.
• 함의: 대각선 트로트에서 진자 방향은 종종 마찰 원뿔 밖으로 기울어집니다 (예: $\mu=0.6$인 경우 $40.6^\circ$ 대 $31^\circ$). 이로 인해 진자 해가 불가능해지고 $\alpha$에 무관한 바닥이 생성됩니다.

3. 작업 주도 편차

• 정리 7: 작업 항 ($\lambda > 0$) 이 추가되면 최적 각운동량 변화율은 부드러운 보간이 됩니다: $\dot{H}_G^* = \frac{\lambda}{\alpha + \lambda} \dot{H}_{task}^*$.
• 상도: 시스템은 '균형 주도' 영역 (진자, $\lambda/\alpha \ll 1$) 에서 '작업 주도' 영역 (비진자, $\lambda/\alpha \gg 1$) 으로 부드럽게 전환됩니다. 이는 동일한 로봇이 정적 균형에서는 LIPM 과 유사하게 보이지만 민첩한 기동 중에는 제어기를 전환하지 않고도 비진자적으로 행동할 수 있는 이유를 설명합니다.

실험 결과

저자들은 Unitree Go1 과 점질량 모델을 사용하여 이러한 주장을 검증했습니다:

• 속도 검증 (점질량): $\alpha$가 증가함에 따라 진자 힘 패턴과의 편차가 227 배 감소했으며, 로그 - 로그 기울기는 -1 로, $O(1/\alpha)$ 속도 (정리 1) 를 확인했습니다.
• 기하학적 상수 (Go1 N=4): 네 발 지지 자세에서 경험적 스케일링 상수 ($K_e = 8.4$) 는 발간격에서 유도된 분석적 예측 ($K_a = 9.8$) 과 약 16% 이내에서 일치했습니다. 모멘트 야코비안의 특이값은 물리적 발간격과 소수점 넷째 자리까지 일치했습니다.
• 바닥 검증 (Go1 N=2): 대각선 트로트에서 $\|\dot{H}_G\|$는 예측된 기하학적 바닥인 $0.2835 \, \text{m}^2/\text{s}^2$에서 포화되었으며, 분석적 하한과 네 자리 숫자까지 일치했습니다.
• 폐루프 현실: 폐루프 토크 수준 MPC 에서 바닥은 유지되었으나 정적 예측 ($\approx 1.71 \, \text{m}^2/\text{s}^2$) 보다 높았습니다. 저자들은 이 간격이 준정적 분석에서 포착되지 않은 동적 효과 (스윙 충격, 접촉 전환, 서보 오차) 에 기인한다고 주장합니다.
• ZMP 동일성: 수치 실험은 진자 다양체에서 힘 기반 ZMP 가 가상 피벗으로 수렴함을 확인하여 LIPM 과 OCP 최적 해 간의 이론적 연결을 검증했습니다.

중요성 및 주장

본 논문은 LIPM 과 최적 제어 간의 연결을 명시적이고 검증 가능하게 만든다고 주장합니다. 다음과 같이 논합니다:

1. LIPM 은 선택이 아닌 결과입니다: 풀랭크 지지 자세에서 각운동량 변화율을 최소화할 때, 가중치 $\alpha$가 충분히 높다면 진자 패턴이 자연스럽게 출현합니다.
2. 한계는 예측 가능합니다: $N=2$ 지지 자세에서 진자 가정의 붕괴는 제어기의 실패가 아니라 근본적인 기하학적 및 마찰적 한계입니다. 본 논문은 하한과 마찰 제약이 작용하는 임계 가속도에 대한 폐쇄형 표현식을 제공합니다.
3. 균형과 작업의 통합적 관점: 비율 $\lambda/\alpha$는 균형 중심 (LIPM 유사) 행동과 작업 중심 (민첩한) 행동 간의 전환을 설명하는 연속 매개변수를 제공합니다.

저자들은 폐루프 결과에 대해 겸손하게 접근하며, 정성적 행동 (바닥의 존재) 은 동적 시뮬레이션에서 유효하지만, 정적 이론과 폐루프 현실 간의 정량적 간격은 크며 모델링되지 않은 동역학 (탄성, 미끄러짐, 유한 제어 속도) 으로 인한 것일 수 있다고 지적합니다. 그들은 이러한 동적 간격을 해결한다고 주장하기보다는, 해당 동역학이 작동하는 정적 경계를 정의하는 데 그칩니다.