哈工大电算实习报告

1、1、 任务已完成如下了任务：1 增加荷载类型（原来程序中只有集中力和均布力）。2 计算指定截面的弯矩和剪力。3 计算指定截面的弯矩影响线。4 计算指定截面的剪力影响线。5 计算指定支座的反力影响线。6 计算均布活荷作用下的弯矩包络图。7 计算均布活荷作用下的剪力包络图。二 添加的子程序框图1计算指定截面的弯矩流程图开始求出只有杆端弯矩时截面的弯矩（左右侧弯矩）I=非结点荷载的个数（NF）荷载是否作用在此单元I=I+1调用SCWX子程序输出左右侧弯矩结束根据不同荷载类型（ID），按简支梁求出截面的左右侧弯矩，再叠加上原有弯矩I=I+1NNYY调出杆端弯矩，所求截面所在单元及其在单元中的位置2.计

2、算指定截面的弯矩流程图开始求出只有杆端弯矩时截面的剪力（左右侧剪力）I=非结点荷载的个数（NF）荷载是否作用在此单元I=I+1调用SCWX子程序输出左右剪力结束根据不同荷载类型（ID），按简支梁求出截面的左右侧剪力，再叠加上原有剪力I=I+1NNYY调出杆端弯矩，所求截面所在单元及其在单元中的位置3.支座反力计算的流程图开始调入所求支座对应的节点号计算支座左侧引起的支座反力FX1计算支座右侧引起的支座反力FX2FX=FX1+FX2输出支座反力结束3.1计算左侧荷载引起的支座反力的流程图 开始I=非结点荷载的个数（NF）荷载是否作用在此单元I=I+1结束根据不同荷载类型（ID），按简支梁求出引起

3、的支座反力，再叠加上原有的支座反力I=I+1NNYY所求支座反力的序号是否为1调出所求支座的序号NYFX1=03.2计算右侧荷载引起的支座反力的流程图开始I=非结点荷载的个数（NF）荷载是否作用在此单元I=I+1结束根据不同荷载类型（ID），按简支梁求出引起的支座反力，再叠加上原有的支座反力I=I+1NNYY所求支座反力的序号是否为NJ调出所求支座的序号NYFX2=04.弯矩影响线的流程图开始调用截面所在单元及其在单元中的位置单位荷载作用的单元NN=1NN=NEI=20单元是否作用在支座处I=I+1Y调用求截面弯矩的程序I=I+1NN=NN+1NYNY输出MX=0X=I*L/20结束N5.剪力

4、影响线流程图开始调用截面所在单元及其在单元中的位置单位荷载作用的单元NN=1NN=NEI=20单元是否作用在支座处I=I+1Y调用求截面左右侧剪力的程序I=I+1NN=NN+1NYNY输出FXZ=FXY=0X=I*L/20结束N6.支座反力的影响线开始调用截面所在单元及其在单元中的位置单位荷载作用的单元NN=1NN=NEI=20单元是否作用在支座处I=I+1Y调用求支座反力的程序I=I+1NN=NN+1NYNY输出FX=0X=I*L/20结束N是否作用在所求支座处NYI=I+1输出FX=17.弯矩包络图的流程图（注：对于弯矩的求解运用叠加原理，对每一跨荷载引起连续梁的内力进行叠加）开始计算隔跨

