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 V1 contours
  msg
  cont V1 z 1 fil;.0001
  pause
 
  cl
  msg
  msg U1 contours
  msg
  cont U1 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
  msg
  msg STRZ contours
  msg
  cont STRZ z 1 fil;.0001
  pause
 
  cl
  msg
  msg STRX contours
  msg
  cont STRX z 1 fil;.0001
  pause
 
  cl
  msg
  msg STRY contours
  msg
  cont STRY z 1 fil;.0001
  pause
 
  cl
  msg
  msg STXY contours
  msg
  cont STXY z 1 fil;.0001
  pause
 
  cl
  msg Displacements fields
  gr ou z 1
  cont DELT z 1 fil;.0001
  vec z 1 col 0
  pause
 
  cl
  AUTOPLOT
  file
  phi 5
  EDZZ\S201_1.txt 2
  EDZZ\S201_2.txt 2
  EDZZ\S201_3.txt 2
 
  data 1
  STRX X 11 Y 1 28
  col1 1
  data 2
  y SXX
  colA 2
  scale
  msg STRX(Y) for X=45 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2
  pause
 
  clear
  data 1
  STRY X 11 Y 1 28
  col1 1
  data 2
  y SYY
  colA 2
  scale
  msg STRY(Y) for X=45 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2
  pause
 
  clear
  data 1
  STXY X 11 Y 1 28
  col1 1
  data 2
  y TXY
  colA 2
  scale
  msg STXY(Y) for X=45 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2
  pause
 
 
  clear
  data 1
  U1 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
  V1 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
  EPSZ X 11 Y 1 28
  col1 1
  data 3
  y EZZ
  colA 2
  scale
  msg EpsZ(Y) for X=45 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2
  pause
 
  clear
  data 1
  U1 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
  V1 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 =         201
 ************************************************************
  Group 1. Run Title
 TEXT(2D xy Plate with square hole; s201      )
 ************************************************************
  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) =DELT ; NAME(141) =DELY
 NAME(142) =DELX ; NAME(143) =STXY
 NAME(144) =EPSZ ; NAME(145) =EPSX
 NAME(146) =EPSY ; NAME(147) =STRZ
 NAME(148) =STRY ; NAME(149) =STRX
 NAME(150) =PRPS
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  )
    * Stored variables list
 STORE(PRPS,STRX,STRY,STRZ,EPSY,EPSX,EPSZ,STXY)
 STORE(DELX,DELY,DELT)
 
 ************************************************************
  Echo InForm settings for Group  7
  INFORM7BEGIN
 
(STORED VAR DelX IS 0.5*(U1+U1[-1]) with imat>100)
(STORED VAR DelY IS 0.5*(V1+V1[&-1]) with imat>100)
(STORED VAR Delt IS SQRT(DELX^2+DELY^2) with imat>100)
  INFORM7END
 ************************************************************
  Group 8. Terms & Devices
 ************************************************************
  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(P1  ) =  0.000000E+00 ;FIINIT(U1  ) =  0.000000E+00
 FIINIT(V1  ) =  0.000000E+00 ;FIINIT(DELT) =  0.000000E+00
 FIINIT(DELY) =  1.001000E-10 ;FIINIT(DELX) =  1.001000E-10
 FIINIT(STXY) =  1.001000E-10 ;FIINIT(EPSZ) =  1.001000E-10
 FIINIT(EPSX) =  1.001000E-10 ;FIINIT(EPSY) =  1.001000E-10
 FIINIT(STRZ) =  1.001000E-10 ;FIINIT(STRY) =  1.001000E-10
 FIINIT(STRX) =  1.001000E-10 ;FIINIT(PRPS) =  1.600000E+02
   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      ,V1  , 1.000000E+00, 0.000000E+00)
 
 PATCH (AXESZZ  ,WWALL ,4,0,0,0,0,0,1,1)
 COVAL (AXESZZ  ,U1  , 1.000000E+00, 0.000000E+00)
 
 PATCH (FORZ01  ,NORTH ,3,0,0,0,0,0,1,1)
 COVAL (FORZ01  ,V1  , FIXFLU      ,-4.000000E+07)
 
 EGWF    =    T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =      2000
 RESREF(P1  ) = -1.000000E+00 ;RESREF(U1  ) = -1.000000E+00
 RESREF(V1  ) = -1.000000E+00
 SELREF  =    F
 ************************************************************
  Group 16. Terminate Iterations
 LITER (U1  ) =       50 ;LITER (V1  ) =       50
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 1.000000E+00)
 ************************************************************
  Group 18. Limits
 VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06
 VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 STRA    =    T
 ASAP    =    T
 PARSOL  =    T
 ISG21   =      2000
 ISG50   =         1
 ISG52   =         1
 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,0.000000E+00,1.000000E+00,0.000000E+00)
 
> DOM,    SIZE,        9.000000E-02, 1.200000E-01, 1.000000E-03
> DOM,    MONIT,       7.928572E-02, 6.428572E-03, 5.000000E-04
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    INCREMENT,   1.000000E-02, 1.000000E-02, 1.000000E-02
> GRID,   RSET_X_1,      7, 1.000000E+00
> GRID,   RSET_X_2,     14, 1.000000E+00
> GRID,   RSET_Y_1,      1, 1.000000E+00
> GRID,   RSET_Y_2,     20, 1.000000E+00
> GRID,   RSET_Y_3,      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,    cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        FORZ01
> OBJ,    POSITION,    0.000000E+00, 4.285715E-03, 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,    cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 ************************************************************
STOP