CN109711082B - 大椭圆冻结轨道卫星光照条件和帆板遮挡联合分析方法 - Google Patents
大椭圆冻结轨道卫星光照条件和帆板遮挡联合分析方法 Download PDFInfo
- Publication number
- CN109711082B CN109711082B CN201910017066.4A CN201910017066A CN109711082B CN 109711082 B CN109711082 B CN 109711082B CN 201910017066 A CN201910017066 A CN 201910017066A CN 109711082 B CN109711082 B CN 109711082B
- Authority
- CN
- China
- Prior art keywords
- sailboard
- satellite
- shielding
- calculating
- reference plane
- 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
Landscapes
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
Abstract
本发明提供了一种大椭圆冻结轨道卫星光照条件和卫星帆板遮挡联合分析方法,具体包括步骤:步骤1:根据光照条件对复杂外形的卫星进行包络化抽象处理;步骤2:对帆板进行网格化处理建模,判断帆板网格点和星体遮挡形成的阴影的关系,得到帆板被遮挡点的集合。本发明通过光照条件和遮挡的联合分析,获取卫星帆板被遮挡的计算结果,为卫星确定帆板面积和能源平衡计算提供依据,可实现复杂光照条件下批量计算卫星帆板遮挡。
Description
技术领域
本发明涉及卫星领域,具体地,涉及大椭圆冻结轨道卫星光照条件和帆板遮挡联合分析方法。尤其是涉及卫星总体方案设计中帆板面积确定问题,通过光照条件和遮挡的联合分析,获取卫星帆板被遮挡的计算结果,为卫星确定帆板面积和能源平衡计算提供依据。
背景技术
与常用的太阳同步轨道、地球静止轨道不同,在J2000.0坐标系下,大椭圆冻结轨道卫星轨道面每天自东向西进动。轨道面进动与太阳自西向东的运动复合构成了大椭圆轨道独特而复杂的光照条件。此外,由于卫星的模型复杂,若按照真实模型和全寿命期内光照条件进行帆板遮挡计算,将面临建模困难、计算量大等困难。因此需要一种复杂光照条件下批量计算卫星帆板遮挡的方法。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种大椭圆冻结轨道卫星光照条件和帆板遮挡联合分析方法。
根据本发明提供的一种大椭圆冻结轨道卫星光照条件和帆板遮挡联合分析方法,包括:
步骤1:根据光照条件对卫星进行包络化抽象处理,得到包络模型;
步骤2:对帆板进行网格化处理建模,得到帆板网格点;基于所述包络模型,判断帆板网格点与星体遮挡所形成阴影之间的关系,得到帆板被遮挡点的集合。
优选地,在步骤1中,通过计算全工况下卫星光照条件,并根据光照条件和卫星实际外形对卫星模型进行实时包络化抽象处理。
优选地,所述步骤1包括:
步骤1.1:获取计算帆板遮挡所需的参数;
步骤1.2:在XYZ空间直角坐标系中,以每一个时刻太阳至卫星矢量在XOY面的投影矢量为基准,实时建立一与该投影矢量垂直的基准平面,所述基准平面在XOY面内尺寸为遮光罩对角线长度,并且所述基准平面在XOY面内所处的位置距离原点的尺寸为卫星星体对角线长度;
其中,卫星星体朝向遮光罩的轴向方向为Z轴正向,卫星星体朝向一侧太阳帆板的轴向方向为Y轴正向,X轴正向分别垂直于Z轴正向、Y轴正向;帆板驱动轴安装点为XYZ空间直角坐标系的原点O。
优选地,所述计算帆板遮挡所需的参数,包括:遮光罩尺寸、星体尺寸、帆板尺寸、帆板连接架尺寸、帆板驱动轴安装点与遮光罩定点的距离尺寸、帆板驱动轴安装点与对天面的距离尺寸。
优选地,所述步骤2包括:
步骤2.1:根据太阳至卫星矢量与卫星Z轴方向的夹角,确定对地面遮挡或对天面遮挡,并根据夹角计算所述基准平面在XOY面内投影,包含基准平面的投影四边形各条边在XOY面内的表示方程;
步骤2.2:在XOY面对帆板作网格化处理,获得表征帆板的点阵集合;
步骤2.3:依次进行表征帆板的各点与基准平面投影的关系计算,判定该点是否处于基准平面投影内,获得帆板被遮挡的点阵集合;
步骤2.4:根据全寿命周期内光照条件,计算不同时刻卫星帆板被遮挡的面积。
根据本发明提供的一种存储有计算机程序的计算机可读存储介质,所述计算机程序被处理器执行时实现所述的大椭圆冻结轨道卫星光照条件和帆板遮挡联合分析方法的步骤。
与现有技术相比,本发明具有如下的有益效果:
本发明通过光照条件和遮挡的联合分析,获取卫星帆板被遮挡的计算结果,为卫星确定帆板面积和能源平衡计算提供依据,可实现复杂光照条件下批量计算卫星帆板遮挡。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为卫星模型的平视结构示意图。
图2为卫星模型的俯视结构示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
针对大椭圆冻结轨道卫星星体主要部件对帆板的遮挡阴影分析问题,本发明提供一种大椭圆冻结轨道卫星光照条件和帆板遮挡联合分析方法,包括:
步骤1:根据光照条件对复杂外形的卫星进行包络化抽象处理;
步骤2:对帆板进行网格化处理建模,判断帆板网格点和星体遮挡形成的阴影的关系,得到帆板被遮挡点的集合。
具体地,首先根据光照条件对复杂外形的卫星进行包络化抽象处理:
步骤1.1:获取计算帆板遮挡所需的各项参数,包括遮光罩尺寸S1和S8、星体尺寸S4和S9、帆板尺寸S6和S7、帆板连接架尺寸S5、帆板驱动轴安装点与遮光罩定点的距离尺寸S2、帆板驱动轴安装点与对天面的距离尺寸S3、如附图1和附图2;
步骤1.2:在如附图1和附图2所示的坐标系中,以每一个时刻太阳至卫星矢量在XOY面的投影矢量为基准,实时建立一与该投影矢量垂直的基准平面,该平面Z方向尺寸为S2,XOY面内尺寸为遮光罩对角线长度,并且该面在XOY面内所处的位置距离原点的尺寸为卫星星体对角线长度;其中,卫星星体朝向遮光罩的轴向方向为Z轴正向,卫星星体朝向一侧太阳帆板的轴向方向为Y轴正向,X轴正向分别垂直于Z轴正向、Y轴正向;帆板驱动轴安装点为XYZ空间直角坐标系的原点O。
然后对帆板进行网格化处理建模,判断帆板网格点和星体遮挡形成的阴影的关系,得到帆板被遮挡点的集合:
步骤2.1:根据太阳至卫星矢量与卫星Z方向的夹角确定对地面遮挡或对天面遮挡,并根据夹角计算步骤2.2中基准平面在XOY面内投影,包含投影四边形各条边在XOY面内的表示方程;
步骤2.2:如附图2,在XOY面对帆板作网格化处理,获得表征帆板的点阵集合;
步骤2.3:依次进行表征帆板的各点与基准平面投影的关系计算,判定该点是否处于基准平面投影内,获得帆板被遮挡的点阵集合;
步骤2.4:根据全寿命期内光照条件计算每时每刻卫星帆板被遮挡的面积。
以上对本发明的具体实施例进行了描述,根据实际计算,该方法计算出的遮挡与Pro-E建模计算实际模型的遮挡面积的差别在10%以内,具有较高的工程参考价值。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,如近似模型中尺寸的确定、基准平面的位置确定等,均可按照工程实际进行适应性确定,这并不影响本发明的实质内容。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (3)
1.一种大椭圆冻结轨道卫星光照条件和帆板遮挡联合分析方法,其特征在于:
具体包括如下步骤:
步骤1:根据光照条件对复杂外形的卫星进行包络化抽象处理;
步骤2:对帆板进行网格化处理建模,判断帆板网格点和星体遮挡形成的阴影的关系,得到帆板被遮挡点的集合;
所述步骤1包括:
步骤1.1:获取计算帆板遮挡所需的参数;
步骤1.2:在XYZ空间直角坐标系中,以每一个时刻太阳至卫星矢量在XOY面的投影矢量为基准,实时建立一与该投影矢量垂直的基准平面,所述基准平面在XOY面内尺寸为遮光罩对角线长度,并且所述基准平面在XOY面内所处的位置距离原点的尺寸为卫星星体对角线长度;
其中,卫星星体朝向遮光罩的轴向方向为Z轴正向,卫星星体朝向一侧太阳帆板的轴向方向为Y轴正向,X轴正向分别垂直于Z轴正向、Y轴正向;帆板驱动轴安装点为XYZ空间直角坐标系的原点O;
所述步骤2包括:
步骤2.1:根据太阳至卫星矢量与卫星Z轴方向的夹角,确定对地面遮挡或对天面遮挡,并根据夹角计算所述基准平面在XOY面内投影,包含基准平面的投影四边形各条边在XOY面内的表示方程;
步骤2.2:在XOY面对帆板作网格化处理,获得表征帆板的点阵集合;
步骤2.3:依次进行表征帆板的各点与基准平面投影的关系计算,判定各点是否处于基准平面投影内,获得帆板被遮挡的点阵集合;
步骤2.4:根据全寿命周期内光照条件,计算不同时刻卫星帆板被遮挡的面积。
2.根据权利要求1所述的大椭圆冻结轨道卫星光照条件和帆板遮挡联合分析方法,其特征在于,在步骤1中,通过计算全工况下卫星光照条件,并根据光照条件和卫星实际外形对卫星模型进行实时包络化抽象处理。
3.根据权利要求1所述的大椭圆冻结轨道卫星光照条件和帆板遮挡联合分析方法,其特征在于,所述计算帆板遮挡所需的参数,包括:遮光罩尺寸、星体尺寸、帆板尺寸、帆板连接架尺寸、帆板驱动轴安装点与遮光罩定点的距离尺寸、帆板驱动轴安装点与对天面的距离尺寸。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910017066.4A CN109711082B (zh) | 2019-01-08 | 2019-01-08 | 大椭圆冻结轨道卫星光照条件和帆板遮挡联合分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910017066.4A CN109711082B (zh) | 2019-01-08 | 2019-01-08 | 大椭圆冻结轨道卫星光照条件和帆板遮挡联合分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109711082A CN109711082A (zh) | 2019-05-03 |
CN109711082B true CN109711082B (zh) | 2023-08-08 |
Family
ID=66261098
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910017066.4A Active CN109711082B (zh) | 2019-01-08 | 2019-01-08 | 大椭圆冻结轨道卫星光照条件和帆板遮挡联合分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109711082B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111547273B (zh) * | 2020-05-14 | 2021-05-25 | 中国人民解放军国防科技大学 | 薄膜航天器 |
CN112613197B (zh) * | 2021-01-08 | 2021-08-24 | 中国科学院微小卫星创新研究院 | 一种星体对转动太阳帆板遮挡面积的分析方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102700728A (zh) * | 2012-06-18 | 2012-10-03 | 上海卫星工程研究所 | 一种卫星太阳电池阵在轨被遮挡的确定方法及其应用 |
CN103150414A (zh) * | 2012-06-13 | 2013-06-12 | 中国空间技术研究院 | 基于stl的通信卫星大天线对太阳翼遮挡分析方法 |
CN104615841A (zh) * | 2015-03-05 | 2015-05-13 | 哈尔滨工业大学 | 考虑遮挡效应的航天器太阳能帆板三维动态仿真方法 |
CN106251335A (zh) * | 2016-07-19 | 2016-12-21 | 中国空间技术研究院 | 一种基于stl网格化模型的敏感器视场遮挡区域确定方法 |
CN108920829A (zh) * | 2018-06-29 | 2018-11-30 | 中国空间技术研究院 | 一种带大型网状天线卫星的太阳光压力矩计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107108047A (zh) * | 2014-08-26 | 2017-08-29 | 有效空间解决方案有限公司 | 用于卫星的对接系统和方法 |
-
2019
- 2019-01-08 CN CN201910017066.4A patent/CN109711082B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103150414A (zh) * | 2012-06-13 | 2013-06-12 | 中国空间技术研究院 | 基于stl的通信卫星大天线对太阳翼遮挡分析方法 |
CN102700728A (zh) * | 2012-06-18 | 2012-10-03 | 上海卫星工程研究所 | 一种卫星太阳电池阵在轨被遮挡的确定方法及其应用 |
CN104615841A (zh) * | 2015-03-05 | 2015-05-13 | 哈尔滨工业大学 | 考虑遮挡效应的航天器太阳能帆板三维动态仿真方法 |
CN106251335A (zh) * | 2016-07-19 | 2016-12-21 | 中国空间技术研究院 | 一种基于stl网格化模型的敏感器视场遮挡区域确定方法 |
CN108920829A (zh) * | 2018-06-29 | 2018-11-30 | 中国空间技术研究院 | 一种带大型网状天线卫星的太阳光压力矩计算方法 |
Non-Patent Citations (1)
Title |
---|
一种组合体航天器帆板的光照遮挡计算新方法;黄岸毅等;《载人航天》;20130930;正文第2-4节 * |
Also Published As
Publication number | Publication date |
---|---|
CN109711082A (zh) | 2019-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109711082B (zh) | 大椭圆冻结轨道卫星光照条件和帆板遮挡联合分析方法 | |
CN108920829B (zh) | 一种带大型网状天线卫星的太阳光压力矩计算方法 | |
CN109459043B (zh) | 一种基于生成式重构图像的航天器相对导航方法 | |
CN107985631B (zh) | 低轨微纳卫星及适用于脉冲微弧电推力器的在轨安装方法 | |
Zhang et al. | Task-driven on-Board Real-Time Panchromatic Multispectral Fusion Processing Approach for High-Resolution Optical Remote Sensing Satellite | |
CN114840812A (zh) | 欧拉角的计算方法、装置、设备及介质 | |
CN113091753B (zh) | 用于星敏视场保护的卫星姿态导引方法及其系统 | |
CN109975836A (zh) | 解算ccd图像地面位置方法、装置、电子设备及介质 | |
CN113674391B (zh) | 一种升降轨成像双星地物表观辐亮度差异计算方法 | |
Zhang et al. | Optical image simulation system for space surveillance | |
CN112857306B (zh) | 一种视频卫星任意视向点的连续太阳高度角确定方法 | |
CN103679780A (zh) | 一种空间目标实时模拟方法 | |
Bechini et al. | Spacecraft pose estimation via monocular image processing: Dataset generation and validation | |
CN114485620A (zh) | 融合轨道动力学的小行星探测器自主视觉定位系统及方法 | |
Grigoriev et al. | Unified algorithm of cone optics to compute solar flux on central receiver | |
Fasano et al. | A stereo-vision based system for autonomous navigation of an in-orbit servicing platform | |
CN108520150B (zh) | 一种航天器含金属网布的构件的透光性遮挡的计算方法 | |
CN111366162B (zh) | 基于太阳帆板投影与模板匹配的小天体探测器位姿估计方法 | |
Inose et al. | Refining outdoor photometric stereo based on sky model | |
CN108413903B (zh) | 一种空间目标过境区域位置检测方法及系统 | |
CN111536982B (zh) | 一种光照不良条件下小行星探测器光学导航图像生成方法 | |
Nugnes et al. | A New Analytical Method for Eclipse Entry/Exit Positions Determination Considering a Conical Shadow and an Oblate Earth Surface | |
CN115801108B (zh) | 一种卫星星座分析方法 | |
Chen et al. | Extending Legacy Climate Models by Adaptive Mesh Refinement for Single Component Tracer Transport | |
Chen et al. | Spacecraft Homography Pose Estimation with Single-Stage Deep Convolutional Neural Network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |