2. 模型降阶的基本概念与数学基础

在数字孪生系统中，物理对象往往由有限元、计算流体力学或多体动力学模型描述，其离散自由度可达百万甚至上亿级。若把这样的全阶模型（Full-Order Model, FOM）直接嵌入实时回路，单步求解时间常常远超孪生体所需的毫秒级刷新周期。模型降阶（Model Order Reduction, MOR）要解决的核心问题，就是在保留系统关键输入输出行为的前提下，把高维状态空间压缩到低维子空间，使在线计算量下降几个数量级。

降阶的必要性：从算得准到算得快

以某热-结构耦合数字孪生为例：全阶有限元模型含 120 万个节点自由度，隐式时间积分单步约 8 s，而孪生体要求 20 ms 内完成一次状态更新，两者相差约 400 倍。降阶的目标不是让模型“更简单”，而是让模型在给定精度容差内满足实时性约束。这里存在一个基本权衡：降阶阶数 r 越小，在线求解越快，但逼近误差越大；r 越大，精度越高，但可能重新逼近实时上限。因此降阶必须与误差度量、稳定性判据一起设计，而不是单纯追求低维。

全阶模型与降阶模型的关系

设全阶模型为线性时不变系统：

E ẋ(t) = A x(t) + B u(t)
y(t)   = C x(t) + D u(t)

其中 x ∈ ℝⁿ 为全阶状态，u 为输入，y 为输出，E 可能奇异（描述代数约束）。降阶模型通过投影构造：选取左、右投影基 V, W ∈ ℝⁿˣʳ，令 x ≈ V xᵣ，并施加 Petrov-Galerkin 条件 Wᵀ(E V ẋᵣ − A V xᵣ − B u) = 0，得到：

Eᵣ ẋᵣ(t) = Aᵣ xᵣ(t) + Bᵣ u(t)
yᵣ(t)    = Cᵣ xᵣ(t) + D u(t)

Eᵣ = WᵀEV,  Aᵣ = WᵀAV,  Bᵣ = WᵀB,  Cᵣ = CV

当 W = V 时称为 Galerkin 投影，适合结构动力学等对称系统；当 W ≠ V 时称为 Petrov-Galerkin 投影，常用于非对称或含代数约束的系统。降阶模型不是全阶模型的“近似副本”，而是全阶模型在低维子空间上的受限映射：它继承了输入输出通道 B、C 的物理含义，但内部状态 xᵣ 不再对应具体节点位移，而是模态或子空间坐标。

误差度量：降阶是否“够用”的判据

降阶误差不能只看状态向量，因为 xᵣ 与 x 的物理量纲不同。工程上更关心输出误差和频域误差。常用度量包括：

• 时域相对输出误差：对给定输入 u(t)，计算 ‖y(t) − yᵣ(t)‖₂ / ‖y(t)‖₂ 的最大值或均方根值，要求小于容差 ε（如 1%）。
• H∞ 误差：‖G(s) − Gᵣ(s)‖∞ = sup_ω σ̄(G(jω) − Gᵣ(jω))，衡量最坏频率点的增益误差，适合鲁棒控制场景。
• H₂ 误差：‖G − Gᵣ‖₂，反映能量意义下的平均误差，与脉冲响应误差直接相关。
• 特征值/模态误差：对结构动力学，关注前 r 阶固有频率与阻尼比的相对偏差，通常要求频率误差小于 0.5%。

以热-结构孪生为例，若取 r = 20，时域输出误差为 0.8%，H∞ 误差为 −42 dB，单步耗时 15 ms，则满足 20 ms 实时约束且精度达标。若把 r 降到 10，单步耗时降至 6 ms，但输出误差升至 4.5%，超出容差，说明降阶阶数不能仅由速度决定。

稳定性与保结构降阶的基本要求

降阶模型必须继承全阶模型的关键结构性质，否则在线仿真可能发散或产生非物理结果。基本要求包括：

• 稳定性保持：若全阶系统渐近稳定（A 的特征值均在左半平面），则降阶系统 Aᵣ 也应稳定。Galerkin 投影对对称正定系统可保持稳定性；对一般系统，需采用稳定保持算法（如 IRKA 的稳定化变体）。
• 无源性/耗散性保持：电路、热网络等系统要求降阶后仍无源，否则与外部电路耦合时可能产生能量不守恒。可用 PRIMA、无源保持投影等方法。
• 结构对称性保持：结构动力学中 M、K、C 矩阵对称，降阶后 Mᵣ = VᵀMV、Kᵣ = VᵀKV 仍对称，保证模态实特征值。
• 代数约束保持：若 E 奇异，降阶后 Eᵣ 也应保持相应奇异结构，避免把微分代数系统误化为纯微分系统。
• 输入输出通道一致：Bᵣ、Cᵣ 必须由同一投影基构造，保证降阶模型与原系统在相同输入下输出可比。

一个可操作的降阶检查清单

在实际数字孪生项目中，可按以下顺序判断降阶是否可用：

只有上述条件同时满足，降阶模型才能作为数字孪生在线计算内核。若某一项不达标，应优先调整投影基构造方法或阶数 r，而不是直接降低精度要求。降阶的本质是在精度、速度与结构性质之间寻找可验证的平衡点，而不是单纯压缩维度。