矩阵位移法 结构力学 学习资料

1、1郑州大学土木工程学院 樊友景编制 基本要求：熟练掌握两种坐标系中的单元刚度矩阵、 结构的整体刚度矩阵、等效结点荷载 的形成，已知结点位移求单元杆端力 的计算方法，整体刚度矩阵和结构结 点荷载的集成过程。理解单元刚度矩阵中和整体刚度矩阵中的 元素的物理意义。了解不计轴向变形时矩形刚架的整体分析.MATRIX DISPLACEMENT METHOD2郑州大学土木工程学院 樊友景编制13-1 概述 矩阵位移法以传统的结构力学作为理论基础，以矩阵作为数学表达形式，以电子计算机作为计算手段，三位一体的方法。 手算与电算的不同：手算：怕繁，讨厌重复性的大量运算，追求机灵的计算技巧， 运算次数较少的方法。

2、电算：怕乱，讨厌头绪太多，零敲碎打的算法，追求计算过 程程序化，通用性强的方法。矩阵位移法（有限单元法finite element method）的基本思路是： 先将结构离散成有限个单元，然后再将这些单元按一定条件集合成整体。这样，就使一个复杂结构的计算问题转化为有限个简单单元的分析与集成问题。有限单元法的两个基本环节：1）单元分析：建立单元刚度方程，形成单元刚度矩阵(物理关系)2）整体分析：由单元刚度矩阵形成整体刚度矩阵，建立结构的 位移法基本方程(几何关系、平衡条件)3郑州大学土木工程学院 樊友景编制 单元刚度矩阵是用来表示杆端力与杆端位移之间的物理关系的，不是新东西，但有几点新考虑：重新

