⚡ electrical engineering

B-spline shaping for low-thrust interplanetary rendezvous

이 논문은 궤적을 매개화하기 위해 클램프된 B-스플라인(clamped B-splines)을 활용하고 제어 가속도를 대수적으로 회복함으로써 최적화 문제를 유한 차원 비선형 계획법으로 축소하며, 다양한 목표 시나리오에서 고차 호도그래프 형상화(high-order hodographic shaping)보다 델타-V 측면에서 일관되게 우수한 성능을 보이는 저추력 행성 간 랑데우를 위한 계산 효율적인 형상 기반 방법을 소개한다.

원저자: Julio C. Sanchez

게시일 2026-07-28

📖 6 분 읽기🧠 심층 분석

CC BY 4.0

원저자: Julio C. Sanchez

원본 논문은 CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) 라이선스로 제공됩니다. ✨ 이것은 아래 논문에 대한 AI 생성 설명입니다. 저자가 작성하거나 승인한 것이 아닙니다. 기술적 정확성을 위해서는 원본 논문을 참조하세요. 전체 면책 조항 읽기

당신은 아주 작고 효율적인 로봇 자동차가 광활하고 보이지 않는 우주의 바다를 가로질러 가도록 안내하려고 한다고 상상해 보십시오. 이것은 급가속을 위해 가스 페달을 세게 밟는 일반적인 자동차가 아닙니다. 이것은 부드럽고 지속적인 미풍을 받는 돛단배와 같은 "저추력(low-thrust)" 탈것입니다. 이 자동차는 빠르게 멈추거나 출발할 수 없습니다. 대신, 지구에서 다른 행성까지 가기 위해 몇 달, 혹은 몇 년 동안 부드럽게 밀어붙여야 합니다. 과학자들의 큰 과제는 이 느리고 꾸준한 추진력을 위한 완벽한 경로를 찾아내는 것입니다. 그들은 연료(무겁고 운반하기 비싼)를 가장 적게 사용하면서도 합리적인 시간 내에 임무를 완수할 수 있는 경로를 찾아내야 합니다. 만약 경로가 잘못된다면, 로봇은 도착하기도 전에 연료가 바닥나거나, 임무가 실패할 정도로 너무 오래 걸릴 수도 있습니다. 이 "슬로 모션" 경로를 찾는 것은 마치 움직이는 조각들로 이루어진 거대한 3차원 퍼즐을 푸는 것과 같으며, 엔지니어들이 실제로 임무를 수행할 수 있도록 솔루션을 빠르게 찾아내야 합니다.

이 논문은 "B-스플라인 셰이핑(B-spline shaping)"이라는 기술을 사용하여 이 퍼즐을 해결하는 영리한 새로운 방법을 소개합니다. B-스플라인이 무엇인지 생각해보면, 그것은 엔지니어들이 원하는 대로 어떤 매끄러운 곡선으로든 구부릴 수 있는 유연한 디지털 자와 같습니다. 엔지니어들이 수년간 엔진이 가하는 모든 아주 작은 추진력을 일일이 계산하는 대신, 이들은 이 유연한 자들을 사용하여 우주선의 전체 경로를 먼저 그려냅니다. "미분 평탄성(differential flatness)"이라는 특별한 수학적 트릭 덕분에, 일단 경로의 형태를 결정하면, 그들은 경로를 따라 유지하기 위해 엔진이 매 순간 얼마나 세게 밀어야 하는지를 즉각적으로 계산할 수 있습니다. 연구진은 이 방법이 더 오래되고 복잡한 경로 그리기 방식(호도그래프 셰이핑, hodographic shaping이라 불리는)을 대상으로 화성, 수성, 근처의 소행성, 그리고 혜성을 향한 여정을 테스트했습니다. 결과는 이 새로운 "유연한 자" 방식이 표준 구성에서 계산 속도가 더 빠를 뿐만 아니라, 특히 수성이나 혜성처럼 까다로운 목적지에 대해서는 훨씬 더 적은 연료를 사용하는 경로를 찾아낸다는 것을 시사합니다. 이는 마치 다른 사람들이 계산하느라 너무 바빠서 알아채지 못한 미로 속의 지름길을 찾아내는 것과 같습니다.

문제: 느리고 꾸준한 우주 경주

