공학기초

동역학

2026-08-19 v1 2

동역학

[L1] Part 1: 질점의 동역학 (Kinematics and Kinetics of Particles) [L2] 1. 질점의 직선 운동 (Rectilinear Kinematics of Particles) [L5] * 개념: 물체의 크기를 무시하고 오직 질량만을 가진 점(Particle)으로 가정하여, 단일 직선 축 상에서 발생하는 위치, 속도, 가속도의 시간 의존적 변화를 해석하는 학문이다. [L5] * 지배 방정식 및 수식 정리: [L5] * 위치 $s(t)$에 대하여 속도 $v(t)$와 가속도 $a(t)$는 다음과 같은 미분 관계를 갖는다. $$v = \frac{ds}{dt}$$ $$a = \frac{dv}{dt} = \frac{d^2s}{dt^2}$$ [L5] * 시간 변수 $dt$를 소거하면 위치와 속도, 가속도 사이의 관계식을 도출할 수 있다. $$a \, ds = v \, dv$$ [L5] * 가속도가 일정한 등가속도 운동의 경우, 초기 조건($t=0, s=s_0, v=v_0$) 하에서 아래의 대수 방정식이 성립한다. $$v = v_0 + a_c t$$ $$s = s_0 + v_0 t + \frac{1}{2} a_c t^2$$ $$v^2 = v_0^2 + 2a_c(s - s_0)$$ [L5] * 특징: 변수 분리법을 통한 적분 시, 가속도 $a$가 시간 $t$, 위치 $s$, 속도 $v$ 중 어느 변수의 함수인지 정의하는 것이 선행되어야 한다. [L2] 2. 질점의 곡선 운동 및 좌표계 변환 (Curvilinear Motion and Coordinate Systems) [L5] * 개념: 2차원 또는 3차원 공간 상에서 변화하는 질점의 운동 궤적을 물리적 상황에 부합하는 해석 기하학적 좌표계를 도입하여 정량화하는 기법이다. [L5] * 지배 방정식 및 수식 정리: [L5] * 직교 좌표계 (Rectilinear Coordinates): $$\mathbf{r} = x\mathbf{i} + y\mathbf{j}$$ $$\mathbf{v} = \dot{x}\mathbf{i} + \dot{y}\mathbf{j}$$ $$\mathbf{a} = \ddot{x}\mathbf{i} + \ddot{y}\mathbf{j}$$ [L5] * 법선-접선 좌표계 (Normal-Tangential Coordinates): 궤적의 접선 방향 단위 벡터 $\mathbf{u}_t$와 곡률 중심을 향하는 법선 방향 단위 벡터 $\mathbf{u}_n$을 사용한다. $$\mathbf{v} = v\mathbf{u}_t$$ $$\mathbf{a} = \dot{v}\mathbf{u}_t + \frac{v^2}{\rho}\mathbf{u}_n \quad (\rho = \text{곡률 반경})$$ [L5] * 극좌표계 (Polar Coordinates): 반경 방향 단위 벡터 $\mathbf{u}_r$과 회전 방향 단위 벡터 $\mathbf{u}_\theta$를 사용한다. $$\mathbf{r} = r\mathbf{u}_r$$ $$\mathbf{v} = \dot{r}\mathbf{u}_r + r\dot{\theta}\mathbf{u}_\theta$$ $$\mathbf{a} = (\ddot{r} - r\dot{\theta}^2)\mathbf{u}_r + (r\ddot{\theta} + 2\dot{r}\dot{\theta})\mathbf{u}_\theta$$ [L5] * 특징: 궤적의 기하학적 형상이 기지(Known)일 때는 법선-접선 좌표계가 유효하며, 특정 고정점으로부터의 거리와 각도가 관측 기준일 때는 극좌표계 변환이 강제된다. [L2] 3. 뉴턴의 제2법칙과 운동 방정식 (Kinetics: Newton's Second Law) [L5] * 개념: 질점에 작용하는 외력의 합과 그로 인해 발생하는 선형 운동량의 시간 변화율 사이의 인과 관계를 정립하는 동역학의 기초 거동 법칙이다. [L5] * 지배 방정식 및 수식 정리: [L5] * 벡터 형식의 운동 방정식: $$\sum \mathbf{F} = m\mathbf{a}$$ [L5] * 각 좌표계별 성분 분해: [L5] * 직교 좌표계: $\sum F_x = m\ddot{x}, \; \sum F_y = m\ddot{y}$ [L5] * 법선-접선 좌표계: $\sum F_t = m\dot{v}, \; \sum F_n = m\frac{v^2}{\rho}$ [L5] * 극좌표계: $\sum F_r = m(\ddot{r} - r\dot{\theta}^2), \; \sum F_\theta = m(r\ddot{\theta} + 2\dot{r}\dot{\theta})$ [L5] * 특징: 가속도 성분을 결정하는 원인력(접선력, 구심력 등)을 명확히 분리하기 위해 자유물체도(Free Body Diagram)와 동역학적 가속도 성분도(Kinetic Diagram)의 기하학적 일치가 필수적이다. [L2] 4. 일과 에너지 원리 (Kinetics: Work and Energy) [L5] * 개념: 운동 방정식을 변위에 대해 공간 적분하여, 힘이 가해진 거리와 물체의 속도 변화(운동에너지) 간의 스칼라 관계를 도출하는 해석법이다. [L5] * 지배 방정식 및 수식 정리: [L5] * 힘 $\mathbf{F}$가 변위 $d\mathbf{r}$을 이동하는 동안 행한 일($U_{1\rightarrow2}$): $$U_{1\rightarrow2} = \int_{s_1}^{s_2} F_t \, ds = \int_{1}^{2} \mathbf{F} \cdot d\mathbf{r}$$ [L5] * 일-에너지 원리: $$T_1 + \sum U_{1\rightarrow2} = T_2 \quad \left(T = \frac{1}{2}mv^2\right)$$ [L5] * 보존력장(중력, 탄성력) 내에서의 역학적 에너지 보존 법칙: $$T_1 + V_1 = T_2 + V_2$$ [L5] * 중력 위치에너지: $V_g = mgy$ [L5] * 탄성 위치에너지: $V_e = \frac{1}{2}ks^2$ [L5] * 특징: 시간 변수가 배제된 상태에서 위치 변화에 따른 속도의 크기만을 추적할 때 연산의 복잡성을 최소화할 수 있다. [L2] 5. 충격량과 운동량 원리 (Kinetics: Impulse and Momentum) [L5] * 개념: 운동 방정식을 시간에 대해 적분하여, 시간 흐름에 따라 누적된 힘의 효과(충격량)와 선형 운동량 벡터의 변화를 연계하는 해석법이다. [L5] * 지배 방정식 및 수식 정리: [L5] * 선형 충격량-운동량 원리: $$m\mathbf{v}_1 + \int_{t_1}^{t_2} \mathbf{F} \, dt = m\mathbf{v}_2$$ [L5] * 외력이 없는 상태에서의 선운동량 보존 법칙: $$\sum m\mathbf{v}_1 = \sum m\mathbf{v}_2$$ [L5] * 2체 충돌 현상에서의 반발 계수($e$): $$e = \frac{v'_{B} - v'_{A}}{v_{A} - v_{B}}$$ [L5] * 특징: 아주 짧은 시간 동안 거대한 충격력이 발생하는 충돌, 폭발 메커니즘 해석에 특화되어 있으며, 충돌면의 법선(Impact Line)과 접선(Plane of Contact)에 따라 성분 보존 여부가 결정된다. [L2] 6. 질점계의 동역학 (Kinetics of a System of Particles) [L5] * 개념: 상호작용하는 다수의 질점 그룹을 하나의 시스템으로 확장하여, 계 내부의 내력(Internal Force)을 상쇄시키고 외력에 의한 질량중심(Center of Mass)의 통합적 거동을 해석하는 이론이다. [L5] * 지배 방정식 및 수식 정리: [L5] * 질량중심 위치 벡터 $\mathbf{r}_G$ 정의 (전체 질량 $m_{tot} = \sum m_i$): $$m_{tot}\mathbf{r}_G = \sum m_i \mathbf{r}_i$$ [L5] * 질점계의 운동 방정식: $$\sum \mathbf{F}_{ext} = m_{tot}\mathbf{a}_G$$ [L5] * 임의의 고정점 $O$에 대한 질점계의 각운동량 원리: $$\sum \mathbf{M}_O = \dot{\mathbf{H}}_O \quad \left(\mathbf{H}_O = \sum (\mathbf{r}_i \times m_i \mathbf{v}_i)\right)$$ [L5] * 특징: 내부 결합력이나 파편의 비산 등 계 내부의 복잡한 물리 현상과 무관하게, 전체 외력의 합은 오직 질량중심의 선형 및 각거동만을 변화시킨다. [L1] Part 2: 강체의 평면 동역학 (Planar Kinematics and Kinetics of Rigid Bodies) [L2] 7. 강체의 평면 운동학: 상대속도 해석 (Kinematics: Relative Velocity) [L5] * 개념: 크기와 형상이 고정되어 변형되지 않는 강체(Rigid Body)가 평면 상에서 병진 운동과 회전 운동을 동시에 수행할 때, 강체 내 두 점 사이의 기하학적 구속 조건을 구하는 방법이다. [L5] * 지배 방정식 및 수식 정리: [L5] * 기준점 $A$에 대한 점 $B$의 상대속도 벡터 방정식: $$\mathbf{v}_B = \mathbf{v}_A + \mathbf{v}_{B/A} = \mathbf{v}_A + \boldsymbol{\omega} \times \mathbf{r}_{B/A}$$ [L5] * 2차원 평면(x-y 평면) 회전 각속도 벡터 $\boldsymbol{\omega} = \omega \mathbf{k}$ 도입 시 크기 관계식: $$v_{B/A} = \omega r_{B/A}$$ [L5] * 특징: 강체 내의 임의의 두 점 사이의 거리는 불변하므로, 상대속도 벡터 $\mathbf{v}*{B/A}$는 항상 두 점을 잇는 위치 벡터 $\mathbf{r}*{B/A}$에 수직해야 한다. [L2] 8. 순간중심 및 상대가속도 해석 (Instantaneous Center of Zero Velocity and Relative Acceleration) [L5] * 개념: 평면 운동 중인 강체에서 특정 순간 속도가 0이 되는 점(순간중심, IC)을 정의하여 회전 운동으로 단순화하는 기법과, 속도의 미분을 통해 상대가속도 및 코리올리 가속도를 유도하는 해석법이다. [L5] * 지배 방정식 및 수식 정리: [L5] * 순간중심(IC)을 이용한 속도 산출: $$v_i = \omega r_{i/IC}$$ [L5] * 상대가속도 벡터 방정식: $$\mathbf{a}_B = \mathbf{a}_A + \mathbf{a}_{B/A} = \mathbf{a}_A + \boldsymbol{\alpha} \times \mathbf{r}_{B/A} - \omega^2 \mathbf{r}_{B/A}$$ [L5] * 회전 좌표계에서의 코리올리 가속도(Coriolis Acceleration): 슬라이더가 회전하는 링크 위를 상대 운동할 때 발생한다. $$\mathbf{a}_{cor} = 2\boldsymbol{\omega} \times \mathbf{v}_{rel}$$ [L5] * 특징: 순간중심은 가속도가 0인 지점이 아니므로 가속도 해석에는 직접 사용할 수 없으며, 가속도 해석 시에는 법선 및 접선 방향의 상대 성분을 엄격히 분리해야 한다. [L2] 9. 강체의 평면 운동역학: 힘과 가속도 (Kinetics: Force and Acceleration) [L5] * 개념: 강체에 작용하는 외력 및 모멘트 시스템과, 질량관성모멘트로 대변되는 회전 저항 특성이 결합하여 발생하는 선형 및 각가속도 거동을 규명하는 학문이다. [L5] * 지배 방정식 및 수식 정리: [L5] * 질량관성모멘트(Moment of Inertia) 정의: $$I = \int r^2 \, dm$$ [L5] * 평행축 정리(Parallel-Axis Theorem): $$I_P = I_G + md^2$$ [L5] * 평면 운동역학 방정식: $$\sum F_x = m(\a_G)_x, \quad \sum F_y = m(\a_G)_y$$ $$\sum M_G = I_G \alpha$$ [L5] * 임의의 점 $P$에 대한 모멘트 방정식: $$\sum M_P = I_G \alpha + \mathbf{r}_{G/P} \times m\mathbf{a}_G$$ [L5] * 특징: 회전 중심이 질량중심($G$)이 아니거나 가속도가 존재하는 임의의 점일 경우, 모멘트 방정식 우변에 관성력에 의한 모멘트 항이 추가적으로 발생한다. [L2] 10. 강체의 평면 운동역학: 일과 에너지 (Kinetics: Work and Energy) [L5] * 개념: 평면 거동을 하는 강체 시스템에 대하여 외력의 일, 토크의 일, 그리고 병진 및 회전 에너지를 결합하여 시스템의 에너지 수지를 해석하는 방법이다. [L5] * 지배 방정식 및 수식 정리: [L5] * 강체의 총 운동에너지($T$): $$T = \frac{1}{2}m v_G^2 + \frac{1}{2}I_G \omega^2$$ [L5] * 고정축 $O$에 대해 회전하는 경우의 운동에너지: $$T = \frac{1}{2}I_O \omega^2 \quad (\because I_O = I_G + mr_G^2)$$ [L5] * 모멘트(토크) $M$이 행한 일($U_\theta$): $$U_{\theta} = \int_{\theta_1}^{\theta_2} M \, d\theta$$ [L5] * 강체의 일-에너지 방정식: $$T_1 + \sum U_{1\rightarrow2} = T_2$$ [L5] * 특징: 바퀴가 미끄러짐 없이 구르는 순수 구름 운동(Pure Rolling)의 경우, 지면과 접촉하는 순간중심점의 접지 마찰력은 가상 변위를 발생시키지 않으므로 일을 하지 않는다. [L2] 11. 강체의 평면 운동역학: 충격량과 운동량 (Kinetics: Impulse and Momentum) [L5] * 개념: 강체 시스템의 선형 운동량 벡터와 각운동량 벡터의 시간 적분형을 적용하여, 외력 벡터 및 외력 모멘트의 유한 시간 누적 효과를 정량화하는 기법이다. [L5] * 지배 방정식 및 수식 정리: [L5] * 선형 충격량-운동량 방정식: $$m(\mathbf{v}_G)_1 + \int_{t_1}^{t_2} \sum \mathbf{F} \, dt = m(\mathbf{v}_G)_2$$ [L5] * 각충격량-각운동량 방정식 (질량중심 $G$ 기준): $$I_G \omega_1 + \int_{t_1}^{t_2} \sum M_G \, dt = I_G \omega_2$$ [L5] * 임의의 고정점 $O$ 기준 각운동량 형태: $$\mathbf{H}_O = I_G \omega \mathbf{k} + \mathbf{r}_G \times m\mathbf{v}_G = I_O \omega \mathbf{k}$$ [L5] * 특징: 편심 충돌(Eccentric Impact) 해석 시, 충돌 지점의 기하학적 반발 계수 관계식과 구조물의 각운동량 보존 법칙을 연립하여 전후 속도 필드를 결정한다. [L2] 12. 기계 진동학 기초 (Introduction to Mechanical Vibrations) [L5] * 개념: 동역학적 평형 상태로부터 이탈한 강체 또는 질점 시스템이 복원력(Restoring Force)에 의해 주기적인 거동을 반복하는 현상을 미분방정식으로 정립하는 응용 분야이다. [L5] * 지배 방정식 및 수식 정리: [L5] * 1자유도 비감쇠 자유 진동(Un-damped Free Vibration) 지배 방정식: $$m\ddot{x} + kx = 0 \implies \ddot{x} + \omega_n^2 x = 0$$ [L5] * 고유 진동수(Natural Frequency, $\omega_n$) 및 주기($\tau$): $$\omega_n = \sqrt{\frac{k}{m}} \quad \text{(rad/s)}$$ $$\tau = \frac{2\pi}{\omega_n}$$ [L5] * 일반해 표현식 (초기 변위 $x_0$, 초기 속도 $v_0$): $$x(t) = x_0 \cos(\omega_n t) + \frac{v_0}{\omega_n} \sin(\omega_n t)$$ [L5] * 특징: 시스템의 질량 등가 성분과 강성 등가 성분을 도출하여 고유 진동 특성을 파악하는 것이 진동 제어 및 공진(Resonance) 회피 설계의 핵심이다. [L1] Part 3: 동역학 응용 및 실무 활용 (Advanced Engineering Applications) [L2] 13. 다관절 기계 시스템의 링크 동역학 (Multi-body Link Dynamics) [L5] * 개념: 기하학적으로 연계된 다수의 강체(Links)로 구성된 기계 구조물(예: 건설 중장비의 붐 및 암)에 대하여, 각 조인트의 구속 조건과 작용 유압력을 기반으로 한 다자유도 역학 해석 기법이다. [L5] * 지배 방정식 및 수식 정리: [L5] * 다관절 시스템의 대표적 해석 모델인 뉴턴-오일러(Newton-Euler) 전방향-역방향 재귀 알고리즘의 기초 속도/가속도 전파식: $$\boldsymbol{\omega}_{i+1} = \boldsymbol{\omega}_i + \dot{\theta}_{i+1}\mathbf{z}_{i+1}$$ $$\mathbf{a}_{G, i} = \mathbf{a}_i + \boldsymbol{\alpha}_i \times \mathbf{r}_{G/i} + \boldsymbol{\omega}_i \times (\boldsymbol{\omega}_i \times \mathbf{r}_{G/i})$$ [L5] * 유압 실린더의 선형 압축력 $F_{cyl}$과 조인트 작용 토크 $\tau_{joint}$ 간의 가상일 원리(Principle of Virtual Work) 관계식: $$\delta W = \tau_{joint} \, \delta \theta - F_{cyl} \, \delta s = 0 \implies \tau_{joint} = F_{cyl} \frac{ds}{d\theta}$$ [L5] * 특징: 개별 강체의 동역학 방정식이 조인트 구속력(Constraint Force)에 의해 결합되므로, 시스템 전체의 자유도 확장에 따른 연립 미분방정식 행렬 구성이 요구된다. [L2] 14. 기하 구속 조건을 포함한 물리 엔진 시뮬레이션 원리 (Constraint-based Physics Engine Principles) [L5] * 개념: 컴퓨터 그래픽스 및 CAD 상에서 물체 간의 물리적 접촉, 충돌, 조인트 결합을 수치적으로 실시간 모사하기 위해 기하학적 제한 조건을 대수적 방정식으로 변환하여 연산하는 알고리즘 원리이다. [L5] * 지배 방정식 및 수식 정리: [L5] * 위치 수준의 기하 구속 조건 함수 $C(\mathbf{q}) = 0$에 대하여, 이를 시간 미분한 속도 야코비안(Jacobian) 행렬 $\mathbf{J}$의 정의: $$\dot{C}(\mathbf{q}) = \frac{\partial C}{\partial \mathbf{q}}\dot{\mathbf{q}} = \mathbf{J}\mathbf{v} = 0$$ [L5] * 구속 조건을 만족하기 위한 라그랑주 승수(Lagrange Multiplier, $\boldsymbol{\lambda}$)를 도입한 제한 조건 하의 운동 방정식: $$\mathbf{M}\dot{\mathbf{v}} = \mathbf{F}_{ext} + \mathbf{J}^T \boldsymbol{\lambda}$$ [L5] * 특징: 가상 환경 내 강체들의 관성 행렬 $\mathbf{M}$과 구속 야코비안 행렬 $\mathbf{J}$를 결합한 대수 방정식 시스템을 매 시간 스텝마다 솔버(Solver)가 수치적으로 해결함으로써 물체의 투과 방지 및 힌지 거동을 구현한다. [L2] 15. 수치해석 기법을 이용한 비선형 동역학 방정식의 이산화 (Numerical Integration of Nonlinear Dynamics) [L5] * 개념: 해석적(Analytical) 방법으로 닫힌 형태의 해(Closed-form solution)를 구할 수 없는 비선형·고차 동역학 미분방정식을 시간 영역에서 이산화(Discretization)하여 근사 가속도, 속도, 위치 플롯을 도출하는 수치 연산 기법이다. [L5] * 지배 방정식 및 수식 정리: [L5] * 동역학 미분방정식을 상태 공간 방정식(State-Space Equation)으로 변환: $$\dot{\mathbf{x}} = f(\mathbf{x}, t), \quad \mathbf{x} = \begin{bmatrix} s \\ v \end{bmatrix}$$ [L5] * 룬게-쿠타 4차 기법 (Runge-Kutta 4th Order, RK4) 이산화 알고리즘: 시간 간격 $h$에 대하여, 다음 시간 스텝의 상태 벡터 $\mathbf{x}_{n+1}$ 계산식: $$\mathbf{x}_{n+1} = \mathbf{x}_n + \frac{h}{6}(k_1 + 2k_2 + 2k_3 + k_4)$$ $$\begin{aligned} k_1 &= f(\mathbf{x}_n, t_n) \\ k_2 &= f\left(\mathbf{x}_n + \frac{h}{2}k_1, t_n + \frac{h}{2}\right) \\ k_3 &= f\left(\mathbf{x}_n + \frac{h}{2}k_2, t_n + \frac{h}{2}\right) \\ k_4 &= f(\mathbf{x}_n + hk_3, t_n + h) \end{aligned}$$ [L5] * 특징: 1차 오일러 방법(Euler's Method) 대비 오차 절단율이 $O(h^4)$로 매우 낮아, 비선형 댐핑이나 대변형을 수반하는 실제 기계 시스템의 거동 예측 시 수치적 안정성을 보장한다.