非线性固体计算力学 宋天霞等著.pdf

第 1 页 第 章　　线性有 限元法 线性有 限元法 的基本理论 ，早 已为人们所 熟悉并被广泛 应用 ，因此 ，在这里不再系统地论述 但 与非线性有 限元有 关 的某 些 内容 ，例 如 ，刚度方程 、单元 刚度矩 阵 （简称 单元 刚阵 ）以及 主 、从 结 点换位 关系等方面 的研 究思路 与进 展 ，则结合工程应用 形式予 以提 纲性地讨论 刚 度 方 程 从 单元 到 结构 考虑 到 能量 泛 函极 值 时 ，就可 得 到 刚度方 程 这里 ，仅 以平面三 角元为例来说 明建立 刚度方程 的一般过程 任意平面三角元 ，如 图 所 示 任一结点的位移和结点力 分别为 图 平 面 三 角 元 单元结点位移和结点力分别为 第 2 页 单元 能量泛 函为 式 中 ， 称为弹性体 内任一点的应变列阵； 称为弹性体 内任一点的应力列阵 若单元位移函数为 则 式 中 ， 单元 的形 函数矩 阵 式 为位移 函数与单元 结点位移之 间的关系式 式 中， 单元应变矩 阵 式 （ 为应变与单元结点位移 之 间的关系式 式 中， 应 力 矩 阵； 弹性矩 阵 式 ）为应 力与单 元 结 点位 移关系 将 式 （ 代入 式 （ ，则有 根据 势能原理 ，将上 式对 取 变 分 ，可 得 由于 的任意性 ，故有 第 3 页 从而可 知单元 结 点位移 与结点力之 间的关系为 这 即是单元 的平衡方程 ，且 称为单元刚阵 将单元能量泛函式 （ ）对结构的所有离散单元予以叠加，则 可得到整体能量泛函： 将式 （ ）代 入 上 式 ，可 得 按 照单元 刚阵组装成整体 结构 刚阵 的 “对 号入 座 ”规律 以及 结 点力与结 点载荷在 结 点上 应 力平衡 的条件 ，可 以得 到 式 中， （ 为结点总数 ） 称为整体结构的结点位移； 称为整体结构的已知结点载荷 将 式 （ ）代 入 式 （ ，可 得 第 4 页 式 中 ， 式 （ ）即为 整 体 结构 能量 泛 函 的工程 实用 显 式 根 据 势 能原理 ，将 式 （ ） 对 取变分 ，可得 式 （ ）即是整 体 结构 的刚度 方程 整体 刚阵 具有三个特性 ，即对称 性 、正定性和稀疏性 ， 这里 不做详 细说 明，读 者可参 阅有 限元法 的基础 教材 值 得 指 出 的是： ① 刚度 方 程 式 （ ）尽 管是针 对平面 三 角元提 出来 的 ，但 其 思路 、作 法 与 形 式具有 普遍 意义 ，换 句 话说 ，无论对 哪一类 型