7.45 감쇠 진동과 로봇 관절 응답 해석

1. 감쇠 진동의 지배 방정식

로봇 관절의 감쇠 진동을 기술하는 기본 방정식은 2차 선형 상미분 방정식이다. 관절 각도 $\theta(t)$ 에 대한 자유 진동(free vibration) 방정식은 다음과 같다.

$J\ddot{\theta} + b\dot{\theta} + k\theta = 0$

여기서 $J$ 는 관절 축에 대한 등가 관성 모멘트, $b$ 는 점성 감쇠 계수, $k$ 는 등가 강성 계수이다. 이 방정식을 표준 형태로 정규화하면 다음과 같다.

$\ddot{\theta} + 2\zeta\omega_n\dot{\theta} + \omega_n^2\theta = 0$

여기서 두 매개변수는 시스템의 진동 특성을 완전히 결정한다.

$\omega_n = \sqrt{\frac{k}{J}} \quad \text{(비감쇠 고유 진동수)}$

$\zeta = \frac{b}{2\sqrt{kJ}} = \frac{b}{2J\omega_n} \quad \text{(감쇠비)}$

2. 특성 방정식과 근의 분류

표준형 방정식의 특성 방정식은 다음과 같다.

$s^2 + 2\zeta\omega_n s + \omega_n^2 = 0$

이 이차 방정식의 근은 다음과 같다.

$s_{1,2} = -\zeta\omega_n \pm \omega_n\sqrt{\zeta^2 - 1}$

감쇠비 $\zeta$ 의 값에 따라 근의 성격이 변하며, 이에 따라 시스템의 응답 특성이 근본적으로 달라진다.

3. 과감쇠 응답

3.1 조건과 해의 형태

감쇠비가 $\zeta > 1$ 인 경우를 과감쇠(overdamped)라 한다. 특성 방정식의 두 근은 모두 실수이며 음수이다.

$s_1 = -\zeta\omega_n + \omega_n\sqrt{\zeta^2 - 1}, \quad s_2 = -\zeta\omega_n - \omega_n\sqrt{\zeta^2 - 1}$

일반해는 다음과 같다.

$\theta(t) = c_1 e^{s_1 t} + c_2 e^{s_2 t}$

여기서 $\lvert s_2 \rvert > \lvert s_1 \rvert$ 이므로, $e^{s_2 t}$ 성분이 더 빠르게 소멸한다. 충분한 시간이 경과하면 응답은 느린 모드 $e^{s_1 t}$ 에 의해 지배된다.

3.2 초기 조건에 따른 해

초기 조건 $\theta(0) = \theta_0$ , $\dot{\theta}(0) = \dot{\theta}_0$ 에 대하여 상수는 다음과 같이 결정된다.

$c_1 = \frac{\dot{\theta}_0 - s_2\theta_0}{s_1 - s_2}, \quad c_2 = \frac{s_1\theta_0 - \dot{\theta}_0}{s_1 - s_2}$

과감쇠 시스템은 진동 없이 평형점에 단조적으로 접근하나, 복귀 속도가 느리다. 로봇 관절에서 과감쇠는 위치 결정 속도를 저하시키므로, 일반적으로 바람직하지 않다.

4. 임계 감쇠 응답

4.1 조건과 해의 형태

감쇠비가 정확히 $\zeta = 1$ 인 경우를 임계 감쇠(critically damped)라 한다. 특성 방정식은 중근 $s = -\omega_n$ 을 갖는다. 일반해는 다음과 같다.

$\theta(t) = (c_1 + c_2 t)e^{-\omega_n t}$

초기 조건 $\theta(0) = \theta_0$ , $\dot{\theta}(0) = \dot{\theta}_0$ 에 대하여 다음을 얻는다.

$\theta(t) = \left[\theta_0 + (\dot{\theta}_0 + \omega_n\theta_0)t\right]e^{-\omega_n t}$

4.2 최속 비진동 복귀

임계 감쇠는 진동 없이 평형점으로 복귀하는 응답 중 가장 빠른 복귀를 실현한다. 임계 감쇠 계수는 $b_{cr} = 2\sqrt{kJ} = 2J\omega_n$ 이다.

로봇 관절 제어에서 임계 감쇠는 오버슈트 없이 가능한 한 빠르게 목표 위치에 도달하여야 하는 경우에 이상적인 감쇠 조건이다. 그러나 실제 시스템에서는 매개변수 불확실성으로 인하여 정확한 임계 감쇠를 달성하기 어렵다.

5. 부족 감쇠 응답

5.1 조건과 해의 형태

감쇠비가 $0 < \zeta < 1$ 인 경우를 부족 감쇠(underdamped)라 한다. 특성 방정식의 근은 복소 켤레쌍이다.

$s_{1,2} = -\zeta\omega_n \pm j\omega_d$

여기서 $j = \sqrt{-1}$ 이고, 감쇠 고유 진동수(damped natural frequency)는 다음과 같다.

$\omega_d = \omega_n\sqrt{1 - \zeta^2}$

일반해는 감쇠하는 정현파 형태이다.

