用于有限元分析的MATLAB函数.doc

? PAGE 8 ? MATLAB有限元分析与应用 ? PAGE 9 ? 第1章 引 言 第1章 引 言 这个简短的引言分为两部分，第一部分是对有限元方法步骤的概括介绍，第二部分是MATLAB的简略使用指南。 1.1 有限元方法的步骤 有许多关于有限元分析的优秀教材，比如在参考文献[1-18]中列出的那些书目。因此，本书不准备对有限元理论或有限元方程进行详细地阐述和推导。每一章只总结概括主要的方程，这些章节都附有示例来说明这些方程。此外，全书只讨论线弹性结构力学的问题。 有限元方法用于解决工程问题的数值计算过程。本书假定所有的行为都是线弹性行为。虽然本书的问题都与结构工程相关，但有限元方法也同样适用于工程的其他领域。本书中使用有限元方法解决问题共包括6个步骤。对有限元分析的6个步骤阐述如下： (1) 离散化域—— 这个步骤包括将域分解成单元和节点。对于像桁架和刚架这类离散系统，已经离散化，这一步就不需要了。此处获得的结果应该已经是精确的。然而，对于连续系统，如板壳，这一步就变得至关重要，因为它只能得到近似的结果。因此解决方案的精确度取决于所使用的离散化方法。本书中，我们将手动完成这一步(对连续系统)。 (2) 写出单元刚度矩阵(element stiffness matrices)—— 写出域内每个单元的单元刚度矩阵。在本书中，这个步骤通过MATLAB实现。 (3) 集成整体刚度矩阵(global stiffness matrices)—— 这一步用直接刚度法(direct stiffness approach)实现。在本书中，该步骤借助于MATLAB实现。 (4) 引入边界条件—— 诸如支座(supports)、外加载荷(applied loads)和位移(displacements)等。本书中手动实现这一步骤。 (5) 解方程—— 这一步骤分解整体刚度矩阵并用高斯消去法求解方程组。在本书中，在用高斯消去法实现求解部分的时候需要手动分解矩阵。 (6) 后处理—— 得到额外的信息，如支反力、单元节点力和单元应力。本书中这一步骤通过MATLAB实现。 从以上步骤可以看出，解决问题的过程结合使用了MATLAB和某些有限的手动操作。使用手动操作只是实现离散化域(步骤1)、设置边界条件(步骤4)和分解整体刚度矩阵(步骤5)。可以看出，所有冗长、反复的计算都由MATLAB完成。 1.2 用于有限元分析的MATLAB函数 我们为本书提供的免费资源文件(请登录站点的下载页面)中包含了75个MATLAB函数(m文件)，主要是由本书的作者专门为书中的有限元分析而编写的。它们组成了MATLAB有限元工具箱(MATLAB Finite Element Toolbox)。下面列出了所有这些函数。对于每个函数的详细使用方法，读者可以参考每一章的内容。 SpringElementStiffness(k) SpringAssemble(K，k，i，j) SpringElementForces(k，u) LinearBarElementStiffness(E，A，L) LinearBarAssemble(K，k，i，j) LinearBarElementForces(k，u) LinearBarElementStresses(k，u，A) QuadraticBarElementStiffness(E，A，L) QuadraticBarAssemble(K，k，i，j，m) QuadraticBarElementForces(k，u) QuadraticBarElementStresses(k，u，A) PlaneTrussElementLength(x1，y1，x2，y2) PlaneTrussElementStiffness(E，A，L，theta) PlaneTrussAssemble(K，k，i，j) PlaneTrussElementForce(E，A，L，theta，u) PlaneTrussElementStress(E，L，theta，u) PlaneTrussInclinedSupport(T，i，alpha) SpaceTrussElementLength(x1，y1，z1，x2，y2，z2) SpaceTrussElementStiffness(E，A，L，thetax，thetay，thetaz) SpaceTrussAssemble(K，k，i，j) SpaceTrussElementForce(E，A，L，thetax，thetay，thetaz，u) SpaceTrussElementStress(E，L，thetax，thetay，thetaz，u) BeamElementStiffness(E，I，L)