2-1杆梁问题的有限元法.ppt

有限元法基础Finite Element Method 第二章 杆件静力学问题有限元法 要 点 了解杆件系统有限元法的解是精确解； 杆件系统中的单元都在单元局部座标中讨 论，然后再投影到整体座标中去； 杆件系统中的单元是梁单元，其变形以弯曲 和拉伸变形为主，常忽略剪切变形的影响。 第二章 杆件静力学问题有限元法 引 言 二、组合弹簧的刚度矩阵 第二节 杆件问题有限元法 所研究的对象是细长的杆件，即轴线方向的尺寸远比其他二个方向的尺寸大的多，如轴支柱、加强筋等； 杆件系统可分为平面杆件系统和空间杆件系统，若组成结构物的杆件都在同一个平面内，外力也在这一平面内，则称为平面杆件系统； 若组成结构物的杆件不在同一个平面内，则称为空间杆件系统； 杆件系统中所用的单元主要为梁单元，如果各杆件只受拉压作用则采用桁架单元。 例2.1 一维拉杆 图示阶梯形直杆，各段长度均为l，横截面积分别为3A，2A，A，材料重度为γ，弹性模量E。求结点位移和各段杆中内力。 离散化：将单元划分为3个单元，4个结点。 单元刚度矩阵： 等效结点荷载：按静力等效原则，有: 整体结构平衡方程 座标变换 令λ＝cos(a),μ=sin(a),则矩阵相乘可写出 刚度矩阵的性质 单元刚度矩阵：对称性；奇异性， 即[K]e =0。 总体刚度矩阵：对称性；奇异性；稀疏性；非零元素呈带状性。 2.4 桁架结构的有限元法分析简例 刚度矩阵的求法 2.4 桁架结构的有限元法分析简例 2.4 桁架结构的有限元法分析简例 2.4 桁架结构的有限元法分析简例 2.4 桁架结构的有限元法分析简例 2.4 桁架结构的有限元法分析简例 2.4 桁架结构的有限元法分析简例 2.4 桁架结构的有限元法分析简例 2.4 桁架结构的有限元法分析简例 2.4 桁架结构的有限元法分析简例 2.4 桁架结构的有限元法分析简例 4. 形成方程 有限元法分析过程概述 1. 结构物的离散 3. 分析单元的力学特性 {R}e= [K]e {δ}e {R}e — 单元节点力 [K]e — 单元刚度矩阵 4. 计算等效节点力 5. 组装整体结构物的刚度矩阵 [K ] [K ] =Σ[K]e 6. 形成结构物的求解方程，并解出位移{δ} (6) 单元坐标单元刚度矩阵 对于等截面铰接杆单元， （2-23） 所讨论的[K]e,{R}e,{ε}, {σ} 都是在单元的局部座 标下进行的。由于杆件系统在空间中各有自己的方向，都 建立了各自的座标系，所以在分析中都必须统一在一个整 体座标中进行，一定要进行座标变换。 1、坐标变换矩阵定义 把单元位移从结构坐标系转换到单元坐标系的变换矩阵定义为坐标变换矩阵，用符号[T]表示。有 （2-58） 2、结构坐标单元力 （2-59） 3、结构坐标单元刚度矩阵 （2-60） 三 坐标变换 坐标变换矩阵因单元类型不同而异。 1)、平面铰接杆单元（桁架元） X Y ?X ?Y ?x ?y ? b a c d x y i ? （2-61） ia+ab cd-ca 取 ——i节点在单元坐标系中的位移向量 ——i节点在结构坐标系中的位移向量 x对X、Y的方向余弦 y对X、Y的方向余弦 写成矩阵形式， 有 是在j节点有单位位移时，而在i节点所需施加的力。 图5-1 简例结构图 由两根杆件组成的桁架，杆件的截面积都为A，弹性模量为E，长度为l1及l2，在节点处受有外力Fx1，Fy1，Fx2，Fy2，Fx3，Fy3，求各节点的位移。 分析步骤： 1. 将结构物离散为有限个单元 分为2个单元，第一个单元的节点编号为1和2，第二个单元的节点编号为2和3。对于第一单元，在第1、2节点处的节点力为 ，表示节点施加在单元1上的节点力。相应的位移为 。下标表示节点的编号，上标表示单元的编号。同样第二单元上节点2、3处的节点力和节点位移分别为 和 图2-2 简例单元图 图2-3 单元特性图 2. 研究单元特性 在单元分析中，主要是建立节点位移分量和节点力的关系式。 令u11=1，v11= u12= v12= 0，则杆1的变形量为1×cos?，由材料力学可知，杆1所受的轴向压力P