플러터

시뮬레이션 구조해석 전산유체역학 유체역학 마지막 수정: 2026-07-19 04:36:05

1. 개요[편집]

플러터 Flutter
분야	항공탄성학 (FSI × 구조 동역학)
본질	비정상 공기력이 에너지를 공급하는 자기여기 진동
판정	복소 고유값의 감쇠가 0을 지나는 속도
대표 해법	K법, p-k법, 시간영역 CFD-CSD 연성

플러터(flutter)는 탄성 구조물의 진동이 그 진동 자체가 만들어내는 비정상 공기력으로부터 에너지를 공급받아, 유효 감쇠가 음(-)이 되어 진폭이 스스로 커지는 동적 항공탄성 불안정 현상이다. 핵심 단어는 자기여기(self-excited)다. 돌풍이나 엔진 같은 외부 가진이 전혀 없어도, 일단 살짝 흔들리기 시작하면 유동이 알아서 그 흔들림을 키워준다.

플러터가 무서운 이유는 두 가지다. 첫째, 임계 속도 아래에서는 아무 징조도 없다. 진동이 잘 잦아들던 구조가 속도를 조금 더 올린 순간 발산으로 돌아선다. 둘째, 발산이 지수적이다. 임계 속도를 넘긴 뒤 구조가 부서지기까지 수 초가 걸리는 경우도 흔하다.¹ 그래서 항공기 인증에서 플러터는 “안전율을 얼마 두느냐”가 아니라 “운용 영역 안에 절대 없어야 한다”는 이진 조건으로 다뤄진다.

항공탄성학의 고전적 지도는 콜라(A. R. Collar)의 항공탄성 삼각형이다. 공기력·탄성력·관성력 세 꼭짓점 중 두 개만 얽히면 정적 항공탄성(정적 발산, 조종면 반전)이나 순수 구조 진동이고, 세 개가 모두 얽힐 때 그 자리에 플러터가 있다.

2. 굽힘-비틀림 플러터[편집]

가장 고전적인 형태는 날개의 굽힘 모드와 비틀림 모드가 결합하는 2자유도 플러터다. 메커니즘은 이렇다.

날개가 위로 휘면서 동시에 앞전이 들리는 쪽으로 비틀리면, 받음각이 커져 양력이 증가하고 그 양력이 날개를 더 위로 밀어 올린다. 반대로 내려갈 때 앞전이 내려가면 아래로 미는 힘이 커진다. 즉 굽힘과 비틀림 사이의 위상차가 특정 범위에 들면, 공기력이 한 주기 동안 구조에 순(net) 일을 해준다. 이 공급 에너지가 구조 감쇠가 소산할 수 있는 양을 넘어서는 속도가 플러터 속도 $V_F$ 다.

주파수 관점에서 보면 속도가 오를수록 두 모드의 주파수가 서로 가까워지다가(coalescence) 한 가지의 감쇠가 0을 뚫고 내려간다. 그래서 플러터 대책의 정석은 두 모드를 떼어놓는 것이다. 비틀림 강성을 키우거나, 무게 중심을 앞으로 옮겨(질량 균형, mass balancing) 굽힘-비틀림 관성 결합을 줄인다. 조종면에 카운터웨이트를 다는 이유가 정확히 이것이다.²

3. 감소 진동수와 비정상 공기력[편집]

플러터 해석이 공력해석과 결정적으로 다른 점은 공기력이 준정상(quasi-steady)이 아니라는 것이다. 날개가 흔들리면 그 뒤로 후류 와류가 흘러가고, 이 와류가 다시 날개에 영향을 준다. 결과적으로 공기력은 현재 받음각뿐 아니라 그 변화의 이력에 의존하며, 변위에 대해 위상이 지연된다. 이 위상 지연이 바로 에너지를 주고받는 통로다.

지배 무차원수는 감소 진동수(reduced frequency)다.

k = \frac{\omega b}{U}

