11平面问题的有限元分析及三角形单元的应用.pptVIP

第十一讲 平面问题的有限元分析及三角形单元的应用  如图8-7所示单元的边界ij上受到沿x方向作用的按三角形分布的荷载，在i点的荷载集度为q,则其等效结点荷载列阵为 （8-42） 如上述两种荷载为任意方向，可分别按其在x方向和y方向的分量进行计算。 如单元上同时作用有集中力，分布体力和分布面力，则结点荷载应为三部分荷载之和，即： （8-43） 下举例说明用单元刚度矩阵形成整体刚度矩阵的方法。图8-8表示一弹性体划分为三个单元，并对结点进行局部编号，整体编号为1，2，3，4，55单元结点局部编号为I,j,m按逆时针方向排序。整体刚度矩阵为10阶方阵。 首先将单元e的刚度矩阵 扩大为10阶方阵，成为单位贡献矩阵。单元局部编号应与其所在位置的整体编号相对应。以单元3为例，其结点局部编号I,j,m与整体编号2，4，5对应，形成单元3的贡献矩阵。 整体号 1 2 3 4 5 i j m 局部号 然后将各个单元贡献矩阵叠加起来，形成整体刚度矩阵。 在实际编程中，为节省容量，将各个单元刚度矩阵 中的子矩阵 逐个搬到整体刚度矩阵对应的位置上去，建立整体矩阵如下： 需要指出说明的是，离散弹性体的整体的结点荷载列阵R也是由各单元的结点荷载加以集成得出的。 在形成整体刚度矩阵之后，利用式（8-44）还不能直接求解结点位移，尚需引入唯一边界条件，以保证结构具有确定的位移状态，求解未知的结点位移。 * * 有限元法及程序设计 第一节 概述 第二节 单元分析 第三节 等效结点荷载 第四节 整体刚度矩阵 第五节 平面问题分析举例 第六节 单元网格的划分和计算成果的整理 下面给出作用于单元边界上荷载的两种简单情况。如图8-6所示的边界ij上沿x方向作用一集中力 ，其作用点距i，j两点的距离分别为 和 。其等效结点和在列阵为 （8-41） 用集成刚度法将e单元的单元贡献矩阵表示如下, （8-45） 整体号1 2 3 4 5 体内任一点应力分量可表示为 （8-3） 任一点的应变分量可表示为 （8-4） 任一点的位移分量可表示为 （8-5） 弹性力学平面问题的几何方程的矩阵表达式为 （8-6） 平面应力问题的物理方程的矩阵表达式为 （8-7） 或简写成为 （8-8） 式中 （8-9） 称为弹性矩阵。 平面应变问题的物理方程也可写成式（8-8），但须将式（8-9）中的E换成，μ换成，因此得出 （8-10） 平衡微分方程及边界条件也可以用矩阵表示，但弹性力学有限元位移法中，通常用虚功方程代替平衡微分方程和应力边界条件。虚功方程的矩阵表达式为 （8-11） 式中：，表示虚位移； ，表示与虚位移相对应的虚应变。 为了便于计算，作用于弹性体上的体力和面力替换为作用在结点上的集中力，即等效结点荷载。设作用于各个结点上的外力分量用如下列阵来表示 与这些结点外力分量相对应得结点虚位移分量列阵为 则外力在虚位移上做的虚功为 如平面弹性体的厚度为t，该虚功除以t，即可得出单位厚度薄板上的外力虚功。于是，式（8-11）所示虚功方程可写成 （8-11） 虚功方程不仅仅应用于弹性力学，也可用于塑性力学。其应用条件是：只要变形体的全部外力和应力满足平衡方程；位移是微小的，