3、动量轮工作原理：动量轮结构、动量交换原理、动量轮数学模型

动量轮是卫星姿态控制系统中最为核心的执行机构之一。它通过内部高速旋转的飞轮转子储存角动量，并与卫星本体进行动量交换，从而实现对卫星姿态的精确控制。理解动量轮的结构组成、动量交换的物理本质以及其数学模型，是后续掌握卸载策略、进行故障诊断和工程调优的基础。本章将从这三个维度展开，并结合工程实际，剖析动量轮在轨工作的核心逻辑。

3.1 动量轮结构：从机械本体到控制电子学的完整组成

动量轮在工程上是一个高度集成的机电产品，其结构设计直接决定了动量轮的寿命、控制精度和可靠性。一个典型的动量轮系统由机械结构、电机驱动、测量反馈和润滑密封四大子系统构成。

3.1.1 核心机械部件

• 飞轮转子：通常采用高比强度材料（如钛合金、铍或碳纤维复合材料）制成，设计为圆盘或轮缘结构。其转动惯量 J_w 是动量轮的关键设计参数，决定了在相同转速下储存角动量的能力。例如，一个转动惯量为 0.05 kg·m² 的飞轮，在 6000 rpm（约 628 rad/s）转速下可储存 31.4 N·m·s 的角动量。
• 驱动电机：绝大多数动量轮采用无刷直流电机（BLDC），因其具有高效率、长寿命、低噪声和易于控制等优点。电机定子绕组产生旋转磁场，驱动永磁体转子旋转。电机的力矩常数 K_t 和反电动势常数 K_e 是关键的电气参数。
• 支撑轴承：轴承是动量轮寿命的瓶颈。传统动量轮使用精密角接触球轴承，配合油润滑或脂润滑。长寿命动量轮则采用主动磁悬浮轴承，消除机械接触，但增加了控制复杂度。轴承的摩擦力矩 T_f 直接影响动量轮的功耗和控制精度。
• 测速装置：通常采用光电编码器或旋转变压器，实时测量飞轮转速。转速测量的精度和更新率决定了转速闭环控制的带宽和精度。高精度动量轮要求测速分辨率优于 0.1 rpm。

3.1.2 控制电子学

动量轮内部集成了电机驱动电路和转速控制电路。驱动电路将母线电压转换为三相桥式PWM驱动信号，控制电机电流。控制电路接收上位机（卫星姿控计算机）发送的转速指令或力矩指令，通过电流环和转速环双闭环实现精确控制。典型的动量轮电子学架构包括：电源模块（DC-DC）、数字信号处理器（DSP或MCU）、功率驱动级（MOSFET H桥）、电流采样电路和通信接口（通常为RS-422或CAN总线）。

3.2 动量交换原理：角动量守恒的工程应用

动量轮工作的物理基础是角动量守恒定律。在没有外力矩作用的理想情况下，卫星与动量轮组成的系统总角动量保持恒定。当动量轮电机加速或减速时，飞轮的角动量发生变化，为了维持系统总角动量不变，卫星本体必然产生反向的角速度变化，从而实现姿态控制。

3.2.1 基本原理推导

设卫星本体的转动惯量为 J_s，角速度为 ω；动量轮飞轮的转动惯量为 J_w，相对卫星的转速为 Ω。则系统总角动量为：

H_total = J_s · ω + J_w · (ω + Ω) = (J_s + J_w) · ω + J_w · Ω

当动量轮电机施加力矩 T_w 使飞轮转速变化时，根据角动量守恒（忽略外力矩）：

dH_total/dt = 0 → (J_s + J_w) · dω/dt = -J_w · dΩ/dt

由此可得卫星角加速度与飞轮角加速度的关系：

α_s = -[J_w / (J_s + J_w)] · α_w

其中 α_s 为卫星角加速度，α_w 为飞轮角加速度。可见，飞轮加速时卫星反向转动，飞轮减速时卫星同向转动。动量轮输出的控制力矩为：

T_control = J_w · dΩ/dt

3.2.2 动量交换的工程意义

动量交换机制使得卫星能够在不消耗推进剂的情况下进行姿态机动和稳定控制。相比推力器，动量轮具有以下优势：

• 高精度：通过精确控制电机电流，可实现微牛·米量级的力矩输出，满足高精度对地观测卫星的指向要求（如 0.001° 级）。
• 无污染：不产生羽流，不会污染光学载荷或太阳帆板。
• 长寿命：仅消耗电能，不消耗工质，寿命受限于轴承和电子学。

3.3 动量轮数学模型：从电气方程到控制模型

建立准确的数学模型是进行动量轮控制器设计、故障诊断和性能仿真的基础。动量轮的数学模型包括电气子模型、机械子模型和热模型三个层面。

3.3.1 电气方程

对于无刷直流电机，在dq旋转坐标系下的电压方程为：

U_d = R·I_d + L·dI_d/dt - ω_e·L·I_q
U_q = R·I_q + L·dI_q/dt + ω_e·L·I_d + K_e·Ω

其中 U_d、U_q 为d/q轴电压，I_d、I_q 为d/q轴电流，R 为定子电阻，L 为定子电感，ω_e 为电角速度（ω_e = p·Ω，p 为极对数），K_e 为反电动势常数。在磁场定向控制（FOC）中，通常令 I_d = 0，则电磁力矩简化为：

