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 为实际力/力矩。这样既便于替换不同执行机构模型,也便于在仿真中注入故障(如推力器失效、轮速卡死)。

执行机构仿真接口与信号流 控制器 控制律 / 分配算法 指令 化学推进 延迟/上升/衰减 · 最小开机 电推进 功率-推力映射 · 慢建立 动量轮 角动量守恒 · 轮速饱和 磁力矩器 T = M × B · 最大磁矩 航天器 六自由度 动力学模型 轨道 + 姿态 质量/惯量变化 环境力矩 状态反馈:轮速 / 推进剂余量 / 饱和标志 / 实际力矩 约束注入点:最小开机时间 · 最大轮速 · 最大磁矩 · 功率上限 · 推进剂余量 · 推力方向安装矩阵 地磁场模型 B 轨道位置 → 磁场矢量

七、建模中的常见问题与处理建议

第一,忽略执行机构动态会导致仿真结果过于乐观。例如,化学推进若按理想脉冲建模,小脉冲控制的速度增量误差可能达到 10% 以上。建议在仿真中至少保留一阶延迟和上升时间。

第二,动量轮饱和处理需要与控制器协同。仿真中若只建模轮速上限而不输出饱和标志,控制器可能持续发出无法执行的力矩指令,导致姿态误差累积。正确做法是:轮速达到上限时,模型输出饱和标志,控制器切换至卸载模式或重新分配力矩。

第三,磁力矩器必须与地磁场模型同步。由于 T = M × B,磁场方向随轨道位置变化,仿真中若使用固定磁场矢量,磁力矩器产生的力矩方向将完全错误。建议在每个仿真步长中根据当前轨道位置更新磁场矢量。

第四,电推进的慢建立特性对长时间轨道机动影响显著。若仿真步长大于推力建立时间,可能完全丢失建立过程。建议电推进仿真步长不大于建立时间的 1/5,或在模型内部进行子步积分。

在工程实践中,建议为每类执行机构建立“标称模型 + 偏差模型”两级结构。标称模型用于控制律设计和快速验证,偏差模型用于蒙特卡洛打靶,评估参数不确定性(如推力偏差 ±5%、比冲偏差 ±3%、轮速摩擦 ±10%)对控制精度的影响。

综上,执行机构建模的关键在于:明确输入输出接口、显式表达物理约束、保留影响控制精度的主要动态环节。只有将执行机构模型与控制器、航天器动力学放在同一仿真回路中闭环验证,才能真实反映轨道控制系统的实际性能。