hein_small.dat

海伦模型UDM自定义开发源代码

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 small

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 = 1

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 44956000 rat 0. f 0.
return