【技术实现步骤摘要】
本专利技术实例涉及航天器动力学领域，尤其涉及一类堆栈平板式卫星的刚柔耦合姿态动力学的变分积分方法。
技术介绍
1、近年来卫星互联网星座技术的快速发展，以轻质化、高功能密度化及柔性化为特征的空间系统逐渐成为卫星设计的重要趋势。其中相比于陆基信号基站，用以对地提供网络直连服务的低轨卫星距终端设备较远，立方体卫星通常使用的传统尺寸天线几乎不可能从地面移动终端接收到信号，而模块化设计的超大型相控阵天线阵列可突破这一物理约束。此外卫星天线尺寸与在卫星通信信道中可实现的数据吞吐量直接相关，除增益之外，增加卫星天线尺寸还可提高空间分辨率，从而在给定的服务区域内实现更大的频率复用。因此在同低功率、低信号强度的地面移动终端进行通信时，部署更大的相控阵天线是建设未来太空基站，实现对地移动终端网络直连的必要途径，比如于2022年由美国ast&science研发的bluewalker-3大型手机直连卫星具备近1500kg的质量，其所搭载的相控阵天线达到了64m2。另一方面随着此类卫星电子设备的升级，星载的太阳能阵列在整个卫星的质量、尺寸所占的比重越来越大，因此搭载大功率电子设备的卫星的太阳能阵列通常具有较大的翼展。2019年以来，亚马逊、telesat和spacex等航天企业正致力于部署巨型星座来为全球提供宽带互联网，其中spacex公司所构建的starlink卫星互联网星座中单星的结构设计上采用了平板式的扁平化设计方案，初代的整星产品质量约为240～260kg，太阳翼面积达到24m2，而第二代的产品中也将具备更大的质量和翼展面积。近期我国也开
2、堆栈平板式阵列卫星的大翼展往往具备一定程度的柔性和较低的固有频率。保证卫星动力学高保真度的关键在于对易弯曲或振动的附件进行建模，如相控阵天线阵列和太阳能阵列。此类大尺寸附件在复杂空间环境中易于弯曲，多数情况下，柔性附件的弯曲需要包含在整星的动力学模型中。初始建模过程中，一般假设卫星为刚体，但若存在柔性附件，该假设会降低数值仿真的保真度。尽管现阶段已有多种柔性体离散化方法，包括假设模态法、有限元法、有限段法和无网格法，但此类动力学模型维数通常较高，通常计算量大、耗时长。
3、工程应用中通常采用单位四元数来表征姿态，其数值积分过程中一般在每一积分步中通过归一化来施加单位长度的约束，这种方法并不利于几何结构的保持和动量、能量的守恒。通过基于旋转矩阵的姿态表征方式和变分积分器可以解决这一问题，但同时引入了过多的姿态表征参数，占用了较多的计算资源，算法的实现成本较高。
技术实现思路
1、针对现有堆栈平板式阵列卫星采用旋转矩阵的姿态表征方式存在冗余表征的问题，本专利技术提供一种针对堆栈平板式卫星的刚柔耦合姿态动力学建模方法。该方法使用单位四元数推导了用于堆栈平板式卫星的李群变分积分器，该积分器全局有效，避免了旋转矩阵的冗余表征弊端，此外，这种方式推导的积分法保留了连续系统的许多特性，例如动量和能量守恒。
2、本专利技术所述一种针对堆栈平板式卫星的刚柔耦合姿态动力学建模方法，该方法包括以下步骤：
3、步骤1、建立惯性坐标系卫星体坐标系和板式附件浮动坐标系
4、步骤2、根据各坐标系分别建立卫星各附件质量元素相对卫星体坐标系的位置、速度关系；
5、步骤3、根据整星姿态运动学方程及板式附件的振动姿态运动学方程获取整星姿态四元数变分、振动姿态四元数变分、卫星体坐标系下整星的质心位置向量、整星质心位置向量的一阶导数及整星质心位置向量的变分；
6、步骤4、获取整星的动能、势能及其变分，并根据实际非保守力、力矩获取虚功的表达式；
7、步骤5、根据连续系统的largrange量获取离散系统的largrange量；
8、步骤6、根据离散legendre变换建立整星离散hamilton动力学方程，完成堆栈平板式卫星的刚柔耦合姿态动力学建模。
9、优选地，步骤1中的卫星体坐标系的原点位于卫星轮毂上的任意位置；
10、附件浮动坐标系的原点位于板式附件与轮毂的铰接点处，三轴xqyqzq与附件固连，轴指向翼展伸展方向，轴垂直于翼平面，轴由右手定则给定，附件浮动坐标系用于表征卫星板式附件偏离初始安装姿态的振动姿态。
11、优选地，步骤2包括轮毂上任一质量元素的位置向量、速度和板式附件上任一质量元素的位置向量、速度；
12、其中，轮毂上任一质量元素的位置向量rh、速度vh按下式获取：
13、
14、式中为卫星体坐标系原点至轮毂质心的向量，在系下表示；r为卫星相对惯性坐标系的姿态矩阵，ω为卫星相对惯性坐标系的角速度；ρ为轮毂质心至该质量元素的位置向量；
15、板式附件上任一质量元素的位置向量rsa、速度vsa按下式获取：
16、
17、式中rh为板式附件相对卫星体坐标系的初始安装姿态；d为卫星体坐标系原点至板式附件浮动坐标系原点的向量，在系下表示；s为板式附件浮动坐标系原点至板式附件质心的向量；rsa为板式附件偏离其安装姿态的振动姿态；ρsa为板式附件质心至板式附件上该质量元素的位置向量；θ为板式附件绕轴的振动角；为板式附件绕轴的振动角速度；e2表示轴单位向量，
18、优选地，步骤3具体包括：
19、步骤31、根据整星姿态运动学方程获取整星姿态四元数变分；
20、整星姿态运动学方程为：
21、
22、式中表示四元数乘法，q表示整星的姿态四元数，表示整星的姿态四元数导数，表示整星角速度的增广纯四元数表示，且
23、根据整星姿态运动学方程获取整星四元数变分为：
24、
25、式中δ[·]表示·的变分，ε表示无穷小变化量，表示四元数指数，表示η的增广纯四元数，η表示姿态无穷小变分；
26、步骤32、根据板式附件的振动姿态运动学方程获取振动姿态四元数变分；
27、板式附件的振动姿态运动学方程为：
28、
29、式中qsa表示板式附件相对卫星体坐标系的振动姿态四元数，表示板式附件相对卫星体坐标系的振动姿态四元数导数，表示轴单位向量的增广纯四元数；
30、根据板式附件的振动姿态运动学方程获取振动姿态四元数变分为：
31、
32、根据可推导得到角速度变分：
33、
34、步骤33、获取卫星体坐标系下整星的质心位置向量csc：
35、
36、式中，msa、mh、msc分别为板式附件质量、轮毂质量和整星质量；
37、步骤34、根据板式附件的振动姿态运动学方程计算整星的质心位置一阶导数：
38、
39、[·]^表示3维向量到3×3斜对称矩阵的映射；本文档来自技高网...
【技术保护点】
1.一种针对堆栈平板式卫星的刚柔耦合姿态动力学建模方法，其特征在于，该方法包括以下步骤：
2.根据权利要求1所述一种针对堆栈平板式卫星的刚柔耦合姿态动力学建模方法，其特征在于，步骤1中的卫星体坐标系的原点位于卫星轮毂上的任意位置；
3.根据权利要求2所述一种针对堆栈平板式卫星的刚柔耦合姿态动力学建模方法，其特征在于，步骤2包括轮毂上任一质量元素的位置向量、速度和板式附件上任一质量元素的位置向量、速度；
4.根据权利要求3所述一种针对堆栈平板式卫星的刚柔耦合姿态动力学建模方法，其特征在于，步骤3具体包括：
5.根据权利要求4所述一种针对堆栈平板式卫星的刚柔耦合姿态动力学建模方法，其特征在于，步骤4具体包括：
6.根据权利要求5所述一种针对堆栈平板式卫星的刚柔耦合姿态动力学建模方法，其特征在于，步骤5具体包括：
【技术特征摘要】
1.一种针对堆栈平板式卫星的刚柔耦合姿态动力学建模方法，其特征在于，该方法包括以下步骤：
2.根据权利要求1所述一种针对堆栈平板式卫星的刚柔耦合姿态动力学建模方法，其特征在于，步骤1中的卫星体坐标系的原点位于卫星轮毂上的任意位置；
3.根据权利要求2所述一种针对堆栈平板式卫星的刚柔耦合姿态动力学建模方法，其特征在于，步骤2包括轮毂上任一质量元素的位置向量、速度和...
【专利技术属性】
技术研发人员：李化义，曹芊，李冬柏，吴凡，张冀鹞，何文韬，
申请(专利权)人：哈尔滨工业大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。