電算實習(xí)報告word版-wenkub.com
2025-03-20 03:29 本頁面
【正文】 ,10x,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 ******************************************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)=10 CONTINUE15 CONTINUEC DL—桿長,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分別存主對角元素和付對角元素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,NE,P,PF,F0,GC) DIMENSION P(NJ),PF(1,2),F0(2),GC(NE) DO 10 I=1,NJ P(I)=10 CONTINUE CALL DJH(NE,PF,F0,GC) ND=INT(PF(1,2)+) 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)= P(1)= ZK2(1)= IF() GOTO 30 GOTO 10030 ZK1(NJ)= ZK2(NJ1)= P(NJ)=100 CONTINUE ENDc ******************************c 解方程求節(jié)點位移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 ******************************c 計算單元桿端內(nèi)力c ***************** SUBROUTINE GDL(NE,NJ,E0,P,WY,F,F0,DK,PF,GC,GX,gd) DIMENSION P(NJ),WY(2),F(2),F0(2),DK(2,2),PF(1,2) DIMENSION GC(NE),GX(NE),gd(ne,2) 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(INT(PF(1,2)+).) THEN CALL DJH(NE,PF,F0,GC) DO 40 J=1,2 F(J)=F(J)F0(J)40 CONTINUE END IF50 CONTINUE gd(ne0,1)=f(1) gd(ne0,2)=f(2)200 CONTINUE endc *********************************************c 計算任意截面的彎矩和剪力c ************************ subroutine jmnl(gd,sqjm,njm,qdnl,gc,ne,pf) dimension gd(ne,2),sqjm(njm,2),gc(ne),qdnl1(1,njm,2),pf(1,2),qdnl(, 1njm,2)c qdnl(njm,1)=彎矩 qdnl(njm,2)=剪力c 取左側(cè)梁為隔離體，彎矩逆時針為正，剪力向上為正， do njm0=1,njm 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)) c=pf(1,1) d=gc(dyh)c if((abs(pf(1,2)dyh)).and.(pf(1,1)=sqjm(njm0,2))) then qdnl1(1,njm0,2)=*d
研究報告相關(guān)推薦
鄂ICP備17016276號-1