電算實(shí)習(xí)報(bào)告

1、結(jié)構(gòu)力學(xué)電算實(shí)習(xí)報(bào)告姓名:* 學(xué)號(hào):*摘 要: 現(xiàn)實(shí)生活中往往存在一些很復(fù)雜的連續(xù)梁模型，它們的支撐性質(zhì)和受力特點(diǎn)都比較復(fù)雜，結(jié)構(gòu)力學(xué)電算則是將它們有規(guī)律的進(jìn)行簡(jiǎn)化，最后將其所有已知的性質(zhì)表達(dá)為一組數(shù)據(jù)，將這些數(shù)據(jù)輸入到我們所編好的程序中，既可以得到我們所要求一些量如：桿端力，各截面內(nèi)力，某界面的內(nèi)力影響線等。關(guān)鍵詞：連續(xù)梁、電算程序、截面內(nèi)力程序功能：連續(xù)梁程序可對(duì)連續(xù)梁結(jié)構(gòu)進(jìn)行相關(guān)的計(jì)算。能夠計(jì)算在任意形式均布荷載、跨中集中力荷載、跨中力偶荷載、節(jié)點(diǎn)力偶荷載組合下連續(xù)梁?jiǎn)卧臈U端力，也可計(jì)算連續(xù)梁各截面內(nèi)力，從而繪制連續(xù)梁的彎矩圖。影響線程序可以計(jì)算單位1集中力荷載作用下的指定截面中內(nèi)力在

2、單位1載荷按指定距離移動(dòng)情況下的變化情況，從而提取數(shù)據(jù)繪制影響線。使用說明：（1）連續(xù)梁內(nèi)力程序在連續(xù)梁計(jì)算程序（.exe文件）所在文件夾中建立一新的txt文件，并命名（1.txt）。在文件其中依次輸入單元數(shù)、支承類型、節(jié)點(diǎn)荷載個(gè)數(shù)、非節(jié)點(diǎn)荷載個(gè)數(shù)、指定截面?zhèn)€數(shù)、彈性模量、桿長(zhǎng)&慣性矩（依次按組輸入）、節(jié)點(diǎn)荷載大小&對(duì)應(yīng)節(jié)點(diǎn)編號(hào)（依次按組輸入）、非節(jié)點(diǎn)荷載值&距該單元左端距離&單元號(hào)&荷載類型號(hào)（依次按組輸入）、所求截面所在單元號(hào)&距該單元左端距離（依次按組輸入），然后保存該文件。運(yùn)行連續(xù)梁內(nèi)力程序，按提示輸入起始文件即初始數(shù)據(jù)的文件（要求全名）

3、，再輸入結(jié)果文件名，在相應(yīng)文件夾下打開結(jié)果文件即可查看結(jié)果。（2）影響線程序在影響線程序（.exe文件）所在文件夾中建立一新的txt文件，并命名。在其中一次輸入單元數(shù)、支承類型、要求影響線的數(shù)量、彈性模量、單位載荷每次移動(dòng)的距離、桿長(zhǎng)&慣性矩、截面所在單元號(hào)&距該單元左端距離，然后保存該文件。運(yùn)行影響線程序，輸入起始文件即初始數(shù)據(jù)的文件，再輸入結(jié)果文件名，在相應(yīng)文件夾下打開結(jié)果文件即可查看結(jié)果。提取期中的數(shù)據(jù)，利用excel或其它工具繪制影響線即可。程序注釋：支承類型：兩端鉸支為1，左固右鉸為2，左鉸右固為3，兩端固定為4。非節(jié)點(diǎn)荷載類型：均布荷載為1，集中力為2，集中力偶為3

