Rayleigh 감쇠행렬

왜 배우는가: 레일리 감쇠는 다자유도 구조물의 시간이력 해석에서 감쇠행렬을 가장 간단하게 만드는 방법이다. 두 모드의 감쇠비만 정하면 전체 감쇠행렬이 결정되므로 실무 해석에서 자주 쓰인다. 시험에서는 감쇠행렬 구성 절차를 묻는 문항으로 반복해서 출제된다.

출제 이력: 제112회 1교시 9번, 제90회 3교시 2번

1. 한 줄 정의와 핵심 키워드

레일리 감쇠는 감쇠행렬을 질량행렬과 강성행렬의 선형 결합으로 놓는 비례감쇠 모델이다.

핵심 키워드

• 레일리 감쇠(Rayleigh damping): 감쇠행렬을 C = aM + bK로 두는 모델
• 비례감쇠: 감쇠가 질량과 강성에 비례한다고 가정하는 방식
• 질량 비례 계수 a: 낮은 진동수 모드의 감쇠에 크게 기여하는 계수
• 강성 비례 계수 b: 높은 진동수 모드의 감쇠에 크게 기여하는 계수
• 모드 감쇠비 ζn: n차 모드의 감쇠 정도를 나타내는 비
• 모드 직교성: 서로 다른 모드 사이에서 행렬 곱이 0이 되는 성질
• 고유진동수 ωn: n차 모드의 원주파수

2. 개념 풀이

쉽게 말하면, 구조물이 흔들릴 때 에너지가 줄어드는 정도를 질량과 강성의 크기에 맞춰 정하는 방법이다. 감쇠를 실험으로 행렬 전체에 대해 측정하기는 어려우므로, 이 방식은 해석에 필요한 감쇠행렬을 두 개의 계수로 만든다.

정식 정의는 다음과 같다. 다자유도계 운동방정식 M\ddot{u} + C\dot{u} + Ku = 0 에서 감쇠행렬 C를 C = aM + bK로 놓는 것을 레일리 감쇠라 한다. 여기서 M은 질량행렬, K는 강성행렬, a와 b는 비례 계수다. 단위는 a가 1/s, b가 s이며, 두 항의 단위가 모두 kg/s가 되어 감쇠행렬의 단위와 맞는다.

질량 비례 계수 a는 질량에 비례하는 감쇠 항의 계수다. 이 항은 진동수가 낮을수록 감쇠비를 키운다. 강성 비례 계수 b는 강성에 비례하는 감쇠 항의 계수로, 진동수가 높을수록 감쇠비를 키운다.

모드별 감쇠비와 고유진동수의 관계는 다음과 같다. $\zeta_n = \frac{a}{2\omega_n} + \frac{b\,\omega_n}{2}$ 여기서 \zeta_n은 n차 모드 감쇠비, \omega_n은 n차 고유진동수(rad/s)다. 이 식에서 두 계수를 구하면 감쇠행렬이 정해진다.

3. 등장 배경과 필요성

실제 구조물의 감쇠는 재료 내부 마찰, 접합부 미끄럼, 지반과 기초의 상호작용 등이 섞여 있다. 이 원인들을 감쇠행렬의 각 원소로 직접 나누어 측정하기 어렵다. 그래서 해석에서는 감쇠행렬을 단순한 식으로 가정할 필요가 생겼다.

모드해석은 감쇠를 무시해도 고유진동수와 모드형상을 얻을 수 있다. 그러나 지진이나 바람에 대한 시간이력 해석에서는 감쇠가 응답의 크기를 좌우하므로 감쇠행렬이 반드시 필요하다. 감쇠가 있는 문제는 고유치 문제가 이차식 형태로 바뀌어 직접 풀기가 복잡해진다.

레일리 감쇠가 널리 쓰이는 이유는 모드 좌표 분리가 가능하기 때문이다. 비감쇠 모드형상은 질량과 강성에 대해 직교한다. 따라서 C = aM + bK로 두면 감쇠행렬도 서로 다른 모드 사이에서 직교가 성립한다. 그 결과 운동방정식이 모드별 1자유도 방정식으로 분리된다.

