응답진폭연산자

1. 개요[편집]

응답진폭연산자 Response Amplitude Operator (RAO)
정체	선형 시스템의 주파수응답 함수
정의	단위 파 진폭당 응답 진폭과 위상
취득	규칙파 시험 · 주파수영역 포텐셜 해석
사용법	응답 스펙트럼 = RAO² × 파 스펙트럼
필수 변환	조우 주파수 (선속 · 조우각)
전제	선형성 · 정상성 · 좁은 대역 근사

배가 파도에 어떻게 반응하는지 알고 싶다고? 한 주파수씩 때려 보고 표를 만들면 된다. 라플라스가 200년 전에 알려 준 그 방법 그대로.

응답진폭연산자(Response Amplitude Operator, RAO)는 파랑을 입력, 부유체의 운동을 출력으로 보는 선형 시스템의 주파수응답 함수다. 정의는 한 줄이다 — 진폭 $A$ 인 규칙파에서 어떤 자유도의 응답 진폭이 $\zeta_a$ 라면

$$\mathrm{RAO}(\omega) = \frac{\zeta_a(\omega)}{A}, \qquad \angle H(\omega) = \varepsilon(\omega)$$

이고, 진폭비와 위상을 묶어 복소 전달함수 $H(\omega)$ 로 다룬다.¹ 이름에 “진폭”이 들어가 있어 크기만 뜻하는 것처럼 보이지만, 위상을 버리면 여러 자유도를 합쳐 한 지점의 상대 운동이나 가속도를 만들 수 없다. 상대 운동·국부 가속도가 필요한 순간 위상은 필수다.

이 문서는 도구를 다룬다. RAO 로 판정하는 현상과 설계 결과는 내항성이, 수평면 조종 운동은 선박 조종성이 담당한다. 개념적으로는 주파수 응답 해석과 같은 물건이고, 지진공학의 응답 스펙트럼 해석과는 사촌이되 결정적으로 다르다 — 저쪽은 최대응답 곡선을 곱하고 이쪽은 파워 스펙트럼을 곱한다.

2. 무엇의 비인가 — 정규화 관례[편집]

RAO 는 무차원화 방식이 자유도마다 다르고, 이걸 틀리면 그래프 전체가 무의미해진다.

자유도	표준 정규화	저주파 극한
상하동요	파 진폭 $A$	1
종동요 · 횡동요	파 경사 $kA$	1
가속도	$\omega_e^2 \times$ 변위 RAO	0
상대 운동	파 진폭 $A$	0

종동요를 파 경사 $kA$ 로 나누는 이유는 장파 극한에서 배가 파면에 얹혀 파면 기울기를 그대로 따라가기 때문이다. 그렇게 정규화해야 저주파 극한이 깔끔하게 1 이 되고, 서로 다른 배의 곡선을 겹쳐 볼 수 있다. 반면 실무 보고서에는 도/미터(deg/m) 단위로 그대로 찍은 곡선도 흔하니, 남의 RAO 를 받으면 단위와 정규화 기준부터 확인하는 것이 예의이자 생존 기술이다.

3. 어떻게 얻나[편집]

3.1. 규칙파 시험[편집]

수조에서 조파기로 한 주파수의 규칙파를 만들고, 정상 상태 응답의 진폭과 파고계 신호에 대한 위상차를 잰다. 주파수를 바꿔 가며 이 짓을 20~40번 반복하면 곡선 하나가 나오고, 조우각과 선속 조합마다 이걸 또 반복한다. 시험 계획이 조합 폭발과의 싸움이 되는 이유이고, 예인수조와 해양공학수조의 일정이 늘 빡빡한 이유다. 신호 처리는 고속 푸리에 변환으로 조우 주파수 성분만 뽑아내는 것이 표준이고, 고조파 성분의 크기가 곧 선형 가정이 얼마나 살아 있는가의 척도가 된다.

과도 시험으로 지름길을 낼 수도 있다. 백색에 가까운 짧은 파군을 때리고 입출력 신호를 함께 변환해 $H(\omega) = S_{xy}(\omega)/S_{xx}(\omega)$ 로 한 번에 전 주파수를 얻는 방식인데, 신호 대 잡음비가 나빠 검증용으로 더 많이 쓴다.

3.2. 주파수영역 해석[편집]

훨씬 흔한 경로는 계산이다. 포텐셜 유동 가정 아래 라디에이션·디프랙션 문제를 풀어 부가질량 $A_{jk}(\omega)$, 방사감쇠 $B_{jk}(\omega)$, 기진력 $F_j(\omega)$ 를 얻고 6×6 복소 연립방정식을 푼다.

$$\left[-\omega_e^{2}\left(\mathbf{M}+\mathbf{A}\right)+i\,\omega_e\mathbf{B}+\mathbf{C}\right]\boldsymbol{\xi} = \mathbf{F}$$

$\boldsymbol{\xi}$ 를 파 진폭으로 나누면 그대로 RAO 다. 스트립 이론이면 초 단위, 3차원 경계요소법이면 분~시간 단위로 곡선 한 벌이 나온다. 계보와 각 기법의 한계는 내항성 문서 참고.

4. 곡선의 세 영역[편집]

RAO 곡선은 어느 부유체든 대체로 세 구간으로 읽힌다. 이 구조를 알면 남의 그래프에서 계산 실수를 눈으로 잡아낼 수 있다.

저주파(장파) 영역. $\lambda \gg L$ 이면 배는 파면 위에 놓인 부표나 다름없다. 관성항 $\omega^2(\mathbf{M}+\mathbf{A})$ 와 감쇠항이 복원 강성 $\mathbf{C}$ 에 비해 무시할 만해지고, 응답은 정수압 평형만으로 결정되는 준정적 거동이 된다. 그래서 상하동요 RAO 는 1 로, 종동요는 파 경사비로 1 로 수렴한다. 여기서 1 로 안 가는 곡선은 대체로 계산이 틀렸다는 신호다.

공진 영역. $\omega_e$ 가 고유진동수에 접근하면 관성과 복원이 상쇄되어 응답이 감쇠에만 지배된다. 봉우리의 높이는 곧 감쇠의 역수다. 상하·종동요는 조파 방사로 에너지를 많이 버려 봉우리가 완만한 편(대개 2 를 크게 넘지 않는다)이지만, 횡동요는 방사감쇠가 거의 없어 봉우리가 뾰족하고 높다. 포텐셜 이론만으로 계산한 횡동요 RAO 가 현실의 몇 배로 솟는 것이 이 때문이고, 점성 감쇠 모형이나 자유횡동요 감쇠 시험이 반드시 따라붙는다.

고주파(단파) 영역. $\lambda \ll L$ 이면 파의 압력 변동이 선체 길이 방향으로 여러 번 부호를 바꿔 적분하면 서로 상쇄된다. 배는 잔물결을 그냥 무시하고, RAO 는 0 으로 떨어진다. 중간에 나타나는 상쇄 골(cancellation notch)은 여기서 한 발 더 나간 현상이다 — 특정 파장에서 선체 앞뒤 부분이 받는 기진력이 정확히 반대 위상이 되어 합이 0 에 가까워지면 RAO 가 깊게 파인다. 골의 위치는 선체 길이 방향 부력 분포가 결정하므로, 선형(hull form)마다 다르고 설계로 옮길 수 있다. 운항 파랑의 지배 주기 근처에 상쇄 골이 오도록 배치하는 것이 내항 성능 설계의 고전적인 수법이다.²

5. 조우 주파수 — 배가 움직이면 파도 주파수가 달라진다[편집]

정지한 부유체가 아니면 이야기가 하나 복잡해진다. 배가 속도 $U$ 로 달리면 배가 느끼는 파의 주파수는 절대 주파수 $\omega$ 가 아니라 조우 주파수 $\omega_e$ 다. 심해파에서 파수가 $k=\omega^2/g$ 이므로

$$\omega_e = \omega - \frac{\omega^{2}U}{g}\cos\mu$$

여기서 $\mu$ 는 조우각으로, $\mu=180°$ 가 선수파(head sea), $\mu=0°$ 가 추파(following sea), $\mu=90°$ 가 횡파(beam sea)다. 선수파에서는 $\cos\mu=-1$ 이라 $\omega_e>\omega$ — 파도를 더 자주 만난다. 추파에서는 반대로 $\omega_e<\omega$ 이고, $U\cos\mu$ 가 충분히 크면 파를 뒤에서 추월해 $\omega_e$ 가 음수가 되기도 한다. 배가 파보다 빨리 가면 파가 뒤로 흘러가는 것으로 보인다는, 물리적으로는 당연하지만 부호 관리에서는 성가신 상황이다.

여기서 실무 사고가 가장 자주 난다. RAO 는 조우 주파수의 함수인데 파 스펙트럼은 절대 주파수의 함수다. 순진하게 스펙트럼을 조우 주파수 축으로 변환하면

$$S(\omega_e) = \frac{S(\omega)}{\left|\,1 - \dfrac{2\omega U\cos\mu}{g}\,\right|}$$

인데, 추파 쪽에서 분모가 0 이 되는 주파수가 존재해 스펙트럼이 발산한다. 이건 물리적 특이점이 아니라 좌표 변환의 야코비안이 죽는 자리다. 해법은 단순하다 — 절대 주파수 $\omega$ 로 적분하고, RAO 만 $\omega_e(\omega)$ 에서 값을 읽어 온다. 변환을 아예 하지 않으면 특이점도 없다. 남의 코드에서 추파 조건의 응답이 터무니없이 크게 나오면 십중팔구 여기다.

6. 응답 스펙트럼과 통계[편집]

RAO 가 있으면 불규칙 해상의 응답 통계는 곱셈 한 번과 적분 몇 번으로 끝난다. 선형 시스템에 정상 가우스 입력이 들어가면 출력도 정상 가우스이고, 파워 스펙트럼은 전달함수 크기의 제곱만큼 배로 곱해진다.

$$S_R(\omega) = \left|H\!\left(\omega_e(\omega)\right)\right|^{2} S_\zeta(\omega)$$

여기서 $S_\zeta$ 는 유의파고와 평균주기 두 모수로 주어지는 브레치나이더/ITTC 스펙트럼이나 JONSWAP 스펙트럼이다. 방향 분산이 있으면 방향 확산 함수를 곱해 조우각에 대해 한 번 더 적분한다.

그 다음은 스펙트럼 모멘트다.

$$m_n = \int_0^{\infty} \omega_e^{\,n}\, S_R(\omega)\, d\omega$$

$n$ 승에 붙는 것이 조우 주파수라는 점이 중요하다. $m_0$ 은 응답의 분산이라 어느 축으로 재도 같지만, $m_2$ 는 “배가 겪는 진동 횟수”와 직결되므로 반드시 조우 기준이어야 한다. 여기서 나오는 값들:

• RMS(표준편차) $\sigma = \sqrt{m_0}$
• 유의 진폭 $= 2\sqrt{m_0}$ — 큰 쪽 1/3 의 평균 진폭. 유의 파고와 비교하려면 이것의 두 배인 유의 양진폭 $4\sqrt{m_0}$ 를 쓴다.
• 평균 조우 주기 $T_z = 2\pi\sqrt{m_0/m_2}$
• 대역폭 모수 $\varepsilon = \sqrt{1 - m_2^2/(m_0 m_4)}$ — 0 에 가까우면 좁은 대역.

좁은 대역 가정이 서면 응답 진폭은 레일리 분포를 따른다.

$$P(a > a_0) = \exp\!\left(-\frac{a_0^{2}}{2 m_0}\right)$$

이 한 줄이 슬래밍 발생 확률, 갑판 침수 확률, 특정 가속도 초과 확률을 전부 만들어 낸다. 예컨대 선수 상대 운동의 분산 $m_0$ 을 구했으면 건현 $f$ 를 넘을 확률이 $\exp(-f^2/2m_0)$ 이고, 그 값을 내항성 판정 기준의 한계치와 비교하는 식이다.

노출 시간 동안 $N$ 회 진동한다면 최대값의 기댓값은

$$E[a_{\max}] \approx \sqrt{2 m_0}\left(\sqrt{\ln N} + \frac{\gamma}{2\sqrt{\ln N}}\right)$$

로 $\gamma\approx0.5772$ 는 오일러-마스케로니 상수다. $\ln N$ 의 제곱근이라 관측 시간을 열 배 늘려도 최대값은 얼마 안 커진다 — 극값 통계 특유의 무딘 성장이다.³ 더 긴 재현주기(설계 수명 20년, 100년 재현 파랑)로 외삽하려면 이 단기 통계를 해상 상태의 장기 분포와 합성해야 하고, 그 영역은 극값 통계의 몫이다.

7. 어디까지 믿을 수 있나[편집]

RAO 절차 전체는 선형성 + 정상성이라는 두 기둥 위에 서 있다. 둘 중 하나가 무너지면 곱셈 한 번으로 끝나는 이 우아함도 함께 무너진다.

