feap ** Wave propagation in elastic bar ** 0,0,1,2,2,4 param h=0.5 l=10 g=0.6 ! 0.5 No numerical dissipation, 0.6 gives some dissipation block bar cart 200 1 1,0,0 2,l,0 3,l,h 4,0,h ebou 1,0,1,0 1,l,1,0 2,0,0,1 2,h,0,1 boun 1,,1,1 edisp 1,l,1,0 mate,1 solid plane strain elastic,isotropic,200.d9,0.3 density body 7850.0 end opti batch ! Load Table prop end 2 0.0 0.0 .5d-4 1.d-2 1.0d-4 0.0 batch ! Setup Printing Turn on newmark plot,mesh transient,newmark,g/2,g plot,pers,1 plot,hide plot,dofs,0,0,1 plot,defo,50,1 end 0 40 -10 5 0 0 1 batch ! Turn off printing to prevent slow down for output nopr print,off end batch dt,,1.0d-5 time tang,,1 plot,defa,off ! Turn off defaults for ranges plot,cont,1,,1 ! contour axial displacment and set ranges plot,prompt,off ! Turn off plot prompting loop,,100 ! 100 sets of five time steps loop,,5 ! Five time steps before replot time form solv next plot,wipe ! Clear plot and redraw plot,cont,1,,-1 ! Last minus sign suppresses mesh plus uses previous range jpgd next end -5.d-3 5.d-3 y inte stop