교수님께서 평면 탄성 유한요소 해석을 위해 fortran 77로 짜여진 파일을 주셨는데요, 생소한 언어인데다 인풋 파일을 어떤식으로 만들어야 하는지
감도 안잡혀서 이렇게 질문 하게 됐습니다. 아무리 찾아보고 공부하려 해봤지만 어째서 [10 FORMAT(12I5)인지] 감도 안잡히네요 제가 이해하기론
12개의 변수를 5자리의 공간을 주고 받는 다는 것인데, 실제로 read해야할 변수는 NODT,NELT,KOX,KOY,NF 5개가 아닌가요? ㅠ
data input 에 관한 코드는 아래와 같습니다.
C
C ****************************************************************
C **** STRESS ANALYSIS OF IN-PLANE ELASTICITY PROBLEM ****
C ****************************************************************
C
DIMENSION JIJOM(100,3),X(50),Y(50),T(100),F(100),SM(6.6)
1 TSM(100,100),INDEX(100),DIS(100),STR(3,6),REAC(100),NOKX(20)
2 NOKY(20),WD(6),WS(6)
C
OPEN(5,FILE='STRESS.DAT')
OPEN(6,FILE='STRESS.OUT')
C
C ******** DATA INPUT ********
C
READ(5,10) NSTRES
IF(NSTRES.EQ.0) GO TO 5
WRITE(6,6)
6 FORMAT(1H0,32H********* PLANE STRESS ********* ) GO TO 7
5 WRITE(6,8)
8 FORMAT(1H0,32H********* PLANE STRAIN ********* ) 7 READ(5,10) NODT,NELT,KOX,KOY,NF
10 FORMAT(12I5)
C
DO 15 I=1,NELT
15 READ(5,11) (JIJOM(I,J),J=1,3),T(I)
11 FORMAT(3I5,5X,F10.0)
DO 20 I=1,NODT
20 READ(5,25) X(I),Y(I)
25 FORMAT(2F10.0)
READ(5,25) E,PO
READ(5,10) (NOKX(I),I=1,KOX)
READ(5,10) (NOKY(I),I=1,KOY)
NT=NODT*2
DO 30 I=1,NT
30 F(I)=0
DO 35 I=1,NF
READ(5,40) NO,FX,FY
40
FORMAT(I5,5X,2F10.0)
F(2*NO-1)=FX
35 F(2*NO)=FY
C
나름대로 작성해본 input file은 다음과 같습니다.
******************************************************
*** STRESS ANALYSIS OF IN-PLANE ELASTICITY PROBLEM ***
******************************************************
1
12 12 3 6 3
2 1 3 45.
1 4 3 45.
2.1 0.25
0 96
60 96
0 0
0 0
60 0
60 96
2 0. 1500
하지만 역시나 컴파일이 되질 않네요...
어떠한 형식으로 STRESS.DAT 을 작성해야 할지 도움좀 주시면 정말 감사하겠습니다!.. ㅠ