3、规定正负规则，以矩阵的形式表示，讨论杆件单元的一般情况。杆端局部编码与局部坐标系杆端局部编码与局部坐标系eE，A，Il局部坐标系中的杆端位移分量1u2u1v2v2q1q2M2Y2X1X1M1Y局部坐标系中的杆端力分量 = F222111MYXMYXee222111qqvuvue= De12yx13-2 单元刚度矩阵（局部坐标系） (element stiffness matrix)4郑州大学土木工程学院 樊友景编制单元刚度方程单元刚度方程FD方程1q 1u1v2v2q 2u1X1Y1M2X2Y2M)(,1221uulNXNX-=D=-=llEAND=)(),(212211uulEAXuulEA

4、X-=-=由虎克定律：由转角位移方程，并考虑：2QYBA=1,QYAB-=12,vv -=D()212212212)(6vvlEIlEIY-+-=qq()212212112)(6vvlEIlEIY-+=qq()212212642vvlEIlEIlEIM-+=qq()212211624vvlEIlEIlEIM-+=qq5郑州大学土木工程学院 樊友景编制222111MYXMYXe222111qqvuvue-lEIlEIlEIlEIlEIlEIlEIlEIlEAlEAlEIlEIlEIlEIlEIlEIlEIlEIlEAlEA4602606120612000002604606120612000002

5、22323222323= kee（135）单元刚度方程（136）单元刚度矩阵单元刚度矩阵的性质单元刚度矩阵的性质 1）单元刚度矩阵是杆端力用杆端位移来表达的联系矩阵。 2）其中每个元素称为单元刚度系数，表示由于单位杆端位移引起的杆端力。63k如 第 个杆端位移分量 =1时引起的第 个杆端力2M1q三六ijkji jDjiijkk=?反力互等定理212221121121DD=kkkkFFD= kF6郑州大学土木工程学院 樊友景编制 3）单元刚度矩阵是对称矩阵。 4）第k列元素分别表示当第k个杆端位移=1时引起的六个杆端力分量。 5）一般单元的单元刚度矩阵是奇异矩阵。 不存在逆矩阵0=k De F

6、e正问题力学模型将单元视为“两端有六个人工控制的附加约束的杆件” De控制附加约束加以指定。解的性质为任何值时， De Fe 都有唯一的解答。且总是一个平衡力系，不可能是不平 衡力系。 De Fe反问题将单元视为“两端自由的杆件”。 Fe直接加在自由端作为指定的杆端力 Fe 为不平衡力系时 没有静力解。 De Fe 为平衡力系时 有无穷多组解。 De1X1Y1M2X2Y2M1X1Y1M2X2Y2MD= kF7郑州大学土木工程学院 樊友景编制2v=10312lEI312lEI-EI26l-EI26l-01v=10312lEI312lEI-EI26lEI26l0-=lEIlEIlEIlEIlEIl

7、EIlEIlEIlEAlEAlEIlEIlEIlEIlEIlEIlEIlEIlEAlEAk460260612061200000260460612061200000222323222323e 第二列元素第二列元素变符号即第五列变符号即第五列 第一列元素第一列元素变符号即第四列变符号即第四列 第二行元素第二行元素变符号即第五行变符号即第五行 第一行元素第一行元素变符号即第四行变符号即第四行8郑州大学土木工程学院 樊友景编制特殊单元特殊单元单元的某个或某些杆端位移的值已知为零。如梁单元、柱单元。特殊单元的单元刚度矩阵，可由一般单元的单元刚度矩阵删除 与零杆端位移对应的行和列得到。02211=vuvu

8、-lEIlEIlEIlEIlEIlEIlEIlEIlEAlEAlEIlEIlEIlEIlEIlEIlEIlEIlEAlEA460260612061200000260460612061200000222323222323 =lEIlEIlEIlEIk4224为了使计算过程程序化、标准化、自动化，只采用一般单元 的刚度矩阵作为标准形式。各种特殊单元的刚度矩阵有计算 机程序去自动形成。某些特殊单元的刚度矩阵是可逆的。121q2q122M1M9郑州大学土木工程学院 樊友景编制选局部坐标系推导单元刚度矩阵方便且单元刚度矩阵的形式简单。选整体坐标系是为进行整体分析。按一个统一的坐标系来建立各 单元的刚度矩

9、阵单元坐标转换矩阵单元坐标转换矩阵1X1Y1M2X2Y2Myx1X1Y1M2X2Y2Myxyxyx11111111cossinsincosMMYXYYXX=+-=+=13-3 单元刚度矩阵（整体坐标系）22222222cossinsincosMMYXYYXX=+-=+=局部坐标系中的杆端力整体坐标系中的杆端力10郑州大学土木工程学院 樊友景编制=222111MYXMYX222111MYXMYX-1000000cossin0000sincos0000001000000cossin0000sincos单元坐标转换矩阵T是一正交矩阵。TTT1=-同理：整体坐标系中的单元刚度矩阵整体坐标系中的单元刚度

10、矩阵D= kF设：将（a）、（b）代入（a）（b）D=TkFTD=TkTFTTTTFTF=FTFT=D=DTD=DTTD= kF =T11郑州大学土木工程学院 樊友景编制1）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 TkTkT=k，k 同阶,性质类似：=000022112221121

11、1221122211211TTkkkkTTkkkkTT12郑州大学土木工程学院 樊友景编制解：（1）求kek1k2=503003012000300-1003003012000300503003012000300-1003003012000300-104（2）求 ke kkIT=,11单元 =9021单元 =0 -=100000001000010000000100000001000010Tk2=500300300030012-1000300300030012-500300300030012-100030030003012104-=100001010100300301200030010000101

12、011111111TkTkT13郑州大学土木工程学院 樊友景编制结点力、结点位移、形成总刚度矩阵结点力、结点位移、形成总刚度矩阵(传统位移法传统位移法)123F1F2F3123F1F2F31=11K11K21K31 K12 K22K32 2=12K13 K23K33 3=13333232131332322212123132121111D+D+D=D+D+D=D+D+D=KKKFKKKFKKKFDDD=321333231232221131211321KKKKKKKKKFFFF = KK 为整体刚度矩阵，简称总刚。 =333231232221131211KKKKKKKKKK13-4 连续梁的整体分

13、析assembly analysis of continuous beam14郑州大学土木工程学院 樊友景编制整体刚度矩阵(assembled stiffniss matrix)的性质 1）总刚是结点力用结点位移来表达的联系矩阵。 2）K中的元素Kij表示第j个结点位移分量j=1（其它结点 位移分量=0）时所产生的第i个结点力。 3）K是对称矩阵。 4）如果引入支承条件，K是可逆矩阵。形成整体刚度矩阵形成整体刚度矩阵2=1K12 K22K32 112=D112k121k221211=D2k112k212结点发生单位位移杆端发生单位位移变形协调条件产生附加约束中约束力(总刚元素)产生杆端力(单刚

14、元素)平衡条件总刚元素是由单刚元素集合而成K22k221k112k212K32 综上所述，直接刚度法是根据单元的结点位移分量的直接刚度法是根据单元的结点位移分量的局部码和总码之间的对应关系，由单元刚度矩阵集成结构局部码和总码之间的对应关系，由单元刚度矩阵集成结构整体刚度矩阵。整体刚度矩阵。15郑州大学土木工程学院 樊友景编制直接刚度法形成总刚(刚度集成法) 首先要注意同一个结点位移在整体中与在各单元中两种编码不同。单元结点位移总码按局部码顺序排列而成的向量称为“单元定位向量” (element localization vector) 。e单元对应关系：局部码总码单元定位向量e12A （1）

15、1B （2） 2 =112B （1） 2C （2） 3 =223 将各单元的单刚的行列局部码（i）、（j）换成对应的结点位移总码i、 j，按此行列总码将单刚元素送入总刚。即:k(i)(j)jiK2112213ABC(1)(2) (1)(2)整体分析用总码，单元分析用局部码整体分析用总码，单元分析用局部码16郑州大学土木工程学院 樊友景编制例13-2 试求图示连续梁的整体刚度矩阵K。i1i2i31230123解：1）编码凡给定 为零的结点位移分 量,其总码均编为零。 =112 =2232）单元定位向量 =3303）求单刚并集成总刚(1) (2)1 2 =K1 2 31 23在给节点位移编码时已经

16、考虑了支承条件。(先处理法) k=4i12i12i14i14i12i12i14i1 k=4i22i22i24i2(1) (2)2 34i22i22i24i2+ k=4i32i32i34i3(1) (2)3 04i3+17郑州大学土木工程学院 樊友景编制1n2312n+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+ = =+=-

17、+njiKKnjiiKiKiKjjjjjjjjjnnn1n2318郑州大学土木工程学院 樊友景编制1）结点位移分量增加到三个；2）各杆方向不尽相同，要进行坐标变换；3）除了刚结点，还要考虑铰结点等其它情况。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 刚架的整体刚度矩阵(assembly analysis of frame

18、)情况复杂：19郑州大学土木工程学院 樊友景编制3、单元集成过程K=1 2 3 4104503003012000300-1003003012000300503003012000300-1003003012000300-104 1 2 3 0 0 430012301003010050305030k =11 2 3 0 0 0 k2=500300300030012-1000300300030012-500300300030012-100030030003012104+1230+30030+10020郑州大学土木工程学院 樊友景编制1）结点位移分量的统一编码总码 铰结点处的两杆端结点应看 作半独立的

19、两个结点（C1和 C2） 它们的线位移相同， 角位移不同，00012321A C1B D000456475C234、铰结点、铰结点(hinge joint)的处理的处理线位移采用同码，角位移采用异码。2）单元定位向量： =11 2 3 4 5 6T =21 2 3 0 0 0T =34 5 7 0 0 0T3）按次序进行单元集成：21郑州大学土木工程学院 樊友景编制503003012000300-1003003012000300503003012000300-1003003012000300-104k =1 1 2 3 4 5 61 2 3 4 5 6 7300500-3010030000-3

20、0000 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=500300300030012-1000300300030012-500300300030012-1000300300030121044 5 7 0 0 0 +12+030+0+300+030+0+10030?22郑州大学土木工程学院 樊友景编制1、整体刚度方程、整体

21、刚度方程 F=K （a） 表示由表示由F结点力的关系式。反映了结构的刚度性质，结点力的关系式。反映了结构的刚度性质， 不涉及结构上的实际荷载。不涉及结构上的实际荷载。2、位移法基本方程、位移法基本方程 在给结点位移分量编总码时，已考虑了结构的支承连接情况，在给结点位移分量编总码时，已考虑了结构的支承连接情况， K是非奇异矩阵。是非奇异矩阵。 如果已知结构上的结点荷载如果已知结构上的结点荷载P，（，（a）就是求结点位移）就是求结点位移 的位移法基本方程的位移法基本方程。P=K（b）注：结点力与结点荷载的不同。结点力是发生给定的结点位移，注：结点力与结点荷载的不同。结点力是发生给定的结点位移， 在

22、结点上所需施加的力，它与体系的刚度有关，由刚度方在结点上所需施加的力，它与体系的刚度有关，由刚度方 程确定。而结点荷载是给定的程确定。而结点荷载是给定的,与体系无关。由结点荷载与体系无关。由结点荷载 产产 生的未知结点位移由位移法基本方程求解。生的未知结点位移由位移法基本方程求解。3、等效结点荷载、等效结点荷载 平衡方程的荷载平衡方程的荷载P是作用在结点上的集中荷载，当荷载不是作用在结点上的集中荷载，当荷载不是结点集中荷载时，应化成等效结点荷载。是结点集中荷载时，应化成等效结点荷载。13-6 等效结点荷载(equivalent joint load)23郑州大学土木工程学院 樊友景编制123F

23、P2FP1FP3P2P1P3结点约束力FP =-FP2 =-FP1-FP3=1323=2=1等效结点荷载P=FPv除了结点位移外，等效结点荷载与原荷载产生的其它位移除了结点位移外，等效结点荷载与原荷载产生的其它位移 和内力并不相同。和内力并不相同。v等效结点荷载为位移法基本体系附加约束中约束力的负值。等效结点荷载为位移法基本体系附加约束中约束力的负值。 而约束力为各固端力之和。所以求结构等效结点荷载应该而约束力为各固端力之和。所以求结构等效结点荷载应该 先求出单元的等效结点荷载，它是单元固端力的负值。先求出单元的等效结点荷载，它是单元固端力的负值。0=+PFP位移法方程：位移法方程：由由 求得

24、求得D= KP注意：注意：v非结点荷载与等效结点荷载等效的条件是，两者产生相同非结点荷载与等效结点荷载等效的条件是，两者产生相同 结点位移。在基本体系上产生相同的结点约束力。结点位移。在基本体系上产生相同的结点约束力。24郑州大学土木工程学院 樊友景编制4、按单元集成法求整体结构的等效结点荷载、按单元集成法求整体结构的等效结点荷载 局部坐标系中的单元固端约束力PFe整体坐标系中的单元等效结点荷载 PTFTP-=ee整体结构的等效结点荷载P 由各单元P 中的元素按在P中进行定位并累加。e等效结点荷载与直接结点荷载叠加，即得结构的结点荷载。4.8kN/m2.5m2.5m5m8kN12yx例13-3

25、 求图示结构的等效结点荷载P. PFe解:1)求单元,.10120111mkNMkNYXPPP-=-=mkNMkNYXPPP.10120222=-=单元,.540111mkNMkNYXPPP=mkNMkNYXPPP.540222-= TPF1012010120-= TPF540540-=25郑州大学土木工程学院 樊友景编制2)求 PTFTP-=单元的倾角1=01 2 3 0 0 4单元的倾角2=90123000P=123401210-10+4 +05 -=-=-=504504540540100000001000010000000100000001000010PTFTP TPPTFIFTP101

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

