TALK=F;RUN(1,1)
DISPLAY
2D (x,y) Y-shape beam is free for Time <= 0.
:----- fixed (DisX=DisY=0)
:
: **************************
V * : Top :
/***** : :
/* ***********************
/* *
/* ***********************
/***** : :
* : Bottom :
**************************
Initial condition Time = 0:
caseno = 1: in Top VDSY = 0.1
caseno = 2: in Top and Bottom VDSY = 0.1
caseno = 3: in Top VDSY = 0.1, in Bottom VDSY = -0.1
ENDDIS
PHOTON USE
p;parphi
1 1 200
clear
set prop off
view y
gr ou y 30
con disy y 30 fi;0.0001
msg DISY (X,Z=Time) on axis of top beam
pause
clear
view 1 1 1
gr ou z 10 x 1 12 y 13 24
gr ou z 10 x 9 m y 1 12
gr ou z 10 x 9 m y 25 m
SET VEC COMP - - -
SET VEC COMP DISX DISY -
VECT IZ 10 SH DASH 0
gr ou z 100 x 1 12 y 13 24
gr ou z 100 x 9 m y 1 12
gr ou z 100 x 9 m y 25 m
SET VEC COMP - - -
SET VEC COMP DISX DISY -
VECT IZ 100 SH DASH 0
gr ou z 190 x 1 12 y 13 24
gr ou z 190 x 9 m y 1 12
gr ou z 190 x 9 m y 25 m
SET VEC COMP - - -
SET VEC COMP DISX DISY -
VECT IZ 190 SH DASH 0
msg Displacement vector (Z=Time=10, 100, 190)
pause
autoplot
FILE
parphi 5
clear
data 1
DISY X m Y m Z 1 m
col1 1
data 1
DISY X m Y 1 Z 1 m
colA 2
scale
plot 1 2
msg DISY(Time) on end of top(W) and bottom(Y) beam
pause
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
TEXT(Free wave in 2D "tuning-fork" model; S673
libref=673
TITLE
integer(caseno)
mesg(caseno 1 : Positive initial Velocity in TOP beam
mesg(caseno 2 : Positive initial Velocity in BOTH beam
mesg(caseno 3 : Symmetrical initial conditions
caseno=1
mesgm(caseno = :caseno: Enter another if not OK
readvdu(caseno, int, 1)
caseno
Declarations and settings
REAL(LZ,LX,LY,POISSON,YOUNG)
LX = 120.e-3
LY = 36.e-3
LZ = 10.e-3
YOUNG = 1/0.5E-11 ! Young's modulus used for exact solutions
POISSON=0.3 ! Poisson's ratio used for exact solutions
! The above must correspond to the data
! below MATFLG
INTEGER(NYBODY,NXBODY)
************************************************************
Group 2. Time dependence
STEADY = F
RSET(U,0,30.e-4,200)
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NXBODY = 32
GRDPWR(X,NXBODY,LX,1)
************************************************************
Group 4. Y-Direction Grid Spacing
NYBODY = 36
GRDPWR(Y,NYBODY,LY,1)
************************************************************
Group 5. Z-Direction Grid Spacing
NZ=1
ZWLAST = LZ
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
SOLVE(DISX,DISY)
SOLUTN(DISY ,Y,Y,Y,N,N,Y)
SOLUTN(DISX ,Y,Y,Y,N,N,Y)
TERMS (DISY,N,N,Y,Y,Y,N)
TERMS (DISX,N,N,Y,Y,Y,N)
STORE(VDSY,VDSX)
STORE(PRPS,DEN1,ENUL,DVO1,DRH1)
STORE(STRX,STRY)
STORE(EPSY,EPSX,EPSZ)
************************************************************
GROUP 8. ITERATION NUMBERS ETC
************************************************************
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
************************************************************
GROUP 11. INITIAL VALUES
**** CALCULATE initial conditions ***
FIINIT(disx)=0
FIINIT(disy)=0
FIINIT(VDSX)=0.0
FIINIT(VDSY)=0.0
FIINIT(PRPS)=160
PATCH(FLUID1,INIVAL,1,8,1,12,1,1,1,1)
INIT(FLUID1,PRPS,FIXVAL,0)
PATCH(FLUID2,INIVAL,1,8,ny-11,ny,1,1,1,1)
INIT(FLUID2,PRPS,FIXVAL,0)
PATCH(FLUID3,INIVAL,13,nx,13,ny-12,1,1,1,1)
INIT(FLUID3,PRPS,FIXVAL,0)
**** Initial Velocity ---
PATCH(INITVEL,INIVAL,26,nx,NY-11,NY,1,1,1,1)
INIT(INITVEL,VDSY,FIXVAL,0.1)
if(caseno.eq.2) then
PATCH(INITVEL1,INIVAL,26,nx,1,12,1,1,1,1)
INIT(INITVEL1,VDSY,FIXVAL, 0.1)
endif
if(caseno.eq.3) then
PATCH(INITVEL1,INIVAL,26,nx,1,12,1,1,1,1)
INIT(INITVEL1,VDSY,FIXVAL, -0.1)
endif
RELAX(DISY,LINRLX, 1.000000E+00)
RELAX(DISX,LINRLX, 1.000000E+00)
************************************************************
GROUP 13. BOUNDARY & SPECIAL SOURCES
PATCH(Left,WWALL,1,1,13,24,1,1,1,LSTEP) ! west edge fixed
COVAL(Left,DISY,1,0.0)
COVAL(Left,DISX,1,0.0)
************************************************************
GROUP 15. TERMINATE SWEEPS
LSWEEP = 120
ISG21=LSWEEP
************************************************************
GROUP 17. RELAXATION
RELAX(P1 ,LINRLX, 1.000000E+00)
************************************************************
GROUP 19. DATA TRANSMITTED TO GROUND
STRA = T
PARSOL = F
ISG52 = 3 ! probe & res
************************************************************
GROUP 23.FIELD PRINT-OUT & PLOT CONTROL
TSTSWP = - 1 ! graphic-mode
NYPRIN = 1
NXPRIN = 1
IXMON = NX
IYMON = NY
IZMON = 1
IDISPA = 1
inform7begin
inform7end
************************************************************
STOP