3. 轨道控制执行机构建模
轨道控制系统的仿真精度,很大程度上取决于执行机构模型是否真实反映了其动力学特性、物理约束以及与控制器之间的信号接口。本节围绕化学推进、电推进、动量轮和磁力矩器四类典型执行机构,说明它们在仿真中应如何建模、需要暴露哪些接口、以及哪些约束条件必须显式表达。
一、执行机构在仿真中的角色定位
在轨道控制仿真回路中,执行机构位于控制器与航天器动力学之间,承担“控制指令→物理作用力/力矩”的转换。仿真中执行机构模型需要完成三件事:接收控制指令(推力大小、力矩指令、转速指令等),根据自身动力学特性计算实际输出,并将输出以力/力矩或角动量的形式施加到航天器六自由度模型上。同时,执行机构模型还需向控制器反馈状态量(如剩余推进剂、轮速、磁力矩器电流等),供控制律进行饱和判断和分配。
因此,执行机构建模的核心不是简单乘以一个增益,而是要在指令与输出之间建立包含动态响应、饱和约束、效率损失和状态反馈的完整映射。
二、化学推进建模
化学推进(如单组元肼推力器、双组元发动机)的仿真模型通常采用“指令—延迟—上升—稳态—关机”的时序描述。关键参数包括:额定推力 F、最小开机时间 tmin、关机后残余推力衰减时间、推力上升/下降斜率、比冲 Isp 以及推力方向安装矩阵。
仿真中,控制器给出开机指令后,推力并非立即达到额定值。一个常用的简化模型为:
F(t) = 0, t < t_delay
F(t) = F_rated * (t - t_delay)/t_rise, t_delay ≤ t < t_delay + t_rise
F(t) = F_rated, t_delay + t_rise ≤ t < t_off
F(t) = F_rated * (1 - (t - t_off)/t_fall), t_off ≤ t < t_off + t_fall
F(t) = 0, t ≥ t_off + t_fall
其中 tdelay 为指令到推力建立的延迟,trise 为上升时间,tfall 为关机衰减时间。这些参数直接影响小脉冲控制的精度,尤其在脉冲宽度接近最小开机时间时,实际产生的速度增量会明显偏离理想值。
约束方面,化学推进模型必须包含:最小开机时间(低于该值推力器可能无法正常点火或重复性极差)、最大占空比(热约束)、推进剂剩余量(积分消耗,dm/dt = −F/(Ispg0))、以及推力器安装方向与质心位置决定的力矩臂。仿真接口通常以“开机指令+期望脉宽”为输入,以“实际推力矢量+推进剂消耗率+累计消耗量”为输出。
三、电推进建模
电推进(如霍尔推力器、离子推力器)与化学推进的最大区别在于:推力小但比冲高,且推力建立过程更复杂。仿真中需要关注:推力与功率的耦合关系、束流方向、羽流发散角、以及推力器对供电功率和热状态的依赖。
一个典型的电推进简化模型为:
T = η * 2 * P / (Isp * g0)
其中 η 为总效率,P 为输入电功率,Isp 为比冲,g0 为标准重力加速度。
实际仿真中,推力并非仅由功率决定,还受到阳极流量、磁场强度、阴极状态等影响。因此更合理的做法是建立“指令推力→功率需求→实际推力”的映射,并加入推力建立时间(通常数十秒到数分钟)和最小稳定工作功率约束。
电推进的仿真接口通常包括:推力指令、功率上限、最小/最大推力、推力方向安装矩阵、以及累计工作时间。输出为实际推力矢量、功率消耗、推进剂消耗率(氙气或氪气)。由于电推进推力小,仿真步长需要足够小以捕捉长时间连续工作的累积效果,同时避免因步长过大导致推力积分误差。
四、动量轮建模
动量轮通过改变自身角动量来对航天器产生反作用力矩。仿真中,动量轮模型的核心是角动量守恒关系:
T_wheel = -dH_wheel/dt
H_wheel = J_wheel * ω_wheel
其中 Jwheel 为轮体转动惯量,ωwheel 为轮速。控制器给出力矩指令后,动量轮通过电机加速或减速,产生反作用力矩作用于星体。仿真中需要建模的约束包括:最大轮速(饱和)、最大力矩(电机能力)、摩擦力矩(低速时非线性明显)、以及轮速过零时的死区效应。
动量轮仿真接口通常以力矩指令为输入,以实际输出力矩、当前轮速、饱和标志为输出。当轮速接近上限时,模型应输出饱和标志,控制器据此进行力矩重新分配或启动卸载。此外,多轮配置(如三正交+一斜装)需要提供安装矩阵,将轮系力矩映射到星体坐标系。
五、磁力矩器建模
磁力矩器通过与地磁场相互作用产生力矩:
T_mag = M × B
其中 M 为磁矩矢量,B 为当地地磁场矢量。仿真中,磁力矩器模型需要输入当地磁场强度(由轨道位置和地磁模型计算),输出磁矩指令对应的实际力矩。约束包括:最大磁矩、电流响应时间、以及磁力矩器只能产生垂直于磁场方向的分量这一固有特性。
磁力矩器的仿真接口通常包括:磁矩指令(三轴)、当地磁场矢量、最大磁矩限制。输出为实际磁矩和产生的力矩。由于磁力矩器效率低、力矩小,通常用于动量轮卸载或低精度姿态控制,仿真中需要与轨道位置和地磁模型同步更新。
六、四类执行机构对比与仿真接口汇总
下表汇总了四类执行机构在仿真建模中的关键差异:
| 执行机构 | 输入指令 | 核心动力学 | 主要约束 | 输出反馈 |
|---|---|---|---|---|
| 化学推进 | 开机/脉宽 | 延迟+上升+稳态+衰减 | 最小开机时间、推进剂余量 | 实际推力、消耗量 |
| 电推进 | 推力/功率 | 功率-推力映射、慢建立 | 最小/最大功率、建立时间 | 实际推力、功率、消耗率 |
| 动量轮 | 力矩指令 | 角动量守恒、电机动态 | 最大轮速、最大力矩 | 实际力矩、轮速、饱和标志 |
| 磁力矩器 | 磁矩指令 | T = M × B | 最大磁矩、只能垂直分量 | 实际磁矩、力矩 |
在仿真实现中,建议将每类执行机构封装为独立模块,统一接口形式为:update(dt, command, state) → (output, state_new)。其中 command 为控制器输出,state 包含内部状态(如轮速、推进剂余量、温度等),output 为实际力/力矩。这样既便于替换不同执行机构模型,也便于在仿真中注入故障(如推力器失效、轮速卡死)。
七、建模中的常见问题与处理建议
第一,忽略执行机构动态会导致仿真结果过于乐观。例如,化学推进若按理想脉冲建模,小脉冲控制的速度增量误差可能达到 10% 以上。建议在仿真中至少保留一阶延迟和上升时间。
第二,动量轮饱和处理需要与控制器协同。仿真中若只建模轮速上限而不输出饱和标志,控制器可能持续发出无法执行的力矩指令,导致姿态误差累积。正确做法是:轮速达到上限时,模型输出饱和标志,控制器切换至卸载模式或重新分配力矩。
第三,磁力矩器必须与地磁场模型同步。由于 T = M × B,磁场方向随轨道位置变化,仿真中若使用固定磁场矢量,磁力矩器产生的力矩方向将完全错误。建议在每个仿真步长中根据当前轨道位置更新磁场矢量。
第四,电推进的慢建立特性对长时间轨道机动影响显著。若仿真步长大于推力建立时间,可能完全丢失建立过程。建议电推进仿真步长不大于建立时间的 1/5,或在模型内部进行子步积分。
综上,执行机构建模的关键在于:明确输入输出接口、显式表达物理约束、保留影响控制精度的主要动态环节。只有将执行机构模型与控制器、航天器动力学放在同一仿真回路中闭环验证,才能真实反映轨道控制系统的实际性能。