TALK=F;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2007, Date 26/09/07
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
DISPLAY
A rectangular plate with a centrally-placed square hole is
extended by uniform equal forces on its upper and lower
surfaces, which are allowed to distort.
By reason of symmetry, only one quarter of the plate is analysed.
Displacements in the third (z) dimension are prevented. The
problem is therefore one of plane strain.
No analytical solution exists; but the computed displacements,
strains and stresses are in close agreement with those computed
by the finite-element code ELCUT.
ENDDIS
PHOTON USE
p;;;;
set prop off
msg
msg DISX contours
msg
cont DISX z 1 fil;.0001
SET VEC COMP - - -
SET VEC COMP DISX DISY -
VECT IZ 1 SH DASH 0
pause
cl
msg
msg DISY contours
msg
cont DISY z 1 fil;.0001
pause
cl
msg
msg EpsX contours
msg
cont EPSX z 1 fil;.0001
pause
cl
msg
msg EpsY contours
msg
cont EPSY z 1 fil;.0001
pause
cl
AUTOPLOT
file
phi 5
EDZZ\S201_1.txt 2
EDZZ\S201_2.txt 2
EDZZ\S201_3.txt 2
clear
data 1
DISX X 11 Y 1 28
col1 1
data 3
y DX
colA 2
scale
msg U1(Y) for X=45 mm (WHITE-Phoenics, YELLOW - Elcut) ...
plot 1 2
pause
clear
data 1
DISY X 11 Y 1 28
col1 1
data 3
y DY
colA 2
scale
msg V1(Y) for X=45 mm (WHITE-Phoenics, YELLOW - Elcut) ...
plot 1 2
pause
clear
data 1
EPSX X 11 Y 1 28
col1 1
data 3
y EXX
colA 2
scale
msg EpsX(Y) for X=45 mm (WHITE-Phoenics, YELLOW - Elcut) ...
plot 1 2
pause
clear
data 1
EPSY X 11 Y 1 28
col1 1
data 3
y EYY
colA 2
scale
msg EpsY(Y) for X=45 mm (WHITE-Phoenics, YELLOW - Elcut) ...
plot 1 2
pause
clear
data 1
DISX Y 14 X 1 21
col1 1
data 4
x DX
colA 2
scale
msg U1(X) for Y=59 mm (WHITE-Phoenics, YELLOW - Elcut) ...
plot 1 2
pause
clear
data 1
DISY Y 14 X 1 21
col1 1
data 4
x DY
colA 2
scale
msg V1(X) for Y=59 mm (WHITE-Phoenics, YELLOW - Elcut) ...
plot 1 2
pause
ENDUSE
************************************************************
IRUNN = 1 ;LIBREF = 3
************************************************************
Group 1. Run Title
TEXT(2D xy Plate with square hole; S601 )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,21,28,1,1.000000E-05)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(140) =VISL ; NAME(141) =DRH1
NAME(142) =STXY ; NAME(143) =EPSZ
NAME(144) =EPSX ; NAME(145) =EPSY
NAME(146) =STRY ; NAME(147) =STRX
NAME(148) =PRPS ; NAME(149) =DISY
NAME(150) =DISX
* Solved variables list
SOLVE(DISY,DISX)
* Stored variables list
STORE(PRPS,STRX,STRY,EPSY,EPSX,EPSZ,STXY,DRH1)
STORE(VISL)
* Additional solver options
SOLUTN(DISY,Y,Y,Y,N,N,Y)
SOLUTN(DISX,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
TERMS (DISY,N,N,Y,N,Y,N)
TERMS (DISX,N,N,Y,N,Y,N)
************************************************************
Group 9. Properties
* Domain material index is 160 signifying:
* NOTSET
SETPRPS(1,160,Q1)
ENUT = 0.000000E+00
DRH1DP = 5.000000E-12
DVO1DT = 1.000000E-05
* List of user-defined materials to be read by EARTH
MATFLG=T;IMAT=1
* Name
*Ind. Dens. Viscos. Spec.heat Conduct. Expans. Compr.
*
160 7800.0 0.3 473.0 43.0 1.0E-5 0.5E-11
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(VISL) = 1.001000E-10 ;FIINIT(DRH1) = 1.001000E-10
FIINIT(STXY) = 1.001000E-10 ;FIINIT(EPSZ) = 1.001000E-10
FIINIT(EPSX) = 1.001000E-10 ;FIINIT(EPSY) = 1.001000E-10
FIINIT(STRY) = 1.001000E-10 ;FIINIT(STRX) = 1.001000E-10
FIINIT(PRPS) = 1.600000E+02 ;FIINIT(DISY) = 0.000000E+00
FIINIT(DISX) = 0.000000E+00
No PATCHes used for this Group
INIADD = F
************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources
PATCH (UP ,NWALL ,2,0,0,0,0,0,1,1)
COVAL (UP ,DISY, 1.000000E+00, 0.000000E+00)
PATCH (FORZ01 ,SOUTH ,3,0,0,0,0,0,1,1)
COVAL (FORZ01 ,DISY, FIXFLU ,-4.000000E+07)
PATCH (AXESZZ ,WWALL ,4,0,0,0,0,0,1,1)
COVAL (AXESZZ ,DISX, 1.000000E+00, 0.000000E+00)
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 100
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
LITER (DISY) = 50 ;LITER (DISX) = 50
************************************************************
Group 17. Relaxation
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
STRA = T
ASAP = T
PARSOL = T
ISG21 = 100
ISG50 = 1
ISG52 = 3
CSG10 ='q1'
SPEDAT(SET,BOUNDARY,ZCONST,R,1.00000E+20)
************************************************************
Group 20. Preliminary Printout
ECHO = T
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 19 ;IYMON = 2 ;IZMON = 1
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+00)
GVIEW(UP,4.999979E-03,9.999875E-01,0.000000E+00)
> DOM, SIZE, 9.000000E-02, 1.200000E-01, 1.000000E-03
> DOM, MONIT, 7.928570E-02, 6.428570E-03, 5.000000E-04
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> GRID, RSET_X_1, 7, 1.000000E+00
> GRID, RSET_X_2, 14, 1.000000E+00
> GRID, RSET_Y_1, 21, 1.000000E+00
> GRID, RSET_Y_2, 7, 1.000000E+00
> GRID, RSET_Z_1, 1, 1.000000E+00
> OBJ, NAME, FLOW1
> OBJ, POSITION, 0.000000E+00, 9.000000E-02, 0.000000E+00
> OBJ, SIZE, 3.000000E-02, 3.000000E-02, 1.000000E-03
> OBJ, GEOMETRY, cubet
> OBJ, ROTATION24, 1
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 0, Air at 20 deg C, 1 atm, treated as i
> OBJ, SCAL_FIXF, Per-unit-volume
> OBJ, NAME, UP
> OBJ, POSITION, 3.000000E-02, 1.200000E-01, 0.000000E+00
> OBJ, SIZE, 6.000000E-02, 0.000000E+00, 1.000000E-03
> OBJ, GEOMETRY, default
> OBJ, ROTATION24, 1
> OBJ, TYPE, USER_DEFINED
> OBJ, NAME, FORZ01
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 9.000000E-02, 0.000000E+00, 1.000000E-03
> OBJ, GEOMETRY, default
> OBJ, ROTATION24, 1
> OBJ, TYPE, USER_DEFINED
> OBJ, NAME, AXESZZ
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 9.000000E-02, 1.000000E-03
> OBJ, GEOMETRY, default
> OBJ, ROTATION24, 1
> OBJ, TYPE, USER_DEFINED
************************************************************
STOP