4、；節(jié)點(diǎn)荷載只有外力偶。非節(jié)點(diǎn)荷載作用距離：指到單元左端的距離，其中均布荷載指均布荷載結(jié)束位置到單元左端的距離。荷載正負(fù)號(hào)：力偶逆時(shí)針為正，力向上為正。算例1：求圖示結(jié)構(gòu)的彎矩圖及指定截面內(nèi)力?！咀⑨尅繄D中形如 的粗線所示斷面為任意斷面，程序會(huì)求出任意斷面的剪力及彎矩根據(jù)此連續(xù)梁的計(jì)算簡(jiǎn)圖，將其性質(zhì)及受力特征表現(xiàn)為一組數(shù)據(jù)，并輸入運(yùn)行程序進(jìn)行計(jì)算，如下：輸入數(shù)據(jù):4 4 2 5 7 13 6 5 10 4 8 6 1210 3 4 4-6 3 1 1-5 3 2 2-4 3 3 14 1 3 1-12 3 4 31 1.5 2 1 2 3 3 1 3 3 4 2.99 4 3.01得出如下數(shù)據(jù)結(jié)

5、果輸出結(jié)果: *連續(xù)梁內(nèi)力計(jì)算* 單元數(shù)= 4 支承類型= 4 節(jié)點(diǎn)荷載個(gè)數(shù)= 2 非節(jié)點(diǎn)荷載個(gè)數(shù)= 5 指定截面?zhèn)€數(shù)= 7 彈性模量= 1.0000 桿長(zhǎng),慣性矩 GC(NE),GX(NE) 3.000 6.000 5.000 10.000 4.000 8.000 6.000 12.000 節(jié)點(diǎn)荷載大小,對(duì)應(yīng)未知數(shù)序號(hào) PJ(I,1),PJ(I,2) 10.000 3.000 4.000 4.000 非結(jié)點(diǎn)荷載值,距離,單元號(hào),荷載類型號(hào) -6.000 3.000 1.000 1.000 -5.000 3.000 2.000 2.000 -4.000 3.000 3.000 1.000 4.

6、000 1.000 3.000 1.000 -12.000 3.000 4.000 3.000 截面所在單元號(hào)，距離 1.000 1.500 2.000 1.000 2.000 3.000 3.000 1.000 3.000 3.000 4.000 2.990 4.000 3.010 ：位 移 ：；： 結(jié)點(diǎn)號(hào)= 1 .0000 結(jié)點(diǎn)號(hào)= 2 .0109 結(jié)點(diǎn)號(hào)= 3 .4815 結(jié)點(diǎn)號(hào)= 4 .5463 結(jié)點(diǎn)號(hào)= 5 .0000 .各單元桿端內(nèi)力. 單元號(hào)= 1 左端彎矩= 4.543 右端彎矩= -4.413 單元號(hào)= 2 左端彎矩= 4.413 右端彎矩= .296 單元號(hào)= 3 左端彎矩

7、= 9.704 右端彎矩= 2.630 單元號(hào)= 4 左端彎矩= 1.370 右端彎矩= -.815 ;所求截面內(nèi)力; 截面 1 彎矩= 2.272 剪力= -.043 截面 2 彎矩= -1.471 剪力= -2.942 截面 3 彎矩= 4.412 剪力= -2.942 截面 4 彎矩= -2.621 剪力= -7.083 截面 5 彎矩= 3.546 剪力= .917 截面 6 彎矩= -7.073 剪力= 1.907 截面 7 彎矩= 4.888 剪力= 1.907= 計(jì)算結(jié)束 =根據(jù)此得出數(shù)據(jù)，繪制連續(xù)梁彎矩圖：如下：算例2：求圖示結(jié)構(gòu)1、2截面的剪力影響線。輸入數(shù)據(jù)：3 1 2 1

8、 14 4 5 5 6 61 2 2 1輸出結(jié)果： *連續(xù)梁彎矩剪力影響線計(jì)算* 單元數(shù)= 3 支承類型= 1 要求影響線數(shù)量= 2 彈性模量= 1.0000 單位荷載每次移動(dòng)距離 1.0000 桿長(zhǎng),慣性矩 GC(NE),GX(NE) 4.000 4.000 5.000 5.000 6.000 6.000 截面所在單元號(hào)，距離 1.000 2.000 2.000 1.000 單位力所在單元號(hào) 1 距離 1.0000 截面 1 彎矩= .375 剪力= .313 截面 2 彎矩= -.187 剪力= -.062 單位力所在單元號(hào) 1 距離 2.0000 截面 1 彎矩= .800 剪力= -.

