OpenModelica电气组件因果适配器实现与多FMU仿真故障问询

一、因果适配器的实现与工作原理 #

1. OpenModelica中因果适配器的实现代码 #

你可以直接基于Modelica标准库的逻辑编写适配器，以下是通用的电压转电流、电流转电压适配器代码：

package ElectricalAdapters
  model GeneralVoltageToCurrentAdaptor "无因果电压转因果电流适配器"
    Modelica.Electrical.Analog.Interfaces.PositivePin p;
    Modelica.Electrical.Analog.Interfaces.NegativePin n;
    output Modelica.SIunits.Current I;
    Modelica.SIunits.Voltage V = p.v - n.v;
  equation
    0 = p.i + n.i;
    I = p.i;
  end GeneralVoltageToCurrentAdaptor;

  model GeneralCurrentToVoltageAdaptor "无因果电流转因果电压适配器"
    Modelica.Electrical.Analog.Interfaces.PositivePin p;
    Modelica.Electrical.Analog.Interfaces.NegativePin n;
    input Modelica.SIunits.Current I;
    Modelica.SIunits.Voltage V = p.v - n.v;
  equation
    p.i = I;
    n.i = -I;
  end GeneralCurrentToVoltageAdaptor;
end ElectricalAdapters;

将上述代码保存为.mo文件后，即可在OpenModelica中直接调用。

2. 因果适配器的工作原理 #

• 为何提取电压(V)和电流(I)：Modelica电气连接器是**无因果（acausal）的，端口的电压、电流没有固定的输入输出关系，完全通过联立方程求解；而FMU的接口必须是因果（causal）**的，必须明确区分输入、输出变量。提取V和I是为了将无因果的电气端口转换为符合FMU要求的因果接口。
• 无因果到因果的转换逻辑：
  • 电压转电流适配器：通过测量无因果端口的电压差（V=p.v-n.v），同时将端口电流作为输出（I=p.i），并满足基尔霍夫电流定律（p.i + n.i = 0），最终实现“输入电压、输出电流”的因果接口。
  • 电流转电压适配器：接收外部输入的电流（I），将其施加到无因果端口（p.i=I，n.i=-I），同时输出端口的电压差（V=p.v-n.v），实现“输入电流、输出电压”的因果转换。

二、多FMU电气电路仿真问题 #

1. 是否可用多个FMU表示电气电路中的各组件？ #

完全可以。通过因果适配器正确转换接口后，多个FMU可以组合成完整的电气电路，这正是FMU模块化复用的核心价值之一。

2. 多FMU连接仿真失败的原因？ #

报错涉及kinsol求解器错误码（-5、-7、-11），核心原因通常是：

• 代数环问题：多个FMU连接后形成了无法直接求解的代数环。每个FMU的因果接口与外部连接形成循环依赖，求解器找不到稳定的迭代收敛路径。
• 求解器参数不匹配：默认的求解器配置（如密度阈值、迭代次数）无法处理多FMU交互带来的刚性或病态方程组。
• 接口变量连续性问题：FMU间的输入输出变量在仿真初始时刻或迭代过程中出现不连续，导致求解器无法收敛。

3. 如何解决该问题？ #

• 调整求解器参数：你使用的-nlssMaxDensity=0.05是通过降低非线性方程组的密度阈值，让求解器采用更适合稀疏矩阵的求解策略，从而缓解代数环或病态方程组的收敛问题。还可以尝试：
  • 增加迭代次数：-nlssMaxIterations=1000
  • 调整收敛容差：-nlssTol=1e-6
• 优化FMU接口设计：合理分配各FMU的输入输出逻辑，减少不必要的因果转换。比如串联电路中，电压源FMU输出电压，电阻FMU输入电压、输出电流，负载FMU输入电流，避免形成循环依赖。
• 打破代数环：在FMU接口间添加极小值的电容（如pF级）或电感（如nH级），将纯代数环转换为微分代数方程（DAE）系统，降低求解难度。
• 检查FMU导出配置：确保每个FMU导出时选择了正确的因果接口，未保留无因果连接器；同时确认FMU版本（如FMU 2.0 for Co-Simulation）符合仿真环境要求。

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