有限元法基本原理.ppt

（ 1.13 ）和（ 1.14 ）用矩阵表示： 即： — 为单元 的节点力向量 ； — 为单元 的节点位移向量。 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? j i jj ji ij ii j i j i u u K K K K u u l EA U U 1 1 1 1 ? ? 15 . 1 ? ? ? ? ? ? e e e K F ? ? ? ? 16 . 1 e ? ? ? ? ? ? ? ? ? j i e u u ? e — 为单元 的刚度矩阵。 刚度系数：单位节点位移分量所引起的节点力分量。 单元：节点力为 ，内力为 ，节点虚位移 ， 单元虚位移 ，单元虚应变 。 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1 1 1 1 l EA K K K K K jj ji ij ii e ? ? 17 . 1 e ? ? e F ? ? ? ? ? e * ? ? ? e f * ? ? e * ? 单元的外力虚功： 单元的虚应变能： 虚功原理 ： ? ? ? ? ? ? e T e F V * ? ? ? ? ? ? ? ? ? ??? ? V T e dV U ? ? ? * U V ? ? ? 将 代入上式，得： ? ? ? ? ? ? ? ? ? ?? ? ? ? ? ? ? ? ? ? ? e e B D B ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?? ? ? ? ??? ? e T T e e T e dV B D B F ? ? ? * * ? ? ? ? ? ? ? ? ? ?? ? ? ? ? ? 0 * ? ? ??? e T e T e dV B D B F ? ? 由虚位移的任意性： 与（ 1.16 ）式对比，得： 称为单元的刚度矩阵，简称单刚。 ? ? ? ? ? ?? ? ? ? 0 ? ? ??? e T e dV B D B F ? ? ? ? ? ? ?? ? ??? ? dV B D B K T e § 1.6 整体刚度矩阵的集成 整体结构 → 离散 → 有限个单元 → 单元分析 → 单元特性集合 台阶轴： 台阶轴 离散化模型 单元模型 1 1 2 2 2 3 3 ① ① ② ② 1 P 3 P ? ? a ? ? b ? ? c 1 1 u P 1 1 u P 3 3 u P 3 3 u P 2 u 2 2 u P 2 2 u P ? ? a ? ? b ? ? c 则单元（ 1 ）的刚度方程，考虑（ 1.15 ）式： 单元（ 2 ）的刚度方程： ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1 2 1 1 1 22 1 21 1 12 1 11 1 2 1 1 u u K K K K P P ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1 2 1 22 1 1 1 21 1 2 1 2 1 12 1 1 1 11 1 1 u K u K P u K u K P ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 2 3 2 2 2 33 2 32 2 23 2 22 2 3 2 2 u u K K K K P P 台阶轴在外力 和 的作用下发生变形时，在节点处 的变形必须是连续的，即： 。 三个节点的位移： ， ， ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 2 3 2 33 2 2 2 32 2 3 2 3 2 23 2 2 2 22 2 2 u K u K P u K u K P ? ? ? ? 1 P 3 P ? ? ? ? 2 2 1 2 u u ? ? ? 1 1 1 u u ? ? ? ? ? 2 2 1 2 2 u u u ? ? ? ? 2 3 3 u u ? 台阶轴处于平衡状态，各个节点的力是平衡的。 ? 节点 1 ：外力 ，节点力 ； ? 节点 2 ：外力为 0 ，节点力 ； ? 节点 3 ：外力 ，节点力 。 1 P 3 P ? ? ? ? 0 2 2 1 2 ? ? P P ? ? 1 1 P ? ? 2 3 P 农业机械有限元软件方法 吉林大学生物与农业工程学院 韩志武 第一章 有限元法基本原理 § 1.1 有限元法方法及其历史 § 1.1 有限元法简介及其历史 ? 精确解：少数方程性质简单，形状规则； ? 复杂问题解：简化假设 + 数值