3. 航空航天对象特性分析

半实物仿真的可信度,最终取决于被仿对象模型是否忠实反映了真实飞行器的动态特性。在把飞控计算机、舵机、惯组接入仿真回路之前,必须先回答一个基础问题:飞行器本体、气动、动力、执行与测量这五类对象,各自以什么物理机理工作,又该以什么形式进入实时模型。本节按“机理—建模需求—实时约束”三条线索逐一展开,并在最后给出一个可复用的对象建模清单。

一、六自由度运动:刚体动力学的核心状态

飞行器六自由度运动指质心三轴平动(沿机体轴 x、y、z 的速度 u、v、w)与绕质心三轴转动(滚转 p、俯仰 q、偏航 r)。建模时通常采用机体坐标系:x 轴指向机头,y 轴指向右翼,z 轴指向机腹,角速度与力矩均按右手定则定义。这一约定必须在全课程、全模型、全接口中保持一致,否则气动力矩符号会整体反号。

动力学方程由牛顿—欧拉方程给出,平动部分为 m(V̇ + ω×V) = F,转动部分为 I ω̇ + ω×(Iω) = M。其中 F、M 是气动力/力矩、发动机推力/推力偏心矩、重力在机体系投影的合力。姿态更新常用四元数积分,避免欧拉角在俯仰 ±90° 附近的奇异:q̇ = 0.5 · q ⊗ [0, p, q, r],再归一化。位置更新则需把机体系速度经方向余弦矩阵转到地面系积分。

实时建模的关键需求有三点:一是积分步长与稳定性,显式欧拉在 1 kHz 下对小型飞行器通常够用,但大惯量比对象需改用四阶龙格库塔或半隐式方法;二是状态量必须可被外部(飞控、视景、测试设备)按固定周期读取,不能出现半更新状态;三是重力、地球自转、大气密度等环境项要作为独立可配置模块,便于在不同飞行包线下替换。

二、气动面:从迎角侧滑角到气动系数

气动面包括升降舵、副翼、方向舵,部分飞行器还有襟翼、扰流板、矢量喷管。其工作机理是:舵面偏转改变局部气流方向,从而改变升力与力矩。工程上不直接算流场,而是用气动系数表建模:给定马赫数 Ma、迎角 α、侧滑角 β 和舵偏 δ,查表得到升力系数 CL、阻力系数 CD、俯仰力矩系数 Cm 等,再乘以动压 q̄ = 0.5 ρ V² 和参考面积/弦长得到力与力矩。

例如某样例无人机在 Ma=0.3、α=6°、升降舵 δe=-5° 时,查表得 CL=0.62、Cm=-0.08,取 ρ=1.225 kg/m³、V=60 m/s、S=0.8 m²、c̄=0.35 m,则 q̄≈2205 Pa,升力 L≈1094 N,俯仰力矩 M≈-61.7 N·m。这些数值属于教学示例,用于说明查表—乘动压—得力的完整链路。

建模需求上,气动数据通常以多维插值表形式存在,实时模型必须支持线性插值与边界外推保护;舵面偏转速率限制、舵面失效(卡死、反向)也应作为可注入故障,供半实物仿真验证飞控重构能力。

三、发动机:推力与动态响应

发动机建模分两层:稳态推力特性和动态响应。稳态推力一般表示为油门 δt、马赫数、高度 H 的函数 T = f(δt, Ma, H),以二维或三维表存储。动态响应则用一阶惯性环节近似:Ṫ = (Tcmd - T)/τ,涡喷发动机 τ 约 0.5~2 s,电动机 τ 可小到 0.05 s。若考虑螺旋桨,还需把推力与转速、前进比关联。

实时建模中,发动机模型要输出推力大小、推力作用点与推力方向(矢量喷管时含偏转角),并把这些量送入六自由度方程。推力偏心会直接产生俯仰或偏航力矩,是半实物仿真中容易被忽略却影响配平的重要项。

四、执行机构与传感器:闭环回路的物理接口

执行机构主要是舵机与油门伺服。舵机模型需包含:指令到实际偏转的带宽(常用二阶或一阶加延迟)、速率限幅、位置限幅、死区与回中特性。半实物仿真中,真实舵机可能直接接入回路,此时模型只保留限幅与故障注入,其余由实物体现。

传感器包括惯性测量单元(IMU)、大气数据系统、磁力计、GPS。IMU 建模要点是:输出比力与角速度,叠加零偏、随机游走、标度因数误差和量化噪声;大气数据需由位置与速度反算静压、动压、迎角侧滑角,再叠加噪声。传感器模型的采样率、延迟和噪声谱必须与真实器件一致,否则飞控滤波器在半实物回路中会表现出与真实飞行不同的相位裕度。

飞行器对象建模信号链与实时接口关系 执行机构 舵机 / 油门伺服 带宽·限幅·速率 气动面 α β Ma δ → C_L C_D C_m q̄ = 0.5ρV² 发动机 T = f(δ_t, Ma, H) 一阶惯性 τ 六自由度刚体动力学 m(V̇ + ω×V) = F I ω̇ + ω×(Iω) = M 状态:u v w p q r · 四元数积分 传感器 IMU 比力/角速度 大气数据 α β V 零偏·噪声·延迟 飞控计算机 控制律 / 滤波器 半实物接入点 实时约束:定步长积分 · 状态原子更新 · 接口周期一致 · 故障可注入 气动/发动机查表需边界保护;传感器噪声谱与真实器件一致;执行机构限幅与实物一致

五、建模需求汇总与可复用清单

把上述五类对象放在同一条信号链上看:飞控输出舵偏与油门指令,执行机构把它变成实际偏转与推力,气动面与发动机据此产生力与力矩,六自由度方程积分出运动状态,传感器再把状态变成带噪声的测量值送回飞控。半实物仿真的价值,就在于把其中一段或多段换成真实硬件,检验闭环在真实接口特性下是否仍然稳定。

下表给出一份可直接用于对象建模评审的清单,字段均为本节专业语境所需,不含无关管理项。

对象 关键输入 关键输出 建模形式 实时约束
六自由度运动 F、M(机体系) u v w p q r、位置、姿态 牛顿—欧拉 + 四元数积分 定步长 ≥1 kHz,状态原子更新
气动面 α β Ma δ C_L C_D C_m 及力/力矩 多维插值表 + 动压计算 插值边界保护,舵面故障可注入
发动机 δ_t Ma H 推力 T、推力方向 稳态表 + 一阶惯性 τ 与真实发动机一致
执行机构 舵偏/油门指令 实际偏转、实际推力 二阶/一阶 + 延迟 + 限幅 与实物带宽、速率一致
传感器 真实运动状态 比力、角速度、气压、α β 理想量 + 零偏 + 噪声 + 延迟 采样率与噪声谱匹配器件
评审时优先核对三件事:坐标系与符号约定是否全文一致;气动/发动机查表在包线边界是否做了外推保护;传感器延迟是否与真实器件标定值一致。这三项出错,半实物仿真结果往往表现为“飞控在地面正常、接入实物后发散”,而问题根源并不在飞控本身。