第十三章矩阵位移法

1、131443523v基本概念v单元刚度矩阵(局部坐标系)(整体坐标系)v整体刚度矩阵(连续梁)(平面刚架)v等效结点荷载v计算步骤和算例v忽略轴向变形的矩形刚架的整体分析v上 机 作 业 （ 连 续 梁 程 序 设 计 ）213-1 概述 矩阵位移法以传统的结构力学作为理论基础，以矩阵作为数学表达形式，以电子计算机作为计算手段，三位一体的方法。 手算与电算的不同：手算：怕繁，讨厌重复性的大量运算，追求机灵的计算技巧， 运算次数较少的方法。电算：怕乱，讨厌头绪太多，零敲碎打的算法，追求计算过 程程序化，通用性强的方法。 矩阵位移法（有限单元法）的基本思路是： 先将结构离散成有限个单元，然后再将这

2、些单元按一定条件集合成整体。这样，就使一个复杂结构的计算问题转化为有限个简单单元的分析与集成问题。有限单元法的两个基本环节：1）单元分析：建立单元刚度方程，形成单元刚度矩阵(物理关系)2）整体分析：由单元刚度矩阵形集成整体刚度矩阵，建立结构的 位移法基本方程(几何关系、平衡条件)3 单元刚度矩阵是用来表示杆端力与杆端位移之间的物理关系的，不是新东西，但有几点新考虑：重新规定正负规则，以矩阵的形式表示，讨论杆件单元的一般情况。杆端局部编码与局部坐标系杆端局部编码与局部坐标系eE，A，Il局部坐标系中的杆端位移分量1u2u1v2v2q1q2M2Y2X1X1M1Y局部坐标系中的杆端力分量 = F22

