TALK=T;RUN(1,1)
PHOTON USE
p
gr ou z 1
con thex z 1 fil;.001
gr ou z 1 x 1 3 y 6 m
gr ou z 1 x 8 m y 6 m
msg Linear thermal expansion coefficient
pause
con cl;red
con stif z 1 fil;.001
gr ou z 1 x 1 3 y 6 m
gr ou z 1 x 8 m y 6 m
msg Stiffnes ( Young's module ) contours
pause
con cl;red
vec z 1 sh
msg Velocity and displacement vectors
pause
vec cl;red
con u1 z 1 fil;.001
gr ou z 1 x 1 3 y 6 m
gr ou z 1 x 8 m y 6 m
msg X-component displacement contours
pause
con cl;red
con v1 z 1 fil;.001
gr ou z 1 x 1 3 y 6 m
gr ou z 1 x 8 m y 6 m
msg Y-component displacement contours
pause
con cl;red
con dil z 1 fil;.001
gr ou z 1 x 1 3 y 6 m
gr ou z 1 x 8 m y 6 m
msg Dilatation contours
pause
con cl;red
con epsx z 1 fil;.001
gr ou z 1 x 1 3 y 6 m
gr ou z 1 x 8 m y 6 m
msg X-strain contours
pause
con cl;red
con epsy z 1 fil;.001
gr ou z 1 x 1 3 y 6 m
gr ou z 1 x 8 m y 6 m
msg Y-strain contours
pause
con cl;red
con epst z 1 fil;.001
gr ou z 1 x 1 3 y 6 m
gr ou z 1 x 8 m y 6 m
msg Thermal strain contours
pause
con cl;red
con strx z 1 fil;.001
gr ou z 1 x 1 3 y 6 m
gr ou z 1 x 8 m y 6 m
msg X-normal stress contours
pause
con cl;red
con stry z 1 fil;.001
gr ou z 1 x 1 3 y 6 m
gr ou z 1 x 8 m y 6 m
msg Y-normal stress contours
pause
con cl;red
con taxy z 1 fil;.001
gr ou z 1 x 1 3 y 6 m
gr ou z 1 x 8 m y 6 m
msg Shear stress contours
pause
msg Press E to end
enduse
DISPLAY
This case, deals with the calculations of
the displacements and stresses
in uniformly heated block composed the materials of
different mechanical properties.
It is used here as an example of extensive processing
the calculation results with the aim to get the
distribution of normal and shear stresses by the
numerical differentiation of displacement field.
The outcome is the stress fields which may be viewed by
PHOTON. The use file for the latter is also supplied.
ENDDIS
PLANTBEGIN
PATCH(STELC,INIVAL,1,NX,1,NY,1,1,1,1)
VAL=(1.-RG(4))/((1.+RG(4))*(1-2.*RG(4)))
INIT (STELC,CON1,0.0,GRND)
VAL=RG(4)/((1.+RG(4))*(1-2.*RG(4)))
INIT (STELC,CON2,0.0,GRND)
VAL=1./(1.-2.*RG(4))
INIT (STELC,CON3,0.0,GRND)
VAL=1./(2.*(1+RG(4)))
INIT (STELC,CON4,0.0,GRND)
Stress - in -Solid BC
---------------- WEST ----------------------------
** Zero displacement at the west boundary
PATCH(WESZD,WWALL,1,1,1,NY,1,1,1,1)
COVAL(WESZD,U1,1.,0.0)
---------------- EAST ----------------------------
** Normal stress at the east solid boundary
PATCH(EASZN,EAST,NX-1,NX-1,1,NY,1,1,1,1)
VAL=U1[-1,,]+0.333*(P1+4.*(THEX*TEM+RG(3)/STIF))*DXU2D
COVAL(EASZN,U1,FIXVAL,GRND)
** Zero shear stress at the east boundary
PATCH(EASZS,EAST,NX,NX,1,NY-1,1,1,1,1)
VAL=-(U1[-1,+1,]-U1[-1,,])/DYG2D
COVAL(EASZS,V1,FIXFLU,GRND)
** Dilatation at the east solid boundary
PATCH(FXDEAS,EAST,NX,NX,2,NY-1,1,1,1,1)
VAL=2.*(THEX*TEM-RG(3)/STIF)-1.5*(V1-V1[,-1,])/DYV2D
COVAL(FXDEAS,P1,FIXVAL, grnd)
PATCH(FXDEAC,EAST,NX,NX,1,1,1,1,1,1)
VAL=2.*(THEX*TEM-RG(3)/STIF)- 1.5*(V1[,+1,]-$
V1)/DYV2D[,+1,]
COVAL(FXDEAC,P1,FIXVAL, grnd)
---------------- NORTH ----------------------------
** Normal stress at the free north boundary
PATCH(NORZN,CELL,1,NX,NY-1,NY-1,1,1,1,1)
VAL=V1[,-1,]+0.333*(P1+4.*(THEX*TEM+RG(2)/STIF))*DYV2D
COVAL(NORZN,V1,FIXVAL,grnd)
** Zero shear stress at the free north boundary
PATCH(NORZS,NORTH,1,NX-1,NY,NY,1,1,1,1)
VAL=-(V1[+1,-1,]-V1[,-1,])/DXG2D
COVAL(NORZS,U1,FIXFLU,GRND)
** Dilatation at the free north boundary
PATCH(FDNW,NORTH,1,1,NY,NY,1,1,1,1)
VAL=2.*(THEX*TEM-RG(2)/STIF)- 1.5*(U1[+1,,]-$
U1)/DXU2D[+1,,]
COVAL(FDNW,P1,FIXVAL, grnd)
PATCH(FDNOR,NORTH,2,NX-1,NY,NY,1,1,1,1)
VAL=2.*(THEX*TEM-RG(2)/STIF)-1.5*(U1-U1[-1,,])/DXU2D
COVAL(FDNOR,P1,FIXVAL, grnd)
PATCH(FDNEC,NORTH,NX,NX,NY,NY,1,1,1,1)
VAL=2.*(THEX*TEM-RG(2)/STIF)- 1.5*(U1[-1,,]-U1[-2$
,,])/DXU2D[-1,,]
COVAL(FDNEC,P1,FIXVAL, grnd)
---------------- SOUTH ----------------------------
** Zero displacement at the south symmetry
PATCH(SOUZD,SWALL,1,NX,1,1,1,1,1,1)
COVAL(SOUZD,V1,1.,0.0)
==================================================
** Volumetric mass source for dilatation
PATCH(DILAT,VOLUME,1,nx,1,ny,1,1,1,1)
VAL=-P1/CON3+THEX*TEM/CON4
COVAL(DILAT,P1,FIXFLU,GRND)
===================================================
OUTPUT INFORMATION
** Dilatation
STORE(DIL)
DIL=-P1/CON3+THEX*TEM/CON4
** Stress and strain calculations
STORE(STRX,STRY,TAXY)
STORE(EPSX,EPSY,EPST)
* Thermal strains
EPST=THEX*TEM
REGION(1,NX,1,NY,1,1)
* Normal strains
EPSX=(U1-WEST(U1))/DXU2D
REGION(2,NX-1,1,NY,1,1)
EPSX=(U1[+1,,]-U1)/DXU2D[+1,,]
REGION(1,1,1,NY,1,1)
EPSX=(U1[-1,,]-U1[-2,,])/DXU2D[-1,,]
REGION(NX,NX,1,NY,1,1)
EPSY=(V1-SOUTH(V1))/DYV2D
REGION(1,NX,2,NY-1,1,1)
EPSY=(V1[,+1,]-V1)/DYV2D[,+1,]
REGION(1,NX,1,1,1,1)
EPSY=(V1[,-1,]-V1[,-2,])/DYV2D[,-1,]
REGION(1,NX,NY,NY,1,1)
* Normal stresses
STRX=STIF*(CON1*EPSX+CON2*EPSY-CON3*EPST)
REGION(1,NX,1,NY,1,1)
STRY=STIF*(CON1*EPSY+CON2*EPSX-CON3*EPST)
REGION(1,NX,1,NY,1,1)
* Shear stresses
STORE(DDYU,DDXV)
== ddy.u ==
DDYU=0.25*((U1-U1[,-1,])/DYG2D[,-1,]+$
(U1[,+1,]-U1)/DYG2D+(U1[-1,,]- U1[-1,-1,])$
/DYG2D[,-1,]+ (U1[-1,+1,]-U1[-1,,])/DYG2D)
REGION(2,NX-1,2,NY-1,1,1)
DDYU=0.5*((U1-U1[,-1,])/DYG2D[,-1,]+ (U1[$
,+1,]-U1)/DYG2D)
REGION(1,1,2,NY-1,1,1)
DDYU=(U1[,+1,]-U1)/DYG2D
REGION(1,1,1,1,1,1)
DDYU=(U1-U1[,-1,])/DYG2D[,-1,]
REGION(1,1,NY,NY,1,1)
DDYU=0.5*((U1[-1,,]-U1[-1,-1,])/DYG2D[,-1,]+$
(U1[-1,+1,]-U1[-1,,])/DYG2D)
REGION(NX,NX,2,NY-1,1,1)
DDYU=(U1[-1,+1,]-U1[-1,,])/DYG2D
REGION(NX,NX,1,1,1,1)
DDYU=(U1[-1,,]-U1[-1,-1,])/DYG2D[,-1,]
REGION(NX,NX,NY,NY,1,1)
DDYU=0.5*((U1-U1[,-1,])/DYG2D[,-1,]+ (U1[-1,$
,]-U1[-1,-1,])/DYG2D[,-1,])
REGION(2,NX-1,NY,NY,1,1)
DDYU=0.5*((U1[,+1,]-U1)/DYG2D+ (U1[-1,+1,]-U$
1[-1,,])/DYG2D)
REGION(2,NX-1,1,1,1,1)
== ddx.v ==
DDXV=0.25*((V1[+1,-1,]-V1[,-1,] )/DXG2D+$
(V1[+1,,]-V1)/DXG2D+(V1[,-1,]- V1$
[-1,-1,])/DXG2D[-1,,]+ (V1-V1[-1,,])/DXG2D$
[-1,,])
REGION(2,NX-1,2,NY-1,1,1)
DDXV=(V1[+1,,] -V1)/DXG2D
REGION(1,1,1,1,1,1)
DDXV=0.5*((V1[+1,-1,]-V1[,-1,])/DXG2D+$
(V1[+1,,]-V1)/DXG2D)
REGION(1,1,2,NY-1,1,1)
DDXV=(V1[+1,-1,]-V1[,-1,] )/DXG2D
REGION(1,1,NY,NY,1,1)
DDXV=(V1-V1[-1,,])/DXG2D[-1,,]
REGION(NX,NX,1,1,1,1)
DDXV=0.5*((V1[,-1,]-V1[-1,-1,])/DXG2D[-1,,]+$
(V1-V1[-1,,])/DXG2D[-1,,])
REGION(NX,NX,2,NY-1,1,1)
DDXV=(V1[,-1,] -V1[-1,-1,])/DXG2D[-1,,]
REGION(NX,NX,NY,NY,1,1)
DDXV=0.5*((V1[+1,,]-V1)/DXG2D+ (V1-V1[-1,,$
])/DXG2D[-1,,])
REGION(2,NX-1,1,1,1,1)
DDXV=0.5*((V1[+1,-1,]-V1[,-1,])/DXG2D+ (V1$
[,-1,]-V1[-1,-1,])/DXG2D[-1,,])
REGION(2,NX-1,NY,NY,1,1)
TAXY=STIF*CON4*(DDYU+DDXV)
REGION(1,NX,1,NY,1,1)
PLANTEND
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(2D Stress in solid, XY plane )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 611
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 10
XULAST =1.
XFRAC(1)=0.1 ;XFRAC(2)=0.2
XFRAC(3)=0.3 ;XFRAC(4)=0.4
XFRAC(5)=0.5 ;XFRAC(6)=0.6
XFRAC(7)=0.7 ;XFRAC(8)=0.8
XFRAC(9)=0.9 ;XFRAC(10)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 10
YVLAST =1.
YFRAC(1)=0.1 ;YFRAC(2)=0.2
YFRAC(3)=0.3 ;YFRAC(4)=0.4
YFRAC(5)=0.5 ;YFRAC(6)=0.6
YFRAC(7)=0.7 ;YFRAC(8)=0.8
YFRAC(9)=0.9 ;YFRAC(10)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 1
ZWLAST =1.
ZFRAC(1)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(3)=U1
NAME(5)=V1 ;NAME(134)=DDXV
NAME(135)=DDYU ;NAME(136)=EPST
NAME(137)=EPSY ;NAME(138)=EPSX
NAME(139)=TAXY ;NAME(140)=STRY
NAME(141)=STRX ;NAME(142)=DIL
NAME(143)=CON4 ;NAME(144)=CON3
NAME(145)=CON2 ;NAME(146)=CON1
NAME(147)=POIS ;NAME(148)=STIF
NAME(149)=THEX ;NAME(150)=TEM
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(P1,Y,Y,N,N,N,Y)
SOLUTN(U1,Y,Y,N,N,N,Y)
SOLUTN(V1,Y,Y,N,N,N,Y)
SOLUTN(DDXV,Y,N,N,N,N,Y)
SOLUTN(DDYU,Y,N,N,N,N,Y)
SOLUTN(EPST,Y,N,N,N,N,Y)
SOLUTN(EPSY,Y,N,N,N,N,Y)
SOLUTN(EPSX,Y,N,N,N,N,Y)
SOLUTN(TAXY,Y,N,N,N,N,Y)
SOLUTN(STRY,Y,N,N,N,N,Y)
SOLUTN(STRX,Y,N,N,N,N,Y)
SOLUTN(DIL,Y,N,N,N,N,Y)
SOLUTN(CON4,Y,N,N,N,N,Y)
SOLUTN(CON3,Y,N,N,N,N,Y)
SOLUTN(CON2,Y,N,N,N,N,Y)
SOLUTN(CON1,Y,N,N,N,N,Y)
SOLUTN(POIS,Y,N,N,N,N,Y)
SOLUTN(STIF,Y,N,N,N,N,Y)
SOLUTN(THEX,Y,N,N,N,N,Y)
SOLUTN(TEM,Y,Y,N,N,N,Y)
************************************************************
Group 8. Terms & Devices
* Y in TERMS argument list denotes:
* 1-built-in source 2-convection 3-diffusion 4-transient
* 5-first phase variable 6-interphase transport
TERMS(P1,Y,Y,Y,N,Y,Y)
TERMS(U1,Y,N,Y,Y,Y,Y)
TERMS(V1,Y,N,Y,Y,Y,Y)
TERMS(TEM,N,N,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 =1. ;TMP1 =0. ;EL1 =0.
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =1. ;ENUT =0.
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRNDTL(TEM)=1.
PRT(U1)=1. ;PRT(V1)=1.
PRT(TEM)=1.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.0E-10
FIINIT(V1)=1.0E-10 ;FIINIT(DDXV)=1.0E-10
FIINIT(DDYU)=1.0E-10 ;FIINIT(EPST)=1.0E-10
FIINIT(EPSY)=1.0E-10 ;FIINIT(EPSX)=1.0E-10
FIINIT(TAXY)=1.0E-10 ;FIINIT(STRY)=1.0E-10
FIINIT(STRX)=1.0E-10 ;FIINIT(DIL)=1.0E-10
FIINIT(CON4)=1.0E-10 ;FIINIT(CON3)=1.0E-10
FIINIT(CON2)=1.0E-10 ;FIINIT(CON1)=1.0E-10
FIINIT(POIS)=1.0E-10 ;FIINIT(STIF)=1.0E-10
FIINIT(THEX)=1.0E-10 ;FIINIT(TEM)=100.
PATCH(STELP ,INIVAL, 1, 10, 1, 10, 1, 1, 1, 1)
INIT(STELP ,POIS,0. ,0.3333 )
PATCH(STELT ,INIVAL, 1, 10, 1, 10, 1, 1, 1, 1)
INIT(STELT ,THEX,0. ,1.0E-05 )
PATCH(STELS ,INIVAL, 1, 10, 1, 10, 1, 1, 1, 1)
INIT(STELS ,STIF,0. ,2.0E+11 )
PATCH(ALMT1 ,INIVAL, 1, 3, 6, 10, 1, 1, 1, 1)
INIT(ALMT1 ,STIF,0. , FIXVAL )
INIT(ALMT1 ,THEX,0. ,1.0E-04 )
PATCH(ALMT2 ,INIVAL, 8, 10, 6, 10, 1, 1, 1, 1)
INIT(ALMT2 ,STIF,0. , FIXVAL )
INIT(ALMT2 ,THEX,0. ,1.0E-04 )
PATCH(STELC ,INIVAL, 1, 10, 1, 10, 1, 1, 1, 1)
INIT(STELC ,CON4,0. , GRND )
INIT(STELC ,CON3,0. , GRND )
INIT(STELC ,CON2,0. , GRND )
INIT(STELC ,CON1,0. , GRND )
INIADD = F
FSWEEP = 1
NAMFI =CHAM
************************************************************
Group 12. Patchwise adjustment of terms
Patches for this group are printed with those
for Group 13.
Their names begin either with GP12 or &
************************************************************
Group 13. Boundary & Special Sources
PATCH(WESZD ,WWALL , 1, 1, 1, 10, 1, 1, 1, 1)
COVAL(WESZD ,U1 ,1. ,0. )
PATCH(EASZN ,EAST , 9, 9, 1, 10, 1, 1, 1, 1)
COVAL(EASZN ,U1 , FIXVAL , GRND )
PATCH(EASZS ,EAST , 10, 10, 1, 9, 1, 1, 1, 1)
COVAL(EASZS ,V1 , FIXFLU , GRND )
PATCH(FXDEAS ,EAST , 10, 10, 2, 9, 1, 1, 1, 1)
COVAL(FXDEAS ,P1 , FIXVAL , GRND )
PATCH(FXDEAC ,EAST , 10, 10, 1, 1, 1, 1, 1, 1)
COVAL(FXDEAC ,P1 , FIXVAL , GRND )
PATCH(NORZN ,CELL , 1, 10, 9, 9, 1, 1, 1, 1)
COVAL(NORZN ,V1 , FIXVAL , GRND )
PATCH(NORZS ,NORTH , 1, 9, 10, 10, 1, 1, 1, 1)
COVAL(NORZS ,U1 , FIXFLU , GRND )
PATCH(FDNW ,NORTH , 1, 1, 10, 10, 1, 1, 1, 1)
COVAL(FDNW ,P1 , FIXVAL , GRND )
PATCH(FDNOR ,NORTH , 2, 9, 10, 10, 1, 1, 1, 1)
COVAL(FDNOR ,P1 , FIXVAL , GRND )
PATCH(FDNEC ,NORTH , 10, 10, 10, 10, 1, 1, 1, 1)
COVAL(FDNEC ,P1 , FIXVAL , GRND )
PATCH(SOUZD ,SWALL , 1, 10, 1, 1, 1, 1, 1, 1)
COVAL(SOUZD ,V1 ,1. ,0. )
PATCH(DILAT ,VOLUME, 1, 10, 1, 10, 1, 1, 1, 1)
COVAL(DILAT ,P1 , FIXFLU , GRND )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 3000 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = F
RESREF(P1)=1.0E-08 ;RESREF(U1)=1.0E-08
RESREF(V1)=1.0E-08 ;RESREF(TEM)=1.0E-08
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(U1)=10
LITER(V1)=10 ;LITER(TEM)=20
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=1.0E-03 ;ENDIT(TEM)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,5.0E-03)
RELAX(U1,FALSDT,5.0E-04)
RELAX(V1,FALSDT,5.0E-04)
RELAX(DDXV,LINRLX,1.)
RELAX(DDYU,LINRLX,1.)
RELAX(EPST,LINRLX,1.)
RELAX(EPSY,LINRLX,1.)
RELAX(EPSX,LINRLX,1.)
RELAX(TAXY,LINRLX,1.)
RELAX(STRY,LINRLX,1.)
RELAX(STRX,LINRLX,1.)
RELAX(DIL,LINRLX,1.)
RELAX(CON4,LINRLX,1.)
RELAX(CON3,LINRLX,1.)
RELAX(CON2,LINRLX,1.)
RELAX(CON1,LINRLX,1.)
RELAX(POIS,LINRLX,1.)
RELAX(STIF,LINRLX,1.)
RELAX(THEX,LINRLX,1.)
RELAX(TEM,FALSDT,1.0E+09)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(DDXV)=1.0E+10 ;VARMIN(DDXV)=-1.0E+10
VARMAX(DDYU)=1.0E+10 ;VARMIN(DDYU)=-1.0E+10
VARMAX(EPST)=1.0E+10 ;VARMIN(EPST)=-1.0E+10
VARMAX(EPSY)=1.0E+10 ;VARMIN(EPSY)=-1.0E+10
VARMAX(EPSX)=1.0E+10 ;VARMIN(EPSX)=-1.0E+10
VARMAX(TAXY)=1.0E+10 ;VARMIN(TAXY)=-1.0E+10
VARMAX(STRY)=1.0E+10 ;VARMIN(STRY)=-1.0E+10
VARMAX(STRX)=1.0E+10 ;VARMIN(STRX)=-1.0E+10
VARMAX(DIL)=1.0E+10 ;VARMIN(DIL)=-1.0E+10
VARMAX(CON4)=1.0E+10 ;VARMIN(CON4)=-1.0E+10
VARMAX(CON3)=1.0E+10 ;VARMIN(CON3)=-1.0E+10
VARMAX(CON2)=1.0E+10 ;VARMIN(CON2)=-1.0E+10
VARMAX(CON1)=1.0E+10 ;VARMIN(CON1)=-1.0E+10
VARMAX(POIS)=1.0E+10 ;VARMIN(POIS)=-1.0E+10
VARMAX(STIF)=1.0E+10 ;VARMIN(STIF)=-1.0E+10
VARMAX(THEX)=1.0E+10 ;VARMIN(THEX)=-1.0E+10
VARMAX(TEM)=1.0E+10 ;VARMIN(TEM)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
NAMSAT =MOSG
PARSOL = F
ISG62 = 1
POISSN =0.3333
SPEDAT(SET,GXMONI,PLOTALL,L,T)
RG( 1) =1.0E-05
RG( 2) =-1.0E+08
RG( 4) =0.3333
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=4.671E-03 ;EX(U1)=3.079E-03
EX(V1)=6.59E-04 ;EX(DDXV)=2.347E-03
EX(DDYU)=3.255E-03 ;EX(EPST)=3.7E-03
EX(EPSY)=3.955E-03 ;EX(EPSX)=6.382E-03
EX(TAXY)=1.362E+08 ;EX(STRY)=3.664E+08
EX(STRX)=1.83E+08 ;EX(DIL)=8.618E-03
EX(CON4)=0.375 ;EX(CON3)=2.999
EX(CON2)=0.7498 ;EX(CON1)=1.5
EX(POIS)=0.3333 ;EX(STIF)=1.46E+11
EX(THEX)=3.7E-05 ;EX(TEM)=100.
************************************************************
Group 21. Print-out of Variables
INIFLD = F ;SUBWGR = F
* Y in OUTPUT argument list denotes:
* 1-field 2-correction-eq. monitor 3-selective dumping
* 4-whole-field residual 5-spot-value table 6-residual table
OUTPUT(P1,Y,N,Y,Y,Y,Y)
OUTPUT(U1,Y,N,Y,Y,Y,Y)
OUTPUT(V1,Y,N,Y,Y,Y,Y)
OUTPUT(DDXV,Y,N,Y,N,N,N)
OUTPUT(DDYU,Y,N,Y,N,N,N)
OUTPUT(EPST,Y,N,Y,N,N,N)
OUTPUT(EPSY,Y,N,Y,N,N,N)
OUTPUT(EPSX,Y,N,Y,N,N,N)
OUTPUT(TAXY,Y,N,Y,N,N,N)
OUTPUT(STRY,Y,N,Y,N,N,N)
OUTPUT(STRX,Y,N,Y,N,N,N)
OUTPUT(DIL,Y,N,Y,N,N,N)
OUTPUT(CON4,Y,N,Y,N,N,N)
OUTPUT(CON3,Y,N,Y,N,N,N)
OUTPUT(CON2,Y,N,Y,N,N,N)
OUTPUT(CON1,Y,N,Y,N,N,N)
OUTPUT(POIS,Y,N,Y,N,N,N)
OUTPUT(STIF,Y,N,Y,N,N,N)
OUTPUT(THEX,Y,N,Y,N,N,N)
OUTPUT(TEM,Y,N,Y,Y,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 9 ;IZMON = 1
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NXPRIN = 1 ;IXPRF = 1 ;IXPRL = 10000
NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 10000
IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
No PATCHes yet used for this Group
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP