TALK=f;RUN(1,1)
DISPLAY
A hollow, thick-walled cylinder, immersed in an outer fluid,
contains a second fluid having a different pressure.
Recirculating square-wave cavity is made in internal surface of the cylinder.
Cylinder is cooled to temperature -30.
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 x 1 y 1 21 z 2 19 fil;.0001
vec x 1 y 2 21 z 2 19 col 1
dump
S225_V
pause
cl
msg
msg W1 contours
msg
cont W1 x 1 y 2 21 y 1 19 fil;.0001
vec x 1 y 2 21 z 2 19 col 1
dump
S225_W
pause
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
READQ1=T
begin_readq1
TEXT(Pressurised cylinder [R,Z] with Hole and not zero TEM
Declarations and settings
REAL(PI,PO,RI,RO,LZ,POISSON,YOUNG)
PI=1.e7 ! inner pressure
PO=1.e6 ! outer pressure
RI=0.2 ! inner radius
RO=0.8 ! outer radius
LZ = 1
YOUNG = 1/(0.5E-11) ! Young's modulus
POISSON=0.3 ! Poisson's ratio
integer(nycyl)
Grid settings
NX=1
NYCYL=20
NZ=20
**** !!! ***
ISG59 = 1
************************************************************
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = F
XULAST = 0.5
GRDPWR(X,NX,XULAST,1)
************************************************************
Group 4. Y-Direction Grid Spacing
NREGY=3 ! 3 regions
IREGY=1;GRDPWR(Y,1,RI,1.0) ! single inner fluid cell
IREGY=2;GRDPWR(Y,NYCYL,RO-RI,1) ! ncyl solid cells
IREGY=3;GRDPWR(Y,1,0.01*RO,1.0) ! single outer fluid cell
************************************************************
Group 5. Z-Direction Grid Spacing
ZWLAST = LZ
GRDPWR(Z,NZ,LZ,1.0)
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
SOLVE(P1,V1,W1)
STORE(P1TH,PRPS,DEN1,ENUL,DVO1,DRH1) ! TH means theoretical
STORE(STRX,EPSX,STRY)
STORE(EPSY,STRZ,EPSZ)
STORE(TEM1,EPST)
************************************************************
GROUP 8. ITERATION NUMBERS ETC
RESFAC=1.e-7
RESREF(V1)=0.0 ! to prevent premature exit
LITER(V1) = 100 ! from solver
RESREF(W1)=0.0; LITER(W1) = 100
************************************************************
GROUP 9. PROPERTIES
CSG10='Q1' ! materials with various POISSON ratios
MATFLG=T;NMAT=2
160 7800.0 0.3 473.0 43.0 1.0e-5 0.5E-11
68 1.0 1.E+20 1000.0 0.0 0 0
************************************************************
GROUP 11. INITIAL VALUES
FIINIT(PRPS)=160
PATCH(POUT,INIVAL,1,1,NY,NY,1,NZ,1,1)
INIT(POUT,PRPS,FIXVAL,0)
PATCH(PIN,INIVAL,1,1,1,1,1,NZ,1,1)
INIT(PIN,PRPS,FIXVAL,0)
PATCH(PHOLE,INIVAL,1,1,2,11,5,10,1,1)
INIT(PHOLE,PRPS,FIXVAL,0)
(initial of TEM1 at BLOCK is :formt:)
real(DTE_R,TE_R)
DTE_R = 0.0
TE_R = -30.0
(STORED VAR TEM1 IS (:TE_R:+:DTE_R:*(YG-0.2)) with imat>100)
-----------------------------------
PATCH(INNER1,CELL,1,1,1,1,1,NZ,1,1)
COVAL(INNER1,P1,FIXVAL,PI)
COVAL(INNER1,W1,FIXVAL,0)
PATCH(INNER2,CELL,1,1,2,11,5,9,1,1)
COVAL(INNER2,W1,FIXVAL,0)
PATCH(INNER3,CELL,1,1,1,10,5,10,1,1)
COVAL(INNER3,V1,FIXVAL,0)
PATCH(INNER,CELL,1,1,2,11,5,10,1,1)
COVAL(INNER,P1,FIXVAL,PI)
PATCH(OUTER,CELL,1,1,NY,NY,1,NZ,1,1)
COVAL(OUTER,P1,FIXVAL,PO)
COVAL(OUTER,W1,FIXVAL,0)
************************************************************
GROUP 13. BOUNDARY & SPECIAL SOURCES
PATCH(Z01,LOW,1,1,2,NY-1,1,1,1,1)
COVAL(Z01,W1,FIXVAL,0.0)
PATCH(Z02,HIGH,1,1,2,NY-1,NZ-1,NZ-1,1,1)
COVAL(Z02,W1,FIXVAL,0.0)
SPEDAT(BOUNDARY,XCONST,R,1.E20) ! because total angle is fixed
************************************************************
GROUP 15. TERMINATE SWEEPS
LSWEEP = 400
ISG21 = LSWEEP
ENDIT(V1) =0.0
RESREF(V1)=0.0
ENDIT(u1) =0.0
RESREF(u1)=0.0
************************************************************
GROUP 17. RELAXATION
#CONPROM
************************************************************
GROUP 19. DATA TRANSMITTED TO GROUND
STRA = T
************************************************************
GROUP 23.FIELD PRINT-OUT & PLOT CONTROL
TSTSWP = - 1 ! graphic-mode
ISG52 = 3 ! probe & res
output(den1,n,n,n,n,n,n)
output(drh1,n,n,n,n,n,n)
output(dvo1,n,n,n,n,n,n)
output(enul,n,n,n,n,n,n)
output(prps,n,n,n,n,n,n)
NYPRIN=1
NZPRIN=1
IXMON = 1
IYMON = NY-1
IZMON = NZ/2
ISG52 = 3 ! probe & res
spedat(rlxfac,rlxw1d,r,0.5)
spedat(rlxfac,rlxv1d,r,0.5)
(stored var v1an is anco(v1))
(stored var v1as is asco(v1))
(stored var v1ah is ahco(v1))
(stored var v1al is alco(v1))
(stored var v1rs is resi(v1))
(stored var w1an is anco(w1))
(stored var w1as is asco(w1))
(stored var w1ah is ahco(w1))
(stored var w1al is alco(w1))
(stored var w1rs is resi(w1))
STOP