항공서보탄성학

1. 개요[편집]

항공서보탄성학(aeroservoelasticity, ASE)은 공기력·구조 탄성·비행제어계 셋이 서로를 밀고 당기는 3자 연성 문제를 다루는 분야다. 항공탄성학의 콜라 삼각형이 공기력(A)·탄성력(E)·관성력(I)의 세 꼭짓점이었다면, 여기에 제어(S)가 네 번째 꼭짓점으로 들어와 사면체가 된 것이라고 보면 된다.

핵심은 이거다. 자이로와 가속도계는 강체 운동과 탄성 진동을 구분하지 못한다. 센서가 붙어 있는 그 지점의 운동을 그냥 다 읽을 뿐이다. 그래서 안정성 증강을 위해 설계된 제어 법칙이 기체의 굽힘 모드를 강체 자세 변화로 착각하고 반응하면, 조종면이 그 모드의 진동수로 흔들리고, 그 조종면이 만든 비정상 공기력이 다시 같은 모드를 가진한다. 되먹임 고리가 닫히는 순간 제어계는 안정화 장치가 아니라 가진기가 된다.¹

플라이-바이-와이어가 표준이 된 이후 이건 예외 상황이 아니라 기본 상황이 됐다. 유연한 대형기, 고기동 전투기, 고고도 장기체공 무인기, 트러스 지지 초고세장비 날개 전부가 이 문제를 안고 설계된다. 개별 공탄성 현상 자체(플러터, 정적 발산, 버피팅, 조종면 반전)는 각 문서에 있으니 여기서는 제어계가 끼어들었을 때 무엇이 달라지는가만 본다.

2. 되먹임 고리의 해부[편집]

ASE 불안정이 성립하려면 다음 다섯 칸이 모두 채워져야 한다.

1. 센서 감지: 관성 센서가 자기 위치에서의 탄성 모드 성분을 읽는다. 센서를 그 모드의 절점(node line) 근처에 놓으면 이 항이 작아지므로, 센서 위치 선정은 그 자체로 ASE 설계 변수다.
2. 제어 법칙 통과: 그 신호가 이득을 받고 제어 법칙을 지난다. 대역폭이 넓을수록, 미분 보상이 많을수록 고주파가 잘 통과한다.
3. 작동기 응답: 조종면이 실제로 그 진동수에서 움직인다. 유압 작동기의 대역폭·유격·율제한이 여기 개입한다.
4. 공기력 발생: 움직인 조종면이 비정상 공기력을 만든다.
5. 모드 가진: 그 공기력이 원래 모드를 다시 밀어 준다. 위상이 맞으면 진폭이 커진다.

전형적으로 문제가 되는 것은 날개 1차 굽힘, 동체 1차 굽힘, 수직미익 굽힘, 외부 장착물이 붙은 파일런 모드처럼 5~30 Hz 대역이다. 강체 조종 대역(1 Hz 내외)과 충분히 떨어져 있어 보이지만, 제어 대역폭이 넓어질수록 그 여유가 좁아진다.

디지털 제어계에서는 함정이 하나 더 있다. 표본화 주파수의 절반(나이퀴스트)을 넘는 구조 모드는 앨리어싱되어 낮은 주파수로 접혀 들어온다. 60 Hz 모드가 표본화 주파수 80 Hz 계에서 20 Hz로 나타나 제어 대역 근처에서 말썽을 부리는 식이다. 그래서 안티앨리어싱 필터가 필수인데, 그 필터가 또 위상 지연을 만든다. 여기에 연산 지연까지 더해지면 상황은 산술적으로 나빠진다 — 20 ms 지연은 20 Hz에서 $360^\circ \times 20 \times 0.02 = 144^\circ$ 의 위상 지연이다. 고주파에서 위상은 순식간에 무너진다.

3. 구조 필터, 그리고 이득 안정화 대 위상 안정화[편집]

가장 고전적인 대책은 노치 필터(대역저지 필터)다. 문제가 되는 모드 진동수에 2차 노치를 걸어 그 주파수 성분만 잘라내면 루프가 닫히지 않는다.

간단해 보이지만 공짜가 아니다. 노치는 저지 대역 아래쪽 주파수에도 위상 지연을 남기고, 그 지연은 그대로 강체 제어 루프의 위상 여유를 깎아먹는다. 조종성을 위해 대역폭을 넓히고 싶은 제어 설계자와, 구조 모드를 확실히 죽이고 싶은 ASE 담당자가 같은 필터를 놓고 싸우는 이유다. 게다가 구조 모드 진동수는 연료 잔량, 외부 장착물 조합, 페이로드에 따라 움직인다. 한 점이 아니라 대역을 덮어야 하고, 대역을 넓히면 위상 손실은 더 커진다. 전투기가 장착물 형상마다 별도의 ASE 해석과 인증을 요구받는 근본 원인이 여기 있다.²

여기서 두 가지 전략이 갈린다.

