회전체 동역학

시뮬레이션 구조해석 물리 마지막 수정: 2026-07-16 04:41:05

1. 개요[편집]

회전체 동역학 Rotordynamics
대상	터빈 · 압축기 · 모터 · 터보차저 축계
핵심 개념	자이로스코프 효과, 임계속도, 세차
대표 도구	캠벨 선도, 불평형 응답 곡선
고전 모델	젭코트(Jeffcott) 로터

회전체 동역학(rotordynamics)은 축을 중심으로 회전하는 구조물 — 터빈 로터, 압축기 축, 모터 샤프트, 터보차저 — 의 진동을 다루는 구조해석의 한 분과이다. 일반 구조동역학과 갈라지는 지점은 회전 자체가 새로운 힘과 새로운 방정식 항을 만든다는 것이다. 가만히 있으면 안 나타나던 자이로스코프 모멘트, 불평형 원심력, 베어링 유막의 비대칭 강성이 회전수와 함께 살아난다.

핵심 명제 하나만 기억하면 된다: 회전체의 고유진동수는 회전수에 따라 변한다. 이 한 줄 때문에 모드 해석 한 번으로 끝나던 일이 회전수 전체에 걸친 지도 그리기가 된다.

2. 젭코트 로터[편집]

젭코트 로터(Jeffcott rotor)는 무질량 탄성축 한가운데에 강체 디스크 하나를 얹고 양 끝을 단순 지지한 모델로, 회전체 동역학의 조화진동자다.¹ 디스크의 질량중심이 기하중심에서 편심 $e$ 만큼 어긋나 있으면, 회전수 $\Omega$ 에서 원심력 $m e \Omega^2$ 이 회전과 함께 도는 가진력으로 작용한다. 정상상태 진폭은

r = \frac{e \, (\Omega/\omega_n)^2}{\sqrt{\left(1 - (\Omega/\omega_n)^2\right)^2 + \left(2\zeta \Omega/\omega_n\right)^2}}

이 식이 알려주는 게 셋이다. (1) $\Omega = \omega_n$ 에서 진폭이 폭발한다 — 임계속도. (2) $\Omega \gg \omega_n$ 이면 $r \to e$ 로 수렴한다. 즉 임계속도를 넘겨서 운전하면 축은 질량중심을 축으로 스스로 돌기 시작한다. 이것이 자동 중심 잡기(self-centering)이며, 고속 회전기계가 임계속도를 굳이 통과해 초임계(supercritical) 영역에서 운전하는 이유다. (3) 진폭이 아무리 커도 편심 $e$ 에 비례하므로, 밸런싱이 전부다.

3. 자이로스코프 효과[편집]

디스크가 축에 대해 기울어진 채로 회전하면, 그 기울어짐의 변화율(세차 운동)과 스핀 각운동량이 곱해져 자이로스코프 모멘트가 생긴다. 운동방정식에 이 항이 들어가면 이렇게 된다.

\mathbf{M}\ddot{\mathbf{u}} + (\mathbf{C} + \Omega \mathbf{G})\dot{\mathbf{u}} + \mathbf{K}\mathbf{u} = \mathbf{F}(t)

$\mathbf{G}$ 가 자이로스코프 행렬이고, 결정적으로 반대칭( $\mathbf{G}^T = -\mathbf{G}$ )이다. 감쇠 행렬 자리에 앉아 있지만 에너지를 소산시키지 않는다 — 대칭이 아니기 때문이다. 자이로 항은 에너지를 빼앗는 대신 모드를 회전 방향에 따라 갈라놓는다.

계수에 $\Omega$ 가 붙어 있다는 점을 보라. 회전수가 시스템 행렬 자체를 바꾼다. 이 순간부터 고유값 문제는 $\Omega$ 에 의존하고, 고유값은 실수가 아니라 복소수가 되며, 모드 형상은 복소 벡터가 된다. 일반 고유값 문제 풀이가 아니라 비대칭 복소 고유값 문제를 풀어야 한다.

3.1. 전방 세차와 후방 세차[편집]

세차(whirl)는 축 중심이 그리는 궤도 운동이다. 방향이 스핀과 같으면 전방 세차(forward whirl), 반대면 후방 세차(backward whirl).

자이로스코프 효과는 전방 세차 모드를 뻣뻣하게(강성 증가), 후방 세차 모드를 무르게 만든다. 그래서 정지 상태에서 하나였던 고유진동수가 회전을 시작하면 두 갈래로 갈라진다. 위로 뻗는 가지가 전방, 아래로 처지는 가지가 후방이다.

실무적으로는 전방 세차가 문제다. 불평형 원심력은 항상 스핀 방향으로 회전하는 가진이므로 전방 모드만 때린다. 등방성 베어링이면 후방 모드는 불평형으로는 아예 여기되지 않는다. 후방 세차는 베어링 강성이 비등방이거나 마찰 접촉이 있을 때 등장하며, 나타나면 뭔가 잘못됐다는 신호에 가깝다.

