[논문 리뷰] Attitude Maneuvers of a Rigid Spacecraft in a Circular Orbit

Taeyoung Lee, Melvin Leok외 1명|ArXiv.org|2005. 09. 14.

Spacecraft Dynamics and Control참고 문헌 2인용 수 8

한 줄 요약

이 논문은 원형 궤도를 도는 강체 위성의 대칭각 자세 제어를 위해 SO(3)의 회전 행렬과 리 군 변분 적분기를 사용하여 전역적으로 유효하고 특이점이 없는 기하학적 모델을 제시한다. 비선형 경계값 문제 및 최소 임펄스 최적 제어 문제를 정확하게 해결하기 위한 계산 방법(근 찾기 및 순차적 목적함수 최적화)을 개발하였으며, 리 대수 so(3)에서 유도된 민감도 도함수를 포함한다. 주요 기여는 수치적으로 안정적이고 기하학적으로 정확한 접근 방식으로, 위성 자세 역학의 본질적 구조를 유지한다는 점이다.

ABSTRACT

A global model is presented that can be used to study attitude maneuvers of a rigid spacecraft in a circular orbit about a large central body. The model includes gravity gradient effects that arise from the non-uniform gravity field and characterizes the spacecraft attitude with respect to the uniformly rotating local vertical local horizontal coordinate frame. An accurate computational approach for solving a nonlinear boundary value problem is proposed, assuming that control torque impulses can be applied at initiation and at termination of the maneuver. If the terminal attitude condition is relaxed, then an accurate computational approach for solving the minimal impulse optimal control problem is presented. Since the attitude is represented by a rotation matrix, this approach avoids any singularity or ambiguity arising from other attitude representations such as Euler angles or quaternions.

연구 동기 및 목표

• 원형 궤도를 도는 강체 위성의 대칭각 자세 제어를 위한 전역적으로 유효하고 기하학적인 모델을 개발한다. 이때 중력 기울기 효과를 고려한다.
• 기존 자세 표현 방식(예: 오일러 각도, 허니우스)의 한계를 극복하기 위해 SO(3) 내의 회전 행렬을 사용하여 특이점과 모호함을 제거한다.
• 두 가지 유형의 자세 제어 문제를 수립하고 해결한다: (1) 경계값 문제를 통한 고정된 종료 자세, (2) 완화된 종료 조건을 갖는 최소 임펄스 제어.
• 리 대수 so(3)에서 분석적으로 유도된 민감도 도함수를 통해 계산 정확도를 확보하고, 역학의 기하학적 구조를 유지한다.
• 반지름에서 궤도면 수직 방향으로의 1/4 궤도 제어를 수치적으로 해결하여 방법의 효과성을 입증한다.

제안 방법

• 위성 자세를 몸체 기준에서 국지 수직-지역 수평(지역 수평 기준, LVLH) 기준으로 변환하는 데 사용되는 $ R^{bl} \in \mathrm{SO(3)} $ 인 회전 행렬로 표현한다.
• SO(3)의 기하학적 구조를 유지하고 관측 가능한 양을 보존하는 리 군 변분 적분기를 사용하여 연속적인 자세 역학을 모델링한다.
• 중력 잠재에너지에서 유도된 중력 기울기 모멘트를 포함하여, LVLH 기준에서 비균일한 중력장 영향을 기록한다.
• 상태 전이 행렬을 통해 분석적으로 유도된 민감도 도함수를 사용한 근 찾기 알고리즘을 통해 비선형 경계값 문제를 해결한다.
• 성능 지수와 제약 조건의 기울기를 so(3)에서 분석적으로 유도하여 순차적 목적함수 최적화(SQP)를 사용해 최소 임펄스 최적 제어 문제를 해결한다.
• 민감도 도함수를 리 대수 $ \mathfrak{so}(3) $의 형태로 표현하여 기하학적 일致성과 오일러 각도 또는 허니우스에서 발생하는 특이점 회피를 보장한다.

실험 결과

연구 질문

• RQ1 원형 궤도를 도는 강체 위성의 대칭각 자세 제어를 오일러 각도나 허니우스에서 기인하는 특이점과 모호함 없이 어떻게 모델링할 수 있는가?
• RQ2 위성 자세 제어에서 비선형 경계값 문제를 해결하기 위한 가장 정확하고 기하학적으로 일致하는 계산적 접근은 무엇인가?
• RQ3 주어진 자세 제어를 위해 필요한 최소 총 임펄스는 어떻게 계산할 수 있으며, 종료 제약 조건을 충족할 수 있는가?
• RQ4 성능 지수 및 제약 조건의 민감도 도함수는 최적화 과정의 수렴성과 정확도에 어떻게 기여하는가?
• RQ5 리 군 변분 적분기를 사용할 경우, SO(3)의 본질적 기하학적 구조가 장기 시뮬레이션에서 얼마나 잘 유지되는가?

주요 결과

• 계산 방법은 제약 위반 값 $ 4.80 \times 10^{-14} $ 이내로 최소 임펄스 최적 제어 문제를 성공적으로 해결하여 높은 수치 정확도를 입증하였다.
• 최적화된 초기 각운동량은 $ \Pi_0 = [-2.915, -2.347, -2.734]^T $ 로, 종료 시 각운동량 $ -\Pi_N = [-0.343, 2.686, 2.734]^T $ 를 얻어 원하는 자세를 달성하였다.
• 최종 회전 행렬 $ R^{bl}_N $ 은 종료 방향 벡터 $ \Lambda_N = [0, 0, -1]^T $ 를 생성하였으며, 기준 프레임 방향에 따라 보정된 후 목표값 $ \Lambda_{N_d} = [0, -1, 0]^T $ 와 일치하였다.
• 성능 지수 $ \mathcal{J} = 6.771 $ 을 달성하여 제어 동안의 총 임펄스 크기를 나타내었다.
• 이 방법은 기하학적으로 정확하고 수치적으로 효율적이며, 민감도 도함수를 $ \mathfrak{so}(3) $ 에서 분석적으로 유도하여 특이점 없이 작동한다.
• 모든 초기 각속도 $ \omega_0 = 0 $ 인 경우 자유 강체 운동으로 수렴하여 고전 역학과의 일致성을 확인하였다.