TALK=f;RUN(1,1) Square channel & TEM; s254

  DISPLAY
  Problem: One quarter of a square pressurised channel.
           Inner pressure = 4.0e8, outer pressure = 0
   1. 2D 'PLANE-STRAIN' : z-surface fixed
   2. Comparison - a calculation on code ELCUT (www.tor.ru)
   3. Thin fluid cells are provided at IX=1 and IY=NY
  ENDDIS
  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 
  EDZZ\S254_1.txt 2
  EDZZ\S254_2.txt 2

  data 1                                                                          
  STRX X 14 Y 2 21                                                                   
  col1 1          
  data 2                                                                          
  y SXX 
  colA 2          
  scale
  msg STRX(Y) for X=40 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2                                                                      
  pause 

  clear
  data 1                                                                          
  STRY X 14 Y 2 21                                                                   
  col1 1          
  data 2                                                                          
  y SYY 
  colA 2          
  scale
  msg STRY(Y) for X=40 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2                                                                      
  pause 

  clear
  data 1                                                                          
  STXY X 14 Y 2 21                                                                   
  col1 1          
  data 2                                                                          
  y TXY 
  colA 2          
  scale
  msg STXY(Y) for X=40 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2                                                                      
  pause 

  clear  
  data 1                                                                          
  DELX X 14 Y 2 21                                                                   
  col1 1          
  data 3                                                                          
  y dx 
  colA 2          
  scale
  msg U1(Y) for X=40 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2                                                                      
  pause 

  clear  
  data 1                                                                          
  DELY X 14 Y 2 21                                                                   
  col1 1          
  data 3                                                                          
  y dy 
  colA 2          
  scale
  msg V1(Y) for X=40 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2                                                                      
  pause 

  clear
  data 1                                                                          
  EPSX X 14 Y 2 21                                                                   
  col1 1          
  data 3                                                                          
  y exx 
  colA 2          
  scale
  msg EpsX(Y) for X=40 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2                                                                      
  pause 

  clear
  data 1                                                                          
  EPSY X 14 Y 2 21                                                                   
  col1 1          
  data 3                                                                          
  y eyy 
  colA 2          
  scale
  msg EpsY(Y) for X=40 mm (WHITE-Phoenics, YELLOW - Elcut) ...
  plot 1 2                                                                      
  pause 


  ENDUSE

  
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 TEXT(square channel & TEM1;2D[xy])
 
  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(BODY,INIVAL,1,NX,1,NY,1,1,1,1)
 INIT(BODY,PRPS,FIXVAL,160)

 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
 #maxmin

 ************************************************************
  GROUP 23.FIELD PRINT-OUT & PLOT CONTROL
TSTSWP = - 1   ! graphic-mode
NYPRIN = 1
NXPRIN = 1
IXMON = NX-2
IYMON = 2
IZMON = 1

  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

STOP