第四章 平面问题有限单元法.pptVIP

计算实例 弹性模量为E， ， 不计自重。 求各角点的位移及应力。 如图所示一块薄板，长度和宽度分别为2米、1米，厚度为t，一端固定，一端受均布拉力q， 1、离散结构物 为了计算简单，划分为2个单元，单元号和节点编号如图所示。 对节点进行编号时，应使同一单元内的节点编号间差值最小。 2、选择单元的位移模式 因为采用的是三节点的三角形单元，位移模式已确定。 3、计算单元刚度矩阵 r = i, j, m s = i, j, m 对于平面应力问题，三角形单元的刚度矩阵 其中，? — 三角形单元的面积，=1/2*2*1=1 t — 三角形单元的厚度 对于单元? 单元?的刚度矩阵 表示为 所以单元?的刚度矩阵 代入上式，可得 单元?的刚度矩阵 表示为 对于单元? 可得 所以单元?的刚度矩阵 4、计算等效节点力 面积坐标对直角坐标的积分 表面分布力的等效节点力 所以，各节点的等效节点力 总体刚度矩阵 因为结构物共有4个节点，所以 分别将[K]1、[K]2扩展为4?4矩阵 5、组装整体结构物的刚度矩阵 [K ] 所以总体刚度矩阵[K] 6、形成结构物的求解方程，并解出位移{δ} 因为节点3、4的约束边界是固定约束，即 u3 = v3 = 0， u4 = v4 = 0，经过划行划列降阶处理后，得 解得 计算单元应力 单元? 得 单元? 得 学习动物精神 11、机智应变的猴子：工作的流程有时往往是一成不变的，新人的优势在于不了解既有的做法，而能创造出新的创意与点子。一味 地接受工作的交付， 只能学到工作方法 的皮毛，能思考应 变的人，才会学到 方法的精髓。 学习动物精神 12、善解人意的海豚：常常问自己：我是主管该怎么办才能有助于更好的处理事情的方法。在工作上善解人意， 会减轻主管、共 事者的负担，也 让你更具人缘。 单元内有虚位移，引起单元内的虚应变 根据虚位移原理 单元刚度矩阵 由于节点虚位移和节点位移都是常量，可提到积分号外面， 令 由于虚位移是任意的，所以欲使上式成立，必须有 则 — 单元刚度方程 Ke — 单元刚度矩阵 刚度矩阵可分为三行三列的子矩阵 3节点三角形单元刚度矩阵的显式 对于平面应力问题 r =i, j, m s =i, j, m 对于平面应变问题 单元刚度矩阵的性质 单元刚度矩阵是对称矩阵 单元刚度矩阵的主对角元素恒为正值 单元刚度矩阵是奇异矩阵 单元刚度矩阵仅与单元本身有关 单元所受的非节点载荷向节点移置，移置到节点后的载荷称为等效节点载荷。 载荷移置必须按照静力等效的原则进行。 载荷移置的方法： §4-4 载荷移置 普遍公式法 直接法 用虚功原理确定等效节点载荷 若三角形单元上作用有集中力P、分布力q(力/面积)和体积力g(力/体积)，则根据静力等效原理，节点力所做的虚功等于三种力所做的虚功。 载荷移置的普遍公式法 代入上式，得 由此可知 由体积力引起的等效节点载荷 由表面力引起的等效节点载荷 由集中力引起的等效节点载荷 由于虚位移是任意的，所以 集中力的等效节点载荷计算 由于 表面分布力的等效节点载荷 由于 体积力的等效节点载荷 由于 单元特性分析 → 建立单元刚度矩阵 单元载荷移植 → 建立节点载荷列阵 根据结构各节点的静力平衡条件，得出结构总体刚度方程。 §4-5 结构总体刚度方程 结构总体刚度方程的建立 ① ② 节点载荷列阵 节点位移列阵 单元的节点力 ① ② 单元①的节点力列阵 单元②的节点力列阵 单元节点力与节点位移之间的关系 对于单元① 对于单元② 由节点的静力平衡条件 将单元① 和②的节点力与节点位移之间的关系代入上式，得 用矩阵形式表示，有 简写为 K 称为总体刚度矩阵 1、按单元形成结构总体刚度矩阵 形成结构总体刚度矩阵的方法 a、将存放结构总体刚度矩阵的数组清零 b、计算第一个单元的刚度矩阵，叠加到相应位置 c、按顺序依次叠加各个单元的刚度矩阵，直至最后一个单元 d、最终得到的即为结构总体刚度矩阵 2、按节点形成结构总体刚度矩阵 a、将存放结构总体刚度矩阵的数组清零 b、从第一个节点开始，检查该节点与哪些节点相邻，是哪些单元的公共节点，写出第一行的子矩阵 c、按照节点编号顺序，对每个节点重复以上步骤，到最后一个节点 d、最终得到的即为结构总体刚度矩阵 结构总体刚度矩阵的性质 结构总体刚度矩阵是一个对称矩阵 结构结构总体刚度矩阵是一个奇异矩阵 结构总体刚度矩阵