结构有限元分析1章.pptVIP

绪论 有限元方法是求解一般连续域问题的数值解法 特点： 1、适用于计算机处理，可充分发挥计算机的 优势 2、可解决多学科、多物理问题 3、概念清晰、易于掌握 有限元法的基本思路 将连续的求解域离散成由节点互相联系的有限个单元，用节点的未知参数表征单元特征，用适当方法组合全部单元形成包含节点未知参数的线性方程组，求解得到待求的节点参数。 单元节点数是有限的，单元大小是有限的，单元数是有限的。 求解格式 1、求解对象划分成单元，离散化； 2、单元特征分析（建立单元的刚度矩阵和刚度方程）； 3、单元集合成整体模型（大型方程组） 4、数值求解，得节点未知量。 第一章 杆系结构的有限元 1、杆系结构是工程结构中常见形式； 2、每一杆件可视为一个单元，受力与位移关系明确 。 1.1直梁问题 一变截面梁，外载Z、M作用 在B点，求解节点位移， 计算单元内力、应力。 建立有限元模型 节点编号：1，2，3，4 单元编号：（1），（2），（3） 直梁结构的有限元分析 直梁有限元模型包括n个节点，任一节点有位移分量分量 称为节点位移{δi},节点的外载Zi,Mi, 称为节点载荷{Qi}。 对于直梁结构，要先求出全部节点载荷作用下引起的全部节点位移，进而再求出单元内的位移和内力。 全部节点位移记为{δ}，全部节点载荷记为{Q}，{δ}和{Q}各有2n个分量，称该结构具有2n个自由度 直梁的单元分析 是e单元的刚度矩阵，此处为4X4的方阵， 每个元素都是常数。 可用材料力学方法解出 的各元素： = 容易看出1、 为对称矩阵； 2、 元素仅与E,J, l 有关。 分块表示： = = = (1.1)式可表示为： 展开得矢量方程： 每个矢量方程包含2个代数方程。 单元的集合与刚度矩阵叠加 原结构共有3个单元，每个单元分别有 有限元方程依据节点平衡条件建立 节点2 的平衡条件为： ,简写为 代入 , 得 同理，由节点3的平衡方程，有 梁端点 将各式合并得 简写为 （1.6） 该矩阵式即结构的有限元方程（位移方程） 其中 ——整个结构的刚度矩阵（总刚） 矩阵的集合方法： 1、单元刚度矩阵 的子块对号入座； 2、 扩大至结构总刚阵同阶； （共m个单元） 对于m个单元，n个节点，式（1.6）应有2n个方程 为 的对称矩阵 单元刚度矩阵 与整体刚度矩阵 ： 表示单元抵抗变形的能力 表示整个结构抵抗变形的能力。 引入边界约束与求解 图1.1所示结构中有4个节点，8个自由度， 8个位移分量中 =0 8个节点载荷中 为约束反力，其余5项已知 必须将已知位移条件引入（1.6）中，方程组才是可解的 1.2 平面刚架问题 1、2、1 单元与节点 平面内，每个节点有3项位移分量 每一节点的节点载荷也有3项分量 每一节点有3个自由度，如刚架有n个节点，则有3n个自由度 表示全部节点位移， 为全部节点载荷 都具有3n个分量。 1、2、2 单元刚阵（相对局部坐标） 为分析方便， 建立单元的局部坐标系 ，在局部坐标系下，e单元的节点位移表示为 e单元的节点力为 在小变形下， 仅与T有关，f, 仅与q,m有关 分析单元的位移与受力的关系，可以分别建立轴向变形和弯曲变形的单元刚度矩阵。 对于弯曲变形，与直梁的弯曲相同 对于轴向变形，有 将两部分合并，得 仍可简写为 这里 为平面刚架单元在局部坐标系下的单元刚阵，为 方阵，对称矩阵。 1.2.3 单元刚阵的坐标变换 各单元在局部坐标系下建立的 需要经坐标变换至公共的坐标系下才能叠加形成结构的整体刚度矩阵 设 为统一坐标系 e单元的局部坐标系为 以位移为例：节点在坐标系 中位移为 在坐