• 응답이 파 진폭에 비례해야 한다. 상하·종동요는 어지간히 큰 파에서도 선형성이 잘 유지되는 편이지만, 횡동요는 점성 감쇠가 진폭에 따라 변하고 대각도에서 복원정 곡선이 휘어 파고를 두 배 키우면 응답이 두 배가 안 된다. 그래서 횡동요 RAO 는 “어느 파고에서 잰 것인가”를 반드시 명시한다.
• 문턱 사건은 원리적으로 밖에 있다. 슬래밍·갑판 침수·프로펠러 노출은 상대 운동이 문턱을 넘어야 발생하는 사건이고, 발생한 뒤의 하중은 선형 이론이 다루는 물리가 아니다. RAO 가 줄 수 있는 것은 발생 확률까지다.
• 2차 효과는 안 나온다. 파랑 중 부가저항, 계류된 부유체의 저주파 표류 운동은 파 진폭의 제곱에 비례하는 2차 효과라, 1차 RAO 만으로는 절대 계산되지 않는다. 이쪽은 2차 전달함수(QTF)라는 별도의 이중 주파수 행렬이 필요하다.
• 매개변수 불안정은 아예 존재하지 않는다. 파라메트릭 횡동요처럼 계수가 시간에 따라 변해 생기는 불안정은 진폭이 0 에서 자라나는 현상이라 주파수응답이라는 개념 자체가 정의되지 않는다. 시간영역 해석이나 전산유체역학으로 넘어가는 수밖에 없다.
• 입력이 정상 가우스여야 한다. 태풍이 지나가는 몇 시간을 하나의 스펙트럼으로 뭉뚱그릴 수는 없다. 실무는 3시간을 정상 구간으로 보는 관례를 쓴다.

그럼에도 이 절차가 반세기 넘게 표준으로 살아남은 것은, 대안이 한 조건마다 수백 파 주기를 시간영역으로 돌리는 것이기 때문이다. 선속 5개 × 조우각 13개 × 해상 상태 20개면 1300 케이스인데, 주파수영역이면 노트북으로 점심 전에 끝나고 CFD 면 몇 달이다. 실무의 분업은 여기서 결정된다 — RAO 로 전 범위를 훑고, 걸리는 몇 케이스만 비선형 해석으로 확인한다.

8. 관련 문서[편집]

• 내항성 · 선박 조종성 · 슬래밍 · 모형실험 · 예인수조
• 주파수 응답 해석 · 응답 스펙트럼 해석 · 모드 해석 · 감쇠
• 파 스펙트럼 · 레일리 분포 · 확률과정 · 극값 통계
• 푸리에 변환 · 고속 푸리에 변환 · 중심극한정리
• 포텐셜 유동 · 경계요소법 · 전산유체역학
• 피로 해석 · 불확실성 정량화 · 검증 및 확인

9. Footnotes[편집]

1. “Operator” 라는 이름은 이 물건이 원래 입력 함수를 출력 함수로 보내는 선형 연산자의 커널로 도입됐기 때문이다. 시간영역에서 보면 임펄스 응답과의 합성곱이고, 그 푸리에 변환이 RAO 다. 이름만 보고 무슨 특별한 조선공학 개념인 줄 알았다면, 미안하지만 그냥 보드 선도의 크기 곡선이다. ↩

2. 상쇄 골의 존재는 “특정 파장에서만 유난히 안 흔들리는 배”를 가능하게 하지만, 반대로 설계 파장이 살짝만 어긋나면 그 이득이 통째로 사라진다는 뜻이기도 하다. 좁은 대역에 최적화된 설계가 실제 바다에서 배신당하는 전형적인 구조라, 요즘은 골의 깊이보다 곡선 전체를 파 스펙트럼과 겹쳐 적분한 값을 목적함수로 삼는다. ↩

3. $\sqrt{\ln N}$ 성장을 처음 보면 안심되지만, 방심하면 안 되는 것은 이게 분포의 꼬리가 지수적으로 얇다는 가정에서 나온 결과라는 점이다. 슬래밍 압력처럼 꼬리가 두꺼운 양에 같은 공식을 쓰면 최대값을 태연하게 과소평가한다. “레일리를 쓸 수 있는가”는 계산 편의의 문제가 아니라 물리의 문제다. ↩