请教关于谐响应分析结果与理论值不相符的问题

最近遇到一个问题，就是利用谐响应分析得到某结构在谐振情况下（简谐力频率与固有频率相等）的谐振位移。但是此位移与理论计算结果不吻合。
结构如图所示，由一个质量块和四根支撑梁组成。质量块的质量m=1.92e-7。支撑梁的尺寸为：梁长L=1000微米，梁宽h=15微米，厚度b=80微米。

这样我们可以根据简单的材料力学知识，得到此系统的一阶频率为4907.2Hz。对该结构进行模态分析，得到一阶频率为：4849.3，基本相符。

然后对该结构进行谐响应分析，频率范围为4844-4854，分100步进行。设置恒定阻尼比的大小为5e-5。

简谐力施加在质量块所有的node上，每node的大小为1e-6N（实部，虚部为0），沿X轴向。

分析完成后，在时间历程后处理中，选择质量块中心node的X位移Ux为变量，绘制曲线。

结果中，谐振幅值为0.023。根据振动理论，在谐振状态下，位移的幅值公式为X=F0/K/(2*DMPRAT)。其中DMPRAT为阻尼比，F0为简谐力幅值，这里就等于1e-6*858（质量块上的节点数）。K是结构刚度，这里可以用频率和质量来反推。这样得到的理论结果X=0.047。与仿真结果有较大差别。此问题困扰多天，跪求大神解救。
另外，按道理说，谐响应分析时，不指定阻尼时，求得的谐振位移应该无穷大。但是我没有设定阻尼时，得到的结果如下图。跪求！