Pyomo DAE仿真报错：声明单个DerivativeVar却提示存在多个

问题核心 #

你在使用Pyomo Simulator仿真DAE模型时，明明只声明了1个DerivativeVar（m.da对应m.a的导数），却触发了Cannot simulate a differential equation with multiple DerivativeVars异常；而给第二个变量m.b添加DerivativeVar后，异常反而消失，这一现象的根源在于你的代数方程存在奇异表达式，导致Pyomo对DAE结构的分析出现误判。

具体原因分析 #

你的代数方程eq1定义为：

def eq1(m,t):
    return (m.a[t]-m.b[t])/(m.a[t]**(-4)) == (m.b[t]-m.a[t])

这个表达式存在两个关键问题：

1. 奇异点：当m.a[t] = 0（你的初始值就是0）时，m.a[t]**(-4)会导致分母为0，触发数值奇异。
2. 冗余化简：对等式整理后可发现，左边等价于(m.a[t]-m.b[t]) * m.a[t]**4，移项后得到(m.a[t]-m.b[t])*(m.a[t]**4 + 1) = 0。由于m.a[t]**4 + 1恒大于0，等式实际等价于m.a[t] == m.b[t]，原表达式完全是冗余且易出错的写法。

Pyomo在分析DAE结构时，遇到这个奇异的代数方程，无法正确识别其为纯代数约束，反而误判系统需要额外的DerivativeVar来满足方程，因此抛出“多个DerivativeVar”的错误。而当你手动给m.b添加DerivativeVar后，Pyomo的结构分析逻辑被改变，不再触发这个误判，但这属于“治标不治本”的做法——你的模型根本不需要m.b的导数。

正确解决方法 #

直接化简代数方程，消除奇异表达式，让DAE结构清晰可辨：

import pyomo.environ as pyo
import pyomo.dae as dae

m = pyo.ConcreteModel()
m.t = dae.ContinuousSet(bounds=(0,1))

m.a = pyo.Var(m.t, initialize=0)
m.b = pyo.Var(m.t, initialize=0)
m.da = dae.DerivativeVar(m.a, initialize=0)

# 化简后的代数方程，彻底消除奇异点
def eq1(m,t):
    return m.a[t] == m.b[t]
m.eq1_con = pyo.Constraint(m.t, rule=eq1)

def eq2(m,t):
    return (m.da[t] * 10) + m.b[t] == 1
m.eq2_con = pyo.Constraint(m.t, rule=eq2)

sim = dae.Simulator(m, package='casadi')
tsim, profiles = sim.simulate(integrator='idas', numpoints=1000)

此时模型结构明确：1个微分方程（含m.da）、1个代数方程（a=b）、1个DerivativeVar，完全符合DAE仿真的要求，不会再触发异常。

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