有限元方法概论01绪论.ppt

有限元方法概论 第一章 绪 论 1.1 课程介绍 科学研究三大方法： 什么是有限元法（FEM，Finite Element Method）？ 1.2有限元的发展概况 1.3 有限元方法发展现状 由单一结构场求解发展到耦合场问题的求解 1.3 有限元方法发展现状 1.3 有限元方法发展现状 1.3 有限元方法发展现状 1.3 有限元方法发展现状 作用：实现产品效益最大化目标（降低开发成本、缩短研制周期） 举例1：对概念设计的应力分析、方案评价 举例2：瞬态应力分析 举例3：整车碰撞分析 举例4：安全气囊展开过程分析 举例5：汽车的流体动力学分析 举例6：发动机缸体热应力分析 举例7：汽车覆盖件冲压成型分析 举例8：复杂钣金件结构拓扑优化 有限元法的基本思路： 有限元法是将连续体理想化为有限个单元集合而成，单元之间仅在有限个节点上相连接，亦即用有限个单元的集合来代替原来具有无限个自由度的连续体。 1.5 有限元法的基本思路及基本分析过程 1.5 有限元法的基本思路基本分析过程 思考题 1.5 有限元法的基本思路基本分析过程 思考题 1. 有限元法的基本思想？ 2. 有限元分析的基本步骤？ 3. 你希望将来从事什么工作？会不会用到有限元理论？ 设有一根受自重作用的等截面直杆，上端 固定，下端自由，单位杆长的重力为q，杆长为L，横截面积为A，材料弹性模量为E。试求该直杆各横截面上的应力。其材料力学的精确解答： 解得 下面我们用虚位移原理把每个单元的重力等效地移置到节点上去。 由于?i*＝1，?j* ＝0，所以得单元虚位移为 于是，单元重力所做虚功为 Pie所做虚功为1? Pie ，Pje而不做功。根据虚功原理，移置前单元荷载所做虚功应该等于移置后节点荷载所做虚功，即有 Pie＝(xj-xi) ?q/2 同理，可求得 Pje＝(xj-xi) ?q/2 若记 则单元荷载列矩阵为 由此可见，只要把单元重力q(xj-xi)平均移置到节点上即可。 对于本例，有 P11 = P21 = P22 = P32 = P33= P43＝ qL/6 再考虑处于固定端的节点1还受有约束反力 P=-qL，于是，单元荷载移置后，各节点的节点荷载分别为 各节点的节点荷载分别为 第二步：分别用几何方程、物理方程以及虚功方程分析单元的应变、应力以及节点受力与节点位移的关系。 其中 1.7 一个简单的例子 一维弹性问题的物理方程为 ? ＝E? 可写为 {?}= [D]{?} 其中，[D]=E 为弹性矩阵。 或 其中 [S] = [D][B] = [Si Sj] = [[D][Bi], [D][Bj]] 即 称为应力矩阵，它反映了单元应力与节点位移之间的关 系。 下面利用虚功原理分析单元的节点受力与节点位移之间的关系。 直杆分割成3个单元之后，相邻单元之间的作用力通过节点来传递。我们把节点对单元和单元对节点的作用力都称为节点力，单元e的节点i、j处的节点力分别用Rie、Rje表示（右图），并规定节点对单元的节点力以沿x轴正向为正。对单元来讲，节点力为外力，若记 则外力所做的虚功为 而内力所做虚功为 根据虚功原理，单元虚功方程为 由于 而单元虚应变与节点位移也满足单元应变与节点位移间的关系式，即 所以，由单元虚功方程可得 单元虚功方程 其中[K]e为单元刚度矩阵，它反映了单元的节点力与节点位移之间的关系。 根据应变矩阵[B]和应力矩阵[S]，可计算单元刚度矩阵[K]e ，由单元刚度矩阵表达式，可得 其中 第三步：由节点平衡方程建立并求解以节点位移为未知量的线性代数方程组。 根据各节点的受力情况，可得 引入 {P}=[P1 P2 P3 P4]T 称为总体荷载矩阵。 要使单元节点力与节点荷载相等（静力平衡），须将单元刚度矩阵升阶： 1.7 一个简单的例子 将它们代入静力平衡式中，使得到以节点位移为未知量的线性方程组，即 若记 其中 [K]——总体刚度矩阵，它反映总体节点位移与总体 荷载之间的关系 {δ}——总体节点位移列矩阵 {P} ——总体节点荷载列矩阵 根据各单元刚度矩阵可计算总体刚度矩阵。由于各单元的