27、制000123000456例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郑州大学土木工程学院 樊友景编制.

28、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. 631. 20003 .83-8 .2794. 609 .1394. 6094.

29、 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)形成k29郑州大学土木工程学院 樊友景编制3) 形成k 单元、(=90o)坐标

30、转换矩阵为 -=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=00032112=654321=0006543+11k+11k=2

31、2211211kkkkK30郑州大学土木工程学院 樊友景编制-=6 .5547. 394. 647. 388.83094. 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-03P31郑州大学土木工程学院 樊友景编制6) 解基本方程-=DDDDDD-0003036 .55

32、47. 394. 647. 388.83094. 6081.549 .1347. 3047. 358. 00005 .529 .1347. 3047. 358. 00005 .526 .5547. 394. 647. 388.83094. 6081.54106543213-=DDDDDD5 .9613. 58244 .2813. 584765432132郑州大学土木工程学院 樊友景编制=00032112=654321=0006543-=DDDDDD5 .9613. 58244 .2813. 5847654321-=D5 .9613. 58244 .2813. 58477) 求杆端力 单元T00

33、04 .2813. 5847-=D1=+D= PFTkF-=-+-49. 876. 443. 009. 224. 143. 03303300004 .2813. 58471000000010000100000001000000010000108 .2794. 6094. 631. 20003 .839 .1394. 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 .2

