4. 动量轮工作原理与理想模型

动量轮是航天器姿态控制中最重要的执行机构之一。它本质上是一个由电机驱动的飞轮，通过改变飞轮自身的角动量，与航天器本体进行角动量交换，从而产生作用于星体的反作用力矩。理解动量轮，需要抓住三个层次：角动量交换的物理本质、反作用轮与偏置动量轮在构型与工作点上的差异，以及理想模型下力矩—转速关系与多轮构型矩阵的数学表达。

一、角动量交换原理

设飞轮绕其自转轴的转动惯量为 J，相对星体的角速度为 Ω，则飞轮相对星体的角动量为 h = JΩ。当飞轮与星体固连、星体本身以角速度 ω 相对惯性空间转动时，飞轮相对惯性空间的角动量为 J(Ω + ω)。对整个航天器（星体 + 飞轮）应用角动量定理，若忽略外部力矩，系统总角动量守恒。飞轮电机施加控制力矩 T 改变 Ω 时，星体必然受到大小相等、方向相反的力矩，即

T_body = −T_wheel = −J dΩ/dt

这就是角动量交换的核心：飞轮加速，星体获得反向力矩；飞轮减速，星体获得正向力矩。飞轮并不“消耗”角动量，而是把角动量在飞轮与星体之间来回搬运。因此动量轮属于内力矩执行机构，它不能改变系统总角动量，只能改变角动量在内部的分配。这一点与喷气推力器有本质区别：推力器通过喷出质量产生外力矩，可以真正卸载系统角动量。

由此也带来动量轮的根本约束：飞轮转速不能无限增加。当飞轮接近饱和转速 Ω_max 时，它已无法继续沿该方向加速，必须借助磁力矩器或推力器等外部力矩进行“卸载”，把积累的角动量排出去，飞轮才能回到工作区间。这是动量轮系统设计中必须处理的饱和与卸载问题。

二、反作用轮与偏置动量轮的差异

按标称工作点不同，动量轮分为两类。反作用轮（Reaction Wheel, RW）的标称转速接近零，转速可正可负，输出力矩方向可正可反，因此能提供双向控制力矩，适合三轴零动量控制。偏置动量轮（Momentum Wheel, MW）的标称转速远离零，通常维持一个较大的恒定偏置角动量，飞轮主要在该偏置点附近小幅调速，利用陀螺效应提供姿态稳定，同时通过转速变化产生控制力矩。

两者的差异可以从工作点、力矩能力、饱和特性和适用场景几个方面对比。反作用轮工作点低，转速穿越零点时摩擦力矩和死区会影响低速性能，但控制灵活、可双向出力；偏置动量轮工作点高，飞轮始终单向旋转，避开了零速死区，且偏置角动量能提供被动陀螺刚度，但只能在小范围内双向调速，饱和方向不对称，且偏置角动量会与星体运动耦合。工程上，三轴稳定卫星常用三个正交反作用轮加一个斜装备份轮，而某些对成本敏感或需要长期稳定的平台会采用偏置动量轮加偏置动量管理方案。

三、理想力矩—转速关系

在理想模型中，忽略摩擦、风阻、电机电感与电流环动态，只保留飞轮转动惯量与电机力矩的关系。设飞轮角动量为 h = JΩ，电机输出力矩 T 满足

T = dh/dt = J dΩ/dt

对星体而言，其受到的反作用力矩为 −T。若把飞轮视为一个理想力矩源，则其输入是力矩指令 T_cmd，输出是飞轮角动量变化率。理想模型下，飞轮可以瞬时输出任意大小（在饱和范围内）的力矩，转速则是对力矩的积分：

Ω(t) = Ω(0) + (1/J) ∫₀^t T(τ) dτ

这个积分关系说明：飞轮转速是力矩的时间累积。持续输出同向力矩，转速会单调上升直至饱和；要维持转速不漂移，力矩指令的长期均值必须为零。实际工程中，飞轮还存在最大力矩 T_max、最大转速 Ω_max、最大角动量 h_max = JΩ_max 等约束，理想模型是这些约束下进行控制分配的基础。

四、多轮构型矩阵

单轴飞轮只能产生绕自身转轴的力矩。要实现三轴姿态控制，需要多个飞轮按一定几何构型安装。设第 i 个飞轮的转轴在星体坐标系中的单位方向矢量为 a_i = [a_ix, a_iy, a_iz]^T，飞轮角动量为 h_i，则所有飞轮对星体的总角动量为

H_w = Σ h_i a_i = A h

其中 h = [h₁, h₂, …, h_n]^T 为各轮角动量组成的列向量，A 为 3×n 的构型矩阵（也称安装矩阵），其第 i 列为 a_i。对时间求导，得到力矩关系：

T_body = −A T_w

式中 T_w = [T₁, …, T_n]^T 为各轮电机力矩。控制分配问题就是已知期望的三轴控制力矩 T_cmd，求解各轮力矩指令 T_w，使 A T_w = −T_cmd。当 n > 3 时方程欠定，通常取最小范数解 T_w = −A^T(AA^T)⁻¹T_cmd，在满足指令的同时使各轮力矩分配最均衡。

常见的构型包括：三正交轮（A = I₃，结构简单但无冗余）、三正交加一斜装轮（第四个轮轴与三个正交轴等角，提供故障冗余）、四轮金字塔构型（四个轮轴指向正四面体顶点，任意三轮可覆盖三轴，冗余度高）、以及偏置动量轮加反作用轮的混合构型。构型矩阵的秩决定了可控自由度：rank(A) = 3 时三轴完全可控；若某方向出现奇异，则该方向失去控制能力，需要靠其他执行机构补足。

以一个四轮金字塔构型为例，设四个轮轴与星体 +Z 轴夹角均为 θ，方位角分别为 0°、90°、180°、270°，则第 i 个轮轴方向为 a_i = [sinθ cosφ_i, sinθ sinφ_i, cosθ]^T。取 θ = 54.7356°（即 cosθ = 1/√3），四个轮轴恰好指向正四面体顶点，此时 AA^T = (4/3)I₃，条件数为 1，力矩分配最均衡，且任意三轮组合仍可张成三维空间，具备单轮失效后的重构能力。这一构型在工程中应用广泛，也是理解构型矩阵与冗余设计关系的典型例子。

综上，动量轮建模的核心链条是：角动量交换给出内力矩本质，工作点选择区分反作用轮与偏置动量轮，理想力矩—转速关系给出单轮动态，构型矩阵把多个单轮统一到三轴控制分配框架中。掌握这四步，就为后续加入摩擦、饱和、死区等真实约束的建模打下了基础。