3、2111MYXMYXee222111qqvuvue= De12yx13-2 单元刚度矩阵（局部坐标系）4单元刚度方程单元刚度方程FD方程1q 1u1v2v2q 2u1X1Y1M2X2Y2M)(,1221uulNXNX-=D=-=llEAND=)(),(212211uulEAXuulEAX-=-=由虎克定律：由转角位移方程，并考虑：2QYBA=1,QYAB-=12,vv -=D()212212212)(6vvlEIlEIY-+-=qq()212212112)(6vvlEIlEIY-+=qq()212212642vvlEIlEIlEIM-+=qq()212211624vvlEIlEIlEIM-+=

4、qq5222111MYXMYXe222111qqvuvue-lEIlEIlEIlEIlEIlEIlEIlEIlEAlEAlEIlEIlEIlEIlEIlEIlEIlEIlEAlEA460260612061200000260460612061200000222323222323= kee（135）单元刚度方程（136）单元刚度矩阵单元刚度矩阵的性质单元刚度矩阵的性质 1）单元刚度矩阵是杆端力用杆端位移来表达的联系矩阵。 2）其中每个元素称为单元刚度系数，表示由于单位杆端位移引起的杆端力。63k如 第 个杆端位移分量 =1时引起的第 个杆端力2M1q三六ijkji jDjiijkk=?反力互等定理

5、212221121121DD=kkkkFFD= kF6 3）单元刚度矩阵是对称矩阵。 4）第k列元素分别表示当第k个杆端位移=1时引起的六个杆端力分量。 5）一般单元的单元刚度矩阵是奇异矩阵。 不存在逆矩阵0=k De Fe正问题力学模型将单元视为“两端有六个人工控制的附加约束的杆件” De控制附加约束加以指定。解的性质为任何值时， De Fe 都有唯一的解答。且总是一个平衡力系，不可能是不平 衡力系。 De Fe反问题将单元视为“两端自由的杆件”。 Fe直接加在自由端作为指定的杆端力 Fe 为不平衡力系时 没有静力解。 De Fe 为平衡力系时 有无穷多组解。 De1X1Y1M2X2Y2M1

6、X1Y1M2X2Y2MD= kF72v=10312lEI312lEI-EI26l-EI26l-01v=10312lEI312lEI-EI26lEI26l0-=lEIlEIlEIlEIlEIlEIlEIlEIlEAlEAlEIlEIlEIlEIlEIlEIlEIlEIlEAlEAk460260612061200000260460612061200000222323222323e 第二列元素第二列元素变符号即第五列变符号即第五列 第一列元素第一列元素变符号即第四列变符号即第四列 第二行元素第二行元素变符号即第五行变符号即第五行 第一行元素第一行元素变符号即第四行变符号即第四行8特殊单元特殊单元单元

7、的某个或某些杆端位移的值已知为零。如梁单元、柱单元。特殊单元的单元刚度矩阵，可由一般单元的单元刚度矩阵删除 与零杆端位移对应的行和列得到。02211=vuvu-lEIlEIlEIlEIlEIlEIlEIlEIlEAlEAlEIlEIlEIlEIlEIlEIlEIlEIlEAlEA460260612061200000260460612061200000222323222323 =lEIlEIlEIlEIk4224为了使计算过程程序化、标准化、自动化，只采用一般单元 的刚度矩阵作为标准形式。各种特殊单元的刚度矩阵有计算 机程序去自动形成。某些特殊单元的刚度矩阵是可逆的。121q2q122M1M9选

8、局部坐标系推导单元刚度矩阵方便且单元刚度矩阵的形式简单。选整体坐标系是为进行整体分析。按一个统一的坐标系来建立各 单元的刚度矩阵单元坐标转换矩阵单元坐标转换矩阵1X1Y1M2X2Y2Myx1X1Y1M2X2Y2Myxyxyx1111111111111111cossinsincoscossinsincosMMYXYYXXMMYXYYXX=+-=+=+-=+=13-3 单元刚度矩阵（整体坐标系）10=222111MYXMYX222111MYXMYX-1000000cossin0000sincos0000001000000cossin0000sincos =T单元坐标转换矩阵T是一正交矩阵。TTT1

9、=-同理：整体坐标系中的单元刚度矩阵整体坐标系中的单元刚度矩阵D= kF设：将（a）、（b）代入（a）（b）D=TkFTD=TkTFTTTTFTF=FTFT=D=DTD=DTTD= kF111）ijk 表示在整体坐标系第j个杆端位移分量=1时引起的第i个杆端力分量。 2）k 是对称矩阵。3）一般单元的k是奇异矩阵。 例13-1 求图示刚架中各单元在整体标系中的单元刚度矩阵。设各杆的几何尺寸相同。l=5m，A=0.5m2, I=1/24m4E=3107kN/m2441025,10300=lEIlEA21与 比较 D= kFD=TkTFTTTT I k，k 同阶,性质类似：=00002211222

10、11211221122211211TTkkkkTTkkkkTT12解：（1）求kek1k2=503003012000300-1003003012000300503003012000300-1003003012000300-104（2）求 ke kkIT=,11单元 =9021单元 =0 -=100000001000010000000100000001000010Tk2=500300300030012-1000300300030012-500300300030012-100030030003012104-=100001010100300301200030010000101011111111TkT

11、kT13结点力、结点位移、形成总刚度矩阵结点力、结点位移、形成总刚度矩阵(传统位移法传统位移法)123F1F2F3123F1F2F31=11K11K21K31 K12 K22K32 2=12K13 K23K33 3=13333232131332322212123132121111D+D+D=D+D+D=D+D+D=KKKFKKKFKKKFF = KK 为整体刚度矩阵，简称总刚。13-4 连续梁的整体刚度矩阵14整体刚度矩阵的性质 1）总刚是结点力用结点位移来表达的联系矩阵。 2）K中的元素Kij表示第j个结点位移分量j=1（其它结点 位移分量=0）时所产生的第i个结点力。 3）K是对称矩阵。

12、4）如果引入支承条件，K是可逆矩阵。形成整体刚度矩阵形成整体刚度矩阵2=1K12 K22K32 112=D112k121k221211=D2k112k212结点发生单位位移杆端发生单位位移变形协调条件产生附加约束中约束力(总刚元素)产生杆端力(单刚元素)平衡条件总刚元素是由单刚元素集合而成K22k221k112k212K3215直接刚度法形成总刚(刚度集成法) 首先要注意同一个结点位移在整体中与在各单元中编码不同。单元结点位移总码按局部码顺序排列而成的向量称为“单元定位向量”。e单元对应关系：局部码总码单元定位向量e12A （1） 1B （2） 2 =112B （1） 2B （2） 3 =22

13、3 将各单元的单刚的行列局部码（i）、（j）换成对应的结点位移总码i、 j，按此行列总码将单刚元素送入总刚。即:k(i)(j)jiK2112213ABC(1)(2) (1)(2)16例13-2 试求图示连续梁的整体刚度矩阵K。i1i2i31230123解：1）编码凡给定 为零的结点位移分 量,其总码均编为零。 =112 =2232）单元定位向量 =3303）求单刚并集成总刚k =14i1 2i12i1 4i1(1) (2) 1 2 =K4i1 2i12i1 4i1k =24i2 2i22i2 4i2(1) (2) 2 3+ 4i2 2i22i2 4i2k =34i3 2i32i3 4i3(1)

14、 (2) 3 0+ 4i31 2 31 2300在给节点位移编码时已经考虑了支承条件。(先处理法)171n2312n+1对于n跨连续梁，有n+1个节点，不难导出整体刚度矩阵如下：4i12i12i14(i1+ i2)02i24(i1+ i2)02i22i3002in-14(In-1+ in) 2in2in4in000K=Kn+1，n+1是稀疏矩阵和带状矩阵。) 1, 3 , 2(2), 3 , 2(44,4,411111, 1111+ = =+=-+njiKKnjiiKiKiKjjjjjjjjjnnn1n2318情况复杂：1）结点位移分量增加到三个；2）各杆方向不尽相同，要进行坐标变换；3）除了

15、刚结点，还要考虑铰结点等其它情况。1、结点位移分量的统一编码总码yx000123040结点位移列阵:=1 2 3 4T =uA vA A CT结点力列阵:F=F1 F2 F3 F4T2、单元定位向量211(1)(2)(3)(4)(6) (5)2(1)(2)(3)(5)(4)(6) =11 2 3 0 0 4T =21 2 3 0 0 0TACB13-5 刚架的整体刚架矩阵193、单元集成过程503003012000300-1003003012000300503003012000300-1003003012000300-104k =1 1 2 3 0 0 4K=1 2 3 4300 000012

16、30100030100500305030104k2=500300300030012-1000300300030012-500300300030012-1000300300030121041 2 3 0 0 0 +12+030+0+300+030+0+100求单刚201）结点位移分量的统一编码总码 铰结点处的两杆端结点应看 作半独立的两个结点（C1和 C2） 它们的线位移相同， 角位移不同，00012321A C1B D000456475C234、铰结点的处理线位移采用同码，角位移采用异码。2）单元定位向量： =11 2 3 4 5 6T =21 2 3 0 0 0T =34 5 7 0 0 0

17、T3）按次序进行单元集成：21503003012000300-1003003012000300503003012000300-1003003012000300-104k =1 1 2 3 4 5 61 2 3 4 5 6 7300500-3010030000-30000 01230 0-12300301000-30500-12-30012-30K=104 -30000300001 2 3 0 0 0k2=500300300030012-1000300300030012-500300300030012-100030030003012104+12+030+0+300+030+0+100k3=500

18、300300030012-1000300300030012-500300300030012-1000300300030121044 5 7 0 0 0 +12+030+0+300+030+0+10030?221、整体刚度方程、整体刚度方程 F=K （a） 表示由表示由F结点力结点力的关系式。反映了结点的刚度性质，的关系式。反映了结点的刚度性质， 不涉及结构上的实际荷载。不涉及结构上的实际荷载。2、位移法基本方程、位移法基本方程 在给结点位移分量编总码时，已考虑了结构的支承连接情况，在给结点位移分量编总码时，已考虑了结构的支承连接情况， K是非奇异矩阵。是非奇异矩阵。 如果已知结构上的结点荷载如

19、果已知结构上的结点荷载P，（，（a）就是求结点位移）就是求结点位移 的位移法基本方程的位移法基本方程。P=K（b）注：结点力与结点荷载的不同。结点力是发生给定的结点位移，注：结点力与结点荷载的不同。结点力是发生给定的结点位移， 在结点上所需施加的力，它与体系的刚度有关，由刚度方在结点上所需施加的力，它与体系的刚度有关，由刚度方 程确定。而结点荷载是给定的与体系无关。由结点荷载程确定。而结点荷载是给定的与体系无关。由结点荷载 产产 生的未知结点位移由位移法基本方程求解。生的未知结点位移由位移法基本方程求解。3、等效结点荷载、等效结点荷载 平衡方程的荷载平衡方程的荷载P是作用在结点上的集中荷载，当

20、荷载不是作用在结点上的集中荷载，当荷载不是结点集中荷载时，应化成等效结点荷载。是结点集中荷载时，应化成等效结点荷载。13-6 13-6 等效结点荷载等效结点荷载23123FP2FP1FP3P2P1P3结点约束力FP =-FP2 =-FP1-FP3=1323=2=1等效结点荷载P=FP可由位移法基本方程可由位移法基本方程(b)求得求得.注意：注意：j非结点荷载与等效结点荷载等效的条件是，两者产生相同非结点荷载与等效结点荷载等效的条件是，两者产生相同 结点位移。结点位移。j除了结点位移外，等效结点荷载与原荷载产生的其它位移除了结点位移外，等效结点荷载与原荷载产生的其它位移 和内力并不相同。和内力并

21、不相同。j等效结点荷载为位移法基本体系附加约束中约束力的负值。等效结点荷载为位移法基本体系附加约束中约束力的负值。 而约束力为各固端力之和。所以求结构等效结点荷载应该而约束力为各固端力之和。所以求结构等效结点荷载应该 先求出单元的等效结点荷载，它是单元固端力的负值。先求出单元的等效结点荷载，它是单元固端力的负值。0=+PFP位移法方程：位移法方程：D= KP244、按单元集成法求整体结构的等效结点荷载 局部坐标系中的单元固端约束力PFe整体坐标系中的单元等效结点荷载 PTFTP-=ee整体结构的等效结点荷载P 由各单元P 中的元素按在P中进行定位并累加。e等效结点荷载与直接结点荷载叠加，即得结

22、构的结点荷载。4.8kN/m2.5m2.5m5m8kN12yx例13-3 求图示结构的等效结点荷载P. PFe解:1)求单元,.10120111mkNMkNYXPPP-=-=mkNMkNYXPPP.10120222=-=单元,.540111mkNMkNYXPPP=mkNMkNYXPPP.540222-= TPF1012010120-= TPF540540-=252)求 PTFTP-=单元的倾角1=01 2 3 0 0 4单元的倾角2=90123000P=123401210-10+4 +05 -=-=-=5045045405401000000010000100000001000000010000

