第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列。 按照节点号，将单元刚度矩阵叠加成总体刚度矩阵 引入约束条件 节点载荷 则位移和力向量为 从而线性代数方程组为 划去位移为零的行和列，方程组成为 解此方程组得未知的位移分量 计算单元应力 由于只有6个非零位移分量，所以单元（1）（2）和（4）的应力矩阵为 对于单元（3）有 各单元应力为 从计算结果可以看出，单元内的应力为常数，相邻单元的应力值会产生突变。只有在单元尺寸很小时，应力的突变值才会减小。这是三角形单元的主要缺点。 课后作业 1. 如图（a）所示一高深悬臂梁，在右端部受集中力F作用，材料弹性模量E、泊松比v＝1/3，悬臂梁的厚度（板厚）为t，如图（b）所示有限元模型，试按平面应力问题，求结构总体刚度矩阵和总体载荷向量，引入边界条件，求解节点1，2的位移。 2. 下图为一固定端梁受集中力P作用，试用图所示单元划分，求出节点位移（按平面应力问题计算，取E为常量，μ=0，t=1)。 3. 三角形单元i j m的j m边作用有右图所示线性分布面载荷，求节点载荷向量。 * * * * * * * * * * * * 如以无量纲坐标系来表达，则上式可以写成 其中： 单元应变场 　　根据单元的位移场函数式，由几何方程可以得到单元的应变场表达式， 记为： 　　这里，[B] 矩阵称为几何矩阵。 [B] 矩阵可以表示为分块矩阵的形式 其中 单元应力场 　　由物理方程及应变矩阵，可以得到单元的应力场表达式， 其中 　 为应力矩阵，[D]称为弹性矩阵，对于平面应力问题， 将应力矩阵表示为分块矩阵的形式 其中： 对于平面应变问题，只需将 E 换为　　　，　 换为　　　。 单元刚度矩阵 　　和三角形单元一样，可以根据虚功理导出节点位移向量和节点力向量之间关系，即单元的刚度矩阵　　，可以将其写成分块的形式。 其中 　　对于平面应力问题，如果单元厚度 t 为常数，则可得刚度矩阵的显式形式： 例题4 图（a）所示的深梁，在跨中受集中力F的作用，若取 试用有限单元法求解跨中的位移。 解 1.将图（a）划分网格，化为离散化结构，如图（b）所示。由于结构具有对称性，可取 部分进行分析，如 图（c） 所示。 2. 在图（c）中，只有两个未知结点位移 其余的结点位移均为零。 未知的结点位移列阵是 对应的结点荷载列阵是 3.下面我们直接来建立对应于未知结点 位移的平衡方程式， 4.对于三角形单元，按照结点的局部编号 结点力一般公式是 当 且结点的局部编号如图（d)、(e) 时， 对于I 单元，结点的局部编号与整体编号的关系是 则有 其中 由上式，得出 I单元中 不存在，而 对于II 单元，结点的局部编号与整体编号的关系是 。 得 其中 由上式，可得II单