10矩阵法3.ppt

1 2 3 4 5 A B C D 1 2 3 4 5 6 7 8 （3）建立整体座标系单刚 e [k] 对称 1 2 5 6 5 6 7 8 对称 5 6 7 8 1 2 5 6 对称 3 4 7 8 3 4 7 8 对称 3 4 5 6 3 4 5 6 1 2 7 8 1 2 7 8 对称 * [K]{?} = {F} 总刚度方程所表示的内容与应用 1：已知所有结点位移，计算所有约束力 2：已知所有结点约束力（仅有的），计算结点位移和内力。（如习题3-1） 已知任意载荷计算位移、内力的问题如何应用？ 位移法的做法 Δ2 F2 F1=0 F2=0 解: l EI P l l EI EI 2EI F1 Δ1 P k21 k11 Δ1=1 Δ2=1 k22 k12 F2P F1P P 3i/l 12i/l 12i/l 3i/l M1 8i 4i 3i M2 MP §10-6 等效结点荷载 {F}= [K]{?} ………………(1) 结构体系刚度方程： 一、位移法基本方程 [K]{?} +{FP} ={0} …………...………(2) {F} +{FP} ={0} …………..………(3) 将(1)式代入(2)式： 表示结点位移{?}和结点力{F}之间的关系，反映了结构的刚度性质，而不涉及原结构上作用的实际荷载，并不是原结构的位移法基本方程。 基本体系在荷载单独作用下产生的结点约束力。 基本体系在结点位移单独作用下产生的结点约束力。 k11? 1+ k12? 2+ · · · · · · · · · ·+ k1n? n+F1P=0 k21? 1+ k22? 2 +· · · · · · · · · ·+ k2n? n+F2P=0 · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · kn1? 1+ kn2? 2+ · · · · · · · · · ·+ knn? n+FnP=0 二、 等效结点荷载的概念 结点结束力——{FP} 结点结束力——{FP} 等效结点荷载{P} 原荷载 显然 {P}= –{FP}………解决了计算等效结点荷载的问题 等效原则是两种荷载在基本体系中产生相同的结点约束力； 也是结点位移等效。 [K]{?} = {F} {FP} + = 已知条件:固端内力 表10-1 三、按单元集成法求整体结构的等效结点荷载{P} (1)局部座标单元的等效结点荷载{P}， 表10-1 e x e e {P} e e (2)整体座标单元的等效结点荷载{P} e e e (3) 结构的等效结点荷载{P} x y 1 1 1 2 x y 1 2 3 4 8kN 4.8kN/m A B C 5m 2.5m 2.5m 单元1: 单元2: 1 2 12 10 -10 +4 +0 -5 1 2 2 2 2 思考：如果有作用在结点上的载荷如何计算？ 结构杆端内力的计算： e e e e e e e e e e e e {F} =[k] {?} e e [K]{?} = {F} {FP} + = 作业 ： 计算习题10-3， 1、整理原始数据，对单元和结构进行局部和总体编码 §10-7 计算步骤和算例 2、形成局部坐标系中的单元刚度矩阵 3、形成整体坐标系中的单元刚度矩阵 4、集成总刚度矩阵 5、形成结构等效结点荷载 6、解总刚度方程，得结点位移 7、求杆端力 建坐标系，结点坐标，杆件E、A、I等参数，编码，定位向量， 转换矩阵 结点载荷，加单元非结点载荷的等效载荷 与结点位移对应的杆端力，加固端力。 [K] 求单元常数 {?} [T] {P} 原始数据、局部码、总码 解方程[K]{?}={P} 求出结点位移 {?} 开始 单元刚度 矩阵 k e 单元固 端力 e 结束 [K]{?} = {F} {FP} + = 程序设计框图 求杆端力 e e e e §10-7 计算步骤和算例 { for(ne=1;ne=dyzs;ne++) { dw(dwxl,dyjd,jdwy,ne) ; /*调用定位向量函数*/ } dg(ne,dyjd,xy,eai,dygd); /*调用计算单刚*/ /*下面组集总刚度矩阵*/ for(j=1;j=6;j++) { zh=dwxl[j]; for(k=1;k=6;k++) { zl=dwxl[k];