22.44 양력과 항력의 비행 시뮬레이션 모델링

1. 비행 시뮬레이션의 목적과 구성

비행 시뮬레이션(flight simulation)은 비행체의 운동 방정식, 공력 모형, 추진 모형, 환경 모형을 결합하여 실제 비행 상태를 수치적으로 재현하는 체계이다. 양력과 항력은 공력 모형의 핵심 출력이며, 기체 6자유도(6-DOF) 운동 방정식의 외력으로 입력된다. 시뮬레이션 루프는 일반적으로 상태 추정→공력·추진력 계산→운동 방정식 적분→환경 상호 작용→계기·센서 모형의 순서로 구성되며, 시간 적분기는 Runge-Kutta 4차 또는 가변 스텝 적분법을 주로 사용한다. 이 표준 구조는 Stevens, Lewis, Johnson이 Aircraft Control and Simulation(3rd ed., Wiley, 2015)에서 체계적으로 정리하였다.

2. 공력 모형의 수학적 구조

6자유도 모형에서 기체에 작용하는 공력 $\mathbf{F}_A$ 와 공력 모멘트 $\mathbf{M}_A$ 는 바람 축 또는 기체 축 좌표계에서 다음과 같이 표현된다.

$\mathbf{F}_A = \frac{1}{2} \rho V^2 S \begin{bmatrix} -C_D \\ C_Y \\ -C_L \end{bmatrix}, \quad \mathbf{M}_A = \frac{1}{2} \rho V^2 S \begin{bmatrix} b C_l \\ c C_m \\ b C_n \end{bmatrix}$

여기서 $V$ 는 상대 대기 속도, $S$ 는 기준 날개 면적, $b$ 는 날개폭, $c$ 는 평균 공력 시위이다. 각 계수는 상태 변수와 조종면 입력의 함수로 모형화된다. 양력과 항력의 비행 시뮬레이션 모델링은 곧 $C_L$ 과 $C_D$ 의 다변수 모형 구조 결정과 매개변수 식별 과정이다.

3. 테이블 기반 공력 모형

비행 시뮬레이션에서 가장 널리 사용되는 공력 모형은 다변수 룩업 테이블(lookup table) 형태이다. $C_L$ 과 $C_D$ 는 다음과 같은 일반 형식으로 정의된다.

$C_L = C_L(\alpha, \beta, M, Re, \delta_e, \delta_a, p, q, r, \dotsc)$

이때 테이블은 풍동 시험, CFD, 비행 시험 자료로부터 구성되며, 다차원 보간(multilinear, spline interpolation)으로 임의 상태에서의 값을 산출한다. 보간 성능과 정밀도의 균형을 위해 변수 축을 물리적으로 의미 있는 구간으로 구성하며, 경계 외(off-nominal) 영역에 대한 외삽 정책을 명확히 설정한다. 이 접근은 NASA의 F-16 모형, JSBSim, MATLAB/Simulink Aerospace Blockset이 공통적으로 채택한 표준 형식이다.

4. 해석적·준해석적 모형

설계 초기 단계 또는 제어기 합성 단계에서는 테이블 대신 해석적·준해석적 공력 모형이 사용된다. 대표적 형식은 선형 유도 계수 모형(stability derivative model)이다.

$\begin{aligned} C_L &\approx C_{L_0} + C_{L_\alpha} \alpha + C_{L_{\delta_e}} \delta_e \\ C_D &\approx C_{D_0} + k C_L^2 \end{aligned}$

여기서 $C_{L_0}$ 는 영 받음각 양력 계수, $C_{L_\alpha}$ 는 받음각 기울기, $C_{L_{\delta_e}}$ 는 승강키 효율, $C_{D_0}$ 는 영 양력 항력, $k$ 는 유도 항력 계수이다. 유도 항력 계수는 $k = 1 / (\pi e \mathrm{AR})$ 로 표현되며, $e$ 는 Oswald 효율이다. 이러한 선형 모형은 Etkin과 Reid의 Dynamics of Flight: Stability and Control(3rd ed., Wiley, 1996)에서 정식화된 비행 동역학 표준 표기에 부합하며, 선형화된 상태 공간 모형을 구성하는 기반이 된다.

5. 비정상 공력 효과의 모형화

급격한 받음각 변화, 조종면의 급속 편위, 돌풍(gust)이 존재하는 조건에서는 비정상 공력(unsteady aerodynamics)이 중요해진다. 대표적 모형은 다음과 같다. 첫째, Theodorsen 함수 기반의 선형 비정상 공력 모형은 조화(harmonic) 운동에 대한 2차원 익형의 양력 응답을 복소 전달 함수로 표현한다. 둘째, Wagner 함수는 단위 받음각 계단 입력에 대한 시간 영역 비정상 양력 응답을 기술한다. 셋째, Leishman-Beddoes 동적 실속 모형은 공격각이 정적 실속각을 초과하는 영역에서의 비정상 양력과 모멘트 응답을 상태 공간 형태로 모형화한다. Leishman이 Principles of Helicopter Aerodynamics(2nd ed., Cambridge University Press, 2006)에서 이러한 비정상 모형을 정리하였다.

6. 멀티로터 기체의 간이 공력 모형

쿼드로터 및 소형 멀티로터 시뮬레이션에서는 개별 로터의 추력을 $T_i = k_T \omega_i^2$ 로 단순화하고, 전진 비행 시의 유도 항력 성분을 $-k_D \mathbf{v}$ 의 선형 공력 저항 항으로 근사하는 경우가 많다. Mahony, Kumar, Corke가 Multirotor Aerial Vehicles: Modeling, Estimation, and Control of Quadrotor(IEEE Robotics and Automation Magazine, vol. 19, no. 3, 2012)에서 제시한 모형이 대표적 예시이다. 또한 기체 동체의 양력은 통상 무시되지만, 동체의 상대 유동에 의한 기생 항력과 프로펠러 간섭에 의한 추력 변동은 고속 비행 조건에서 반드시 포함된다.

