多自由度参数振动受迫响应的三角级数逼近

针对多自由度的参数振动系统,研究了参数周期与激励力周期不相同情况下受迫振动稳态响应的三角级数解。首先根据调制反馈原理将受迫振动响应表示为各谐波成分线性组合的三角级数形式;然后运用谐波平衡及分块带状矩阵求逆算法,求解得到各个谐波项的系数向量;最后将本文提出的级数逼近法与标准的Runge-Kutta法进行了对比。结果表明:两种方法得到的相图结果高度一致,而本文方法在分析和计算上更具有优势:1)用三角级数来表达振动受迫响应,可以直接得到各频率成分的谐波系数,利于参数振动的时域、频域分析;2)所有的谐波系数向量都可以通过计算机进行数值计算,采用该方法得到的逼近误差与级数项数有关,随着级数项数的增加精度会进一步提高;当级数项数取15项时,逼近误差为8.12×10-3;3)对于同样的10万数据点的计算量,Runge-Kutta法耗时69.51s,而本文采用了大型稀疏矩阵快速算法的级数逼近法计算耗时仅0.1317s,计算效率显著提高。