求助：关于单自由度弹簧阻尼系统扫频振动响应的疑惑

xdcyy

紧急求助，请各位好友帮忙分析可能的原因及解决方案！
问题描述：对单自由度弹簧阻尼系统进行扫频振动分析，如图所示，k=60,m=1,c=7.5,激振力f=10sin(2\pi(0.1+4t)t)，初速度和位移均为零。
目标是用MATLAB软件采用直接积分算法Newmark计算位移响应。然而，newmark算法得到的响应与simulink得到的结果完全一致，却与ADAMS仿真结果截然不同。
尝试了很多种可能的方法，包括用ode45求解等，都无法得到与ADAMS或motionview一致的仿真结果。
对于高阶系统差别更大。
各种结果图片如下图所示。

锄头雨

Mathematica用NDSolve函数，龙哥-库塔四阶算出来是这样的，看起也有点不一样。

xdcyy

锄头雨 发表于 2018-12-26 17:20
Mathematica用NDSolve函数，龙哥-库塔四阶算出来是这样的，看起也有点不一样。

谢谢您的帮助，问题已经解决，是我自己出来的问题，在motionview和ADAMS中选择的单位制有点问题，更正之后问题解决。再次感谢您的关注