34、813. 58478 .2747. 309 .1347. 3047. 358. 0047. 358. 00005 .52005 .529 .1347. 308 .2747. 3047. 358. 0047. 358. 00005 .52005 .5210333郑州大学土木工程学院 樊友景编制=00032112=654321=0006543-=DDDDDD5 .9613. 58244 .2813. 5847654321=D0005 .9613. 5824=+D= PFTkF-=-38. 424. 143. 004. 324. 143. 00005 .9613. 58241000000010000

35、100000001000000010000108 .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. 143. 0-=8.492.093.044.38M图(kN.m)4.761.240.43 1.241.24Q图(kN)N=0.43

36、N=1.24N=0.43N图(kN)34郑州大学土木工程学院 樊友景编制1）结点位移分量的统一编码总码 在刚结点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 忽略轴向变形时矩形刚架的整体分析35郑州大学土木工程学院 樊友景编制503003012000300-1003003012000300503003012000300-10

37、03003012000300-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-1000300300030121041 0 4 0 0 0104000+12303

38、0+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-=lEIlEIlEIlEIlEAlEIlIA梁解:1)原始数据

39、及编码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. 6

40、31. 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)

41、形成k38郑州大学土木工程学院 樊友景编制3) 形成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

42、) 形成K=00020112=301201=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 .83