23、10PTFTP TPPTFIFTP101010120-=-=-=4 125-10 TPF1012010120-= TPF540540-=261）整理原始数据，进行局部编码和整体编码。2）用式（13-6）形成局部坐标系中的单元刚度矩阵3）用式（13-21）形成整体坐标系中的单元刚度矩阵4）用单元集成法形成整体刚度矩阵K5）形成整体结构的等效结点荷载6）解方程K=P， 求出结点位移。7）求各杆杆端力123FP2FP1FP3P2P1P3 =-FP2 =-FP1-FP3=1323=2=1固端力FP杆端位移产生的杆端力PkkFDkPFkF+D=计算步骤13-7 计算步骤和算例27000123000456

24、例13-4:求内力。横梁b1h1=0.5m 1.26m,立柱b2h2=0.5m 1m.6m12m1kN/m213xy.1031. 212,1094. 66,108 .274,109 .132,103 .83,1094. 6, 6,241, 5 . 0:33323333-=lEIlEIlEIlEIlEAlEIlIA柱.1058. 012,1047. 36,108 .274,109 .132,105 .52,1094. 6,12,121,63. 0:33323333-=lEIlEIlEIlEIlEAlEIlIA梁解:1)原始数据及编码28.1031. 212,1094. 66,108 .274,1

