所有知识点均参考《航天器轨道力学理论与方法》. 张洪波著
2021.7.29
辐照度:投射到单位面积上的辐射能量。
太阳常数:地球大气上界处太阳辐照度
岁差:地轴长期运动。
章动:地轴短周期运动
质心定义: r c = m 1 r 1 + m 2 r 2 m 1 + m 2 \boldsymbol r_c = \frac{m_1\boldsymbol r_1 + m_2 \boldsymbol r_2}{m_1+m_2} rc=m1+m2m1r1+m2r2
质心性质: m 1 ( r 1 − r c ) + m 2 ( r 2 − r c ) = 0 m_1(\boldsymbol r_1-\boldsymbol r_c)+m_2(\boldsymbol r_2-\boldsymbol r_c) = 0 m1(r1−rc)+m2(r2−rc)=0
第三章 3.1
二体系统的质心无加速度(从系统的角度来看,二体系统不受外力作用,质心加速度必然为0),即 r ¨ c = 0 \ddot \boldsymbol r_c=0 r¨c=0。
二体问题基本运动方程: r ¨ + μ r 3 r = 0 \ddot \boldsymbol r+\frac{\mu}{r^3}\boldsymbol r=0 r¨+r3μr=0
第三章 3.2.1
动量矩为 h = r × r ˙ \boldsymbol h=\boldsymbol r \times \dot \boldsymbol r h=r×r˙,在二体系统中, d d t ( r × r ˙ ) = 0 \frac{d}{dt}(\boldsymbol r \times \dot \boldsymbol r)=0 dtd(r×r˙)=0,所以二体问题中动量矩守恒。动量矩守恒说明航天器位置和速度矢量始终处于惯性空间中某固定平面中,所以二体运动是平面运动。
h = ∣ h ∣ = r 2 u ˙ h=|\boldsymbol h|=r^2\dot u h=∣h∣=r2u˙, u u u为掠面速率, h = h [ s i n i s i n Ω , − s i n i c o s Ω , c o s i ] T \boldsymbol h=h[sini sin\Omega,-sini cos\Omega, cosi]^T h=h[sinisinΩ,−sinicosΩ,cosi]T,这个可以看作 h \boldsymbol h h从椭圆轨道坐标系转换到地心惯性坐标系。 h \boldsymbol h h在椭圆轨道坐标系表示为: [ 0 , 0 , h ] T [0,0,h]^T [0,0,h]T
混合积: ( A × B ) ⋅ C = ( B × C ) ⋅ A = ( C × A ) ⋅ B (A\times B) \cdot C=(B \times C)\cdot A=(C\times A)\cdot B (A×B)⋅C=(B×C)⋅A=(C×A)⋅B
A × ( B × C ) = ( A ⋅ C ) B − ( A ⋅ B ) C A\times(B\times C)=(A\cdot C)B-(A\cdot B)C A×(B×C)=(A⋅C)B−(A⋅B)C
轨道积分方程: r = h 2 / μ 1 + e c o s f r=\frac{h^2/\mu}{1+ecosf} r=1+ecosfh2/μ, f f f为真近点角, e = ∣ e ∣ e=|\boldsymbol e| e=∣e∣,而 e \boldsymbol e e称为偏心率矢量。 e = v × h μ − r / r \boldsymbol e = \frac{\boldsymbol v\times \boldsymbol h}{\mu}-\boldsymbol r/r e=μv×h−r/r, e \boldsymbol e e在轨道面内方向固定,f=0时指向近心点方向。
半通径 p = h 2 μ = a ( 1 − e 2 ) p=\frac{h^2}{\mu}=a(1-e^2) p=μh2=a(1−e2),且由 r = p 1 + e c o s f r=\frac{p}{1+ecosf} r=1+ecosfp得 p = r + r ⋅ f p=r+\boldsymbol r\cdot \boldsymbol f p=r+r⋅f
椭圆轨道面积微元 d σ = 1 2 r ⋅ r d f d\sigma=\frac{1}{2}r\cdot rdf dσ=21r⋅rdf两边同除以dt可得 σ ˙ = 1 2 r 2 f ˙ \dot \sigma=\frac{1}{2}r^2\dot f σ˙=21r2f˙,则椭圆轨道周期为 T = S σ ˙ = 2 π a 3 μ , S = π a b , b = a 1 − e 2 T=\frac{S}{\dot \sigma}=2\pi \sqrt{\frac{a^3}{\mu}},S=\pi ab,b=a\sqrt {1-e^2} T=σ˙S=2πμa3,S=πab,b=a1−e2
平均角速度: n = 2 π / T = μ a 3 n=2\pi /T=\sqrt{\frac{\mu}{a^3}} n=2π/T=a3μ
时间积分: n ( t − τ ) = E − e s i n E , E n(t-\tau)=E-esinE,E n(t−τ)=E−esinE,E表示偏近点角, M = n ( t − τ ) M=n(t-\tau) M=n(t−τ)为平近点角。
能量守恒: v 2 2 − μ r = ε \frac{v^2}{2}-\frac{\mu}{r}=\varepsilon 2v2−rμ=ε
r ⋅ r ˙ = r ⋅ r ˙ \boldsymbol r \cdot \dot \boldsymbol r=r\cdot \dot r r⋅r˙=r⋅r˙,这是因为 r ˙ \dot \boldsymbol r r˙包括切向分量和横向分量(见《航天器轨道力学理论与方法》 P 57 P_{57} P57)其中切向分量就是 r \boldsymbol r r所在方向,且切向分量大小为 r ˙ \dot r r˙。
ε = − μ 2 a \varepsilon=-\frac{\mu}{2a} ε=−2aμ
活力公式(能量积分): v 2 = μ ( 2 r − 1 a ) v^2=\mu(\frac{2}{r}-\frac{1}{a}) v2=μ(r2−a1)
惠特克定理:飞行速度分解到垂直于长轴分量和垂直于矢径分量时为常数