平面刚架问题的有限元法-公开课件.pptVIP

1 有限元法基础Finite Element Method 四 总刚矩阵集成 组装各单元在总体坐标系中的刚度矩阵 [K ] =Σ[K]e 五 节点载荷列阵 各类载荷移置为单元等效节点载荷后，叠加得节点载荷列阵{F} 。 约束处理和求解线性方程组 经约束处理消除[K]的奇异后，求得位移{δ}e。从而求得局部坐标下的位移 例2-2 变截面梁 有一变截面梁，一端固定，另一端铰支。梁长为2l，固支端的截面尺寸为b×1.6h，铰支端的截面尺寸为b×h。梁上作用均布载荷p0。求梁端的约束反力。 离散化 将梁划分成2个单元，3个结点。每个单元 长度为,截面取平均截面。 对号入座，组合整体刚度矩阵 荷载等效结点力向量 引入边界条件 有限元法分析过程概述 1. 结构物的离散 3. 分析单元的力学特性 {R}e= [K]e {δ}e {R}e — 单元节点力 [K]e — 单元刚度矩阵 4. 计算等效节点力 5. 组装整体结构物的刚度矩阵 [K ] [K ] =Σ[K]e 6. 形成结构物的求解方程，并解出位移{δ} (6) 等价节点力 四 总刚矩阵集成 组装各单元在总体坐标系中的刚度矩阵 [K ] =Σ[K]e 五 节点载荷列阵 各类载荷移置为单元等效节点载荷后，叠加得节点载荷列阵{F} 。 约束处理和求解线性方程组 经约束处理消除[K]的奇异后，求得位移{δ}e。从而求得局部坐标下的位移 (5) 单元坐标单元刚度矩阵 梁单元刚度矩阵公式为 （2-48） x y i j l 图2-9 qy(x) qx 将式弯剪梁（2-36）、一维杆（2-11）膨胀成6×1矩阵后相加，并注意到式（2-43），有 （2-36） （2-11） 一维杆 弯剪梁 最后得等价节点力矩阵 （2-47） Qj Qi Mj Nj Mi Ni 荷载分布 表2－2 几种横向分布荷载等价节点力 i j qy qx qy i j qx qy i j qx 1、坐标变换矩阵定义 把单元位移从结构坐标系转换到单元坐标系的变换矩阵定义为坐标变换矩阵，用符号[T]表示。有 （2-58） 2、结构坐标单元力 （2-59） 3、结构坐标单元刚度矩阵 （2-60） 三 坐标变换 i j x y ?2 ?3 ?5 ?6 l ?1 ?4 于是得到： X Y 由于?1 、 ?2 、?4、 ?5的性质和平面铰接杆相同，因而有相同的[T]矩阵。又因单元坐标系xy平面和结构坐标系XY平面在同一平面上， 因而单元坐标系z轴和结构坐标系的Z轴总有相同指向，所以恒有： 3)轴剪弯平面梁单元（刚架） （2-63） Institute of Mechanical Engineering and Automation Long yuhong, Ph D School of Mechantronic Engineering, Guilin University of Electronic Technology, Guilin, 541004, P.R.C. TelE-mail: longyuhong02@ Institute of Mechanical Engineering and Automation IMEA §2.2 平面刚架问题的有限元法 梁：两端刚性连接的杆件 固定连接 拉压、弯曲、扭转 与截面大小、形状、方位有关 刚架：梁单元组成的系统 平面刚架：所有杆轴线处于同一平面的系统 一 结构离散 梁单元离散，各杆有局部坐标系，方向为杆的轴线方向。结构的刚度方程是在统一的坐标系（总体坐标系）中建立并求解，因此需将单元局部坐标各量转换到总体坐标系中。 取杆件与杆件交点、集中力作用点、杆件与支承的交点为节点。相邻两节点间的杆件段是单元。节点编号时力求单元两端点号差最小。 二 单元分析 1、两端承受剪力、弯矩的平面梁单元 i j x y i j x y ?1 ?2 ?3 ?4 l F1 F2 F3 F4 l （1）局部坐标下单元位移和单元力 ① 单元位移 （2-24） 其中， v——y方向位移，即挠度。