梁自由横振动方程的有限差分方法-重庆交通大学学报.PDF

第 卷第 期重庆 交通 大 学学报 自然科 学版 年 月 梁自由横振动方程的有限差分方法 杜 绍 洪 重庆交通大学 理学院重庆 摘要以简支梁的自由横振动问题为背景求解一个非稳态四阶线性偏微分方程的初边值问题 通过引进辅助函数 组将四阶问题转化为二阶混合初边值问题 对两个辅助函数和二阶混合初边值问题进行离散并消掉中间变量对 由简支梁两端挠度为零得来的二阶偏微分边界条件进行近似处理构造出求解四阶非稳态线性偏微分方程的差分 隐格式 数值实验表明构造的隐格式绝对稳定并且具有很高的精度阶 关键词梁横振动方程有限差分方法差分格式的绝对稳定 中图分类号 文献标志码文章编号 随着新技术和新材料的发展梁以其高强度低 其中为梁的长度 是与梁的体密度梁的横截面 质量的特点在土木桥梁船舶和航空航天等领域逐 积以及抗弯刚度有关的一个正常数 渐得到了广泛的应用 为了实现结构安全和长时间 由于工程实践中初始位移和初始速度通常获得 稳定工作 的是其离散形式 必须对该类构建的横振动进行研究以避 给出解析解的形式 是不可能的 免重大事故的发生 梁发生横振动时由于出现了 即使能给出形式解对级数截取会使得精度低计算 切应力和挠矩方程出现四阶导数项边界条件也要 量大等困难 运用有限元求解可以通过选用相互 相应的增加 若梁的两端固定且端点的挠矩为 正交的高阶形函数可以使单元矩阵成为带宽很窄 即端点与支架是铰接性联接致使端点可以自由旋 的稀疏矩阵从而在实现显著提高计算精度同时计 转但没有位移 设初始位移和初始速度分别为 算量又大为减少的目的但是有限元法待求未知数 和 则梁横振动问题归结为求解如下非 多且要求稀疏矩阵必须保持对称正定性需要求解 稳态四阶线性偏微分方程的初边值问题 的方程规模大导致输入数据多计算的准备工作量 大 尽管有限元方法法适宜处理复杂区域精度可  选但缺点在于内存和计算量巨大不易编程实   现   有限差分方法由于有着简单灵活等优点所以   经常被人们采用 但通常求解最多计算到二阶对   于高阶方程有限差分方法具有相当的困难其表现