TALK=T;RUN(1,1) DISPLAY Experiment Q1 ==> VR 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 U1 z 1 fil;.0001 vec z 1 co 1 dump S204_U pause cl msg y-displacement fields gr ou z 1 cont V1 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; s204 ) ************************************************************ ************************************************************ 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(1)=P1 ;NAME(3)=U1 NAME(5)=V1 ;NAME(135)=KOND NAME(136)=SPH1 ;NAME(137)=ENUL NAME(138)=DEN1 ;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(P1,Y,Y,Y,N,N,N) SOLUTN(U1,Y,Y,Y,N,N,Y) SOLUTN(V1,Y,Y,Y,N,N,Y) SOLUTN(KOND,Y,N,N,N,N,Y) SOLUTN(SPH1,Y,N,N,N,N,Y) SOLUTN(ENUL,Y,N,N,N,N,Y) SOLUTN(DEN1,Y,N,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) DEN1 = 138 VISL = 137 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 =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(U1)=1. ;PRNDTL(V1)=1. PRT(U1)=1. ;PRT(V1)=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(P1)=0. ;FIINIT(U1)=0. FIINIT(V1)=0. ;FIINIT(KOND)=1.0E-10 FIINIT(SPH1)=1.0E-10 ;FIINIT(ENUL)=1.0E-10 FIINIT(DEN1)=1.0E-10 ;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(INI1 ,INIVAL, 1, 16, 1, 8, 1, 1, 1, 1) INIT(INI1 ,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 ,V1 ,1. ,0. ) PATCH(FORS01 ,SOUTH , 1, 16, 1, 1, 1, 1, 1, 1) COVAL(FORS01 ,V1 , FIXFLU ,-4.0E+07 ) PATCH(AXESZZ ,WWALL , 1, 1, 1, 16, 1, 1, 1, 1) COVAL(AXESZZ ,U1 ,1. ,0. ) XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 100 ;ISWC1 = 1 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 SELREF = T RESFAC =1.0E-07 ************************************************************ Group 16. Terminate Iterations LITER(P1)=50 ;LITER(U1)=50 LITER(V1)=50 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03 ENDIT(V1)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,1.) RELAX(U1,FALSDT,1.) RELAX(V1,FALSDT,1.) RELAX(KOND,LINRLX,1.) RELAX(SPH1,LINRLX,1.) RELAX(ENUL,LINRLX,1.) RELAX(DEN1,LINRLX,1.) 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(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(KOND)=1.0E+10 ;VARMIN(KOND)=-1.0E+10 VARMAX(SPH1)=1.0E+10 ;VARMIN(SPH1)=-1.0E+10 VARMAX(ENUL)=1.0E+10 ;VARMIN(ENUL)=-1.0E+10 VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-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 = 100 ISG52 = 3 ISG62 = 1 CSG10 ='q1' SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,160) SPEDAT(SET,DOMAIN,PROPS_FILE_1,C,Q1) SPEDAT(SET,MATERIAL,162,L,T) SPEDAT(SET,BOUNDARY,ZCONST,R,0.) SPEDAT(SET,STORED,U1T,C,=-6.000001E-05*XU!ZSLFIN) SPEDAT(SET,STORED,V1T,C,=2.0E-04*(YV-0.12)!ZSLFIN) SPEDAT(SET,STORED,U1/T,C,=U1/(U1T+1.E-20)!ZSLFIN) SPEDAT(SET,STORED,V1/T,C,=V1/(V1T+1.E-20)!ZSLFIN) SPEDAT(SET,INITIAL,PRPS!2NDMAT,C,=162) 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,OBJNAM,!INI1,C,2NDMAT) SPEDAT(SET,OBJTYP,!INI1,C,USER_DEFINED) 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(P1,Y,N,Y,Y,Y,Y) OUTPUT(U1,Y,N,Y,Y,Y,Y) OUTPUT(V1,Y,N,Y,Y,Y,Y) OUTPUT(KOND,Y,N,Y,N,N,N) OUTPUT(SPH1,Y,N,Y,N,N,N) OUTPUT(ENUL,Y,N,Y,N,N,N) OUTPUT(DEN1,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(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