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

************************************************************
Q1 created by VDI menu, Version 2007, Date 16/01/07
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
photon use
p;;;;

set prop off
con disx z 1 fi;.001
pause
con disy z 1 fi;.001
pause
con u1/t z 1 fi;.001
vec z 1 co 1
pause
con v1/t z 1 fi;.001
vec z 1 co 1
pause
enduse
DISPLAY
This file presents data-input sets which correspond to 2D (xy)
situations capable of analytical representation in terms of
polynomial stress functions, as described by Timoshenko and
Goodier, 'Theory of Elasticity', Chapter 3, Article 18.

PLate with linear variation of shear stress and
correspondingly non-uniform direct stresses
Analytical solution:
1.   StrY = B3 * Y ; StrXY = - B3 * X; StrX = StrZ = 0

U = - P/E * B3 * X * Y;
V = B3/(2*E)* [ Y^2 - (2+P)*X^2]

2. Boundary Condition for U:
LEFT  : U = 0;
RIGHT : free
UP    : StrXY = - B3*X
DOWN  : StrXY =   B3*X

3. Boundary Condition for V:
LEFT  : V(Y=LY/2) = 0;
RIGHT : StrXY = - B3*LX
UP    : StrY  =   B3*LY/2
DOWN  : StrY  =   B3*LY/2

ENDDIS
************************************************************
IRUNN   =         1 ;LIBREF =       600
************************************************************
Group 1. Run Title
TEXT(2D Plate with polynomial stress function)
************************************************************
Group 2. Transience
************************************************************
Groups 3, 4, 5  Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,10,10,1,1.000000E-05)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS  =    T
* Non-default variable names
NAME(139) =V1/T ; NAME(140) =U1/T
NAME(141) =V1T  ; NAME(142) =U1T
NAME(143) =EPSZ ; NAME(144) =EPSX
NAME(145) =EPSY ; NAME(146) =DRH1
NAME(147) =ENUL ; NAME(148) =PRPS
NAME(149) =DISY ; NAME(150) =DISX
* Solved variables list
SOLVE(DISY,DISX)
* Stored variables list
STORE(PRPS,ENUL,DRH1,EPSY,EPSX,EPSZ,U1T ,V1T )
STORE(U1/T,V1/T)
SOLUTN(DISY,Y,Y,Y,N,N,Y)
SOLUTN(DISX,Y,Y,Y,N,N,Y)

************************************************************
Echo InForm settings for Group  7
inform7begin
REAL(FY,LX,LY,POISSON,YOUNG)
FY= 40.0e6  ! H/m^2 = 40 N/mm^2
LX=120.e-3
LY=20.e-3
REAL(B3)
B3 = FY/(LY/2)

YOUNG   = 1/0.5E-11   ! Young's modulus
POISSON=0.3           ! Poisson's ratio

char(formU,formV)
REAL(DLY)
DLY = LY/2
real(CO1,CO2,CO3)
CO1 = -POISSON/YOUNG*B3
CO2 = 0.5/YOUNG*B3
CO3 = (2+POISSON)
formU=:CO1:*XG*(YG-:DLY:)
formV=:CO2:*((YG-:DLY:)^2-:CO3:*XG^2)

(STORED VAR U1T IS :FORMU:  with swps)
(STORED VAR V1T IS :FORMV:  with swps)

(STORED VAR U1/T IS DISX/(U1T+1.e-20))
(STORED VAR V1/T IS DISY/(V1T+1.e-20))
inform7end
************************************************************
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(V1/T) =  1.001000E-10 ;FIINIT(U1/T) =  1.001000E-10
FIINIT(V1T ) =  1.001000E-10 ;FIINIT(U1T ) =  1.001000E-10
FIINIT(EPSZ) =  1.001000E-10 ;FIINIT(EPSX) =  1.001000E-10
FIINIT(EPSY) =  1.001000E-10 ;FIINIT(DRH1) =  1.001000E-10
FIINIT(ENUL) =  1.001000E-10 ;FIINIT(PRPS) =  1.600000E+02
FIINIT(DISY) =  0.000000E+00 ;FIINIT(DISX) =  0.000000E+00
No PATCHes used for this Group

************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources

PATCH (UPV     ,NORTH ,2,0,0,0,0,0,1,1)

PATCH (DOWNV   ,SOUTH ,3,0,0,0,0,0,1,1)

PATCH (RIGHTV  ,EAST  ,4,0,0,0,0,0,1,1)

PATCH (UPU     ,NORTH ,2,0,0,0,0,0,1,1)

PATCH (DOWNU   ,SOUTH ,3,0,0,0,0,0,1,1)

EGWF    =    T
************************************************************
Echo InForm settings for Group 13
inform13begin
(SOURCE of DISX at LEFTU is COVAL(1.e4,:formU:))
(SOURCE of DISY at LEFTU is COVAL(1.e4,:formV:))
(SOURCE of DISY at UPV is COVAL(FIXFLU,:FY:))
(SOURCE of DISY at DOWNV is COVAL(FIXFLU,:FY:))
(SOURCE of DISY at RIGHTV is COVAL(FIXFLU,-:B3:*:LX:))
(SOURCE of DISX at UPU is COVAL(FIXFLU,:-B3:*XG))
(SOURCE of DISX at DOWNU is COVAL(FIXFLU,:B3:*XG))
inform13end
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP  =       500
RESFAC  = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD  =    T  ;USEGRX =    T
STRA    =    T
ASAP    =    T
PARSOL  =    F
ISG21   =       500
ISG50   =         1
ISG52   =         3
CSG10   ='q1'
SPEDAT(SET,BOUNDARY,ZCONST,R,0.00000E+00)
************************************************************
Group 20. Preliminary Printout
ECHO    =    T
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON   =         8 ;IYMON  =         2 ;IZMON  =         1
NPRMON  =    100000
NPRMNT  =         1
TSTSWP  =        -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT  =    100000
NXPRIN  =         1
NYPRIN  =         1
ISWPRF  =         1 ;ISWPRL =    100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts

GVIEW(P,-4.999979E-03,0.000000E+00,9.999875E-01)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)

> DOM,    SIZE,        1.200000E-01, 2.000000E-02, 1.000000E-03
> DOM,    MONIT,       9.000000E-02, 3.000000E-03, 5.000000E-04
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> GRID,   RSET_X_1,      1, 1.000000E+00
> GRID,   RSET_X_2,      9, 1.000000E+00
> GRID,   RSET_Y_1,     10, 1.000000E+00
> GRID,   RSET_Z_1,      1, 1.000000E+00

> OBJ,    NAME,        LEFTU
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.200000E-02, 2.000000E-02, 1.000000E-03
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED

> OBJ,    NAME,        UPV
> OBJ,    POSITION,    0.000000E+00, 2.000000E-02, 0.000000E+00
> OBJ,    SIZE,        1.200000E-01, 0.000000E+00, 1.000000E-03
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED

> OBJ,    NAME,        DOWNV
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.200000E-01, 0.000000E+00, 1.000000E-03
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED

> OBJ,    NAME,        RIGHTV
> OBJ,    POSITION,    1.200000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 2.000000E-02, 1.000000E-03
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
************************************************************
STOP
```