求解壳体结构的频率特征值

江山美人

我在求解一个壳体结构的频率特征值，先假定了位移函数，位移函数采用多项式求和的形式，理论利用能量法求得了刚度和质量矩阵，再用MATLAB编程求解频率特征值，但是遇到一个问题：当项数取得不是很大的时候，随着项数的增大，频率特征值在慢慢收敛，但是仍未收敛到较精确的结果，而且项数再继续增大，频率特征值开始有点发散，不再收敛。    请问各位大侠，这是什么问题，根据理论，应该是项数越多越收敛才对，难道是因为malab程序计算过程中误差的原因吗？请大家指教，谢谢