基于增广拉格朗日乘子法的全电推进卫星多学科优化.docxVIP

基于增广拉格朗日乘子法的全电推进卫星多学科优化 0 全电推进卫星系统设计方法的提出 与传统的化学推进卫星相比，全球电推高压电气发射卫星能够执行地球同步轨迹的转移、地球静态轨迹的保护和生命末的离轨，从而显著减少施工量。“一箭双星”的发射对于降低卫星的开发成本和提高卫星的有效载荷比具有重要意义。 由于电推进系统推力较小 (几十至几百毫牛) , 全电推进卫星地球同步轨道转移时间较传统化学推进卫星大幅延长 (一般需3至6个月) , 较长的变轨周期还导致卫星长期运行在辐射带中, 造成太阳翼性能退化与输出功率降低, 对卫星抗辐射能力与供配电系统等方面的设计提出了新的技术挑战。由于全电推进卫星系统设计是一个多学科耦合的复杂工程问题, 有必要采用多学科设计优化 (Multidisciplinary design optimization, MDO) 方法 然而, 在航天器系统实际设计中需要调用高精度仿真分析模型 (如有限元模型等) 以提高设计可信度, 导致了优化成本急剧增加。此外, 航天器MDO问题往往需要通过多学科分析 (Multidisciplinary analysis, MDA) 过程进行迭代求解, 从而进一步加剧计算复杂性 常见的SBAO方法包括高效全局优化方法 1 基于mdo模型的城市轨道系统设计 考虑到全电推进卫星推力小、变轨时间长以及辐射环境恶劣等特点, 本文参考文献[17], 主要考虑轨道转移、位置保持、空间环境、供配电、结构及质量学科建立系统MDO模型, 其设计结构矩阵 (DSM) 如图1所示, 图中各参数定义如表1所示。下面对各学科建模方法进行简要说明。 1.1 小推力轨道转移动力学方程 轨道转移学科分析的目的是获取转移时间T 由于电推进系统推力较小, 与卫星所受到的摄动加速度处于同一量级, 可将电推进加速度作为摄动加速度处理, 建立小推力轨道转移动力学方程如下所示 式中:各参数定义可参考文献[18]。转移过程中保持推力大小不变。 卫星初始位于近地点几百千米, 远地点位于地球同步轨道的大椭圆GTO轨道。本文采用推力方向固定的两阶段变轨策略 在轨道学科需满足的约束条件如下所示: 1.2 不点面和推力器安装位置 位置保持学科分析的目的是获取东西方向位保精度l 本文基于文献[19]中的位保策略, 建立全电推进卫星完备状态下的小推力位置保持模型。卫星上共装有四个推力器, 关于质心对称安装于背地板上, 且推力方向通过质心, 其位置定义在TNR坐标系中。TNR坐标系及推力器安装位置如图3所示。 由于GEO轨道倾角和偏心率均为零, 为避免动力学方程奇异, 本文采用无奇点轨道要素描述GEO轨道, GEO轨道脉冲推力控制方程参考文献[19]。 本文以2天为一个位保小周期, 7个小周期构成一个完整的位保周期, 总位保时间为15年。每个小周期中每个推力器仅开机一次, 顺序为推力器1-推力器4-推力器2-推力器3, 其中推力器1和推力器2开机赤经相同, 推力器3和推力器4开机赤经相同, 两组开机赤经相差180°, 根据脉冲推力控制方程和推力矢量余弦, 获得小控制周期内各推力器产生的速度增量, 继而得到该小控制周期内卫星的位移。 位置保持学科需满足的约束条件如下所示: 1.3 空间环境学科 空间环境学科分析的目的是获取太阳翼功率下降系数p 式中:D 1.4 供能源科学 供配电学科分析的目的分为三部分。第一部分是根据太阳翼功率下降系数, 计算太阳翼寿命初期输出功率P 供配电学科需满足的约束条件如下所示: 1.5 卫星结构分析 结构学科分析的目的是获取卫星结构质量M 本文研究的全电推进卫星采用中心承力筒结构。为减小模型复杂程度, 提高设计效率, 本学科采用文献[22]中的经验公式进行卫星结构分析。卫星轴向频率f 式中:E为杨氏模量, 取值为7.06×10 本文中卫星采用复合材料, 包括一层芯子及4层铺层, 芯子密度为30.95 kg/m 结构学科需满足的约束条件如下所示: 1.6 质量分析 质量分析学科主要目的是获取卫星整星质量M 式中:M 2 alm-ego算法流程 针对现有大部分SBAO方法 (如EGO, MPS等) 无法处理复杂工程问题中的高耗时约束的问题, 本文提出一种基于增广拉格朗日乘子法的高效全局优化方法 (Augmented Lagrange multiplier based efficient global optimization, ALM-EGO) 来实现高效求解包含高耗时约束的全电推进卫星MDO问题。EGO以Kriging为代理模型, 并引入期望改善度 (Expected improvement, EI) 作为采样准则, 通过在EI最大值处新增样本点实现对代理模型的更新, 从而引导算法快速收敛, 具体流程参考文献[13]。 传