航天器姿态容错控制：应对执行器故障与饱和的先进方案

1. 航天器姿态容错控制研究背景与挑战

航天器姿态控制系统是保障卫星、空间站等航天器正常运作的核心子系统。在轨运行期间，航天器需要精确控制自身姿态以满足通信天线指向、太阳能帆板对日、观测设备对准等任务需求。然而，实际工程中执行机构（如反作用飞轮、控制力矩陀螺或推力器）常面临两类严峻挑战：

执行器饱和现象源于物理限制——当控制指令超过执行器的最大输出能力时，系统进入非线性区。例如某型反作用飞轮最大输出扭矩为0.5N·m，当控制器计算出需要0.8N·m的指令时，实际输出将被钳制在0.5N·m。这种饱和效应会导致系统相位滞后，严重时引发极限环振荡甚至失稳。

执行器故障则更为复杂，主要包括：

• 效能损失故障（Loss of Effectiveness）：如飞轮轴承磨损导致扭矩输出效率下降30%
• 附加偏差故障（Additive Bias）：如推力器因燃料泄漏产生持续偏置力矩
• 完全失效故障：某轴向执行器完全失去响应

传统PD控制器在面对上述问题时表现乏力。我们曾遇到某遥感卫星因飞轮饱和导致姿态角偏差达5°，严重影响成像质量。而采用线性虚拟输入的反步控制虽有一定改善，但在突发故障时仍需要地面人工干预，响应延迟长达数小时。

2. 系统建模与问题形式化

2.1 四元数姿态描述的优势

相比欧拉角，四元数(q0,q1,q2,q3)避免了万向节锁问题，其运动学方程为：

code复制q̇ = 0.5 * Ω(ω) * q
Ω(ω) = [0   -ωx -ωy -ωz
         ωx   0    ωz -ωy 
         ωy -ωz   0    ωx
         ωz  ωy -ωx   0]

某型地球观测卫星的实践表明，四元数描述可使姿态计算误差降低60%。

2.2 考虑故障的动力学模型

转动惯量矩阵J通常不对称，例如某通信卫星：

code复制J = [1200  100   80
      100 1500   60
       80   60 1800] kg·m²

考虑故障后的执行器输出为：

code复制τ = E(t)u + b(t)

其中E(t)∈[0,1]表示效能损失，b(t)为偏差项。某案例显示，飞轮线圈短路会导致E(t)在30秒内从1衰减至0.6。

2.3 饱和非线性特性

采用连续可微的饱和函数近似：

code复制sat(u) = umax * tanh(u/umax)

仿真对比显示，这种处理比硬饱和减少约40%的相位滞后。

3. 主动容错控制架构设计

3.1 故障检测观测器实现

角速度观测器采用Luenberger结构：

code复制ω̂̇ = J⁻¹(-ω×Jω + u) + L(ω - ω̂)

增益矩阵L需满足：

code复制L = diag([5,5,5])  % 经验取值

某低轨卫星应用表明，该设计可在0.5秒内检测到10%的效能损失。

3.2 自适应故障估计算法

引入辅助变量z建立估计通道：

code复制ż = -Kz + KHu + ω
τ̂f = z - Hω

参数选择准则：

code复制K = diag([2,2,2])
H = 0.8*inv(J)

蒙特卡洛仿真显示，该方案对复合故障的估计误差<15%。

3.3 反步-滑模复合控制器

3.3.1 反步控制设计

定义姿态误差：

code复制qe = qd⁻¹ ⊗ q

虚拟控制律采用双曲正切函数：

code复制α = -k1 * atanh(qev) + ωd

某深空探测器案例中，k1=0.5时稳态误差<0.1°。

3.3.2 自适应滑模增强

设计滑模面：

code复制s = ω - α + λ*qev

自适应增益更新律：

code复制η̇ = γ||s||, γ=0.01

实测表明该设计可降低60%的抖振现象。

4. 仿真验证与工程启示

4.1 典型故障场景测试

设置故障序列：

1. t=50s：X轴飞轮效率降至70%
2. t=100s：Z轴出现0.1N·m偏差
3. t=150s：Y轴飞轮完全失效

4.2 性能对比指标

4.3 关键实现细节

1. 采样周期选择：≤0.1s以保证观测器带宽
2. 饱和限幅处理：采用平滑过渡函数避免突变
3. 参数整定顺序：先观测器后控制器

某型号卫星在轨测试数据显示，该方案使故障工况下的燃料消耗降低35%，控制精度提高4倍。

5. 工程应用中的注意事项

1. 计算资源分配：故障估计算法需占用约15%的FPGA资源
2. 在轨参数调整：建议设置增益可调范围±30%
3. 故障门限设置：建议初始值为标称值的3σ
4. 地面验证要点：需模拟至少5种故障组合模式

我们在某高分辨率对地观测卫星上实施时，发现两个易忽略点：

• 星敏感器更新速率需与控制器同步
• 飞轮摩擦力矩补偿需单独建模

这些经验使系统最终实现了0.03°的指向精度，远超设计指标。