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 CIRCULAR hole is
  extended by uniform equal forces on its left and right
  surfaces, which are allowed to distort.
  By reason of symmetry, only one quarter of the plate is analysed.
 
  Comparison - a analytical solution (Timoshenko and Goodier).
  See also: I. Demirdzic, S. Musafija
  "Finite Volume method for stress analysis in Complex Domains"
  Int. J. for numerical methods in engineering",
  v. 37, 3751-3756 (1994).)
 
  ENDDIS
 ************************************************************
 IRUNN   =         1 ;LIBREF =       202
 ************************************************************
  Group 1. Run Title
 TEXT(USP:2D xy Plate with CIRCULAR hole   )
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  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(133) =XY/T ; NAME(134) =SXYT
 NAME(135) =SY/T ; NAME(136) =SYTH
 NAME(137) =SX/T ; NAME(138) =SXTH
 NAME(139) =TET7 ; NAME(140) =R7
 NAME(141) =EPSZ ; NAME(142) =EPSX
 NAME(143) =EPSY ; NAME(144) =STXY
 NAME(145) =STRZ ; NAME(146) =STRY
 NAME(147) =STRX ; NAME(148) =PRPS
 NAME(149) =DISX ; NAME(150) =DISY
    * Solved variables list
 SOLVE(DISX,DISY)
    * Stored variables list
 STORE(PRPS,STRX,STRY,STRZ,STXY,EPSY,EPSX,EPSZ)
 STORE(R7  ,TET7,SXTH,SX/T,SYTH,SY/T,SXYT,XY/T)
 STORE(ZZZ7)

 ************************************************************
  Echo save-block settings for Group  7
   SAVE7BEGIN
 
REAL(FX,LX,LY,R0,POISSON,YOUNG)
FX= 1.0e4
LX=2.0
LY=LX
R0=0.5
YOUNG   = 1/1.0E-7   ! Young's modulus
POISSON=0.3           ! Poisson's ratio

char(form1,form2)

(STORED VAR R7 IS SQRT(XG^2+YG^2))
(STORED VAR TET7 IS ATAN(YG/(XG+1.e-10)))

form1=:FX:*(1.0-(:R0:/R7)^2*(1.5*COS(2*TET7)+
form2=COS(4*TET7))+1.50*(:R0:/R7)^4*COS(4*TET7))
(STORED VAR SXTH IS :form1::form2: with imat>100)
(STORED VAR SX/T IS STRX/SXTH with imat>100)

char(frm3,frm4)
frm3=:FX:*(-(:R0:/R7)^2*(0.5*COS(2*TET7)-
frm4=COS(4*TET7))-1.5*(:R0:/R7)^4*COS(4*TET7))
(STORED SYTH IS :frm3::frm4: with imat>100)    
(STORED SY/T IS STRY/SYTH with imat>100)
 

char(ff1,ff2,ff3,Tf1,Rf1)
real(rf2)
rf2=R0**2
Rf1=(:Rf2:/(YG^2+XG^2))
Tf1=ATAN(YG/XG)
ff1=:FX:*(-:Rf1:*(0.5*SIN(2*:Tf1:)+
ff2=SIN(4*:Tf1:))+
ff3=1.5*:Rf1:^2*SIN(4*:Tf1:))
(STORED VAR SXYT IS :ff1::ff2::ff3: with imat>100)
(STORED VAR XY/T IS STXY/SXYT with imat>100)

   SAVE7END
 ************************************************************
  Group 8. Terms & Devices
 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  =         4 ; UMXCLS=     10000
 UNOLAY  =         2 ; UIOLEV=         2
 UWALEV  =         2 ; UGRITR=         2
 UMXCRS  =        16
 USPREL  = 7.000000E-01 ; UMNCEL= 1.000000E+00
 UDVREL  = 1.000000E+00
 ************************************************************
  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 1.0E-7
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(XY/T) =  1.001000E-10 ;FIINIT(SXYT) =  1.001000E-10
 FIINIT(SY/T) =  1.001000E-10 ;FIINIT(SYTH) =  1.001000E-10
 FIINIT(SX/T) =  1.001000E-10 ;FIINIT(SXTH) =  1.001000E-10
 FIINIT(TET7) =  1.001000E-10 ;FIINIT(R7  ) =  1.001000E-10
 FIINIT(EPSZ) =  1.001000E-10 ;FIINIT(EPSX) =  1.001000E-10
 FIINIT(EPSY) =  1.001000E-10 ;FIINIT(STXY) =  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
 ************************************************************
  Echo save-block settings for Group 11
 ************************************************************
  Group 12. Convection and diffusion adjustments
   No PATCHes used for this Group
 ************************************************************
  Group 13. Boundary & Special Sources
 
 EGWF    =    T
 ************************************************************
  Echo save-block settings for Group 13
   SAVE13BEGIN
 real(r02)
 r02 = R0**2
 ***** TOP : V - normal ******
 char(fV1,fV2,TV1,RV1)
 RV1=(:R02:/(:LY:^2+XG^2))
 TV1=ATAN(:LY:/XG)
 fV1=:FX:*(-:RV1:*(0.5*COS(2*:TV1:)-
 fV2=COS(4*:TV1:))-1.5*:RV1:^2*COS(4*:TV1:))
(source of DISY at TOP is COVAL(FIXFLU,:fV1::fV2:)!VOLU)
 
 ***** TOP : U - tangential ******
 char(fU1,fU2,TU1,RU1)
 RU1=(:R02:/(:LY:^2+XG^2))
 TU1=ATAN(:LY:/XG)
 fU1=:FX:*(-:RU1:*(0.5*SIN(2*:TU1:)+
 fU2=SIN(4*:TU1:))+1.5*:RU1:^2*SIN(4*:TU1:))
(source of DISX at TOP is COVAL(FIXFLU,:fU1::fU2:)!VOLU)

     (STORED VAR zzz7 IS :fU1::fU2:)

 
 ***** RIGHT : U - normal ******
 char(fU3,fU4,TU2,RU2,fZZ1,fZZ2)
 RU2=(:R02:/(:LX:^2+YG^2))
 TU2=ATAN(YG/:LX:)
 fU3=:FX:*(1-:RU2:*(1.5*COS(2*:TU2:)+
 fU4=COS(4*:TU2:))+1.5*:RU2:^2*COS(4*:TU2:))
(source of DISX at RIGHT is COVAL(FIXFLU,:fU3::fU4:)!VOLU)
       (source of DISX at RIGHT is COVAL(FIXFLU,-4.e7*YG)!VOLU)

(STORED VAR zzz7 IS :fU3::fU4:)


     fZZ1=:FX:*(1-:RU2:*(1.5*COS(2*:TU2:)
     fZZ2=+COS(4*:TU2:)))
    fZZ1=:RU2:*(1.5*COS(2*:TU2:)
     fZZ2=+COS(4*:TU2:))     
        
    (STORED VAR ZZZ7 IS :FZZ1::FZZ2:)
    (STORED VAR ZZZ7 IS :RU2:)   

 
 ***** RIGHT : V - tangential ******
 char(fV3,fV4,TV2,RV2)
 RV2=(:R02:/(:LX:^2+YG^2))
 TV2=ATAN(YG/:LX:)
 fV3=:FX:*(-:RV2:*(0.5*SIN(2*:TV2:)+
 fV4=SIN(4*:TV2:))+1.5*:RV2:^2*SIN(4*:TV2:))
(source of DISY at RIGHT is COVAL(FIXFLU,:fV3::fV4:)!VOLU)
 
 (source of DISX at LEFT is COVAL(1,0)!VOLU)
 (source of DISY at BOT is COVAL(1,0)!VOLU)
  
 
   SAVE13END
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =       800
 RESFAC  = 1.000000E-07
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 ************************************************************
  Group 18. Limits
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 STRA    =    T
 PARSOL  =    F
 ISG21   =       400
         ISG50   =         1
 ISG52   =         3
 CSG10   ='q1'
 SPEDAT(SET,MATERIAL,160,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 DISTIL=T
 EX(DISY)=   3.649E-04
 EX(DISX)=   1.184E-03
 EX(SXYT)=   1.296E+03
 EX(SY/T)=   1.234E+00
 EX(SYTH)=   1.445E+03
 EX(SX/T)=   1.055E+00
 EX(SXTH)=   9.872E+03
 EX(TET7)=   7.854E-01
 EX(R7  )=   1.351E+00
 EX(EPSZ)=   3.253E-04
 EX(EPSX)=   9.710E-04
 EX(EPSY)=   2.903E-04
 EX(STXY)=   1.327E+03
 EX(STRZ)=   2.444E-13
 EX(STRY)=   1.477E+03
 EX(STRX)=   9.801E+03
 EX(XY/T)=   1.206E+00
 EX(ZZZ7)=   9.338E+03

 ************************************************************
  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
 NXPRIN  =         3
 NYPRIN  =         4
 ISWPRF  =         1 ;ISWPRL =    100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 
 GVIEW(P,-1.445024E-01,2.191776E-01,9.649254E-01)
 GVIEW(UP,-3.081209E-02,9.736905E-01,-2.257833E-01)
 
> DOM,    SIZE,        2.000000E+00, 2.000000E+00, 5.000000E-01
> DOM,    MONIT,       1.300000E+00, 1.000000E-01, 2.500000E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    INCREMENT,   1.000000E-02, 1.000000E-02, 5.000000E-01
> GRID,   RSET_X_1,     10, 0.000000E+00
> GRID,   RSET_Y_1,     10, 1.000000E+00
  > GRID,   RSET_Z_1,      1, 1.000000E+00
 
> OBJ,    NAME,        LEFT
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 2.000000E+00, 5.000000E-01
> OBJ,    GEOMETRY,    cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
 
> OBJ,    NAME,        BOT
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        2.000000E+00, 0.000000E+00, 5.000000E-01
> OBJ,    GEOMETRY,    cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
 
> OBJ,    NAME,        TOP
> OBJ,    POSITION,    0.000000E+00, 2.000000E+00, 0.000000E+00
> OBJ,    SIZE,        2.000000E+00, 0.000000E+00, 5.000000E-01
> OBJ,    GEOMETRY,    cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
 
> OBJ,    NAME,        RIGHT
> OBJ,    POSITION,    2.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 2.000000E+00, 5.000000E-01
> OBJ,    GEOMETRY,    cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
 
> OBJ,    NAME,        HOLE
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        5.000000E-01, 5.000000E-01, 5.000000E-01
> OBJ,    GEOMETRY,    quarter-cylinder
> OBJ,    ROTATION24,        1
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    198,Solid with smooth-wall friction
 ************************************************************
STOP