Equivalent Mechanical Models for Sloshing

기술 요약: 슬로싱(Sloshing)에 대한 등가 기계 모델

문제 정의
추진제 슬로싱은 우주선의 안정성과 제어에 있어 매우 중요한 과제로 남아 있으며, 특히 두 가지 뚜렷한 영역에서 나타난다: 높은 지향 안정성이 요구되는 미세 중력 환경(예: 행성 간 탐사 임무)과 발사, 재진입 및 착륙 시 발생하는 고가속도($g$) 조건이다. 제어 분석가들이 슬로싱 유발 섭동을 평가하기 위해 펜듈럼(pendulum) 기반 모델과 질량-스프링-댐퍼(mass-spring-damper) 모델을 표준 도구로 사용하고 있지만, 특히 고-$g$ 영역에서 명목 종방향 힘(nominal longitudinal forces)의 효과를 고려하여 이러한 표현 방식들 사이를 연결하는 엄밀한 수학적 정식화가 견고한 유도 항법 및 제어(GNC) 설계를 위해 필요하다.

방법론
본 논문은 슬로싱에 대한 등가 기계 모델을 비선형 동역학에서 제어 분석에 적합한 선형화된 표현으로 발전시키며 엄밀한 유도를 제시한다.

1. 비선형 정식화:

   • 본 연구는 질량 $m_P$와 길이 $\ell$을 가진 단일 구형 펜듈럼이 질량 $m_B$인 강체 $B$에 부착된 모델로부터 시작한다.
   • 저자들은 뉴턴 역학적 접근 방식을 사용하여 8자유도 시스템(강체 병진/회전 + 펜듈럼 각도)에 대한 결합 운동 방정식을 유도한다.
   • 유도 과정에서는 장력 $N$을 명시적으로 다루고 이를 제거하여 시스템 변수를 축소하며, 그 결과 $\mathcal{M}\ddot{\mathbf{x}} = \mathbf{u}$ 형태의 조밀한 행렬 형식을 얻는다.
   • 이 정식화는 여러 개의 펜듈럼(여러 탱크 또는 슬로싱 모드를 나타냄)을 수용할 수 있도록 확장되었으며, 효율적인 수치 해법을 위해 희소성(sparsity)을 활용하는 블록 구조 행렬을 사용한다.
   • 점성 댐핑 모델은 각속도에 비례하는 항을 추가함으로써 통합되었으며, 이는 강체 방정식 자체를 직접적으로 변화시키지 않으면서 펜듈럼 동역학에 영향을 미친다.

2. 선형화 및 모드 분석:

   • 탑재된 제어 시스템이 일반적으로 추진제 분포가 강체라고 가정하고 무게 중심을 추정한다는 점을 인식하여, 저자들은 시스템을 강체의 무게 중심이 아닌 명목 전역 무게 중심($\bar{G}$)을 기준으로 재정식화한다.
   • 시스템은 관성 좌표계와 정렬된 일정한 자세와 규정된 명목 종방향 힘 $\bar{F}_z$에 의해 특징지어지는 평형 상태를 중심으로 선형화된다.
   • 두 가지 경우가 분석된다: 관성 좌표계에서 정의된 힘과 물체 고정 좌표계(body-fixed frame)에서 정의된 힘이다.
   • 좌표 변환과 대수적 조작을 통해, 선형화된 동역학은 "모드 형태(modal form)"로 변환된다. 이는 펜듈럼 동역학을 두 개의 독립적인 평면 진동자로 분리하는 새로운 모드 좌표($\eta_\theta, \eta_\phi$)를 정의하는 과정을 포함한다.
   • 결과적인 시스템은 2차 선형 MIMO 시스템($M\delta\ddot{\mathbf{x}} + K\delta\mathbf{x} = \delta\mathbf{u}_{ext}$)으로 표현되며, 여기서 모드 참여 행렬(modal participation matrix)과 명목 힘에 의해 유도된 강성(stiffness) 항을 명시적으로 식별한다.

3. 질량-스프링-댐퍼 모델과의 등가성:

   • 본 논문은 슬로싱 질량이 강성 $k$를 가진 스프링에 의해 횡방향 이동이 제한되는 질량-스프링-댐퍼 시스템의 선형화된 방정식을 유도한다.
   • 스프링 강성을 $k = m_P \omega_0^2$ (여기서 $\omega_0$는 펜듈럼 모델에서 유도된 고유 진동수)로 설정함으로써, 저자들은 선형화된 질량-스프링-댐퍼 동역학이 선형화된 펜듈럼 동역학과 수학적으로 동일함을 입증한다.
   • 분석 결과, 명목 종방향 힘 $\bar{F}_z$의 존재는 펜듈럼이 스프링과 유사한 진동 거동을 보이기 위한 필수 요소임을 강조한다. 즉, 이 힘이 없는 경우($\bar{F}_z=0$), 펜듈럼 모델은 0-g 슬로싱을 표현하기 위해 수정(예: 세 번째 모드 자유도 추가)이 필요하다.

주요 결과

• 수학적 등가성: 본 연구는 고-$g$ 조건에서 선형화된 펜듈럼 모델과 고전적인 질량-스프링-댐퍼 모델 사이의 정확한 등가성을 공식적으로 확립하였으며, 이를 통해 명목 힘이 펜듈럼 정식화에서 어떻게 복원 메커니즘(restoring mechanism)으로 작용하는지 명확히 하였다.
• 모드 표현: 유도된 모드 표현은 강성 행렬에서 슬로싱 모드가 강체 동역학으로부터 디커플링(decoupling)되도록 하며, 결합은 모드 참여 행렬을 통해 질량 행렬에서만 나타나도록 한다.
• 검증:
  • 비선형 검증: 비선형 모델은 MATLAB Simscape Multibody 구현(4개의 펜듈럼이 달린 강체 허브)을 통해 검증되었다. 시계열 시뮬레이션 결과, 응답 진폭의 0.002% 미만 수준에서 매우 우수한 일치도를 보였으며, 불일치는 오직 수치 적분 방식의 차이에서 기인한 것으로 분석되었다.
  • 선형 검증: 해석적으로 선형화된 모델과 선형화된 Simscape 플랜트를 비교한 Bode plot은 관성 및 물체 고정 힘 입력 모두에 대해 전 주파수 범위에서 매우 우수한 일치도를 보였다. 관성 좌표계 케이스에서의 미세한 저주파 불일치는 Simscape 선형화 과정의 수치적 정확도 한계로 판단된다.

의의 및 주장
본 논문은 고가속도 환경에서의 추진제 슬로싱 모델링을 위한 엄밀한 수학적 기초를 제공한다고 주장한다. 주요 기여는 다중 펜듈럼 시스템에 대한 비선형 운동 방정식의 명시적 유도와 명목 종방향 힘을 고려한 후속 선형화이다. 펜듈럼 모델과 질량-스프링-댐퍼 모델 사이의 수학적 등가성을 증명함으로써, 본 연구는 후자의 제어 설계 활용을 정당화하는 동시에 이를 더 근본적인 펜듈럼 동역학에 근거하게 한다. 저자들은 제안된 모델이 슬로싱 각도가 합리적인 범위 내에 있는 실제 우주선 응용 분야에 적합함을 강조하며, 검증을 통해 시계열 시뮬레이션 및 주파수 영역 분석 모두에서 도출된 정식화의 정확성을 확인하였다.