7. 고정익 무인기 시뮬레이션의 구조

고정익 무인기에서는 날개, 미익, 동체가 분리된 공력 구성을 이루며, 각 구성 요소의 양력·항력을 합산하여 전체 공력을 산출한다. 전형적 분해식은 다음과 같다.

$\begin{aligned} C_L &= C_{L,\text{wing}} + C_{L,\text{tail}} + C_{L,\text{fuselage}} \\ C_D &= C_{D,0} + \frac{C_L^2}{\pi e \mathrm{AR}} + C_{D,\text{trim}} + C_{D,\text{other}} \end{aligned}$

여기서 $C_{D,\text{trim}}$ 은 수평 안정판의 트림 항력, $C_{D,\text{other}}$ 는 랜딩기어, 안테나, 센서 등의 부가 항력이다. 이러한 분해 모형은 JSBSim의 XML 기반 구성 파일 및 X-Plane의 Blade Element Theory 기반 내부 모형에서 실제 구현에 반영된다.

8. 환경 요소의 결합

비행 시뮬레이션의 공력 모형은 환경 모형과 결합된다. 주요 환경 요소는 대기 밀도 $\rho$ , 온도 $T$ , 음속 $a$ , 점성 계수 $\mu$ , 풍속·난류 프로파일 등이다. 국제 표준 대기(International Standard Atmosphere, ISA)가 기본 기준이며, 실내 비행 시뮬레이션은 실내 환경 상수 또는 국지 관측값을 사용한다. 난류 및 돌풍 모형으로는 Dryden 난류 모형, von Karman 난류 모형, Kaimal 스펙트럼 모형이 사용되며, MIL-F-8785C 및 MIL-HDBK-1797과 같은 군사 규격 지침에 따라 구현된다. 이러한 환경 모형이 양력과 항력에 주는 영향은 상대 대기 속도 $\mathbf{V}_{\text{rel}} = \mathbf{V}_{\text{aircraft}} - \mathbf{V}_{\text{wind}}$ 의 동압과 기하 받음각·측풍각 변화를 통해 반영된다.

9. 시뮬레이션의 실시간성과 계산 효율

실시간 비행 시뮬레이터(Hardware-in-the-Loop, HIL)는 공력 모형이 주어진 시간 스텝 내에 반드시 계산 완료되어야 한다. 이를 위해 다차원 테이블 보간, 선형 유도 계수 모형, 적은 수의 미분 방정식으로 구성된 비정상 모형이 선호된다. 이와 반대로 오프라인 분석용 시뮬레이션은 시간 지연을 허용하므로, 실시간 Large Eddy Simulation과 같은 고해상도 모형도 부분적으로 결합될 수 있다.

10. 시뮬레이션의 검증과 매개변수 식별

시뮬레이션의 공력 모형은 풍동 시험, 비행 시험, 정지 시험 자료를 이용하여 매개변수 식별(system identification)과 검증(validation)을 거친다. 자주 사용되는 식별 기법은 극대우도법(Maximum Likelihood Estimation, MLE), 최소자승법, 확장 칼만 필터, 부분 공간 식별(subspace identification) 등이다. Jategaonkar가 Flight Vehicle System Identification: A Time Domain Methodology(2nd ed., AIAA Progress in Astronautics and Aeronautics, 2015)에서 이러한 시간 영역 식별 방법을 상세히 기술하였다. 식별 결과는 신뢰 구간과 함께 문서화되어 모형의 적용 범위를 명시한다.

11. 로봇공학적 의의

자율 비행 로봇과 실내 비행 로봇의 연구 개발에서 시뮬레이션은 필수적이다. 공력 모형의 정확도는 제어기 설계의 강건성, 학습 기반 제어의 수렴성, 임무 계획의 현실성에 직접 영향을 미친다. Gazebo, RotorS, AirSim, Isaac Sim과 같은 로봇 시뮬레이터는 공력 모형을 모듈식으로 구성하며, 다양한 수준의 공력 모형을 교체 적용할 수 있도록 설계되어 있다. 이러한 환경에서 실측 자료를 반영한 고신뢰 공력 모형의 구축은 자율 비행 기술의 실환경 이식성(sim-to-real) 향상에 기여한다.

12. 출처

Stevens, B. L., Lewis, F. L., and Johnson, E. N. Aircraft Control and Simulation, 3rd ed. Wiley, 2015.
Etkin, B., and Reid, L. D. Dynamics of Flight: Stability and Control, 3rd ed. Wiley, 1996.
Leishman, J. G. Principles of Helicopter Aerodynamics, 2nd ed. Cambridge University Press, 2006.
Mahony, R., Kumar, V., and Corke, P. “Multirotor Aerial Vehicles: Modeling, Estimation, and Control of Quadrotor.” IEEE Robotics and Automation Magazine, vol. 19, no. 3, 2012.
Jategaonkar, R. V. Flight Vehicle System Identification: A Time Domain Methodology, 2nd ed. AIAA Progress in Astronautics and Aeronautics, vol. 245, 2015.
U.S. Department of Defense. Flying Qualities of Piloted Aircraft. MIL-HDBK-1797, 1997.
U.S. Air Force. Military Specification: Flying Qualities of Piloted Airplanes. MIL-F-8785C, 1980.

13. 버전

v1.0 (2026-04-17)