计算结构力学第四章.ppt

第四章 单元刚度矩阵 二、选择位移函数应遵守的准则 1、允许发生刚体位移(不产生单元自应变)； 2、能反应常应变。 令 则单刚便可写成简洁的形式，以利编程。 3. Kii0 * 首页 上页 返回 下页 * 计算结构力学 　　形成单元刚度矩阵是整个结构分析中的一个重要环节。 　　静力法推导利用了结构力学中的转角位移方程，也是采用了Euler梁理论的结果。 Euler梁：简单梁 　　有限元分析的计算精度在很大程度上取决于单元刚度矩阵，也就是取决于单元形状函数(位移函数)的选择。 4-1 概 述 单元位移函数： 　　指单元位移场，或称形状函数。 　　在形成单元刚度矩阵前需事先假设，如梁单元可假设为三次多项式等，位移函数一般选取代数多项式。 一、单元位移函数 Kij的定义 　　[K]的某一行的元素表示{δ}中六个分量分别发生单位位移时所引起{F}中某一分量的值； 　　 [K]的某一列的元素表示{δ}中某一分量发生单位位移时所引起{F}中六个分量的值。 三、单元刚度矩阵中行列元素的物理意义 4-2 平面刚架的单元刚度矩阵(静力法推导) 　　分别采用静力法、能量法(虚位移原理和势能驻值原理)推导局部坐标系下的单元刚度矩阵。为方便计，均省略“ ”的记号。 1、单刚是对称的，Kij=Kji ① ② 反力互等定理 讨　论： 2、单刚是奇异的，不存在逆矩阵，如第一行(列)加第四行(列)等于零，或第二行(列)加第五行(列)等于零。 　 　　奇异性的说明：没有约束的自由体是不能求解位移的，故要引入足够的约束条件消除刚体位移，这样也就消除了矩阵的奇异性。 4-3 利用能量原理推导轴力杆单元的刚度矩阵(积分法) 1、选择位移函数 　　考虑图示桁架单元(等截面直杆)在杆端力作用下杆内应力应变均为常量，根据轴力杆单元的几何条件(应变-位移关系)，故可将单元位移场选为： 　u=a1+a2x 　　　(1) 其中：u=u(x)即为单元位移函数(单元位移模式) 式(1)中： 　　a1——刚体位移项； 　　a2——常应变项。 　　满足单元位移模式的选择准则。将(1)式写成矩阵形式 　　　　　　u={X}T{a} (2) 式中： {X}T={1 x} 　　 {a}T={a1 a2} x i j x y y x i j 1 0 (11) (12) (13) 5．①采用虚位移原理推导 任给单元一个虚位移δu*(x)，并设结点虚位移为δ{δ*}，于是虚应变能可写为 (14) (15) (16) (17) (18) (19) (20) 4-4 利用能量原理推导梁单元的刚度矩阵(积分法) 　　仍采用虚位移原理和势能驻值原理进 行推导，现讨论较复杂的等截面直梁单元 (如图)，采用Euler梁理论。 1、选择位移函数 因不考虑分布荷载(均简化为等效结点荷载处理)，故有EI’’’’=q(x)=0，由此得 ∴ v=v(x)=a1+a2x+a3x2+a4x3 ={X}T{a}　　　　　　 (1) v(x)满足刚体位移及常应变准则。 (2) (3)