$\theta(t) = e^{-\zeta\omega_n t}\left(A\cos\omega_d t + B\sin\omega_d t\right)$

이를 진폭-위상 형태로 다시 쓰면 다음과 같다.

$\theta(t) = Ce^{-\zeta\omega_n t}\cos(\omega_d t - \phi)$

여기서 진폭 $C$ 와 위상 $\phi$ 는 초기 조건에 의해 결정된다.

$C = \sqrt{A^2 + B^2}, \quad \phi = \arctan\left(\frac{B}{A}\right)$

초기 조건 $\theta(0) = \theta_0$ , $\dot{\theta}(0) = \dot{\theta}_0$ 에 대하여 다음이 성립한다.

$A = \theta_0, \quad B = \frac{\dot{\theta}_0 + \zeta\omega_n\theta_0}{\omega_d}$

5.2 포락선과 감쇠율

응답의 진폭은 포락선(envelope) $\pm Ce^{-\zeta\omega_n t}$ 에 의해 지수적으로 감쇠한다. 시간 상수는 $\tau = 1/(\zeta\omega_n)$ 이며, 이는 진폭이 초기값의 $1/e \approx 36.8\%$ 로 감소하는 데 소요되는 시간이다. 실용적으로 정착 시간(settling time)은 시간 상수의 배수로 정의된다.

6. 과도 응답의 성능 지표

로봇 관절의 계단 응답(step response)에서 과도 응답 성능을 정량화하는 주요 지표는 다음과 같다.

6.1 상승 시간

상승 시간(rise time) $t_r$ 은 응답이 최종값의 0%에서 100%까지(또는 10%에서 90%까지) 도달하는 데 소요되는 시간이다. 부족 감쇠 시스템에서 0-100% 상승 시간은 다음과 같이 근사된다.

$t_r = \frac{\pi - \phi}{\omega_d}$

여기서 $\phi = \arctan(\omega_d / (\zeta\omega_n))$ 이다.

6.2 최대 오버슈트

최대 오버슈트(maximum overshoot) $M_p$ 는 응답이 최종값을 초과하는 최대 비율이다.

$M_p = e^{-\pi\zeta / \sqrt{1 - \zeta^2}} \times 100\%$

이 공식에서 오버슈트는 감쇠비 $\zeta$ 에만 의존한다. 오버슈트를 5% 이하로 제한하려면 $\zeta \geq 0.69$ 이어야 한다.

6.3 최대 오버슈트 시간

최대 오버슈트에 도달하는 시간(peak time) $t_p$ 는 다음과 같다.

$t_p = \frac{\pi}{\omega_d} = \frac{\pi}{\omega_n\sqrt{1 - \zeta^2}}$

6.4 정착 시간

정착 시간(settling time) $t_s$ 는 응답이 최종값의 허용 오차 범위 내에 진입하여 더 이상 벗어나지 않는 데 소요되는 시간이다. 2% 기준과 5% 기준의 근사 공식은 각각 다음과 같다.

$t_s \approx \frac{4}{\zeta\omega_n} \quad \text{(2\% 기준)}, \quad t_s \approx \frac{3}{\zeta\omega_n} \quad \text{(5\% 기준)}$

정착 시간은 $\zeta\omega_n$ 의 곱에 반비례하므로, 빠른 정착을 위해서는 감쇠비와 고유 진동수 모두를 높여야 한다.

7. 감쇠비에 따른 응답 특성 비교

다음 표는 감쇠비에 따른 시스템 응답 특성을 요약한 것이다.

감쇠비 $\zeta$	분류	응답 형태	오버슈트	복귀 속도
$\zeta = 0$	비감쇠	지속 정현파 진동	100%	복귀하지 않음
$0 < \zeta < 1$	부족 감쇠	감쇠 진동	존재	중간
$\zeta = 1$	임계 감쇠	비진동 감쇠	없음	최속 비진동
$\zeta > 1$	과감쇠	비진동 감쇠	없음	느림

8. 강제 진동과 정상 상태 응답

8.1 정현파 강제 입력

정현파 토크 입력 $\tau(t) = \tau_0 \sin(\omega t)$ 가 인가되는 경우, 정상 상태 응답은 다음과 같다.

$\theta_{ss}(t) = \frac{\tau_0 / k}{\sqrt{(1 - r^2)^2 + (2\zeta r)^2}} \sin(\omega t - \phi)$

여기서 주파수비 $r = \omega / \omega_n$ 이고, 위상 지연은 다음과 같다.

$\phi = \arctan\left(\frac{2\zeta r}{1 - r^2}\right)$

8.2 진폭비와 공진

정상 상태 진폭을 정적 변위 $\theta_{st} = \tau_0 / k$ 로 정규화한 진폭비(magnification factor) $M$ 은 다음과 같다.

$M(r) = \frac{1}{\sqrt{(1 - r^2)^2 + (2\zeta r)^2}}$

공진 주파수(resonance frequency)는 진폭비가 최대가 되는 주파수이며, $\zeta < 1/\sqrt{2}$ 일 때 다음과 같다.

