電算實(shí)習(xí)報(bào)告word版(編輯修改稿)
2025-04-19 03:29 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 NUEC DL—桿長(zhǎng),DI慣性矩,S線剛度 DL=GC(NE0) DI=GX(NE0) S=E0*DI/DL DK(1,1)=*S DK(1,2)=*S DK(2,1)=*S DK(2,2)=*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)= ZK2(I)=10 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 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)=10 CONTINUE IF() THEN DO 20 I=1,NP J=INT(PJ(I,2)+) P(J)=PJ(I,1)20 CONTINUE END IF IF() THEN DO 100 NF0=1,NF CALL DJH(NF0,NE,NF,PF,F0,GC) ND=INT(PF(NF0,3)+) 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)= P(1)= ZK2(1)= IF() GOTO 30 GOTO 10030 ZK1(NJ)= ZK2(NJ1)= P(NJ)=100 CONTINUE ENDc ******************************c 解方程求節(jié)點(diǎn)位移Pc ****************************** SUBROUTINE FCQJ(NJ,ZK1,ZK2,P) DIMENSION ZK1(NJ),ZK2(NJ),P(NJ) DO 10 I=1,NJ1 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,NJ1 P(NJI)=(P(NJI)ZK2(NJI)*P(NJI+1))/ZK1(NJI)20 CONTINUE ENDc *******************************c 輸出位移c ******************************* SUBROUTINE SCWY(NJ,P) DIMENSION P(NJ) WRITE(9,10)10 FORMAT(1X,39。：：：：：：：：：位 移 ：；：：：：：：：：39。) DO 100 I=1,NJ WRITE(9,20)I,P(I)20 FORMAT(1X,39。結(jié)點(diǎn)號(hào)=39。,I2,5X,)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(NE),gd(ne,2) WRITE(9,2)2 FORMAT(1X,39。.................各單元桿端內(nèi)力....................39。) 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)= DO 30 J=1,2 F(I)=F(I)+DK(I,J)*WY(J)30 CONTINUE35 CONTINUE IF() THEN DO 50 I=1,NF IF(INT(PF(I,3)+).) THEN CALL DJH(I,NE,NF,PF,F0,GC) DO 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,39。單元號(hào)=39。,I2,5X,39。左端彎矩=39。,2X,39。右端彎矩=39。,)200 CONTINUE endc *********************************************c 計(jì)算任意截面的彎矩和剪力 （利用疊加法）c ************************ subroutine jmnl(gd,sqjm,njm,qdnl,gc,ne,nf,pf) dimension 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)20 format(1x,39。。。所求截面內(nèi)力。。。39。) do njm0=1,njm !按要求的截面順序循環(huán) dyh=sqjm(njm0,1) qdnl(njm0,2)=(gd(dyh,2)+gd(dyh,1))/gc(dyh) !單元上無(wú)荷載情況 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 i
研究報(bào)告相關(guān)推薦
鄂ICP備17016276號(hào)-1