5、布置均布活荷载所引起的截面内力，用矩阵WX(NE,21)来存储。计算另一种隔跨布置均布活荷载所引起的截面内力，用矩阵WXX(NE,21)来存储。比较WX(NE,21)与WXX(NE,21)，如果WXX对应的元素比WX大，则将WXX对应的元素赋值给WX对于引起支座最大负弯矩的均布活荷载引起的各截面弯矩用WXX(NE,21)来存储，再与WX比较，如果WXX对应的元素比WX大，则将WXX对应的元素赋值给WX（WX即为各截面的最大弯矩）。按行输出WX(NE,21)最大弯矩计算隔跨布置均布活荷载所引起的截面内力，用矩阵WX(NE,21)来存储。计算另一种隔跨布置均布活荷载所引起的截面内力，用矩阵WXX(

6、NE,21)来存储。比较WX(NE,21)与WXX(NE,21)，如果WXX对应的元素比WX小，则将WXX对应的元素赋值给WX对于引起支座最大负弯矩的均布活荷载引起的各截面弯矩用WXX(NE,21)来存储，再与WX比较，如果WXX对应的元素比WX小，则将WXX对应的元素赋值给WX（WX即为各截面的最小弯矩）。按行输出WX(NE,21)最小弯矩结束8.剪力包络图的流程图(此时所加的荷载与求弯矩包络图时的荷载相同,对于剪力的求解运用叠加原理，对每一跨荷载引起连续梁的内力进行叠加)：开始计算隔跨布置均布活荷载所引起的截面内力，用矩阵FX(NE,21)来存储。计算另一种隔跨布置均布活荷载所引起的截面内

7、力，用矩阵FXX(NE,21)来存储。比较FX(NE,21)与FXX(NE,21)，如果FXX对应的元素比FX大，则将FXX对应的元素赋值给FX对于引起支座截面最大剪力的均布活荷载引起的各截面剪力用FXX(NE,21)来存储，再与FX比较，如果FXX对应的元素比FX大，则将FXX对应的元素赋值给WX（WX即为各截面的最大弯矩）。按行输出FX(NE,21)最大剪力计算隔跨布置均布活荷载所引起的截面内力，用矩阵FX(NE,21)来存储。计算另一种隔跨布置均布活荷载所引起的截面内力，用矩阵FXX(NE,21)来存储。比较FX(NE,21)与FXX(NE,21)，如果FXX对应的元素比FX小，则将FX

8、X对应的元素赋值给FX对于引起支座最大剪力的均布活荷载引起的各截面剪力用FXX(NE,21)来存储，再与FX比较，如果FXX对应的元素比FX小，则将FXX对应的元素赋值给FX（FX即为各截面的最小剪力）。按行输出FX(NE,21)最小剪力结束三 程序功能本程序可以解决任意跨数连续梁在集中力偶，集中力以及任意分布的均布荷载作用下连续梁的内力（包括剪力和弯矩）以及支座反力的求解问题；可以解决不同跨的截面惯性矩不同的问题，但要求连续梁的弹性模量相同以及任意一跨内惯性矩保持不变；也可以计算在单位荷载作用下指定截面的弯矩，剪力以及引起某支座的反力的数据并结合Excel求出对应的弯矩，剪力和支座反力影响线

9、；还可以用于求解在均布荷载作用下各截面最大，最小的剪力与弯矩的数据并结合Excel求出弯矩和剪力包络图。四 使用说明4.1建立文本文档用于存储初始数据，文本文档应按顺序输入连续梁以及荷载的具体情况，应分别按以下顺序输入：支承类型，节点荷载个数，非节点荷载个数，弹性模量；杆长,惯性矩；节点荷载大小,对应未知数序号（如没有节点荷载可不用输入）；非结点荷载值,起点位置,终点位置,单元号,荷载类型号（如没有非节点荷载可不用输入）；所求弯矩对应截面所在单元以及在单元中的位置；所求剪力对应截面所在单元以及在单元中的位置；所求支座反力的节点号；弯矩影响线对应截面所在单元以及对应截面在单元中的位置；剪力影响线

10、对应截面所在单元以及对应截面在单元中的位置；支座反力影响线对应的节点号。具体数据如下：支承类型1：两端铰接支承类型2：左端固结，右端铰接支承类型3：左端铰接，右端固结支承类型4：两端固结对于集中荷载（如集中力，集中力偶），起点位置,终点位置输入相同的数据（即作用点的位置），对于均匀分布荷载对应起始点以及终止点的位置。4.2力的正负号规定：外荷载：竖向荷载以向上为正，向下为负；弯矩以逆时针为正，顺时针为负。杆端弯矩以顺时针为正，以逆时针为负。各截面的内力：弯矩以下侧受拉为正，上侧受拉为负；剪力使隔离体产生顺时针旋转为正，产生逆时针旋转为负。4.3由于数组定义数组大小的限制，本程序最多可以解决20

