CN103942404B - 导航geo卫星太阳辐射压力矩计算方法 - Google Patents
导航geo卫星太阳辐射压力矩计算方法 Download PDFInfo
- Publication number
- CN103942404B CN103942404B CN201310026185.9A CN201310026185A CN103942404B CN 103942404 B CN103942404 B CN 103942404B CN 201310026185 A CN201310026185 A CN 201310026185A CN 103942404 B CN103942404 B CN 103942404B
- Authority
- CN
- China
- Prior art keywords
- antenna
- radiation pressure
- satellite
- pressure moment
- panel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 80
- 238000000205 computational method Methods 0.000 title abstract 2
- 238000004364 calculation method Methods 0.000 claims abstract description 13
- 238000000034 method Methods 0.000 claims description 18
- 230000011218 segmentation Effects 0.000 claims description 5
- 230000007774 longterm Effects 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 abstract description 5
- 230000000903 blocking effect Effects 0.000 abstract 1
- 238000000638 solvent extraction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Details Of Aerials (AREA)
Abstract
本发明公开了一种导航GEO卫星太阳辐射压力矩计算方法,包括以下步骤,(1)获得输入条;(2)计算太阳对星本体辐射压力矩;(3)计算太阳对太阳翼辐射压力矩;(4)计算太阳对S/L天线辐射压力矩;(5)对所述星本体辐射压力矩、所述太阳翼辐射压力矩和所述天线辐射压力矩进行傅立叶拟合得到太阳辐射压力矩的傅立叶拟合系数。利用该太阳辐射压力矩的傅立叶拟合系数进行卫星总体设计,本发明针对导航GEO卫星的形状、表面材料的光学特性,考虑了部件之间的遮挡等因素并考虑了S/L天线向北偏置的角度,同时针对网孔状S/L天线建立模型,对导航GEO卫星太阳辐射压力矩进行了计算,保证了导航GEO卫星太阳辐射压力矩计算结果的精度。
Description
技术领域
本发明属于航天飞行器设计领域,涉及导航GEO卫星太阳辐射压力矩计算方法。
背景技术
由于卫星导航具有覆盖区域大、作用距离远以及受自然灾害破坏小等优点,因此卫星导航得到了广泛应用。由于GEO卫星能覆盖固定区域,因此我国在区域星座中得到了应用。太阳辐射压力是导航GEO卫星定点以后长期工作期间的主要干扰力,其作用力使卫星轨道偏心率呈长周期变化。当太阳辐射压力相对应的压心和卫星质心不一致时,会产生太阳辐射压力力矩,这是影响卫星姿态运动的最主要干扰力矩。太阳辐射压力力矩与卫星的形状、表面材料的光学特性,部件之间的遮挡以及卫星位置等有关。
导航GEO卫星因其使命特点,采用动量轮控制卫星姿态,动量轮会不断吸收干扰力矩,到一定程度达到转速极限就会饱和,从而丧失姿态控制能力。因此需要定期对动量轮进行卸载。为了使得动量轮卸载周期满足用户要求,卫星设计之初,需要准确知道太阳干扰力矩的大小。
从以上分析可知,有必要提供导航GEO卫星太阳辐射压力矩计算方法来支持导航卫星总体设计,以支导航星座组网建设工作。
一代导航GEO卫星太阳辐射压计算时考虑了卫星的形状、表面材料的光学特性,部件之间的遮挡等因素,但没有考虑S/L天线向北偏置的角度。本算法不但考虑了卫星的形状、表面材料的光学特性,部件之间的遮挡等因素,同时考虑了S/L天线向北偏置的角度(即天线反射面按右手系绕星体X轴转过的角度),也针对网孔状S/L天线建立模型,计算结果更接近实际。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提出导航GEO卫星太阳辐射压力矩计算方法,该方法针对导航GEO卫星的形状、表面材料的光学特性,考虑了部件之间的遮挡等因素并考虑了S/L天线向北偏置的角度,同时针对网孔状S/L天线建立模型,对导航GEO卫星太阳辐射压力矩进行了计算,保证了导航GEO卫星太阳辐射压力矩计算结果的精度。
为解决上述技术问题,本发明的技术方案,通过以下过程实现:
(1)获得输入条件,所述输入条件包括:卫星构型;卫星几何尺寸;卫星表面特性参数;卫星质心位置;太阳翼位置、尺寸和表面特性参数;S/L天线位置几何尺寸和表面特性参数;
(2)根据卫星构型、几何尺寸、卫星质心以及卫星表面特性参数计算太阳对星本体辐射压力矩;星本体辐射压力矩包括星本体各个面板的辐射压力矩;所述各面板包括西面板、东面板、对地板、背地板、南面板和北面板;
(3)根据卫星质心、太阳翼位置、太阳翼尺寸和表面特性参数计算太阳对太阳翼辐射压力矩;
(4)根据卫星质心、S/L天线位置、S/L天线几何尺寸和表面特性参数计算太阳对S/L天线辐射压力矩,具体过程如下:
对所述S/L天线进行微元划分;
判断所述S/L天线受东面板的遮挡:
ya=yij+|AC|sinδ;
za=zij+|AC|cosδcosα,
若满足:-l2<ya<l1,-l5<za<l6则存在遮挡,否则无遮挡,
式中:
|AC|是天线上某微元到东板上某微元的距离;
xij是天线上某微元X坐标;
yij是天线上某微元Y坐标;
zij是天线上某微元Z坐标;
ya是东板上某微元的Y坐标;
za是东板上某微元的Z坐标;
力矩计算:
对于天线正面:
存在遮挡时,Mx8=0,My8=0,Mz8=0;
无遮挡时,光压力矩:
式中:
Mx8、My8、Mz8分别是天线正面X、Y、Z轴太阳辐射压力矩;
Mx(i,j)、My(i,j)、Mz(i,j)分别是微元X、Y、Z轴太阳辐射压力矩,按下式计算:
Mx(i,j)=[(P+τcosβ)(yijcosβNcosβ′+zijsinβNcosβ′)-τ(yijcosδcosα-zijsinδ)]KdΔs;
My(i,j)=[-(zijsinβ′+xijcosβ′cosβN)(P+τcosβ)-τ(zijcosδsinα-xijcosδcosα)]KdΔs;
Mz(i,j)=[(P+τcosβ)(yijsinβ′-xijsinβNcosβ′)-τ(xijsinδ-yijcosδsinα)]KdΔs,
式中:
Δs-微元的面积;
A1、A2、A3分别是微元的X、Y、Z坐标;
β′是天线平面与星体系yox平面的夹角;
βN是天线向北偏置的角度;
Kd是网状天线面积比例因子;
对于天线背面:
天线背面三轴力矩分别为:
Mx9=[-(P+τcosβ)(a2cosβ′cosβN+a3cosβ′sinβN)-τ(a2cosδcosα-a3sinδ)]KdπR2;
My9=πR2Kd[(P+τcosβ)(a3sinβ′+a1cosβ′cosβN)+τcosδ(a1cosα-a3sinα)];
Mz9=[(P+τcosβ)(a1sinβNcosβ′-a2sinβ′)-τ(a1sinδ-a2cosδsinα)]KdπR2;
(5)对所述星本体辐射压力矩、所述太阳翼辐射压力矩和所述天线辐射压力矩进行傅立叶拟合得到太阳辐射压力矩的傅立叶拟合系数。
利用上述太阳辐射压力矩的傅立叶拟合系数进行卫星总体设计。
进一步的,在上述步骤(2)中,采用微元分割法计算所述东面板的辐射压力矩,具体步骤为:
以所述东面板一角为坐标原点,对所述东面板进行微元划分;
根据所述微元确定所述东面板受所述S/L天线的遮挡情况;
分别计算所述东面板的有遮挡力矩和无遮挡的光压力矩。
进一步的,在上述步骤(5)中,所述傅立叶拟合以太阳在形体坐标系中的位置为自变量,将所述太阳辐射压力矩表示成傅立叶级数形式,其中包含长期项和以天为周期的项。
本发明与现有技术相比的有益效果是:
本发明分别星本体辐射压力矩、太阳翼辐射压力矩和S/L天线辐射压力矩分别进行计算,并最终进行傅里叶拟合,从而使得获得的太阳辐射压力矩的傅立叶拟合系数可更精确的指导卫星的总体设计。
本发明针对S/L天线的网状特点及于东面板相互遮挡的情况,采用微元划分方法分别计算S/L天线的正反两面的力矩,且在整个步骤中考虑了S/L天线向北偏置的角度,从而可保证获得更加精确的S/L天线的辐射压力矩。
进一步的,针对受S/L天线遮挡的东面板,采用微元划分方法对东面板进行划分,从而可对受遮挡部分和不受遮挡部分分别计算力矩,因此,同样保证了结果的精度。本发明考虑了S/L网状天线的独特构型,并建立了网孔状S/L天线的辐射压力矩模型。
附图说明
图1为本发明流程图。
具体实施方式
参考图1所示流程图,本发明方法的具体步骤如下:
第一步获得输入条件。计算所需输入条件主要包括:卫星构型、几何尺寸和卫星表面特性参数;卫星质心位置;太阳翼位置、尺寸和表面特性参数;S/L天线位置几何尺寸和表面特性参数,为计算太阳辐射压力矩做准备;
第二步根据卫星构型、几何尺寸、卫星质心以及卫星表面特性参数计算太阳对星本体各个面板的辐射压力矩,星体东面板与S/L天线存在遮挡,为了考虑遮挡的影响,采用微元分割法计算西面板辐射压力矩。
针对各面板的具体计算方法如下:
(1).西面板
X轴力矩:Mx1=-τ(J2cosδcosα-J1sinδ)
Y轴力矩:My1=-PJ1-τl4J3cosδcosα
Z轴力矩:Mz1=J2P+τl4J3sinδ
式中:
Mx1、My1、Mz1分别是西面板X、Y、Z轴太阳辐射压力矩,X、Y、Z轴分别为卫星本体坐标系的三个方向;
δ为太阳光矢量和星体系XOZ面的夹角;
是法向调节常数;
τ=kcosβ[1-ρ(1-B)];
β-面元法向和太阳光方向间的夹角;
ρ-反射系数;
τ-切向系数;
A,B-法向、切向调节系数;
κ-温度调节系数;
l4-卫星质心到西面板距离;
计算常数
计算常数;
J3=(l1+l2)(l5+l6),计算常数;
l1-卫星质心到南面板距离;
l2-卫星质心到北面板距离;
l5-卫星质心到背地板距离;
l6-卫星质心到对地板距离;
α-太阳在星体在XOZ平面内投影与星体Z轴夹角。
(2).东面板
由于东面板和天线相互遮挡,因此采用微元分割法计算,具体方法如下:
A、微元划分
选东面板一角为坐标原点,该点的坐标为(l3,-l2,-l5)
将Y方向分为m等份,Z方向分为n等份,每份的长度为:
取i=1,2,3,...,n,j=1,2,3,...,m;
第(i,j)微元中心坐标为:(l3,-l2+(j-1/2)ΔY,-l5+(i-1/2)ΔZ)
式中:
l3是卫星质心到东面板距离;
ΔY是微元Y轴方向长度;
ΔZ是微元Z轴方向长度;
n是东面板Z方向等份数;
m是东面板Y方向等份数。
B、微元被天线遮挡判断
x=lcosδsinα+A1
y=lsinδ+A2
z=lcosδcosα+A3
|A′O′|=[(x-a1)2+(y-a2)2+(z-a3)2]1/2,
当|A′O′|≥R时不遮挡,|A′O′|<R时遮挡,
式中:
l是比例因子;
(x,y,z)是天线反射面上任一点坐标;
(a1,a2,a3)为天线中心点的坐标;
R是天线口径;
|A′O′|是天线反射面上任一点到天线中心点距离;
(Nx,Ny,Nz)为天线法向坐标。
C、力矩计算
存在遮挡时,Mx2=0,My2=0,Mz2=0;
无遮挡时,东面板光压力矩:
式中:
Mx2、My2、Mz2分别是东面板X、Y、Z轴太阳辐射压力矩;
Mx(i,j)、My(i,j)、Mz(i,j)分别是微元X、Y、Z轴太阳辐射压力矩,按下式计算:
Mx(i,j)=-τ(A2cosδcosα-A3sinδ)ΔS
My(i,j)=(PA3+τA1cosδcosα)ΔS
Mz(i,j)=(-PA2-τA1sinδ)ΔS
式中:
Δs-微元的面积;
A1、A2、A3分别是微元的X、Y、Z坐标。
(3).对地板
X轴力矩:Mx3=PW1+τl6W3sinδ;
Y轴力矩:My3=-PW2-τl6W3cosδsinα;
Z轴力矩:Mz3=τ(W1cosδsinα-W2sinδ),
式中:
Mx3、My3、Mz3分别是对地板X、Y、Z轴太阳辐射压力矩;
计算常数;
计算常数;
W3=(l3+l4)(l1+l2),计算常数。
(4).背地板
X轴力矩:Mx4=-PW1-τl5W3sinδ;
Y轴力矩:My4=PW2+τl5W3cosδsinα;
Z轴力矩:Mz4=τ(W1cosδsinα-W2sinδ),
式中:
Mx4、My4、Mz4分别是背地板X、Y、Z轴太阳辐射压力矩。
(5).南面板
X轴力矩:Mx5=-PQ1-τl1Q3cosδcosα;
Y轴力矩:My5=-τ(Q1cosδsinα-Q2cosδcosα);
Z轴力矩:Mz5=PQ2+τl1Q3cosδsinα,
式中:
Mx5、My5、Mz5分别是南面板X、Y、Z轴太阳辐射压力矩;
计算常数;
计算常数;
Q3=(l3+l4)(l5+l6),计算常数。
(6).北面板
X轴力矩:Mx6=PQ1+τl2Q3cosδcosα;
Y轴力矩:My6=-τ(Q1cosδsinα-Q2cosδcosα);
Z轴力矩:Mz6=-PQ2-τl2Q3cosδsinα,
式中:
Mx6、My6、Mz6分别是北面板X、Y、Z轴太阳辐射压力矩。
第三步根据卫星质心、太阳翼位置、太阳翼尺寸和表面特性参数计算太阳对太阳翼辐射压力矩。
X轴力矩:Mx7=PI1cosα+τI3D2sinδ;
Y轴力矩:My7=I3[P(D2sinα-D1cosα)];
Z轴力矩:Mz7=-PI1sinα-τI3D1sinδ,
式中:
Mx7、My7、Mz7分别是太阳翼X、Y、Z轴太阳辐射压力矩;
计算常数;
I3=2HL,计算常数;
D1、D2分别是太阳翼转轴中心在星体XOZ面投影点的X、Z坐标;
L是太阳帆板宽;
H是太阳帆板长。
第四步根据卫星质心、S/L天线位置、S/L天线几何尺寸和表面特性参数计算太阳对S/L天线辐射压力矩,采用微元分割法分别针对天线正面和反面进行计算。
天线正面可能被东面板遮挡,用微元法求解,计算每一个微元力矩,然后叠加。具体方法如下:
A、微元的划分
将圆分为n个同心圆,最里面的圆分为两个半圆,以后从里到外按倍数增加。
B、遮挡判断
ya=yij+|AC|sinδ;
za=zij+|AC|cosδcosα,
若满足:-l2<ya<l1,-l5<za<l6则存在遮挡,否则无遮挡,
式中:
|AC|是天线上某微元到东板上某微元的距离;
xij是天线上某微元X坐标;
yij是天线上某微元Y坐标;
zij是天线上某微元Z坐标;
ya是东板上某微元的Y坐标;
za是东板上某微元的Z坐标。
C、力矩计算
对于天线正面:
存在遮挡时,Mx8=0,My8=0,Mz8=0;
无遮挡时,光压力矩:
式中:
Mx8、My8、Mz8分别是天线正面X、Y、Z轴太阳辐射压力矩;
Mx(i,j)、My(i,j)、Mz(i,j)分别是微元X、Y、Z轴太阳辐射压力矩,按下式计算:
Mx(i,j)=[(P+τcosβ)(yijcosβNcosβ′+zijsinβNcosβ′)-τ(yijcosδcosα-zijsinδ)]KdΔs;
My(i,j)=[-(zijsinβ′+xijcosβ′cosβN)(P+τcosβ)-τ(zijcosδsinα-xijcosδcosα)]KdΔs;
Mz(i,j)=[(P+τcosβ)(yijsinβ′-xijsinβNcosβ′)-τ(xijsinδ-yijcosδsinα)]KdΔs,
式中:
Δs-微元的面积;
A1、A2、A3分别是微元的X、Y、Z坐标;
β′是天线平面与星体系yox平面的夹角;
βN是天线向北偏置的角度;
Kd是网状天线面积比例因子。
对于天线背面:
天线背面三轴力矩分别为:
Mx9=[-(P+τcosβ)(a2cosβ′cosβN+a3cosβ′sinβN)-τ(a2cosδcosα-a3sinδ)]KdπR2;
My9=πR2Kd[(P+τcosβ)(a3sinβ′+a1cosβ′cosβN)+τcosδ(a1cosα-a3sinα)];
Mz9=[(P+τcosβ)(a1sinβN cosβ′-a2sinβ′)-τ(a1sinδ-a2cosδsinα)]KdπR2。
第五步根据前面计算结果得到最终太阳辐射压力矩,经过傅立叶拟合得到太阳辐射压力矩傅立叶拟合系数。
从前面计算可得到卫星不同轨道位置的三轴干扰力矩为:
Mx=Mx1+Mx2+Mx3+Mx4+Mx5+Mx6+Mx7+Mx8+Mx9
My=My1+My2+My3+My4+My5+My6+My7+My8+My9
Mz=Mz1+Mz2+Mz3+Mz4+Mz5+Mz6+Mz7+Mz8+Mz9
经过傅立叶拟合,卫星三轴太阳辐射压力矩的傅立叶级数形式如下:
式中:
N是傅立叶拟合阶数;
θ是太阳在ZOX平面内投影与Z轴的夹角,以Y轴为右手系度量,角度从0~360;
ax0是X轴力矩常值系数;
ay0是Y轴力矩常值系数;
az0是Z轴力矩常值系数;
axk、bxk、ayk、byk、azk、bzk是周期性力矩系数。
利用上述获得的结果制导卫星的总体设计,从而为获得更优的卫星总体设计精度。我国某型号卫星用了本发明计算了导航GEO卫星太阳辐射压力矩,经过卫星在轨飞行验证,导航GEO卫星太阳辐射压力矩计算结果满足实际工程需求。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
Claims (4)
1.一种导航GEO卫星太阳辐射压力矩计算方法,包括以下步骤,
(1)获得输入条件,所述输入条件包括:卫星构型;卫星几何尺寸;卫星表面特性参数;卫星质心位置;太阳翼位置、尺寸和表面特性参数;S/L天线位置几何尺寸和表面特性参数;
其特征在于,还包括以下步骤:
(2)根据卫星构型、几何尺寸、卫星质心以及卫星表面特性参数计算太阳对星本体辐射压力矩;星本体辐射压力矩包括星本体各个面板的辐射压力矩;所述各面板包括西面板、东面板、对地板、背地板、南面板和北面板;
(3)根据卫星质心、太阳翼位置、太阳翼尺寸和表面特性参数计算太阳对太阳翼辐射压力矩;
(4)根据卫星质心、S/L天线位置、S/L天线几何尺寸和表面特性参数计算太阳对S/L天线辐射压力矩,具体过程如下:
对所述S/L天线进行微元划分;
判断所述S/L天线受东面板的遮挡:
ya=yij+|AC|sinδ;
za=zij+|AC|cosδcosα,
若满足:-l2<ya<l1,-l5<za<l6则存在遮挡,否则无遮挡,
式中:
|AC|是天线上某微元到东板上某微元的距离;
xij是天线上某微元X坐标;
yij是天线上某微元Y坐标;
zij是天线上某微元Z坐标;
ya是东板上某微元的Y坐标;
za是东板上某微元的Z坐标;
l1是质心到南面板距离;
l2是质心到北面板距离;
l3是质心到东面板距离;
l5是质心到背地板距离;
l6是质心到对地板距离;
δ是太阳光矢量和星体系XOZ面的夹角;
α是太阳光矢量在星体坐标系XOZ平面内投影与Z轴的夹角;
力矩计算:
对于天线正面:
存在遮挡时,Mx8=0,My8=0,Mz8=0;
无遮挡时,光压力矩:
式中:
Mx8、My8、Mz8分别是天线正面X、Y、Z轴太阳辐射压力矩;
Mx(i,j)、My(i,j)、Mz(i,j)分别是微元X、Y、Z轴太阳辐射压力矩,按下式计算:
Mx(i,j)=[(P+τcosβ)(yij cosβN cosβ'+zij sinβN cosβ')-τ(yij cosδcosα-zij sinδ)]KdΔs;
My(i,j)=[-(zij sinβ'+xij cosβ'cosβN)(P+τcosβ)-τ(zij cosδsinα-xij cosδcosα)]KdΔs;
Mz(i,j)=[(P+τcosβ)(yij sinβ'-xij sinβN cosβ')-τ(xij sinδ-yij cosδsinα)]KdΔs,
式中:
Δs-微元的面积;
β'是天线平面与星体系yox平面的夹角;
βN是天线向北偏置的角度;
Kd是网状天线面积比例因子;
是法向调节常数;
τ=kcosβ[1-ρ(1-B)];
对于天线背面:
天线背面三轴力矩分别为:
Mx9=[-(P+τcosβ)(a2 cosβ'cosβN+a3 cosβ'sinβN)-τ(a2 cosδcosα-a3 sinδ)]KdπR2;
My9=πR2Kd[(P+τcosβ)(a3 sinβ'+a1 cosβ'cosβN)+τcosδ(a1 cosα-a3 sinα)];
Mz9=[(P+τcosβ)(a1 sinβN cosβ'-a2 sinβ')-τ(a1 sinδ-a2 cosδsinα)]KdπR2;
其中,
R是天线的半径,天线中心点的坐标为(a1,a2,a3);
a1天线中心点X坐标;
a2天线中心点Y坐标;
a3天线中心点Z坐标;
(5)对所述星本体辐射压力矩、所述太阳翼辐射压力矩和所述天线辐射压力矩进行傅立叶拟合得到太阳辐射压力矩的傅立叶拟合系数。
2.根据权利要求1所述导航GEO卫星太阳辐射压力矩计算方法,其特征在于:在所述步骤(2)中,采用微元分割法计算所述东面板的辐射压力矩,具体步骤为:
以所述东面板一角为坐标原点,对所述东面板进行微元划分;
根据所述微元确定所述东面板受所述S/L天线的遮挡情况;
分别计算所述东面板的有遮挡力矩和无遮挡的光压力矩。
3.根据权利要求2所述导航GEO卫星太阳辐射压力矩计算方法,其特征在于:所述步骤(5)中,所述傅立叶拟合以太阳在形体坐标系中的位置为自变量,将所述太阳辐射压力矩表示成傅立叶级数形式,其中包含长期项和以天为周期的项。
4.根据权利要求1~3中任一所述导航GEO卫星太阳辐射压力矩计算方法,其特征在于:所述太阳辐射压力矩的傅立叶拟合系数用于卫星总体设计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310026185.9A CN103942404B (zh) | 2013-01-21 | 2013-01-21 | 导航geo卫星太阳辐射压力矩计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310026185.9A CN103942404B (zh) | 2013-01-21 | 2013-01-21 | 导航geo卫星太阳辐射压力矩计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103942404A CN103942404A (zh) | 2014-07-23 |
CN103942404B true CN103942404B (zh) | 2017-03-15 |
Family
ID=51190072
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310026185.9A Active CN103942404B (zh) | 2013-01-21 | 2013-01-21 | 导航geo卫星太阳辐射压力矩计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103942404B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108959734B (zh) * | 2018-06-15 | 2023-08-29 | 北京控制工程研究所 | 一种基于实时递推太阳光压力矩辨识方法及*** |
CN108920829B (zh) * | 2018-06-29 | 2022-05-24 | 中国空间技术研究院 | 一种带大型网状天线卫星的太阳光压力矩计算方法 |
CN110132582B (zh) * | 2019-05-21 | 2021-11-09 | 北京空间飞行器总体设计部 | 真空环境下无摩擦的力矩加载及测量装置 |
CN112163305B (zh) * | 2020-09-07 | 2023-07-25 | 北京空间飞行器总体设计部 | 一种网状天线卫星的光压计算方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0949144A2 (en) * | 1998-04-09 | 1999-10-13 | Hughes Electronics Corporation | Methods for reducing solar array power variations while managing the system influences of operating with off-pointed solar wings on spacecraft |
CN101226561A (zh) * | 2007-12-28 | 2008-07-23 | 南京航空航天大学 | 用于微型航天器姿态轨道控制***的微型仿真支持***及工作方法 |
-
2013
- 2013-01-21 CN CN201310026185.9A patent/CN103942404B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0949144A2 (en) * | 1998-04-09 | 1999-10-13 | Hughes Electronics Corporation | Methods for reducing solar array power variations while managing the system influences of operating with off-pointed solar wings on spacecraft |
CN101226561A (zh) * | 2007-12-28 | 2008-07-23 | 南京航空航天大学 | 用于微型航天器姿态轨道控制***的微型仿真支持***及工作方法 |
Non-Patent Citations (5)
Title |
---|
Feedback Linearization and Solar Pressure Satellite Attitude Control;S.N.Singh 等;《IEEE TRANSACTIONS ON AEROAPACE AND ELECTRONIC SYSTEMS》;19960430;第32卷(第2期);第732-741页 * |
SOLAR PRESSURE EVALUATION ON LARGE REFLECTORS FOR SPACE APPLICATIONS;P.BRUNELLO;《COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING》;19931231;第9卷(第10期);第787-795页 * |
卫星太阳光压力矩计算中有效作用面积的计算;刘暾;《南京理工大学学报》;20071231;第31卷(第6期);第684-688页 * |
双翼箱式航天器太阳光压干扰力矩模型;李连军;《航天控制》;20050831;第23卷(第4期);第13-17页 * |
导航卫星太阳光压任务分析;高照照;《第二届中国卫星导航学术年会》;20111231;第1-3页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103942404A (zh) | 2014-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103942404B (zh) | 导航geo卫星太阳辐射压力矩计算方法 | |
WO2017113567A1 (zh) | 火星探测器自主导航方法 | |
CN104462776B (zh) | 一种低轨道地球观测卫星对月球绝对辐射定标方法 | |
CN107655485B (zh) | 一种巡航段自主导航位置偏差修正方法 | |
CN103017760B (zh) | 一种大椭圆轨道火星探测器自主对火定向方法 | |
CN108181916B (zh) | 小卫星相对姿态的控制方法及装置 | |
CN102322862B (zh) | 一种编队飞行卫星绝对和相对轨道确定方法 | |
Takahashi et al. | Spin state and moment of inertia characterization of 4179 Toutatis | |
CN107797130A (zh) | 低轨航天器多点多参数轨道上行数据计算方法 | |
CN104960674A (zh) | 一种运动目标的指向跟踪控制方法 | |
CN104457705A (zh) | 基于天基自主光学观测的深空目标天体初定轨方法 | |
CN106197425A (zh) | 基于卫星姿态角的地面目标点位置的计算方法 | |
CN108959734B (zh) | 一种基于实时递推太阳光压力矩辨识方法及*** | |
CN102519455B (zh) | 基于紫外敏感器的自主导航半物理仿真试验*** | |
Früh et al. | Coupled orbit–attitude motion of high area-to-mass ratio (hamr) objects including efficient self-shadowing | |
CN104021241A (zh) | 一种导航卫星复杂模型的地球反照光压摄动建模方法 | |
CN112783183B (zh) | 一种太阳同步圆回归轨道的轨道规划方法 | |
Godard et al. | Orbit determination of Rosetta around comet 67P/Churyumov-Gerasimenko | |
CN112629543A (zh) | 一种大椭圆轨道及小倾角圆轨道的轨道规划方法 | |
Linares et al. | Space object mass-specific inertia matrix estimation from photometric data | |
Feng et al. | Autonomous orbit determination and its error analysis for deep space using X-ray pulsar | |
CN108680171A (zh) | 一种对日定向卫星星敏感器安装指向获取方法和*** | |
CN113589832B (zh) | 对地表固定区域目标稳定观测覆盖的星座快速设计方法 | |
Huang et al. | Precise orbit determination of a maneuvered GEO satellite using CAPS ranging data | |
CN109269508A (zh) | 一种卫星相对小行星视觉自主导航方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |