課程學(xué)習(xí)-有限元作業(yè)

平面剛架靜力分析程序例平面剛架，已知各豎直桿截面A=0.2m2，I=4.1667×10-3m4，各水平桿截面A=0.16m2，I=2.1333×103m4，各桿E=3.65×104MPa，求該剛架的各桿端1-1計(jì)算分對該平面剛架進(jìn)行單元?jiǎng)澐?，然后對各結(jié)點(diǎn)、單元進(jìn)行。共13個(gè)單元，13個(gè)節(jié)點(diǎn)坐標(biāo)軸及單元與節(jié)點(diǎn)的如圖將均不荷載分配到各相應(yīng)的節(jié)點(diǎn)上，1-2輸入文 124.166700E-234.166700E-252.133300E-362.133300E-454.166700E-564.166700E-582.133300E-692.133300E-784.166700E-894.166700E-82.133300E-2.133300E-4.166700E-00000000000000003200300-00-734-42-572-584-933輸出文TheInputDataTeralInformationE1TheInformationofAI1122.000000E-4.166700E-2232.000000E-4.166700E-3251.600000E-2.133300E-4361.600000E-2.133300E-5452.000000E-4.166700E-6562.000000E-4.166700E-7581.600000E-2.133300E-8691.600000E-2.133300E-9782.000000E-4.166700E-892.000000E-4.166700E-81.600000E-2.133300E-1.600000E-2.133300E-2.000000E-4.166700E-TheJointXY123456789TheInformationof(258(927TheLoadingsat323--TheLoadingsat 34-42-72-84-933TheResultsofCalculationTheJointDisplacementsuv19.727862E--7.663514E--5.946894E-28.890688E--1.049796E--7.120790E-31.470285E--2.323802E--7.465713E-49.683458E--3.536778E--5.930754E-58.886289E--4.844902E--7.160653E-61.469823E--7.581843E-71.058868E--2.196870E--6.233836E-88.890698E--3.009411E--6.177454E-91.470137E--3.792985E--1.977170E-8.890698E--3.520155E--7.782789E-8.890698E--7.411550E--7.782790E-TheTerminal1--2----3----4--5--6--7---8--9----(258(951變量說Exy受載類型號(hào)（1為水平集中荷載；23為水平均布荷載；4為豎向均布荷載!PF.FOR(Aprogramforysisofplaneframe)MAINPROGRAM!主程序讀入控制數(shù)據(jù)，分配動(dòng)態(tài)數(shù)組地址，并調(diào)用各子程序來調(diào)用整個(gè)計(jì)DIMENSIONW(80000)CHARACTERIDFN*20,TITLE(5)*72FORMAT(1X,'InputDataName:')FORMAT(1X,'OutputFileName:')OPEN(4,FILE=OUTFN,FORMAT(1X,'InputDataName:',1X,A12)FORMAT(1X,'OutputFileName:',1X,A12)FORMAT(/13X,'TheInputData'//10X,'Teral//6X,'E',9X,'NM',5X,'NJ',5X,'NS',5X,FORMAT(/10X,'TheInformationof//1X,'member',2X,'start',2X,'end',9X,'A',15X,CALLCALLLCVCTCALLFORMACALLDO100LC=1,NLCREAD(3,*)NLJCALLB0(LIF(NLJ.NE.0)CALLIOLJB(N,NLJ,W(L61),W(L62),&CALLIF(NLM.NE.0)CALLIOLMFB(NM,NJ,N,NLM,W(L81),&CALLCALLOJD(NJ,N,W(L54))CALL100CONTINUESUBROUTINE!讀入桿件、節(jié)點(diǎn)及支承條件的數(shù) CHARACTERTITLE(5)*72 READSUBROUTINE!確定單元定位向DO100M=1,NMLV(1,M)=I-LV(2,M)=I-LV(4,M)=J-LV(5,M)=J-100CONTINUESUBROUTINE!確定總綱ADO100I=1,N100CONTINUEDO400M=1,NMDO200I=2,6200CONTINUEDO300300CONTINUE400CONTINUEDO500500CONTINUEDO600I=2,N600CONTINUESUBROUTINE 計(jì)算桿件的桿長、余弦和正DIMENSIONIST(NM),IEN(NM),X(NJ),Y(NJ)X1=X(J)-SUBROUTINEKEBAR 計(jì)算局部坐標(biāo)DIMENSIONIST(NM),IEN(NM),X(NJ),Y(NJ),AR(NM),RI(NM)CALLRLCS(M,NM,NJ,IST,IEN,X,Y,RL,C,S)SUBROUTINEKE 計(jì)算整體坐標(biāo)DIMENSIONAE(4,1)=-AE(5,2)=-SUBROUTINEFORMA 形成總剛DIMENSIONIST(NM),IEN(NM),AR(NM),RI(NM),X(NJ),Y(NJ),&DOUBLEPRECISIONA(NA)DO300M=1,NMDO200I=1,6DO100A(LD(LV(I,M))-A(LD(LV(J,M))-ENDIF100200300CONTINUESUBROUTINEAS!總剛ADIMENSIONIS(NS),LD(N)DOUBLEPRECISIONA(NA)DO100M=1,NS100CONTINUE!解方程組之一總綱分DIMENSIONDOUBLEPRECISIONA(NA),T(N),SUMDO300I=2,NDO200J=I1,I-1IF(I1.GT.J1)J1=I1DO100K=J1,J-1100CONTINUE200300CONTINUE!解方程組之二向前和向后回DIMENSIONDOUBLEPRECISIONA(NA),B(N)DO200I=2,NDO100J=I1,I-1100CONTINUE200DO300I=1,N300DO500I=N-1,1,-1DO400J=I+1,IMINJ1=J-LDJ+LD(J-IF(I.GE.J1)B(I)=B(I)-A(LDJ-J+I)*B(J)400CONTINUE500CONTINUESUBROUTINEB0(L!節(jié)點(diǎn)荷載列陳B置零DOUBLEPRECISIONB(N)WRITE(*,1)LC,NLJ//10X,'TheLoadingsatDO100I=1,N100SUBROUTINEIOLJB!讀入并輸出節(jié)點(diǎn)荷載，形成節(jié)點(diǎn)荷載列陣BDIMENSIONILJ(NLJ),PX(NLJ),PY(NLJ),PM(NLJ)DOUBLEPRECISION WRITE(*,1) WRITE DO100M=1,NLJB(I-B(I-100SUBROUTINE!桿端力置DIMENSIONF(6,NM)WRITE RMAT(/10X,'TheLoadingsatMembers'//17X,'NLM=',I4)DO200J=1,NMDO100I=1,6100200CONTINUESUBROUTINE!考慮非節(jié)點(diǎn)荷載，計(jì)算固端力，計(jì)算等效節(jié)點(diǎn)荷載并疊加到節(jié)點(diǎn)荷載列陣BDIMENSIONILM(NLM),ITL(NLM),PV(NLM),DST(NLM),IST(NM),&DOUBLEPRECISIONB(N)WRITE(*,1) DO100M=1,NLMD2=RL-IF(ITL(M).EQ.1.OR.ITL(M).EQ.3)THENF4=-P1-B(LV(2,L))=B(LV(2,L))-F1*S-B(LV(3,L))=B(LV(3,L))-B(LV(4,L))=B(LV(4,L))-B(LV(5,L))=B(LV(5,L))-F4*S-100CONTINUESUBROUTINE!引入支承條DIMENSIONIS(NS),VS(NS)DOUBLEPRECISIONB(N)DO100M=1,NS100CONTINUESUBROUTINE!輸出節(jié)點(diǎn)位DOUBLEPRECISIONB(N)WRITE(*,1)//10X,'TheJointWRITE(*,2)(M,B(3*M-2),B(3*M- SUBROUTINE!計(jì)算并輸出桿端DIMENSIONIST(NM),IEN(NM),AR(NM),RI(NM),X(NJ),Y(NJ),&DOUBLEPRECISIONB(N),U1,U2,U3,U4,U5,U6WRITE(*,1)FORMAT(/10X,'TheTerminalDO100M=1,NMF(2,M)=F(2,M)+E2*(U2-F(3,M)=F(3,M)+E3*(U2-F(4,M)=F(4,M)+E1*(U4-F(5,M)=F(5,M)+E2*(U5-U2)-100CONTINUE平面三角形三節(jié)點(diǎn)單元有限元分析程序例q1Nmm2作用。E2.1105Nmm2=0.3，厚度h10mm2-1算例分81080個(gè)單元，51155個(gè)節(jié)點(diǎn)，坐標(biāo)軸及單元與節(jié)點(diǎn)的如圖。將均51、52、53、54、55視作固定2-1所示的離散化計(jì)算模型。輸入文400.,80.,400.,60.,400.,40.,0.,-200.,0.,0.,0.,0.,0.,0.,0.,0.,0.,-0.,0.,0.,0.,0.,0.,0.,0.,0.,-0.,0.,0.,0.,0.,0.,0.,0.,0.,-0.,0.,0.,0.,0.,0.,0.,0.,0.,-0.,0.,0.,0.,0.,0.,0.,0.,0.,-0.,0.,0.,0.,0.,0.,0.,0.,0.,-0.,0.,0.,0.,0.,0.,0.,0.,0.,-0.,0.,0.,0.,0.,0.,0.,0.,0.,-0.,0.,0.,0.,0.,0.,0.,0.,0.,-0.,0.,0.,0.,0.,0.,0.,0.,0.,-輸出文 T2 172283394456--------------------------------------6.27E-------------------8.84E--6.76E--7.74E--2.63E--4.72E--2.43E--7.39E--8.58E--1.60E-2.5變量說E節(jié)點(diǎn)位移度T結(jié)構(gòu)位移度總單元位移度總被固定的節(jié)點(diǎn)位移度序各單元節(jié)點(diǎn)2.5PROGRAMDIMENSIONSK(300,30),EK(12,12),Q(300),MC(55),XY(3,100),XYE(3,4),& REWIND7 FORMAT(10X,'NFNENNMBNDLELSEUMREWIND9 STOP1000DIMENSIONXY(LS,NN),Q(MS),NX(NE,LE),MC(ND) DO35DO35 DO200DO40J=1,NEDO40I=1,LS DO50DO50 DO200DO200II=1,NFDO200J=1,NEDO200JJ=1,NF DO220I=1,ND SUBROUTINESOLVE(SK,Q,MS,MB)DIMENSIONSK(MS,MB),Q(MS)DO125IF(K+MB-1- N=K+MB-GOTO DO125I=I1,NDO122 SK(I,J)=SK(I,J)- Q(I)=Q(I)-DO145I=MS-IF(MS-I+1- N=MS-GOTO DO142J=2,N Q(I)=Q(I)- DIMENSIONQ(MS),QE(ME),NX(NE,LE),XY(LS,NN),XYE(LS,NE)DO400DO160DO160J=1,NF DO200J=1,NEDO200I=1,LS CALL SUBROUTINESTIFE(EK,XYE,ME,NE,NF,LS,E,UM,T)DIMENSIONEK(ME,ME),XYE(LS,NE),B(3),C(3)B(2)=XYE(2,3)-B(3)=XYE(2,1)-C(1)=XYE(1,3)-C(2)=XYE(1,1)-C(3)=XYE(1,2)-DO30DO30J=1,3N=2*(J-EK(M,N)=D*(B(I)*B(J)+C(I)*C(J)*(1-EK(M,N+1)=D*(UM*B(I)*C(J)+C(I)*B(J)*(1-EK(M+1,N)=D*(UM*C(I)*B(J)+B(I)*C(J)*(1- SUBROUTINEDIMENSIONXYE(LS,NE),QE(ME),SG(3),B(3),C(3),S(3,6)B(2)=XYE(2,3)-B(3)=XYE(2,1)-C(1)=XYE(1,3)-C(2)=XYE(1,1)-C(3)=XYE(1,2)-DO220I=1,3 S(3,2*I)=D*(1-DO250I=1,3DO250 平面四節(jié)點(diǎn)矩形薄板單元有限元分析程序例題及分1/4E=2.1×105N/mm2=0.3h=20mm3-1將它劃分為4×5的矩形網(wǎng)格，其節(jié)點(diǎn)及單元如圖3-1，建立相應(yīng)的離散輸入文-125000,-10416670,20833330,-250000,--250000,-20833330,0,-250000,--125000,-10416670,-20833330,--500000,0,0,--500000,0,0,-250000,0,--250000,0,41666670,--500000,0,0,--250000,0,-41666670,--500000,0,0,--500000,0,0,-250000,0,--250000,0,41666670,-500000,0,--500000,0,0,-500000,0,--250000,0,-41666670,--250000,20833330,0,--250000,20833330,0,-125000,10416670,-輸出文-1.92E--1.89E-1.32E--6.07E--6.27E--6.08E--6.64E--6.75E--1.96E--2.51E--1.19E--5.01E--4.21E--1.04E---5.28E---9.52E-1.54E--6.12E--1.20E--3.26E--8.94E-1.31E--8.81E--8.26E-4.14E---6.13E--1.13E--2.66E--5.57E--4.79E-2.39E--3.31E--2.55E--1.44E--1.11E--2.81E--2.25E-1.23E--2.90E--1.38E--2.62E--7.91E--1.10E--3.58E--5.04E-1.49E--3.52E--5.25E--9.48E-3.17E--1.15E----------------------------------------------------------------------------------------------------------------------------變量說E節(jié)點(diǎn)位移度T結(jié)構(gòu)位移度總單元位移度總被固定的節(jié)點(diǎn)位移度序各單元節(jié)點(diǎn)PROGRAMDIMENSIONSK(300,30),EK(12,12),Q(300),MC(55),XY(3,100),XYE(3,4),& REWIND7 FORMAT(10X,'NFNENNMBNDLELSEUM REWIND9 STOP1000DIMENSIONXY(LS,NN),Q(MS),NX(NE,LE),MC