11、跨连续梁的结构，对于一般情况都能满足要求。4.4单位的规定：力的单位：KN 力偶单位：KN*m 长度单位：m4.5荷载类型的规定：1.集中力 2.集中力偶 3.任意均布荷载五源程序如下 连续梁静力计算程序 PROGRAM LXL DIMENSION GC(20),GX(20),PJ(20,2),PF(10,5),DK(2,2),P(45) DIMENSION F0(2),WY(2),F(2),ZK1(45),ZK2(45),GDM(20,2) character*12 dat0,dat1 write(*,*)输入:初始数据文件名 read(*,10)dat0 write(*,*)输入:计算结果

12、数据文件名 read(*,10)dat110 format(A12) open(8,file=dat0,status=old) open(9,file=dat1,status=new) write(9,11)11 format(8X,*连续梁内力计算*) READ(8,*)NE,NZ,NP,NF,E0 WRITE(9,12)NE,NZ,NP,NF,E0 12 FORMAT(1X,单元数=,6x,I3,5x,支承类型=,I3,5x,节点荷载个数=I3,/, 1 1x,非节点荷载个数=,5x,I3,弹性模量=,F12.4) NJ=NE+1 CALL SRSJ(NE,NP,NF,GC,GX,PJ,P

13、F,N,X,NN,XX,K,M,XL,MM,XLL,KK) CALL XCP(NJ,NP,NF,NE,P,PJ,PF,F0,GC) CALL JCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX) CALL ZCCL(NZ,NJ,ZK1,ZK2,P) CALL FCQJ(NJ,ZK1,ZK2,P)CALL SCWY(NJ,P) CALL GDL(NE,NJ,NF,E0,P,WY,F,F0,DK,PF,GC,GX,GDM)CALL SCGDL(GDM,NE)CALL JMM(NF,GC,GDM,N,X,PF,NE,WXZ,WXY)CALL SCWX(WXZ,WXY)CALL JMF(NF

14、,GC,GDM,NN,XX,PF,NE,FXZ,FXY)CALL SCFX(FXZ,FXY)CALL ZZFL(NF,GC,GDM,K,PF,NE) CALL MIL(GC,NE,NJ,E0,GX,M,XL,NZ)CALL FIL(GC,NE,NJ,E0,GX,MM,XLL,NZ)CALL ZZFLIL(GC,NE,NJ,E0,GX,KK,NZ)CALL WBLT(GC,NE,NJ,E0,GX,NZ)CALL FBLT(GC,NE,NJ,E0,GX,NZ) ENDC *c 输入初始数据SRSJ子程序C * SUBROUTINE SRSJ(NE,NP,NF,GC,GX,PJ,PF,N,X,NN,

15、XX,K,M,XL,MM,XLL,KK) DIMENSION GC(NE),GX(NE),PJ(NP,2),PF(NF,5)C GC(NE)：杆长 GX(NE)：惯性矩C 结点荷载:PJ(NP,1)=荷载大小 PJ(NP,2)=对应未知数序号 C 非结点荷载 : PF(NE,1)=荷载大小 PF(NE,2)=荷载起始点位置 PF(NE,3)=荷载终点位置 PF(NE,4)=单元号 PF(NE,5)=荷载类型号C 截面的弯矩，剪力：N=单元号 X=所求截面在单元中的位置 WRITE(9,11)11 FORMAT(1X,杆长,惯性矩 GC(NE),GX(NE) READ(8,*)(GC(I),GX

16、(I),I=1,NE)write(9,12)(GC(I),GX(I),I=1,NE)12FORMAT(10X,2F8.3) IF(NP.GT.0)THEN WRITE(9,22)22 FORMAT(1X,节点荷载大小,对应未知数序号 PJ(I,1),PJ(I,2) READ(8,*)(PJ(I,1),PJ(I,2),I=1,NP)write(9,23)(PJ(I,1),PJ(I,2),I=1,NP)23FORMAT(1X,2F8.3) END IF IF(NF.GT.0) THEN WRITE(9,33) 33 FORMAT(1X,非结点荷载值,起点位置,终点位置,单元号,荷载类型号) REA

17、D(8,*)(PF(I,J),J=1,5),I=1,NF) WRITE(9,34)(PF(I,J),J=1,5),I=1,NF) 34 FORMAT(1X,5F8.3) END IFwrite(9,35)35format(1x,所求弯矩对应截面所在单元,在单元中的位置)read(8,*) N,Xwrite(9,36) N,X36 FORMAT(1X,I2,15x,F5.2) write(9,45)45format(1x,所求剪力对应截面所在单元,在单元中的位置)read(8,*)NN,XXwrite(9,46)NN,XX46 FORMAT(1X,I2,15x,F5.2)write(9,37)3

