C 主程序
C
DIMENSION C(500000),IA(50000),EK(36)
CHARACTER*12 IN,OUT
C IN和OUT为输入文件和输出文件的文件名
WRITE(*,*)' '
WRITE(*,*)' PLEASE INPUT THE INPUT-FILE NAME (A<12)'
WRITE(*,*)' '
READ(*,5) IN
C 输入输入文件的文件名
WRITE(*,*)' '
WRITE(*,*)' PLEASE INPUT THE OUTPUT-FILE NAME (A<12)'
WRITE(*,*)' '
READ(*,5) OUT
C 输入输出文件的文件名
5 FORMAT(A12)
8 OPEN(5,FILE=IN, STATUS='OLD')
OPEN(6,FILE=OUT,STATUS='UNKNOWN')
C 打开对应的输入和输出文件
10 READ(5,*) NG
IF(NG.EQ.0) STOP
write(6,15)
15 FORMAT(//,5X,
& '******************* PLANE PROBLEM *******************',/)
C 输入结构结点数;如果结点数为0则停止运行
READ(5,*)NE,MC,NX,NB,ND,EO,VO,T
WRITE(6,20)NG,NE,MC,NX,NB,ND,EO,VO,T
20 FORMAT(//,10X,'NG=',I4,5X,'NE=',I4,5X,'MC=',I4,5X,'NX=',I4,/,
& 10X,'NB=',I4,5X,'ND=',I4,/,10X,'EO=',E10.4,5X,'VO=',F7.2,
& 5X,'T=',F7.2)
C 按顺序输入结构单元数,问题类型参数,载荷组数,给定位移个数
C 结构刚度阵的半带宽,弹性模量,泊松比和结构厚度
READ(5,*)NWA,NWE,NWK,NWP,NWD
WRITE(6,25)NWA,NWE,NWK,NWP,NWD
25 FORMAT(/10X,'NWA=',I2,3X,'NWE=',I2,3X,'NWK=',I2,3X,
& 'NWP=',I2,3X,'NWD=',I2)
C 按顺序输入各输出控制参数
NT=2*NG
C 确定总刚度矩阵阶数NT
C
C 计算变界数组的下限
C
M1=3*NE+1
M2=M1+2*NB
N1=2*NG+1
N2=N1+NB
N3=7*NE+N2
N4=N3+NT*ND
N5=N4+NT
C 得到各变界数组在一维大数组中的起始元素编号
C
C 检验实型数组C的下限
C
NEND=N5
IF(NEND.LE.500000) GOTO 35
WRITE(*,*)'*** EXCEED THE LIMIT OF ARRAY C(IN THE MIDDLE)!! ***'
