非齐次动力学方程的一种精细积分单步方法

来源：计算力学学报杂志2020年第2期北京时间：

作者：王永,马骏,李靖翔,陈汝斯,许杨,郝跃东,胡鹏

单位：国网上海市电力公司特高压换流站分公司, 上海 201413,国网湖北省电力有限公司, 武汉 430070,南方电网超高压公司广州局, 广州 510000,国网湖北省电力有限公司电力科学研究院, 武汉 430070,国网上海市电力公司特高压换流站分公司, 上海 201413,国网上海市电力公司特高压换流站分公司, 上海 201413,国网上海市电力公司特高压换流站分公司, 上海 201413

摘要：针对非齐次动力学方程v=Hv+f（v，t），结合精细积分法和微分求积法，利用同阶的显式龙格-库塔法对计算过程中待求的vk+i/s（i=1，2，…，s）进行预估，提出了一种避免状态矩阵求逆的高效精细积分单步方法。该方法采用精细积分法计算eHt，而Duhamel积分项采用s级s阶的时域微分求积法，计算格式统一且易于编程，可灵活实现变阶变步长。仿真结果表明，与其他单步法及预估校正-辛时间子域法进行数值比较，该方法具有高精度、高效率及良好的稳定性，在求解大规模动力系统时间响应问题中具有较大的优势。

关键词：非齐次动力学方程,精细积分法,微分求积法,变阶,单步法,