25、09 .132,103 .83,1094. 6, 6,241, 5 . 0:33323333-=lEIlEIlEIlEIlEAlEIlIA柱-=lEIlEIlEIlEIlEIlEIlEIlEIlEAlEAlEIlEIlEIlEIlEIlEIlEIlEIlEAlEAk460260612061200000260460612061200000222323222323e=8 .2794. 6094. 631. 20003 .83kk139 .1394. 6094. 631. 20003 .83-8 .2794. 609 .1394. 6094. 631. 2094. 631. 20003 .83003

26、 .83-.1058. 012,1047. 36,108 .274,109 .132,105 .52,1094. 6,12,121,63. 0:33323333-=lEIlEIlEIlEIlEAlEIlIA梁103-=8 .2747. 3047. 358. 00005 .529 .1347. 3047. 358. 00005 .529 .1347. 3047. 358. 00005 .528 .2747. 3047. 358. 00005 .52k21032)形成k293) 形成k 单元、(=90o)坐标转换矩阵为 -=100000001000010000000100000001000010T

27、-=1000010108 .2794. 6094. 631. 20003 .8310000101011111111TkTkT-=8 .27094. 603 .83094. 6031. 2kk131039 .13094. 603 .83094. 6031. 2-8 .27094. 69 .13094. 603 .83003 .83094. 6031. 294. 6031. 2-单元(=0o)坐标转换矩阵为单位矩阵所以:kk=224) 形成K=00032112=654321=0006543+11k+11k=22211211kkkkK30-=6 .5547. 394. 647. 388.83094.

