26.19 이륙 과정에서의 지면 효과 천이 해석

1. 이륙 과정에서의 지면 효과 천이 개요

이륙 과정은 기체가 정지 혹은 저속 호버링 상태에서 자유류 고도(out-of-ground-effect, OGE)로 이행하는 시간 의존적 비행 단계이다. 이 단계에서 기체는 수직 속도 $w(t)$ 와 수평 속도 $u(t)$ 를 동시에 증가시키며 지면과의 상대 고도 $h(t)$ 가 급격히 변화한다. 지면 효과(in-ground-effect, IGE)는 유효 고도와 로터 또는 익형의 기하학적 치수비에 의해 결정되므로, 이륙 중의 공력 계수는 준정상(quasi-steady) 가정으로 단순 해석하기 어려운 동적 천이(transient) 특성을 보인다. 본 절에서는 이 천이 구간에서 나타나는 유도 속도의 재분배, 후류 축소, 지면 와류 붕괴 현상을 정량적으로 서술한다.

2. 천이 구간의 정의와 비차원 파라미터

지면 효과 천이 해석에서 중심적으로 사용되는 비차원 고도는 로터 반경 $R$ 에 대한 상대 고도 $z/R$ 또는 주익 반폭 $b/2$ 대비 $h/b$ 이다. Cheeseman과 Bennett의 “The Effect of the Ground on a Helicopter Rotor in Forward Flight”(Aeronautical Research Council Reports and Memoranda, 1955)에 의하면 로터 추력 증가량 $T/T_\infty$ 는 $z/R$ 의 역함수 형태로 근사된다. 이륙 중 $z/R$ 이 1 미만에서 2를 초과하는 영역으로 변화할 때 해당 식은 시간 의존 입력을 갖는 비선형 상미분 방정식이 되며, 공력 응답 지연을 갖는 1차 또는 2차 동적 모델로 표현할 필요가 있다.

또 다른 핵심 파라미터는 유도 속도 시간 상수 $\tau_w$ 이다. 이는 로터 원판 하부의 유도 유동장이 새로운 평형 상태에 도달하는 데 소요되는 시간 척도이며, 일반적으로 $\tau_w \approx R/v_h$ 로 근사된다. 여기서 $v_h$ 는 호버 유도 속도이다. 이륙 중 기체의 상승률이 $\tau_w$ 의 역수보다 충분히 작을 때는 준정상 해석이 유효하지만, 드론의 급가속 이륙과 같이 상승률이 크면 비정상(unsteady) 항을 포함해야 한다.

3. 호버 이륙에서의 추력 잉여와 가속 동역학

수직 이륙 상태에서 로터는 IGE 조건으로 인해 동일 집합 피치 입력에 대해 추가 추력을 발생한다. 이를 추력 잉여(thrust surplus) $\Delta T(h) = T_{IGE}(h) - T_{OGE}$ 로 정의한다. 기체 질량 $m$ 의 수직 가속도는 다음과 같이 표현된다.

$m\ddot{h} = T(h, \theta_{col}) - mg$

여기서 $\theta_{col}$ 은 집합 피치 지령이다. 지면 근처에서는 동일 $\theta_{col}$ 에 대해 $T$ 가 증가하므로 기체가 비선형적으로 가속된다. 이륙 초기에 집합 피치가 호버 트림 값으로 설정되어 있으면 기체는 예상보다 큰 상승률로 이륙하게 되며, 상승률이 증가함에 따라 $z/R$ 이 증가하고 지면 효과가 감소한다. 이 과정이 진행됨에 따라 추력은 OGE 값으로 점근하고 상승률은 안정화된다. 제어기는 이 추력 천이를 감지하고 집합 피치를 적시에 감소시켜 과도한 상승을 방지해야 한다.

로터 후류 구조의 천이

Curtiss, Sun, Putman, and Hanker의 “Rotor Aerodynamics in Ground Effect at Low Advance Ratios”(Journal of the American Helicopter Society, 1984)에 따르면 로터가 지면 근처에서 이륙하는 과정에서 후류는 세 가지 뚜렷한 단계를 거친다. 첫째, 낮은 고도에서 후류는 지면에 충돌한 후 반경 방향 분사(radial jet) 형태로 전개된다. 둘째, 고도가 증가함에 따라 후류는 수축(contraction)하며 지면 충돌 강도가 감소한다. 셋째, $z/R > 2$ 근방에서 후류는 완전한 원추형 수축을 회복하며 지면의 영향은 측정 범위 이하가 된다. 이러한 후류 천이는 로터 원판에 유도되는 속도 분포를 시간 의존적으로 변화시키며 결과적으로 로터 추력의 시간 변동을 야기한다.

고정익 이륙 활주에서의 지면 효과 천이

고정익 비행체의 이륙 활주는 공력 양력 계수의 점진적 변화 과정으로 해석한다. 정지 상태에서 회전 속도(rotation speed) $V_R$ 에 도달하기까지 기체는 지면에 접촉한 채 가속하며, 받음각은 지상 접촉 각도로 고정된다. $V_R$ 에서 조종사 또는 자동 비행 제어기가 피치업 명령을 입력하면 양력이 중력을 초과하여 기체가 지면을 벗어난다. 이후 상승각 $\gamma$ 로 상승하는 동안 $h/b$ 가 급격히 증가한다. 유도 항력 계수의 지면 효과 보정 인자 $\sigma$ 는 Wieselsberger의 고전식

$\sigma(h/b) = \frac{1 - 1.32(h/b)}{1.05 + 7.4(h/b)}$

