17.37 링크 유연성 모델과 탄성 변형 해석

1. 개요

링크 유연성(link flexibility)은 로봇의 물리적 구조를 구성하는 링크 자체가 강체가 아닌 탄성체로서 변형되는 현상을 말한다. 대부분의 산업용 매니퓰레이터에서는 강체 가정이 타당하지만, 우주용 긴 매니퓰레이터, 경량 협동 로봇, 고속 경량 구조 로봇, 고정밀 반도체 장비의 이송 암에서는 링크의 굽힘과 비틀림 변형이 제어 정확성에 직접적인 영향을 미친다. 본 절은 링크 유연성을 수학적으로 기술하기 위한 대표적 접근을 정리한다. 연속체 기반의 부분 미분 방정식 모델, 정상 모드 해석, 가정 모드 방법, 유한 요소 기반 이산화, 부동 좌표계(floating frame) 정식화, 그리고 강체-유연체 결합 동역학의 기본 구조를 체계적으로 제시한다.

2. 연속체 관점의 링크 유연성

링크의 기본 기하학이 가늘고 긴 보 구조라면, 오일러-베르누이(Euler-Bernoulli) 또는 티모셴코(Timoshenko) 보 이론을 기반으로 한 연속체 모델이 적합하다. 균일 보의 횡방향 변위 $w(s,t)$ 에 대한 운동 방정식은 다음과 같이 표현된다.

$\rho A\,\frac{\partial^2 w}{\partial t^2} + EI\,\frac{\partial^4 w}{\partial s^4} = f(s,t)$

여기서 $\rho$ 는 밀도, $A$ 는 단면적, $E$ 는 영률, $I$ 는 관성 모멘트, $f(s,t)$ 는 단위 길이당 분포 하중이다. 이 편미분 방정식은 무한 차원의 자유도를 가진 연속체 모델을 기술하며, 로봇의 동역학과 결합하려면 공간 이산화를 통해 유한 차원으로 축소해야 한다. 티모셴코 보는 전단 변형과 회전 관성을 포함하여 두꺼운 보의 경우에 더 정확하지만, 수학적 복잡성이 증가한다.

모드 해석과 정상 모드

보의 경계 조건이 주어졌을 때 자유 진동 해는 정상 모드(normal mode) $\phi_k(s)$ 와 고유 주파수 $\omega_k$ 의 집합으로 표현된다. 변위는 모드 좌표 $\eta_k(t)$ 의 선형 결합으로 다음과 같이 전개된다.

$w(s,t) = \sum_{k=1}^{\infty}\phi_k(s)\,\eta_k(t)$

유한 개의 모드로 근사하는 것은 가정 모드 방법(assumed modes method)으로 불리며, 저차 모드의 에너지가 지배적임을 활용한다. 로봇 매니퓰레이터의 경우 처음 2~4개의 횡방향 굽힘 모드가 대개 충분한 정확성을 제공하며, 모드 좌표는 일반화 좌표의 일부로 동역학 방정식에 통합된다.

3. 가정 모드 방법과 라그랑지언 접근

가정 모드 방법은 각 링크의 변형을 고정된 공간 형상 함수와 시간 의존 일반화 좌표의 곱으로 표현한다.

$\mathbf{u}_i(s,t) = \sum_{k=1}^{m_i}\phi_{ik}(s)\,\eta_{ik}(t)$

링크의 운동에너지와 변형에너지는 일반화 좌표 $\eta_{ik}$ 에 대해 이차 형식으로 표현되며, 오일러-라그랑주 방정식을 적용하면 강체 좌표 $q$ 와 탄성 좌표 $\eta$ 가 결합된 확장 동역학 방정식이 얻어진다.

$\begin{bmatrix}M_{qq} & M_{q\eta}\\ M_{\eta q} & M_{\eta\eta}\end{bmatrix}\begin{bmatrix}\ddot{q}\\ \ddot{\eta}\end{bmatrix} + \begin{bmatrix}C_{qq}\\ C_{\eta q}\end{bmatrix} + \begin{bmatrix}g_q\\ K_\eta \eta\end{bmatrix} = \begin{bmatrix}\tau\\ 0\end{bmatrix}$

여기서 $K_\eta$ 는 탄성 변형의 복원력을 기술하는 강성 행렬이다. 이 형식은 강체-유연체 결합 동역학의 전형적 구조를 보여준다.

4. 부동 좌표계 정식화

유연 링크의 대변형을 다루는 데는 부동 좌표계(floating frame of reference) 방법이 널리 사용된다. 이 접근은 각 링크에 강체 운동을 기술하는 부동 좌표계를 부여하고, 링크의 탄성 변형을 이 좌표계 위에서의 작은 변형으로 표현한다. 이로써 강체 운동과 탄성 변형의 수학적 분리가 가능해지며, 대회전을 포함하는 매니퓰레이터에서도 탄성 모델이 직관적으로 결합된다. 부동 좌표계 방법은 ANCF(absolute nodal coordinate formulation)와 함께 현대 다체 동역학 소프트웨어에서 유연체 해석의 표준 틀을 이룬다.

5. 유한 요소 기반 이산화

보다 일반적인 기하의 링크에는 유한 요소 방법(FEM)을 이용한 공간 이산화가 사용된다. 각 요소에 대해 변위와 회전의 절점 값으로 정의된 형상 함수가 도입되며, 전체 링크의 운동에너지와 변형에너지가 절점 좌표 함수의 이차 형식으로 표현된다. 이로써 다음과 같은 표준 형식이 얻어진다.

