Dymola与OpenModelica运行Modelica代码结果不一致问题咨询

在Dymola 2022和OpenModelica v1.19.2中运行以下Modelica测试代码时，得到了不同的结果。实验仅运行至1秒，if time > 10代码块不会执行，预期b取a的初始值{0,0,0,0}，但实际输出存在明显差异：

model Test
  Real a[2,2](each start=0);
  Real b[2,2];
algorithm 
  b := a;
  a := [1,2; 2,4];
  if time > 10 then
    a := b;
  end if;
end Test;

结果差异对比 #

• Dymola输出
  • b变量值：
  • a变量值：
• OpenModelica输出
  • b变量值：
  • a变量值：

差异原因 #

1. Modelica规范的执行逻辑：Modelica标准明确规定，algorithm段内的语句是严格按顺序执行的——先执行b := a（此时a为初始值0），再执行a := [1,2; 2,4]，因此b应保留a的初始值，Dymola的行为符合规范。
2. OpenModelica旧版本的实现缺陷：OpenModelica v1.19.2属于较旧版本，在处理algorithm段的顺序赋值时存在逻辑错误，错误地将两个赋值语句当成联立方程处理，导致b被赋值为a更新后的值，违背了顺序执行的要求。

解决方案 #

• 升级OpenModelica版本：更新至最新稳定版（如v1.21及以上），新版本已修复此类顺序执行的逻辑问题，能严格按照algorithm段的语句顺序执行赋值。
• 兼容旧版本的代码调整：若需适配旧版本，可通过显式初始化或引入临时变量确保逻辑正确，示例代码如下：

model Test
  Real a[2,2](each start=0);
  Real b[2,2](each start=0);
  Real temp_a[2,2];
algorithm 
  temp_a := a;
  b := temp_a;
  a := [1,2; 2,4];
  if time > 10 then
    a := b;
  end if;
end Test;

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