로 근사되며, 0 < $h/b$ < 1.5 범위에서 널리 인용된다. 이륙 직후 유도 항력이 급격히 증가하는 천이 구간에서는 추력 잉여가 감소하므로 상승률 손실이 발생할 수 있다. 이 현상은 지면 효과 감소에 의한 초기 상승 정체로 알려져 있으며, 단거리 이륙 활주 성능 해석에서 주요 고려 사항이다.

4. 멀티로터의 이륙 천이 동역학

멀티로터 무인기는 정지 상태에서 스로틀이 증가함에 따라 네 개 또는 그 이상의 로터가 동시에 추력을 발생시킨다. 저고도에서 각 로터는 독립적인 IGE 상태에 놓이나, 로터 간격이 좁은 경우 로터 상호 간섭으로 비선형 결합이 발생한다. Sharf, Nahon, Harmat, Khan, Michini, Speal, Trentini, Tsadok, and Wang의 “Ground Effect Experiments and Model Validation with Draganflyer X8 Rotorcraft”(International Conference on Unmanned Aircraft Systems, 2014)는 쿼드로터와 옥토로터의 이륙 천이 중 추력 계수의 시간 응답을 실험적으로 측정하였다. 이 연구에 따르면 이륙 천이 시간은 로터 직경, 기체 질량, 상승 속도 지령에 의해 결정되며, 일반적으로 0.5초에서 2초 범위의 시간 상수를 갖는다.

5. 지면 와류의 생성과 소멸

이륙 초기 단계에서 후류가 지면에 충돌하면 반경 방향 분사 유동이 주변 정지 공기와 상호작용하여 도넛 형태의 지면 와류(ground vortex)를 형성한다. 이 와류는 로터 전진 방향에서 특히 강하게 형성되며, 로터 원판의 전방 영역으로 이류(advection)되어 양력 비대칭과 피칭 모멘트를 유발한다. 기체가 상승함에 따라 지면 와류의 강도는 감소하며, 임계 고도 $z/R \approx 1.5$ 근방에서 유의한 영향이 소멸한다. Light의 “Tip Vortex Geometry of a Hovering Helicopter Rotor in Ground Effect”(Journal of the American Helicopter Society, 1993)는 광학 측정을 통해 이 와류의 공간 구조를 규명하였다.

6. 천이 해석을 위한 수학적 모델링

지면 효과 천이를 해석하기 위한 전통적 접근은 다음 세 가지로 분류된다. 첫째, 준정상 공력 계수 내삽(quasi-steady coefficient interpolation) 방식은 고도별 정상 상태 계수를 시간에 따라 내삽한다. 둘째, Peters와 He의 유한상태 유도 유동 이론에 지면 효과 이미지 항을 추가하는 방식은 유도 속도의 과도 응답을 적절히 포착한다. 셋째, 격자 기반 전산유체역학 해석에서 이동 격자(dynamic mesh) 또는 중첩 격자(overset mesh) 기법을 사용하여 기체의 실시간 이동을 모사하는 방식은 가장 엄밀하나 계산 비용이 매우 크다. 무인기 제어 시스템 설계에서는 일반적으로 첫째 또는 둘째 접근법을 채택한다.

7. 제어기 설계 관점에서의 시사점

이륙 과정의 지면 효과 천이는 개방 루프(open-loop) 비행에서 추력 과잉 및 상승률 폭주를 초래할 수 있다. 폐쇄 루프 고도 제어기는 다음 세 가지 관점에서 보상이 요구된다. 먼저 피드포워드(feed-forward) 경로에 지면 효과 보상항을 추가하여 예측 가능한 추력 변화를 상쇄한다. 다음으로 고도 측정 기반 이득 스케줄링(gain scheduling)을 도입하여 IGE와 OGE 영역의 이득을 분리 설계한다. 마지막으로 모델 예측 제어(Model Predictive Control)에서는 지면 효과를 포함하는 비선형 공력 모델을 예측기로 사용하여 상승 궤적을 최적화한다. Powers, Mellinger, Kushleyev, Kothmann, and Kumar의 “Influence of Aerodynamics and Proximity Effects in Quadrotor Flight”(International Symposium on Experimental Robotics, 2012)는 이러한 제어기가 이륙 중 과도 응답을 유의미하게 감소시킴을 실증하였다.

8. 출처

Cheeseman, I. C., and Bennett, W. E., “The Effect of the Ground on a Helicopter Rotor in Forward Flight,” Aeronautical Research Council Reports and Memoranda No. 3021, 1955.
Curtiss, H. C., Sun, M., Putman, W. F., and Hanker, E. J., “Rotor Aerodynamics in Ground Effect at Low Advance Ratios,” Journal of the American Helicopter Society, Vol. 29, No. 1, 1984.
Light, J. S., “Tip Vortex Geometry of a Hovering Helicopter Rotor in Ground Effect,” Journal of the American Helicopter Society, Vol. 38, No. 2, 1993.
Sharf, I., Nahon, M., Harmat, A., Khan, W., Michini, M., Speal, N., Trentini, M., Tsadok, T., and Wang, T., “Ground Effect Experiments and Model Validation with Draganflyer X8 Rotorcraft,” International Conference on Unmanned Aircraft Systems, 2014.
Powers, C., Mellinger, D., Kushleyev, A., Kothmann, B., and Kumar, V., “Influence of Aerodynamics and Proximity Effects in Quadrotor Flight,” International Symposium on Experimental Robotics, 2012.
Leishman, J. G., “Principles of Helicopter Aerodynamics,” 2nd ed., Cambridge University Press, 2006.

9. 버전

문서 버전: v1.0
작성 기준일: 2026-04-17