单个移动力下winkler梁振动问题

topmny

本人想研究单个移动力作用下的winkler梁振动问题，如下图所示

图1
设梁的E为210Gpa，I为6110*10^-8m⁴，单位长度质量为119.87kg/m，地基支撑刚度k=427kN/m²，荷载P=166.8kN

这样计算得到共振波速为197.6m/s

针对上述参数其他学者取原点左右各100m，计算力以188.65m/s运行至原点左侧1m位置时得到梁上个点位移如下图

图2
他的计算方法为解析公式推导

而我采用有限元法，编制程序将winkler梁左右2端各取100m，假设支撑条件为简支，其他条件相同，计算结果为

图3
各位高人可以指点这是什么原因造成的差异吗？

感觉图2更像静态解，也就是力作用在-1m处引起的位移，而图3和图2差别很大啊？

另外波在介质中传播有提到subcritical, critical,和supercritical，这三种不同波速下，winkler梁都可以用有限元来求解吗？
谢谢