提升 AMESim 仿真速度——刚性问题 原创

刚性（stiffness）是指同一仿真模型中同时存在时间尺度相差悬殊的动态。液压系统里，小容腔中微秒级的压力波（|λ|≈1e5）和秒级的活塞运动（|λ|≈1）共存，刚性比 S = |λmax|/|λmin| 可达 1e6。刚性是系统的物理属性，不是建模错误；它拖慢仿真的机制，是让积分步长被最快的模态绑架，导致总步数成数量级增加。

什么是刚性

模型编译后得到一组微分方程，每个状态对应一个模态，时间常数 τ = 1/|λ|。刚性系统的特征是快慢模态共存且跨度极大：

物理域	快模态来源	慢模态来源
液压	小容腔 + 高体积模量（油液 B≈1.7e10 Pa）挤出高频压力波	活塞/缸体大质量运动
热	小热容节点	大热容体温度演化
机械	高刚度弹簧短轴	大惯量回转
电气	低电感支路	大时间常数 RC

判断刚性的量化指标是刚性比 S = |λmax|/|λmin|。S 在 1e3 以上开始影响效率，1e6 量级是典型刚性系统。

刚性如何影响仿真速度

积分步长同时受两道闸约束：

• 精度闸由截断误差决定：h_精度 ∝ 容差^(1/(p+1))，容差越严步长越小。
• 稳定性闸由误差传播决定：对显式法（Adams/RK），h_稳定性 ≤ 2/|λmax|。

刚性系统慢的根因，是稳定性闸把步长压到了与精度无关的量级。液压里 |λmax|≈1e5，h_稳定性 ≤ 2e-5 s；容腔再小一点直接 2e-8。显式法被迫用微秒级步长去算秒级仿真，刚性比 1e6 时步数多 1e6 倍——不崩也算到天荒地老，多数情况求解器在步长触底时报 step size unacceptably small 中止。

BDF/DASSL 属于隐式法。向后欧拉作用在测试方程 y'=λy 上，误差传播因子为 1/(1−hλ)；对衰减模态（Re(λ)<0）无论 h 多大恒小于 1，这叫 A-稳定。所以 BDF 可以用大步长跑，只追踪慢动态包络，快模态一带而过。代价是每步要做 Newton 迭代、单步更贵，但刚性系统上步数大减的收益远超单步成本。

如何判断刚性

1. Run statistics 中 log10(min step) < -8（尤其 -10），是刚性最直观的信号。
2. Frequencies 面板点 Start 做线性化，出现红/橙低阻尼（damping ratio < 0.01）的高频模态，即刚性来源。
3. BDF steps by order 若长期停留在 1 阶，说明系统极刚。

最可靠的对照法：看 Frequencies 面板执行线性化，确认存在红/橙低阻尼高频模态，且 min step 极低，即可坐实为刚性。不可通过"手动切 Adams"做对照——AMESim 标准积分器不支持手动切换求解器。

需排除另外两类常伴生问题：代数环（看编译的 Generated implicits、State contributions 中的蓝色行）和不连续（看 Total discontinuities、Discontinuities 面板占比）。

如何优化

优化有两条方法，优先级严格：

方法 A：求解器配置（首选，零物理代价）

• 求解器交给标准积分器自动选，不需手动指定：含代数环（隐式变量）走 DASSL，纯 ODE 走 LSODA（内部自动在 Adams 与 BDF 间切换）
• 启用稳态初始化，从平衡点起步
• 最大步长 ≈ 0.1 × 快模态时间常数
• 分级容差：RelTol 1e-3～1e-4；AbsTol 按量纲给（压力 1e2～1e3 Pa、位移 1e-6～1e-8）

方法 B：物理降刚（次选，有代价）

• 体积模量从 1.7e10 调到 1.3e9（约降 13 倍），或设 1～5% 含气率
• 合并 / 增大相邻小容腔
• 步长约放大 3～4 倍，但会改变波速与瞬态响应

两条红线：① 放宽容差对刚性无效——精度闸管不到稳定性闸，必须先确认标准积分器自动走 DASSL（含代数环）或 LSODA（纯 ODE），容差才重新生效；② 给零质量阀芯加质量是解代数环，不是降刚性。

例程

某液压回路仿真在 0.31 s 卡死，日志反复 step size unacceptably small，Performance Analyzer 显示 log10(min step) = -10.4、BDF 1 阶步占 94%。

处理：求解器已是 BDF，启用稳态初始化，最大步长设为 0.1 × 快模态时间常数，压力 AbsTol 提至 1e3，物理参数未改动。仿真由 2 小时缩至 18 分钟，峰值压力与切换时刻均一致。

方法 B 未启用——该回路分析目标含瞬态冲击，不宜为提速牺牲响应精度。若方法 A 仍不满足速度要求，且物理上允许近似，再考虑方法 B。

三个易混淆点

• 提速 ≠ 降刚。多数情况方法 A 足够，不必动物理参数。物理降刚会改变动态响应，做错会得到"算得快但错"的结果。
• 刚性 ≠ 代数环。刚性是数值尺度问题（特征值跨度大），代数环是拓扑缺陷（因果链缺储能件）。小容腔既招高频模态又易触发代数环，但解法完全不同。
• 加质量 ≠ 降刚。给零质量阀芯加质量是在解代数环（引入惯性积分缓冲），不是降刚性。

#Amesim #液压仿真 #CAE工程师 #仿真提速