beamgn平面梁大挠度分析程序

《beamgn平面梁大挠度分析程序》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、BEAMGN平面梁大挠度分析程序BEAMGN——平面梁大挠度分析程序****************************************************THENONLINEARFINITEELEMENTANALYSIS****PROGRAMFORBEAMT.L.ANDU.L.METHOD**************************************************PROGRAMMAIN**************************************************

2、**变量说明*FILE1:输入文件名*FILE2:输出文件名*FILE3:图形结果数据文件名*NTU:计算方法指示数，1为T.L.，2为U.L.*LOAD：载荷因子数*EF:拉伸刚度*EI:抗弯刚度*NE:单元总数*IME:增量步总数*ACC:计算精度*MM:迭代信息数，1为等刚度，2为变刚度*NRM:广义外力作用个数*NU:初始位移约束个数*XX1,YY1,XX2,YY2:梁两端点坐标*PNM:广义力信息，PNM(I,1),PNM(I,2),PNM(I,3)分别为作用点号，方向和数值*UXY:结点位移约束信息，UXY(I,

3、1),UXY(I,2),UXY(I,3)为作用点号，方向和数值*M:结构自由度数*NE:单元总数*BO：单元两结点连线与坐标方向夹角*PP:总结点载荷数组*U:总位移数组*P:结点增量载荷数组*************************************************IMPLICITREAL*8(A-H,O-Z)DIMENSIONWK(63,63),U(63),DU(63),RS(63,1),PP(63,1),RN(20),&XY(21,2),DLL(20),BO(20),XYO(21,2),T

4、XY(20,2),TL1(20),&PNM(10,4),UXY(10,3),FP(63,1),P(63,1),UI(63),TL2(20),&PP1(63,1),UU(63)***INPUTDATA*CHARACTERFILE1*12,FILE2*12WRITE(*,*)'PLEASEINPUTTHEINPUTFILENAME'READ(*,'(A)')FILE1OPEN(1,FILE=FILE1,STATUS='OLD')OUT=2IF(OUT.EQ.2)TH

5、ENWRITE(*,*)'PLEASEINPUTTHEOUTPUTFILENAME'READ(*,'(A)')FILE2OPEN(2,FILE=FILE2,STATUS='NEW')WRITE(2,*)'------THERESULTSOFCALCULATION------'ENDIFREAD(1,*)NTU,LOADWRITE(*,*)'------THERESULTSOFCALCULATION------'READ(1,*)EF,EIREAD

6、(1,*)NE,IME,ACC,MM,NRM,NUNP=NE+1READ(1,*)XX1,YY1READ(1,*)XX2,YY2XXX=(XX2-XX1)/NEYYY=(YY2-YY1)/NEDO10I=1,NPXY(I,1)=XXX*(I-1)XY(I,2)=YYY*(I-1)10CONTINUEDO11I=1,NRMREAD(1,*)PNM(I,1),PNM(I,2),PNM(I,3)11CONTINUEDO20I=1,NUREAD(1,*)UXY(I,1),UXY(I,2),UXY(I,3)20CONTINUEM=3*

7、NPM1=M+1PI=3.1415926DO50I=1,NE50BO(I)=0.0DO60I=1,MPP(I,1)=0.0U(I)=0.060P(I,1)=0.0***FORMTOTALLOADCOLUMNMATRIX*DO70K=1,NRMII=3*(PNM(K,1)-1)+PNM(K,2)70P(II,1)=PNM(K,3)/IMEDO75I=1,NPDO75J=1,275XYO(I,J)=XY(I,J)NO=0***FORMLOCALCOORDINATE*110CALLFLC(NE,M,RS,UI,DU,RN,XY,B

8、O,DLL)1500NO=NO+1***FORMINCREMENTLOADPP*IF(LOAD.GT.0)THENIF(NTU.EQ.2)THENDO90I=1,M90PP(I,1)=P(I,1)ELSEDO95I=1,M95PP(I,1)=P(I,1)*NOENDIFELSEE=M-1F