TALK=F;RUN(1,1)                                           
  DISPLAY
  A rectangular beam of dimensions Lx 
  for Time = 0 spraining force FX. For Time > 0 force FX=0 
  and beam makes the longitudal fluctuations.
  ENDDIS

  PHOTON USE 
  p;parphi 
  1 1 2000
  
  clear
  set prop off
  view y
  gr ou y 1
  con disx y 1 fi;0.0001
  msg DISX(X,Z=Time)
  pause
  
  
  autoplot 
  FILE 
  parphi 5 

  clear
  data 1                                                                          
  DISX Z 10 X 1 m                                                                  
  col1 1          
  data 1                                                                          
  DISX Z 50 X 1 m                                                                  
  colA 2          
  data 1                                                                          
  DISX Z 100 X 1 m                                                                  
  colF 3          
  scale
  plot 1 2 3                                                                      
  msg DISX(X) for Time=10(W),50(Y),100(R)
  pause

  clear
  data 1                                                                          
  DISX X 80 Z 1 m                                                                  
  col1 1          
  scale
  plot 1                                                                      
  msg DISX(Time) for X=LX
  pause


 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
TEXT(Wave in 1DX free beam; S670
libref=670
TITLE
 
  Declarations and settings
REAL(FX,LYZ,LX,POISSON,YOUNG) 
FX= 40.0e6  ! H/m^2 = 40 N/mm^2
LYZ=40.e-3
LX=120.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(NXBODY)

 ************************************************************
  Group 2. Time dependence
 STEADY  =    F
 RSET(U,0,2.e-4,800)

 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES  =    T
 NY = 1
 YVLAST  = LYZ

 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NXBODY = 80
 GRDPWR(X,NXBODY,LX,1)  
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 NZ=1
 ZWLAST  = LYZ
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
 SOLVE(DISX)
 SOLUTN(DISX  ,Y,Y,Y,N,N,Y)

 TERMS (DISX,N,N,Y,Y,Y,N)

 STORE(VDSX)
 
 STORE(PRPS,DEN1,ENUL,DVO1,DRH1)
 STORE(STRY,VTHP,VTHM)
 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
 real(CEPSX,CEPSY,CEPS)
 CEPSX = FX/YOUNG
 
   **** CALCULATE analytical solution ***
 PATCH(BODY,INIVAL,1,NX,1,1,1,1,1,1)
 (INITIAL of disX at BODY is :CEPSX:*XG)

FIINIT(VDSX)=0.0
FIINIT(PRPS)=160

RELAX(DISX,LINRLX, 1.000000E+00)



 ************************************************************
  GROUP 13. BOUNDARY & SPECIAL SOURCES
 
PATCH(Left,WWALL,1,1,1,1,1,1,1,LSTEP)    ! left edge fixed 
COVAL(Left,DISX,1,0.0)

 ************************************************************
  GROUP 15. TERMINATE SWEEPS
 LSWEEP  =      10
 ISG21=LSWEEP
  
      RESREF(disx)=0.0  ! to prevent premature exit
      LITER(disx) = 40  ! from solver

 
 ************************************************************
  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
IXMON = NX
IYMON = 1
IZMON = 1

IDISPA = 1


  inform7begin
(STORE VAR VTHP is :CEPSX:*XG)
(STORE VAR VTHM is -:CEPSX:*XG)
 
  inform7end
  
 ************************************************************
  
  
STOP