无意中看见的
有弹性力学平衡式子（位移场）、流体连续性方程（压力场）两个场
如果将位移场和压力场分开写pde文件，对于第2式中左端的时间求导项该怎么处理？

如果将位移场和压力场放在一起直接处理，也可以进行计算，但按理说两场分开考虑也没问题，但分开之后第一个方程乘以虚位移处理，第二项则乘以虚压力处理，对于第二式左端的体应变对时间导数项不知该如何处理，高手指点啊，急
如果需要分开写方程，\epsilon_v这项就不能当做自由度了。只能通过coef引入进来。
因为有求导，所以需要读入两个时间步的值。

只不过我对那个前一时间步的引入知道，在前一步那？

上一步引入应该是un，vn，在前一步那？能详细点吗？
比如，位移场有2个自由度，现在需要引入位移场的当前时间步的值和上个时间步的值，
所以PDE文件中就需要引入四个即 coef un vn un1 vn1
然后对时间导数的这项就可以放到方程右端项，当作已知的。
而nfe里面，需要
coef un un1
vect un,un1
read(s,unod) un
read(s,unod1) un1
于是在gcn中
在:1后面，bft前加一行， copy unod unod1
即，首先，一开始unod1和unod一样，以后的时刻步里面，这句话的作用是，
在求解当前步位移场之前，保存上一步的unod结果文件为unod1。
当求解完当前步位移场，我们就有2个时刻步的位移值。