[发明专利]一种基于二次反射的高精度卫星太阳光压摄动力建模方法

权利要求书

1.一种基于二次反射的卫星太阳光压摄动力建模方法，其特征在于，包括如下步骤：

步骤1、获取待建模卫星的表面三维模型，确定卫星的所有表面部件，将每个表面部件的表面划分成多个面元；

步骤2、用面积等于卫星外形最大包络的方形像元阵列作为模拟太阳光源，每个像元作为一个小光源并发出一条垂直于方形像元阵列的太阳光线，所述模拟太阳光源与卫星的距离为一个天文单位；

步骤3、针对所述卫星的各个表面部件，记录所述模拟太阳光源的每个像元发出的光线与表面部件最短距离的交点和相应的入射角，统计第i个表面部件上有所述交点的面元数量，即为初次照射第i个表面部件的有效被照射面元数量A_ieff1；i＝1,2,...N，其中N为表面部件数量；

根据有效被照射面元数量A_ieff1，得到第i个表面部件由模拟太阳光源初次照射所受的法向光压摄动力f_ni1以及切向光压摄动力f_si1；

步骤4、首先，针对各个表面部件，根据步骤3获得的所有入射光线的交点和入射角，继而确定入射光线在所在面元上的反射光线；

其次，统计第i个表面部件与所述反射光线有最短距离交点的面元个数A_ieff2，则第i个表面部件被二次照射中漫反射的有效面元数为A_ieff2·v_i(1-μ_i)，镜面反射的有效面元数为A_ieff2·v_iμ_i，其中，ν_i表示第i个表面部件的反射率、μ_i第i个表面部件的镜面系数；

然后，根据A_ieff2·v_i(1-μ_i)，获得第i个表面部件在二次照射中由漫反射引起的法向光压摄动力f_ni2d；根据A_ieff2·v_iμ_i，获得第i个表面部件在二次照射中由镜面反射引起的法向光压摄动力f_ni2m和切向光压摄动力f_si2m；

最后，得到第i个部件所受的总的法向光压摄动力f_ni为：

f_ni＝f_ni1+f_ni2d+f_ni2m；

总的切向光压摄动力f_si为：

f_si＝f_si1+f_si2m；

步骤5、根据各个表面部件的安装位置和法线方向，将所有表面部件的法向光压摄动力f_ni和切向光压摄动力f_si分别分解到卫星本体坐标系的三个轴方向，并得到各轴方向的合力，最后得到各轴方向的光压摄动加速度；

步骤6、采用步骤2-5的方法，在一个卫星运动周期内，以设定角度为步长，计算设定时间内每个卫星运动周期不同时刻的太阳光压摄动加速度；针对计算结果，通过傅里叶多项式拟合得到太阳光压摄动加速度的数学模型。

2.如权利要求1所述的一种基于二次反射的卫星太阳光压摄动力建模方法，其特征在于，所述步骤6中的设定时间为半年。

3.如权利要求1所述的一种基于二次反射的卫星太阳光压摄动力建模方法，其特征在于，所述像元大小为1mm×1mm。

4.如权利要求1所述的一种基于二次反射的卫星太阳光压摄动力建模方法，其特征在于，所述面元大小为1mm×1mm。