第九章 矩阵位移法

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

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

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

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

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

6、D= kF82v=10312lEI312lEI-EI26l-EI26l-01v=10312lEI312lEI-EI26lEI26l0-=lEIlEIlEIlEIlEIlEIlEIlEIlEAlEAlEIlEIlEIlEIlEIlEIlEIlEIlEAlEAk460260612061200000260460612061200000222323222323e 第二列元素第二列元素变符号即第五列变符号即第五列 第一列元素第一列元素变符号即第四列变符号即第四列 第二行元素第二行元素变符号即第五行变符号即第五行 第一行元素第一行元素变符号即第四行变符号即第四行9特殊单元特殊单元单元的某个或某些杆端位移的

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

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

9、FTF=FTFT=D=DTD=DTT12整体坐标系中的整体坐标系中的单元刚度矩阵单元刚度矩阵D= kF设：将（a）、（b）代入D=TkFTD=TkTFTTTTD= kFFTFT=D=DTT（a）（b）与 比较 D= kFD=TkTFTTTTk，k 同阶,性质类似：131）ijk 表示在整体坐标系第j个杆端位移分量=1时引起的第i个杆端力分量。 2）k 是对称矩阵。3）一般单元的k是奇异矩阵。 I k，k 同阶,性质类似：=0000221122211211221122211211TTkkkkTTkkkkTT14 例13-1 求图示刚架中各单元在整体标系中的单元刚度矩阵。设各杆的几何尺寸相同。l=

10、5m，A=0.5m2, I=1/24m4E=3107kN/m2441025,10300=lEIlEA21 I 15解：（1）求kek1k2=503003012000300-1003003012000300503003012000300-1003003012000300-104（2）求 ke kkIT=,11单元 =9021单元 =0 -=100000001000010000000100000001000010Tk2=500300300030012-1000300300030012-500300300030012-100030030003012104-=100001010100300301200

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

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

13、） 2 =112B （1） 2B （2） 3 =223 将各单元的单刚的行列局部码（i）、（j）换成对应的结点位移总码i、 j，按此行列总码将单刚元素送入总刚。即:k(i)(j)jiK2112213ABC(1)(2) (1)(2)19（2）单刚）单刚k eKe和单元贡献和单元贡献中元素的对应关系中元素的对应关系单元单元单元单元k =4i12i14i12i11(1)(2)(1)(2) 1=21K =11230000000004i12i12i14i1123k =4i22i24i22i22(1)(2)(1)(2) 2=32K =20000000004i22i24i22i2123123单元贡献矩阵是单

14、元刚单元贡献矩阵是单元刚度矩阵，利用度矩阵，利用“单元定单元定位向量位向量”进行进行“换码重换码重排位排位”。2112213ABC(1)(2) (1)(2)20（3 3） 单元集成法的实施单元集成法的实施（定位（定位 累加）累加）K123123000000000k 1 10000000004i12i12i14i1123123k 2 24i12i14i12i1000002i22i24i24i1+4i2123123（1）将）将K置零，得置零，得K=0；（2）将）将k的元素在的元素在K中按中按 定位并进行累加，得定位并进行累加，得K=K；（3）将）将k的元素在的元素在K中按中按 定位并进行累加，得定

15、位并进行累加，得K=K+K；按此作法对所有单元循环一遍，最后即得整体刚度矩阵按此作法对所有单元循环一遍，最后即得整体刚度矩阵K。21例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) (2) 3 0+ 4i31 2 31 2300在

16、给节点位移编码时已经考虑了支承条件。(先处理法)22四、整体刚度矩阵四、整体刚度矩阵 K 的性质的性质（1）整体刚度系数的意义）整体刚度系数的意义: KijD Dj=1 (其余其余D D=0)时产生的结点力时产生的结点力Fi（2）K是对称矩阵是对称矩阵（3）对几何不变体系，）对几何不变体系，K是可逆矩阵，如连续梁是可逆矩阵，如连续梁i1i2D1D2D3F1F2F3F=KD DD=K-1F（4）K是稀疏矩阵和带状矩阵，如连续梁是稀疏矩阵和带状矩阵，如连续梁D1D2D3F1F2F3123nDnFnDn+1Fn+1000000000000000000000000000000000000=+1321n

17、FFFFDDDD+1321n4i12i12i12i22i24i2+4i34i1+4i24in2i32in231n2312n+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+ = =+=-+njiKKnjiiKiKiKjjjjjjjjjnnn1n2324情况复杂：1）结点位移分量增加到三个；2）各杆方向不尽相同，要进行坐标变换；3

18、）除了刚结点，还要考虑铰结点等其它情况。9-5 刚架的整体刚架矩阵思路要点思路要点：（：（1）设各单元已形成了整体座标系下的单元刚度矩阵）设各单元已形成了整体座标系下的单元刚度矩阵；ek（2）各）各 经由经由e 进行累加集成进行累加集成K。与连续梁相比：与连续梁相比： (1)各单元考虑轴向变形各单元考虑轴向变形；(2)每个刚结点有三个位移每个刚结点有三个位移； (3)要采用整体座标要采用整体座标；(4)要处理非刚结点的特殊情况要处理非刚结点的特殊情况。25情况复杂：1）结点位移分量增加到三个；2）各杆方向不尽相同，要进行坐标变换；3）除了刚结点，还要考虑铰结点等其它情况。1、结点位移分量的统一

