我正在尝试解决滞后问题。考虑到滞后效应，水分含量取决于历史：考虑一块木材，其初始含水量w0，与phi0的初始环境平衡;然后你把块放到phi1的环境中;现在你需要（w0，phi0和phi1）来确定对应于phi1的w1;然后你把块放到phi2的环境中，确定对应phi2的w2，你需要（w1，phi1和phi2）; .... 在Modelica中，您只能在当前时间步使用状态变量;我怎样才能同时使用phi0和phi1？不是phi，der（phi）的时间导数。当确定w1时，w1和phi1的值将被转移到w0和phi0。 http://vbn.aau.dk/ws/files/13648994/Selected_Constitutive_Models_for_Simulating_the_Hygromechanical_Response_of_Wood 抱歉，我无法上传图片，上传了论文的链接，请查看图6和第65页的公式17（90/135）
对不起，我不知道是否清楚。我想在Modelica中实现吸附滞后。滞后是历史依赖的，我需要先前的状态作为输入来计算当前状态变量。下连续线是吸附曲线，上连续线是解吸曲线，这两条曲线是相对湿度h的函数，它们是已知的; S（H）=（W [H] -wads并[h]）/（wdes并[h] -wads并[h]）。 z0（h0，m0）是初始/先前的水分状态（可以通过初始条件或前一时间步的计算结果给出）。如果当前h大于h0，则zd0（相对湿度hd0）的水分状态将从水分状态z0，那么沿zd0-z0-右箭头的曲线将被定义为hd0的函数;现在我想在定义的曲线zd0-z0-右箭头上计算当前h（h> h0）的水分含量;但我无法在Modelica代码中同时提供h0和h。我将h定义为当前状态，我需要前一步骤h0的水分状态。并且每一个步骤，当计算h处的水分状态时（m为h，m和h的值将转移到m0和h0，以便在下一步计算。）