이 분리 덕분에 해석이 단순해지고 계산 비용이 줄어든다. 다만 감쇠가 질량과 강성 분포에 비례한다는 가정이므로, 실제 감쇠 분포와 차이가 날 수 있다. 이 한계는 6절과 7절에서 다룬다. 두 모드의 감쇠비만 정하면 전체 감쇠행렬이 결정된다. 즉 두 개의 모드 감쇠비 입력이 해석에서 가장 중요한 가정이다.

4. 구성요소·원리·절차

Rayleigh 감쇠행렬(Rayleigh Damping Matrix)은 질량행렬과 강성행렬의 선형 결합으로 만든다. 감쇠행렬 C는 다음 식으로 정한다.

$C = a_0 M + a_1 K$

여기서 M은 질량행렬(kg), K는 강성행렬(N/m), a0는 질량비례 계수(1/s), a1은 강성비례 계수(s)다. 두 계수를 정하면 감쇠행렬 C가 결정된다. 이 방식을 비례감쇠(Proportional Damping)라고 하며, 감쇠행렬이 M과 K에 비례하므로 모드 좌표계에서 대각 형태를 유지한다.

위 표는 계수의 역할을 보인다. a0 항은 저진동수 모드의 감쇠를 키우고, a1 항은 고진동수 모드의 감쇠를 키운다. 모드 n의 감쇠비는 다음과 같이 정해진다.

$\zeta_n = \frac{1}{2}\left(\frac{a_0}{\omega_n} + a_1\omega_n\right)$

여기서 ζn은 n차 모드 감쇠비, ωn은 n차 고유원진동수(rad/s)다. 이 식은 질량항과 강성항이 각 모드 감쇠비에 기여하는 방식을 보여준다.

구성 절차는 다음과 같다.

• 해석 대상 구조물의 질량행렬 M과 강성행렬 K를 구한다.
• 고유치 해석으로 대상 두 모드의 고유원진동수 ω1, ω2를 구한다. 보통 지배 모드인 1차 모드를 포함한다.
• 두 모드의 목표 감쇠비 ζ1, ζ2를 정한다.
• 다음 식으로 a0와 a1을 구한다.

$a_0 = \frac{2\omega_1\omega_2(\zeta_1\omega_2 - \zeta_2\omega_1)}{\omega_2^2 - \omega_1^2}, \qquad a_1 = \frac{2(\zeta_2\omega_2 - \zeta_1\omega_1)}{\omega_2^2 - \omega_1^2}$

• C = a0M + a1K로 감쇠행렬을 구성하고, 다른 모드의 감쇠비 ζn을 확인한다.

계산 예제는 다음과 같다. 주어진 값은 ω1 = 10rad/s, ω2 = 30rad/s, ζ1 = ζ2 = 0.05이다. 분자는 2·10·30·(0.05·30 − 0.05·10) = 600·1 = 600이고, 분모는 30² − 10² = 800이다. 따라서 a0 = 600/800 = 0.75 1/s이다. a1은 2(0.05·30 − 0.05·10)/800 = 2/800 = 0.0025 s이다. 검산으로 ζ1 = (0.75/10 + 0.0025·10)/2 = 0.05이고 ζ2 = (0.75/30 + 0.0025·30)/2 = 0.05이다. 두 값이 목표 감쇠비를 재현하므로 결과는 상식 범위에 있다.

5. 관련 법령·기준·표준

Rayleigh 감쇠 자체는 법령 조항으로 정해지지 않는다. 국내 교량과 건축물의 내진·동적 해석에서는 설계기준이 해석 조건으로 감쇠비를 요구한다. 아래 표는 기준값과 적용 조건을 정리한 것이다. 조항 번호는 확실한 것만 쓰고, 나머지는 기준 이름만 적는다.

외국 기준과의 비교는 감쇠비 정의와 계수 산정 방식이 같은지를 보는 수준으로 한정한다.

6. 유형과 비교

Rayleigh 감쇠는 질량과 강성에 비례하는 감쇠행렬이다. 실무에서는 감쇠를 직접 행렬로 정의하기 어려워 이 모델을 많이 쓴다. 시험에서는 다른 감쇠 모델과 구분하는 문제로 자주 나온다.

위 표에서 Rayleigh 감쇠와 모달 감쇠는 모두 고유벡터 기저에서 대각행렬이 된다. 차이는 Rayleigh가 두 계수로 여러 모드의 감쇠비를 간접적으로 정한다는 점이다. 따라서 중간 모드의 감쇠가 실제와 다를 수 있다.