用于非线性动力分析的一种高效精细积分单步法

针对非线性动力状态方程v=H·v+f(v,t),结合精细积分法和Romberg数值积分,对计算过程中待求的v_(k+j/m)(j=1,2,…,m),利用当前时刻vk,通过二阶龙格库塔法进行预估,提出了一种避免状态矩阵求逆的高效精细积分单步法。该方法计算格式统一,易于编程,通过选取m值,可进行不同计算精度的计算。与两种单步法、一次预-校法及预估校正-辛时间子域法进行数值比较,计算结果表明,该方法具有高精度、高效率及较好的稳定性。在求解多自由度、强非线性动力响应问题中具有较大优势。