config creep
set large
gr 2,1

model salzer

def _emodul
_emod=7e3
_ny=0.27
_bulk=_emod/(3*(1-2*_ny))
_shear=_emod/(2*(1+_ny))
end
_emodul

prop bulk _bulk shear _shear dens 0.0024
prop a_p 2.97e-27 n_p 5.26 m_u 4.26
prop t_null 10 e0 1e-5

fix y j=1
fix x i=2
ini syy -20
app press 20 j 2
hist unbal
hist nstep=20 syy sxx i=1 j=1
hist ydis i 1 j 2
hist xdisp i 1 j 2
def seff
    seff= syy(1,1)-sxx(1,1)
end
hist seff

hist crtime
hist unbal

set fobl  1e-7
set fobu  1e-5
set maxdt 1e-4
set mindt 1e-10
set lmul  1.01
set umul  0.5
prop ssfa 30

set crdt auto
solve force 0.0 age 100

ret