有限元法课件.pptVIP

对于平面应力问题: 工程有限单元法 2.2.4 单元刚度矩阵 单元节点力与单元位移的关系式，称为单元刚度方程组。 工程有限单元法 单元刚度矩阵的性质： （1）单元刚度矩阵中每个元素有明确的物理意义； （2）刚度矩阵是对称矩阵； （3）刚度矩阵是奇异矩阵； 另外，单元刚度矩阵取决于： （1）单元的位移函数； （2）单元的几何参数； （3）单元的材料性质。 工程有限单元法 2.2.5 单元等效节点载荷 连续弹性体离散为单元组合体时，为简化受力情况，需把弹性体承受的任意分布的载荷都向节点移置(分解)，而成为结点载荷。如果弹性体受承受的载荷全都是集中力，则将所有集中力的作用点取为节点，就不存在移置的问题，集中力就是节点载荷。但实际问题往往受有分布的面力和体力，都不可能只作用在节点上。因此，必须进行载荷移置。如果集中力的作用点未被取为节点，该集中力也要向结节移置。 将载荷移置到节点上，必须遵循静力等效的原则。静力等效是指原载荷与节点载荷在任意虚位移上做的虚功相等。在一定的位移模式下，移置结果是唯一的，且总能符合静力等效原则。 工程有限单元法 在线性位移模式下，对于常见的一些载荷，可以通过简单的虚功计算，得出所需的载荷列矩阵。 均质等厚度的三角形单元所受的重力，把1/3的重力移到每个节点 工程有限单元法 例： 总载荷的2/3移置到节点i，1/3移置到节点j，与原载荷同向 工程有限单元法 载荷向节点的移置，可以用普遍公式来表示。 体力的移置 分布面力的移置 在线性位移模式下，用直接计算法简单；非线性模式下，要用普遍公式计算。 工程有限单元法 2.2.6 总刚度矩阵 [K]为总刚度矩阵，{R}为节点力分量矩阵。 为节点位移分量矩阵。 总刚度矩阵性质： （1）总刚度矩阵也是对称矩阵； （2）总刚度矩阵呈稀疏带状分布； （3）总刚度矩阵奇异矩阵。 工程有限单元法 2.2.7 边界约束条件 有限元法中通常采用两种方法: 划行划行法和乘大数法. 其中前者适用于简单的手算练习,后者适合于实际问题的计算机处理. 工程有限单元法 2.2.8 解题步骤与算例 有限元法的一般分析步骤如下: (1)首先绘出结构的几何简图,在此基础上将结构离散; (2)其次进行单元分析; (3)组集总刚度矩阵; (4)最终求单元应力和节点应力. 工程有限单元法 算例讲解 P27 工程有限单元法 2.2.9 计算结果处理 有限元中计算结果主要包括位移和应力两方面. 其中位移可根据计算结果中的节点位移分量画出结构的位移图. 而对于应力计算结果必须进行整理.方法有: (1) 形心法; (2) 绕节点法; (3) 二单元法. 工程有限单元法 2.2.10 平面高阶单元 为了提高有限元法计算结果的精度，除了增加单元数目外，还常采用具有较高次位移函数的单元，即高阶单元。 常用的四节点矩形单元和六节点三角形单元。 工程有限单元法 1、四节点任意四边形等参数单元 任意四结点四边形单元 四结点正方形单元 工程有限单元法 1、八节点任意四边形等参数单元 四边形八结点单元 八结点基本单元 工程有限单元法 3、应用等参单元应注意以下几点问题 1）各向长度的相对大小：单元长度之比不宜相差太大，接近正方形的单元误差最小，长宽比很大，误差也很大。 2）棱边的曲折：应使单元边上没有折点，如边上不可避免有折点，应使棱边只有凸出的折点。 3）棱边的夹角：尽量接近90度。 4）棱边上节点的间距：尽量均匀。 工程有限单元法 对于简单的线性、边界规则问题，可采用常规的数学方法进行解析求解，得到问题的精确解。 而对于非线性、边界不规则的复杂问题，在满足精度预期条件下，可采用整体离散化、局部连续化方法建立有限元模型，通过数值解得到原问题的近似解。 * * * * * * * * * * * * 平面应力问题 应力矩阵(1-2) 可以简化为： 工程有限单元法 物理方程(1-10)中后两式可见，这时的剪应变： 由物理方程(1-10)中的第三式可见： 一般 ， 并不一定等于零，但可由 及 求得，在分析问题时不必考虑。于是只需要考虑 三个应变分量即可，于是应变矩阵(1-3-2)简化为： 工程有限单元法 平面应力问题 物理方程(1-10)简化为： 转化成应力分量用应变分量表示的形式： 工程有限单元法 平面应力问题 将(1-21)式用矩阵方程表示： 它仍然可以简写为： 弹性矩阵[D]则简化为： 工程有限单元法 平面应力问题 只有