항공우주/동역학29 수정 로드리게스 파라미터 (MRP) 수정 로드리게스 파라미터(MRP, modified Rodrigues parameters)는 1962년에 T. F. Wiener에 의해서 고안되었다. MRP \(\mathbf{e}_b^a\) 의 정의는 다음과 같다. \[ \begin{align} \mathbf{e}_b^a= \frac{ \mathbf{q}_{1:3}}{1+q_0 } \tag{1} \end{align} \] 앞선 게시글 (https://pasus.tistory.com/377) 에 있는 고전(classical) 로드리게스 파라미터의 정의와 비교해 보면 분모에 1을 더한 것을 볼 수 있다. 식 (1)과 쿼터니언의 특성에 의하면 다음식이 성립한다. \[ \begin{align} (1+q_0 )^2 (\mathbf{e}_b^a )^T \ma.. 2025. 6. 3. 로드리게스 파라미터 로드리게스 파라미터(Rodrigues parameters) 벡터 또는 깁스(Gibbs) 벡터는 1840년 로드리게스(Olinde Rodrigues)에 의해서 처음 도입되었다. 로드리게스 파라미터 벡터 \(\mathbf{g}_b^a\)는 쿼터니언으로부터 다음과 같이 정의된다. \[ \begin{align} \mathbf{g}_b^a= \begin{bmatrix} g_1 \\ g_2 \\ g_3 \end{bmatrix} = \frac{ \mathbf{q}_{1:3} }{ q_0} \tag{1} \end{align} \] 여기서 \(\mathbf{q}_{1:3}\) 는 쿼터니언의 벡터부이고 \(q_0\) 는 스칼라부다. 참고로 좌표계 \(\{a\}\) 를 회전축 \(\hat{p}\) 축을 중심으로 \(\.. 2025. 6. 1. 회전벡터 (Rotation Vector)의 시간 변화율 모든 좌표변환은 어떤 회전축과 그 회전축을 중심으로 하는 한번의 회전을 통해서 가능하다. 쿼터니언(quaternions)의 정의도 이 회전축과 회전각에 기반을 두고 있다. 단위벡터는 크기가 1인 벡터이기 때문에 방향을 표시하는데 자주 쓰인다. 여기서도 회전축 방향을 정하는데 단위벡터를 이용하기로 하고 기호로 \(\hat{p}\) 으로 표시한다. 회전벡터(rotation vector)는 회전축과 회전각을 간결하게 표현한 벡터로서 그 크기는 회전각이며 방향은 회전축과 같은 방향으로서 다음과 같이 정의한다. \[ \begin{align}\vec{\beta} = \beta \hat{p} \tag{1} \end{align} \] 여기서 \(\beta\) 는 회전각, \(\vec{\beta}\) 는 .. 2025. 4. 9. 오일러 운동방정식 (Euler’s Equation of Motion) 질량중심을 기준으로 한 강체의 운동방정식은 다음과 같았다 (https://pasus.tistory.com/191). \[ \begin{align} \vec{M}_G &= \frac{ ^bd \vec{H}_G}{dt}+ \ ^i\vec{\omega} ^b \times \vec{H}_G \tag{1} \\ \\ &= \bar{I}_G \cdot \frac{ ^b d \ ^i\vec{\omega} ^b }{dt} + \ ^i\vec{\omega} ^b \times (\bar{I}_G \cdot \ ^i\vec{\omega} ^b ) \end{align} \] 여기서 \( ^i \vec{\omega} ^b\) 는 관성 좌표계 \(\{i\}\) 에 대한 강체 좌표계 \(\{b\}\) 의 각속도벡터, \(\b.. 2024. 3. 22. 관성 주축 (Principal Axes of Inertia) 강체에 고정된 임의의 점 B를 원점으로 하고 강체에 고정된 좌표계 \(\{a\}\) 에 관한 관성행렬(inertia matrix)은 다음 식으로 계산할 수 있다 (https://pasus.tistory.com/191). \[ [I_B^a ]= \begin{bmatrix} I_{xx} & I_{xy} & I_{xz} \\ I_{yx} & I_{yy} & I_{yz} \\ I_{zx} & I_{zy} & I_{zz} \end{bmatrix} \tag{1} \] 여기서 \(x, y, z\) 는 좌표계 \(\{a\}\) 의 원점에서 강체 내의 임의의 점까지의 위치 좌표이고, 즉 \(\vec{r}=x\hat{a}_1+y\hat{a}_2+z\hat{a}_z\), 행렬 \([I_B^a ]\) 의 구성 성분은 다음과 .. 2023. 2. 19. 좌표변환과 관성행렬 (Inertia Matrix) 관성 다이아딕(inertia dyadic)은 특정 좌표계와 무관하지만 관성 다이아딕을 특정 좌표계로 표현한 관성행렬(inertia matrix)은 좌표계에 따라 달라진다. 어떤 강체의 질량중심 \(G\) 를 원점으로 하고 강체에 고정된 좌표계 \(\{a\}\) 와 좌표계 \(\{b\}\) 가 있다고 하자. 질량중심 \(G\) 에 관한 관성 다이아딕 \(\bar{I}_G\) 를 좌표계 \(\{a\}\) 와 좌표계 \(\{b\}\) 로 각각 표현하면 다음과 같다 (https://pasus.tistory.com/191). \[ \begin{align} \bar{I}_G &= \sum_{i=1}^3 \sum_{j=1}^3 I_{ij}^a \ \hat{a}_i \hat{a}_j \tag{1} \\ \\ &= \s.. 2023. 2. 17. 오일러의 회전 정리 (Euler’s Rotation Theorem) 오일러각 좌표변환 방법에서 알아본 회전축은 좌표계의 \(x\) 축, \(y\) 축, \(z\) 축이었다. 하지만 좌표계를 구성하는 좌표축만이 아니라 임의의 축, 즉 임의의 방향을 중심으로 좌표계를 회전시킬 수도 있다. 단위벡터는 크기가 \(1\) 인 벡터이기 때문에 방향을 표시하는데 자주 쓰인다. 여기서도 회전축 방향을 정하는데 단위벡터를 이용하기로 하고 기호로 \(\hat{p}\) 으로 표시하기로 한다. 좌표계 \(\{a\}\) 를 회전축 \(\hat{p}\) 축을 중심으로 \(\beta\) 만큼 회전시키면 새로운 좌표계로 변환되는데 이 좌표계를 \(\{b\}\) 라고 하자. 그러면 그림에서 보듯이 좌표계 \(\{a\}\) 의 좌표축과 회전축 사이의 각도는 좌표계 \(\{b\}\) 의 좌표축과 회전축 .. 2022. 3. 22. 라그랑지 방정식을 이용한 강체 운동방정식 유도 강체(rigid body)의 다양한 지점에 가해지는 모든 외력(external force)은 질량중심(center of mass)에 가해지는 총 외력으로 합산할 수 있고 질량중심은 마치 강체의 모든 질량이 그 중심에 집중되어 있는 질점(point mass)처럼 운동한다. 또한 외력은 강체의 다양한 지점에서 작용하기 때문에 질량중심에 대해서 모멘트를 만들고 이 모멘트는 질량중심에 대한 회전운동을 생성한다. 이와 같이 강체의 운동은 질량중심의 병진운동과 질량중심에 대한 회전운동으로 분리할 수 있다. 이제 강체 운동방정식을 라그랑지 방정식(Lagrange's Equation)을 이용하여 유도해 보도록 하겠다. 강체의 운동에너지도 질량중심의 병진 운동에너지와 질량중심에 대한 회전 운동에너지의 합으로 표현할 수 .. 2022. 2. 14. 강체의 운동방정식 - 4 지금까지 질량중심을 기준으로 강체(rigid body)의 운동방정식을 유도하였다. 이번에는 강체에 고정되어 있는 임의의 점 \(A\) 에 대해서 강체의 운동방정식을 유도해 보도록 하겠다. 임의의 점 \(A\) 에 대한 파티클 시스템(systems of particles)의 운동방정식은 다음과 같았다. \[ \begin{align} & \sum_{j=1}^n \vec{F}_j = m \frac{^i d^2 \vec{r}_G}{dt^2} = m \frac{^i d \vec{v}_G }{dt} \tag{1} \\ \\ & \frac{^i d \vec{H}_A}{dt} = m \frac{^i d \vec{r}_{G/A}}{dt} \times \vec{v}_G + \sum_{j=1}^n \vec{M}_{jA} \.. 2022. 2. 7. 강체의 운동방정식 - 3 지금까지 파티클 시스템(systems of particles)에 대해서 다음과 같은 운동방정식을 얻었다. \[ \begin{align} & \sum_{j=1}^n \vec{F}_j =m \frac{^id^2 \vec{r}_G }{dt^2}= m \frac{^id \vec{v}_G }{dt} \tag{1} \\ \\ & \sum_{j=1}^n \vec{M}_{jG} = \frac{^id \vec{H}_G }{dt} \tag{2} \\ \\ & \vec{H}_G= \sum_{j=1}^n \vec{r}_{j/G} \times m_j \frac{^id \vec{r}_j}{dt} \\ \\ & T= \frac{1}{2} m \vec{v}_G \cdot \vec{v}_G + \frac{1}{2} \sum_{j=1}.. 2022. 2. 6. 이전 1 2 3 다음
