feap ** Resonator test with Lysmer-Kuhlemeyer dampers in axis-symm **
0 0 0 2 2 4


parameters
R = 10
T = R/3
E = 10
nu = 0.3
n = 20
p = 4*atan(1)
d = 1
mu = E/2/(1+nu)
la = E*nu/(1+nu)/(1-2*nu)


snode
1 0 0
2 0 -T 
3 T -T 
4 T 0
5 0 -R
6 R*cos(-p/4) R*sin(-p/4)
7 R 0

side
 polar 5 6 1
 polar 6 7 1

blend
  surface n n 
  1 2 3 4

blend
  surface n n
  2 5 6 3

blend
  surface n n
  6 7 4 3

blend
  line n 0 0 0 2
  5 6

blend
  line n 0 0 0 2
  6 7


disp
  1 0 0 -1e-2

boun
  1 0 0  1

eboun
  1 0 1 0

mate 1
solid
 elastic isotropic E nu
 density material d
 axis symm


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


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
  tplot
end
disp 21 2
disp 412 2
show


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


inte

stop