우주 여행은 보통 느리다는 오명을 쓰곤 하지만, 저추력 엔진은 사실 효율성의 챔피언입니다. 거대한 불꽃 폭발과 함께 튀어 나가는 화학 로켓과 달리, 전기 추진기(Dawn이나 Hayabusa 임무에 사용된 것과 같은)는 아주 작은 지속적인 추진력을 제공합니다. 이것은 폭발적으로 출발하는 단거리 선수와 완벽하고 꾸준한 페이스를 유지하는 마라톤 선수의 차이와 같습니다. 문제는 추진력이 매우 부드럽기 때문에, 우주선이 목적지에 도달하기 위해 태양 주위를 여러 번 나선형으로 돌아야 한다는 점입니다. 이로 인해 임무 설계는 수학적 악몽이 됩니다. 당신은 연료가 떨어지기 전에 이 느리게 움직이는 물체를 몇 년 동안 어떻게 조종할지를 정확히 계산해야 합니다.

전통적으로 과학자들은 여 journey를 아주 작은 단계로 나누어 각 단계에서의 힘을 계산하거나, 정답을 추측하고 계속해서 다듬어가는 복잡한 "간접적" 수학을 사용하여 이를 해결하려고 노력해 왔습니다. 이러한 방법들은 종종 느리고 까다로우며, 초기 추측이 완벽하지 않으면 막힐 수 있습니다. 이 논문의 저자들은 이러한 저추력 엔진을 위한 궤적을 더 빠르고 유연하게 설계할 수 있는 방법을 원했으며, 슈퍼컴퓨터가 며칠 동안 연산에 매달리지 않고도 다양한 목표에 대해 빠르게 좋은 솔루션을 생성할 수 있는 방법을 원했습니다.

해결책: 디지털 자로 경로 그리기

저자들은 B-스플라인 셰이핑이라고 불리는 방법을 제안합니다. 이를 이해하기 위해, 당신이 종이 위에 매끄러운 곡선을 그리고 있다고 상상해 보십시오. 당신은 선 전체를 한 번에 그리는 것이 아니라, 몇 개의 "제어점(control points)"(판 위의 핀과 같은 것)을 배치하고 그 사이를 통과하는 유연한 띠(스플라인)를 늘립니다. 핀을 움직임으로써 곡선의 모양을 바꿀 수 있습니다.

이 논문에서 연구진은 **클램프드 B-스플라인(clamped B-spline)**이라는 특정 유형의 디지털 곡선을 사용합니다. 그들은 이 곡선들을 사용하여 우주선의 경로, 즉 태양으로부터의 거리, 각도, 그리고 태양계의 평면을 기준으로 한 높낮이를 3차원으로 그립니다.

여기에는 마법 같은 트릭이 있습니다. 저추력 엔진의 물리 법칙은 미분 평탄성이라는 특별한 성질을 가지고 있습니다. 이것은 경로의 형태와 변화 속도를 알고 있다면, 그 경로를 유지하기 위해 엔진이 적용해야 할 정확한 힘을 즉시 계산할 수 있다는 것을 의미하는 멋진 표현입니다. 매 초마다 비행을 시뮬레이션할 필요가 없습니다. 그냥 경로를 그리면, 수학이 엔진 설정을 알려줍니다.

연구진은 곡선의 시작과 끝이 지구와 목표 행성의 위치 및 속도와 완벽하게 일치하도록 "핀"들을 설정했습니다. 그런 다음, 가장 적은 연료를 사용하는 경로를 찾기 위해 컴퓨터가 중간 핀들의 위치를 최적화하도록 했습니다. 시작과 끝은 수학에 의해 고정되어 있기 때문에, 컴퓨터는 중간 부분에만 집중하면 되며, 따라서 문제는 훨씬 작아지고 빠르게 해결됩니다.

테스트 드라이브: 화성, 수성, 그리고 그 너머

이 새로운 방법이 실제로 작동하는지 확인하기 위해, 팀은 네 가지 다른 임무에 대한 시뮬레이션을 실행했습니다:

지구에서 화성까지: 고전적인 붉은 행성 여행.
지구에서 수성까지: 수성이 태양의 중력 우물 깊숙이 있기 때문에 매우 어려운 여행.
지구에서 소행성 1989 ML까지: 근지구 소행성.
지구에서 혜성 템펠 1까지: 기묘하게 기울어진 궤도를 가진 혜성.

