항공탄성학

1. 개요[편집]

항공탄성학(aeroelasticity, 공탄성학)은 탄성 구조물이 유동 속에서 변형될 때, 그 변형이 다시 공기력을 바꾸고 바뀐 공기력이 또 변형을 바꾸는 되먹임 고리를 다루는 학문이다. 유체와 구조가 서로를 밀고 당기는 문제라는 점에서 유체-구조 연성의 한 갈래이지만, 항공탄성학은 그중에서도 “이 되먹임이 발산하는가, 수렴하는가”라는 안정성 판정에 초점이 있다는 점에서 특유의 성격을 갖는다.

날개를 무한히 단단하게 만들면 이 문제는 사라진다. 그런데 그러면 비행기가 못 뜬다. 구조를 가볍게 만들수록 성능이 좋아지고, 가벼워질수록 유연해지고, 유연해질수록 공탄성이 물어뜯는다 — 항공기 설계가 영원히 이 삼각관계에 붙잡혀 있는 이유다.¹ 이 문서는 항공탄성학 전체 지도를 그리는 허브이며, 개별 현상의 상세는 플러터, 정적 발산, 버피팅 문서에 위임한다.

2. 콜라의 공탄성 삼각형[편집]

1946년 A. R. 콜라가 제안한 삼각형은 이 분야의 지도로 지금도 그대로 쓰인다. 꼭짓점 셋은 공기력(A), 탄성력(E), 관성력(I) 이고, 각 변과 내부가 서로 다른 문제 영역이 된다.

조합	영역	대표 현상
A + E	정적 공탄성	발산, 조종면 반전, 하중 재분포
A + E + I	동적 공탄성	플러터, 버피팅, 돌풍 응답
E + I	구조 동역학	고유진동(모드 해석)
A + I	비행 역학	강체 안정성, 조종 응답

즉 관성이 빠지면 정적 문제, 세 힘이 다 들어오면 동적 문제다. 여기에 제어계까지 끼어들면 항공서보탄성학(aeroservoelasticity)이 되는데, 플라이-바이-와이어 항공기에서는 오히려 이쪽이 표준 상황이다. 제어기가 구조 모드를 잘못 밟으면 안정화 장치가 오히려 진동을 키우는 사태가 벌어진다.

3. 정적 공탄성[편집]

관성을 무시할 수 있을 만큼 천천히 벌어지는 문제들. 지배 구도는 “공기력 모멘트가 구조 강성보다 커지는가”라는 고유값 문제 하나로 요약된다.

• 정적 발산(divergence): 날개가 비틀리면 받음각이 커지고, 받음각이 커지면 양력과 비틀림 모멘트가 더 커진다. 동압이 어떤 임계값 $q_{div}$ 를 넘으면 이 되먹임이 강성을 이겨 변형이 무한대로 발산한다. 후퇴익은 굽힘이 앞전을 내려주는 방향이라 발산에 유리하고, 전진익은 반대라 불리하다 — X-29 같은 전진익기가 복합재 굽힘-비틀림 커플링(aeroelastic tailoring)을 반드시 필요로 했던 이유다(복합재 해석).
• 조종면 반전(control reversal): 에일러론을 아래로 꺾으면 그 부분 양력은 커지지만 동시에 날개를 앞전 아래로 비트는 모멘트가 생긴다. 고속에서는 이 비틀림이 만드는 양력 손실이 에일러론 효과를 역전시켜, 롤을 오른쪽으로 주려던 조작이 왼쪽 롤을 만든다. 실기에서 실제로 발생했던 사고 모드이며, 현대 여객기가 고속에서 외측 에일러론을 잠그고 내측만 쓰는 이유이기도 하다.
• 하중 재분포: 발산까지 가지 않아도, 유연한 날개는 변형하면서 양력 분포를 뿌리 쪽으로 옮긴다. 강체 가정으로 계산한 굽힘 모멘트와 실제 값이 다르다는 뜻이라, 피로 해석과 하중 설계에서 반드시 보정한다.

4. 동적 공탄성[편집]

관성이 들어오면 안정성이 아니라 진동의 문제가 된다.

• 플러터: 굽힘 모드와 비틀림 모드가 공기력을 매개로 결합해, 유동에서 에너지를 뽑아내는 자려진동이 되는 현상. 유속이 오르면 두 모드의 진동수가 접근하다가 어느 지점에서 한쪽 모드의 감쇠가 음으로 넘어간다. 그 속도가 플러터 속도 $V_f$ 이며, 넘으면 진폭이 지수적으로 커져 수 초 안에 구조가 파괴된다. 항공탄성학에서 가장 무서운 현상이자 가장 많은 연구가 붙은 주제.
• 버피팅: 박리 난류나 앞쪽 구조의 후류가 만드는 강제 진동. 플러터와 달리 자려진동이 아니라 외부 랜덤 하중에 대한 응답이므로, 안정성이 아니라 주파수 응답 해석과 확률적 피로 문제로 다룬다.
• LCO(limit cycle oscillation): 선형 이론상 불안정한데 진폭이 무한히 커지지 않고 일정 크기의 진동으로 정착하는 경우. 자유 유격(freeplay) 같은 구조 비선형이나 천음속 충격파 운동의 비선형성이 원인이며, 선형 플러터 해석으로는 원리적으로 예측되지 않는다(비선형 구조해석).
• 돌풍·착륙 응답: 대기 난류와 이산 돌풍에 대한 동적 응답. 인증 하중 케이스의 핵심 항목이다.
• 위르 플러터(whirl flutter): 프로펠러·틸트로터의 자이로스코픽 효과가 나셀 지지 강성과 결합해 일으키는 불안정. 회전체 동역학과 겹치는 영역.

