Modelica纯方程化分压电路建模报错问题排查求助

你遇到的结构奇异报错，核心原因是冗余的KCL方程和缺失关键电流约束，以下是具体分析和修正方案：

原代码问题分析 #

1. 冗余的KCL方程 #

FourPin接口已经自带两个TwoPin的电流守恒方程：

• 第一个端口（p1-n1）：0 = p1.i + n1.i
• 第二个端口（p2-n2）：0 = p2.i + n2.i

你手动添加的两个KCL方程：

0 = p1.i + p2.i - iR2;
0 = n1.i + n2.i + iR2;

可以通过接口方程推导得出，属于冗余方程，直接导致系统结构奇异。

2. 缺失电流约束 #

流过R1的电流iR1同时关联p1.i和p2.i，你只定义了p1.i = iR1，但没有约束p2.i与iR1、iR2的关系，导致变量约束不足。

修正后的代码 #

model mathBox_resistorBridge
    extends Modelica.Electrical.Analog.Interfaces.FourPin;

  parameter Modelica.Units.SI.Resistance R1 = 1;
  parameter Modelica.Units.SI.Resistance R2 = 1;

  Modelica.Units.SI.Current iR1 "Current through resistor R1";
  Modelica.Units.SI.Current iR2 "Current through resistor R2";

/*     p1       R1          p2
i1 ->  x------|||||---------x  <- i2
                      |
                     |||
                     |||  R2                    
                     |||
       n1             |     n2
i1 <-  x--------------------x  -> i2
*/
equation 
  // 短接n1和n2，电位相等
  n1.v = n2.v;

  // 单个电阻的电压-电流关系，清晰对应电路结构
  p1.v - p2.v = iR1 * R1;  // R1两端电压降
  p2.v - n1.v = iR2 * R2;  // R2两端电压降（n1与n2等电位）

  // 流入p1的电流等于R1的电流
  p1.i = iR1;

  // p2节点KCL：R1流入的电流 + 外部流入p2的电流 = R2流出的电流
  iR1 + p2.i = iR2;

  annotation (uses(Modelica(version="4.0.0")));
end mathBox_resistorBridge;

修改说明 #

1. 移除了与接口重复的冗余KCL方程，解决结构奇异问题
2. 拆分电压方程为单个电阻的独立VCR，更贴合电路物理意义
3. 补充p2节点的KCL方程，完善变量约束，确保方程数与变量数匹配

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