레일리 감쇠모델

왜 배우는가: 구조물 동적해석에서는 감쇠행렬을 직접 정의하기 어려워 감쇠 모델을 가정해야 한다. 레일리 감쇠는 두 개의 계수만으로 감쇠행렬을 만들어 주파수에 따른 감쇠비를 쉽게 맞출 수 있어, 시간이력해석과 내진 해석에서 반복적으로 쓰이고 시험에서도 계수 산정과 감쇠비 계산 문제로 자주 나온다.

출제 이력: 제128회 1교시 1번, 제117회 3교시 1번

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

답안 첫 줄 정의: 레일리 감쇠는 질량행렬과 강성행렬의 선형 결합으로 감쇠행렬을 가정하는 비례감쇠 모델이다.

핵심 키워드

• 비례감쇠(Proportional Damping): 감쇠행렬이 질량행렬과 강성행렬에 비례하는 감쇠.
• 감쇠행렬(Damping Matrix, C): 구조물의 감쇠를 행렬로 나타낸 것.
• 질량 비례항: 감쇠행렬 중 질량행렬에 비례하는 항(계수 a0).
• 강성 비례항: 감쇠행렬 중 강성행렬에 비례하는 항(계수 a1).
• 감쇠비(Damping Ratio, ζ): 임계감쇠 대비 실제 감쇠의 비율.
• 고유진동수(Natural Frequency, ω): 구조물의 고유 원진동수.
• 모드 직교성: 감쇠행렬이 모드형상에 대해 대각화되어 모드별 방정식이 분리되는 성질.

2. 개념 풀이

쉬운 말로 하면 레일리 감쇠는 "질량에 비례하는 감쇠와 강성에 비례하는 감쇠를 더해 구조물 전체의 감쇠를 만드는 방법"이다. 질량이 크면 저주파에서 감쇠가 크게 작용하고, 강성이 크면 고주파에서 감쇠가 크게 작용한다는 직관에 맞춘 모델이다.

정식 정의: 레일리 감쇠는 감쇠행렬 C를 질량행렬 M과 강성행렬 K의 선형 결합으로 가정한다. 계수 a0와 a1은 두 개의 기준 고유진동수에서 목표 감쇠비를 맞추도록 정한다. 이때 감쇠행렬은 각 모드의 고유벡터에 대해 대각화되므로, 다자유도계의 운동방정식이 모드별 단자유도계 방정식으로 분리된다.

$\mathbf{C}=a_0\mathbf{M}+a_1\mathbf{K}$

여기서 C는 감쇠행렬(N/(m/s)), M은 질량행렬(kg), K는 강성행렬(N/m), a0는 질량 비례계수(1/s), a1는 강성 비례계수(s)이다.

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

여기서 ζn은 n차 모드의 감쇠비, ωn은 n차 모드의 원진동수(rad/s)이다. 이 식을 보면 a0 항은 주기가 긴 저주파 모드의 감쇠비를 키우고, a1 항은 주기가 짧은 고주파 모드의 감쇠비를 키운다. 따라서 두 계수는 주기가 다른 두 모드의 감쇠비를 목표값으로 맞추는 데 쓰인다.

두 기준 모드 ω1, ω2에서 감쇠비가 같은 값 두 고유진동수에서 감쇠비 5%로 주어지면, 계수는 다음과 같이 구한다.

$a_0=\frac{2\zeta\omega_1\omega_2}{\omega_1+\omega_2},\qquad a_1=\frac{2\zeta}{\omega_1+\omega_2}$

이 식에서 ζ는 기준 감쇠비, ω1·ω2는 기준 원진동수(rad/s)다. 두 기준 모드 사이의 주기에서는 감쇠비가 목표값보다 작게 나오고, 기준 모드 바깥의 주기에서는 감쇠비가 커지는 것이 레일리 감쇠의 특징이다.

핵심 용어 풀이

• 비례감쇠: 감쇠행렬을 질량과 강성의 조합으로 만든 감쇠 모델이다. 모드 해석이 가능하도록 하는 것이 목적이다.
• 모드 직교성: 질량행렬과 강성행렬의 고유벡터에 대해 감쇠행렬도 직교하므로, 모드 좌표계에서 방정식이 독립적으로 풀린다.
• 감쇠비: 감쇠계수 c를 임계감쇠계수 cc=2mω로 나눈 값이다. 일반적으로 강구조는 약 2~5%를 쓴다.

3. 등장 배경과 필요성

구조물의 감쇠는 재료 내부 마찰, 접합부 마찰, 비구조요소와의 상호작용 등 여러 원인이 섞여 있어 실험으로 측정할 수 있는 것은 모드별 감쇠비뿐이다. 다자유도계 구조물에서 감쇠행렬을 요소 단위로 정의하려면 감쇠 메커니즘을 모두 알아야 하므로 현실적으로 불가능하다. 레일리 감쇠는 이 문제를 해결하기 위해 두 개의 계수만 정하면 되는 간편한 모델로 제안됐다.

