基于Matlab语言的有限元法及其应用.PDF

维普资讯 烟台师范学院学报(自然科学版) YantaiNormalUniversityJournal(NaturalScience) 2004。20(2)I10一10 厂 、‘ t综述t - + 基于 Matlab语言的有 限元法及其应用 宋克志 一，刘智儒 (1．烟台师范学院交通学院，山东 烟台 264025；2．北京交通大学隧道及地下工程试验研究中心 ，北京 100044) 摘要：Or绍了有限元法的本质特征及用变分试函数法和残值试函数法导出有限元法的过程，给出了在Matlab 语言环境下实现有限元法的步骤．利用Matlab语言中的PDE工具箱求解偏微分方程具有简便、快速、可视化 程度高等优点，能满足精度要求，并以一个工程实例说明了利用有限元法求解偏微分方程从而解决实际问题 的方法 ． 关键词：有限元法；偏微分方程；Matlab语言 中图分类号：O241．82；TB115 文献标识码：A 文章编号2004)02．0100-03 在解决各类工程问题时，有限元法发挥着重 中水头与流量之间的关系、电势与电量之间的关 要作用 ，而许多工程问题 的解决是利用微分方程 系等都可称为刚度矩阵． 来求解的，大多数微分方程求解析解较困难 ，而求 有限元法可以由多种方式推导出来 ，若忽略 数值解较容易，有限元法是求解微分方程的方法， 推导过程 ，此方法的结果就是微分方程的数值解 ． 它是数值分析中Rayleigh．Ritz方法的具体应用． 有限元法可以由变分试函数法 (也称 Ritz法)和残 本文在 阐明有 限元法与微分方程关系的基础上， 值试函数法 (也称 Galerkin法)，这两种方法是有 结合工程实例，给出了在 Matlab语言环境下实现 限元法的基础 ．这两种方法开始都要假设试函数 ， 有限元法的步骤 ． 基本原理如图 1所示 ，从 中可以看出有限元法与 微分方程的关系 ．这两种方法在数学上是等价的， 1 有限元法及其推导 可以证明，Galerkin法和变分法可导出相等的公 有限元法是将物体划分成有限个单元，这些 单元之间通过有限个节点相互连接，单元看作是 不可变形的刚体，单元之间的力通过节点传递，然 后利用能量原理建立各单元矩阵，然后组成总刚 矩阵，求解线性方程组．结构分析中经典的刚度方 法可以在不考虑控制微分方程的情况下导出来 ， 但是推导刚度方法 的基本关系是 以微分方程的解 为基础的，也即能量原理与控制微分方程之间是 有必然联系的，使用时此原始含义很容易被忽 略…．实际上 ，微分方程渗透到工程领域的方方面 面，常见的控制微分方程有弹性问题、热传导、位 I 势流、质量输运、电流等．结构分析中的刚度矩阵 匝圃 用以描述力和位移之间的矩阵关系，但是现在提 I 到刚度矩阵时可 以不考虑其具体应用，如热传导 匾亟巫 固 问题中温度与热通量之间的矩阵关系、渗流问题 图1 有限元求解微分方程的两种方法 收稿 E1期 ：20∞．11．13 作者简介 ：宋克志(197O一)，男 ，讲师 ，在读博士研究生，主要从事隧道及地下工程方面的研究 ． 维普资讯 第 2期 宋克志，等：基于Matlab语言的有限元法及其应用 101 式…．它们之间的差异主要在于满足边界条件的 还未铺设 防水层)及隧道边界 3(开挖面)为透水 要求不同．Galerkin法