上海交大计算结构力学课件ppt平面问题的有限元法01.doc

CHAPTER02 平面问题的有限元 2.1 引言 平面应力和平面应变问题的有限元分析 位移有限元建立的一般方法 三角形单元、矩形单元 线性单元、高阶单元 轴对称单元 2.2 弹性力学平面问题基本方程 平面应力和平面应变 杆系结构: 长度远大于两个方向的尺度 变形在轴向采用了平面假设 只有一个未知函数 非杆系: 比如深梁,平面假定不再成立 图2-1 平面应力 平面应力: 厚度很小, 载荷(包括位移边界)和平面平行,沿均匀分布,因此可以近似地认为沿向所有的应力分量为零. 图2-2 平面应变 平面应变: 横向尺度远小于纵向尺度, 载荷(包括位移边界)和平面平行,沿均匀分布,认为沿位移分量为零. （1）平面应力: 沿厚度很薄，可以认为沿厚度均匀 各个分量和无关，只和有关 板的两个面上没有外力， 因此两个面上有： 也即对薄片所有点有： 并进一步有: 于是： 注意: （2）平面应变: 沿没有位移， 只有和不为零. 但 ， 平面应力和平面应变的关系 以下仅研究平面应力 2.2 平面问题的三角形单元格式 2.2.1 三角形单元 任意的一个平面采用三角形单元离散 节点编号 单元节点位移： 如果单元节点位移已知，如何给出单元内部任意一点的位移， 其中 FEM的思想：单元内部插值 2.2.2 单元分析 若节点位移已知，如何求单元内部任意一点的位移，称为位移模式或位移函数 一般取为多项式，多项式选取：由低向高 为待定 在节点上应满足节点位移 依照方程组求解的行列式方法，可以求出，它应该是节点坐标和节点位移的函数。 其中A是单元三角形的面积 上述相同 同理可以求出 其中A是单元三角形的面积 上述相同 代入假定的位移表达式 经过演算： 也可以同上。 最后给出位移的形式为： 其中： 或： 或 由于N是线性函数，因此任意一点的位移就表现为单元节点位移（取决于）的线性组合（取决于N）或就是插值，因此N也叫做插值函数。 整个有限元的核心就是构造N和 N所具有的性质（也是今后构造N所必须满足的条件） （1） （2） 任意点插值函数的和为1 意味着如果单元三个节点位移相同（具有刚体位移），则内部任意一点也有这个位移。 从而有： 由于单元受相邻单元的影响，必然产生刚体位移，因此我们在构造N的时候必须使得位移模式包含刚体位移。 （3） 单元与单元边界上的位移是唯一的。 保证唯一的连续性，否则单元之间出现间隙或重叠。（验证留做作业） 将位移、应变、应力表示为节点位移的函数。 注意由于是和的线性函数，因此是常数 因此我们得到的常应变单元和常应力单元。 作业 2．1 验证上述三角形单元边界上位移的连续性。应变是否连续？ 2.2.3三角形单元有限元平衡方程-最小势能原理 其中： 单元的刚度矩阵： 可以看出，此时为矩阵 体积力等效节点载荷列阵： 表面分布力等效节点载荷列阵： 根据最小势能原理： 于是得到总体平衡方程为： 或： 讨论： (1) 关于等效节点力概念的解释： 外力做功= 因此可以理解为假想的作用在节点上的力， 其效果是做功的大小不变， 所以称为等效节点载荷。 (2) 关于单元刚度矩阵的力学含义 只考虑一个单元 如果该单元产生的位移为： （ B1） 上述平衡方程有： 这给出了刚矩阵各元素的力学意义。由于产生（B2）位移的载荷应该平衡， 于是进一步有： （3）关于对称性 由于的对称性， 因此也对称。 （4）刚度矩阵的奇异性 行列式值为零， 力学解释：外力下平衡的单元可以有任意的刚体位移， 从而解不唯一。因此要想真正给出解答，必须将位移边界约束条件反映到平衡方程里。 2.2.5 等效节点载荷列阵 （1） 体积力 如果存在重力（一类体积力）的作用，则 计算后 （2）分布力 作用在垂直于轴的三角形边界（边长为）上。 讨论：如果沿长度为线性分布， 则等效节点力会怎样？ 该问题为作业。 2.2.6 整体刚度矩阵的合成 举例说明： （1） 预备知识-关于矩阵分块 其中 ……. （2）单元刚度矩阵的扩维 首先将单元刚度矩阵和单元节点位移列阵的由单元维数扩充为总体维数。 先将每个单元的刚度阵表示为如下的形式， 其中上标为单元号，下标为单元对应的节点。后面将看到，这个下标也是该单元刚度矩阵在这个位置的元素组装到总体刚度矩阵时的位置。 此外还引入矩阵的分块表示方法。 上述矩阵中的各元素均为，0实际表示的是。 （3）总体刚度矩