有限元方法课件 第四章 杆系结构有限元.ppt

第四章 杆系结构有限元 有限元法主要优点： （1） 概念浅显，容易掌握。（离散、插值、能量原理、数学分析） （2）适用性强，应用范围广，几乎适用于所有连续体和场问题的分析。（结构、热、流体、电磁场和声学等问题） （3）计算规格化（采用矩阵表示），便于计算机编程。 例4-1 一维拉杆 图示阶梯形直杆，各段长度均为，横截面积分别为3A，2A，A，材料重度为γ，弹性模量E。求结点位移和各段杆中内力。 离散化：将单元划分为3个单元，4个结点。 单元刚度矩阵： 等效结点荷载：按静力等效原则，有: 整体结构平衡方程 例4-2 变截面梁 有一变截面梁，一端固定，另一端铰支。梁长为2l，固支端的截面尽寸为b×1.6h，铰支端的截面尺寸为b×h。梁上作用均布载荷p0。求梁端的约束反力。 离散化 将梁划分成2个单元，3个结点。每个单元 长度为,截面取平均截面。 对号入座，组合整体刚度矩阵 荷载等效结点力向量 引入边界条件 3? 用式（4-77）算A1、A2、A3 和A。 本部分内容结束 只需要将式（5-48）中的Iz换成Iy，并注意编号次序。同时考虑到式（5-51），即得 （4-52） 4、空间梁单元(BEAM4) 空间梁单元，每个节点有6个自由度，单元自由度为12。下图给出了空间梁单元节点位移分量的正方向及其编号。单元力的正向及其编号与单元位移相同。 i j l x y z ui ?xi vi ?yi wi ?zi vj ?yj wj ?zj uj ?xj i j l x y z 1 4 2 5 3 6 8 11 9 12 7 10 综合前述结果，得空间梁单元单元坐标单元刚度矩阵(式4-53）～（式4-57）。 （4-53） 轴弯剪(4-48）、弯剪扭（4-52） （4-54） （5-55） （5-56） （5-57） 在前面两节中，单元位移和单元力都是按单元坐标系的坐标轴分量定义的，由此建立的单元刚度矩阵属于单元坐标单元刚度矩阵。 进行系统分析时，需要把单元力按统一的结构坐标轴的分量表示出来，以便建立结点平衡方程。因此，在进行系统分析之前，必须把单元坐标系中的单元力以及单元刚度矩阵都转换到结构坐标系中去。此外，还需要把结构坐标系中的节点位移转换到单元坐标系中去，以计算结构内力。这一转换过程称为坐标变换。(一维杆和弯剪梁单元不需要坐标变换，因两种坐标系统方向一致) §4-5 坐标变换 结构坐标符号约定： ???——结构坐标单元位移 ?F?——结构坐标单元力 [k]——结构坐标单元刚度矩阵 单元坐标中的符号约定： ???e——单元坐标单元位移 ?F?e——单元坐标单元力 [k]e——单元坐标单元刚度矩阵 1、坐标变换矩阵定义 把单元位移从结构坐标系转换到单元坐标系的变换矩阵定义为坐标变换矩阵，用符号[T]表示。有 式（4-58）给出了结构坐标单元位移转换为单元坐标单元位移的转换式，同时是坐标变换矩阵[T]的定义式。 （4-58） 2、结构坐标单元力 单元力在单元位移上作的功，不因其坐标系的改变而变。则有 将式（4-58）代入， 消去???，得 对上式两端进行转置，注意到 即得 （4-59） 式（4-59）表明：结构坐标单元力等于单元坐标单元力前乘坐标变换矩阵的转置。 在单元坐标系中，有 3、结构坐标单元刚度矩阵 上式两端左乘[T]T， 注意到式（4-58）、（4-59），有 （4-58） （4-59） [k]——结构坐标单元刚度矩阵。 得 （4-60） 式（4-60）给出了把单元坐标单元刚度矩阵转换为结构坐标单元刚度矩阵的转换式。 引入 坐标变换矩阵因单元类型不同而异。 1、平面铰接杆单元（桁架元） 设OXY为结构坐标，oxy为单元坐标。?为从单元 i 端出发的任一矢量。它在结构坐标系中的分量为?X、?Y, 在单元坐标系中的分量为?x、?y。结构坐标系中的分量?X、?Y 在单元坐标x轴上投影的代数和给出?x 。同理， ?X、?Y 在单元坐标y轴上投影的代数和给出?y 。 §4-6 坐标变换矩阵 （4-61） 写成矩阵形式， ia+ab x y i ? X Y ?X ?Y ?x ?y ? b a c d cd-ca 取 ——i节点在单元坐标系中的位移向量 ——i节点在结构坐标系中的位移向量 x对X、Y的方向余弦 y对X、Y的方向余弦 同理可得单元j节点在单元坐标系和结构坐标系中的位移向量： 有 组合上述结果，得平面铰接杆单元的单元坐标单元位移和结构坐标单元位移之间关系： i、j两节点间的位移变换