主题:fortran解平面梁问题
fortran读取数据到单元5就结束了,forrt1:severe(24):end-of-file during read,unit5,
我怀疑是读数据的问题,开头部分数据不能全读入,到单元5就终止了。
我输入,输出文件都有,但是输出文件只有部分输出
我把程序上传了,大神有时间给看看,万分感激
http://image.keyan.cc/data/bcs/2015/0630/w145h3869417
http://image.keyan.cc/data/bcs/2015/0630/w146h3869417_1435648114_346.png
http://image.keyan.cc/data/bcs/2015/0630/w166h3869417_1435648129_349.png
http://image.keyan.cc/data/bcs/2015/0630/w110h3869417_1435648984_742.png
SUBROUTINE CST(LOC,IFIX,AREA,FINT,CX,CY,F,FR,FE) DIMENSION LOC(NE,2),IFIX(NFIX),AREA(NE),FINT(NE),CX(NN),CY(NN), 1 F(ND),FR(ND),FE(NE,6), 2 P1(3,6),P2(3,6),GK(12,12),GK1(12,12),AL(3),TT(3,6,6),BKL(3,6,6), 3 T(6,6),CSS(3),SNN(3),EK(6,6),EKL(6,6),EKK(3,6,6),XX(6),BA(6,6) COMMON NN,NE,ND,NFIX,E C------------------------------------------------------- DO 5 I=1,NE I1=LOC(I,1) I2=LOC(I,2) X12=CX(I2)-CX(I1) Y12=CY(I2)-CY(I1) AL(I)=SQRT(X12* *2+Y12* *2) CSS(I)=X12/AL(I) SNN(I)=Y12/AL(I) 5 CONTINUE C FOR CONCENTRAL FORCE,SET TO A NODE READ (5,*) NP,NVD IF (NP.EQ.0) GO TO 15 DO 10 II=1,NP READ (5,*) (I,F(3*I-2),F(3*I-1),F(3*I),I=2,4) 10 CONTINUE 15 CONTINUE C FOR VERTICAL DISTRIBUTED LOADING OF AN ELEMENT IF(NVD.EQ.0) GO TO 25 DO 20 II=1,NVD READ(5,*)IE,Q I1=LOC(IE,1) I2=LOC(IE,2) P1(IE,3*I1-1)=Q*AL(IE)/2. P1(IE,3*I1)=Q*AL(IE)* *2/12. P1(IE,3*I2-1)=Q*AL(IE)/2. P1(IE,3*I2)=-Q*AL(IE)* *2/12. 20 CONTINUE
http://image.keyan.cc/data/bcs/2015/0630/w149h3869417_1435649483_667.png