两点边值问题的Fortran-Matlab程序实现与误差分析

需积分: 19 9 下载量 66 浏览量 更新于2024-08-08 1 收藏 787KB PDF 举报
本文主要探讨了两点边值问题在有限元方法中的应用,并提供了详细的程序实现与计算分析。该研究基于山东省优秀中青年科学家科研奖励基金和自然科学基金的支持,针对的是一个典型的边界值问题,即Lu=-d的求解。论文首先介绍了两点边值问题作为有限元方法学习中的基础模型,强调了理论分析的重要性,包括变分形式、有限元空间构建、误差估计和收敛性。 然而,作者指出在实际的编程和计算展示方面,相关文献相对较少,且缺乏完整的源代码供读者参考。因此,本文的主要贡献在于提供了一种实用的方法,通过Fortran 77和Matlab语言编写程序,实现了线性有限元方法,并展示了有限元解与真解在L2和H1范数下的误差和收敛阶。这种程序设计不仅展示了理论知识的实践应用,而且对于研究生和本科生来说,具有重要的指导意义,能提升他们使用计算机程序实现数值方法的技能和处理计算结果的能力。 在程序实现部分,作者详细解释了模型选择、子程序设计、函数定义与调用、数据管理和计算结果的可视化技术。通过这些技术,读者不仅能理解有限元方法的实施过程,还能学会如何优化代码和有效展示计算结果,这对于教学和科研人员来说具有很高的参考价值。 通过对比有限元解和真解的误差及收敛性,文章验证了理论分析的有效性,证明了有限元方法在解决两点边值问题上的有效性。这不仅增强了读者对有限元方法的理解,也为实际工程问题的数值模拟提供了实用工具。 本文是一篇结合理论和实践的重要研究,旨在填补两点边值问题有限元方法程序实现方面的研究空白,为相关领域的专业人士提供了宝贵的学习资源。
2019-03-19 上传
本书的第1 章简要地介绍了 MATLAB 的基本知识和编程中常用的语句及函数 , 使者能够阅读本书各章节中的程序。 第2 章系统地介绍了有限元的理论基础———微分方程的近似解法。 这部分内容在一的有限元书籍中是很少介绍的 , 它不仅可以使我们了解有限元的发展过程 , 也能够使读者加深对有限元方法的理解。 第3 章介绍了广义坐标有限元方法。 它是物理坐标下的直接方法 , 读者可以通过该章的学习了解和掌握有限元方法的一般步骤。 第4 章简要介绍了有限元编程方法。 与大多数有限元书籍不同的是 , 用其他高级语言编写有限元程序时所需的一般编程技巧在 MATLAB 中不再需要 , 因此 , 本书不再赘述。 第5 章详细讨论了构造单元和插值函数的原则和方法 , 并着重讨论了在实际中有着广泛应用的等参数单元的构造方法和表达格式 , 以及与广义坐标单元的变换方法。 第6 章和第7 章讨论了杆系结构有限元问题。由于杆系结构与一般的二维和三维弹性体结构有较大的区别 , 因此 , 杆系结构的单元及其插值函数是区别于一般二维和三维单元的特殊单元 , 同时 , 桁架的杆单元和框架的梁单元也是完全不同的两类单元。 第8 章详细讨论了一般弹性力学问题的有限元方法 , 包括稳定问题和动力学问题。 第9 章讨论了板问题的有限元方法 , 其中介绍了多种类型和不同位移模式的板单元 ,包括用于复合材料结构的层状单元。 第10 章介绍了系统建模、 线性系统分析及结构振动控制的基础知识 , 并详细地介绍了如何用 MATLAB 来实现。