salzer_example_5.dat

马尔克斯模型UDM自定义开发源代码

; deformation of a rock salt pillar
;
config creep cppudm
model load salzer.dll
set large
;
grid 11 34
m salzer i=1,11    j=1,34       ; - Na2
;
; Na2
prop bu=16000 sh=10000 den=0.0025 i=1,11 j=1,34
prop a_p=1.7e-22 n_p=10.0 m_u=1.6 t_null=87600 e0=1e-7 i=1,11 j=1,34
;
; Na2
gen   0,-535   0,-450  27.5,-450  27.5,-535  i=1,12 j=1,35
;
fix x i=1
fix x i=12
fix x y j=1
;
set grav 10.
ini syy -13.375 var 0,2.125  j=1,35
ini sxx -13.375 var 0,2.125  j=1,35
ini szz -13.375 var 0,2.125  j=1,35
app press 11.25 i=1,12 j=35
;
def jahre
  jahre = crtime/8760
end
;
his unbal
his crtime
his jahre
his ydis i=12 j=35
his yvel i=12 j=35
his xdis i=8 j=12
his xvel i=8 j=12
;
set fobl=0.01 fobu=5 min=1e-25 max=10000 lmul=1.1 umul=0.9 lat=25
;set fobl=0.1 fobu=10 min=1e-25 max=1000 lmul=1.1 umul=0.8 lat=25
set crdt 0
set force 1e-30
set srat  1e-30
set step 9000000
;
; model settlement
step 500
save bep0.aus
ini xdisp=0
ini ydisp=0
;
;
;  0 years
; drivage of chamber
m n i=8,11 j=5,18
int 1 a from 8,19 to 12,19 b from 8,5 to 8,19
int 1 kn 1000 ks 1000 fric 10 coh 0
int 2 a from 8,5 to 12,5 b from 8,5 to 8,19
int 2 kn 1000 ks 1000 fric 10 coh 0
;
set crdt auto
;  10 years
solve age 87600
save bep11.aus
;
;  50 years
solve age 438000
save bep12.aus
;
; 100 years
solve age 876000
save bep13.aus
;
ret