18、7format(1X,所求支座反力的节点号)read(8,*) Kwrite(9,38) K 38 FORMAT(1X,I2)write(9,39)39 format(1X,弯矩影响线对应截面所在单元,对应截面在单元中的位置)read(8,*)M,XLwrite(9,40)M,XL40 FORMAT(1X,I2,15X,F5.2)write(9,43)43 format(1X,剪力影响线对应截面所在单元,对应截面在单元中的位置)read(8,*)MM,XLLwrite(9,44)MM,XLL44 FORMAT(1X,I2,15X,F5.2)write(9,41)41 format(1X,支座反

19、力影响线对应的节点号)read(8,*) KKwrite(9,42) KK42 FORMAT(1X,I2) END C *c 计算第NHF个非结点荷截引起的等效结点荷截列阵F0C * SUBROUTINE DJH(NHF,NE,NF,PF,F0,GC) DIMENSION PF(NF,5),F0(2),GC(NE)C G荷载值,Q-荷载起始点位置,Z-荷载终点位置，NT-单元号,ID荷载类型号 G=PF(NHF,1)Q=PF(NHF,2) Z=PF(NHF,3) NT=INT(PF(NHF,4)+0.1) ID=INT(PF(NHF,5)+0.1) BL=GC(NT) GOTO(10,20,3

20、0),IDC 集中力作用下的等效荷截列阵P0(2)10 D=BL-QC1=Q/BLC2=C1*C1 F0(1)=G*Q*D*D/BL/BL F0(2)=-G*D*C2GOTO 200C 集中力偶作用下的等效荷截列阵P0(2)20D=BL-QC1=Q/BLC2=C1*C1 F0(1)=G*D*(-2+3*D/BL)/BLF0(2)=G*Q*(-2+3*Q/BL)/BLGOTO 200C 均布荷载作用下的等效荷载列阵F0(2)30D=BL-QC1=Q/BL C2=C1*C1C3=Z/BLC4=C3*C3 F0(1)=(G*Z*Z*(6.0-8.0*C3+3.0*C4)-G*Q*Q*(6.0-8.0

21、*C1+3.0*C2)/12.0 F0(2)=(-G*Z*Z*Z*(4.0-3.0*C3)+G*Q*Q*Q*(4.0-3.0*C1)/12.0/BL200 RETURN ENDC *c 计算第NE0个单元的单元刚度矩阵DK(2,2)C * SUBROUTINE DG(NE0,NE,E0,GC,GX,DK) DIMENSION GC(NE),GX(NE),DK(2,2) DO 15 I=1,2 DO 10 J=1,2 DK(I,J)=0.010 CONTINUE15 CONTINUEC BL杆长,DI-惯性矩,S-线刚度 BL=GC(NE0) DI=GX(NE0) S=E0*DI/BL DK(1

22、,1)=4.0*S DK(1,2)=2.0*S DK(2,1)=2.0*S DK(2,2)=4.0*S ENDc *C 集成总体刚度矩阵,ZKl,ZK2分别存主对角元素和付对角元素c * SUBROUTINE JCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX) DIMENSION ZK1(NJ),ZK2(NJ),DK(2,2),GC(NE),GX(NE) DO 10 I=1,NJ ZK1(I)=0.0 ZK2(I)=0.010 CONTINUE DO 100 I=1,NE CALL DG(I,NE,E0,GC,GX,DK) ZK1(I)=ZK1(I)+DK(1,1) ZK2(I)=

23、ZK2(I)+DK(1,2) ZK1(I+1)=ZK1(I+1)+DK(2,2)100 CONTINUE ENDC *C 形成总荷裁矩阵C * SUBROUTINE XCP(NJ,NP,NF,NE,P,PJ,PF,F0,GC) DIMENSION P(NJ),PJ(NP,2),PF(NF,5),F0(2),GC(NE) DO 10 I=1,NJ P(I)=0.010 CONTINUE IF(NP.GT.0) THEN DO 20 I=1,NP J=INT(PJ(I,2)+0.1) P(J)=PJ(I,1)20 CONTINUE END IF IF(NF.GT.0) THEN DO 100 NF

24、0=1,NF CALL DJH(NF0,NE,NF,PF,F0,GC) ND=INT(PF(NF0,4)+0.1) P(ND)=P(ND)+F0(1) P(ND+1)=P(ND+1)+F0(2)100 CONTINUE END IF ENDc *c 进行支承条件处理c * SUBROUTINE ZCCL(NZ,NJ,ZK1,ZK2,P) DIMENSION ZK1(NJ),ZK2(NJ),P(NJ) GOTO(10,20,30,20),NZ10 GOTO 10020 ZK1(1)=1.0 P(1)=0.0 ZK2(1)=0.0 IF(NZ.EQ.4) GOTO 30 GOTO 10030 ZK

25、1(NJ)=1.0 ZK2(NJ-1)=0.0 P(NJ)=0.0100 CONTINUE ENDc *c 解方程求节点位移Pc * SUBROUTINE FCQJ(NJ,ZK1,ZK2,P) DIMENSION ZK1(NJ),ZK2(NJ),P(NJ) DO 10 I=1,NJ-1 ZK1(I+1)=ZK1(I+1)-ZK2(I)*ZK2(I)/ZK1(I) P(I+1)=P(I+1)-ZK2(I)*P(I)/ZK1(I)10 CONTINUE P(NJ)=P(NJ)/ZK1(NJ) DO 20 I=1,NJ-1 P(NJ-I)=(P(NJ-I)-ZK2(NJ-I)*P(NJ-I+1)/Z

26、K1(NJ-I)20 CONTINUE ENDc *c 输出位移c *SUBROUTINE SCWY(NJ,P) DIMENSION P(NJ) WRITE(9,10)10 FORMAT(1X,：位 移 ：；：) DO 100 I=1,NJ WRITE(9,20)I,P(I)20 FORMAT(1X,结点号=,I2,5X,F10.4)100 CONTINUE ENDc *c 计算单元杆端内力c * SUBROUTINE GDL(NE,NJ,NF,E0,P,WY,F,F0,DK,PF,GC,GX,GDM) DIMENSION P(NJ),WY(2),F(2),F0(2),DK(2,2),PF(N

27、F,5),GDM(NE,2) DIMENSION GC(NE),GX(NE) DO 200 NE0=1,NE CALL DG (NE0,NE,E0,GC,GX,DK) WY(1)=P(NE0) WY(2)=P(NE0+1) DO 35 I=1,2 F(I)=0.0 DO 30 J=1,2 F(I)=F(I)+DK(I,J)*WY(J)30 CONTINUE35 CONTINUE IF(NF.GT.0) THEN DO 50 I=1,NF IF(INT(PF(I,4)+0.1).EQ.NE0) THEN CALL DJH(I,NE,NF,PF,F0,GC) DO 40 J=1,2 F(J)=F(

28、J)-F0(J)40 CONTINUE END IF50 CONTINUE END IFGDM(NE0,1)=F(1)GDM(NE0,2)=F(2)200 CONTINUE ENDc *c 输出杆端力c *SUBROUTINE SCGDL(GDM,NE) DIMENSION GDM(NE,2)DO 210 I=1,NE WRITE(9,150) I,GDM(I,1),GDM(I,2) 150 FORMAT(1X,单元号=,I2,5X,左端弯矩=,F9.3,2X,右端弯矩=,F9.3)210CONTINUE ENDC *c 计算指定截面的弯矩C *subroutine JMM(NF,GC,GDM

29、,N,X,PF,NE,WXZ,WXY)dimension GC(NE),GDM(NE,2),PF(NF,5)BL=GC(N)A=-GDM(N,1)B=GDM(N,2)WXZ=(B-A)*X/BL+AWXY=(B-A)*X/BL+A DO 50 I=1,NFIF(INT(PF(I,4)+0.1).EQ.N) THENBL=GC(N)G=PF(I,1)Q=PF(I,2) Z=PF(I,3)ID=INT(PF(I,5)+0.1)goto(10,20,30),IDC 集中力作用下的截面弯矩 10 IF(X=Q)THENWXZ=WXZ-G*X*(1-Q/BL)WXY=WXY-G*X*(1-Q/BL)el

30、seWXZ=WXZ-G*Q*(BL-X)/BLWXY=WXY-G*Q*(BL-X)/BLEND IFGOTO 50C 集中力偶作用下的截面弯矩20IF(XQ)THENWXZ=WXZ-G*(BL-X)/BLWXY=WXY-G*(BL-X)/BLELSEWXZ=WXZ+G*X/BLWXY=WXY-G*(BL-X)/BLEND IFEND IFGOTO 50C 均布荷载作用下的截面弯矩30IF(X=Q)THEN WXZ=WXZ-(G*Z-G*Z*Z/(2*BL)*X+(G*Q-G*Q*Q/2/BL)*XWXY=WXY-(G*Z-G*Z*Z/(2*BL)*X+(G*Q-G*Q*Q/2/BL)*Xels

31、eIF(XZ)THENWXZ=WXZ+G*X*X/2-(G*Z-G*Z*Z/2/BL)*X+G*Q*Q*(BL-X)/2/BLWXY=WXY+G*X*X/2-(G*Z-G*Z*Z/2/BL)*X+G*Q*Q*(BL-X)/2/BLelseWXZ=WXZ+G*Z*Z*(X-BL)/2/BL+G*Q*Q*(BL-X)/2/BLWXY=WXY+G*Z*Z*(X-BL)/2/BL+G*Q*Q*(BL-X)/2/BLend ifend if end if50 CONTINUEENDc *c 输出截面弯矩c *SUBROUTINE SCWX(WXZ,WXY)write(9,250) WXZ,WXY250

32、FORMAT(1X,所求截面左侧弯矩=,F7.3,所求截面的右侧弯矩=,F7.3) ENDC *c 计算指定截面的剪力C *subroutine JMF(NF,GC,GDM,N,X,PF,NE,FXZ,FXY)dimension GC(NE),GDM(NE,2),PF(NF,5)BL=GC(N)A=-GDM(N,1)B=GDM(N,2)FXZ=(B-A)/BLFXY=(B-A)/BL DO 50 I=1,NFIF(INT(PF(I,4)+0.1).EQ.N) THENBL=GC(N)G=PF(I,1)Q=PF(I,2) Z=PF(I,3)ID=INT(PF(I,5)+0.1)goto(10,2

33、0,30),IDC 集中力作用下的截面剪力 10 IF(XQ)THENFXZ=FXZ+G*Q/BLFXY=FXY+G*Q/BL ELSE FXZ=FXZ-G*(1-Q/BL) FXY=FXY+G*Q/BL END IFEND IFGOTO 50C 集中力偶作用下的截面剪力 20FXZ=FXZ+G/BL FXY=FXY+G/BLGOTO 50C 均布荷载作用下的截面剪力30IF(X=Q)THENFXZ=FXZ+G*Z*Z/2/BL-G*Z-G*Q*Q/2/BL+G*QFXY=FXY+G*Z*Z/2/BL-G*Z-G*Q*Q/2/BL+G*QelseIF(XNE)THEN FX2=0 RETURN

34、ELSEBL=GC(N)A=-GDM(N,1)B=GDM(N,2)FX2=(B-A)/BL DO 50 I=1,NFIF(INT(PF(I,4)+0.1).EQ.N) THENBL=GC(N)G=PF(I,1)Q=PF(I,2) Z=PF(I,3)ID=INT(PF(I,5)+0.1)goto(10,20,30),IDC 集中力作用下的截面剪力 10 FX2=FX2-G*(1-Q/BL)GOTO 50C 集中力偶作用下的截面剪力 20FX2=FX2+G/BLGOTO 50C 均布荷载作用下的截面剪力30FX2=FX2-G*Z+G*Z*Z/2/BL+G*Q-G*Q*Q/2/BL end if50

35、 CONTINUE END IF ENDC *c 计算支座反力C *subroutine ZZFL(NF,GC,GDM,K,PF,NE)dimension GC(NE),GDM(NE,2),PF(NF,5)CALL ZZZ(NF,GC,GDM,K,PF,NE,FX1) CALL ZZY(NF,GC,GDM,K,PF,NE,FX2) FX=FX1+FX2WRITE(9,150)FX150FORMAT(1X,所求支座反力,F9.3) ENDC *c 计算指定截面弯矩的影响线 C *subroutineMIL(GC,NE,NJ,E0,GX,M,XL,NZ)dimension GC(NE),PF1(5

36、),GX(NE),PJ(20,2),DK(2,2),P(45),F0(2),WY(2)dimension F(2),ZK1(45),ZK2(45),GDM(20,2)DO 70 NN=1,NE write(9,240) NN240 FORMAT(1X,单位荷载作用的单元号为,I2) DO 60 I=0,20 IF(I=0) THEN WX=0 write(9,250) WX250 FORMAT(1X,所求截面的弯矩=,F9.3) ELSE IF(I=20) THEN WX=0 write(9,260) WX260 FORMAT(1X,所求截面的弯矩=,F9.3) ELSE NP=0 NF=1

37、BL=GC(NN) XX=I*BL/20 PF1(1)=-1 PF1(2)=XX PF1(3)=XX PF1(4)=NN PF1(5)=1 CALL XCP(NJ,NP,NF,NE,P,PJ,PF1,F0,GC) CALL JCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX) CALL ZCCL(NZ,NJ,ZK1,ZK2,P) CALL FCQJ(NJ,ZK1,ZK2,P) CALL GDL(NE,NJ,NF,E0,P,WY,F,F0,DK,PF1,GC,GX,GDM) CALL JMM(NF,GC,GDM,M,XL,PF1,NE,WXZ,WXY) write(9,900) WXZ

38、900 FORMAT(1X,所求截面弯矩=,F7.3) END IF END IF 60 CONTINUE 70CONTINUEEND C *c 计算指定截面剪力的影响线C *subroutineFIL(GC,NE,NJ,E0,GX,M,XL,NZ) dimension GC(NE),PF1(5),GX(NE),PJ(20,2),DK(2,2),P(45),F0(2),WY(2)dimension F(2),ZK1(45),ZK2(45),GDM(20,2)DO 70 NN=1,NE write(9,240) NN240 FORMAT(1X,单位荷载作用的单元号为,I2) DO 60 I=0,

39、20 IF(I=0) THEN WX=0 write(9,250) WX250 FORMAT(1X,所求截面的剪力=,F9.3) ELSE IF(I=20) THEN WX=0 write(9,260) WX260 FORMAT(1X,所求截面的剪力=,F9.3) ELSE NP=0 NF=1 BL=GC(NN) XX=I*BL/20 PF1(1)=-1 PF1(2)=XX PF1(3)=XX PF1(4)=NN PF1(5)=1 CALL XCP(NJ,NP,NF,NE,P,PJ,PF1,F0,GC) CALL JCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX) CALL ZCCL(NZ,NJ,ZK1,ZK2,P) CALL FCQJ(NJ,ZK1,ZK2,P) CALL GDL(NE,NJ,NF,E0,P,WY,F,F0,DK,PF1,GC,G