그들은 자신들의 B-스플라인 방법을 다른 수학적 구성 요소를 사용하는 잘 알려진 고품질 기술인 "고차 호도그래프 셰이핑(high-order hodographic shaping)" 방법과 비교했습니다. 그들은 수천 번의 시뮬레이션을 실행하여 서로 다른 출발 날짜와 여행 기간을 테스트했으며, 이는 "포크숍 플롯(porkchop plots)"(모든 가능한 여행에 대한 연료 비용을 보여주는 지도)이라 불리는 것을 생성했습니다.

결과: 더 빠르고 더 가볍게 (단, 주의사항 있음)

결과는 인상적이었지만, 몇 가지 중요한 뉘앙스가 있었습니다. 전반적으로 B-스플라인 방식은 전통적인 방식보다 일반적으로 더 적은 연료를 사용하는 경로를 찾아냈지만, 그 개선 정도는 목적지에 따라 달랐습니다.

연료 절감: 새로운 방식은 더 어려운 목표들에 대해 필요한 연료( $\Delta V$ , 즉 속도 변화량으로 측정됨)를 크게 줄였습니다. 수성 여행의 경우, "중간값" 연료 비용이 약 45.4% 감소했습니다. 소행성 1989 ML의 경우 37.9% 감소했으며, 혜성 템펠 1의 경우 31.6% 감소했습니다. 화성의 경우 개선 효과는 더 완만했습니다. 중간값 연료 비용은 21.8% 감소했지만, 기존 방식이 찾아낼 수 있었던 가장 좋은(최소) 경로와 비교했을 때 절대적인 최적 경로는 단 1.4% 더 나은 수준이었습니다. 이는 새로운 방식이 평균적으로는 일관되게 더 낫지만, 화성을 위한 "완벽한" 지름길은 이미 기존 방식이 찾을 수 있는 수준에 매우 근접해 있었음을 강조합니다.
속도: 방법의 속도는 "자"의 복잡성이 어떻게 설정되느냐에 따라 크게 달라집니다. 연구진은 자신들의 표준 구성(10개의 제어점 사용)이 기존 방식보다 계산 시간이 약 30~33% 적게 걸려 훨씬 빠르다는 것을 발견했습니다. 기존 방식이 0.158초 걸렸던 여행에 대해 새로운 방식은 0.109초밖에 걸리지 않았습니다. 그러나 추가적인 연료 절감을 위해 복잡성을 40개의 제어점으로 늘리면, 이 방식은 표준 버전보다 11~13배 더 느려집니다. 즉, 트레이드오프가 존재합니다: 단순한 버전은 속도의 제왕이지만, 복잡한 버전은 시간이 더 오래 걸리는 연료 절약형입니다.
"스위트 스팟(Sweet Spot)": 연구진은 디지털 자에 있는 다양한 수의 "핀"(제어점)을 테스트했습니다. 그들은 10개의 제어점과 5차(quintic) 곡선을 사용하는 것이 완벽한 균형임을 발견했습니다. 이 방식은 빠르고, 연료를 매우 적게 사용하며, 너무 많은 컴퓨팅 파워를 요구하지 않았습니다. 만약 40개의 제어점을 사용한다면 (특히 최상의 시나리오에서) 아주 약간의 연료 절감을 더 얻어낼 수 있지만, 컴퓨터는 이를 위해 약 10배 더 오래 작동해야 했습니다.

이것이 미래에 의미하는 바

이 논문은 B-스플라인 접근 방식이 임무 설계의 초기 단계에서 강력한 도구가 될 수 있음을 시사합니다. 이를 통해 엔지니어들은 느리고 복잡한 계산에 발이 묶이지 않고도 수천 가지의 가능한 발사 날짜를 빠르게 훑어보고 최적의 날짜를 찾을 수 있습니다.

하지만 저자들은 이 연구가 하지 않은 것에 대해서도 주의를 기울였습니다. 그들은 주요 계산에 연료가 타서 없어지는 무게(질량 감소)를 포함하지 않았으며, 실제 우주선에 대해 테스트하지도 않았습니다. 그들이 찾은 솔루션은 "국소적 최적해(locally optimal)"입니다. 즉, 그들이 테스트한 특정 수학적 범주 내에서는 최고이지만, 수학의 전체 우주에서 가능한 절대적인 최적의 경로는 아닐 수 있습니다. 또한, 이 방법은 엔진이 어떤 방향으로든 어떤 강도로도 밀 수 있다고 가정하는데, 이는 실제 하드웨어에서는 항상 사실이 아닙니다.

