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; us202 )
************************************************************
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
************************************************************
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