***behavior gen_evp
  **elasticity  isotropic
    young 2e5
    poisson 0.3
***return