第4章 平面问题有限单元法分析.ppt

[例]有一正方形板，沿对角承受压力作用，板厚 ，载荷 ，为简化计算，设泊松 比 ，材料弹性模量为 ，求它的应力分布。 解：由于结构几何对称，所受载荷也对称，所以可取板的1/4做为计算对象。单元分划如图。 因为对称性，节点1，2，4没有水平位移；节点4，5，6没有垂直位移。 单元（1），（2），（4）为一种形式，单元（3）为另一种形式。 模型为平面应力问题。 由于 ，所以 单元（1）的节点坐标为 面积为 代入单元刚度矩阵的一般表达式可得 或者 同样可以得到单元（2），（4）得刚度矩阵，结果与单元（1）一样。 单元（3）的节点号为2，5，3，节点坐标为 单元面积仍为 。刚度矩阵为 或者 单元刚度矩阵中的子刚度矩阵 ，角标 表示节点号，即实际在总体刚度矩阵中的位置，叠加总体刚度矩阵时一定要严格按节点号叠加。 例如， 中的4个元素，要写在第3、第4行和第9、第10列。 按照节点号，将单元刚度矩阵叠加成总体刚度矩阵 引入约束条件 节点载荷 则位移和力向量为 从而线性代数方程组为 划去位移为零的行和列，方程组成为 解此方程组得未知的位移分量 计算单