[논문 리뷰] Convex Optimization of Linear Impulsive Rendezvous

Boris Benedikter, Alessandro Zavoli|arXiv (Cornell University)|2019. 12. 17.

Spacecraft Dynamics and Control참고 문헌 26인용 수 5

한 줄 요약

이 논문은 타임-픽스드 선형 인パル시브 레이던수스 문제를 풀기 위해 콘벡스 최적화 방법을 제안한다. 타우셔너-헤임플 역학을 사용하며, 충격 위치를 격자에 따라 공간 이산화함으로써 비콘벡스 최적 제어 문제를 콘벡스 문제로 변환한다. 이 방법은 최소한의 계산 시간으로 전역 수렴을 보장하며, 최적의 충격 수와 시기를 정확히 파악한다. 기존 기준값과의 검증 결과 근사적으로 완벽한 일致를 보였다.

ABSTRACT

This note aims at providing a concise and self-contained document that describes a clear and easy-to-understand method, that could be useful for a reader that is approaching the linear-impulsive rendezvous topic for the first time, but that is also flexible enough to accommodate problem variations, such as additional constraints like bounded-magnitude DV, or direction constraints with minor modifications.

연구 동기 및 목표

• 타원 궤도에서 시간-픽스드 선형 인파르시브 레이던수스 문제를 해결하기 위한 신뢰성 있고 계산 효율적인 방법을 개발한다.
• 사전에 충격 수와 위치에 대한 가정이 필요한 전통적 방법의 한계를 극복한다.
• 초기화에 의존하지 않고 추진제 소비(Δv)를 전역 최적화할 수 있도록 한다.
• 자율 레이던수스와 같은 시간이 중요한 우주 응용 분야에서 높은 정확도와 빠른 수렴성을 확보한다.
• 콘벡스 전환에서 격자 해상도, 해의 정확도, 수치적 안정성 간의 상호 상충 관계를 조사한다.

제안 방법

• 독립 변수(진공각)를 유한한 격자 노드 집합으로 이산화함으로써 원래의 최적 제어 문제를 콘벡스 최적화 문제로 변환한다.
• 충격 발생 위치는 오직 이 격자 노드에서만 허용되며, 상태 전이 행렬을 사용해 연속 시간 역학을 대체로 대수적 상태 전이 시퀀스로 변환한다.
• 수치 조건과 안정성을 향상시키기 위해 변환된 변수(예: 스케일링된 위치 및 속도)를 사용해 문제를 재구성한다.
• 결과적으로 유도된 콘벡스 문제를 내부점 방법으로 해결하며, 초기화에 관계없이 전역 최적해로 수렴함을 보장한다.
• 최적화 과정에서 충격 수가 암묵적으로 결정되므로 반복적 조정이 필요 없어진다.
• 반올림 오차와 절삭 오차를 줄이기 위해 문제 스케일링과 무차원화를 적용한다. 특히 정밀한 격자에서 효과적이다.

실험 결과

연구 질문

• RQ1 타우셔너-헤임플 역학을 사용한 타원 궤도에서의 시간-픽스드 선형 인파르시브 레이던수스 문제에 대해 콘벡스 공식화를 유도할 수 있는가?
• RQ2 격자 기반 충격 위치 이산화가 초기화에 의존하지 않고 최적해로의 전역 수렴을 가능하게 하는가?
• RQ3 격자 해상도의 선택이 Δv 비용과 충격 시점 정확도에 어떤 영향을 미치는가?
• RQ4 격자 밀도가 증가함에 따라 계산 비용과 해의 품질 간의 상충 관계는 어떻게 변화하는가?
• RQ5 문제 스케일링과 허용 오차 제어를 통해 수치적 문제(예: 충격 분산)를 완화할 수 있는가?

주요 결과

• 콘벡스 최적화 방법은 원형 궤도에서 원형 궤도로의 레이던수스에 대해 전역 최적의 이중 충격 해를 정확히 복원하였으며, Arzelier 등이 제시한 기준해와 100% 일致를 보였다.
• 전반적인 Δv 비용은 조밀한 격자조차도 빠르게 최적값으로 수렴하며, 8~16개의 격자 노드로도 거의 최적 성능을 달성한다.
• 계산 시간은 격자 노드 수에 거의 선형적으로 증가하지만, 콘벡스 솔버의 효율성 덕분에 대규모 격자에서도 매우 낮은 수준을 유지한다.
• 미세한 격자에서는 유한 정밀도와 솔버 허용 오차로 인해 수치적 불안정성이 발생하여 충격이 여러 노드에 걸쳐 퍼지는 현상이 발생한다.
• 중간에서 조밀한 격자 수준이 정확도와 계산 효율성의 최적 균형을 이룹니다. 오차 누적이 방지되며 고해상도 해의 정밀도를 유지합니다.
• 보장된 수렴성과 낮은 계산 오버헤드 덕분에 이 방법은 실시간 응용에 적합하며 강건하다.