19、编码总码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 0TACB9-5 刚架的整体刚架矩阵263、单元集成过程503003012000300-1003003012000300503003012000300-1003003012000300-104k =1 1 2 3 0 0 4K=1 2 3 4300 00001230100030100500305030104k2=500

20、300300030012-1000300300030012-500300300030012-1000300300030121041 2 3 0 0 0 +12+030+0+300+030+0+100求单刚27四、铰结点的处理四、铰结点的处理K 求单元常数T单元刚度矩阵单元刚度矩阵程序设计框图（局部：集成整体刚度矩阵）程序设计框图（局部：集成整体刚度矩阵） 1122刚结点刚结点：变形连续，：变形连续，截面截面1 1和截面和截面2 2具有相同具有相同的结点位移。的结点位移。铰结点铰结点：部分变形连续，：部分变形连续，截面截面1 1和截面和截面2 2具有相同的结点具有相同的结点线位移；而其角位移不相

21、等。线位移；而其角位移不相等。281）结点位移分量的统一编码总码 铰结点处的两杆端结点应看 作半独立的两个结点（C1和 C2） 它们的线位移相同， 角位移不同，00012321A C1B D000456475C234、铰结点的处理线位移采用同码，角位移采用异码。2）单元定位向量： =11 2 3 4 5 6T =21 2 3 0 0 0T =34 5 7 0 0 0T3）按次序进行单元集成：29503003012000300-1003003012000300503003012000300-1003003012000300-104k =1 1 2 3 4 5 61 2 3 4 5 6 73005

22、00-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=500300300030012-1000300300030012-500300300030012-1000300300030121044 5 7 0 0 0 +12+030+0+300+030+0+10030?30F= KD D (1)结

23、构体系刚度方程：结构体系刚度方程：一、位移法基本方程一、位移法基本方程k11D 1+ k12D 2+ + k1nD n+F1P=0 k21D 1+ k22D 2 + + k2nD n+F2P=0 kn1D 1+ kn2D 2+ + knnD n+FnP=0 KD +FP =0 .(2)F +FP =0 .(3)将将(1)式代入式代入(2)式：式： 表示结点位移表示结点位移D D和结点力和结点力F之间的关系，反映了结构的刚度性质，而之间的关系，反映了结构的刚度性质，而不涉及原结构上作用的不涉及原结构上作用的实际荷载实际荷载，并不是原结构的位移法基本方程。，并不是原结构的位移法基本方程。基本体系在

24、荷载单独作用下基本体系在荷载单独作用下产生的结点约束力。产生的结点约束力。基本体系在结点位移单独作基本体系在结点位移单独作用下产生的结点约束力。用下产生的结点约束力。9-6 9-6 等效结点荷载等效结点荷载311、整体刚度方程、整体刚度方程 F=K （a） 表示由表示由F结点力的关系式。反映了结点的刚度性质，结点力的关系式。反映了结点的刚度性质， 不涉及结构上的实际荷载。不涉及结构上的实际荷载。2、位移法基本方程、位移法基本方程 在给结点位移分量编总码时，已考虑了结构的支承连接情况，在给结点位移分量编总码时，已考虑了结构的支承连接情况， K是非奇异矩阵。是非奇异矩阵。 如果已知结构上的结点荷载

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

26、载不是作用在结点上的集中荷载，当荷载不是结点集中荷载时，应化成等效结点荷载。是结点集中荷载时，应化成等效结点荷载。9-6 9-6 等效结点荷载等效结点荷载32（1）设荷载单独作用（结点位移为零）设荷载单独作用（结点位移为零）此时在基本结构此时在基本结构中引起的结点约束力，记为中引起的结点约束力，记为FP。（2）设结点位移单独作用（荷载设为零）设结点位移单独作用（荷载设为零）此时在基本结此时在基本结构中引起的结点约束力为构中引起的结点约束力为F=K。 位移法基本方程应为位移法基本方程应为 F+FP=0即即 K+FP=0 （9-47）9-6 9-6 等效结点荷载等效结点荷载33123FP2FP1F

27、P3P2P1P3结点约束力FP =-FP2 =-FP1-FP3=1323=2=1等效结点荷载P=FP可由位移法基本方程可由位移法基本方程(b)求得求得.注意：注意：（1）非结点荷载与等效结点荷载等效的条件是，两者产生相同）非结点荷载与等效结点荷载等效的条件是，两者产生相同 结点位移。结点位移。（2）除了结点位移外，等效结点荷载与原荷载产生的其它位移）除了结点位移外，等效结点荷载与原荷载产生的其它位移 和内力并不相同。和内力并不相同。（3）等效结点荷载为位移法基本体系附加约束中约束力的负值）等效结点荷载为位移法基本体系附加约束中约束力的负值 而约束力为各固端力之和。所以求结构等效结点荷载应该而约

