现代设计方法基础(有限元法基础).pptVIP

1943年，Courant提出有限元法概念 1956年，Turner和Clough第一次用三角形单元离散飞机机翼，借助有限元法概念研究机翼的强度及刚度 1960年，Clough正式提出有限元法（FEM） 20世纪60年代以后，由于数学界的参与，FEM得到蓬勃发展，并且扩大了应用 有限元分析是一种工程物理问题的数值分析方法，根据近似分割和能量极值原理，把求解区域离散为有限个单元的组合，研究每个单元的特性，组装各单元，通过变分原理，把问题化成线性代数方程组求解。 化整为零，裁弯取直，变难为易，先拆后搭 杆系单元 定义 杆系结构中的杆件、梁、柱等称为杆系单元。连接的点称为节点。 杆系单元为一维单元。 结构离散 一般原则： 杆系的交叉点、边界点、集中力作用点、杆件截面尺寸突变处等都应该设置节点，节点之间的杆件即构成单元。 2 .求单元刚度矩阵 注意 集成[K]的步骤为： （2）将角标相同的系数相加，并按总码的顺序排列，则总体刚度矩阵为 解：（1）结构的离散化 （2）各单元的刚度矩阵 集成总体刚度矩阵 （4）由总体平衡方程求节点位移 式中 第二章 有限元法基础 3. 整体分析 ? 把各个单元的刚度矩阵集成为总体刚度矩阵，以及将各单元的节点力向量集成总的力向量，求得整体平衡方程。 集成总体刚度矩阵[K]并写出总体平衡方程： [K]是由整体节点位移向量 求整体节点力向量 的转移矩阵，其关系式为 ，这就是总体 平衡方程； 第二章 有限元法基础 确定总体刚度矩阵[K]的办法 1）直接利用总体刚度系数的定义 在求出整体结构中各节点力与节点位移关系的基础上获得总体刚度矩阵。此方法只在简单情况下才能采用。 2）集成法 将整体坐标下的单元刚度矩阵进行迭加而得。这里所说的迭加不是简单的相加，而是将下角标相同的刚度系数相加，然后按总码的顺序对号入座。 第二章 有限元法基础 3）利用节点间的刚度系数直接写出总体刚度矩阵 总体刚度矩阵对角线上的刚度系数 等于在节点i汇交的几个单元的刚度系数 之和；非对角线上的刚度系数 等于联结节点i与节点j间几个单元的刚度系数 之和。 第二章 有限元法基础 4. 确定约束条件 　　由上述所形成的整体平衡方程是一组线性代数方程，在求解之前，必修根据具体情况分析，确定求解对象问题的边界约束条件，并对这些方程进行适当修正。 第二章 有限元法基础 5. 有限元方程求解 　　通过求解整体平衡方程，即可求得各节点的位移进而根据位移可计算单元的应力及应变。 6. 结果分析与讨论 第二章 有限元法基础 网架杆件 节点位移 单元刚度矩阵 总刚度矩阵 总刚度方程 节点位移值 杆件内力 单元内力与节点位移间关系 引入边界条件 节点平衡及变形协调条件 基本单元 基本未知量 第二章 有限元法基础 有限元求解实例分析 【例1】一根由两段组成的阶梯轴，一端固定，另一端承受一个轴向载荷F3。这两段的横截面积分别为A(1)和A(2)，长度分别为L(1)和L(2)，弹性模量分别为E(1)和E(2) ，求出这两段的应力和应变。已知数据分别为F3=100N， 第二章 有限元法基础 1 ① 2 ② 3 A(1) E(1) A(2) E(2) L(1) L(2) ① ② 2 Φ1 F2 F3 Φ2 Φ3 F1 F3 第二章 有限元法基础 【解】 1.离散化 把这根阶梯轴看成是由两个单元组成的，节点选在截面积突变处，两个单元的连接处是一个节点，该阶梯轴的两端视为另外两个节点，所以整个结构共有三个节点。这根轴是一维结构，并只受轴向载荷，因此各单元内只有轴向位移。三个节点位置的位移量分别记为 、 、 。在整个结构中节点载荷及节点位移均用大写字母标记，其角标为节点在总体结构中的编码，简称总码。 第二章 有限元法基础 下面分析某等截面单元（e）。当两端分别承受两个轴向力 和 作用时的位移情况。根据材料力学的知识可知，在两端节点i、j处的位移量 和 与轴向力 和 的关系式为 第二章 有限元法基础 在分析单元刚度矩阵时，载荷F和位移 等参数的上角标为该单元的编码，下角标为该单元内节点的局部编码。上两式可写成： 或简写为： 第二章 有限元法基础 式中 —为单元刚度矩阵或单元特性矩阵，其阶数等于单