그럼에도 불구하고, 이 연구는 전체 여정의 더 나은 그림을 먼저 그리는 데 유연한 수학적 "자"를 사용함으로써, 우리가 태양계의 먼 구석으로 로봇 탐사선을 보내는 더 좋고, 더 빠르고, 더 저렴한 방법을 찾을 수 있음을 보여줍니다. 이는 때때로 복잡한 문제를 해결하는 최선의 방법은 모든 작은 세부 사항을 계산하는 것이 아니라, 전체 여정에 대한 더 나은 그림을 그리는 것임을 상기시켜 줍니다.

기술 요약: 저추력 행성 간 랑데부(Rendezvous)를 위한 B-스플라인 형상 설계

문제 정의
본 논문은 고정된 시간 내의 최소 $\Delta V$ 저추력 행성 간 랑데부 궤도의 예비 설계를 다룬다. 목적은 정해진 수의 궤도 회전( $N$ ) 동안 두 천체 사이를 이동하는 우주선의 제어 가속도 크기의 시간 적분( $\int_{t_0}^{t_f} \|u\| dt$ )을 최소화하는 것이다. 문제는 태양 중심 원통 좌표계( $r, \theta, z$ )를 사용하여 정식화되었으며, 경계 조건(위치 및 속도)은 출발 및 도착 시간의 에페메리스(ephemerides)로부터 도출된다. 복잡한 이점-이점 경계값 문제(two-point boundary-value problems)를 해결하거나 동역학적 제약 조건이 포함된 대규모 비선형 계획법을 요구하는 전통적인 최적 제어 방식과 달리, 본 연구는 신속한 예비 설계를 위해 계산 효율성과 충분한 해의 품질을 결합한 파라미터화 기법을 추구한다.

방법론
제안된 방법은 저추력 원통형 동역학의 **미분 평탄성(differential flatness)**을 활용한 형상 기반 접근 방식을 사용한다. 핵심 혁신은 태양 중심 원통 좌표를 **클램프된 B-스플라인(clamped B-splines)**으로 직접 파라미터화하는 것이다.

미분 평산성: 원통형 저추력 동역학은 위치 좌표 $q = [r, \theta, z]^T$ 에 대해 미분 평탄하다. 이 특성을 통해 제어 가속도 $u$ 를 궤적과 그 1차 및 2차 미분으로부터 대수적으로 복원할 수 있다( $u \equiv u(q, \dot{q}, \ddot{q})$ ). 결과적으로 운동 방정식은 수치적 전파(numerical propagation)를 거치지 않고도 구조적으로 만족된다.
B-스플라인 파라미터화: 궤적 $q(\tau)$ 는 B-스플라인 기저 함수와 제어점(control points)의 선형 결합으로 표현된다. 독립 변수는 정규화된 시간 $\tau \in [0, 1]$ 이다.
클램프된 경계 조건: 양 끝점( $\tau=0$ 및 $\tau=1$ )에서 매듭(knots)을 반복함으로써, 본 방법은 위치 및 속도 경계 조건을 해석적으로 만족시킨다. 처음 두 개와 마지막 두 개의 제어점은 초기 및 최종 상태와 일치하도록 고정되며, 오직 내부 제어점들만이 최적화 변수로 남는다.
최적화 정식화: 연속 최적 제어 문제는 유한 차원 비선형 계획법(NLP) 문제로 축소된다. 결정 변수는 자유로운 내부 B-스플라인 제어점들이다. 목적 함수( $\Delta V$ )는 대수적으로 도출된 제어 크기의 수치 적분을 통해 평가된다. 동역학 콜로케이션(collocation) 제약 조건은 필요하지 않다.
구현: NLP는 심볼릭 CasADi 포뮬레이션을 갖춘 IPOPT 알고리즘을 사용하여 해결된다. 초기 추측값은 수정된 등각 요소(modified equinoctial elements)에서의 5차 다항식 프로파일에 B-스플라인 제어점을 피팅하여 생성된다.