$b$ 는 반현(semi-chord), $U$ 는 자유류 속도다. $k$ 는 “공기 입자가 익현을 지나가는 동안 구조가 몇 번 흔들리는가”를 재는 값으로, $k \ll 1$ 이면 준정상 근사가 통하고 $k$ 가 커질수록 비정상 효과가 지배한다. 플러터는 대개 $k$ 가 0.1~1 정도인 영역에서 일어난다.

비압축·조화진동하는 얇은 익형에 대해서는 1935년 테오도르센(Theodorsen)이 공기력을 닫힌 형태로 풀었다. 그 결과에 등장하는 복소 함수 $C(k)$ (테오도르센 함수)가 후류의 위상 지연을 요약한다. 3차원 아음속 형상에는 더블릿 격자법(Doublet Lattice Method, DLM)이 산업 표준으로, 양력면을 패널로 나눠 주파수 영역에서 공기력 영향계수 행렬을 만든다. NASTRAN의 항공탄성 해석 체계가 이 위에 서 있다.

4. 플러터 방정식과 풀이법[편집]

구조를 모드 해석으로 축약한 뒤 일반화 좌표 $\mathbf{q}$ 로 쓰면, 플러터 방정식은 다음 꼴의 복소 고유값 문제가 된다.

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

$q_\infty$ 는 동압, $\mathbf{Q}$ 는 감소 진동수와 마하수에 의존하는 비정상 공기력 행렬이다. 골치 아픈 점은 $\mathbf{Q}$ 가 미지수인 $\omega$ 에 의존한다는 것 — 답을 알아야 행렬을 만들 수 있다. 그래서 반복이 필요하고, 역사적으로 두 갈래 방법이 자리 잡았다.

K법(V-g법) — 방정식이 조화 진동한다고 가정하고 인위적 구조 감쇠 $g$ 를 미지수로 넣어, 해가 존재하는 $g$ 를 찾는다. $g = 0$ 이 되는 속도가 플러터 속도다. 빠르고 견고하지만, 임계점 밖에서 나오는 감쇠 값은 물리적 의미가 없다.
p-k법 — 복소 고유값의 실수부를 실제 감쇠로 해석하며, 속도를 훑으면서 각 모드의 감쇠 추이를 그린다. 감쇠 곡선이 물리적으로 읽히기 때문에 비행 시험 데이터와 직접 비교할 수 있어 현업에서 선호된다.

결과는 언제나 V-g / V-f 선도(속도 대 감쇠, 속도 대 진동수)로 정리된다. 어느 가지의 감쇠 곡선이 0선을 아래로 뚫는 지점이 플러터 개시점이며, 그때의 모드 형상을 보면 어떤 모드끼리 결합했는지가 드러난다.

5. 천음속과 비선형 — CFD가 필요한 이유[편집]

DLM은 선형 포텐셜 이론이라 충격파를 모른다. 문제는 천음속 영역에서 플러터 경계가 뚝 떨어지는 현상(transonic dip)이 실제로 존재하고, 하필 그 골짜기가 여객기 순항 마하수 근처에 있다는 것이다. 원인은 익면 위 충격파가 구조 진동에 따라 앞뒤로 움직이며 공기력의 위상을 크게 바꾸는 데 있다. 선형 이론은 이 구간을 낙관적으로 예측하는데, 안전 쪽으로 틀리는 게 아니라 위험한 쪽으로 틀린다.

그래서 천음속 플러터는 압축성 유동 CFD와 구조 솔버를 시간 영역에서 묶는 양방향 유체-구조 연성 해석으로 간다. 구조를 모드 좌표로 축약하고 유체 격자를 ALE 기법으로 변형시키며 시간 적분한 뒤, 응답이 커지는지 잦아드는지로 안정성을 판정하는 방식이다. 다만 속도·마하수·고도 조합마다 해석을 다시 돌려야 해서 비용이 무섭다. 그래서 CFD로 소수의 조건만 계산하고 축소차수모델이나 대리 모델로 경계 전체를 메우는 절충이 널리 쓰인다.