43、094. 6031. 294. 6031. 2103k 1 0 3 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

44、 2 3 123K=10340郑州大学土木工程学院 樊友景编制5) 求等效节点荷载P=3-30330F: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. 2

45、25. 103303300001 .2608381000000010000100000001000000010000108 .2794. 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. 3

46、08 .2747. 3047. 358. 0047. 358. 00005 .52005 .5210342郑州大学土木工程学院 樊友景编制=D0009 .970838=+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.

47、 325. 10-TF09. 343. 0009. 243. 00-=TF41. 875. 4009. 225. 10-=8.412.093.094.47M图(kN.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郑州大学土木工程学院 樊友景编制单元

48、的刚度方程(局部坐标）u1u2X1X2e12)(),(212211uulEAXuulEAX-=-=-=21211111uulEAXXX1X2Y2Y1X1X2yxx-=221122110000010100000101vuvulEAYXYX -=cossin00sincos0000cossin00sincosT注意：桁架单元的结点转角不是基本未知量。无须求等效结点荷载。杆端力是由结点位移产生的。TkTkT=D=TkF13-9 桁架及组合结构的整体分析坐标转换矩阵单元的刚度方程(整体坐标）44郑州大学土木工程学院 樊友景编制 3ll 10kN10kN例13-6 求图示桁架内力(EA=常数)。解：1、

49、编码如图 ；124002、形成k6 , 5,000001010000010124 , 3 , 2 , 1,0000010100000101=-=-=lEAklEAk3、形成kk=k=k=k单元=90-=011011Tk=k=lEA-=01100001011011111111lEATkTkT1000-1000-1000100045郑州大学土木工程学院 樊友景编制单元=45-=11112111Tk=lEA22- -=1111000111112211111111lEATkTkT1111-11111111-1111单元=135-=11112111Tk=lEA22-=111100011111221111

50、1111lEATkTkT1111-1111-1111-1111- 3ll 10kN10kN124004、集成总刚KTT4321 ,0021 =TT00000043=TT00430021 =lEAK-=35. 135. 00035. 035. 1010035. 135. 00135. 035. 15、节点荷载T0010-10=P46郑州大学土木工程学院 樊友景编制-=-00101035. 135. 00035. 035. 1010035. 135. 00135. 035. 1DDCCvuvulEA6、解基本方程lEAvuvuDDCC-=58. 536.2142.1494.267、杆端力计算=D.

51、-=-=D=042.14042.140042.1494.2601001000000100100000010100000101TkF-=-=D=058. 5058. 558. 536.2142.1494.2610000100001000010000010100000101TkF0021lEA-=D0042.1494.26 4321.=D lEA-=D58. 536.2142.1494.2647郑州大学土木工程学院 樊友景编制节点位移分量自由节点位移分量(基本未知量，相应的节点荷载已知)受约束的位移分量(已知量，相应的约束反力未知)先处理法1、节点位移分量中不含受约束的支座位移，节点力分量中不 含

52、未知的支座反力。2、由单刚考虑边界条件K3、对于具有非刚性连接、支承节点较多且分散、不考虑轴向 变形的结构最为方便。可减少内存，提高计算速度。但要 对各节位移进行统一编码，形成各单元的定位向量。后处理法1、节点位移分量中含有受约束的支座位移，节点力分量中含 有未知的支座反力。2、由单刚考虑边界条件KK(原始刚度矩阵,奇异)3、每个节点位移分量数相同， 的阶数是节点总数乘节点位 移 分量数，整个分析过程便于编制通用程序。适用于节点 多支座约束少，考虑轴向变形的结构。但占用内存大。K48郑州大学土木工程学院 樊友景编制后处理法边界条件的处理 在后处理法中，由于没有考虑边界条件，由k集成的 是奇异矩

53、阵，由单元集合成的体系是自由体，具有刚体位移。 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该法虽为近似处理，程序设计容易实现，故被广泛应用。该法虽为近似处理，程序设计容易实现，故被广泛应用。RRd049郑州大学