武科大现代设计方法4–2杆有限元分析.ppt

§4.2 弹性力学平面问题的有限元法 有 限 元 求 解 基 本 原 理 §4.2 弹性力学平面问题的有限元法 弹性力学的有限元分析计算可分为三个步骤: 1 结构离散化 这是有限元法的基础，用由有限个方位不同但几何性质及物理性质均相似的单元组成的集合体来代替原来的连续体或结构。每个单元仅在节点处和其他单元及外部有联系。对于不同的问题，根据自身的特点，可选用不同类型的单元。对同一问题也可以分别或同时选用多种单元。 单元：即原始结构离散后，满足一定几何特性和物理特性的最小结构域 节点：单元与单元间的连接点。 节点力：单元与单元间通过节点的相互作用力。 节点载荷：作用于节点上的外载。 节点自由度(DOFs) ：用于描述一个物理场的响应特性。 分离但节点重叠的单元A和B之间没有信息传递（需进行节点合并处理） 非法结构离散 2.单元分析 主要内容：由节点位移求内部任意点的位移，由节点位移求单元应变，应力和节点力。 3.整体分析 (1) 由节点平衡方程，建立以整体刚度矩阵[K]为系数的，整体节点位移{d }和外载{R}的关系式——整体平衡方程。 (2) 考虑几何边界条件，修改总体刚度矩阵，求解全部未知位移分量。 二. 有 限 元 求 解 基 本 原 理(一维问题) 引例:用有限元法求图1所示受拉阶梯杆的位移和应力。已知杆截面面积A(1)=2×10-4m2，A(2)=1×10-4m2，各段杆长l(1)=l(2)=0.1m；材料弹性模量E(1)=E(2)=2×107Pa，作用于杆端的拉力F3=10N。 1.单元划分 根据材料力学的平面假设，等截面受拉杆的同一截面可认为具有相同的位移和应力，即位移只与截面的轴向坐标(x) 有关，所以可将阶梯杆看作由两个“一维单元”组成，同一个单元内截面面积及材料特性不变。最简单的情况是，每一个单元有两个节点，他们分别位于单元两端。相邻两单元靠公共节点联结。受拉阶梯杆就简化为由两个一维单元和三个节点构成的有限单元模型。图中①和②是单元号，1,2,3是节点号。取节点位移作为基本未知量，应力由求得的节点位移算出。 2.确定单元插值函数(形函数) 有限元法将整个求解域离散为一系列仅靠公共节点联结的单元，而每一个单元本身却视为光滑连续体。单元内任一点的场变量(如位移)可由本单元的节点值根据场变量在单元中的假定分布规律(插值函数)插值求得。 本例中，每单元有两个节点，采用线性插值。图c是一典型单元图，两节点分别为i和j，节点场变量值分别记为ui和uj 。设单元中坐标为x处的场变量为u(x) 。 3.单元方程(单元节点位移与节点力的关系) 由等截面杆变形与拉力的关系(虎克定律)得到 (3) 式中， Pi和Pj分别为作用于单元e的节点i和节点j的节点力。 式(3)写成矩阵形式为： (4) 或简记为： [k]e qe = Pe (5) [k]e常称为单元刚度矩阵(stiffness matrix of element)，简称单元刚阵: P e={Pi Pj}T 称为单元节点力列阵(nodal force vector)。 式(5)称为单元方程。 到目前为止，单元方程(4)或(5)尚不能求解，因为节点力列阵Pe尚属未知。 Pe的分量Pi和Pj为相邻单元作用于单元e的节点i和j的力，即属于单元之间的作用力。只有将具有公共节点的单元“组 集”在一起才能确定上述节点力和节点外载荷之间的关系。 4.单元组集 建立总体方程组为获得总体方程组，必须先将单元方程按照局部自由度(ui和uj)和总体自由度(u1、u2和u3)的对应关系进行扩展。 具体来说，单元1的扩展方程为： (6) 式中，各项上角码表示单元序号；下角码表示自由 度总体序号。 单元2的扩展方程为： (7) 由于相邻两单元公共节点上的基本场变量（位移）相同，所以可将扩展后的各单元方程相加。 将式(6)和式(7)相加得： 上述组集过程可记为： (9) 组集后的结果简记为：Kq = P 式中，K称为总体特性矩阵(常称为总体刚度矩阵和总刚阵)，P称为总体节点载荷列阵。需指出的是，对单元的一个公共节点而言，除了有相邻单元作用于该节点的力之外，还可能有做用于该节点的外载荷。若一节点上无外载荷作用(如本例