9、400 截面 2 彎矩= -.300 剪力= -.100 單位力所在單元號(hào) 1 距離 3.0000 截面 1 彎矩= .325 剪力= -.163 截面 2 彎矩= -.262 剪力= -.087 單位力所在單元號(hào) 1 距離 4.0000 截面 1 彎矩= .000 剪力= .000 截面 2 彎矩= .000 剪力= .000 單位力所在單元號(hào) 2 距離 1.0000 截面 1 彎矩= -.160 剪力= .080 截面 2 彎矩= .512 剪力= -.832 單位力所在單元號(hào) 2 距離 2.0000 截面 1 彎矩= -.200 剪力= .100 截面 2 彎矩= .216 剪力= -.

10、616 單位力所在單元號(hào) 2 距離 3.0000 截面 1 彎矩= -.160 剪力= .080 截面 2 彎矩= .064 剪力= -.384 單位力所在單元號(hào) 2 距離 4.0000 截面 1 彎矩= -.080 剪力= .040 截面 2 彎矩= .008 剪力= -.168 單位力所在單元號(hào) 2 距離 5.0000 截面 1 彎矩= .000 剪力= .000 截面 2 彎矩= .000 剪力= .000 單位力所在單元號(hào) 3 距離 1.0000 截面 1 彎矩= .051 剪力= -.025 截面 2 彎矩= .000 剪力= .102 單位力所在單元號(hào) 3 距離 2.0000 截面

11、 1 彎矩= .074 剪力= -.037 截面 2 彎矩= .000 剪力= .148 單位力所在單元號(hào) 3 距離 3.0000 截面 1 彎矩= .075 剪力= -.038 截面 2 彎矩= .000 剪力= .150 單位力所在單元號(hào) 3 距離 4.0000 截面 1 彎矩= .059 剪力= -.030 截面 2 彎矩= .000 剪力= .119 單位力所在單元號(hào) 3 距離 5.0000 截面 1 彎矩= .032 剪力= -.016 截面 2 彎矩= .000 剪力= .065 單位力所在單元號(hào) 3 距離 6.0000 截面 1 彎矩= .000 剪力= .000 截面 2 彎矩

12、= .000 剪力= .000= 計(jì)算結(jié)束 =1截面彎矩影響線1截面剪力影響線截面2的彎矩剪力影響線同理，不再做出。附：程序清單（程序中含有注釋）（1）連續(xù)梁內(nèi)力程序 PROGRAM LXL DIMENSION GC(20),GX(20),PJ(20,2),PF(10,4),sqjm(20,2),DK(2,2),P(45), 1,F0(2),WY(2),F(2),ZK1(45),ZK2(45),gd(20,2),qdnl(20,2) character*12 dat0,dat1 write(*,*)'輸入:初始數(shù)據(jù)文件名' read(*,10)dat0 write(*,*)&#

13、39;輸入:計(jì)算結(jié)果數(shù)據(jù)文件名' read(*,10)dat110 format(A12) open(8,'old') open(9,'new') write(9,11)11 format(8x,'*連續(xù)梁內(nèi)力計(jì)算*') READ(8,*)NE,NZ,NP,NF,njm,E0 WRITE(9,12)NE,NZ,NP,NF,njm,E0 12 FORMAT(1X,'單元數(shù)=',6x,I3,5x,'支承類型=',I3,5x,'節(jié)點(diǎn)荷載個(gè)數(shù)='I3,/, 1 1x,'非節(jié)點(diǎn)荷載個(gè)數(shù)=

14、9;,I3,5x,'指定截面?zhèn)€數(shù)=',i3,5x,'彈性模量=',F12., 14) NJ=NE+1 CALL SRSJ(NE,NP,NF,GC,GX,PJ,PF,sqjm,njm) 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,gd) call

15、jmnl(gd,sqjm,njm,qdnl,gc,ne,nf,pf)ENDC *c 輸入初始數(shù)據(jù)SRSJ子程序C * SUBROUTINE SRSJ(NE,NP,NF,GC,GX,PJ,PF,sqjm,njm) DIMENSION GC(NE),GX(NE),PJ(NP,2),PF(NF,4),sqjm(njm,2)csqjm:所求截面 njm:截面數(shù)目 sqjm(njm,1)=所屬單元 sqjm(njm,2)=距左端距離C GC(NE)：桿長(zhǎng) GX(NE)：慣性矩C 結(jié)點(diǎn)荷載:PJ(NP,1)=荷載大小 PJ(NP,2)=對(duì)應(yīng)未知數(shù)序號(hào) C 非結(jié)點(diǎn)荷載 : PF(NE,1)=荷載大小 PF(

16、NE,2)=距離 PF(NE,3)=單元號(hào) PF(NE,4)=荷載類型號(hào) WRITE(9,11)11 FORMAT(1X,'桿長(zhǎng),慣性矩 GC(NE),GX(NE)') READ(8,*)(GC(I),GX(I),I=1,NE) write(9,12)(GC(I),GX(I),I=1,NE)12 FORMAT(10X,2F8.3) IF(NP.GT.0)THEN WRITE(9,22)22 FORMAT(1X,'節(jié)點(diǎn)荷載大小,對(duì)應(yīng)未知數(shù)序號(hào) PJ(I,1),PJ(I,2)') READ(8,*)(PJ(I,1),PJ(I,2),I=1,NP) write(9,2

17、3)(PJ(I,1),PJ(I,2),I=1,NP)23 FORMAT(1X,2F8.3) END IF IF(NF.GT.0) THEN WRITE(9,33) 33 FORMAT(1X,'非結(jié)點(diǎn)荷載值,距離,單元號(hào),荷載類型號(hào)') READ(8,*)(PF(I,J),J=1,4),I=1,NF) WRITE(9,34)(PF(I,J),J=1,4),I=1,NF) 34 FORMAT(1X,4F8.3) END IFif(njm.gt.0) thenwrite(9,35)35 format(1x,'截面所在單元號(hào)，距離') read(8,*)(sqjm(i,

18、j),j=1,2),i=1,njm)write(9,36)(sqjm(i,j),j=1,2),i=1,njm)36 format(1x,2f8.3) end if END C *#*c 計(jì)算第NHF個(gè)非結(jié)點(diǎn)荷截引起的等效結(jié)點(diǎn)荷截列陣F0C * SUBROUTINE DJH(NHF,NE,NF,PF,F0,GC) DIMENSION PF(NF,4),F0(2),GC(NE)C G荷載值,C作用距離,NT-單元號(hào),ID荷載類型號(hào) G=PF(NHF,1) C=PF(NHF,2) NT=INT(PF(NHF,3)+0.1) ID=INT(PF(NHF,4)+0.1) BL=GC(NT) D=BL-C

19、 C1=C/BL C2=C1*C1 GOTO(10,20,30),IDC 均布荷載作用下的等效荷載列陣F0(2)10 F0(1)=G*C*C*(6.0-8.0*C1+3.0*C2)/12.0 F0(2)=-G*C*C*C*(4.0-3.0*C1)/12.0/BL GOTO 200C 集中力作用下的等效荷截列陣P0(2)20 F0(1)=G*C*D*D/BL/BL F0(2)=-G*D*C2goto 200C跨中力偶作用下的等效荷載列陣Q0(2) 30F0(1)=+(6.0*G*c*c/bl/bl-6.0*g*c/bl+(2.0*g*c*bl-3.0*g*c*c)/bl/bl+g)F0(2)=-

20、(2.0*g*c*bl-3.0*g*c*c)/bl/bl200 RETURN ENDC *c 計(jì)算第NE0個(gè)單元的單元?jiǎng)偠染仃嘍K(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 DL桿長(zhǎng),DI-慣性矩,S-線剛度 DL=GC(NE0) DI=GX(NE0) S=E0*DI/DL DK(1,1)=4.0*S DK(1,2)=2.0*S DK(2,1)=2.0*S DK(2,2)=

21、4.0*S ENDc *C 集成總體剛度矩陣,ZKl,ZK2分別存主對(duì)角元素和付對(duì)角元素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)=ZK2(I)+DK(1,2) ZK1(I+1)=ZK1(I+1)+DK(2,2)100 CO

22、NTINUE ENDC *C 形成總荷裁矩陣C * SUBROUTINE XCP(NJ,NP,NF,NE,P,PJ,PF,F0,GC) DIMENSION P(NJ),PJ(NP,2),PF(NF,4),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 NF0=1,NF CALL DJH(NF0,NE,NF,PF,F0,GC) ND=INT(PF(

23、NF0,3)+0.1) P(ND)=P(ND)+F0(1) P(ND+1)=P(ND+1)+F0(2)100 CONTINUE END IF ENDc *c 進(jìn)行支承條件處理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 ZK1(NJ)=1.0 ZK2(NJ-1)=0.0 P(NJ)=0.0100 CONTINUE

24、ENDc *c 解方程求節(jié)點(diǎn)位移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)/ZK1(NJ-I)20 CONTINUE ENDc *c 輸出位移c * SUBROUTINE

25、 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,'結(jié)點(diǎn)號(hào)=',I2,5X,F10.4)100 CONTINUE ENDc *c 計(jì)算單元桿端內(nèi)力c * SUBROUTINE GDL(NE,NJ,NF,E0,P,WY,F,F0,DK,PF,GC,GX,gd) DIMENSION P(NJ),WY(2),F(2),F0(2),DK(2,2),PF(NF,4) DIMENSION GC(NE),GX(N

26、E),gd(ne,2) WRITE(9,2)2 FORMAT(1X,'.各單元桿端內(nèi)力.') 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,3)+0.1).EQ.NE0) THEN CALL DJH(I,NE,NF,PF,F0,GC) DO

27、40 J=1,2 F(J)=F(J)-F0(J)40 CONTINUE END IF50 CONTINUE END IF WRITE(9,150) NE0,F(1),F(2)gd(ne0,1)=f(1)!將桿端力數(shù)據(jù)存入內(nèi)存 以備調(diào)用gd(ne0,2)=f(2)150 FORMAT(1X,'單元號(hào)=',I2,5X,'左端彎矩=',F9.3,2X,'右端彎矩=',F9.3)200 CONTINUEendc *c計(jì)算任意截面的彎矩和剪力（利用疊加法）c *subroutine jmnl(gd,sqjm,njm,qdnl,gc,ne,nf,pf)dim

28、ension gd(ne,2),sqjm(njm,2),gc(ne),qdnl1(nf,njm,2),pf(nf,4),qdn，1l(njm,2)c qdnl(njm,1)=彎矩 qdnl(njm,2)=剪力c 取左側(cè)梁為隔離體，彎矩逆時(shí)針為正，剪力向上為正，write(9,20)20format(1x,'所求截面內(nèi)力;')do njm0=1,njm!按要求的截面順序循環(huán)dyh=sqjm(njm0,1)qdnl(njm0,2)=-(gd(dyh,2)+gd(dyh,1)/gc(dyh)!單元上無荷載情況qdnl(njm0,1)=-(gd(dyh,1)+qdnl(njm0,2)*

29、sqjm(njm0,2)do nf0=1,nfc=pf(nf0,2)d=gc(dyh)-cif(abs(pf(nf0,3)-dyh)<1.0e-6).and.(pf(nf0,2)>=sqjm(njm0,2) thenif(pf(nf0,4)=1)then !單元上有均布荷載且荷載結(jié)束位置在所求截面左側(cè)qdnl1(nf0,njm0,2)=pf(nf0,1)*c*(2.0*gc(dyh)-c)/2.0/gc(dyh)-pf(nf0,1,1)*sqjm(njm0,2)qdnl1(nf0,njm0,1)=-(pf(nf0,1)*c*(2.0*gc(dyh)-c)/2.0/gc(dyh)*s

30、qjm(n,1jm0,2)-pf(nf0,1)*sqjm(njm0,2)*sqjm(njm0,2)/2.0)else if(pf(nf0,4)=2)then !單元上有集中力荷載且位置在所求截面左側(cè)qdnl1(nf0,njm0,2)=pf(nf0,1)*d/gc(dyh)qdnl1(nf0,njm0,1)=-qdnl1(nf0,njm0,2)*sqjm(njm0,2)else if(pf(nf0,4)=3)then !單元上有集中力偶且位置在所求截面左側(cè)qdnl1(nf0,njm0,2)=-pf(nf0,1)/gc(dyh)qdnl1(nf0,njm0,1)=-qdnl1(nf0,njm0,2

31、)*sqjm(njm0,2)end if else if(abs(pf(nf0,3)-dyh)<1.0e-6).and.(pf(nf0,2)<sqjm(njm0,2) &&thenif(pf(nf0,4)=1)then !單元上有均布荷載且荷載結(jié)束位置在所求截面右側(cè)qdnl1(nf0,njm0,2)=-pf(nf0,1)*c*c/2.0/gc(dyh)qdnl1(nf0,njm0,1)=-pf(nf0,1)*c*c/2.0/gc(dyh)*(gc(dyh)-sqjm(njm0,2),1)else if(pf(nf0,4)=2)then !單元上有集中力荷載且位置在所

32、求截面左側(cè)qdnl1(nf0,njm0,2)=-pf(nf0,1)*c/gc(dyh)qdnl1(nf0,njm0,1)=-pf(nf0,1)*d*c/gc(dyh)-qdnl1(nf0,njm0,2)*(sqjm(n,1jm0,2)-c)else if(pf(nf0,4)=3)then !單元上有集中力偶且位置在所求截面左側(cè)qdnl1(nf0,njm0,2)=-pf(nf0,1)/gc(dyh)qdnl1(nf0,njm0,1)=pf(nf0,1)/gc(dyh)*(sqjm(njm0,2)-gc(dyh)end ifelseqdnl1(nf0,njm0,2)=0 !其他情況下 需疊加的內(nèi)力

33、為0qdnl1(nf0,njm0,1)=0end ifqdnl(njm0,1)=qdnl(njm0,1)+qdnl1(nf0,njm0,1)!疊加求內(nèi)力qdnl(njm0,2)=qdnl(njm0,2)+qdnl1(nf0,njm0,2)end do write(9,30) njm0,qdnl(njm0,1),qdnl(njm0,2)30format(1x,'截面',i3,2x,'彎矩=',f9.3,2x,'剪力=',f9.3)end do WRITE(9,11)11 FORMAT('= 計(jì)算結(jié)束 =')end（2）影響線程序 P

34、ROGRAM LXLyxx DIMENSION GC(20),GX(20),PF(1,2),sqjm(20,2),DK(2,2),P(45),F0(2), WY,1(2),F(2),ZK1(45),ZK2(45),gd(20,2),qdnl(20,2),fs(20),fs0(20) character*12 dat0,dat1 write(*,*)'輸入:初始數(shù)據(jù)文件名' read(*,10)dat0 write(*,*)'輸入:計(jì)算結(jié)果數(shù)據(jù)文件名' read(*,10)dat110 format(A12) open(8,'old') open(

35、9,'new') write(9,11)11 format(8x,'*連續(xù)梁彎矩剪力影響線計(jì)算*') READ(8,*)NE,NZ,njm,E0,hz WRITE(9,12)NE,NZ,njm,E0,hz 12 FORMAT(1X,'單元數(shù)=',6x,I3,5x,'支承類型=',I3,5x,'要求影響線數(shù)量=', 1i3,5x,'彈性模量=',F12.4,6x,'單位荷載每次移動(dòng)距離',f12.4,6x) NJ=NE+1 CALL SRSJ(NE,GC,GX,sqjm,njm) do

36、 ne0=1,ne!讓單位載荷按指定距離移動(dòng)的循環(huán)fs(ne0)=gc(ne0)/hzfs0(ne0)=int(fs(ne0)+0.1)end done0=0do ne0=1,nepf(1,2)=ne0do fs00=1,fs0(ne0)pf(1,1)=hz*fs00write(9,14) ne0,pf(1,1)14format(/8x,'單位力所在單元號(hào)',i3,12x,'距離',f8.4) CALL XCP(NJ,NE,P,PF,F0,GC) CALL JCZK(NE,NJ,E0,ZK1,ZK2,DK,GC,GX) CALL ZCCL(NZ,NJ,ZK1,Z

37、K2,P) CALL FCQJ(NJ,ZK1,ZK2,P) CALL GDL(NE,NJ,E0,P,WY,F,F0,DK,PF,GC,GX,gd) call jmnl(gd,sqjm,njm,qdnl,gc,ne,pf)end doend do WRITE(9,13)13 FORMAT('= 計(jì)算結(jié)束 =')ENDC *c 輸入初始數(shù)據(jù)SRSJ子程序C * SUBROUTINE SRSJ(NE,GC,GX,sqjm,njm) DIMENSION GC(NE),GX(NE),sqjm(njm,2)csqjm:所求截面 njm:截面數(shù)目 sqjm(njm,1)=所屬單元 sqjm(

38、njm,2)=距左端距離C GC(NE)：桿長(zhǎng) GX(NE)：慣性矩C 非結(jié)點(diǎn)荷載單位1 : PF(1,1)=距離 PF(1,2)=單元號(hào) WRITE(9,11)11 FORMAT(1X,'桿長(zhǎng),慣性矩 GC(NE),GX(NE)') READ(8,*)(GC(I),GX(I),I=1,NE) write(9,12)(GC(I),GX(I),I=1,NE)12 FORMAT(10X,2F8.3)if(njm.gt.0) thenwrite(9,35)35 format(1x,'截面所在單元號(hào)，距離') read(8,*)(sqjm(i,j),j=1,2),i=1

39、,njm)write(9,36)(sqjm(i,j),j=1,2),i=1,njm)36 format(1x,2f8.3) end if END C *#*c 計(jì)算單位1力引起的等效結(jié)點(diǎn)荷截列陣F0C * SUBROUTINE DJH(NE,PF,F0,GC) DIMENSION PF(1,2),F0(2),GC(NE)C C作用距離,NT-單元號(hào) C=PF(1,1) NT=INT(PF(1,2)+0.1) BL=GC(NT) D=BL-C C1=C/BL C2=C1*C1C 集中力作用下的等效荷截列陣P0(2)20 F0(1)=-1.0*C*D*D/BL/BL F0(2)=1.0*D*C2

40、ENDC *c 計(jì)算第NE0個(gè)單元的單元?jiǎng)偠染仃嘍K(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 DL桿長(zhǎng),DI-慣性矩,S-線剛度 DL=GC(NE0) DI=GX(NE0) S=E0*DI/DL DK(1,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分別存主對(duì)角元素和

41、付對(duì)角元素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)=ZK2(I)+DK(1,2) ZK1(I+1)=ZK1(I+1)+DK(2,2)100 CONTINUE ENDC *C 形成總荷裁矩陣C * SUBROUTINE XC

42、P(NJ,NE,P,PF,F0,GC) DIMENSION P(NJ),PF(1,2),F0(2),GC(NE) DO 10 I=1,NJ P(I)=0.010 CONTINUE CALL DJH(NE,PF,F0,GC) ND=INT(PF(1,2)+0.1) P(ND)=P(ND)+F0(1) P(ND+1)=P(ND+1)+F0(2) ENDc *c 進(jìn)行支承條件處理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)=