Scribd
Upload
15 views2 pages

GRID

The document describes a program for generating a finite element grid. It defines various dimensions and reads in title, control card and coordinate data. It then loops through regions, generating element nodal coordinates and connectivity data, and outputs this information.

Uploaded by

vaskomane
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
15 views2 pages

GRID

The document describes a program for generating a finite element grid. It defines various dimensions and reads in title, control card and coordinate data. It then loops through regions, generating element nodal coordinates and connectivity data, and outputs this information.

Uploaded by

vaskomane
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
You are on page 1/ 2

PROGRAM GRID

DIMENSION TITLE(10),XP(100),YP(100),XRG(9),YRG(9)010),NON(8)
DIMENSION NN(21,21),YC(21,21),XC(21,21),NNRB(20,4,21),JT(20,4)
DIMENSION LB(3),NE(400),XE(400),YE(400),NR(4),ICOMP(4,4)
REAL N
DATA ICOMP/-1,1,1,-1,1,-1,-1,1,1,-1,-1,1,-1,1,1,-1/
DATA IN/60/,IO/61/,IP/62/,NBW/0/,NB/0/,NEL/8/
READ(IN,17) TITLE
17 FORMAT(A8)
READ(IN,1) INRG,INBP,IPCH
1 FORNAF(3I3)
READ(IN,3)(XP(1),I=1,INBP)
READ(IN,3)(YF(1),I=1,INBP)
3 FORMAT(8F10,5)
DO 2 I = 1,INRG
2 READ(IN,8) NRG,(JT(NRG,J),J = 1,4)
FORMAT(5I3)
WRITE(I0,36) TITLE
36 FORMAT(1H1////lX,10A8//1X,18HGLOBAL COOROINATES //1X,30HNUMBER
1 X COORD Y COORO )
WRITE(10,30) (I,XP(I),YP(1),I = 1,INBP)
30 FORMAT(2X,I3,7X,F7.2,5X,F7.21)
WRITE(10,21)
21 FORMAT(//1X,17HCONNECTIVITY DATA/1X,41HREGION SIDE 1
1 3 4 )
D0 26 I=1,INRG
26 WRITE(10,22) I,(JT(I,J)),J = 1,4)
22 FORMAT(2X,I3,14X,4(I2,5X))
C*******************
C LOOP ON THE REGIONS TO GENERATE THE ELEMENTS
C*******************
DO 16 KK = 1, INRG
READ(IN,4) NRG,NROWS,NCOL,NON
4 FORMAT(11I3)
WRITE(I0,18) NRG,NROWS,NCOL,(NON(I),I = 1,8)
10 FORMAT(1H1///1X,12H*** REGION ,I2,6H), ****//10X,I2,5H ROWS, 10X, I2
1 ,8H COLUMNS//10X,21HBOUNDARY NODE NUMBERS, 10X,8I5)
C
C GENERATION OF THE ELEMENT NODAL COOROINATES
C
DO 5 I=1,0
II=NON(I)
XRG(I)=XP(II)
5 YRG(I) = - YP(II)
XRG(9) = XRG(1)
YRG(9) = YRG(1)
TR=NROWS-1
OETA = 2./TR
I =NC01 - 1
DSI = 2./TP
DO 12 I =1,NROWS
TR=I-1
ETA=10.- TR*OETA
DO 12 J=1,NCOL
TR = J-1
SI=-1.+TR*OSI
N(1)=-0.25*(1.-SI)*(1.-ETA)*(SI+ETA+1.)
N(2)=0.50*(1.-SI"2)
N(3),10.25*(1.+SI10(1
N(4)=0.50*(1.+SI))*(1
N(5)=0.25*(1.+SI)*(1+
N(6)=0.50*(1.,,,S1**2) W1m0,25*(1,-SI)*(1 N(0)=0,50441,�SI)*(1

70 XC(I,J)=0,0
YC(I,J)3,0,0
0012X=1,6
XC(I,J)=XC(I,J),XRG(K)*N(A)
12 YCII,J/sYC(1,J)0YRG(K)*N(K)
75
C GENERATION OF THE REGION NODE NUMBERS
C
KN1=1
KS1 1
50 KN2=NROWS
KS2= NCOL
00501:1,4
NRI-JT(NRG,I)
ITORT,E0,0,0RORT,GT,NRG) GOT050
$5 0056,,1=1,4
56 IF(JT(NRI,J),E0,NRG) NRTS=,0
K=NCOL KiNROWS
IF(I,E0,2,0R,I,E0,4)
JL-1
90 JK-ICOMP(I,NRIS1
IF(JK,E0,-1) JOIK
0044J:1,K
5 10 15 20 25 30 35 40 45 50 55
C
C INPUT ANO OUTPUT OF TITLE, CONTROL CARO, GLOBAL COORDINATES ANO
C CONNECTIVITY DATA
C ***
1,-ETA)*(SI+ETA,1,)
'(1,-ETA)
,-ETA)*(SI-ETA-1,)
,-ETA*02)
,,,ETA1*(SI*ETA-1,)
�(1,,,ETA)
,*ETA)*(ETA-S1-1,)
,-ETA**21

You might also like