28、束力为各固端力之和。所以求结构等效结点荷载应该 先求出单元的等效结点荷载，它是先求出单元的等效结点荷载，它是单元固端力单元固端力的负值。的负值。0=+PFP位移法方程：位移法方程：D= KP344、按单元集成法求整体结构的、按单元集成法求整体结构的等效结点荷载等效结点荷载局部坐标系中在单元两端加六个附加约束，在给定荷载作用下，求出单元固端约束力向量PFe整体坐标系中的单元等效结点荷载 PTFTP-=ee整体结构的等效结点荷载P 由各单元P 中的元素按在P中进行定位并累加。e等效结点荷载与直接结点荷载叠加，即得结构的结点荷载。=eP单元等效结点荷载（局部坐标系）PFe354.8kN/m2.5m2

29、.5m5m8kN12yx例9-3 求图示结构的等效结点荷载P.PFe解:1)求单元,.10120111mkNMkNYXPPP-=-=mkNMkNYXPPP.10120222=-=单元,.540111mkNMkNYXPPP=mkNMkNYXPPP.540222-= TPF1012010120-= TPF540540-=362)求 PTFTP-=单元的倾角1=01 2 3 0 0 4单元的倾角2=90123000P=123401210-10+4 +05 -=-=-=504504540540100000001000010000000100000001000010PTFTP TPPTFIFTP1012

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

31、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)原始数据及编码39.1031. 212,1094. 66,108 .274,109 .132,103 .83,1094.

32、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 .83-.1058. 012,1047.

33、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)形成k403) 形成k 单元、(=90o)坐标转换矩阵为 -=100000001000010000000100000001000010T-=1000010108 .2794. 60

34、94. 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=22211211kkkkK41-=6 .5547. 394. 647. 388.83094. 6081.549 .1347. 3047.

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

36、394. 647. 388.83094. 6081.54106543213-=DDDDDD5 .9613. 58244 .2813. 584765432143-=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. 6094. 631. 20003 .83

37、9 .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 .5210344=D0005 .9613. 5824=+D= PFTkF-=-38. 424. 143. 004.

38、 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. 143. 0-=8.492.093.044.38

39、M图(kN.m)4.761.240.43 1.241.24Q图(kN)N=0.43N=1.24N=0.43N图(kN)451）结点位移分量的统一编码总码 在刚结点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 忽略轴向变形时矩形刚架的整体分析46503003012000300-1003003012000300503003012

40、000300-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-1000300300030121041 0 4 0 0 0104

41、000+123030+100K=104 47例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)原始数据及编码0001

42、02000103213xy48.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.

43、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)形成k493) 形成k 单元、(=90o)坐标

44、转换矩阵为 -=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=301201=000301350 1 0 2 0

45、 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 0 0 0103000+2.316.946

46、.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=103515) 求等效节点荷载P=3-30330F:P固端力1=-=3033-03PPTF

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

48、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 .5210353=D0009 .970838=+D= PFTkF-=-47. 425

49、. 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图(kN.m)4.751.250.43 1.251

50、.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-= 轴向变形影响不大。54单元的刚度方程(局部坐标）u1u2X1X2e12)(),(212211uulEAXuulEAX-=-=-=21211111uulEAXXX1X2Y2Y1X1X2yxx-=221122110000010100000101vuvulE

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

52、11lEATkTkT1000-1000-1000100056单元=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. 135. 00135. 035. 15、节

53、点荷载T0010-10=P57-=-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-=D0042.1494.26 4321.=

54、D lEA-=D58. 536.2142.1494.2658节点位移分量自由节点位移分量(基本未知量，相应的节点荷载已知)受约束的位移分量(已知量，相应的约束反力未知)先处理法1、节点位移分量中不含受约束的支座位移，节点力分量中不 含未知的支座反力。2、由单刚考虑边界条件K3、对于具有非刚性连接、支承节点较多且分散、不考虑轴向 变形的结构最为方便。可减少内存，提高计算速度。但要 对各节位移进行统一编码，形成各单元的定位向量。后处理法1、节点位移分量中含有受约束的支座位移，节点力分量中含 有未知的支座反力。2、由单刚考虑边界条件KK(原始刚度矩阵,奇异)3、每个节点位移分量数相同， 的阶数是节点总数乘节点位 移 分量数，整个分析过程便于编制通用程序。适用于节点 多支座约束少，考虑轴向变形的结构。但占用内存大。K59后处理法边界条件的处理 在后处理法中，由于没有考虑边界条件，由k集成的 是奇异矩阵，由单元集合成的体系是自由体，具有刚体位移。 KPK=D没有确定的位移解。位移边界条件处理的三种方法：1、划行划列法 编制程序较复杂，不常采用。2、主对角元置大数法设第 i 个节点位移分量 (已知)0di=D =D D DD nininnninn