비선형이 개입하면 진폭이 무한정 커지지 않고 일정 크기로 안착하는 리미트 사이클 진동(LCO)이 나타나기도 한다. 조종면 유격이나 충격파 운동이 원인이며, 당장 부서지지는 않지만 피로 수명을 갉아먹는다. 진폭에 의존하는 현상이라 선형 고유값 해석으로는 원리적으로 잡히지 않는다.

6. 플러터의 친척들[편집]

정적 발산(divergence) — 비틀림 강성이 공기력의 비틀림 모멘트에 지는 정적 불안정. 진동 없이 그냥 비틀려 부러진다. 관성이 개입하지 않는다는 점에서 플러터와 다르다.
휠 플러터(whirl flutter) — 프로펠러/로터의 자이로스코프 효과와 나셀 지지 강성이 결합한 불안정. 1959~60년 록히드 일렉트라의 연속 추락 원인으로 밝혀지며 유명해졌다.
버피팅(buffeting) — 박리 유동이나 후류의 난류가 구조를 때리는 강제 진동. 자기여기가 아니라는 점에서 플러터와 근본적으로 다르지만, 현장에서 자주 혼용된다.
터보기계 블레이드 플러터 — 압축기·팬 블레이드 열이 특정 노드 직경(nodal diameter) 패턴으로 함께 진동하며 불안정해진다. 터보기계 해석에서는 인접 블레이드 간 위상차(IBPA)를 훑으며 공력 감쇠를 계산하는 것이 표준 절차다.
교량 플러터 — 1940년 타코마 내로스 교량 붕괴는 단순 공진이 아니라 1자유도 비틀림 플러터에 가까운 자기여기 불안정으로 보는 것이 오늘날의 정설이다. 이후 장대교 설계에서 내풍 안정성 검토가 필수가 되었다.³

7. 인증과 실무[편집]

민항기는 감항 기준(FAR/CS 25.629)에 따라 급강하 속도 $V_D$ 의 1.15배까지 플러터가 없어야 한다. 이 여유는 해석 오차와 구조 특성 산포를 흡수하기 위한 것이다.

실무 절차는 대체로 이렇다. 유한요소 모델을 만들고 → 지상 진동 시험(GVT) 으로 실제 모드 주파수·형상을 측정해 모델을 갱신하고 → 공기력 모델과 결합해 플러터 해석을 돌리고 → 비행 플러터 시험에서 속도를 단계적으로 올리며 각 모드의 감쇠를 추출해 안전 여유를 확인한다. 비행 시험에서는 감쇠 추세를 외삽해 다음 속도 단계로 갈지 결정하는데, 감쇠가 예상보다 빠르게 줄면 그 자리에서 시험을 중단한다. 항공기 시험 중 가장 신경이 곤두서는 항목으로 꼽히는 이유다.

8. 관련 문서[편집]

플러터 시험 영상이 늘 짧은 이유가 있다. 임계점을 넘으면 진폭이 지수적으로 커지므로, 위험을 인지하고 속도를 줄이는 시간 자체가 없는 경우가 많다. “천천히 나빠지는” 여유를 주지 않는다는 점에서 좌굴이나 피로와 성격이 다르다. ↩
조종면 앞쪽에 붙은 정체불명의 쇳덩어리는 무게를 늘리려고 단 게 아니라, 조종면 무게중심을 힌지 축 앞으로 당겨 조종면 플러터를 막으려고 단 것이다. 항공기에서 일부러 무게를 더하는 몇 안 되는 사례이며, 그만큼 플러터가 무섭다는 방증이기도 하다. ↩
타코마 다리를 “공진으로 무너진 사례”로 소개하는 물리 교과서가 아직도 많은데, 바람은 다리의 고유 진동수로 주기 가진을 한 적이 없다. 구조의 비틀림 운동이 스스로 유동을 조직해 에너지를 받아낸 자기여기 현상이라는 것이 항공탄성학계의 공통된 설명이다. 유명한 사례일수록 잘못된 설명이 오래 산다. ↩

플러터