• 이득 안정화(gain stabilization): 해당 모드 주파수에서 루프 이득을 0 dB 훨씬 아래로 눌러 버린다. 위상이 어떻든 불안정이 성립하지 않으므로 모드 진동수 추정 오차나 감쇠 불확실성에 강건하다. 대부분의 구조 모드는 이렇게 처리한다.
• 위상 안정화(phase stabilization): 이득을 살려 둔 채 위상을 맞춰 오히려 감쇠를 더해 준다. 모드가 제어 대역과 너무 가까워 눌러 버릴 수 없을 때, 그리고 능동 플러터 억제처럼 일부러 그 모드에 개입해야 할 때 쓴다. 대신 모드 특성 추정이 틀리면 그대로 사고로 이어지므로 불확실성 관리가 훨씬 빡세다.

4. 제어계를 아군으로 쓰기 — AFS·GLA·승차감[편집]

방어만 하는 분야는 아니다. 같은 되먹임 고리를 반대 방향으로 쓰면 구조가 얻는 것이 많다.

• 능동 플러터 억제(AFS): 센서가 읽은 굽힘·비틀림 운동에 위상을 맞춰 조종면을 움직여 인위적인 공력 감쇠를 만든다. 플러터 속도를 높이거나, 같은 속도에서 날개를 더 가볍게 만들 수 있다. 1970년대 이후 실험기와 풍동 모델에서 반복 실증됐고, NASA의 X-56A처럼 일부러 유연하게 만든 무인 시험기로 비행 실증까지 갔다.
• 돌풍 하중 경감(GLA)·기동 하중 경감(MLA): 돌풍이나 급기동으로 날개에 하중이 실릴 때 스포일러·플래퍼론·에일러론을 대칭 조작해 양력 중심을 안쪽으로 당긴다. 날개 뿌리 굽힘 모멘트의 첨두값이 줄면 구조 중량이 줄고, 결과적으로 더 큰 세장비의 날개를 쓸 수 있게 된다. 피로 해석 관점에서도 하중 스펙트럼의 진폭이 줄어드는 직접적 이득이 있다.
• 승차감·모드 안정화: 저고도 고속 침투 비행에서 동체 굽힘 진동은 조종사의 시계와 작업 능력을 갉아먹는다. B-52의 LAMS, B-1의 SMCS(기수의 작은 베인으로 동체 굽힘 모드를 감쇠)가 고전적 사례이고, 현대 여객기에도 난류 중 승차감 개선 기능이 들어간다.

그런데 현업에서 AFS가 민항기에 잘 안 실리는 이유가 있다. 감항 규정은 제어계가 고장 난 상태에서도 플러터 자유를 요구하는 쪽으로 설계돼 있어서, AFS에 기대어 기본 플러터 여유를 확보하는 설계는 고장 확률이 극히 희박함을 별도로 입증해야 한다. 그래서 GLA·MLA처럼 하중을 줄이는 기능은 널리 쓰이는 반면, 안정성 자체를 제어계에 맡기는 AFS는 여전히 조심스럽다.

5. 비정상 공기력을 상태공간에 집어넣기[편집]

ASE 해석에는 구조적인 난관이 하나 있다. 고전 플러터 해석의 공기력 영향계수 행렬 $\mathbf{Q}(ik, M_\infty)$ 는 조화 진동을 전제로 감소진동수 $k$ 의 함수로 표로 주어진다. 반면 제어계는 라플라스 영역의 상태공간·전달함수로 기술된다. 표로 된 주파수 응답과 상태공간을 한 루프에 넣을 수는 없으므로, $\mathbf{Q}$ 를 $s$ 의 유리함수로 근사해야 한다. 이것이 유리함수 근사(RFA) 다.

가장 널리 쓰이는 로저(Roger) 형식은 다음과 같다.

$$\mathbf{Q}(\bar{s}) \approx \mathbf{A}_0 + \mathbf{A}_1 \bar{s} + \mathbf{A}_2 \bar{s}^2 + \sum_{j=1}^{n_L} \mathbf{A}_{j+2} \, \frac{\bar{s}}{\bar{s} + \beta_j}$$

여기서 $\bar{s} = sb/V$ 는 무차원 라플라스 변수이고, $\beta_j$ 는 후류의 기억 효과를 나타내는 지연 극점(lag pole)이다. 앞의 세 항은 각각 준정상 강성·감쇠·부가질량에 대응하고, 합 기호 안의 항들이 비정상성(후류가 과거를 기억하는 효과)을 담당한다. 계수 행렬 $\mathbf{A}_j$ 는 표로 주어진 $\mathbf{Q}(ik)$ 에 최소자승법으로 맞춘다.

문제는 상태 수다. 지연 극점 하나마다 구조 모드 수만큼 상태가 추가되므로, 모드 50개에 극점 6개면 공력만으로 300개의 추가 상태가 생긴다. 제어 설계용 모델로는 감당이 안 된다. 그래서 카펠(Karpel)의 최소상태(minimum-state) 근사가 나왔다.

