第三讲平面问题有限元分析.pptVIP

2-7 单元刚度矩阵的物理意义及其性质 单元刚度矩阵的性质： 1)对称性： 是对称矩阵 2)奇异性： 是奇异矩阵， 单元刚度矩阵所有奇数行的对应元素之和为零，所有偶数行的对应元素之和也为零。由此可见，单元刚度矩阵各列元素的总和为零。由对称性可知，各行元素的总和也为零。 2-3 单元位移函数 选择单元位移函数时，应当保证有限元法解答的收敛性，即当网格逐渐加密时，有限元法的解答应当收敛于问题的正确解答。因此，选用的位移模式应当满足下列两方面的条件： (1) 必须能反映单元的刚体位移和常量应变。 6个参数 到 反映了三个刚体位移和三个常量应变。 (2) 必须保证相邻单元在公共边界处的位移连续性。 (线性函数的特性) 形态函数Ni具有以下性质： 1) 在单元结点上形态函数的值为1或为0。 2）在单元中的任意一点上，三个形态函数之和等于1。 用 来计算三角形面积时，要注意单元结点的排列 顺序，当三个结点i，j，m取逆时针顺序时， 当三个结点i，j，m取顺时针顺序时， 图 4-3 2-3 单元位移函数 作业：图示等腰三角形单元，求其形态矩阵[N],应变矩阵,应力矩阵,单元刚度矩阵。 2-3 单元位移函数 由三角形的面积 2-4 单元载荷移置 连续弹性体离散为单元组合体时，为简化受力情况，需把弹性体承受的任意分布的载荷都向结点移置(分解)，而成为结点载荷。如果弹性体受承受的载荷全都是集中力，则将所有集中力的作用点取为结点，就不存在移置的问题，集中力就是结点载荷。但实际问题往往受有分布的面力和体力，都不可能只作用在结点上。因此，必须进行载荷移置。如果集中力的作用点未被取为结点，该集中力也要向结点移置。 将载荷移置到结点上，必须遵循静力等效的原则。静力等效是指原载荷与结点载荷在任意虚位移上做的虚功相等。在一定的位移模式下，移置结果是唯一的，且总能符合静力等效原则。 单元的虚位移可以用结点的虚位移 表示为 表示为， （3-15） 令结点载荷为 令结点载荷为 集中力的移置 如图所示，在单元内任意一点作用集中力 由虚功相等可得， 由于虚位移是任意的，则 例题1：在均质、等厚的三角形单元ijm的任意一点p（xp，yp） 上作用有集中载荷。 体力的移置 令单元所受的均匀分布体力为 由虚功相等可得， 分布面力的移置 设在单元的边上分布有面力 ，同样可以得到结点载荷， 例题：设有均质、等厚的三角形单元ijm，受到沿y方向的 重力载荷qy的作用。求均布体力移置到各结点的载荷。 同理， 例题：在均质、等厚的三角形单元ijm的ij边上作用有 沿x方向按三角形分布的载荷，求移置后的结点载荷。 取局部坐标s，在i点s=0，在j点s=l，L为ij边的长度。 在ij边上，以局部坐标表示的插值函数为， ， ， 载荷为 例： 总载荷的2/3移置到结点i，1/3移置到结点j，与原载荷同向 2-5 单元应力矩阵 本节利用几何方程、物理方程，实现用结点位移表示单元的应变和单元的应力。 用结点位移表示单元的应变的表达式为 ，[B]矩阵称为几何矩阵。 2-5 单元应力矩阵 由物理方程，可以得到单元的应力表达式 为应力矩阵 2-6 单元刚度矩阵 讨论单元内部的应力与单元的结点力的关系，导出用结点位移表示结点力的表达式。 由应力推算结点力，需要利用平衡方程。第一章中已经用虚功方程表示出平衡方程。 2-6 单元刚度矩阵 考虑上图三角形单元的实际受力，结点力和内部应力为： 任意虚设位移，结点位移与内部应变为 2-6 单元刚度矩阵 令实际受力状态在虚设位移上作虚功，外力虚功为 2-6 单元刚度矩阵 计算内力虚功时，从弹性体中截取微小矩形，边长为dx和dy，厚度为t，图示微小矩形的实际应力和虚设变形。 2-6 单元刚度矩阵 微小矩形的内力虚功为 整个弹性体的内力虚功为 2-6 单元刚度矩阵 根据虚功原理，得 这就是弹性平面问题的虚功方程，实质是外力与应力之间的平衡方程。 虚应变可以由结点虚位移求出： 代入虚功方程 2-6 单元刚度矩阵 接上式，将应力用结点位移表示出 有 令