OpenModelica中Delay模块奇异矩阵与全主元求解器切换问题咨询

问题现象 #

开发大型模型时持续触发以下报错，求解器自动切换为全主元求解器后仿真可完成：

Failed to solve linear system of equations (no. 5) at time 0.000000, system is singular for U[2, 2].
The default linear solver fails, the fallback solver with total pivoting is started at time 0.000000. That might raise performance issues, for more information use -lv LOG_LS.
Matrix singular!

The initialization finished successfully without homotopy method.
### STATISTICS ###
The simulation finished successfully.

最小可复现模型（MWE） #

model mwe 
  
  Modelica.Blocks.Nonlinear.FixedDelay fixedDelay(delayTime = 1) annotation(
    Placement(visible = true, transformation(origin = {10, -80}, extent = {{10, -10}, {-10, 10}}, rotation = 0)));
  Modelica.Blocks.Sources.BooleanConstant booleanConstant(k = false) annotation(
    Placement(visible = true, transformation(origin = {-128, -40}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
  Modelica.Blocks.Sources.Constant const7(k = 1) annotation(
    Placement(visible = true, transformation(origin = {-88, -22}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
  Modelica.Blocks.Logical.Switch switch18 annotation(
    Placement(visible = true, transformation(origin = {-30, -40}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));

equation
  connect(switch18.u2, booleanConstant.y) annotation(
    Line(points = {{-42, -40}, {-117, -40}}, color = {255, 0, 255}));
  connect(switch18.y, fixedDelay.u) annotation(
    Line(points = {{-19, -40}, {60, -40}, {60, -80}, {22, -80}}, color = {0, 0, 127}));
  connect(fixedDelay.y, switch18.u3) annotation(
    Line(points = {{-1, -80}, {-60, -80}, {-60, -48}, {-42, -48}}, color = {0, 0, 127}));
  connect(const7.y, switch18.u1) annotation(
    Line(points = {{-77, -22}, {-60, -22}, {-60, -32}, {-42, -32}}, color = {0, 0, 127}));
  
annotation(
    preferredView = "diagram",
    Icon(graphics = {Rectangle(fillColor = {255, 255, 255}, fillPattern = FillPattern.Solid, extent = {{-100, 100}, {100, -100}}), Text(origin = {6, 25}, extent = {{-52, 59}, {42, -47}}, textString = "Delay"), Text(origin = {-2, -31}, extent = {{-38, 45}, {42, -47}}, textString = "Flag")}, coordinateSystem(initialScale = 0.1)));
end mwe;

关键疑问 #

• 仅将FixedDelay输入输出直接连接时，触发奇异矩阵报错且仿真失败；加入Switch后的MWE切换求解器后可运行，二者差异原因是什么？
• 尝试添加初始方程fixedDelay.y = 1;未解决问题，如何让模型正常编译运行并消除报错？

问题根源与解决方案 #

根源分析 #

1. FixedDelay的初始化缺陷：FixedDelay在t=0时的输出默认无明确初始值，直接闭环连接输入输出会形成无约束代数环，线性系统因无法唯一确定初始值而奇异（矩阵不可逆），导致默认求解器崩溃。
2. 带Switch模型的运行逻辑：初始时booleanConstant.k=false，Switch会选择u1的固定值（const7.k=1）作为输出传递给FixedDelay输入，此时FixedDelay的输入在t=0时有明确约束，代数环被打破。全主元求解器可处理剩余弱耦合结构，因此仿真能完成，但默认求解器仍会检测到潜在奇异结构触发报错。

解决方案 #

1. 正确初始化FixedDelay状态

不要直接赋值fixedDelay.y，而是约束初始时刻的输入输出关系，打破无约束代数环：

initial equation
  // 让FixedDelay初始输出等于初始输入，明确约束
  fixedDelay.y = fixedDelay.u;

若模型支持，也可通过FixedDelay的initialOutput参数直接指定初始输出值。

2. 避免初始时刻的无约束闭环

对于直接闭环的FixedDelay模型，需给初始输入明确赋值：

model DirectDelayLoop
  Modelica.Blocks.Nonlinear.FixedDelay fixedDelay(delayTime=1);
equation
  connect(fixedDelay.y, fixedDelay.u);
initial equation
  fixedDelay.u = 1; // 明确初始输入值，消除奇异环
end DirectDelayLoop;

3. 强制使用全主元求解器（临时方案）

在OpenModelica仿真设置中添加启动参数 -ls total，强制使用全主元线性求解器，避免切换报错。但该方案可能降低大型模型仿真性能，仅作临时规避手段。

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