feap ** Bar test with Lysmer-Kuhlemeyer dampers **
0 0 0 2 2 4


parameters
L = 10
T = L/4
E = 10
nu = 0.3
p = 4*atan(1)
d = 1
mu = E/2/(1+nu)
la = E*nu/(1+nu)/(1-2*nu)
n = 20
m = 1
f = 0
s = f*L/n


snode
 1 0 0
 2 L 0
 3 L+s T
 4 0 T


blend
  surface n m 
  1 2 3 4

blend
  line m 0 0 0 2
  2 3

edisp
  1 0 -1e-2 0

eboun
  1 0 1 0
  2 0 0 1
  2 T 0 1
  



mate 1
solid
 elastic isotropic E nu
 density material d



mate 2
user 40
 Cp d*sqrt(mu/d+la/d)
 Cs d*sqrt(mu/d)


end

tie



batch
  prop,,1
end
1 1 0 1e5 0 0 1 2*p/.6 0        ! Sine drive


! Step load
! 2 1
! 0 1
! 0.06 1
! 0.09 1
! 1e6  1


batch
 tplo
end
disp 10 1
show


batch
  dt,,0.03/3
  tran,explict
  mass,lump
  form accel
  loop,,600*4
    time        ! explicit
    form explicit
  next
  plot,stre,1
end


inte

stop