Fortran Coder

查看: 263|回复: 4

[文件读写] 小白码完之后出现断点,按书上码的,求救

[复制链接]

10

帖子

4

主题

0

精华

入门

F 币
51 元
贡献
30 点
发表于 2018-5-13 22:41:00 | 显示全部楼层 |阅读模式
[Fortran] 纯文本查看 复制代码
PROGRAM CBSAP
! STRUCTURAL ANALYSIS PROGRAM FOR CONTINUOUS BEAM
INTEGER NJ,IL,IR,NNE,I,J,K
REAL GC(20),GX(20),AMF(20,2),AM(20,2),P(21),AK(21,21),XG(20)
CHARACTER*80 TITLE
OPEN(2,FILE='CBSAP.IN')
OPEN(6,FILE='CBSAP.OUT')
READ(2,*)TITLE
WRITE(6,*)TITLE
READ(2,20)NJ,IL,IR
20 FORMAT(3I5)
WRITE(6,30)NJ,IL,IR
30 FORMAT(/3X,'NUMBER OF NODE='I2,5X,'BIND(L)='I1,2X,'BIND(R)='I1,3X,'1=FIXED 0=SIMPLY SUPPORTED')
NNE=NJ-1
READ(2,40)(GC(I),GX(I),(AMF(I,J),J=1,2),I=1,NNE)
40 FORMAT(4F10.3)
WRITE(6,50)(I,GC(I),GX(I),(AMF(I,J),J=1,2),I=1,NNE)
50 FORMAT(/3X,'NO.E',7X,'LENGHT',6X,'STIFFNESS',6X,'MOMENTENT(L)',6X,'MOMENT(R)'/(I5,4F15.4))
DO I=1,NNE
   XG(I)=GX(I)/GC(I)
END DO
P(1)=-AMF(1,1)
P(NJ)=-AMF(NNE,2)
DO J=2,NNE
   P(J)=-AMF(J-1,2)-AMF(J,1)
END DO
DO I=1,NJ
  DO J=1,NJ
     AK(I,J)=0
  END DO
END DO
AK(1,1)=4.0*XG(1)
AK(NJ,NJ)=4.0*XG(NNE)
DO J=2,NNE
   AK(J,J)=4.0*(XG(J-1)+XG(J))
END DO
DO J=2,NJ
   AK(J,J-1)=2.0*XG(J-1)
   AK(J-1,J)=2.0*XG(J-1)
END DO
IF(IL.EQ.1)THEN
   AK(1,1)=1.0
   AK(1,2)=0.0
   AK(2,1)=0.0
   P(1)=0.0
ELSE
END IF
IF(IR.EQ.0)THEN
   AK(NJ,NJ)=1.0
   AK(NNE,NJ)=0.0
   AK(NJ,NNE)=0.0
   P(NJ)=0.0
ELSE
END IF
DO K=1,NNE
   DO I=K+1,NJ
      C=AK(I,K)/AK(K,K)
      P(I)=P(I)-C*P(K)
      DO J=1,NJ
         AK(I,J)=AK(I,J)-C*AK(K,J)
      END DO
    END DO
END DO
P(NJ)=P(NJ)/AK(NJ,NJ)
DO I=NNE,1,-1
   DO J=I+1,NJ
      P(I)=P(I)-AK(I,J)*P(J)
   END DO
   P(I)=P(I)/AK(I,I)
END DO
WRITE(6,60) (I,P(I),I=1,NJ)
60 FORMAT(/10X,'ANGLE(RADIAN) OF ROTATION'//(5(I3,F10.4)))
DO I=1,NNE
   AM(I,1)=4.0*XG(I)*P(I)+2.0*XG(I)*P(I+1)+AMF(I,1)
   AM(I,2)=2.0*XG(I)*P(I)+4.0*XG(I)*P(I+1)+AMF(I,2)
END DO
WRITE(6,70)(I,(AM(I,J),J=1,2),I=1,NNE)
70 FORMAT(/17X,'MOMENTS OF BEAM ENDS'//3X,'NO.BEAM',6X,'1(LEFT)',5X,'2(RIGHT)'/(I8,2F16.4))
END PROGRAM CBSAP          

U%7_4`$H198QH107T]@_1CK.png BA4RE%~J5ENB1A1{`GXQ[52.png
回复

使用道具 举报

27

帖子

0

主题

0

精华

熟手

F 币
212 元
贡献
101 点
发表于 2018-5-13 22:59:31 | 显示全部楼层
请上传文件CBSAP.IN;无此文件无法重现您所遇到的问题。

10

帖子

4

主题

0

精华

入门

F 币
51 元
贡献
30 点
 楼主| 发表于 2018-5-13 23:11:39 | 显示全部楼层
吉大渣渣-固体 发表于 2018-5-13 22:59
请上传文件CBSAP.IN;无此文件无法重现您所遇到的问题。

没这个文件,我看书上也没有,书上是程序运行直接输入的,不知道理解的偏差大不?

27

帖子

0

主题

0

精华

熟手

F 币
212 元
贡献
101 点
发表于 2018-5-13 23:58:53 | 显示全部楼层
本帖最后由 吉大渣渣-固体 于 2018-5-14 00:00 编辑

那您需要自己制作一个文件读这个文件,或者在程序内使用语句赋值,或者在屏幕输入;总之这里的变量TITLE,NJ,IL,IR、数组GC,GX,AMF需要你初始化。

10

帖子

4

主题

0

精华

入门

F 币
51 元
贡献
30 点
 楼主| 发表于 2018-5-14 09:37:53 | 显示全部楼层
吉大渣渣-固体 发表于 2018-5-13 23:58
那您需要自己制作一个文件读这个文件,或者在程序内使用语句赋值,或者在屏幕输入;总之这里的变量TITLE,NJ ...

好的,谢谢你,我试试
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|群聊|Fcode

GMT+8, 2018-12-19 11:06

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表