```
TALK=f;RUN(1,1)
DISPLAY
Stresses are computed for a rotating disk with zero
axial-direction (z) stress
Omega   - Angular velocity of the rotation,
DenBody - density,
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
************************************************************
Group 1. Run Title and Number
************************************************************
TEXT(Stress in a rotating disk; 1Dy
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

************************************************************
Group 2. Time dependence
************************************************************
Group 3. X-Direction Grid Spacing
CARTES  =    F
NX=1
XULAST  = 1.000000E-03
XFRAC (   1) =  1.000000E+00
************************************************************
Group 4. Y-Direction Grid Spacing
GRDPWR(Y,NYCYL,RO,1)

************************************************************
Group 5. Z-Direction Grid Spacing
NZ=1
ZWLAST  = 1.000000E+00

**** !!! ***
ISG59 = 1

************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS  =    T
SOLVE(P1,V1)
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,1,1,NY-1,1,1,1,1)
(SOURCE of V1 at BODY is COVAL(FIXFLU,:FVY:*RV))

SPEDAT(BOUNDARY,XCONST,R,1.E20)
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)
NYPRIN=1
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
```