T_e = 1.5 · p · K_e · I_q = K_t · I_q

3.3.2 机械方程

飞轮转子的力矩平衡方程为：

J_w · dΩ/dt = T_e - B·Ω - T_f - T_load

其中 B 为粘性阻尼系数，T_f 为库仑摩擦力矩（与转速方向相反），T_load 为外部负载力矩（通常为0，因为飞轮在真空中运行）。在稳态时（dΩ/dt = 0），电机力矩等于摩擦力矩与阻尼力矩之和。

3.3.3 完整控制模型

将电气方程和机械方程联立，并考虑转速反馈，可得到动量轮的传递函数模型。在工程中，通常将动量轮简化为一个一阶惯性环节加积分环节：

G(s) = K_t / (J_w · s + B) · 1/s

其中 s 为拉普拉斯算子。该模型用于姿控系统设计时，可近似描述动量轮从力矩指令到转速响应的动态特性。实际动量轮的带宽通常在 10~50 Hz 之间，响应延迟约 5~20 ms。

3.4 案例拆解：某遥感卫星动量轮转速异常波动分析

某低轨遥感卫星在轨运行第 14 个月时，姿控分系统遥测数据显示，动量轮转速在指令恒定的情况下出现 ±30 rpm 的周期性波动，波动周期约为 90 分钟（一个轨道周期）。同时，卫星姿态角出现 0.002° 的微小振荡，虽未超出任务指标，但已引起地面测控团队的关注。

3.4.1 初始现象与数据

• 动量轮转速指令：4000 rpm（恒定）
• 实际转速遥测：3970~4030 rpm 波动
• 电机电流：0.35 A ± 0.05 A 波动
• 轴承温度：22°C ± 1.5°C 波动
• 卫星姿态：滚动轴 0.0015° 振荡，俯仰轴 0.002° 振荡

3.4.2 分析判断过程

地面团队首先排除了控制指令异常的可能性，因为姿控计算机输出的指令码流校验无误。随后，利用动量轮数学模型进行仿真复现：

仿真条件：J_w = 0.08 kg·m²，K_t = 0.05 N·m/A，B = 2e-4 N·m·s/rad，T_f = 5e-3 N·m
输入：恒定电流 0.35 A
输出：转速波动 ±25 rpm（与遥测接近）

仿真结果表明，仅凭恒定的电机电流无法产生如此大的转速波动，说明存在外部干扰力矩。进一步分析轨道环境，发现卫星在光照区和阴影区交替时，太阳帆板产生的热弹性变形导致卫星质心发生微小偏移，从而产生周期性干扰力矩。该干扰力矩通过动量交换传递到飞轮，引起转速波动。

3.4.3 实际操作步骤

1. 数据比对：将动量轮转速波动曲线与轨道光照周期进行相关性分析，确认波动周期与轨道周期一致。
2. 模型修正：在动量轮数学模型中增加周期性干扰力矩项 T_dist = A·sin(2πt/T_orbit)，其中 A = 0.002 N·m，T_orbit = 5400 s。
3. 仿真验证：修正后模型输出的转速波动幅值与遥测数据误差小于 5%。
4. 控制策略调整：在姿控软件中增加一个陷波滤波器，频率设置为轨道频率（约 0.185 mHz），衰减系数 0.7，以抑制该频段的干扰传递。
5. 在轨验证：上注新控制参数后，动量轮转速波动降至 ±5 rpm 以内，姿态振荡减小至 0.0005°。

3.4.4 结果与复盘

通过本次排障，团队认识到动量轮转速波动不仅是动量轮自身的问题，还可能是卫星平台其他子系统（如太阳帆板、热控）的干扰耦合。复盘经验如下：

• 动量轮数学模型必须包含干扰力矩项，才能准确反映在轨行为。
• 转速波动分析应结合轨道环境数据（光照、地磁、大气阻力）进行联合诊断。
• 陷波滤波器是抑制特定频率干扰的有效手段，但需注意不要引入相位滞后影响姿控稳定性。

3.5 常见误区与设计要点

在动量轮原理理解和工程应用中，工程师常陷入以下误区：

• 误区一：将动量轮视为理想力矩源。实际动量轮存在摩擦力矩、电流环带宽限制和转速测量噪声，其输出力矩并非理想的阶跃响应。设计姿控律时，应使用包含电气和机械动态的完整模型。
• 误区二：忽略动量轮的饱和特性。动量轮只能在有限转速范围内工作。当飞轮接近饱和转速时，其输出力矩能力急剧下降，甚至无法提供所需力矩。姿控系统设计必须包含卸载逻辑。
• 误区三：混淆角动量与力矩的概念。角动量是动量轮储存的能量状态，力矩是角动量的变化率。卸载的本质是改变飞轮的角动量，而非直接产生控制力矩。
• 误区四：忽视温度对参数的影响。轴承摩擦力矩 T_f 随温度变化显著，低温时润滑剂黏度增大，T_f 可增大数倍。电机电阻 R 也随温度变化，影响电流环增益。高精度应用需对动量轮进行温度补偿。