弹性力学平面问题的有限单元法02综述.ppt

讲 授：陈得良 TelQQ：416501065 Email：deliang_chen@ ;*; 严格地说，实际的弹性结构都是空间结构，并处于空间受力状态，属于空间问题，然而，对于某些特定问题，根据其结构和外力特点可以简化为平面问题来处理。这种近似，可大大减少计算工作工作量，为有限元分析提供方便。弹性力学平面问题可分为两类： 如图柱形管道和长柱形坝体，具有如下特点：a纵向尺寸远大于横向尺寸，且各横截面尺寸都相同；b 载荷和约束沿纵向不变，因此可以认为，沿纵向的位移分量 等于零。;; 1.平衡方程;若令; 在边界;设变形体处于平衡受力状态：体积力为;其中;分析一个典型三角形单元的力学特性。首先建立以单元结点位移表示单元内各点位移的关系式。设单元的结点号码 ，其结点的坐标分别 ; 由于单元体也是一个二维的弹性体，单元内各点的位移分量是坐标;设结点;其中;将(d)式代入(b)式中的第一式，并稍加整理得;如令;*;有了单元的位移模式，就可以利用平面问题的几何方程;求得应变分量，将(e)、(d)两式代入上式即得 ;其中; 便可导出以结点位移表示应力的关系式。把(2-7)式代入上式，得到;对于平面应力问题，;对于平面应变问题;如果注意到(4-1)式，则(4-11)式可写成;上节讨论常应变三角形单元时，曾提出形函数;由(4-3)式可知，常数;类似地有（参见《结构及弹性力学有限元》，刘怀恒主编;2 在单元任一点上三个形函数之和等于1;3 在三角形单元; 利用这一性质，很容易证明相邻两个单元的位移，分别进行线性插值之后 ，在公共边是连续的，例如图中单元 和 具有公共边;;1．试证：在三结点三角形单元内的任意一点都有; ;*;*;*;在3结点等厚三角形单元中，如果单元的材料是均质的，矩阵 中的元素是常量，而且在三角形常应变单元情况下，矩阵 中的元素也是常量，当单元的厚度 t 也是常量时，则有 ，; 物理意义：单元刚度矩阵 中的任一列 的元素分别等于该单元的某个结点 沿坐标 j方向发生单位位移时，在各结点上所引起的结点力，它决定于该单元的形状、大小、方位和弹性常数，而与单元的位置无关。即不随单元或坐标轴的平行移动而改变。 将表达式(4-8)代入(4-22)式，即该平面应力问题中三角形单元的刚度矩阵，写成分块形式如下;其中; 用有限元解题，需把弹性体离散化后的每个单元所受的体力、面力、集中力都移置到有关的受载结点上，形成单元等效结点载荷列阵 ，这种移置时按静力等效原理进行的。;*;1 集中力 的移置：; ;2 体积力 的移置。;3 面力 的移置; 关于刚体静力等效原理这里不在赘述；现将四种特殊载荷简化结果说明如下：; 3. 线性分布载荷，在结点 处强度为 ，在结点 处强度为0，则合力大小为 只需将合力的 移置到结点 ， 移置到 结点。 ;*; 假说弹性体被划分为 个单元和 个结点，并对每一个单元都进行上述运算，则得到 组形如 式的方程，把这些方程集合起来，便可得到表征整个弹性体平衡的表达式。为此目的，我们首先引进整个弹性体的结点位移列阵 ，它是由各结点位移按给定的号码从小到大顺序排列组成的，即; 再引进整个弹性体的载荷列阵 ，它是移置到结点上的等效结点载荷，按照点的号码从小到大顺序排列组成的，即;将各单元的结点力列阵 加以扩大，使之成为 阶列阵; （4-24）式中，矩阵号上面的 表示在分块矩阵意义下 所在 的列的位置。这里以假定的次序恰和结点号码的次序从小到大的排列是一致的。各单元的结点力列阵经过这样扩大以后便可相加，将全部单元的结点力列阵叠加在一起，就得到下式表示的弹性体的载荷列阵，即;将(4-24)式确定的六阶方阵 加以扩大，使成为 阶的方阵;(4-32)式中的2×2阶子矩阵 被放到(4-35)式中的第 双行、第 双列中。; 上式左边 是弹性体所有单元刚度矩阵的总