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 1.55e-32 n_p 6.67 m_u 5.667
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-4
set fobu 1e-2
set maxdt 1e-3
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