주요 기여
본 논문은 두 가지 주요 기여를 한다.

새로운 형상 기반 정식화: 원통 좌표의 클램프된 B-스플라인 파라미터화와 동역학의 미분 평탄성을 결합하였다. 이를 통해 제어 가속도를 대수적으로 복원하고 경계 조건을 해석적으로 만족시킬 수 있어, 운동 방정식의 수치적 전파 없이 문제를 정적인 NLP로 축소한다. 대상별로 특화된 기저 함수를 직접 구축해야 하는 기존의 형상 기반 방식과 달리, 이 방식은 표준 B-스플라인 하이퍼파라미터(차수 및 제어점 수)를 통해 궤적의 유연성을 체계적으로 조정할 수 있다.
체계적인 벤치마킹: 본 방법은 네 가지 서로 다른 랑데부 시나리오(지구-화성, 수성, 근지구 소행성 1989 ML, 혜성 Tempel 1)에 대해 고차 호도그래프(hodographic) 형상화 기법(선도적인 형상 기반 기법)을 대상으로 엄격하게 테스트되었다.

결과
출발 날짜와 전이 시간의 미세 격자에 대해 수치 실험을 수행하였다. 결과에 따르면, B-스플라인 방식은 모든 테스트 대상 및 격자 캠페인에서 호도그래프 벤치마크에 비해 해의 품질 측면에서 일관되게 우수한 성능을 보였으며, 특정 구성에서는 계산 효율성 또한 높았다.

해의 품질: 테스트된 모든 대상에 대해, 모든 B-스플라인 구성은 호도그래프 벤치마크 대비 평균, 중앙값 및 최소 누적 $\Delta V$ $Δ V$ 를 감소시켰다.
- 화성: 10개의 제어점을 가진 5차 B-스플라인은 중앙값 $\Delta V$ 를 21.8%, 최소값을 1.4% 감소시켰다.
- 수성: 10개의 제어점을 가진 5차 B-스플라인은 중앙값 $\Delta V$ 를 45.4%, 최소값을 4.8% 감소시켰다.
- 1989 ML: 10개의 제어점을 가진 5차 B-스플라인은 중앙값 $\Delta V$ 를 37.9%, 최소값을 9.3% 감소시켰다.
- Tempel 1: 10개의 제어점을 가진 5차 B-스플라인은 중앙값 $\Delta V$ 를 31.6%, 최소값을 10.3% 감소시켰다.
계산 효율성: 10개의 제어점을 가진 5차 B-스플라인은 호도그래프 방식보다 최적화 시도당 약 30~~33% 더 빨랐다. 40개의 제어점을 가진 구성은 더 낮은 최소 $\Delta V$ 값을 나타냈으나, 10개 제어점 구성보다 계산 비용이 현저히 높았다(7~~13배 느림).
제어 프로파일: 더 높은 차수의 스플라인(5차 vs 3차)은 더 부드러운 제어 전이를 생성하였고, 더 많은 제어점 수(40개 vs 10개)는 각각 더 국부적인 추력 아크(충격-휴지 형태의 구조)를 생성하였다.

의의 및 주장
본 논문은 제안된 방법이 저추력 행성 간 전이를 위한 예비 설계에 있어 견고하고 효율적인 도구를 제공한다고 주장한다. 구체적으로, 10개의 제어점을 가진 5차 B-스플라인이 해의 품질과 계산 효율성 사이의 유리한 절충안임을 확인하였으며, 이는 광범위한 예비 설계 캠페인에 적합하다. 대상별로 특화된 기저 함수를 수동으로 구축할 필요 없이 다양한 타겟을 처리할 수 있는 능력은 기존의 형상 기반 기술에 비해 중요한 장점으로 강조된다.

저자들은 본 결과가 유한 차원 형상 공간(finite-dimensional shape spaces) 및 무제약 가속도 문제(질량 변화 및 엄격한 추력 제한이 포함되지 않음)로 제한된다는 점을 겸허히 언급하였다. 도출된 해는 선택된 B-스플라인 공간 내에서 국부 최적해이며, 원래의 연속 문제에 대한 전역 최적성을 입증하는 것은 아니다. 그러나 형상화된 궤적은 고충실도 직접 또는 간접 최적 제어 방법을 위한 효과적인 초기 추측값으로서 제시된다.