19.20 회전 날개와 로터 공기역학
회전 날개(rotary wing)는 중심축을 기준으로 회전 운동을 수행하며 공기에 대한 상대 속도를 생성함으로써 양력과 추력을 발생시키는 양력면(lifting surface)이다. 헬리콥터, 멀티로터 무인 항공기(multirotor UAV), 틸트로터(tiltrotor) 등 회전익 항공 로봇 플랫폼의 공력 성능은 로터(rotor)의 유동장 해석을 기반으로 결정되므로, 본 절에서는 회전 날개의 공기역학적 원리를 운동량 이론, 블레이드 요소 이론, 와류 이론의 순으로 체계적으로 기술한다.
1. 회전 날개의 정의와 공력 좌표계
회전 날개는 고정 날개(fixed wing)와 달리 자유류(freestream)가 존재하지 않더라도 회전에 의한 상대 유속을 스스로 생성한다. 로터 블레이드의 임의 단면에서의 상대 속도 성분은 회전 속도에 의한 접선 성분 U_T = \Omega r과, 수직 유입 속도(induced inflow velocity) 및 축 방향 자유류의 합으로 이루어지는 수직 성분 U_P로 구성된다. 여기서 \Omega는 로터의 각속도이며 r은 회전축으로부터의 반경 방향 거리이다. 이때 블레이드 단면에서의 유효 받음각(effective angle of attack)은 다음과 같이 정의된다.
\alpha = \theta - \phi, \qquad \phi = \tan^{-1}\!\left(\frac{U_P}{U_T}\right)
\theta는 기하학적 피치 각(pitch angle)이며 \phi는 유입 각(inflow angle)이다. 로터 디스크(rotor disk)는 블레이드 끝단이 그리는 원판 형상의 가상 면을 의미하며, 로터 공기역학 해석의 기준 좌표계로 사용된다.
19.20.2 운동량 이론에 의한 호버링 해석
운동량 이론(momentum theory)은 Rankine과 Froude에 의해 선박 프로펠러 해석에 처음 적용되었으며, 이후 Glauert에 의해 회전 날개로 확장되었다. 운동량 이론은 로터를 무한히 많은 블레이드로 구성된 등가의 작동 원판(actuator disk)으로 이상화하고, 디스크를 통과하는 공기에 대해 질량, 운동량, 에너지 보존 법칙을 적용한다.
호버링(hovering) 상태에서 디스크 면적을 A, 공기 밀도를 \rho, 디스크에서의 유도 속도를 v_h라 하면, 디스크를 통과하는 질량 유량은 \dot m = \rho A v_h로 주어진다. 운동량 변화에 의한 추력은
T = \dot m \, w = 2 \rho A v_h^2
이 되며, 여기서 후류 최종 속도 w가 2 v_h임은 베르누이 방정식을 디스크 상류와 하류에 개별적으로 적용하여 유도된다. 위 식을 유도 속도에 대해 정리하면
v_h = \sqrt{\dfrac{T}{2 \rho A}}
를 얻는다. 호버링 시 로터가 공기에 전달하는 이상 유도 동력(ideal induced power)은 P_i = T v_h = T^{3/2} / \sqrt{2 \rho A}로 나타난다. 이 관계식은 동일 추력 조건에서 디스크 면적이 클수록 유도 동력이 감소함을 보이며, 대형 헬리콥터 로터의 공력 효율이 소형 멀티로터보다 높게 나타나는 근본적 이유가 된다.
19.20.3 전진 비행에서의 운동량 이론
로터가 자유류 속도 V_\infty로 전진하는 경우, 디스크 면에 경사각 \alpha_D가 존재한다고 가정하면 유도 속도 v_i는 Glauert의 운동량 관계식
T = 2 \rho A v_i \sqrt{(V_\infty \cos \alpha_D)^2 + (V_\infty \sin \alpha_D + v_i)^2}
으로부터 결정된다. 위 식은 비선형 대수 방정식이므로 일반적으로 반복적 수치 해법을 통해 v_i를 구한다. 전진 비율(advance ratio) \mu = V_\infty \cos \alpha_D / (\Omega R)이 증가함에 따라 유도 속도는 감소하며, 로터의 유도 동력 손실 역시 저감된다. 여기서 R은 로터 반경이다.
2. 로터 추력 계수와 동력 계수
로터의 공력 성능은 무차원화된 추력 계수와 동력 계수를 통해 기술된다. 표준 정의는 다음과 같다.
C_T = \dfrac{T}{\rho A (\Omega R)^2}, \qquad C_P = \dfrac{P}{\rho A (\Omega R)^3}
호버링의 경우 이상 유도 동력 계수는 C_{P_i} = C_T^{3/2} / \sqrt{2}로 표현되며, 실제 로터에서는 블레이드 형상의 비균일성, 팁 손실, 점성 항력 등으로 인해 이보다 큰 동력이 요구된다. 이러한 손실을 반영하기 위해 호버링 품질 지수(figure of merit) FM이 도입되며, 다음과 같이 정의된다.
FM = \dfrac{C_T^{3/2} / \sqrt{2}}{C_P}
잘 설계된 헬리콥터 메인 로터의 FM은 약 0.75 내외의 값을 가지는 것으로 알려져 있으며, 이는 이상 운동량 이론이 예측하는 최대 효율 FM = 1에 대한 실제 성능의 비율을 의미한다.
3. 블레이드 요소 이론
블레이드 요소 이론(Blade Element Theory, BET)은 Drzewiecki에 의해 제안된 해석 기법으로, 블레이드를 반경 방향으로 미소 구간으로 분할하여 각 요소에 2차원 익형 공기역학을 국소적으로 적용한다. 반경 r 위치의 폭 dr인 블레이드 요소에 작용하는 미소 양력과 항력은
dL = \tfrac{1}{2} \rho U^2 c \, C_l \, dr, \qquad dD = \tfrac{1}{2} \rho U^2 c \, C_d \, dr
으로 표현되며, 여기서 c는 단면의 코드 길이, C_l과 C_d는 각 익형 단면의 양력 계수와 항력 계수, U = \sqrt{U_T^2 + U_P^2}는 국소 상대 속도의 크기이다. 미소 추력과 미소 토크는 양력과 항력을 디스크 수직 방향 및 접선 방향으로 분해하여
dT = (dL \cos \phi - dD \sin \phi), \qquad dQ = (dL \sin \phi + dD \cos \phi) \, r
로 주어진다. 이를 블레이드 개수 N_b에 걸쳐 합산하고 반경 방향으로 적분함으로써 전체 로터 추력과 토크를 산출한다. 블레이드 요소 이론은 블레이드의 테이퍼, 비틀림, 익형 변화 등 기하학적 세부 사항을 해석에 반영할 수 있다는 장점을 가진다.
4. 결합 블레이드 요소-운동량 이론
블레이드 요소 이론은 유도 속도를 외부에서 부여해야 해석이 가능하다는 한계를 지닌다. 이를 극복하기 위해 블레이드 요소 이론과 운동량 이론을 결합한 BEMT(Blade Element Momentum Theory)가 사용된다. 반경 r, 폭 dr인 환형 고리(annular ring)에 대해 국소 운동량 보존을 적용하면
dT = 4 \pi \rho \, v_i^2 \, r \, dr
가 성립한다. 이를 블레이드 요소로 얻은 dT와 등치하면 각 반경에서의 유도 속도 v_i(r)에 대한 대수 방정식이 유도되며, 이로부터 반경 방향 유입 분포를 구할 수 있다. BEMT는 호버링과 축 상승 비행에서 신뢰할 만한 정확도를 제공하며, 소형 무인 로터 설계와 성능 예측에 널리 활용된다.
19.20.7 와류 이론과 자유 후류 해석
로터는 회전에 따라 블레이드 팁과 루트에서 강한 와류(vortex)를 방출하며, 이들 와류는 나선형 후류(helical wake) 구조를 형성한다. 와류 이론(vortex theory)은 Biot-Savart 법칙에 근거하여 후류 와류가 블레이드에 유도하는 속도를 계산함으로써 유입 분포를 해석적으로 구한다. 후류를 회전축과 평행한 강체 나선으로 가정하는 경직 후류(prescribed wake) 모델과, 후류 요소가 자체 유도 속도에 의해 자유롭게 변형되도록 허용하는 자유 후류(free wake) 모델이 존재한다. 자유 후류 해석은 전진 비행, 와류 고리 상태(vortex ring state) 등 복잡한 비정상 유동 현상을 포착하는 데 적합하다.
19.20.8 팁 손실과 비균일 유입
실제 로터에서는 블레이드 팁 근방에서 하면 고압 영역의 공기가 상면 저압 영역으로 누설되어 팁 와류를 형성하며, 이는 유효 블레이드 길이를 단축시키는 결과를 초래한다. Prandtl은 이러한 현상을 근사적으로 반영하는 팁 손실 계수(tip loss factor) B를 도입하였으며, 유효 반경은 BR로 정의된다. 일반적으로 B는 0.96 내외의 값을 가진다. 또한 실제 로터 디스크 면에서의 유입 속도는 균일하지 않으며, 블레이드 팁으로 갈수록 증가하는 경향을 보인다. 비균일 유입 분포는 동력 손실을 증가시키고 블레이드에 주기적 하중을 유발한다.
19.20.9 전진 비행의 비대칭 유동 환경
로터가 전진 비행을 수행할 때, 전진 측 블레이드(advancing blade)는 회전 속도에 자유류 속도가 더해져 상대 속도가 증가하는 반면, 후퇴 측 블레이드(retreating blade)는 상대 속도가 감소한다. 이로 인해 좌우 양력 비대칭이 발생하며, 후퇴 측에서는 역류 영역(reverse flow region)이 형성된다. 역류 영역은 반경 r < \mu R \, |\sin \psi|를 만족하는 원형 구역으로 나타나며, \psi는 방위각이다. 전진 측 블레이드 팁에서는 마하 수 증가에 따른 압축성 효과와 충격파 형성이, 후퇴 측에서는 받음각 증가로 인한 실속(stall)이 성능 한계를 결정짓는 주요 요인이 된다.
19.20.10 블레이드 플래핑과 힌지 운동
힌지 장착(articulated) 로터에서는 비대칭 양력에 의한 모멘트를 해소하기 위해 블레이드가 플래핑(flapping) 힌지를 중심으로 상하 운동을 수행한다. 플래핑 운동은 일반적으로 푸리에 급수
\beta(\psi) = \beta_0 - \beta_{1c} \cos \psi - \beta_{1s} \sin \psi - \cdots
로 표현되며, \beta_0은 코닝 각(coning angle), \beta_{1c}와 \beta_{1s}는 각각 세로 및 가로 플래핑 성분이다. 플래핑 운동은 블레이드 단면의 유효 받음각을 방위각에 따라 조정함으로써 좌우 양력 균형을 자동적으로 회복시키는 역할을 수행한다. Lock 수(Lock number) \gamma = \rho a c R^4 / I_\beta는 공력 모멘트와 관성 모멘트의 비를 나타내며, 플래핑 운동의 동적 응답 특성을 결정짓는 핵심 무차원 수이다. 여기서 a는 양력 곡선 기울기, I_\beta는 플래핑 축에 대한 블레이드의 질량 관성 모멘트이다.
5. 특이 유동 상태
수직 강하 비행 시 로터는 일련의 특이 유동 상태를 거친다. 정상 작동 상태(normal working state)에서는 후류가 디스크 하방으로 원활하게 배출된다. 강하 속도가 증가하여 유도 속도와 상쇄되기 시작하면 후류가 디스크 주변에 재순환되는 와류 고리 상태(vortex ring state)가 발생하며, 이 영역에서는 추력의 급격한 변동과 진동이 동반된다. 이어 난류 후류 상태(turbulent wake state)를 거쳐 최종적으로 풍차 브레이크 상태(windmill brake state)에 도달하며, 이 상태에서는 로터가 공기로부터 에너지를 흡수하여 자동 회전(autorotation)이 성립된다. 자동 회전은 동력 상실 시 헬리콥터의 비상 착륙을 가능하게 하는 핵심 원리이다.
6. 로봇 공학에서의 응용과 의의
회전 날개 공기역학의 이론 체계는 멀티로터 무인 항공기의 설계와 제어 전반에 직접적으로 적용된다. 소형 로터의 추력과 토크 상수는 BEMT와 실험적 캘리브레이션을 결합하여 산출되며, 이는 멀티로터의 동역학 모델에 입력 파라미터로 사용된다. 로터 간 공력 간섭, 지면 효과(ground effect), 벽면 근접 효과 등은 저고도 자율 비행의 제어 성능에 직접 영향을 미치므로 공기역학적 모델링이 필수적이다. 또한 헬리콥터형 무인기의 경우 플래핑 동역학과 유입 동역학이 비행 제어기 설계에 반영되며, 틸트로터 및 가변 형상 로터 플랫폼의 해석에도 본 절에서 기술한 이론들이 기반을 제공한다.
7. 출처
- Leishman, J. G., Principles of Helicopter Aerodynamics, 2nd ed., Cambridge University Press, 2006.
- Johnson, W., Helicopter Theory, Dover Publications, 1994.
- Bramwell, A. R. S., Done, G., and Balmford, D., Bramwell’s Helicopter Dynamics, 2nd ed., Butterworth-Heinemann, 2001.
- Stepniewski, W. Z., and Keys, C. N., Rotary-Wing Aerodynamics, Dover Publications, 1984.
- Seddon, J., and Newman, S., Basic Helicopter Aerodynamics, 3rd ed., Wiley, 2011.
- Prouty, R. W., Helicopter Performance, Stability, and Control, Krieger Publishing, 1995.
8. 버전
v1.0