TALK=f;RUN(1,1)
DISPLAY
Problem: The quarter of the square channel.
Inner pressure = 4.0e8, outer pressure = 0
1. 2D 'PLANE-STRAIN' : EpsZ = 0
2. Comparison - a calculation on code ELCUT (www.tor.ru)
3. Add Flow Cells for IX=1 and IY=NY
PHOTON USE
p;;;;
msg
msg TEM1 contours
msg
cont TEM1 z 1 x 1 20 y 2 21 fil;.0001
pause
set prop off
msg
msg V1 contours
msg
cont DELY z 1 x 1 20 y 2 21 fil;.0001
pause
cl
msg
msg U1 contours
msg
cont DELX z 1 x 1 20 y 2 21 fil;.0001
pause
cl
msg
msg EpsX contours
msg
cont EPSX z 1 x 1 20 y 2 21 fil;.0001
pause
cl
msg
msg EpsY contours
msg
cont EPSY z 1 x 1 20 y 2 21 fil;.0001
pause
cl
msg
msg EpsZ contours
msg
cont EPSZ z 1 x 1 20 y 2 21 fil;.0001
pause
cl
msg
msg STRX contours
msg
cont STRX z 1 x 1 20 y 2 21 fil;.0001
pause
cl
msg
msg STRY contours
msg
cont STRY z 1 x 1 20 y 2 21 fil;.0001
pause
cl
msg
msg STXY contours
msg
cont STXY z 1 x 1 20 y 2 21 fil;.0001
pause
cl
msg Displacements fields
gr ou z 1
cont DELT z 1 x 1 20 y 2 21 fil;.0001
vec z 1 x 1 20 y 2 21 col 0
pause
cl
AUTOPLOT
file
phi 5
ELDATA\str051.txt 2
ELDATA\dis051.txt 2
data 1
STRX X 14 Y 2 21
col1 1
data 2
y SXX
col2 2
scale
msg STRX(Y) for X=40 mm (WHITE-Phoenics, BLUE - Elcut) ...
plot 1 2
pause
clear
data 1
STRY X 15 Y 2 21
col1 1
data 2
y SYY
col2 2
scale
msg STRY(Y) for X=40 mm (WHITE-Phoenics, BLUE - Elcut) ...
plot 1 2
pause
clear
data 1
STXY X 15 Y 2 21
col1 1
data 2
y TXY
col2 2
scale
msg STXY(Y) for X=40 mm (WHITE-Phoenics, BLUE - Elcut) ...
plot 1 2
pause
clear
data 1
DELX X 15 Y 2 21
col1 1
data 3
y dx
col2 2
scale
msg U1(Y) for X=40 mm (WHITE-Phoenics, BLUE - Elcut) ...
plot 1 2
pause
clear
data 1
DELY X 15 Y 2 21
col1 1
data 3
y dy
col2 2
scale
msg V1(Y) for X=40 mm (WHITE-Phoenics, BLUE - Elcut) ...
plot 1 2
pause
clear
data 1
EPSX X 15 Y 2 21
col1 1
data 3
y exx
col2 2
scale
msg EpsX(Y) for X=40 mm (WHITE-Phoenics, BLUE - Elcut) ...
plot 1 2
pause
clear
data 1
EPSY X 15 Y 2 21
col1 1
data 3
y eyy
col2 2
scale
msg EpsY(Y) for X=40 mm (WHITE-Phoenics, BLUE - Elcut) ...
plot 1 2
pause
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
TEXT(square channel & TEM1;2D[xy])
integer(SolvMod)
mesg(SolvMod = 1 : ModModel361
mesg(SolvMod = 2 : StressModel
SolvMod=1
label ask
mesg(SolvMod=:SolvMod: Enter 1, 2 or blank
readvdu(SolvMod,int,1)
if(SolvMod.lt.1) then
goto ask
endif
if(SolvMod.gt.4) then
goto ask
endif
SolvMod
Declarations and settings
REAL(PIN,LX,LY,R0,POISSON,YOUNG)
PIN=4.0e8 ! H/m^2
LX=60.e-3
LY=60.e-3
R0=30.e-3
YOUNG = 1/0.5E-11 ! Young's modulus
POISSON=0.3 ! Poisson's ratio
INTEGER(NXC,NYC)
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NXC = 10
NREGX=4
IREGX=1;GRDPWR(X,1,0.001*LX,1.0) ! 18.12.05 single outer fluid cell
IREGX=2;GRDPWR(X,NXC,R0,1.0)
IREGX=3;GRDPWR(X,NXC,LX-R0,1)
IREGX=4;GRDPWR(X,1,0.001*LX,1.0) ! single outer fluid cell
************************************************************
Group 4. Y-Direction Grid Spacing
NYC = 10
NREGY=4
IREGY=1;GRDPWR(Y,1,0.001*LY,1.0) ! single outer fluid cell
IREGY=2;GRDPWR(Y,NYC,LY-R0,1.0)
IREGy=3;GRDPWR(Y,NYC,R0,1)
IREGY=4;GRDPWR(Y,1,0.001*LY,1.0) ! 18.12.05 single outer fluid cell
************************************************************
Group 5. Z-Direction Grid Spacing
NZ=1
ZWLAST = 0.001
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
SOLVE(P1,V1,U1,TEM1)
STORE(PRPS,DVO1,DRH1)
STORE(STRX,STRY,STRZ)
STORE(EPSY,EPSX,EPSZ)
STORE(STXY,EPST)
STORE(Delt,DelX,DelY)
************************************************************
GROUP 8. ITERATION NUMBERS ETC
RESFAC=1.e-7
RESREF(V1)=0.0
RESREF(U1)=0.0 ! to prevent premature exit
RESREF(TEM1)=0.0
LITER(V1) = 100 ! from solver
LITER(U1) = 100
LITER(P1) = 2
LITER(TEM1) = 100
************************************************************
GROUP 9. PROPERTIES
CSG10='Q1' ! materials with various POISSON ratios
MATFLG=T;NMAT=3
68 1.189 1.544E-05 1005.0 1000.0 0 0
69 1.189 1.544E-05 1005.0 1.e-20 0 0
160 7800.0 0.3 473.0 10.0 1.0e-4 0.5e-11
************************************************************
FFFFFFFFFFFFF00000000000 ny = 2*nyc+2
FFFFFFFFFFFFFuuuuuuuuuu0 ny-1
FFFFFFFFFFFFFssssssssss0
FFFFFFFFFFFFFssssssssss0
FFFFFFFFFFFFFssssssssss0
FFFFFFFFFFFFFssssssssss0 nyc+2
0lsssssssssssssssssssss0 nyc+1
0lsssssssssssssssssssss0
0lsssssssssssssssssssss0
0lsssssssssssssssssssss0 2
000000000000000000000000 1
1 nxc+1
nxc+2 nx = 2*nxc+2
F : inner ( P1 = PIN)
s : solid
u : solid + UP
l : solid + Left
0 : FU,FD,FL,FR : P1 = 0
************************************************************
GROUP 11. INITIAL VALUES
fiinit(u1)=0.0
fiinit(v1)=0.0
fiinit(TEM1)=0.0
fiinit(Delt)=0.0
==== Initial Value of PROP & P1 =======
FIINIT(PRPS)=160
FIINIT(P1)=0.0
PATCH(INNER,INIVAL,1,NXC+1,NYC+2,NY,1,1,1,1)
INIT(INNER,PRPS,FIXVAL,68)
INIT(INNER,P1,FIXVAL,PIN)
PATCH(FL,INIVAL,1,1,1,NY,1,1,1,1)
INIT(FL,PRPS,FIXVAL,69)
PATCH(FU,INIVAL,1,NX,NY,NY,1,1,1,1)
INIT(FU,PRPS,FIXVAL,69)
PATCH(FD,INIVAL,2,NX,1,1,1,1,1,1)
INIT(FD,PRPS,FIXVAL,68)
PATCH(FR,INIVAL,NX,NX,1,NY-1,1,1,1,1)
INIT(FR,PRPS,FIXVAL,68)
************************************************************
GROUP 13. BOUNDARY & SPECIAL SOURCES
PATCH(UP,north,NXC+2,NX-1,NY-1,NY-1,1,1,1,1) ! UP - fixed V1
COVAL(UP,V1,FIXVAL,0.0)
PATCH(LEFT,WEST,1,1,2,NYC+1,1,1,1,1) ! LEFT - fixed U1
COVAL(LEFT,U1,FIXVAL,0.0)
Fixed P1 in inner-zone
PATCH(PINNER,CELL,1,NXC+1,NYC+2,NY,1,1,1,1) ! pressurizing inner
COVAL(PINNER,P1,FIXVAL,PIN) ! fluid
COVAL(PINNER,TEM1,FIXVAL,0) ! TEM1 Normal Wall
Fixed U1 & V1 in inner-zone
PATCH(UINNER,CELL,1,NXC,NYC+2,NY,1,1,1,1)
COVAL(UINNER,U1,FIXVAL,0.0)
PATCH(VINNER,CELL,1,NXC+1,NYC+2,NY,1,1,1,1)
COVAL(VINNER,V1,FIXVAL,0.0)
Fixed P1 in FL/FD/FR/FU-zones
PATCH(PFL,CELL,1,1,1,NYC+1,1,1,1,1)
INIT(PFL,P1,FIXVAL,0)
PATCH(PFD,CELL,2,NX,1,1,1,1,1,1)
INIT(PFD,P1,FIXVAL,0)
INIT(PFD,TEM1,FIXVAL,100.0) ! TEM1 Hot Wall
PATCH(PFR,CELL,NX,NX,2,NY,1,1,1,1)
INIT(PFR,P1,FIXVAL,0)
COVAL(PFR,TEM1,FIXVAL,0) ! TEM1 Normal Wall
PATCH(PFU,CELL,NXC+2,NX-1,NY,NY,1,1,1,1)
INIT(PFU,P1,FIXVAL,0)
Fixed U1 in FL/FD/FR/FU-zones
PATCH(PFD,CELL,1,NX,1,1,1,1,1,1)
INIT(PFD,U1,FIXVAL,0)
PATCH(PFU,CELL,NXC+1,NX-1,NY,NY,1,1,1,1)
INIT(PFU,U1,FIXVAL,0)
Fixed V1 in FL/FD/FR/FU-zones
PATCH(VFL,CELL,1,1,1,NYC+1,1,1,1,1)
INIT(VFL,V1,FIXVAL,0)
PATCH(VFR,CELL,NX,NX,2,NY-1,1,1,1,1)
INIT(VFR,V1,FIXVAL,0)
! EpsZ = 0
SPEDAT(BOUNDARY,ZCONST,R,1.0e20)
************************************************************
GROUP 15. TERMINATE SWEEPS
LSWEEP = 200
ISG21=LSWEEP
************************************************************
GROUP 17. RELAXATION
#CONPROM
RELAX(P1 ,LINRLX, 1.000000E+00)
relax(U1,linrlx,1.0)
relax(V1,linrlx,1.0)
************************************************************
GROUP 19. DATA TRANSMITTED TO GROUND
STRA = T
PARSOL = F
ISG52 = 3 ! probe & res
************************************************************
GROUP 23.FIELD PRINT-OUT & PLOT CONTROL
TSTSWP = - 1 ! graphic-mode
NYPRIN = 1
IXMON = NX-2
IYMON = 2
IZMON = 1
#conprom
inform7begin
! DelX - U1, DelY - V2 in CENTRE scalar cells
! DelT - displacement in CENTRE scalar cells
(STORED VAR DelX IS 0.5*(U1+U1[-1]) with imat>100)
(STORED VAR DelY IS 0.5*(V1+V1[&-1]) with imat>100)
(STORED VAR Delt IS SQRT(DELX^2+DELY^2) with imat>100)
inform7end
if(SolvMod.EQ.1) then
spedat(rlxfac,rlxu1d,r,0.5)
spedat(rlxfac,rlxv1d,r,0.5)
LG(40) = F
else
LG(40) = T
endif
STOP