关键词:
梁振动偏微分方程
微分–代数方程
L-稳定方法
非线性
稳定性
摘要:
基于高阶非线性梁振动偏微分方程的一般形式,构造了数值求解的L-稳定格式。首先,选取三角插值基函数,基于插值定理进行空间离散,将带有初边值条件的偏微分方程求解问题转化为微分–代数方程求解。然后在时间区间上构造L-稳定求解格式进行求解。以无轴向运动简支梁在外部激励下的强迫振动方程为例进行数值仿真,对梁的位移轨迹、边界条件及系统能量进行探究,并与龙格–库塔法、微分求积法进行对比,结果表明,L-稳定方法可以在较大步长下满足边界,位移轨迹与模型方程一致,在计算精度和稳定性上都有较好的体现。