```
TALK=F;RUN( 1, 1)

************************************************************
Q1 created by VDI menu, Version 2009, Date 17/02/09
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
************************************************************
IRUNN = 1 ;LIBREF = 3
************************************************************
Group 1. Run Title
TEXT(USP: 2D xy Plate with square hole; U505)
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5  Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,7,7,1,1.0E-05)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* 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)
USP = T ; USPDBG= F
UTCPLT = T ; USPIMB= F
UAUTO = T ; UGTGRD= F
USPVTK = T ; USPTER= F
USCANX = T ; USCANY= T
USCANZ = T
MXLEV = 0 ;MYLEV = 0 ;MZLEV = 0
DOMAT = -1
CELLST = 10 ; FACEST= 10
MINPRP = -1 ; MAXPRP= 250
UMXLEV = 3 ; UMXCLS= 10000
UNOLAY = 2 ; UIOLEV= 1
UWALEV = 2 ; UGRITR= 2
UMXCRS = 16
USPREL =0.7 ; UMNCEL=1.
UDVREL =1.
************************************************************
Group 9. Properties
* Domain material index is 160 signifying:
* NOTSET
SETPRPS(1,160,Q1)
ENUT =0.
DRH1DP =5.0E-12
DVO1DT =1.0E-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.001E-10 ;FIINIT(DRH1)=1.001E-10
FIINIT(STXY)=1.001E-10 ;FIINIT(EPSZ)=1.001E-10
FIINIT(EPSX)=1.001E-10 ;FIINIT(EPSY)=1.001E-10
FIINIT(STRY)=1.001E-10 ;FIINIT(STRX)=1.001E-10
FIINIT(PRPS)=160. ;FIINIT(DISY)=0.
FIINIT(DISX)=0.
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
No PATCHes used for this Group
SAVE13BEGIN
(SOURCE DISY at UP is COVAL(1,0)!VOLU)
(SOURCE DISX at AXESZZ is COVAL(1,0)!VOLU)
(SOURCE DISY at FORZ01 is COVAL(FIXFLU,-4.0e+7)!VOLU)
SAVE13END

EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 100
RESFAC =1.0E-03
************************************************************
Group 16. Terminate Iterations
LITER(DISY)=50 ;LITER(DISX)=50
************************************************************
Group 17. Relaxation
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
STRA = T
ISG21 = 100
ISG50 = 1
ISG52 = 3
CSG10   ='q1'
SPEDAT(SET,BOUNDARY,ZCONST,R,1.0E+20)
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 7 ;IYMON = 1 ;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.023095,0.148834,0.988593)
GVIEW(UP,-0.138978,0.979733,-0.144253)

> DOM,    SIZE,        9.000000E-02, 1.200000E-01, 1.000000E-02
> 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-04, 1.000000E-03, 1.000000E-06
> GRID,   RSET_X_1,      7, 1.000000E+00
> GRID,   RSET_Y_1,      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-02
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    198,Solid with smooth-wall friction

> OBJ,    NAME,        UP
> OBJ,    POSITION,    0.000000E+00, 1.200000E-01, 0.000000E+00
> OBJ,    SIZE,        9.000000E-02, 0.000000E+00, 1.000000E-02
> OBJ,    GEOMETRY,    cube11
> OBJ,    TYPE,        PLATE

> OBJ,    NAME,        FORZ01
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        9.000000E-02, 0.000000E+00, 1.000000E-02
> OBJ,    GEOMETRY,    cube11
> OBJ,    TYPE,        PLATE

> OBJ,    NAME,        AXESZZ
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 1.200000E-01, 1.000000E-02
> OBJ,    GEOMETRY,    cube11
> OBJ,    TYPE,        PLATE
************************************************************
STOP
```