Modelica中基于求解器时间步的Variable Clock事件实现问询

要实现不干扰求解器自适应步长、每个时间步都触发事件的需求，核心是不要用固定周期时钟，而是绑定求解器自身的时间步完成逻辑，以下是具体可行的方法：

• 直接利用求解器的内置回调/钩子
  绝大多数仿真和数值求解框架都提供了「每个时间步完成后执行自定义代码」的原生接口，这是最稳妥的方式：

  • MATLAB Simulink：在模型设置的「Solver Callback」中配置PostStepFcn，指定你要运行的算法函数。求解器会在每一步计算完成后自动触发这个函数，完全不会影响它的自适应步长选择逻辑。
  • Python SciPy：使用scipy.integrate.solve_ivp时，设置t_eval=None（让求解器自主选择步长），然后通过callback参数传入你的处理函数，注意回调函数返回False才能让求解继续推进。
  • Modelica/OpenModelica：可以用when time > pre(time)作为事件触发条件——每次求解器推进时间步后，time会更新为新值，这个条件会立刻触发，而且不会强制求解器缩小步长，因为它是被动响应时间变化，而非预设固定触发点。

• 自定义事件触发条件（适配无原生回调的框架）
  如果你的框架没有直接的步长回调，可以手动设计一个基于时间变化的触发逻辑：

  1. 定义一个状态变量last_recorded_time，初始值设为仿真起始时间。
  2. 设置事件触发条件为time > last_recorded_time。
  3. 每次事件触发后，立即更新last_recorded_time = time，确保下一次只有当求解器推进了新的步长才会再次触发。
     这种方式本质是被动检测求解器的时间推进，不会干涉求解器的步长决策。

• 关键注意事项

  • 不要在事件处理函数中修改求解器的步长参数或核心状态变量，否则可能打乱自适应步长的计算逻辑。
  • 区分「求解器的完成步长」和「内部迭代步」：有些隐式求解器会在一个时间步内进行多次迭代，要确保触发的是最终完成的那个步长，优先用框架原生的步长回调来避免这个问题。

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