항공우주102 [CR3BP] 리야프노프 궤도, 헤일로 궤도, 그리고 리사주 궤도 라그랑지 포인트 L1, L2 및 L3에서의 선형화 운동방정식은 다음과 같았다 (https://pasus.tistory.com/272). \[ \begin{align} & \delta \ddot{x}-2 \delta \dot{y}-(1+2c_2 ) \delta x=0 \tag{1} \\ \\ & \delta \ddot{y}+2 \delta \dot{x}+(-1+c_2 ) \delta y=0 \\ \\ & \delta \ddot{z}+c_2 \delta z=0 \end{align} \] 여기서 \[ c_2= \frac{(1-\mu)}{|x_0+\mu|^3 }+ \frac{\mu}{ |x_0+\mu-1|^3 } \tag{2} \] 이다. 식 (1)에서 \(\delta x, \ \delta y\) 운동을 벡터 .. 2023. 6. 27. [CR3BP] L1, L2 및 L3 포인트에서의 궤도 운동 CR3BP의 선형화된 운동방정식을 이용하여 라그랑지 포인트(Lagrange point) L4 및 L5 포인트는 (중립) 안정 평형점이지만, L1, L2 및 L3 포인트는 불안정한 평형점이라는 것을 확인했다 (https://pasus.tistory.com/271). 하지만 L1, L2 및 L3 포인트의 고유값(eigenvalue) 분석에 의하면 평형점 주위에 주기 궤도(periodic orbit)가 존재함을 시사한다. 즉 특정한 초기조건을 설정하면 불안정한 운동 모드를 배제하고 주기 운동을 하는 모드만을 나타나게 할 수가 있다. 라그랑지 포인트에서의 선형화 운동방정식은 다음과 같다. \[ \begin{align} & \delta \ddot{x}-2 \delta \dot{y} = -\bar{U}_{xx} \.. 2023. 6. 25. [CR3BP] 라그랑지 포인트 안정성 해석 CR3BP의 무차원화된 운동방정식은 다음과 같았다 (https://pasus.tistory.com/147). \[ \begin{align} & \ddot{x}-2 \dot{y} = -\bar{U}_x \tag{1} \\ \\ & \ddot{y}+2 \dot{x} = -\bar{U}_y \\ \\ & \ddot{z} = -\bar{U}_z \end{align} \] 여기서 \[ \begin{align} & U_{eff}= -\frac{1}{2} (x^2+y^2 ) - \frac{1-\mu}{r_1} - \frac{\mu}{r_2} - \frac{1}{2} \mu (1-\mu) \\ \\ & r_1= \sqrt{(x+\mu)^2+y^2+z^2 } \\ \\ & r_2= \sqrt{(x+\mu-1)^2+y^2.. 2023. 6. 22. [CR3BP] 힐의 영역 (Hill’s Region) 원궤도 제한 삼체문제(CR3BP)는 질량중심을 중심으로 원궤도 운동을 하는 두 개의 기본 질점에 의해 생성된 중력장에서 제3의 질점의 운동을 기술한다. CR3BP의 무차원화된 운동방정식은 다음과 같았다 (https://pasus.tistory.com/147). \[ \begin{align} & \ddot{x}-2\dot{y}- x= - \frac{(1-\mu)(x+\mu) }{r_1^3 }- \frac{\mu (x+\mu-1)}{ r_2^3 } \tag{1} \\ \\ & \ddot{y}+2 \dot{x}-y= - \frac{(1-\mu)y}{r_1^3 }- \frac{\mu y}{ r_2^3 } \\ \\ & \ddot{z}=- \frac{(1-\mu )z}{r_1^3 }- \frac{\mu z}{ r.. 2023. 6. 19. 항공기의 질점 (Point Mass) 운동 모델 전투기 교전 운동 모델, 비행기 성능(performance) 해석 모델, 그리고 단거리 미사일 운동 모델로서 다음과 같이 평평한 지구(flat Earth) 가정 하에서 유도된 질점(point mass) 운동 모델을 많이 사용한다 (https://pasus.tistory.com/181). \[ \begin{align} & \dot{x} =V \cos \psi \cos \gamma \tag{1} \\ \\ & \dot{y}=V \sin \psi \cos \gamma \\ \\ & \dot{h}=V \sin \gamma \\ \\ & \dot{V}= -\frac{D}{m}+ \frac{T \cos \epsilon }{m} -g \sin \gamma \\ \\ & \dot{\psi} = \frac{ (L+T \.. 2023. 3. 12. 상대 궤도요소의 섭동 (Perturbed Relative Orbital Elements) Clohessy-Wiltshire(CW) 방정식을 \[ \begin{align} & \ddot{x}-3n^2 x-2n \dot{y}=f_1 \tag{1} \\ \\ & \ddot{y}+2n \dot{x}=f_2 \\ \\ & \ddot{z}+n^2 z=f_3 \end{align} \] 벡터 행렬식으로 표현하면 다음과 같다. \[ \begin{align} & \begin{bmatrix} \dot{x} \\ \dot{y} \\ \dot{z} \\ \ddot{x} \\ \ddot{y} \\ \ddot{z} \end{bmatrix} =\begin{bmatrix} 0 & 0 & 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 0 & 0 & 1 \\ 3n^2 & 0 .. 2023. 3. 6. 관성 주축 (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. 상대 궤도요소 (Relative Orbital Elements) - 2 chief 위성의 궤도가 원궤도 또는 근 원궤도(near-circular orbit)일 경우, 시간 \(t=t_0\) 에서 주어진 상대 궤도요소(ROM, relative orbital elements) \(\delta \kappa\) 를 이용하면 Hill 좌표계에서 상대 위치벡터 \(\delta \vec{r}=x\hat{o}_1+y\hat{o}_2+z\hat{o}_3\) 를 다음 식으로 표현할 수 있다 (https://pasus.tistory.com/240). \[ \begin{align} x & \approx a \delta a-a \lVert \delta \vec{e} \rVert_2 \cos (u-\varphi) \tag{1} \\ \\ y & \approx -\frac{3}{2} ua \delta.. 2023. 2. 7. 상대 궤도요소 (Relative Orbital Elements) - 1 식 (1)로 주어지는 Clohessy-Wiltshire(CW) 방정식 \[ \begin{align} & \ddot{x}-3n^2 x-2n \dot{y}=0 \tag{1} \\ \\ & \ddot{y}+2n\dot{x}=0 \\ \\ & \ddot{z}+n^2 z=0 \end{align} \] 의 해는 다음과 같았다 (https://pasus.tistory.com/239). \[ \begin{align} & x(t)= \frac{\dot{x}_0}{n} \sin nt- \left( 3x_0+\frac{2}{n} \dot{y}_0 \right) \cos nt+ \frac{2}{n} (2nx_0+ \dot{y}_0 ) \tag{2} \\ \\ & y(t)=2 \left( 3x_0+ \frac{2}{n} \do.. 2023. 2. 4. 이전 1 2 3 4 5 6 7 ··· 11 다음
