梁方程降阶计算的重心插值配点法-山东建筑大学学报.pdfVIP

第３２卷 第３期 山东 建 筑 大 学 学 报 Ｖｏｌ．３２ Ｎｏ．３ ２０１７年 ６月 　 　　　　　　　ＪＯＵＲＮＡＬ ＯＦＳＨＡＮＤＯＮＧ ＪＩＡＮＺＨＵ ＵＮＩＶＥＲＳＩＴＹ　　　　　　　 　 Ｊｕｎ． ２０１７ 文章编号：１６７３－７６４４（２０１７）０３－０２４５－０６ 梁方程降阶计算的重心插值配点法 １ １， ２ ３ 徐子康 ，王兆清 ，孙浩森 ，李金 （１．山东建筑大学 工程力学研究所，山东 济南２５０１０１；２．山东建筑大学 学报编辑部，山东 济南２５０１０１；３．山东 建筑大学 理学院，山东 济南２５０１０１） 摘要：采用重心插值配点法求解梁方程时，随着计算节点数量的持续增加，其计算精度将逐步下降。通过对降 阶计算重心插值配点法的研究，可为数值求解梁方程提供一种数值稳定性好、计算精度高的新方法。文章基于 重心Ｌａｇｒａｎｇｅ插值及其微分矩阵，推导了梁方程降阶计算重心插值配点法的公式，并通过数值算例验证其有效 －１０ －１２ 性。结果表明：随着计算节点数量的持续增加，降阶法的计算精度仍保持在１０ ～１０ 范围内；求解两端简支 的梁方程时，两步降阶法的计算精度高于一步降阶法；直接法计算矩阵条件数与节点数的７次方是同阶的，而 一步降阶法计算矩阵条件数与节点数的４次方是同阶的，降阶法可以有效地降低计算矩阵的条件数，提高计算 精度；重心插值配点法采用矩阵—向量形式的计算公式，便于程序的编写，提高了计算效率。 关键词：梁方程；降阶法；重心Ｌａｇｒａｎｇｅ插值；配点法 中图分类号：Ｏ３０２　　　　　文献标识码：Ａ Ｂａｒｙｃｅｎｔｒｉｃｉｎｔｅｒｐｏｌａｔｉｏｎｃｏｌｌｏｃａｔｉｏｎｍｅｔｈｏｄｂａｓｅｄ ｏｎｄｅｐｒｅｓｓｉｏｎｏｆｏｒｄｅｒｆｏｒｓｏｌｖｉｎｇｂｅａｍｅｑｕａｔｉｏｎｓ １ １， ２ ＸｕＺｉｋａｎｇ，ＷａｎｇＺｈａｏｑｉｎｇ ，ＳｕｎＨａｏｓｅｎ，ｅｔａｌ． （１．ＩｎｓｔｉｔｕｔｅｏｆＭｅｃｈａｎｉｃｓ，ＳｈａｎｄｏｎｇＪｉａｎｚｈｕＵｎｉｖｅｒｓｉｔｙ，Ｊｉｎａｎ２５０１０１，Ｃｈｉｎａ；２．ＥｄｉｔｏｒｉａｌＤｅｐａｒｔｍｅｎｔｏｆＪｏｕｒｎａｌｏｆ ＳｈａｎｄｏｎｇＪｉａｎｚｈｕＵｎｉｖｅｒｓｉｔｙ，Ｊｉｎａｎ２５０１０１，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｗｈｅｎｔｈｅｂｅａｍｅｑｕａｔｉｏｎｓａｒｅｓｏｌｖｅｄｂｙｂａｒｙｃｅｎｔｒｉｃｉｎｔｅｒｐｏｌａｔｉｏｎｃｏｌｌｏｃａｔｉｏｎｍｅｔｈｏｄ， ｃｏｍｐｕｔａｔｉｏｎａｌａｃｃｕｒａｃｙｗｉｌｌｄｅｃｌｉｎｅｇｒａｄｕａｌｌｙａｓｔｈｅｎｕｍｂｅｒｏｆｎｏｄｅｓｉｎｃｒｅａｓｅｓ．Ｔｈｅｒｅｓｅａｒｃｈｏｎ ｂａｒｙｃｅｎｔｒｉｃｉｎｔｅｒｐｏｌａｔｉｏｎｃｏｌｌｏｃａｔｉｏｎｍｅｔｈｏｄｂａｓｅｄｏｎｄｅｐｒｅｓｓｉｏｎｏｆｏｒｄｅｒ，ｃａｎｐｒｏｖｉｄｅｔｈｅｎｅｗｍｅｔｈｏｄ ｔｈａｔｈａｓａｇｏｏｄｎｕｍｅｒｉｃａｌｓｔａｂｉｌｉｔｙａｎｄｈｉｇｈｃｏｍｐｕｔａｔｉｏｎａｌａｃｃｕｒａｃｙｆｏｒｂｅａｍｅｑｕａｔｉｏｎｓ．Ｂａｓｅｄｏｎ ｂａｒｙｃｅｎｔｒｉｃＬａｇｒａｎｇｅｉｎｔｅｒｐｏｌａｔｉｏｎａｎｄｉｔｓｄｉｆｆｅｒｅｎｔｉａｌｍａｔｒｉｃｅｓ，ｔｈｅｆｏｒｍｕｌａｏｆｂａｒｙｃｅｎｔｒｉｃｉｎｔｅｒｐｏｌａｔｉｏｎ ｃｏｌｌｏｃａｔｉｏｎｍｅｔｈｏｄｂａｓｅｄｏｎｄｅｐｒｅｓｓｉｏｎｏｆｏｒｄｅｒｉｓｄｅｒｉｖｅｄ．Ｎｕｍｅｒｉｃａｌｅｘａｍｐｌｅｓａｒｅｇｉｖｅｎｔｏｖｅｒｉｆｙｔｈｅ ｅｆｆｅｃｔｉｖｅｎｅｓｓｏｆｔｈｅｐｒｏ