4. 캠벨 선도[편집]

캠벨 선도(Campbell diagram)는 가로축 회전수, 세로축 진동수에 고유진동수 가지들을 그린 지도다. 회전체 동역학의 결과물 중 가장 자주 인쇄되는 한 장.

여기에 가진 차수 직선을 겹쳐 그린다. 불평형은 회전당 1회 가진이므로 기울기 1인 직선(1×선, synchronous line), 블레이드 $N$ 장의 통과 가진은 $N$ ×선, 커플링 부정렬은 보통 2×선.

가진선과 고유진동수 가지가 만나는 점이 임계속도다. 설계 규칙은 단순하다: 운전 회전수 대역에 교차점이 없도록, API 표준 기준으로 **분리 마진 15~20%**를 확보한다.² 마진이 안 나오면 베어링 강성을 바꾸거나, 스팬을 줄이거나, 디스크 위치를 옮긴다.

주의할 함정: 캠벨 선도에 그려진 임계속도라고 다 위험한 게 아니다. 그 모드가 실제로 가진되는지(모드 형상과 불평형 분포의 매칭), 감쇠가 충분한지(로그 감쇠율)를 봐야 한다. 그래서 캠벨 선도 옆에는 반드시 불평형 응답 해석이 따라붙는다 — API 610/617이 요구하는 잔류 불평형을 넣고 회전수를 스윕해 실제 진폭과 주파수 응답 해석적 증폭 계수를 뽑는 것.

5. 베어링과 불안정[편집]

로터 해석의 정확도는 대부분 베어링 모델링에서 결정된다. 유체 윤활 저널 베어링은 유막이 강성과 감쇠를 동시에 제공하며, 이를 선형화한 8개 계수로 표현한다.

\begin{pmatrix} f_x \\ f_y \end{pmatrix} = -\begin{pmatrix} k_{xx} & k_{xy} \\ k_{yx} & k_{yy} \end{pmatrix}\begin{pmatrix} x \\ y \end{pmatrix} - \begin{pmatrix} c_{xx} & c_{xy} \\ c_{yx} & c_{yy} \end{pmatrix}\begin{pmatrix} \dot{x} \\ \dot{y} \end{pmatrix}

핵심은 교차 강성 $k_{xy}, k_{yx}$ 다. 이 둘의 부호가 반대이면( $k_{xy} = -k_{yx}$ ) 그 반대칭 성분은 세차 궤도를 따라 일을 해준다 — 즉 에너지를 공급한다. 감쇠보다 이게 커지는 순간 시스템은 자려 진동(self-excited vibration)에 들어간다.

5.1. 오일 휩과 오일 휠[편집]

저널 베어링 유막은 평균적으로 축 회전수의 약 절반으로 끌려 돈다. 이 유막이 축을 대략 0.42~0.48× 회전수로 밀고 다니는 것이 오일 휠(oil whirl)이다. 회전수를 계속 올려서 이 0.5× 성분이 로터의 1차 고유진동수와 맞아떨어지면 — 대략 1차 임계속도의 2배 회전수 — 진동이 그 고유진동수에 고정된 채 회전수를 더 올려도 따라오지 않고 진폭만 커진다. 이것이 오일 휩(oil whip)이며, 방치하면 베어링이 갈려나간다.³

구별법은 명확하다. 오일 휠은 주파수가 회전수에 비례해 따라 올라가고(0.45×가 유지됨), 오일 휩은 주파수가 고유진동수에 고정된다. 폭포수 선도(waterfall plot)에서 사선이면 휠, 수직선이면 휩. 대책은 틸팅 패드 베어링 — 패드가 각자 기울어져 교차 강성을 구조적으로 제거한다. 비싸지만 확실하다.

6. 관련 문서[편집]

1919년 헨리 젭코트(Henry Homan Jeffcott)의 논문에서 왔다. 유럽 문헌에서는 1895년에 먼저 비슷한 분석을 한 August Föppl의 이름을 따 라발-젭코트 또는 Föppl 로터라고도 부른다. 이름 싸움은 회전체 동역학에서도 예외가 아니다. ↩
API 617(원심압축기) 기준. 마진이 안 나오면 보고서에 “운전 대역 내 임계속도 존재”라고 쓰게 되는데, 이 한 줄이 설계 리뷰를 통째로 뒤집는다. 그래서 다들 어떻게든 15%를 만들어낸다. ↩
오일 휩은 1924년 뉴킬크(Newkirk)와 테일러가 보고한 이래 100년 넘게 회전기계 엔지니어를 괴롭히고 있다. 진동 스펙트럼에 0.45× 성분이 뜨는 걸 보는 순간 그날 퇴근은 없다고 보면 된다. ↩