如何对Modelica变量积分且不影响动态状态选择

方案合理性判断 #

你给出的实现是Modelica领域对变量做时间积分的通用常规方案，完全符合物理定义逻辑，量纲匹配、可读性强，仅用于后处理绘图的场景下完全可用。

更优替代方案 #

分使用场景有两类更优选择：

• 仅需积分结果做后处理、不需要参与仿真过程中的逻辑计算：不需要修改模型代码，主流Modelica仿真工具（Dymola、OpenModelica等）均支持仿真结束后在结果后处理模块直接对P_comp变量做时间积分运算，轻量化无负担，不会增加模型的编译和仿真开销。
• 需要把积分结果保留在模型中：可以在原有代码基础上做优化，也可以直接调用标准库模块：
  • 优化基础实现：给积分变量补充初始值定义，避免工具报初始条件缺失警告：

    Modelica.Units.SI.Power P_comp;
    Modelica.Units.SI.Energy E_comp(start=0, fixed=true);
    equation
    der(E_comp) = P_comp;
    

  • 有复位、限幅等额外需求：直接调用Modelica.Blocks.Continuous.Integrator标准库积分模块，内置了各类常用功能，无需自己手写方程，可复用性更高。

对模型动态状态选择的影响 #

默认情况下该实现会把E_comp识别为新增的动态状态变量，但因为这个状态是纯线性的、和其他非线性变量无耦合，对绝大多数中小规模模型的仿真效率影响可以忽略。如果要降低对原有模型状态选择的影响，可以在E_comp的声明中添加stateSelect=StateSelect.avoid修饰符，提示仿真工具优先选择原有模型的状态变量，仅当没有其他可选状态时才将E_comp选为动态状态。

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