제한거리 자율비행을 위한 인버스 시뮬레이션 비행제어 사전계산 알고리즘

초록

본 논문은 6자유도(6‑DOF) 비행역학 방정식을 기반으로, 목표 궤적을 입력하면 필요한 엔진 추력 및 3축 제어면(에일러론·엘리베이터·러더) 각도를 역으로 계산하는 InvSim 알고리즘을 제시한다. 기호 연산, RK4 적분, 유한차분을 결합해 시간 연속적인 제어값을 도출하고, 미라지 III 전투기의 더블 롤(Double‑Roll) 기동을 사례로 시뮬레이션 결과를 검증한다.

상세 분석

이 연구는 고전적인 6‑DOF 비행역학 모델을 확장하여, 항공기의 몸체축, 관성축, 풍향축을 동시에 고려한 통합 좌표계 체계를 구축한다. 특히, 구면 좌표계(방위각·고도각)와 비행각(받음각·측면슬립각)을 명시적으로 포함함으로써, 비대칭 구조나 비정상 비행 상태에서도 정확한 동역학을 기술한다는 점이 돋보인다. 기존의 순방향 시뮬레이션은 주어진 제어 입력에 대해 궤적을 예측하지만, 본 논문이 제안하는 ‘역 시뮬레이션(Inverse Simulation, InvSim)’은 목표 궤적을 사전 정의하고, 이를 만족시키기 위한 최소한의 제어 입력을 역으로 계산한다는 근본적인 패러다임 전환을 시도한다.

InvSim 알고리즘의 핵심은 다음과 같다. 첫째, 목표 궤적을 시간에 대한 3차원 위치와 롤 각도로 정의하고, 이를 미분하여 속도·가속도 정보를 얻는다. 둘째, 6‑DOF 방정식에서 제어 입력(추력·에일러론·엘리베이터·러더)과 상태 변수(속도·각속도·오일러 각) 사이의 관계를 기호적으로 전개한다. 여기서 기호 연산(sympy 등)을 활용해 비선형 항을 명시적으로 분리하고, 제어 입력을 선형화할 수 있는 형태로 재배열한다. 셋째, 재배열된 방정식을 유한차분(FDM) 형태로 이산화하여, 각 시간 스텝에서 제어 입력을 직접 해석적으로 구한다. 넷째, 이 해를 명시적 4차 루젠쿠타(RK4) 적분기에 삽입해 시간 전진을 수행함으로써, 연속적인 제어 신호 시퀀스를 생성한다.

알고리즘 구현 시 주목할 점은 수치적 안정성이다. 비행역학은 강한 비선형성(예: 고속·고각도 비행)과 급격한 제어 입력 변화를 포함하므로, RK4와 FDM의 결합이 적절히 시간 스텝을 선택하지 않으면 발산하거나 물리적으로 불가능한 제어값(예: 과도한 추력·극단적인 면적 변위)이 산출될 위험이 있다. 논문은 이를 완화하기 위해, 제어 입력에 물리적 한계(추력 상한·면적 각도 범위)를 강제하고, 뉴턴-라프슨 방식의 보정 루프를 추가해 수렴성을 확보한다는 점을 강조한다.

시연 사례로 선택된 미라지 III 전투기의 더블 롤(Double‑Roll) 기동은 고속·고각도 회전과 동시에 급격한 피치·요각 변화를 요구한다. 목표 궤적을 3초 구간에 360° 롤을 두 번 수행하도록 설계하고, InvSim이 도출한 엔진 추력·에일러론·엘리베이터·러더 시퀀스는 실제 비행 데이터와 높은 일치도를 보였다. 특히, 롤 각도와 위치 오차가 각각 0.5° 이하, 2 m 이하로 제한된 점은, 역 시뮬레이션이 실시간 제어 설계에 적용 가능함을 시사한다.

이 논문의 의의는 두 가지로 요약할 수 있다. 첫째, 목표 기반 제어 설계 프레임워크를 제공함으로써, 제한된 연료·전력·통신 환경에서 사전 계산된 제어 프로파일을 활용해 자율 비행을 구현할 수 있다. 둘째, 기호 연산·고정밀 적분·유한차분을 결합한 수치 해법이 복잡한 비선형 비행역학 문제를 실용적인 계산 비용으로 해결한다는 점이다. 향후 연구에서는 다중 목표(예: 최소 연료 소모·최소 노출) 최적화와 실시간 재계산을 위한 GPU 가속, 그리고 비행 제어 시스템에의 통합 검증이 필요할 것으로 보인다.