TALK=T;RUN(1,1)
DISPLAY
A rectangular plate of dimensions LX * LY is
* held fixed along its north and west edges
* subjected to uniform tension along its south edge
* allowed to contract freely along its east edge
* in case 1 (plane stress) allowed to contract freely in the
z-direction
* in case 2 (plane strain) prevented from contracting in the
z-direction.
The distributions of displacement in the xy plane are computed,
and also the corresponding stresses and strains in all three
directions.
A uniform cartesian grid is used.
The computed results are compared with the analytically-derived
values.
Cases 3 and 4 represent the same situation but with then lower-IY
part of the plate being composed of a different material from the
uper part. This material has 10 times the compressibility of the
first material.
In case 3, the poisson's ratio is one tenth that of the first
material, so that the lateral extension remains the same, and
the solution remains one-dimensional.
In case 4 the poisson's ratio is the same as that of the first
material; so the lateral extensions differ, and the solution
becomes two-dimensional.
One-dimensionality can however be restored by making LX much
smaller than LY, for example by moving two spaces to the left
(and so activating) the line LX=LY/100 below.
ENDDIS
PHOTON USE
p;;;;
set prop off
cl
msg x-displacement fields
gr ou z 1
cont DISX z 1 fil;.0001
SET VEC COMP - - -
SET VEC COMP DISX DISY -
VECT IZ 1 SH DASH 0
dump
S204_U
pause
cl
msg y-displacement fields
gr ou z 1
cont DISY z 1 fil;.0001
vec z 1 co 1
dump
S204_V
pause
cl
msg x-direction strain fields
msg of interest only for casenos 3 and 4
gr ou z 1
cont epsx z 1 fil;.0001
pause
cl
msg y-direction strain fields
msg of interest only for casenos 3 and 4
gr ou z 1
cont epsy z 1 fil;.0001
pause
cl
msg shear-stress fields
msg of interest only for casenos 3 and 4
gr ou z 1
cont stxy z 1 fil;.0001
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(2D xy 2-material plate in tension; S604 )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 16
XULAST =0.09
XFRAC(1)=0.0625 ;XFRAC(3)=0.1875
XFRAC(5)=0.3125 ;XFRAC(7)=0.4375
XFRAC(9)=0.5625 ;XFRAC(11)=0.6875
XFRAC(13)=0.8125 ;XFRAC(15)=0.9375
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 16
YVLAST =0.12
YFRAC(1)=0.0625 ;YFRAC(3)=0.1875
YFRAC(5)=0.3125 ;YFRAC(7)=0.4375
YFRAC(9)=0.5625 ;YFRAC(11)=0.6875
YFRAC(13)=0.8125 ;YFRAC(15)=0.9375
************************************************************
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(137)=DISY ;NAME(138)=DISX
NAME(139)=V1/T ;NAME(140)=U1/T
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(DISY,Y,Y,N,N,N,Y)
SOLUTN(DISX,Y,Y,N,N,N,Y)
SOLUTN(V1/T,Y,N,N,N,N,Y)
SOLUTN(U1/T,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)
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(DISY,N,Y,Y,Y,N,Y)
TERMS(DISX,N,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 =7800. ;TMP1 =0. ;EL1 =0.
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =1.0E-05 ;DRH1DP =5.0E-12
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =0.3 ;ENUT =0.
PRNDTL(DISY)=1. ;PRNDTL(DISX)=1.
PRT(DISY)=1. ;PRT(DISX)=1.
CP1 =473. ;CP2 =1.
* List of user-defined materials to be read by EARTH
MATFLG=T;IMAT=3
* 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.03 473.0 43.0 1.0E-5 0.5E-10
*
162 7800.0 0.3 473.0 43.0 1.0E-5 0.5E-10
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(DISY)=0. ;FIINIT(DISX)=0.
FIINIT(V1/T)=1.001E-10 ;FIINIT(U1/T)=1.001E-10
FIINIT(V1T)=1.001E-10 ;FIINIT(U1T)=1.001E-10
FIINIT(EPSZ)=1.001E-10 ;FIINIT(EPSX)=1.001E-10
FIINIT(EPSY)=1.001E-10 ;FIINIT(STXY)=1.001E-10
FIINIT(STRZ)=1.001E-10 ;FIINIT(STRY)=1.001E-10
FIINIT(STRX)=1.001E-10 ;FIINIT(PRPS)=160.
PATCH(2NDMAT ,INIVAL, 1, 16, 1, 8, 1, 1, 1, 1)
INIT(2NDMAT ,PRPS,In-Form:initial value - G19)
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(UP ,NWALL , 1, 16, 16, 16, 1, 1, 1, 1)
COVAL(UP ,DISY,1. ,0. )
PATCH(FORS01 ,SOUTH , 1, 16, 1, 1, 1, 1, 1, 1)
COVAL(FORS01 ,DISY,In-Form:source - see Grp 19)
PATCH(AXESZZ ,WWALL , 1, 1, 1, 16, 1, 1, 1, 1)
COVAL(AXESZZ ,DISX,1. ,0. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 400 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-07
************************************************************
Group 16. Terminate Iterations
LITER(DISY)=20 ;LITER(DISX)=20
ENDIT(DISY)=1.0E-03 ;ENDIT(DISX)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(DISY,FALSDT,1.0E+09)
RELAX(DISX,FALSDT,1.0E+09)
RELAX(V1/T,LINRLX,1.)
RELAX(U1/T,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(DISY)=1.0E+10 ;VARMIN(DISY)=-1.0E+10
VARMAX(DISX)=1.0E+10 ;VARMIN(DISX)=-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(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 = 400
ISG50 = 1
ISG52 = 1
ISG62 = 1
CSG10 ='q1'
SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,160)
SPEDAT(SET,DOMAIN,PROPS_FILE_1,C,Q1)
SPEDAT(SET,BOUNDARY,ZCONST,R,1.0E+20)
SPEDAT(SET,STORED,U1T,C,=-7.8E-05*XG!ZSLFIN)
SPEDAT(SET,STORED,V1T,C,=1.82E-04*(YG-0.12)!ZSLFIN)
SPEDAT(SET,STORED,U1/T,C,=DISX/(U1T+1.E-20)!ZSLFIN)
SPEDAT(SET,STORED,V1/T,C,=DISY/(V1T+1.E-20)!ZSLFIN)
SPEDAT(SET,INITIAL,PRPS!2NDMAT,C,=160)
SPEDAT(SET,SOURCE,DISY!FORS01,C,=COVAL(2.0E-10&1.0/ANORTH*(-4.0E+$)
SPEDAT(SET,SOURCE,DISY!FORS01,C,07)))
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,OBJNAM,!2NDMAT,C,2NDMAT)
SPEDAT(SET,OBJTYP,!2NDMAT,C,USER_DEFINED)
SPEDAT(SET,OBJNAM,!UP,C,UP)
SPEDAT(SET,OBJTYP,!UP,C,USER_DEFINED)
SPEDAT(SET,OBJNAM,!FORS01,C,FORS01)
SPEDAT(SET,OBJTYP,!FORS01,C,USER_DEFINED)
SPEDAT(SET,OBJNAM,!AXESZZ,C,AXESZZ)
SPEDAT(SET,OBJTYP,!AXESZZ,C,USER_DEFINED)
SPEDAT(SET,FACETDAT,NUMOBJ,I,4)
SPEDAT(SET,MATERIAL,160,L,T)
SPEDAT(SET,MATERIAL,0,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(DISY,Y,N,Y,Y,Y,Y)
OUTPUT(DISX,Y,N,Y,Y,Y,Y)
OUTPUT(V1/T,Y,N,Y,N,N,N)
OUTPUT(U1/T,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 = 14 ;IYMON = 2 ;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 = 2 ;IXPRF = 1 ;IXPRL = 10000
NYPRIN = 2 ;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