28、 6081.549 .1347. 3047. 358. 00005 .529 .1347. 3047. 358. 00005 .526 .5547. 394. 647. 388.83094. 6081.54K1035) 求等效节点荷载P=3-30330F:P固端力1=-=3033-03PPTFT11单元在整体坐标系中的等效节点荷载000321集 成 等效 节 点荷载=0003-03P316) 解基本方程-=DDDDDD-0003036 .5547. 394. 647. 388.83094. 6081.549 .1347. 3047. 358. 00005 .529 .1347. 3047. 3

29、58. 00005 .526 .5547. 394. 647. 388.83094. 6081.54106543213-=DDDDDD5 .9613. 58244 .2813. 584765432132-=D5 .9613. 58244 .2813. 58477) 求杆端力 单元T0004 .2813. 5847-=D1=+D= PFTkF-=-+-49. 876. 443. 009. 224. 143. 03303300004 .2813. 58471000000010000100000001000000010000108 .2794. 6094. 631. 20003 .839 .1394

30、. 6094. 631. 20003 .839 .1394. 6094. 631. 20003 .838 .2794. 6094. 631. 20003 .83103=+D= PFTkF-=-04. 343. 024. 109. 243. 024. 15 .9613. 58244 .2813. 58478 .2747. 309 .1347. 3047. 358. 0047. 358. 00005 .52005 .529 .1347. 308 .2747. 3047. 358. 0047. 358. 00005 .52005 .5210333=D0005 .9613. 5824=+D= PFTk

31、F-=-38. 424. 143. 004. 324. 143. 00005 .9613. 58241000000010000100000001000000010000108 .2794. 6094. 631. 20003 .839 .1394. 6094. 631. 20003 .839 .1394. 6094. 631. 20003 .838 .2794. 6094. 631. 20003 .83103T38. 424. 143. 004. 324. 143. 0-TF04. 343. 024. 109. 243. 024. 1 -=TF49. 876. 443. 009. 224. 14

32、3. 0-=8.492.093.044.38M图(kN.m)4.761.240.43 1.241.24Q图(kN)N=0.43N=1.24N=0.43N图(kN)341）结点位移分量的统一编码总码 在刚结点A铰结点C1和C2处， 竖向位移均为零，故其编码 也应为零，另外它们的水平 位移分量都相等，因此它们 的水平位移应采用同码。00010221A C1B D000103140C232）单元定位向量： =11 0 2 1 0 3T =21 0 2 0 0 0T =31 0 4 0 0 0T3）按次序进行单元集成：13-8 忽略轴向变形时矩形刚架的整体分析35503003012000300-100

33、3003012000300503003012000300-1003003012000300-104k =1 1 0 2 1 0 31 2 3 4 1234102103300030001000 +050300+0+300+0050+01001 0 2 0 0 0k2=500300300030012-1000300300030012-500300300030012-1000300300030121041020000 0 00 100 500 50 100+123030+100k3=500300300030012-1000300300030012-500300300030012-1000300300

