OpenModelica求解平面摆DAEs时初始条件不一致问题咨询

问题明确 #

使用OpenModelica求解平面摆DAE模型时，出现以下异常：

• 当设置正x方向初始坐标（如x0=0.5, y0=0，满足摆长约束x²+y²=L²=0.25），模型可正常运行；
• 当设置负x方向初始坐标（如x0=-0.5, y0=0，同样满足约束），系统提示初始条件不一致；
• 其他满足约束的负x初始组合（如x0=-0.013437983982246, y0=-0.499819387965786）也触发相同错误。

原因剖析 #

OpenModelica处理DAE初始条件时，会检查约束方程的数值相容性。对于平面摆的约束x²+y²=L²：

1. 若未显式指定速度约束，求解器默认的初始迭代猜测方向可能偏向正x区域，与负x初始值产生数值冲突；
2. 部分版本的求解器在处理符号突变的初始约束时，会因雅可比矩阵的局部奇异性，误判初始条件不一致。

可行解决方案 #

1. 显式添加速度约束 #

对位置约束求导得到速度约束x*dx + y*dy = 0，将其加入模型，帮助求解器明确初始状态的相容性：

model PlanarPendulum
  parameter Real L = 0.5;
  parameter Real g = 9.81;
  Real x(start=-0.5), y(start=0);
  Real dx(start=0), dy(start=0);
equation
  // 位置约束
  x^2 + y^2 = L^2;
  // 速度约束（对位置约束求导）
  x*dx + y*dy = 0;
  // 运动方程
  der(dx) = -x*(g/L);
  der(dy) = -y*(g/L) - g;
end PlanarPendulum;

2. 强制锁定初始值 #

通过fixed=true属性强制求解器使用给定的初始值，避免自动猜测导致的冲突：

Real x(start=-0.5, fixed=true), y(start=0, fixed=true);
Real dx(start=0, fixed=true), dy(start=0, fixed=true);

3. 调整求解器参数 #

在仿真配置中选择对DAE更鲁棒的求解器（如DASSL或IDA），并适当调大初始条件容忍度：

• 打开OpenModelica仿真设置，将Solver设为DASSL；
• 在Initialization选项中，将Tolerance调整为1e-6或更高。

4. 改用角度参数化建模（推荐） #

避开DAE的初始条件问题，直接用角度θ参数化位置，转为ODE建模：

model PlanarPendulum_Angle
  parameter Real L = 0.5;
  parameter Real g = 9.81;
  Real theta(start=3.14159265); // 对应x=-0.5,y=0的初始角度
equation
  der(der(theta)) = - (g/L)*sin(theta);
  // 导出x、y坐标
  Real x = L*cos(theta);
  Real y = L*sin(theta);
end PlanarPendulum_Angle;

验证说明 #

按上述方案调整后，负x方向的初始条件均可正常运行。其中角度参数化建模是最简便的规避方式，显式添加速度约束则更贴合DAE建模的原始需求。

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