$\mathbf{M}_e\,\ddot{\mathbf{u}} + \mathbf{K}_e\,\mathbf{u} = \mathbf{f}_e$

여기서 $\mathbf{u}$ 는 절점 좌표 벡터이다. 유한 요소 모델은 자유도 수가 크므로, 동역학 통합 전에 모드 축소(Craig-Bampton 방법 등)를 통해 차원을 줄이는 것이 일반적이다.

강체-유연체 결합 동역학

유연 링크를 포함하는 매니퓰레이터의 전체 동역학은 강체 좌표, 탄성 좌표, 구동기 좌표가 결합된 확장 방정식으로 기술된다. 이 방정식은 관성 결합 항을 통해 강체 운동이 탄성 변형을 유도하고, 탄성 복원력이 강체 운동에 반작용하는 상호 의존적 구조를 형성한다. 특히 고속 또는 경량 매니퓰레이터에서는 이러한 상호 작용이 두드러지며, 엔드 이펙터의 실제 위치는 강체 모델이 예측하는 위치와 탄성 변형에 의한 편차의 합으로 나타난다.

탄성 변형의 경계 조건

링크의 경계 조건은 링크가 관절에 어떻게 고정되는지와 부하가 어떻게 적용되는지에 따라 결정된다. 전형적인 경계 조건은 다음과 같다. 첫째, 고정-자유 보(cantilever beam)는 한쪽 끝이 관절에 고정되고 다른 쪽 끝이 자유롭거나 부하를 받는 형태이다. 둘째, 고정-고정 보는 양 끝이 모두 관절로 구속된 경우이다. 셋째, 감쇠된 자유 끝 조건은 엔드 이펙터의 동적 질량을 고려한 경우이다. 경계 조건은 정상 모드와 고유 주파수를 결정하므로, 정확한 조건 설정이 해석 결과의 타당성을 좌우한다.

탄성 매개변수와 감쇠

링크의 탄성 특성은 영률, 관성 모멘트, 단면 기하, 재료 밀도 등으로 결정되며, 이 매개변수들은 CAD 설계에서 비교적 정확히 알려진다. 반면 감쇠는 재료 내부 마찰, 구조적 감쇠, 공기 저항 등에 의해 발생하며, 정확한 모델링이 어려워 비례 감쇠 행렬 $\mathbf{C} = \alpha\mathbf{M} + \beta\mathbf{K}$ 와 같은 레일리(Rayleigh) 형식으로 근사된다. 이 근사는 각 모드에 대해 일정한 감쇠 비율을 부여하며, 실험적 식별을 통해 $\alpha, \beta$ 가 결정된다.

관측성과 제어 가능성

유연 링크 시스템의 저차 모드는 일반적으로 관측 가능하고 제어 가능하지만, 고차 모드는 제어 대역폭을 초과하거나 측정 센서의 분해능 이하로 나타나기 때문에 실질적으로 관측되지 않는다. 이러한 미모델링 고차 모드는 제어 시스템에서 스필오버(spillover) 현상을 유발할 수 있으며, 제어기가 의도하지 않은 고차 진동을 여기시킬 수 있다. 이를 방지하기 위해 제어 설계는 저역 통과 롤오프와 고주파 차단 필터를 포함한다.

본 절의 의의

본 절은 링크 유연성의 수학적 기술에 필요한 주요 도구를 체계적으로 제시한다. 연속체 방정식, 정상 모드 해석, 가정 모드 방법, 유한 요소 이산화, 부동 좌표계 접근, 강체-유연체 결합 동역학, 경계 조건의 선택, 감쇠 모델링, 관측 및 제어 가능성은 모두 유연 매니퓰레이터의 정확한 해석과 제어 설계를 위해 요구되는 요소들이다. 이러한 기초는 이후 절의 유연 시스템 진동 해석과 대응 제어 기법의 이론적 출발점이 된다.

학습 권장사항

독자는 단일 유연 링크에 대해 오일러-베르누이 보 방정식의 정상 모드를 해석적으로 구하고, 첫 두 모드의 고유 주파수와 형상 함수를 수치적으로 그려 볼 것을 권장한다. 또한 가정 모드 방법과 유한 요소 방법으로 각각 이산화한 모델의 응답을 비교하고, 차원 축소가 정확성에 미치는 영향을 분석하는 실습이 유익하다. 부동 좌표계 방법을 2자유도 유연 매니퓰레이터에 적용하여 대회전과 탄성 변형이 결합된 시뮬레이션을 구성해 보는 것도 권장된다.

참고 문헌

Shabana, A. A. (2020). Dynamics of Multibody Systems (5th ed.). Cambridge University Press.
Book, W. J. (1984). Recursive Lagrangian dynamics of flexible manipulator arms. International Journal of Robotics Research, 3(3), 87–101.
De Luca, A., & Siciliano, B. (1991). Closed-form dynamic model of planar multilink lightweight robots. IEEE Transactions on Systems, Man, and Cybernetics, 21(4), 826–839.
Meirovitch, L. (2001). Fundamentals of Vibrations. McGraw-Hill.
Craig, R. R., & Bampton, M. C. C. (1968). Coupling of substructures for dynamic analyses. AIAA Journal, 6(7), 1313–1319.
Dwivedy, S. K., & Eberhard, P. (2006). Dynamic analysis of flexible manipulators, a literature review. Mechanism and Machine Theory, 41(7), 749–777.
Shabana, A. A. (1997). Definition of the slopes and the finite element absolute nodal coordinate formulation. Multibody System Dynamics, 1(3), 339–348.

version: 1.0