使用OMPython的ModelicaSystem时线性系统求解器失败问题排查

OMCSessionZMQ与ModelicaSystem的核心差异在于模型加载、参数配置及模拟选项的默认行为不同，具体可能的原因：

1. 参数设置格式错误：setParameters传入的params.xxx前缀可能无效，导致参数未被正确修改，模型状态与OMCSessionZMQ的成功场景不一致，引发线性系统求解失败。
2. 模拟选项默认值差异：ModelicaSystem的默认线性求解器配置、初始化策略可能与OMCSessionZMQ不同，未启用鲁棒性更强的 fallback 求解器逻辑。
3. 编译流程差异：ModelicaSystem可能采用了不同的编译优化或初始化步骤，导致模型的线性系统结构更病态。

1. 修正参数设置格式 #

ModelicaSystem的setParameters方法不需要params.前缀，直接使用顶层模型的参数名即可。修改参数设置代码：

param_changes = []
for param_name, param_value in parameters.items():
    param_changes.append(f'{param_name}={param_value}')
mod.setParameters(param_changes)

2. 强制设置鲁棒性线性求解器 #

显式指定更稳定的线性求解器，并启用全主元策略，覆盖ModelicaSystem的默认配置：

mod.setSimulationOptions([
    "startTime=0.0",
    f"stopTime={simulation_time}.0",
    "tolerance=1e-06",
    "stepSize=0.002",
    "linearSolver='ma57'",  # 选择鲁棒性更强的稀疏线性求解器
    "useTotalPivoting=true",  # 启用全主元消去，应对病态矩阵
    "lv=LOG_LS"  # 可选：启用线性求解器日志，排查具体失败原因
])
mod.simulate()

3. 对齐OMCSessionZMQ的模拟选项 #

先通过OMCSessionZMQ获取默认模拟选项，再同步到ModelicaSystem：

# 用OMCSessionZMQ获取默认选项
omc = omp.OMCSessionZMQ()
default_options = omc.sendExpression("getSimulationOptions()")
# 将选项转换为ModelicaSystem可接受的字符串格式
option_strings = [f"{key}={repr(value)}" for key, value in default_options.items()]
# 在ModelicaSystem中设置这些选项
mod.setSimulationOptions(option_strings)

4. 直接复用OMCSessionZMQ的模拟逻辑 #

如果上述方法无效，可以在ModelicaSystem中直接调用OMC的simulate命令，完全对齐成功场景：

# 先设置参数（修正格式后）
mod.setParameters(param_changes)
# 直接调用OMC模拟命令
result = mod.sendExpression(f"simulate(System, stopTime={simulation_time}.0)")
# 从结果中获取输出文件路径，再读取结果
import pandas as pd
res_df = pd.read_csv(result['resultFile'])

• 启用lv=LOG_LS日志后，查看线性系统961的具体特征（比如矩阵秩、条件数），定位模型中导致病态的组件。
• 先尝试不修改参数，直接用ModelicaSystem模拟，如果成功，说明问题出在参数修改的格式或参数值本身；如果失败，说明是默认配置差异导致。

内容的提问来源于stack exchange，提问作者Andrea