平面钢架程序设计.ppt

第八讲平面钢架程序设计,第一节概述第二节平面刚架内力和位移计算的框图与程序第三节平面刚架程序计算实例,第一节概述,一程序编制说明,平面刚架在荷载作用下的结点位移和杆端力图3-15所示的六种固端梁的等效结点荷载平面和桁架组合结构,本程序建立整体刚度矩阵时，采用直接刚度法中的“先处理法”。各单元为等截面杆。当两结点间的梁自然段抗弯刚度不等时，可将截面突变点作为结点处理。本章根据第一章所述的计算原理和方法，讨论用矩阵位移计算平面刚架结构的程序设计,本程序可用于计算,编制说明,计算模型与计算方法,以杆件联结点、支座结点、截面突变和外伸端点作为计算结点，任意两结点间的杆件作为计算单元。在局部坐标系下，单元两端的杆端力、杆端位移列阵如式（1-15）、式（1-16）所示，即：,1.计算模型,在局部坐标系下，单元刚度矩阵如下式所示，即:,3.支承条件的引入及整体刚度矩阵的组集,整体刚度矩阵的组集采用“直接刚度矩阵法”。整体坐标系下单元刚度矩阵各元素的下标由单元定位数组确定，即在组集整体刚度矩阵之前引入支承条件。,确定单元数组时应注意以下两个问题,（1）支座结点的未知位移分量编号,若单元的某一端与支座相联，则单元支座结点的未知位移分量信息应按表5-1输入,（2）杆件联结点未知位移分量编号,若单元的某一端与其他杆件相联，则应首先根据联结情况确定结点编码，而后再确定与结点相应的单元未知位移分量编码。现将常遇到的几种情况列于表5-2中。,开始,（1）标题及数组说明（读入题目序号NO）,NO是否为零,（2）读入数据并打印,（3）形成整体刚度矩阵,（4）形成结点荷载,（5）解方程并打印杆端位移,（6）计算并打印杆端力,结束,(7)子程序READ,(8)子程序MKE,（9）子程序MR,（10）子程序MAKE,（11）子程序CALM,（12）子程序MK,（13）子程序PE,（14）子程序MULV6,（15）子程序MF,（16）子程序SOLV,（17）子程序MADE,（19）子程序MULV,（18）子程序TRAN,第二节,平面刚架内力和位移计算的框图与程序,一、,程序标识符说明,平面刚架静力分析程序名为,PFSAP(PlaneFarmeStructuralAnalysisProgram),。其,主要标识符说明如下,：,TITLE(2,0),算力标题。实型数组，输入参数。,NJ,结点总数。,，输入参数。,N,结点自由度，即整体刚度矩阵的阶数。,整型变量，输入参数。,NE,单元总数。,整型变量，输入参数。,NM,单元类型总数。,同类型的单元,E,、,A,、,I,相同。,整型变量，输入参数。,NPJ,结点载荷总数。,整型变量，输入参数。,NPF,非结点载荷总数。,整型变量，输入参数。,JN,（,3,，,100,）,结点位移号数组。,整型数组，输入参数。,X,（,100,），,Y,（,100,）,结点坐标数组。,实型数组，输入参数。,JE,（,2,，,100,）,单元两端结点号数组。,整型,数组，输入参数。,JEAI,（,100,）,单元类型信息数组，,JEAI,（,e,）为单元的类型号。,同类型的单元,弹性模量、横截面积及极惯性矩均相同。,整型数组，输入参数。,RT,（,6,，,6,）,单元坐标转换矩阵的转置矩阵。实型数组。,P,（,100,）,综合结点荷载数组。实型数组。,FF,（,6,）,局部坐标系下单元杆端力数组。实型数组。,FE,（,6,）,局部坐标系下单元等效荷载数组。实型数组。,AFE,（,6,）,整体坐标系下单元等效荷载数组。实型数组。,D,（,50,）,整体坐标系下自由解点位移数组。实型数组。,ADE,（,6,）,整体坐标系下单元杆端位移数组。实型数组。,DE,（,6,）,局部坐标系下单元杆端位移数组。实型数组。,F,（,3,）,整体坐标系下结点位移数组。实型数组。,NO,计算题目的符,号，整型变量，输入参数。,SQRT,标准函数，计算非负实数的平方根。,READ,子程序，输入原始数据。,MKE,子程序，计算局部坐标系下单元刚度矩阵。,MR,子程序，计算单元坐标转换矩阵。,平面刚架静力分析原程序（PFSAP.FOR)1:CANALYSISPROGRAMFORPLANEFRAME2:REALK(200,200),KE(6,6),AKE(6,6),X(100),Y(100),AL(100)3:&EAI(3,100),PJ(100),PF(2,100),R(6,6),P(100),FF(6),4:&FE(6),D(100),ADE(6),DE(6),RT(6,6),AFE(6),F(3)5:INTEGERJE(2,100),JN(3,100),JPJ(100),JPF(2,00),M(6),6:&JEAI(100),NO7:OPEN(6,FILE=PFSAP.IN)8:OPEN(8,FILE=PFSAP.OUT)9:READ(6,)NO10:IF(NO.EQ.0)STOP11:WRITE(8,(/A5,I3,A1)(NO.=,NO,12:CALLREAD(NJ,N,NE,NM,NPJ,NPF,JN,X,Y,JE,&JEAI,EAI,JPJ,PJ,JPF,PF)13:DO5I=1,N14:P(I)=015:DO5J=1,N16:K(I,J)=017:DO10IE=1,NE18:CALLMKE(KE,IE,JE,JEAI,EAI,X,Y,AL)19:CALLMR(R,IE,JE,X,Y)20:CALLMAKE(KE,R,AKE),2:6:数组说明7:8:定义输入输出设备9:11:读题目序号。若序号为零，停止运算，程序结束，否则，将序号打印出来。,18:调用MKE，形成局部坐标系下的单元刚度矩阵。19:调用MR，形成单元坐标转换矩阵。20:调用MAKE，形成整体坐标系的单元刚度矩阵。在该子程序中又调用了TRAN和MULV。,21:CALLCALM(M,IE,JN,JE22:CALLMK(K,AKE,M)23:10CONTINUE24:DO20IP=1,NP25:CALLMR(R,JPF(1,IP),JE,X,Y）26:CALLTRAN(R,RT）27:CALLPE(FE,IP,JPF,PF,AL）28:CALLMULV6(RT,FE,AFE)29:CALLCALM(M,JPF(1,IP),JN,JE)30:CALLMF(P,AFE,M）31:20CONTINUE32:DO30I=1,NP33:30P(JPJ(I)=P(JPJ(I)+PJ(I)34:CALLSLOV(K,P,D,N）35:WRITE(8,(/2(26(1H),A)RESULTSOFCALCULATION36:WRITE(8,40)37:40FORMAT(/5X,NO.N,4X,X-DISPLACEMENT,2X,38:Y-DISPLACEMENT,3X,ANG.ROT.(RAD)39:DO60KK=1,NJ40:DO50II=1,341:F(II)=042:I1=JN(II,KK),21:调用CALM,形成单元定位数组。22:调用MK子程序，按单元定位数组给出的编号，“对号入座”组集整体坐标系下的整体结构刚度矩阵K.,25-30:分别调用MR,TRAN（单元坐标转换矩阵的转置矩阵）,PE（等效结点荷载）,MULV6（6阶矩阵与6元素列阵相乘）,CALM（单元定位数组）,MF（整体坐标系下荷载列阵）32:对直接作用在结点上的荷载循环。33：将直接作用在结点上的荷载按位移编号累加入P.,35:38:打印计算结果表名和结点位移表名。39:45:打印各结点的位移值。,43:50IF(I1.GT.0)F(II)=D(I1)44:60WRITE(8,70)KK,F(1),F(2),F(3)45:70FORMAT(I8,2X,3G16.5)46:WRITE(8,80)47:80FORMAT(/NO.E,5X,N(1),8X,Q(1),8X,M(1),48:DO130IE=1,NE50:CALLMADE(IE,JN,JE,D,ADE)51:CALLMKE(KE,IE,JE,JEAI,EAI,X,Y,AL)52:CALLMR(R,IE,JE,X,Y)53:CALLMULV6(R,ADE,DE)54:CALLMULV6(KE,DE,EF)55:DO100IP=1,NPF56:IF(JPF(1,IP).EQ.IE)THEN57:CALLPE(EF,IP,JPF,PF,AL)58:DO90I=1,659:90FE(I)=FF(I)-FE(I)60:ENDIF61:100CONTINUE62:WRITE(8,110)IE,(FF(I),I=1,6)63:110FORMAT(I5,2X,6G12.5)64:130CONTINUE,50:54:分别调用MADE,MKE,MR,MULV6子程序。56:60:若非结点荷载不为零，则将由于非结点荷载引起的固端力取出并存入FF数组。,65:GOTO166:END（7）子程序READ67:SUBROUTINEREAD(NJ,N,NE,NM,NPJ,NPF,JN,X,Y,JE,JEAI,EAI,68:&JPJ,PJ,JPF,PF)69:REALX(100),Y(100),EAI(3,100),PJ(100),PF(2,100)70:INTEGERJE(2,100),JN(3,1000),JPJ(100),JPF(2,100),JEAI(100)71:&TITLE(20)72:READ(6,(20A4)(TITLE(I),I=1,20)73:WRITE(8,(/7X,20A4)TITLE74:READ(6,)NJ,N,NE,NM,NPJ,NPF75:WRITE(8,(/3(5X,A4,1H:I2)NJ=,NJ,76:&N=,N,NE=,NE,NM=,NM,NPJ=,NPJ,NPF=,NPF77:WRITE(8,5)78:5FORMAT(/4X,NO.(1)(2)(3),10X,X,8X,Y)79:READ(6,10)(JN(J,I),J=1,3),X(I),Y(I)，=1,NJ)80:10FORMAT(2,(3I5,2G16.4)81:DO20I=1,NJ82:20WRITE(8,(2X,1H(,I2,1H),3I6,4X,2F10.3)I,JN(1,I),JN(2,I),83:&JN(3,I),X(I),Y(I)84:WRITE(8,30)85:30FORMAT(/10X,ELEMENTNO.NODE-1NODE-2MATERIALS)86:READ(6,40)（JE(1,I),JE(2,I),JEAI(I),I=1,NE）87:40FORMAT(5(3I5),88:DO50I=1,NE89:50WRITE(8,(14X,I2,3(7X,I3)I,JE(1,I),JE(2,I),JEAI(I)90:READ(6,*)(EAI(I,J),I=1,3),J=1,NM)91:WRITE(8,60)(J,(EAI(I,J),I=1,3),J=1.NM)92:60FORMAT(/3X,NO.MAT,6X,ELASTICMODULUS,8X,93:&AREA,5X,MOMENTOFINERTIA/(I6,9X,3G16.6)94:IF(NPJ.EQ.0)GOTO9095:WRITE(8,(20X,16HNODELLOADS)96:WRITE(8,(16XA)NO.DISP.VALUE97:READ(6,70)(JPJ(I),PJ(I),I=1,NPJ)98:70FORMAT(5(I5,G16.4)99:DO80I=1,NP100:80WRITE(8,(14X,I7,F16.3)JPJ(I),PJ(I)101:90CONTINUE102:IF(NPF.EQ.0)GOTO130103:WRITE(8,(/20X,16HNON-NODELLOADS)104:WRITE(8,(11X,A,8X,A,9X,A)NO.E.NO.LOAD.MODEL,A,C105:READ(6,100)(JPF(1,J),JPF(2,I),PF(1,I),PF(2,I),I=1,NPF)106:FORMAT(2(2I5,2G16.4)107:DO110I=1,NPF108:110WRITE(8,120)(JPF(J,I),J=1,2),PF(1,I),PF(2,I)109:120FORMAT(6X,2I8,10X,2F10.3)110:130CONTINUE112END,（8）子程序MKE113:SUBROUTINEMKE(KE,IE,JE,JEAI,EAI,X,Y,AL)114:REALKE(6,6),X(100),Y(100),EAI(3,100),AL(100),L115:INTEGERJE(2,100),JEAI(100)116:II=JE(1,IE)117:JJ=JE(2,IE)118:MT=JEAI(IE)119:L=SQRT(X(JJ)-X(II)*+(Y(JJ)-Y(II)*2)120:AL(IE)=L121:A1=EAI(1,MT)*EAI(2,MT)/L122:A2=EAI(1,MT)*EAI(3,MT)/L*3123:A3=EAI(1,MT)*EAI(3,MT)/L*2124:A4=EAI(1,MT)*EAI(3,MT)/L125:KE(1,1)=A1126:KE(1,4)=-A1127:KE(2,2)=12*A2128:KE(2,3)=6*A3129:KE(2,5)-12*A2130:KE(2,6)=6*A3131:KE(3,3)=4*A4132:KE(3,5)=-6*A3133:KE(3,6)=2*A4134:KE(4,4)=A1,135:KE(5,5)=12*A2136:KE(5,6)=-6*A3137:KE(6,6)=4*A4138:DO10I=1,6139:DO10K=I,6140:10KE(K,I)=KE(I,K)141:RETURN142:END（9）子程序MR143:SUBROUTINEMR(R,IE,JE,X,Y)144:REALR(6,6),X(100),Y(100),L,CX,CY145:INTEGERJE(2,100)146:I=JE(1,IE)147:J=JE(2,IE)148:L=SQRT(X(J)-X(I)*2+(Y(J)-Y(I)*2)149:CX=(X(j)-X(I)/L150:CY=(Y(J)-Y(I)/L151:DO10J=1,6152:DO10I=1,6153:10R(I,J)=CX154:DO20I=1,4,3155:R(I,I)=CX156:R(I,I+1)=CY157:R(I+1,I)=-CY,158:R(I+1,I+1)=CX159:20R(I+2,I+2)=1160:RETURN161:END(10)子程序MAKE162:SUBRUTINEMAKE(KE,R,AKE)163:REALKE(6,6),R(6,6),RT(6,6),TMP(6,6),AKE(6,6)164:CALLTRAN(R,RT)165:CALLMULV(RT,KE,TMP)166:CALLMULV(TMP,R,AKE)167:RETURN168:END(11)子程序CALM169:SUBROUTINECALM(M,JE,JN,JE)170:INTEGERM(6),JN(3,100),JE(2,100),IE171:DO10I=1,3172:M(I)=JN(I,JE(1,IE)173:10M(I+3)=JN(I,JE(2,IE)174:RETURN175:END(12)子程序MK176:SUBROUTINEMK(K,AKE,M)177:REALK(200,200),AKE(6,6),MAKE子程序需用TRAN和MULV调用,178:INTEGERM(6)179:DO10I=1,6180:DO10J=1,6181:IF(M(I).NE.0.AND.M(J).NE.0)182:&K(M(I),M(J)=K(M(I),M(J)+AKE(I,J)183:10CONTINUE184:RETURN185:END（13）子程序PE186:SUBROUTINEPE(FE,IP,JPF,PF,AL)187:REALFE(6),PF(2,100),AL(100),L188:INTEGERJPF(2,100)189:A=PF(1,IP)190:C=PF(2,IP)191:L=AL(JPF(1,IP)192:IND=JPF(2,IP)193:DO5I=1,6194:5FE(I)=0195:GOTO(10,20,30,40,50,60),IND196:10FE(2)=(7*A/20+3*C/20)*L197:FE(3)=(A/20+C/30)*L*2198:FE(5)=(3*A/20+7*C/20)*L199:FE(6)=-(A/30+C/20)*L*2,200:RETURN201:20FE(5)=A*C*3*(2*L-C)/2/L*3202:FE(2)=A*C-FE(5)203:FE(3)=A*C*2*(6*L*L-8*C*L+3*C*C)/12/L/L204:FE(6)=-A*C*3*(4*L-3*C)/12/L/L205:RETURN206:FE(2)=A*(L-C)*2*(L+C*2)/L*3207:FE(3)=A*C*(C-L)*2/L*2208:FE(5)=A-FE(2)209:FE(6)=-A*C*2*(L-C)/L*2210:RETURN211:FE(2)=-6*A*C*(L-C)/L*2212:FE(3)=A*(L-C)*(L-3*C)/L*2213:FE(5)=-FE(2)214:FE(6)=A*C*(3*C-2*L)/L*2215:RETURN216:FE(1)=A*(1-C/L)217:FE(4)=A*C/L218:RETURN219:END220:FE(4)=FE(1)221:RETURN222:END,梁长为L的两端固定梁的等效结点荷载。见书57页。,(14)子程序MULV6223:SUBROUTINEMULV6(A,B,C)224:REALC(6),A(6,6),B(6)225:DO10I=1,6226:C(I)=0.0227:DO10J=1,6228:10C(I)=C(I)+A(I,J)*B(J)229:RETURN230:END(15)子程序MF231:SUBROUTINEMF(P,AFE,M)232:REALP(100),AFE(6)233:INTEGERM(6)234:DO10I=1,6235:IF(M(I).NE.0)P(M(I)=AFE(I)+P(M(i)236:10CONTINUE237:RETURN238:END(16)子程序SLOV239:SUBROUTINESLOV(AK,P,D,N)240:REALAK(200,200),P(100),D(100)241:DO5I=1,100242:5DO(I)=P(I),243:DO10K=1,N-1244:DO10I=K+1,N245:C=-AK(K,I)/AK(K,K)246:DO20J=