杆系结构有限元.pptVIP

第1章 杆系结构有限元 第1章 杆系结构有限元 1.1 有限单元法及其发展概况 1.2? 有限元方法及解题步骤 1.3? 单元分析? 1.4? 整体分析 1.5 按单元定位向量形成总刚度方程 1.6 约束处理及求解 一、对象、任务 对象：研究有限大小的个体(element) 任务：1. 建立应变与结点位移分量之间的关系; 2. 建立应力与结点位移分量之间的关系; 3. 建立结点力与结点位移分量之间的关系; 4. 把作用在单元内的外载转化成结点荷载。 位移函数的一般形式 形函数与形函数矩阵 3、几何矩阵[B] 4、应力矩阵[S] 5、用虚功原理建立单元刚度方程 设单元在杆端力{F} 以及荷载{f(x)}的作用下处于平衡状态，在某种可能的虚位移影响下，杆端虚位移为{d}*，相应的虚应变为{ε}* 。 由虚功原理: 外力虚功的总和为: T={d}*T{F} + {d}*T[N]T{f(x)} 实际应力{σ}在虚应变上所做的内力虚功总和为： 将{ε}和{σ}的表达式代入上式得 由虚功原理（T = W），并考虑到{d}*T为任意已知虚位移，整理可得： 简写成： 式中 为单元刚度矩阵的一般表达式。 6、用最小势能原理建立单元刚度方程 设单元在杆端力{F} 以及荷载{f(x)}的作用下处于平衡状态，在某种可能的位移影响下，杆端位移为{d}*，相应的应变为{ε}* 。 外力势能的总和为: T={d}*T{F} + {d}*T[N]T{f(x)} 应变为{ε}*与应力{σ} *对应的变形能为： 将{ε} *和{σ} *的表达式代入上式得 单元刚度矩阵的一般表达式。 平面梁单元刚度矩阵 等效节点力 单元分析过程总结： 1.设定位移函数，建立单元内的位移和节点位移的关系－形函数[N] 2.建立应变和节点位移的关系－几何矩阵[B] 3.建立应力和节点位移的关系－[S]=[D][B] 4.利用虚功原理－单元刚度矩阵、等效节点力 作业1： 有一个等截面两节点二力杆，杆长为L，截面积为A，材料弹性模量为E。每个节点只考虑一个水平位移，对于图（a）、（b）所示的坐标系统和位移插值函数，分别求相应的[B]矩阵和单元刚度矩阵[K]。 作业2： 有一个等截面三节点轴力杆单元，杆长为L，截面积为A，材料弹性模量为E。每个节点只考虑一个水平位移，求相应的[B]矩阵和单元刚度矩阵[K]。 考虑剪切变形的梁单元 1）在经典梁单元基础上引入剪切变形影响 Timoshenko梁单元 剪切闭锁 减缩积分 假定剪切应变 自学内容： 等效节点力可按下式计算： 其中： 如：满跨横向均布荷载q作用： 距i节点x1处集中力偶M作用： 其他类型的单元也可采取相同的方法进行推导 思考题： 形函数形式简单，而结构受力有时很复杂，简单的形函数能表示复杂的变形吗？ i j x, u L, EA o (a) (b) i j x, u o L/2 L/2 u=α1+α2x2 u=α1x+α2x2 为了保证解收敛，位移函数必须满足： 1）位移函数必须能反映单元的刚性位移 2）位移函数必须能反映单元的常量应变 3）位移函数在单元内部必须是连续函数 收敛的必要条件 作业3： 有一个圆形截面两节点扭转杆单元，杆长为L，截面直径为D，材料弹性模量为E。每个节点只考虑一个扭转位移，求相应的[B]矩阵和单元刚度矩阵[K]。 前面讨论的梁单元基于变形前垂直于中面的截面变形后仍保持垂直的Kirchhoff假设。通常所述的梁单元，即是指这种单元。它在实际中得到广泛的应用，一般情况下也能得到满意的结果。但是应该指出，它是以梁的高度远小于跨度为条件的。因为只有在此条件下，才能忽略横向剪切变形的影响。 理论解： 但是在工程实际中，也常常会遇到需要考虑横向剪切变形影响的情况。例如高度相对跨度不太小的高梁即属此情况。此时梁内的横向剪切力Q所产生的剪切变形将引起梁的附加挠度，并使原来垂直于中面的截面变形后不再和中面垂直，且发生翘曲。但在考虑剪切变形的Timoshenko梁弯曲理论中，仍假设原来垂直于中面的截面变形后仍保持为平面。 图中g 表示截面和中面相交处的剪切应变，并且有如下关系式，即 其中q 是截面的转动。在经典的梁弯曲理论中，忽略剪切应变，即认为g =0，所以 即截面的转动等于挠度