第二次课-2杆的有限元分析摘要.ppt

作业 用有限元法求图示受拉阶梯杆的位移和应力。已知杆截面面积A(1)=4×10-4m2，A(2)=2×10-4m2， ，A(3)=1×10-4m2各段杆长l(1)=l(2) =l(3)=0.1m；材料弹性模量E(1)= E(2)= E(3)=2×107Pa，作用于杆端的拉力F4=10N。建立结构的有限元方程，并求解相应的变形、应变、应力及支反力。 ① u1 u2 u3 E(1), A(1) ② 1 2 3 E(2), A(2) l(1) l(2) R1 u4 ③ F4 4 E(3), A(3) l(3) 弹性力学的有限元分析计算可分为三个步骤: 1 结构离散化 这是有限元法的基础，用由有限个方位不同但几何性质及物理性质均相似的单元组成的集合体来代替原来的连续体或结构。每个单元仅在节点处和其他单元及外部有联系。对于不同的问题，根据自身的特点，可选用不同类型的单元。对同一问题也可以分别或同时选用多种单元。 例：一个受载的悬臂梁和用三角形单元离散化的模型 真实系统 有限元模型 离散 有限元模型由一些简单形状的单元组成，单元之间通过节点连接，并承受一定载荷。 注意：1)节点是有限元法的重要概念，有限元模型中，相邻单元的作用通过节点传递，而单元边界不传递力，这是离散结构与实际结构的重大差别； 2）节点力与节点载荷的差别 单元：即原始结构离散后，满足一定几何特性和物理特性的最小结构域 节点：单元与单元间的连接点。 节点力：单元与单元间通过节点的相互作用力。 节点载荷：作用于节点上的外载。 节点自由度(DOFs) ：用于描述一个物理场的响应特性。 分离但节点重叠的单元A和B之间没有信息传递（需进行节点合并处理） 信息是通过单元之间的公共节点传递的。 具有公共节点的单元之间存在信息传递 非法结构离散 不同材料 节点不合法 典型单元类型 单元类型 单元形状 节点数 节点自由度 杆单元 2 2 梁单元 2 3 平面单元 3 2 平面四边形 4 2 轴对称问题 3 2 板壳单元 4 3 四面体单元 4 3 2.单元分析 主要内容：由节点位移求内部任意点的位移，由节点位移求单元应变，应力和节点力。 3.整体分析 (1) 由节点平衡方程，建立以整体刚度矩阵[K]为系数的，整体节点位移{d }和外载{R}的关系式——整体平衡方程。 (2) 考虑几何边界条件，修改总体刚度矩阵，求解全部未知位移分量。 a)受拉阶梯杆示意图 ① u1 u2 u3 E(1), A(1) ② F3 1 2 3 E(2), A(2) l(1) l(2) R1 x 二. 有 限 元 求 解 基 本 原 理(一维问题) 引例:用有限元法求图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是节点号。取节点位移作为基本未知量，应力由求得的节点位移算出。 c)单元图 b)有限元模型 ① u1 E(1), A(1) Node 1 Node2 l(1) R1 u2 u2 u3 ② F3 Node 3 E(2), A(2) l(2) Node 2 ui Ae Ee Node i Node j le Pi uj Pj x 图5-6 ② ① 1 2 3 2.确定单元插值函数(形函数) 有限元法将整个求解域离散为一系列仅靠公共节点联结的单元，而每一个单元本身却视为光滑连续体。单元内任一点的场变量(如位移)可由本单元的节点值根据场变量在单元中的假定分布规律(插值函数)插值求得。 本例中，每单元有两个节点，采用线性插值。图c是一典型单元图，两节点分别为i和j，节点场变量值分别记为ui和uj 。设单元中坐标为x处的场变量为u(x) 。 单元的位移场为u(x)， 由两个端点的位移来进行线形插值确定，设u(x) 为： (1.a) 单元节点条件： (1.b) 将节点条件(1.b)带入(1.a)，可以求得a0和a1： (1.c) 则 其中N(x)叫做形状函数矩阵(shape function matrix),为 qe叫做节点位移列阵