또한 감쇠행렬이 모드별로 분리되지 않으면, 시간이력해석에서 전체 연립방정식을 직접 풀어야 한다. 비례감쇠 가정은 모드 중첩법(Mode Superposition)을 쓸 수 있게 하여 계산 부담을 줄인다. 설계 실무에서는 1차 모드와 주요 고차 모드의 감쇠비를 목표값으로 맞춘 뒤, 나머지 모드의 감쇠비가 과도하게 작은지 확인하는 순서로 활용된다.

이 모델의 필요성은 다음 세 가지로 정리할 수 있다.

• 감쇠행렬을 실험 감쇠비로부터 역으로 결정할 수 있다.
• 모드별 방정식이 분리되어 해석 효율이 높다.
• 감쇠 메커니즘을 정확히 모르더라도 해석 모델을 만들 수 있다.

한계도 함께 알아야 한다. 레일리 감쇠는 물리적 감쇠 기구를 반영하지 않는 수학적 가정이므로, 기준 모드에서 멀리 떨어진 모드의 감쇠비는 실제와 차이가 날 수 있다. 따라서 기준 주기를 해석에서 중요한 모드의 주기 범위에 맞춰 선택하는 것이 중요하다.

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

레일리 감쇠(Rayleigh damping)는 감쇠행렬을 질량행렬과 강성행렬의 선형 결합으로 둔다. 물리적인 감쇠 장치를 개별로 모델링하지 않고, 해석 모델 전체에 감쇠를 한 번에 부여하는 방법이다.

레일리 감쇠의 원리는 모드별 감쇠비로 설명된다. 모드 i의 고유원진동수를 ω_i라 하면 모드 감쇠비는 다음과 같다.

$\zeta_i=\frac{1}{2}\left(\frac{a_0}{\omega_i}+a_1\omega_i\right)$

여기서 ζ_i는 모드 i의 감쇠비(무차원), ω_i는 모드 i의 고유각진동수(rad/s)다. 질량 비례 항은 진동수가 낮을수록 감쇠에 크게 기여하고, 강성 비례 항은 진동수가 높을수록 크게 기여한다. 따라서 ζ(ω)를 그리면 아래로 볼록한 곡선이 된다. 두 기준 주기 사이 구간은 설계 감쇠비보다 작고, 구간 밖에서는 감쇠비가 커진다. 이 곡선의 모양이 레일리 감쇠의 한계를 결정한다.

두 기준 진동수를 정하고 계수를 구하는 절차는 다음과 같다.

• 해석 대상의 고유치 해석으로 주요 모드의 고유주기 T1, T2를 구한다. 보통 1차 모드와 질량참여계수가 큰 모드를 기준으로 삼는다.
• 각 기준 주기를 고유각진동수로 바꾼다. ω = 2π/T 를 적용한다.
• 두 기준 모드에 목표 감쇠비 ζ를 부여하여 연립방정식을 세운다.
• 연립방정식을 풀어 a0, a1을 구한다.
• C = a0 M + a1 K 로 감쇠행렬을 구성하고 시간이력 또는 모드 중첩 해석에 적용한다.

두 기준 주기 ω1, ω2에서 ζ 값을 같게 두면 다음 해를 얻는다.

$a_0=\frac{2\zeta\,\omega_1\omega_2}{\omega_1+\omega_2},\qquad a_1=\frac{2\zeta}{\omega_1+\omega_2}$

ζ는 목표 감쇠비(무차원), ω1·ω2는 기준 모드의 고유각진동수(rad/s)다. 이 식은 두 식을 ω1², ω2²로 소거하여 얻는다. 주기 단위로 나타내면 계산이 편하다.

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

레일리 감쇠는 특정 법령이 정한 공식이 아니라 구조 해석 기법이다. 따라서 법령·기준은 그 해석 결과에 적용되는 감쇠비와 해석 절차를 정한다. 국내 내진설계기준에서 응답스펙트럼 해석에 쓰는 일반 건축물의 기준 감쇠비는 감쇠비 5%로 보는 것이 일반적이다. 레일리 감쇠의 목표 감쇠비도 이 값을 출발점으로 삼는다.

위 표의 구조형식별 감쇠 기준과 비선형 해석에서의 중복 검토는 해당 설계기준 해설에 따라 적용한다. 기준 조항의 세부 번호는 이 절에서 제시하지 않는다.

레일리 감쇠 계수 산정 시 실무에서 자주 생기는 문제는 기준 주기를 잘못 고르는 것이다. 기준 주기가 실제 응답을 지배하는 모드와 떨어지면 해당 모드의 감쇠비가 목표값과 달라진다. 해석 결과를 제출할 때는 두 기준 주기와 계수 a0, a1의 값을 함께 밝힌다.

국제 기준에서도 레일리 감쇠 계수는 같은 원리로 결정된다. 개념 차이는 감쇠비 목표값과 적용 모드 선택의 기준에 있을 뿐이며, 계산식 자체는 동일하다.

6. 유형과 비교

레일리 감쇠는 질량비례 감쇠와 강성비례 감쇠의 합으로 구성되며, 두 계수 a_0, a_1의 결정 방식에 따라 여러 형태로 나눠 볼 수 있다. 실무에서 흔히 비교하는 대상은 레일리 감쇠, 모드 감쇠(모드별 감쇠비 직접 지정), 그리고 일정 감쇠비 가정(고유모드마다 같은 감쇠비)이다.