请教各位老师：
如下图所示，我用mathcad做了个用矩阵位移法计算连续梁内力的模板。当每根梁划分为12个单元时，计算速度不算慢（可能20多秒），但是梁划分到20个单元时，可能要算100秒；我需要知道离支座距离s的内力；如果不在取值位置划分节点（控制每根梁的单元数，再读取该位置的内力，很准），那么就需要用插值法来得到近似值（优点是不受单元数影响，不影响计算速度，缺点是准确度）；
剪力图是直线，可以用线性插值；弯矩图是长度是二次方，可能需要用二次插值；位移图是长度的四次方，这个要如何插值；

仅仅是划分20个单元的连续梁，用有限元求解基本不需要求解时间，怎么还需要20秒左右？？？？？除非梁的受力采用了非常多移动荷载步（并且移动荷载步还带非常多循环）。
任何插值均可用：vs=pspline（vx,vy）
f（x）=interp（vs，vx，vy，x）

我用了多项式插值。也尝试了另一种方法：控制单元长度，在求截面位置放个节点

不到100个节点，确定要20s？