34、030121041 0 4 0 0 0104000+123030+100K=104 36例13-5:求内力。横梁b1h1=0.5m 1.26m,立柱b2h2=0.5m 1m 忽略轴向变形的影响。.6m12m1kN/m.1031. 212,1094. 66,108 .274,109 .132,103 .83,1094. 6, 6,241, 5 . 0:33323333-=lEIlEIlEIlEIlEAlEIlIA柱.1058. 012,1047. 36,108 .274,109 .132,105 .52,1094. 6,12,121,63. 0:33323333-=lEIlEIlEIlEIlEA

35、lEIlIA梁解:1)原始数据及编码000102000103213xy37.1031. 212,1094. 66,108 .274,109 .132,103 .83,1094. 6, 6,241, 5 . 0:33323333-=lEIlEIlEIlEIlEAlEIlIA柱-=lEIlEIlEIlEIlEIlEIlEIlEIlEAlEAlEIlEIlEIlEIlEIlEIlEIlEIlEAlEAk460260612061200000260460612061200000222323222323e=8 .2794. 6094. 631. 20003 .83kk139 .1394. 6094. 63

36、1. 20003 .83-8 .2794. 609 .1394. 6094. 631. 2094. 631. 20003 .83003 .83-.1058. 012,1047. 36,108 .274,109 .132,105 .52,1094. 6,12,121,63. 0:33323333-=lEIlEIlEIlEIlEAlEIlIA梁103-=8 .2747. 3047. 358. 00005 .529 .1347. 3047. 358. 00005 .529 .1347. 3047. 358. 00005 .528 .2747. 3047. 358. 00005 .52k21032)形

37、成k383) 形成k 单元、(=90o)坐标转换矩阵为 -=100000001000010000000100000001000010T-=1000010108 .2794. 6094. 631. 20003 .8310000101011111111TkTkT-=8 .27094. 603 .83094. 6031. 2kk131039 .13094. 603 .83094. 6031. 2-8 .27094. 69 .13094. 603 .83003 .83094. 6031. 294. 6031. 2-单元(=0o)坐标转换矩阵为单位矩阵所以:kk=224) 形成K=00020112=30

38、1201=000301339 1 0 2 0 0 0-=-8 .27094. 69 .13094. 603 .83003 .83094. 6031. 294. 6031. 29 .13094. 68 .27094. 603 .83003 .83094. 6031. 294. 6031. 2103k1020002.316.9427.86.94-=-8 .27094. 69 .13094. 603 .83003 .83094. 6031. 294. 6031. 29 .13094. 68 .27094. 603 .83003 .83094. 6031. 294. 6031. 2103k 1 0 3

39、 0 0 0103000+2.316.946.9427.8-=-8 .2747. 309 .1347. 3047. 358. 0047. 358. 00005 .52005 .529 .1347. 308 .2747. 3047. 358. 0047. 358. 00005 .52005 .52103k 1 0 2 1 0 3102103+52.5 52.5+0+0+52.552.5+0+27.8+013.9+0+0+013.9+27.8+0 4.62 6.94 6.946.94 55.6 13.96.94 13.9 55.61 2 3 123K=103405) 求等效节点荷载P=3-3033

40、0F:P固端力1=-=3033-03PPTFT11单元在整体坐标系中的等效节点荷载000201集 成 等效 节 点荷载=03-3P6) 解基本方程=-=-9 .971 .26838:0336 .559 .1394. 69 .136 .5594. 694. 694. 662. 4103BAABAAuuqqqq解得41=D9 .9708381 .2608387) 求杆端力 单元T0001 .260838=D1=+D= PFTkF-=-+-41. 875. 4009. 225. 103303300001 .2608381000000010000100000001000000010000108 .27

41、94. 6094. 631. 20003 .839 .1394. 6094. 631. 20003 .839 .1394. 6094. 631. 20003 .838 .2794. 6094. 631. 20003 .83103=+D= PFTkF-=-09. 343. 0009. 243. 009 .9708381 .2608388 .2747. 309 .1347. 3047. 358. 0047. 358. 00005 .52005 .529 .1347. 308 .2747. 3047. 358. 0047. 358. 00005 .52005 .5210342=D0009 .9708

