TALK=f;RUN(1,1)
  DISPLAY
    Stresses are considered in rotating disk with StrZ = 0:
     Omega   - Angular velocity of the rotation,
     DenBody - density,
     RO      - radius of disk. 
    Exact solution (Timoshenko, paragraph 32):
    
     V = DenBody*Omega^2*RO^3 * (1-P)/8/E * { R1*(3+P) - 
         R1^3*(1+P) },
     R1 = R/RO    
  
  ENDDIS
  
    PHOTON USE
  p;;;;
  
  

  set prop off
  cl
  msg PRPS field
  gr ou z 1
  con prps z 1  fi;1.0
  
  pause
  cl
  msg x-displacement field
  gr ou z 1
  con U1 z 1  fi;.0001
  vec z 1  col 0
  
  pause
  cl
  msg y-displacement field
  gr ou z 1
  con V1 z 1 fi;.0001
  
  pause
  cl
  msg epsx field
  gr ou z 1
  con epsx z 1 fi;.0001
  
  pause
  cl
  msg epsy field
  gr ou z 1
  con epsy z 1 fi;.0001

  pause
  cl
  msg strx field
  gr ou z 1
  con strx z 1 fi;.0001
  
  pause
  cl
  msg stry field
  gr ou z 1
  con stry z 1 fi;.0001

  ENDUSE

 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
readq1=t
  READQ1_BEGIN 
 TEXT(Stress in a rotating disk; 2D[x,y]
  Declarations and settings
REAL(DenBody,RO,OMEGA,POISSON,YOUNG) 
DenBody = 7000 
OMEGA =100 * 6.28
RO = 0.1 
YOUNG   = 1/(0.5E-11) ! Young's modulus
POISSON=0.3           ! Poisson's ratio
INTEGER(CASENO,NYCYL) ! nycyl is y-direction cell number in solid
NYCYL=31

real(FVY)
FVY = DenBody*Omega**2


  READQ1_END 
  
  **** !!! ***
ISG59 = 1
  
  
 ************************************************************
  Group 2. Time dependence
 STEADY  =    T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES  =    F
 NX=10
 GRDPWR(X,NX,1.0,1)
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 GRDPWR(Y,NYCYL,RO,1)
 
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 NZ=1
 ZWLAST  = 1.000000E+00
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
 SOLVE(P1,V1,U1)
 STORE(PRPS,DEN1,ENUL,DVO1,DRH1)
 STORE(STRX,EPSX,STRY,EPSY,EPSZ)
 STORE(V1TH,V1/T)
 ************************************************************
  GROUP 8. ITERATION NUMBERS ETC
 RESREF(V1)=0.0   ! to prevent premature exit
 LITER(V1) = 20   ! from solver
 ************************************************************
  GROUP 9. PROPERTIES
  
 CSG10='Q1'                  ! materials with various POISSON ratios
  MATFLG=T;NMAT=1         
  160    7800.0    0.3       473.0   43.0      1.0e-5   0.5E-11 
  161    7800.0    0.2       473.0   43.0      1.0e-5   0.5E-11 
  162    7800.0    0.1       473.0   43.0      1.0e-5   0.5E-11 
  163    7800.0    0.0       473.0   43.0      1.0e-5   0.5E-11 
  164    7800.0    0.4       473.0   43.0      1.0e-5   0.5E-11 
  165    7800.0    0.05      473.0   43.0      1.0e-5   0.5E-11 
  166    7800.0    0.01      473.0   43.0      1.0e-5   0.5E-11 
 
 ************************************************************
  GROUP 11. INITIAL VALUES
FIINIT(PRPS)=160

 ************************************************************
  GROUP 13. BOUNDARY & SPECIAL SOURCES
 
PATCH(BODY,VOLUME,1,NX,1,NY-1,1,1,1,1)
(SOURCE of V1 at BODY is COVAL(FIXFLU,:FVY:*RV))

patch(lowx,wwall,1,1,1,NY,1,1,1,1)
coval(lowx,u1,fIXval,0.0)
patch(highx,ewall,nx-1,nx-1,1,NY,1,1,1,1)
coval(highx,u1,fixval,0.0)

spedat(rlxfac,rlxV1d,r,0.5) 
 
 ************************************************************
  GROUP 15. TERMINATE SWEEPS
 LSWEEP  =  200
 ISG21   =   LSWEEP
 ENDIT(V1) =0.0
 RESFAC=1.e-7 
 RESREF(V1)=0.0                                   
 ************************************************************
  GROUP 17. RELAXATION
#CONPROM
 ************************************************************
  GROUP 19. DATA TRANSMITTED TO GROUND
 STRA    =    T
 ************************************************************
  GROUP 23.FIELD PRINT-OUT & PLOT CONTROL
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)  
TSTSWP=-1  
IYPRF=1
IYPRL=NY

ISG52   =    3   !   probe & res

 
  inform7begin       
  ! Declarations first:
REAL(CC0, CC1, CC2)
CC0 = FVY*RO**3*(1-POISSON)/YOUNG/8
CC1 = (3+POISSON)/RO
CC2 = (1+POISSON)/RO**3
char(formV)
formV=:CC0:*(RV*:CC1:-:CC2:*RV^3)
(STORED VAR V1TH IS :formV:)
(STORED VAR V1/T IS V1/V1TH)
  inform7end
STOP