航天器轨道力学理论与方法笔记

所有知识点均参考《航天器轨道力学理论与方法》. 张洪波著

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(r1rc)+m2(r2rc)=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)=(AC)B(AB)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×hr/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(1e2),且由 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+rf

椭圆轨道面积微元 d σ = 1 2 r ⋅ r d f d\sigma=\frac{1}{2}r\cdot rdf dσ=21rrdf两边同除以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=a1e2

平均角速度: 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τ)=EesinE,E表示偏近点角, M = n ( t − τ ) M=n(t-\tau) M=n(tτ)为平近点角。

能量守恒: v 2 2 − μ r = ε \frac{v^2}{2}-\frac{\mu}{r}=\varepsilon 2v2rμ=ε

r ⋅ r ˙ = r ⋅ r ˙ \boldsymbol r \cdot \dot \boldsymbol r=r\cdot \dot r rr˙=rr˙,这是因为 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=μ(r2a1)

惠特克定理:飞行速度分解到垂直于长轴分量和垂直于矢径分量时为常数

  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值