42、38=+D= PFTkF-=-47. 425. 1009. 325. 100009 .9708381000000010000100000001000000010000108 .2794. 6094. 631. 20003 .839 .1394. 6094. 631. 20003 .839 .1394. 6094. 631. 20003 .838 .2794. 6094. 631. 20003 .83103T47. 425. 1009. 325. 10-TF09. 343. 0009. 243. 00-=TF41. 875. 4009. 225. 10-=8.412.093.094.47M图(k

43、N.m)4.751.250.43 1.251.25Q图(kN)N=0.43N=1.25N=0.43N图(kN) 由单元刚度方程求出的杆端轴 力为零。为什么？ 根据节点平衡由剪力求轴力。TF04. 343. 024. 109. 243. 024. 1 -=TF49. 876. 443. 009. 224. 143. 0-=TF38. 424. 143. 004. 324. 143. 0-= 轴向变形影响不大。43单元的刚度方程(局部坐标）u1u2X1X2e12)(),(212211uulEAXuulEAX-=-=-=21211111uulEAXXX1X2Y2Y1X1X2yxx-=22112211

44、0000010100000101vuvulEAYXYX -=cossin00sincos0000cossin00sincosT注意：桁架单元的结点转角不是基本未知量。无须求等效结点荷载。杆端力是由结点位移产生的。D=TkF13-9 桁架及组合结构的整体分析坐标转换矩阵单元的刚度方程(整体坐标）44 3ll 10kN10kN例13-6 求图示桁架内力(EA=常数)。解：1、编码如图 ；124002、形成k6 , 5,000001010000010124 , 3 , 2 , 1,0000010100000101=-=-=lEAklEAk3、形成kk=k=k=k单元=90-=011011Tk=k=l

45、EA-=01100001011011111111lEATkTkT1000-1000-1000100045单元=45-=11112111Tk=lEA22- -=1111000111112211111111lEATkTkT1111-11111111-1111单元=135-=11112111Tk=lEA22-=1111000111112211111111lEATkTkT1111-1111-1111-1111- 3ll 10kN10kN124004、集成总刚KTT4321 ,0021 =TT00000043=TT00430021 =lEAK-=35. 135. 00035. 035. 1010035.

46、 135. 00135. 035. 15、节点荷载T0010-10=P46-=-00101035. 135. 00035. 035. 1010035. 135. 00135. 035. 1DDCCvuvulEA6、解基本方程lEAvuvuDDCC-=58. 536.2142.1494.267、杆端力计算=D.-=-=D=042.14042.140042.1494.2601001000000100100000010100000101TkF-=-=D=058. 5058. 558. 536.2142.1494.2610000100001000010000010100000101TkF0021lEA

47、-=D0042.1494.26 4321.=D lEA-=D58. 536.2142.1494.2647节点位移分量自由节点位移分量(基本未知量，相应的节点荷载已知)受约束的位移分量(已知量，相应的约束反力未知)先处理法1、节点位移分量中不含受约束的支座位移，节点力分量中不 含未知的支座反力。2、由单刚考虑边界条件K3、对于具有非刚性连接、支承节点较多且分散、不考虑轴向 变形的结构最为方便。可减少内存，提高计算速度。但要 对各节位移进行统一编码，形成各单元的定位向量。后处理法1、节点位移分量中含有受约束的支座位移，节点力分量中含 有未知的支座反力。2、由单刚考虑边界条件KK(原始刚度矩阵,奇异

48、)3、每个节点位移分量数相同， 的阶数是节点总数乘节点位 移 分量数，整个分析过程便于编制通用程序。适用于节点 多支座约束少，考虑轴向变形的结构。但占用内存大。K48后处理法边界条件的处理 在后处理法中，由于没有考虑边界条件，由k集成的 是奇异矩阵，由单元集合成的体系是自由体，具有刚体位移。 KPK=D没有确定的位移解。位移边界条件处理的三种方法：1、划行划列法 编制程序较复杂，不常采用。2、主对角元置大数法设第 i 个节点位移分量 (已知)0di=D =D D DD nininnninniniiiininiPPPPkkkkkkkkkkkkkkkk212121212222211112110di=D为了将第 i 个方程改为:将 kii 置一大数如R=1020 Pi 改为 Rd0第 i 个方程变为:00112211dRdkRkkiiiii=D=D+ +D+ +D+D该法虽为近似处理，程序设计容易实现，故被广泛应用。该法虽为近似处理，程序设计容易实现，故被广泛应用。RRd0492、主对角元置1法 =D