5. 해석 방법[편집]

플러터 해석의 표준 형태는 다음 행렬 방정식의 고유값 문제다.

$$\left[-\omega^2 \mathbf{M} + i\omega \mathbf{C} + \mathbf{K} - q\,\mathbf{Q}(k, M_\infty)\right]\mathbf{q} = \mathbf{0}$$

$\mathbf{Q}$ 는 감소진동수 $k = \omega b / V$ 와 마하수에 의존하는 비정상 공기력 영향계수 행렬이고, 자유도는 보통 모드 해석으로 얻은 구조 고유모드 수십 개로 축소한다(축소차수모델).

• DLM(Doublet Lattice Method, 이중격자법): 1969년 알바노와 로든이 정리한 아음속 진동 양력면 이론. 양력면을 패널로 나누고 각 패널의 1/4 시위선에 가속도 퍼텐셜 더블렛 선을 놓아 $\mathbf{Q}$ 를 만든다. 반세기가 지난 지금도 NASTRAN SOL 145를 비롯한 인증용 해석의 주력이다. 빠르고 검증 이력이 두껍지만, 얇은 양력면·소진폭·아음속 선형이라는 가정이 깔려 있다.
• p-k 법: 위 방정식을 각 유속마다 반복적으로 풀어 감쇠와 진동수를 유속의 함수로 뽑는 표준 절차. 결과를 $V\!-\!g$, $V\!-\!f$ 선도로 그려 감쇠가 0을 가로지르는 지점을 플러터 속도로 읽는다. 유리함수 근사(RFA)로 $\mathbf{Q}$ 를 시간영역 상태공간으로 바꾸면 제어계와 함께 푸는 항공서보탄성 해석이 가능해진다.
• 고충실도 CFD-CSD 연성: 천음속 영역에서는 충격파의 위치가 진동에 따라 움직이면서 플러터 속도가 급격히 떨어지는 천음속 딥(transonic dip) 이 나타나는데, 선형 패널법으로는 잡히지 않는다. 이때는 압축성 유동 솔버와 구조 솔버를 유체-구조 연성으로 묶고 ALE 기법으로 격자를 움직이며 시간적분한다. 비싼 대신 LCO와 박리 효과까지 포함된다.

6. 시험·인증과 응용 분야[편집]

해석만으로 끝나는 분야가 아니다. 감항 규정(FAA/EASA의 25.629 등)은 설계 급강하 속도 $V_D$ 의 1.15배까지 플러터가 없음을 입증하도록 요구한다.

• 지상진동시험(GVT): 완성된 기체를 탄성 지지 위에 올려 놓고 가진해 실제 고유진동수·모드형상·감쇠를 측정, 유한요소 모델을 갱신(model updating)한다. 해석 모델의 신뢰도가 여기서 결정된다(검증 및 확인).
• 플러터 비행시험: 속도를 조금씩 올리면서 각 속도에서 구조를 가진하고 감쇠를 측정, 감쇠-속도 곡선을 외삽해 안전 여유를 확인한다. 가장 조심스럽게 진행되는 시험 중 하나다.²

응용은 항공기 밖으로도 넓다. 가스터빈·압축기 블레이드의 플러터와 강제응답은 터보기계 해석의 핵심 설계 제약이고, 풍력 터빈 블레이드는 길이가 100 m를 넘어가면서 공탄성이 지배 설계 인자가 됐다. 장대교량의 비틀림 플러터는 1940년 타코마 내로스 교 붕괴 이후 교량 설계의 필수 검토 항목이 되었으며³, 최근에는 초고층 건물의 와류유기진동, 심지어 접이식 드론의 유연 날개까지 같은 이론 틀로 다룬다.

7. 관련 문서[편집]

• 플러터 · 정적 발산 · 버피팅
• 유체-구조 연성 · 다물리 연성해석 · ALE 기법
• 모드 해석 · 주파수 응답 해석 · 감쇠 · 축소차수모델
• 공력해석 · 압축성 유동 · 충격파 · 난류
• 복합재 해석 · 피로 해석 · 비선형 구조해석
• NASTRAN · 회전체 동역학 · 터보기계 해석

8. Footnotes[편집]

1. 그래서 항공 구조 엔지니어에게 “무게 좀 줄여 주세요”는 “플러터 마진 좀 깎아 주세요”와 동의어로 들린다. 반대로 공탄성 담당자가 강성을 요구하면 중량 담당자의 표정이 굳는다. 회의실에서 가장 자주 붙는 두 부서다. ↩

2. 감쇠 곡선을 외삽한다는 게 무서운 점이다. 아직 가 보지 않은 속도에서의 안정성을, 가 본 속도들의 추세로 추정해야 한다. 게다가 플러터 직전에는 감쇠가 급격히 무너지는 경우가 많아 외삽이 낙관적으로 나오기 쉽다. ↩

3. 타코마 내로스 교의 붕괴 원인을 “카르만 와열과의 공진”으로 설명하는 교과서가 아직도 많은데, 현재 정설은 비틀림 자려진동(공탄성 플러터)이다. 와류 방출 진동수와 비틀림 고유진동수가 맞지 않았다는 점이 결정적 반증. 이 오해는 유체역학 교과서 오류 목록의 부동의 1위다. ↩