$\omega_r = \omega_n\sqrt{1 - 2\zeta^2}$

공진 시 진폭비의 최댓값은 다음과 같다.

$M_{max} = \frac{1}{2\zeta\sqrt{1 - \zeta^2}}$

감쇠비가 작을수록 공진 피크가 뾰족해지며 진폭이 급격히 증가한다. 로봇 구조물과 제어 시스템에서 공진은 구조적 손상이나 제어 불안정을 초래할 수 있으므로, 작동 주파수 범위에서 공진을 회피하는 설계가 필수적이다.

9. 로봇 관절 응답 해석의 실제적 고려 사항

9.1 PD 제어와 등가 감쇠 시스템

PD(비례-미분) 제어기를 적용한 관절 제어 시스템은 감쇠 진동 모델과 직접적으로 대응된다. 제어 입력 $\tau = -K_p(\theta - \theta_d) - K_d\dot{\theta}$ 를 관절 방정식에 대입하면 다음을 얻는다.

$J\ddot{\theta} + (b + K_d)\dot{\theta} + (k + K_p)\theta = K_p\theta_d$

이 방정식에서 등가 감쇠비와 등가 고유 진동수는 다음과 같다.

$\omega_{n,eq} = \sqrt{\frac{k + K_p}{J}}, \quad \zeta_{eq} = \frac{b + K_d}{2\sqrt{(k + K_p)J}}$

따라서 제어 이득 $K_p$ 와 $K_d$ 를 조정함으로써 원하는 감쇠비와 고유 진동수를 설정할 수 있다. 이는 극배치(pole placement)의 관점에서 폐루프 시스템의 극점을 원하는 위치에 배치하는 것과 동일하다.

9.2 감쇠비 설계 지침

로봇 관절 제어에서 감쇠비의 선택은 오버슈트와 응답 속도 사이의 상충 관계(trade-off)를 수반한다. 일반적으로 다음의 설계 지침이 적용된다.

$\zeta = 0.4 \sim 0.8$ : 적절한 오버슈트와 빠른 응답을 동시에 달성할 수 있는 범위이다.
$\zeta \approx 0.707$ ( $= 1/\sqrt{2}$ ): ITAE(Integral of Time-weighted Absolute Error) 기준에서 최적에 가까운 감쇠비로, 오버슈트 약 4.3%와 합리적인 정착 시간을 제공한다.
$\zeta = 1$ : 오버슈트가 절대적으로 허용되지 않는 안전 관련 작업에 적용한다.

9.3 대수 감쇠율

실험적으로 감쇠비를 측정하는 방법으로 대수 감쇠율(logarithmic decrement) $\delta$ 를 이용한다. 연속하는 두 피크의 비로부터 $\delta$ 를 구할 수 있다.

$\delta = \ln\left(\frac{\theta(t_k)}{\theta(t_{k+1})}\right) = \frac{2\pi\zeta}{\sqrt{1 - \zeta^2}}$

여기서 $t_k$ 와 $t_{k+1}$ 은 연속하는 두 피크의 시각이다. $N$ 개의 주기에 걸쳐 측정하면 정밀도가 향상된다.

$\delta = \frac{1}{N}\ln\left(\frac{\theta(t_k)}{\theta(t_{k+N})}\right)$

감쇠비는 대수 감쇠율로부터 다음과 같이 역산된다.

$\zeta = \frac{\delta}{\sqrt{4\pi^2 + \delta^2}}$

이 방법은 로봇 관절의 마찰 특성 식별 및 제어 이득 튜닝에 실용적으로 활용된다.

10. 요약

감쇠 진동 해석은 로봇 관절의 과도 응답 특성을 이해하고 제어기를 설계하는 핵심 도구이다. 감쇠비에 따른 과감쇠, 임계 감쇠, 부족 감쇠의 세 가지 응답 유형은 각각 고유한 동적 특성을 가지며, 설계 요구 사항에 따라 적절한 감쇠비를 선택하여야 한다. 상승 시간, 오버슈트, 정착 시간 등의 성능 지표는 감쇠비와 고유 진동수로 표현되므로, 이 두 매개변수를 제어 이득으로 조정함으로써 원하는 과도 응답 사양을 달성할 수 있다. 공진 현상의 이해와 회피는 로봇의 구조적 안전성과 제어 안정성을 보장하기 위하여 반드시 필요하다.

참고 문헌

Ogata, K. (2010). Modern Control Engineering (5th ed.). Prentice Hall.
Rao, S. S. (2017). Mechanical Vibrations (6th ed.). Pearson.
Craig, J. J. (2005). Introduction to Robotics: Mechanics and Control (3rd ed.). Pearson.
Spong, M. W., Hutchinson, S., & Vidyasagar, M. (2006). Robot Modeling and Control. Wiley.
Nise, N. S. (2015). Control Systems Engineering (7th ed.). Wiley.
Franklin, G. F., Powell, J. D., & Emami-Naeini, A. (2015). Feedback Control of Dynamic Systems (7th ed.). Pearson.

version: 0.1.0