Modelica带硬限位滑阀位移超限问题及修复方案

问题背景 #

我正在开发一款弹簧压紧式提升阀模型，阀门开启逻辑为：当入口侧压力产生的力大于弹簧力与出口侧压力产生的力之和时，阀门开启。通过编写控制器，以进出口流体压力为输入，输出范围为[0,1]的阀门开度信号。

参考Modelica Standard Library（MSL）中的带限位滑动质量组件，引入了stopped变量（位移s < smin时为-1，s > smax时为1），计划触发限位时重新初始化质量的位置和速度，但运行时发现位移s始终远超设定的smax值，限位逻辑完全失效。

解决方法 #

将部分when语句从equation段移至algorithm段后，限位逻辑正常生效，问题得到解决。

代码对比 #

失效版本（限位逻辑无效） #

model SpringLoadedValve_Broken
  // 变量定义
  Real s(start=0.0); // 阀芯位移
  Real v(start=0.0); // 阀芯速度
  parameter Real smin = 0.0; // 最小位移限位
  parameter Real smax = 0.01; // 最大位移限位
  Integer stopped = 0; // 限位状态标记
  // 输入输出
  Real p_in, p_out; // 进出口压力
  output Real opening(start=0.0); // 阀门开度

  // 物理参数
  parameter Real A = 0.001; // 阀芯受力面积
  parameter Real k = 10000; // 弹簧刚度
  parameter Real c = 10; // 阻尼系数
  parameter Real m = 0.1; // 阀芯质量
equation
  // 阀芯动力学方程
  der(s) = v;
  der(v) = (p_in*A - p_out*A - k*s - c*v)/m;

  // 原限位逻辑：when语句放在equation段导致失效
  when s < smin then
    stopped = -1;
    reinit(s, smin);
    reinit(v, 0.0);
  end when;
  when s > smax then
    stopped = 1;
    reinit(s, smax);
    reinit(v, 0.0);
  end when;

  // 开度计算：映射到[0,1]范围
  opening = clamp(s/smax, 0.0, 1.0);
end SpringLoadedValve_Broken;

修复版本（限位逻辑正常） #

model SpringLoadedValve_Fixed
  // 变量定义
  Real s(start=0.0); // 阀芯位移
  Real v(start=0.0); // 阀芯速度
  parameter Real smin = 0.0; // 最小位移限位
  parameter Real smax = 0.01; // 最大位移限位
  Integer stopped = 0; // 限位状态标记
  // 输入输出
  Real p_in, p_out; // 进出口压力
  output Real opening(start=0.0); // 阀门开度

  // 物理参数
  parameter Real A = 0.001; // 阀芯受力面积
  parameter Real k = 10000; // 弹簧刚度
  parameter Real c = 10; // 阻尼系数
  parameter Real m = 0.1; // 阀芯质量
equation
  // 阀芯动力学方程
  der(s) = v;
  der(v) = (p_in*A - p_out*A - k*s - c*v)/m;

  // 开度计算：映射到[0,1]范围
  opening = clamp(s/smax, 0.0, 1.0);
algorithm
  // 修复：将when语句移至algorithm段，限位逻辑生效
  when s < smin then
    stopped := -1;
    reinit(s, smin);
    reinit(v, 0.0);
  end when;
  when s > smax then
    stopped := 1;
    reinit(s, smax);
    reinit(v, 0.0);
  end when;
end SpringLoadedValve_Fixed;

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