config cppudm creep
model load ..\models\hein
gen zone brick p0=(0,0,0) p1=(1,0,0) p2=(0,1,0) p3=(0,0,1) size 1 1 1
model hein
set large
prop dens 1000
prop K_0=1.812E+04 c_k=9.06544611969
prop G_0=9.843E+03 c_g=9.06544611969
prop eta_0 0.310
prop temperature 373
prop struct_fac = 1.48910222316209E-18
prop gas_c = 1.0
prop act_energy = 17531.4480962282
prop n_stress 8.943640857
prop c_1 2.869698044
prop c_2 -0.153360351
prop c_3 0.59
prop c_4 0.126172031
prop c_5 1.361506092
prop c_6 6.700607645
prop c_7 0.887852607
prop c_8 0.001
prop evol_complete = 0
def time_step
z_eta = z_prop(zone_head, 'poros')
z_evol = z_prop(zone_head, 'evol')
z_bulk = z_prop(zone_head, 'bulk')
z_shear = z_prop(zone_head, 'shear')
z_vol = z_volume(zone_head)
time_step = crdt
days = crtime / 86400.
end
hist crtime
hist days
hist time_step
hist z_eta
hist z_evol
hist z_bulk
hist z_shear
hist z_vol
hist zone szz 1,1,1
hist nstep 200
fix x y z ra z -0.1 0.1
fix x ra x -1.0,0.1
fix x ra x 0.9,1.1
fix y ra y 0.9,1.1
fix y ra y -1.0,0.1
fix zvel -6.9e-9 ra z 0.9 1.1
set creep on
set creep mindt 1e-24 maxdt 100
set creep dt auto on
set creep ufob 1e-2
set creep lfob 5e-2
solve age 39120000 rat 0. f 0.
return
