(精)1.2有限元法的工程应用——最新.pptVIP

1.2 有限元法的工程应用 一 有限元法的求解方法与步骤 有限元是一种以计算机为手段，通过离散化将研究对象变换成一个与原始结构近似的数学模型，再经过一系列规范化的步骤以求解应力、应变、位移等参数的数值计算方法。 什么是离散化呢？ 所谓离散化是将一个连续体分割成若干个通过节点相连的单元，这样一个有无限个自由度的结构就变换成一个具有有限个自由度的近似结构。该过程还包括对单元和节点进行编码以及局部坐标系和整体坐标系的确定。 什么是单元和节点呢？ 在有限元法中，将求解的工程结构看成是由许多小的、彼此用点联接的基本构件，如杆和梁、板和壳组成，这些基本构件称为单元，单元与单元之间的联接点称为节点。 1.2有限元法在结构工程中的应用 根据研究对象的不同，有限元法中采用的单元形式也不相同，常见的有以下5种： 1.桁架杆单元：主要应用于受轴向力作用的杆和杆系，如桁架结构； 2. 刚架杆单元：用于梁及刚架结构分析； 3.三角形平面单元：主要用于弹性力学中平面应力问题和平面应变问题的有限元分析； 4.三棱圆环单元：用于轴对称问题的有限元分析； 5.等参数单元：用于一些具有曲线轮廓的复杂结构。其特点是能简化复杂结构的单元划分工作，又能满足同样精度的要求时大大减少使用的单元数，成功地解决许多二维和三维弹性力学问题。 1.3有限元法求解问题的基本步骤 3.集成总体刚度矩阵[K]并写出总体平衡方程： 确定总体刚度矩阵的方法有三种： 确定总体刚度矩阵的方法有三种： 1.4 有限元法求解实例 解：1）离散化 2）求单元刚度矩阵 注意 3)总体刚度矩阵的集成和总体平衡方程的写出 （1）直接法求总体刚度矩阵[K] 因此，总体刚度矩阵为 （2）用集成法求总体刚度矩阵[K] 注意 该阶梯轴两种编码的对应关系为 集成[K]的步骤为： （2）将角标相同的系数相加，并按总码的顺序排列，则总体刚度矩阵为 （4）引入支撑条件，计算节点位移 （5）求单元中的应力及应变 例2 解：（1）结构的离散化 （2）各单元的刚度矩阵 （3）集成总体刚度矩阵 （4）由总体平衡方程求节点位移 练习 P189 1-3题 上式中的未知量 仍不能求出，因为[K]是一个奇异矩阵，必须引入支撑条件。在本例中支撑条件是节点1的位移为零，即 。这样总体平衡方程简化为 将已知数值代入上面公式解得 单元1中的应变： 单元2中的应变： 单元1中的应力： 单元2中的应力： 在光滑的水平面上有三个小车，他们彼此用四根弹簧相连，其连接方式如图1-2所示。小车1又通过弹簧k1与墙壁固连。每个小车上的作用力分别为F1、F2、F3。求每个小车的位移。 将每根弹簧都看成是一个单元，共有5个单元，将他们按弹簧的编号进行单元编号。三个小车处为三个节点，编码如图。从前面的例子可以看出，对固定节点若无需求出其约束力时，可不予以考虑。 单元1：一端为固定约束，另一端节点总码为1，则 单元2：两端节点总码分别为1和2，则 单元3：两端节点总码也分别为1和2，则 单元4：两端节点总码分别为1和3，则 单元5：两端节点总码分别为2和3，则 * * 1.1有限元的基本思想 1.结构离散化 对整个结构进行离散化，将其分割成若干个单元，单元间彼此通过节点相连； 2.求出各单元的刚度矩阵 是由单元节点位移量 求单元节点力向量 的转移矩阵，其关系式为 总体刚度矩阵[K]是由整体节点位移向量 求整体节点力向量 的转移矩阵，其关系式为 ，此即为总体平衡方程。 1.3有限元法求解问题的基本步骤 1）直接利用总体刚度系数的定义 在求出整体结构中各节点力与节点位移关系的基础上获得总体刚度矩阵。此方法旨在简单情况下才能采用。 2）集成法 将整体坐标下的单元刚度矩阵进行迭加而得。 这里所说的迭加不是简单的相加，而是将下角标相同的总体刚度系数相加，然后按总码的顺序对号入座。 3）利用节点间的刚度系数直接写出总体刚度矩阵 总体刚度矩阵对角线上的刚度系数 等于在节点i汇交的几个单元的刚度系数 之和；非对角线上的刚度系数 等于联结节点i与节点j间几个单元的刚度系数 之和。 4.引入支撑条件，求出各节点的位移 节点的支撑条件有两种：一种是节点n沿某个方向的位移为零，另一种是节点n沿某个方向的位移为一给定值。 5.求出各单元内的应力和应变。 1.3有限元法求解问题的基本步骤 例1 一根由两段组成的阶梯轴，一端固定，另一端承受一个轴向载荷F3。这两段的横截面积分别为 和 ，长度分别为 和 ，弹性模量分别为 和 ，求