$$\mathbf{Q}(\bar{s}) \approx \mathbf{A}_0 + \mathbf{A}_1 \bar{s} + \mathbf{A}_2 \bar{s}^2 + \mathbf{D}\,(\bar{s}\mathbf{I} - \mathbf{R})^{-1}\mathbf{E}\,\bar{s}$$

행렬 $\mathbf{D}, \mathbf{E}$ 를 얇게 잡아 지연 상태를 수십 개 수준으로 줄인다. 대신 맞춤이 비선형 최적화가 되어 초기값 의존성이 생긴다. 이렇게 얻은 공력 상태를 모드 해석으로 축소한 구조 모델(축소차수모델), 작동기 전달함수, 센서 모델, 시간 지연, 그리고 제어 법칙과 하나로 묶으면 비로소 전체 ASE 플랜트가 완성된다. 그다음은 고유값 문제를 풀거나 루프를 끊어 개루프 주파수 응답을 뽑는, 익숙한 제어 해석의 세계다.

6. 안정성 여유와 인증[편집]

ASE의 합격 기준은 결국 여유(margin) 로 표현된다. 루프를 한 지점에서 끊고 개루프 주파수 응답을 그린 뒤, 각 구조 모드 피크에서 이득 여유와 위상 여유를 읽는다. 요구치는 규정·기종·모드 종류에 따라 다르지만 강체 모드에 이득 여유 6 dB·위상 여유 45° 남짓, 탄성 모드에는 그보다 넉넉한 이득 여유를 요구하는 것이 일반적인 그림이다.

이걸 해석만으로 통과시키지는 않는다.

• 지상 구조연성시험(SCT): 지상진동시험으로 모드 특성을 확정한 뒤, 유압을 올린 상태에서 비행제어계에 신호를 주입해 실제 개루프 응답을 측정한다. 해석 모델이 예측한 피크가 실제로 거기 있는지, 여유가 실제로 남아 있는지를 확인하고 필요하면 노치를 다시 판다. 항공기 개발 일정에서 조용히 애를 많이 먹이는 시험이다.
• 형상별 전수 검토: 연료 상태, 외부 장착물 조합, 페이로드, 그리고 고장 상태(센서 1개 고장, 작동기 성능 저하, 제어 법칙 저하 모드)까지 조합이 폭발한다. 조합마다 여유를 확인해야 하므로 실무에서는 강건 제어 쪽 도구 — 파라미터 불확실성을 묶어 한 번에 판정하는 구조적 특이값($\mu$) 해석 — 로 사례 수를 줄이려는 시도가 계속된다.
• 비행시험: 속도를 단계적으로 올리며 각 지점에서 구조를 가진하고 감쇠를 측정한다. 여기서 확인하는 것은 플러터 문서의 감쇠-속도 곡선과 같은 그림이되, 제어계를 켠 상태와 끈 상태 양쪽이다.

요약하면 ASE는 “제어 설계자가 구조를 몰라서, 구조 설계자가 제어를 몰라서” 터지는 사고를 막기 위해 생긴 분야다. 그리고 실제로 이 분야의 가장 흔한 실패 원인은 방정식이 아니라 조직도라는 말이 농담이 아니게 돌아다닌다.³

7. 관련 문서[편집]

• 항공탄성학 · 플러터 · 정적 발산 · 버피팅 · 조종면 반전
• 모드 해석 · 감쇠 · 주파수 응답 해석 · 축소차수모델 · 부분구조법
• 유체-구조 연성 · 다물리 연성해석 · 공력해석
• 고유값 문제 · 최소자승법 · 강건 제어 · 칼만 필터
• 지상진동시험 · 검증 및 확인 · 피로 해석 · NASTRAN

8. Footnotes[편집]

1. 이 현상의 옛 이름이 “구조연성(structural coupling)“인데, 어감이 너무 순하다. 실제로는 기체가 자기 자신의 자세제어기 손에 흔들리는 상황이다. 초기 FBW 개발 시절 지상시험에서 기체가 스스로 떨기 시작해 급히 전원을 내린 일화들이 여러 기종에 하나씩은 있다. ↩

2. 노치를 하나 팔 때마다 저주파 위상이 깎이므로, 노치가 여러 개 겹치면 정작 조종 특성이 둔해진다. 그래서 ASE 회의는 대개 “이 모드를 정말 눌러야 합니까”와 “그럼 조종성은 누가 책임집니까”의 무한 반복이다. 결론은 보통 센서를 옮기자는 쪽으로 난다 — 공짜에 가까운 유일한 해법이라서. ↩

3. 구조팀은 모드 데이터를 주고, 공력팀은 $\mathbf{Q}$ 표를 주고, 제어팀은 그걸 받아 이득을 정한다. 그런데 세 팀이 서로 다른 좌표계·부호 규약·모드 정규화를 쓰고 있다는 사실은 대개 통합 시험 전날 밤에 발견된다. 부호 하나 뒤집히면 감쇠를 더하려던 제어가 정확히 반대로 작동한다는 점이 이 분야를 특히 스릴 있게 만든다. ↩