CN105183948A - 一种基于二次反射的高精度卫星太阳光压摄动力建模方法 - Google Patents
一种基于二次反射的高精度卫星太阳光压摄动力建模方法 Download PDFInfo
- Publication number
- CN105183948A CN105183948A CN201510497095.7A CN201510497095A CN105183948A CN 105183948 A CN105183948 A CN 105183948A CN 201510497095 A CN201510497095 A CN 201510497095A CN 105183948 A CN105183948 A CN 105183948A
- Authority
- CN
- China
- Prior art keywords
- solar radiation
- surface elements
- satellite
- sun
- theta
- 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.)
- Granted
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Photovoltaic Devices (AREA)
Abstract
本发明公开了一种基于二次反射的高精度卫星太阳光压摄动力建模方法,采用模拟太阳光源照射卫星表面部件,研究光线入射后反射再次入射到表面部件后,表面部件的受到的光压摄动力,最终得到太阳光压摄动加速度模型,该建模方法能够建立复杂卫星包含二次反射的高精度太阳光压摄动模型,作为中高轨卫星最主要的非保守力,能够完善并提高高精度卫星动力学模型,进一步提高精密定轨、轨道预报精度。
Description
技术领域
本发明属于飞行器设计领域,涉及一种基于二次反射的高精度卫星太阳光压摄动力建模方法。
背景技术
随着高精度卫星精密定轨与轨道预报精度要求提高,尽可能精确的卫星动力学模型成为工程需要。以导航卫星为例,光压摄动作为目前中高轨卫星最主要的误差源,其精确度直接影响卫星的精密定轨与轨道预报。太阳光压是光压摄动最主要的摄动源,因此,考虑二次反射的高精度卫星太阳光压摄动力建模的研究是必要的。
发明内容
有鉴于此,本发明提供了一种基于二次反射的高精度卫星太阳光压摄动力建模方法,能够完善并提高高精度卫星动力学模型,进一步提高精密定轨、轨道预报精度。
本发明的一种基于二次反射的卫星太阳光压摄动力建模方法,包括如下步骤:
步骤1、获取待建模卫星的表面三维模型,确定卫星的所有表面部件,将每个表面部件的表面划分成多个面元;
步骤2、用面积等于卫星外形最大包络的方形像元阵列作为模拟太阳光源,每个像元作为一个小光源并发出一条垂直于方形像元阵列的太阳光线,所述模拟太阳光源与卫星的距离为一个天文单位;
步骤3、针对所述卫星的各个表面部件,记录所述模拟太阳光源的每个像元发出的光线与表面部件最短距离的交点和相应的入射角,统计第i个表面部件上有所述交点的面元数量,即为初次照射第i个表面部件的有效被照射面元数量Aieff1;i=1,2,...N,其中N为表面部件数量;
根据有效被照射面元数量Aieff1,得到第i个表面部件由模拟太阳光源初次照射所受的法向光压摄动力fni1以及切向光压摄动力fsi1;
步骤4、首先,针对各个表面部件,根据步骤3获得的所有入射光线的交点和入射角,继而确定入射光线在所在面元上的反射光线;
其次,统计第i个表面部件与所述反射光线有最短距离交点的面元个数Aieff2,则第i个表面部件被二次照射中漫反射的有效面元数为Aieff2·v(1-μ),镜面反射的有效面元数为Aieff2·vμ;
然后,根据Aieff2·v(1-μ),获得第i个表面部件在二次照射中由漫反射引起的法向光压摄动力fni2d;根据Aieff2·vμ,获得第i个表面部件在二次照射中由镜面反射引起的法向光压摄动力fni2m和切向光压摄动力fsi2m;
最后,得到第i个部件所受的总的法向光压摄动力fni为:
fni=fni1+fni2d+fni2m;
总的切向光压摄动力fsi为:
fsi=fsi1+fsi2m;
步骤5、根据各个表面部件的安装位置和法线方向,将所有表面部件的法向光压摄动力fni和切向光压摄动力fsi分别分解到卫星本体坐标系的三个轴方向,并得到各轴方向的合力,最后得到各轴方向的光压摄动加速度;
步骤6、采用步骤2-5的方法,在一个卫星运动周期内,以设定角度为步长,计算设定时间内每个卫星运动周期不同时刻的太阳光压摄动加速度;针对计算结果,通过傅里叶多项式拟合得到太阳光压摄动加速度的数学模型。
较佳的,所述步骤6中的设定时间为半年。
较佳的,所述像元大小为1mm×1mm。
较佳的,所述面元大小为1mm×1mm。
较佳的,所述初次照射法向光压摄动力fni1为:
切向光压摄动力fsi1为:
其中,j=1,2,...,M,M为表面部件上面元数量,τ为日蚀因子,c为光速,Esun=1367W/m2为太阳辐射强度;θij为第i个表面部件的第j个面元的光线入射角;νi表示第i个表面部件的反射率、μi第i个表面部件的镜面系数。
较佳的,二次照射中,第i个表面部件由漫反射引起的法向光压摄动力fni2d为:
第i个表面部件由镜面反射引起的法向光压摄动力fni2m和切向光压摄动力fsi2m为:
其中,j=1,2,...,M,M为表面部件上面元数量,τ为日蚀因子,c为光速,Esun=1367W/m2为太阳辐射强度;θij为第i个表面部件的第j个面元的光线入射角;νi表示第i个表面部件的反射率、μi第i个表面部件的镜面系数。
较佳的,所述步骤6中所述设定角度为10°。
本发明具有如下有益效果:
(1)本发明的建模方法能够建立复杂卫星包含二次反射的高精度太阳光压摄动模型,作为中高轨卫星最主要的非保守力,能够完善并提高高精度卫星动力学模型,进一步提高精密定轨、轨道预报精度。
附图说明
图1为卫星本体坐标系示意图;
图2为卫星轨道弧段示意图;
图3为太阳光线像元阵列示意图;
图4为本发明的建模方法流程图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
如图4所示,本发明的建模方法的具体步骤如下:
1)选取卫星本体坐标系作为太阳光压摄动建模的基准坐标系,如图1所示,以卫星的中心的坐标原点,沿太阳翼方向为Y轴,X轴垂直Y同时垂直卫星的一个表面,Z轴同时垂直于Y轴、X轴以及卫星的一个表面;
2)基于卫星本体坐标系,获取待建模研究卫星的表面三维模型,对卫星所有表面部件编号记为i,i为1,…,N,建立包括每个卫星表面部件的卫星表面参数数据库,每个表面部件的信息包括:部件的形状(例如矩形、圆柱、抛物面等)、部件相对本体系的安装位置(rxi,ryi,rzi)T、表面积Ai、反射率vi、镜面系数μi、表面部件的外法线在本体系中的单位矢量Pi=(xi,yi,zi)T;
3)卫星姿态控制模式和太阳矢量运动规律分析。设太阳在J2000惯性系中的矢量为卫星的轨道六根数中升交点经度Ω、轨道倾角i、升交点幅角u,J2000惯性系下的太阳矢量依次经过3-1-2三次坐标转换,可以转到轨道坐标系,旋转矩阵为:
再设卫星按321转序姿态角为:滚动角俯仰角θ、偏航角ψ,则轨道坐标系与卫星本体坐标系之间的转换矩阵为:
则太阳矢量在基准本体系中的矢量为:
本体系中,沿太阳矢量方向的单位向量记为
4)根据卫星在轨运行参数,判断卫星的光照条件,卫星在轨所经历的不同光照弧段如图2所示。
设τ为日蚀因子,日照区τ=1,半影区τ线性减小,本影区减为τ=0。每个表面部件受的摄动力根据所处弧段乘对应的日蚀因子。
5)用面积等于卫星外形最大包络的方形像元阵列模拟太阳光源,像元阵列的指向垂直于太阳矢量,距离卫星本体系为一个天文单位,如图3所示,相邻像元阵列的间隔1mm。定义像元阵列坐标系为:坐标原点在像元阵列的中心点,Z轴沿太阳矢量方向,X轴在像元阵列内,平行于像元阵列某一边长,Z轴与X、Y轴符合右手旋转。设每个像元中心点在像元阵列坐标系中的坐标为Pj,k=(j,k,0),在每个像元中心添加一个太阳到卫星方向的单位矢量,作为模拟光线,记为Rj,k(j表示矢量在阵列中所处的行,k表示矢量在阵列中所处的列)。j、k的范围为-N,…,0,…,N,N的大小由卫星表面积的最大包络确定。
根据像元阵列坐标系的定义和本体系中太阳矢量的指向,本体坐标系与像元阵列坐标系的转换矩阵为:
L=Ry(θy)·Rz(θz)
其中:
每个像元中心点在本体系中的坐标为:
6)计算初次照射引起的太阳光压摄动力。由本体系中太阳指向、像元中心点坐标、以及部件数据库,建立光路跟踪函数,判断每个表面部件的实际受照射情况及有效受照面积。通过每个表面部件的安装位置,表面积,确定部件在本体系中的空间函数及边界。记某一像元中心点在本体系中的坐标为(xsi,ysi,zsi),太阳矢量为(Xs,Ys,Zs)T,则太阳光线的空间方程为:
计算并记录每条光线与卫星表面距离最短的交点(即为光线所照射到的点)和相应的入射角(光线和被照射面元法线的夹角),并对该部件有效受照光线数量加1,对每一条照射光线对应有效被照射面元面积记为Δ=1mm2。通过与卫星部件存在最短距离的光线数,确定初次照射时每个表面部件有效被照射面元数为Aieff1。
对每一个受太阳光初次照射的部件,以i表示卫星表面部件代号,其所受法向光压摄动力fni1为:
切向光压摄动力fsi1为:
式中:Esun=1367W/m2为太阳辐射强度;θij为第i个表面部件的第j个面元的光线入射角。
7)计算二次照射引起的太阳光压摄动力。由于卫星是个构型复杂的结构体,初次入射光线的反射光线会再次照射在某些表面部件上,形成二次太阳光压摄动。
根据卫星表面部件在空间的安装位置、外法向方向Pi=(xi,yi,zi)T,确定固定平面部件面元法向与部件相同、固定曲面部件面元法向为星体坐标系原点指向面元中心,转动部件的面元法向由部件运动情况决定。综合考虑太阳光线在星本体坐标系中的矢量,确定反射光线的方向:漫反射光线沿面元法向,镜面反射光线沿入射光线的对称方向。
由漫反射和镜面反射光线的方向和光线出发点(即首次被照射面元点),利用第6)步中光路跟踪函数法计算反射光线所照射到的表面面元,对每一条照射光线对应有效被照射面元面积记为Δ=1mm2。通过反射光线与部件最短距离的判断,确定每个表面部件由二次照射的有效被照射面元数为Aieff2。根据反射类型,被反射光线照射的有效被照射面元数为Aieff2乘以相应的比例,漫反射的比例为v(1-μ),镜面反射的比例为vμ。因此被二次照射中漫反射的有效面元数为Aieff2·v(1-μ),镜面反射的有效面元数为Aieff2·vμ。
二次照射中,第i个表面部件由漫反射引起的法向光压摄动力fni2d为:
二次照射中,第i个表面部件由镜面反射引起的法向光压摄动力fni2m和切向光压摄动力fsi2m为:
所以,第i个部件所受法向光压摄动力fni为:
fni=fni1+fni2d+fni2m
切向光压摄动力fsi为:
fsi=fsi1+fsi2m
8)根据部件的安装位置、法线方向等,将所有表面部件的太阳光压摄动力(fni,fsi)分解为沿本体系三轴方向的三个分力,分别记为(Fxi,Fyi,Fzi)。
并对所有表面部件的地球反照光压摄动力矢量叠加,计算整星在本体系下的太阳光压摄动力和光压摄动加速度:
9)步骤1)-8)是计算特定太阳矢量方向的太阳光压摄动加速度,在一个卫星运动周期内,以卫星在轨道内运动10°为步长,计算半年内每个卫星运动周期不同时刻的太阳光压摄动加速度。对计算结果,通过傅里叶多项式拟合得到太阳光压摄动加速度的数学模型。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种基于二次反射的卫星太阳光压摄动力建模方法,其特征在于,包括如下步骤:
步骤1、获取待建模卫星的表面三维模型,确定卫星的所有表面部件,将每个表面部件的表面划分成多个面元;
步骤2、用面积等于卫星外形最大包络的方形像元阵列作为模拟太阳光源,每个像元作为一个小光源并发出一条垂直于方形像元阵列的太阳光线,所述模拟太阳光源与卫星的距离为一个天文单位;
步骤3、针对所述卫星的各个表面部件,记录所述模拟太阳光源的每个像元发出的光线与表面部件最短距离的交点和相应的入射角,统计第i个表面部件上有所述交点的面元数量,即为初次照射第i个表面部件的有效被照射面元数量Aieff1;i=1,2,...N,其中N为表面部件数量;
根据有效被照射面元数量Aieff1,得到第i个表面部件由模拟太阳光源初次照射所受的法向光压摄动力fni1以及切向光压摄动力fsi1;
步骤4、首先,针对各个表面部件,根据步骤3获得的所有入射光线的交点和入射角,继而确定入射光线在所在面元上的反射光线;
其次,统计第i个表面部件与所述反射光线有最短距离交点的面元个数Aieff2,则第i个表面部件被二次照射中漫反射的有效面元数为Aieff2·v(1-μ),镜面反射的有效面元数为Aieff2·vμ;
然后,根据Aieff2·v(1-μ),获得第i个表面部件在二次照射中由漫反射引起的法向光压摄动力fni2d;根据Aieff2·vμ,获得第i个表面部件在二次照射中由镜面反射引起的法向光压摄动力fni2m和切向光压摄动力fsi2m;
最后,得到第i个部件所受的总的法向光压摄动力fni为:
fni=fni1+fni2d+fni2m;
总的切向光压摄动力fsi为:
fsi=fsi1+fsi2m;
步骤5、根据各个表面部件的安装位置和法线方向,将所有表面部件的法向光压摄动力fni和切向光压摄动力fsi分别分解到卫星本体坐标系的三个轴方向,并得到各轴方向的合力,最后得到各轴方向的光压摄动加速度;
步骤6、采用步骤2-5的方法,在一个卫星运动周期内,以设定角度为步长,计算设定时间内每个卫星运动周期不同时刻的太阳光压摄动加速度;针对计算结果,通过傅里叶多项式拟合得到太阳光压摄动加速度的数学模型。
2.如权利要求1所述的一种基于二次反射的卫星太阳光压摄动力建模方法,其特征在于,所述步骤6中的设定时间为半年。
3.如权利要求1所述的一种基于二次反射的卫星太阳光压摄动力建模方法,其特征在于,所述像元大小为1mm×1mm。
4.如权利要求1所述的一种基于二次反射的卫星太阳光压摄动力建模方法,其特征在于,所述面元大小为1mm×1mm。
5.如权利要求1所述的一种基于二次反射的卫星太阳光压摄动力建模方法,其特征在于:
所述初次照射法向光压摄动力fni1为:
切向光压摄动力fsi1为:
其中,j=1,2,...,M,M为表面部件上面元数量,τ为日蚀因子,c为光速,Esun=1367W/m2为太阳辐射强度;θij为第i个表面部件的第j个面元的光线入射角;νi表示第i个表面部件的反射率、μi第i个表面部件的镜面系数。
6.如权利要求1所述的一种基于二次反射的卫星太阳光压摄动力建模方法,其特征在于:二次照射中,第i个表面部件由漫反射引起的法向光压摄动力fni2d为:
第i个表面部件由镜面反射引起的法向光压摄动力fni2m和切向光压摄动力fsi2m为:
其中,j=1,2,...,M,M为表面部件上面元数量,τ为日蚀因子,c为光速,Esun=1367W/m2为太阳辐射强度;θij为第i个表面部件的第j个面元的光线入射角;νi表示第i个表面部件的反射率、μi第i个表面部件的镜面系数。
7.如权利要求1所述的一种基于二次反射的卫星太阳光压摄动力建模方法,其特征在于:步骤6中所述设定角度为10°。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510497095.7A CN105183948B (zh) | 2015-08-13 | 2015-08-13 | 一种基于二次反射的高精度卫星太阳光压摄动力建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510497095.7A CN105183948B (zh) | 2015-08-13 | 2015-08-13 | 一种基于二次反射的高精度卫星太阳光压摄动力建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105183948A true CN105183948A (zh) | 2015-12-23 |
CN105183948B CN105183948B (zh) | 2018-04-10 |
Family
ID=54906027
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510497095.7A Active CN105183948B (zh) | 2015-08-13 | 2015-08-13 | 一种基于二次反射的高精度卫星太阳光压摄动力建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105183948B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108460176A (zh) * | 2018-01-02 | 2018-08-28 | 佛山科学技术学院 | 一种提高卫星轨道摄动力模型表示精度的方法 |
CN108490966A (zh) * | 2018-01-31 | 2018-09-04 | 中国人民解放军国防科技大学 | 基于微分代数的静止轨道摄动相对轨迹高阶制导方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5626315A (en) * | 1990-11-30 | 1997-05-06 | Aerospatiale Societe Nationale Industrielle | Apparatus controlling the pitch attitude of a satellite by means of solar radiation pressure |
CN103488833A (zh) * | 2013-09-24 | 2014-01-01 | 北京空间飞行器总体设计部 | 一种导航卫星复杂模型的热辐射力建模方法 |
CN104021241A (zh) * | 2014-05-13 | 2014-09-03 | 北京空间飞行器总体设计部 | 一种导航卫星复杂模型的地球反照光压摄动建模方法 |
-
2015
- 2015-08-13 CN CN201510497095.7A patent/CN105183948B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5626315A (en) * | 1990-11-30 | 1997-05-06 | Aerospatiale Societe Nationale Industrielle | Apparatus controlling the pitch attitude of a satellite by means of solar radiation pressure |
CN103488833A (zh) * | 2013-09-24 | 2014-01-01 | 北京空间飞行器总体设计部 | 一种导航卫星复杂模型的热辐射力建模方法 |
CN104021241A (zh) * | 2014-05-13 | 2014-09-03 | 北京空间飞行器总体设计部 | 一种导航卫星复杂模型的地球反照光压摄动建模方法 |
Non-Patent Citations (1)
Title |
---|
陈秋丽等: "导航卫星高精度太阳光压摄动分析建模及验证", 《中国空间科学技术》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108460176A (zh) * | 2018-01-02 | 2018-08-28 | 佛山科学技术学院 | 一种提高卫星轨道摄动力模型表示精度的方法 |
CN108490966A (zh) * | 2018-01-31 | 2018-09-04 | 中国人民解放军国防科技大学 | 基于微分代数的静止轨道摄动相对轨迹高阶制导方法 |
CN108490966B (zh) * | 2018-01-31 | 2021-02-05 | 中国人民解放军国防科技大学 | 基于微分代数的静止轨道摄动相对轨迹高阶制导方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105183948B (zh) | 2018-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sánchez-González et al. | Solar flux distribution on central receivers: A projection method from analytic function | |
Yu et al. | Modeling and simulation of 1 MWe solar tower plant’s solar flux distribution on the central cavity receiver | |
Wang et al. | A novel integrated simulation approach couples MCRT and Gebhart methods to simulate solar radiation transfer in a solar power tower system with a cavity receiver | |
Lin et al. | Experimental and theoretical analysis on a linear Fresnel reflector solar collector prototype with V-shaped cavity receiver | |
CN108920829B (zh) | 一种带大型网状天线卫星的太阳光压力矩计算方法 | |
CN104021241B (zh) | 一种导航卫星复杂模型的地球反照光压摄动建模方法 | |
CN104573251A (zh) | 一种星载光学遥感器全视场表观光谱辐亮度确定方法 | |
CN114936471B (zh) | 一种基于并行计算的航天器碰撞预警分层快速筛选方法 | |
CN101858747A (zh) | 一种有效利用地球辐照能的卫星帆板对日定向目标姿态的解析确定方法 | |
CN105445738A (zh) | 基于遗传算法的geo星机双基sar接收站飞行参数设计方法 | |
CN104898264A (zh) | 一种大口径旋转对称非成像自由曲面反射镜及其设计方法 | |
CN111680455B (zh) | 基于搭载形式的撞击探测轨道设计方法和系统 | |
CN106096204A (zh) | 一种基于太阳帆推进技术的航天器日心椭圆悬浮轨道设计方法 | |
CN103593551A (zh) | 一种基于高轨光学成像卫星的太阳光压计算方法 | |
CN102564574B (zh) | 一种测量地球反照辐射照度的方法 | |
Lee et al. | Numerical study on optical performances of the first central-receiver solar thermal power plant in Korea | |
CN105486315B (zh) | 遥感卫星对月绝对定标姿态调整方法 | |
Hachicha | Numerical modelling of a parabolic trough solar collector | |
CN105183948A (zh) | 一种基于二次反射的高精度卫星太阳光压摄动力建模方法 | |
CN104408279B (zh) | 一种运载火箭空间外热流计算方法 | |
CN103234552A (zh) | 光学导航目标星模拟仿真图像生成方法 | |
CN105740594A (zh) | 月球轨道飞行器表面红外辐射和反射辐射热流计算方法 | |
CN103488833B (zh) | 一种导航卫星复杂模型的热辐射力建模方法 | |
CN105868464A (zh) | 一种任意大型相控阵太阳能卫星聚光系统模型的构建方法 | |
Guo et al. | Simulations of reflected sun beam traces over a target plane for an azimuth–elevation tracking heliostat with fixed geometric error sources |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |