【正文】 ite(9,20)20 format(1x,39。所求截面內(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)*sqjm(njm0,2)) do nf0=1,nf c=pf(nf0,2) d=gc(dyh)c if((abs(pf(nf0,3)dyh)).and.(pf(nf0,2)=sqjm(njm0,2))) then if(pf(nf0,4)==1)then !單元上有均布荷載且荷載結(jié)束位置在所求截面左側(cè) qdnl1(nf0,njm0,2)=pf(nf0,1)*c*(*gc(dyh)c)/(dyh)pf(nf0,1, 1)*sqjm(njm0,2) qdnl1(nf0,njm0,1)=(pf(nf0,1)*c*(*gc(dyh)c)/(dyh)*sqjm(n, 1jm0,2)pf(nf0,1)*sqjm(njm0,2)*sqjm(njm0,2)/) 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)*sqjm(njm0,2) end if else if((abs(pf(nf0,3)dyh)).and.(pf(nf0,2)sqjm(njm0,2))) amp。then if(pf(nf0,4)==1)then !單元上有均布荷載且荷載結(jié)束位置在所求截面右側(cè) qdnl1(nf0,njm0,2)=pf(nf0,1)*c*c/(dyh) qdnl1(nf0,njm0,1)=pf(nf0,1)*c*c/(dyh)*(gc(dyh)sqjm(njm0,2), 1) else if(pf(nf0,4)==2)then !單元上有集中力荷載且位置在所求截面左側(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 if else qdnl1(nf0,njm0,2)=0 !其他情況下 需疊加的內(nèi)力為0 qdnl1(nf0,njm0,1)=0 end if qdnl(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)30 format(1x,39。,i3,2x,39。,2x,39。,) end do WRITE(9,11)11 FORMAT(39。) end（2）影響線程序 PROGRAM 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(*,*)39。 read(*,10)dat0 write(*,*)39。 read(*,10)dat110 format(A12) open(8,file=dat0,status=39。) open(9,file=dat1,status=39。) write(9,11)11 format(8x,39。) READ(8,*)NE,NZ,njm,E0,hz WRITE(9,12)NE,NZ,njm,E0,hz 12 FORMAT(1X,39。,6x,I3,5x,39。,I3,5x,39。, 1i3,5x,39。,6x,39。,6x) NJ=NE+1 CALL SRSJ(NE,GC,GX,sqjm,njm) do ne0=1,ne !讓單位載荷按指定距離移動的循環(huán) fs(ne0)=gc(ne0)/hz fs0(ne0)=int(fs(ne0)+) end do ne0=0 do ne0=1,ne pf(1,2)=ne0 do fs00=1,fs0(ne0) pf(1,1)=hz*fs00 write(9,14) ne0,pf(1,1)14 format(/8x,39。,i3,12x,39。,) CALL XCP(NJ,NE,P,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 GDL(NE,NJ,E0,P,WY,F,F0,DK,PF,GC,GX,gd) call jmnl(gd,sqjm,njm,qdnl,gc,ne,pf) end do end do WRITE(9,13)13 FORMAT(39。) ENDC ************************************c 輸入初始數(shù)據(jù)SRSJ子程序C ************************************ SUBROUTINE SRSJ(NE,GC,GX,sqjm,njm) DIMENSION GC(NE),GX(NE),sqjm(njm,2)c sqjm:所求截面 njm:截面數(shù)目 sqjm(njm,1)=所屬單元 sqjm(njm,2)=距左端距離C GC(NE)：桿長 GX(NE)：慣性矩C 非結(jié)點荷載單位1 : PF(1,1)=距離 PF(1,2)=單元號 WRITE(9,11)11 FORMAT(1X,39。) READ(8,*)(GC(I),GX(I),I=1,NE) write(9,12)(GC(I),GX(I),I=1,NE)12 FORMAT(10X,) if() then write(9,35)35 format(1x,39。) read(8,*)((sqjm(i,j),j=1,2),i=1,njm) write(9,36)((sqjm(i,j),j=1,2),i=1,njm)36 format(1x,) end if END C **********************************************c 計算單位1力引起的等效結(jié)點荷截列陣F0C *********************************************** SUBROUTINE DJH(NE,PF,F0,GC) DIMENSION PF(1,2),F0(2),GC(NE)C C—作用距離,NT單元號 C=PF(1,1) NT=INT(PF(1,2)+) BL=GC(NT) D=BLC C1=C/BL C2=C1*C1C 集中力作用下的等效荷截列陣P0(2)20 F0(1)=*C*D*D/BL/BL F0(2)=*D*C2 ENDC *********************************