TALK=T;RUN(1,1)
DISPLAY
This is a 2D (xy) Plane-Strain problem, in which a beam is bent
by a torqe applied at its right-hand end.
The problem has an exact analytical solution which is used to
enable the accuracy of the computed values to be assessed by
being divided by the exact values, and printed as U1/T, etc ,
which, if equal to 1.0 signify correctness. [for casenos 1-11
only at present.]
A menu is provided which enables:
* the Poisson's ratio to be varied
* the z-dimension to be varied
* the beam to be made much more slender
* the end torque to be replaced by a tension force
* point-wise and distributed forces to be applied
* temperature gradients to be applied in x and y directions
* the grid-fineness and distribution to be varied.
In order to ensure accuracy, the grid distribution in beam-
thickness is such as to provide one very thin cell near each beam
surface.
ENDDIS
PHOTON USE
p;;;;
set prop off
cl
msg x-direction-displacement field
gr ou z 1
cont U1 z 1 fil;.0001
vec z 1 y 2 22 col 0
pause
cl
msg y-direction-displacement field
gr ou z 1
cont V1 z 1 fil;.0001
vec z 1 col 0
pause;cl
msg v1/v1_exact field
gr ou z 1
cont V1/t z 1 fi;.0001
pause;cl
msg u1/u1_exact field
gr ou z 1
cont u1/t z 1 fil;.0001
pause;cl
msg p1/p1_exact field
gr ou z 1
cont p1/t z 1 fil;.0001
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(bent beam in plane strain; s220 )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 3
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 22
XULAST =0.1224
XFRAC(1)=9.803922E-03 ;XFRAC(2)=0.058824
XFRAC(3)=0.107843 ;XFRAC(4)=0.156863
XFRAC(5)=0.205882 ;XFRAC(6)=0.254902
XFRAC(7)=0.303922 ;XFRAC(8)=0.352941
XFRAC(9)=0.401961 ;XFRAC(10)=0.45098
XFRAC(11)=0.5 ;XFRAC(12)=0.54902
XFRAC(13)=0.598039 ;XFRAC(14)=0.647059
XFRAC(15)=0.696078 ;XFRAC(16)=0.745098
XFRAC(17)=0.794118 ;XFRAC(18)=0.843137
XFRAC(19)=0.892157 ;XFRAC(20)=0.941176
XFRAC(21)=0.990196 ;XFRAC(22)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 22
YVLAST =0.0408
YFRAC(1)=9.803922E-03 ;YFRAC(2)=0.058824
YFRAC(3)=0.107843 ;YFRAC(4)=0.156863
YFRAC(5)=0.205882 ;YFRAC(6)=0.254902
YFRAC(7)=0.303922 ;YFRAC(8)=0.352941
YFRAC(9)=0.401961 ;YFRAC(10)=0.45098
YFRAC(11)=0.5 ;YFRAC(12)=0.54902
YFRAC(13)=0.598039 ;YFRAC(14)=0.647059
YFRAC(15)=0.696078 ;YFRAC(16)=0.745098
YFRAC(17)=0.794118 ;YFRAC(18)=0.843137
YFRAC(19)=0.892157 ;YFRAC(20)=0.941176
YFRAC(21)=0.990196 ;YFRAC(22)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 1
ZWLAST =1.0E-03
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)=P1/T
NAME(135)=V1/T ;NAME(136)=U1/T
NAME(137)=DRH1 ;NAME(138)=DVO1
NAME(139)=VISL ;NAME(140)=DILT
NAME(141)=V1T ;NAME(142)=U1T
NAME(143)=EPSZ ;NAME(144)=EPSX
NAME(145)=EPSY ;NAME(146)=STXY
NAME(147)=STRZ ;NAME(148)=STRY
NAME(149)=STRX ;NAME(150)=PRPS
* 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(P1/T,Y,N,N,N,N,N)
SOLUTN(V1/T,Y,N,N,N,N,N)
SOLUTN(U1/T,Y,N,N,N,N,N)
SOLUTN(DRH1,Y,N,N,N,N,Y)
SOLUTN(DVO1,Y,N,N,N,N,Y)
SOLUTN(VISL,Y,N,N,N,N,Y)
SOLUTN(DILT,Y,N,N,N,N,Y)
SOLUTN(V1T,Y,N,N,N,N,Y)
SOLUTN(U1T,Y,N,N,N,N,Y)
SOLUTN(EPSZ,Y,N,N,N,N,Y)
SOLUTN(EPSX,Y,N,N,N,N,Y)
SOLUTN(EPSY,Y,N,N,N,N,Y)
SOLUTN(STXY,Y,N,N,N,N,Y)
SOLUTN(STRZ,Y,N,N,N,N,Y)
SOLUTN(STRY,Y,N,N,N,N,Y)
SOLUTN(STRX,Y,N,N,N,N,Y)
SOLUTN(PRPS,Y,N,N,N,N,Y)
VISL = 139
PRPS = 150
************************************************************
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,Y,Y,Y,Y,Y)
TERMS(V1,Y,Y,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.0E-05 ;ENUT =0.
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRT(U1)=1. ;PRT(V1)=1.
CP1 =1. ;CP2 =1.
* List of user-defined materials to be read by EARTH
MATFLG=T;IMAT=2
* Name
*Ind. Dens. Viscos. Spec.heat Conduct. Expans. Compr.
*
160 7800.0 0.3 473.0 43.0 1.0E-5 0.5E-11
*
161 7800.0 0.0 473.0 43.0 1.0E-5 0.5E-11
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=0. ;FIINIT(U1)=0.
FIINIT(V1)=0. ;FIINIT(P1/T)=1.0E-10
FIINIT(V1/T)=1.0E-10 ;FIINIT(U1/T)=1.0E-10
FIINIT(DRH1)=1.0E-10 ;FIINIT(DVO1)=1.0E-10
FIINIT(VISL)=1.0E-10 ;FIINIT(DILT)=0.
FIINIT(V1T)=0. ;FIINIT(U1T)=0.
FIINIT(EPSZ)=1.0E-10 ;FIINIT(EPSX)=1.0E-10
FIINIT(EPSY)=1.0E-10 ;FIINIT(STXY)=1.0E-10
FIINIT(STRZ)=1.0E-10 ;FIINIT(STRY)=1.0E-10
FIINIT(STRX)=1.0E-10 ;FIINIT(PRPS)=0.
PATCH(BODY ,INIVAL, 2, 21, 2, 21, 1, 1, 1, 1)
INIT(BODY ,PRPS, FIXVAL ,160. )
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(AXESUU ,WEST , 1, 1, 2, 21, 1, 1, 1, 1)
COVAL(AXESUU ,U1 ,1.0E+06 ,0. )
PATCH(FORC01 ,EAST , 21, 21, 2, 21, 1, 1, 1, 1)
PATCH(RIGHTEND,CELL , 21, 21, 2, 21, 1, 1, 1, 1)
COVAL(RIGHTEND,U1 ,In-Form:source - see Grp 19)
PATCH(RIGHTENV,CELL , 21, 21, 1, 21, 1, 1, 1, 1)
COVAL(RIGHTENV,V1 ,In-Form:source - see Grp 19)
PATCH(AXESVV ,CELL , 2, 2, 1, 21, 1, 1, 1, 1)
COVAL(AXESVV ,V1 ,In-Form:source - see Grp 19)
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1600 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-08
************************************************************
Group 16. Terminate Iterations
LITER(P1)=4 ;LITER(U1)=10
LITER(V1)=100
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=0.
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(U1,FALSDT,1.)
RELAX(V1,FALSDT,1.)
RELAX(P1/T,LINRLX,1.)
RELAX(V1/T,LINRLX,1.)
RELAX(U1/T,LINRLX,1.)
RELAX(DRH1,LINRLX,1.)
RELAX(DVO1,LINRLX,1.)
RELAX(VISL,LINRLX,1.)
RELAX(DILT,LINRLX,1.)
RELAX(V1T,LINRLX,1.)
RELAX(U1T,LINRLX,1.)
RELAX(EPSZ,LINRLX,1.)
RELAX(EPSX,LINRLX,1.)
RELAX(EPSY,LINRLX,1.)
RELAX(STXY,LINRLX,1.)
RELAX(STRZ,LINRLX,1.)
RELAX(STRY,LINRLX,1.)
RELAX(STRX,LINRLX,1.)
RELAX(PRPS,LINRLX,1.)
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(P1/T)=1.0E+10 ;VARMIN(P1/T)=-1.0E+10
VARMAX(V1/T)=1.0E+10 ;VARMIN(V1/T)=-1.0E+10
VARMAX(U1/T)=1.0E+10 ;VARMIN(U1/T)=-1.0E+10
VARMAX(DRH1)=1.0E+10 ;VARMIN(DRH1)=-1.0E+10
VARMAX(DVO1)=1.0E+10 ;VARMIN(DVO1)=-1.0E+10
VARMAX(VISL)=1.0E+10 ;VARMIN(VISL)=-1.0E+10
VARMAX(DILT)=1.0E+10 ;VARMIN(DILT)=-1.0E+10
VARMAX(V1T)=1.0E+10 ;VARMIN(V1T)=-1.0E+10
VARMAX(U1T)=1.0E+10 ;VARMIN(U1T)=-1.0E+10
VARMAX(EPSZ)=1.0E+10 ;VARMIN(EPSZ)=-1.0E+10
VARMAX(EPSX)=1.0E+10 ;VARMIN(EPSX)=-1.0E+10
VARMAX(EPSY)=1.0E+10 ;VARMIN(EPSY)=-1.0E+10
VARMAX(STXY)=1.0E+10 ;VARMIN(STXY)=-1.0E+10
VARMAX(STRZ)=1.0E+10 ;VARMIN(STRZ)=-1.0E+10
VARMAX(STRY)=1.0E+10 ;VARMIN(STRY)=-1.0E+10
VARMAX(STRX)=1.0E+10 ;VARMIN(STRX)=-1.0E+10
VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
STRA = T
PARSOL = F
ISG21 = 1600
ISG50 = 1
ISG52 = 1
ISG62 = 1
CSG10 ='Q1'
SPEDAT(SET,SOURCE,U1!RIGHTEND,C,=COVAL(1.E10&U1T))
SPEDAT(SET,SOURCE,V1!RIGHTENV,C,=COVAL(1.0E10&V1T))
SPEDAT(SET,SOURCE,V1!AXESVV,C,=COVAL(1.E10&V1T))
SPEDAT(SET,BOUNDARY,ZCONST,R,1.0E+20)
SPEDAT(SET,LONGNAME,P1,C,pressure_or_dilatation)
SPEDAT(SET,LONGNAME,PTH,C,dilatation_theoretical)
SPEDAT(SET,LONGNAME,ENUL,C,viscosity_or__LAME's_G*1.E11)
SPEDAT(SET,LONGNAME,VISL,C,viscosity_or__LAME's_G*1.E11)
SPEDAT(SET,LONGNAME,DRH1,C,compressibility_or_LAME's_L+G*1.E11)
SPEDAT(SET,LONGNAME,DVO1,C,th.exp.co_r_or_H*1.E11)
SPEDAT(SET,LONGNAME,U1,C,x-direction_velocity_or_displacement)
SPEDAT(SET,LONGNAME,EPSX,C,x-direction_strain)
SPEDAT(SET,LONGNAME,EXTH,C,x-direction_strain_theoretical)
SPEDAT(SET,LONGNAME,STRX,C,x-direction_stress)
SPEDAT(SET,LONGNAME,SXTH,C,x-direction_stress_theoretical)
SPEDAT(SET,LONGNAME,V1,C,y-direction_velocity_or_displacement)
SPEDAT(SET,LONGNAME,EPSY,C,y-direction_strain)
SPEDAT(SET,LONGNAME,EYTH,C,y-direction_strain_theoretical)
SPEDAT(SET,LONGNAME,STRY,C,y-direction_stress)
SPEDAT(SET,LONGNAME,SYTH,C,y-direction_stress_theoretical)
SPEDAT(SET,LONGNAME,W1,C,z-direction_velocity_or_displacement)
SPEDAT(SET,LONGNAME,EPSZ,C,z-direction_strain)
SPEDAT(SET,LONGNAME,STRZ,C,z-direction_stress)
SPEDAT(SET,LONGNAME,EZTH,C,z-direction_strain_theoretical)
SPEDAT(SET,LONGNAME,SZTH,C,z-direction_stress_theoretical)
SPEDAT(SET,LONGNAME,EPST,C,linearthermalexpansion)
SPEDAT(SET,STORED,V1T,C,=-4.55E-03*((YV-0.51*0.04)^2*0.428571+(XG$)
SPEDAT(SET,STORED,V1T,C,-0.01*0.12)^2)!ZSLFIN)
SPEDAT(SET,STORED,DILT,C,=5.2E-03*(YG-0.51*0.04)!IMAT>100!ZSLFIN)
SPEDAT(SET,STORED,U1T,C,=9.099999E-03*(XU-0.01*0.12)*(YG-0.51*0.0$)
SPEDAT(SET,STORED,U1T,C,4)!IMAT>100!ZSLFIN)
SPEDAT(SET,STORED,U1/T,C,=U1/U1T!IMAT>100!ZSLFIN)
SPEDAT(SET,STORED,V1/T,C,=V1/V1T!ZSLFIN)
SPEDAT(SET,STORED,P1/T,C,=P1/DILT!IMAT>100!ZSLFIN)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,MATERIAL,0,L,T)
SPEDAT(SET,MATERIAL,160,L,T)
************************************************************
Group 20. Preliminary Printout
************************************************************
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(P1/T,Y,N,Y,N,N,N)
OUTPUT(V1/T,Y,N,Y,N,N,N)
OUTPUT(U1/T,Y,N,Y,N,N,N)
OUTPUT(DRH1,N,N,N,N,N,N)
OUTPUT(DVO1,Y,N,Y,N,N,N)
OUTPUT(VISL,N,N,N,N,N,N)
OUTPUT(DILT,Y,N,Y,N,N,N)
OUTPUT(V1T,Y,N,Y,N,N,N)
OUTPUT(U1T,Y,N,Y,N,N,N)
OUTPUT(EPSZ,Y,N,Y,N,N,N)
OUTPUT(EPSX,Y,N,Y,N,N,N)
OUTPUT(EPSY,Y,N,Y,N,N,N)
OUTPUT(STXY,Y,N,Y,N,N,N)
OUTPUT(STRZ,Y,N,Y,N,N,N)
OUTPUT(STRY,Y,N,Y,N,N,N)
OUTPUT(STRX,Y,N,Y,N,N,N)
OUTPUT(PRPS,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 20 ;IYMON = 11 ;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
PATCH(VPROF ,PROFIL, 1, 22, 1, 1, 1, 1, 1, 1)
PLOT(VPROF ,V1 ,0. ,0. )
PLOT(VPROF ,V1T ,0. ,0. )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP