有限元法基础与程序设计.pdf

第 1 页 第 章　　平面 问题 的常应变单元 弹性力学经典解法与有 限元法 的不 同特 点 弹性力学经典解法 弹性力学的任务是研 究弹性体在外力作用下而产生应力 、应变和位移 的规律 。解弹性 力学 问题 ，必须考虑平衡微分方程 、几何方程 、物理方程和边界条件 。问题归结为偏微分 方程 的边值 问题 。 以平面弹性力 学 问题为例 ，弹性力学 的基本方程共 个 ： 个 平 衡 方程 ， 个 几 何 方 程 ， 个物理方程 。这 个基本方程 中包含 个未知函数 ，即： 个 位 移 分量 ， 个 应 变 分 量 个应力分量 。基本方程 的数 目等于未知函数 的数 目。弹性力学的任务 ，就是在适当的 边界条件下，从基本方程 中求解这些未知函数 。 弹性力学 问题 的基本解法有三种 ，即按位移求解 、按应力求解和混合求解 。用弹性力 学经典解法解决实际 问题 的主要 困难在于求解偏微分方程 的复杂性 。区域 内各点的位移 、 应变 、应力都是待求 的，因此未知数有无穷多个 。所求解满足弹性力学基本方程 的位移 函 数 、应变函数 、应力 函数 的表达式要覆盖整个 区域 ，而且还要满足边界条件 ，因此求解这 样的函数形式是十分困难的。 有限元法 有限元法是近 年来 随着 电子计算机 的广泛应用而发展起来 的一种数值方法 。它具 有极大 的通用性和灵活性，可 以用来求解弹性力学中的各种复杂边界 问题 。 用有 限元法分析弹性力学 问题 ，首先是把原来连续 的弹性体离散化 。如 图 所 示 的 悬臂梁 ，采用最简单 的三角形单元对弹性体进行分割 ，形成一个如 图 所示 的单元集合 体 。 图 悬臂梁弹性体和有 限元模型 第 2 页 对于每个三角形单元 ，可选择最简单 的线性函数为位移模式，即分片插值 的方法 。单 元 中任一点的位移可通过 个结点的位移进行插值计算 。因此 ，整个 区域 中无穷多个未知 位移量可 以用有限多个结点位移来表示 。这样就避免 了求解覆盖整个 区域 的位移 函数 的困 难 。 用三角形单元 的结点位移 ，可 以表示单元 中的应变 、应力 、结点力 。将各个单元集合 成离散化 的结构模型进行整体分析 ，问题最后归结为求解 以结点位移为未知量 的线性方程 组 。有限元法 中求解这种线性方程组 比弹性力学经典解法 中求解偏微分方程要容 易得 多。 基本未知量和基本方程 的矩阵表示 在有 限元法 中，为 了简洁、清晰地表示各个基本量 以及它们之 间的关系 ，也为 了便于 应用计算机进行实际计算，广泛采用矩阵表示和矩阵运算。 在平面 问题 中，物体所受 的体积力可用列阵表示为 式 中上标 表示矩阵转置 。 同样 ，物体所受 的表面力可用列阵表示为 一点的位移可用列阵表示为 一点的应变分量可用列阵表示为 一点的应力分量可用列阵表示为 由几何方程 ，式（ 所表示 的应变分量可 以写为 （ 弹性力学 中，平面 问题可划分为平面应力 问题和平面应变 问题 。 对 于 弹 性 力 学 的平 面 应力 问题 ，物理方程可用矩 阵形式表示为 （ ） 式中 、 分别为弹性模量和泊松 比。上式可简写为 （