我也用的是3d呀。下面就是flac3d用户手册里的一个power蠕变模型的例子
也是用 solve age 2000。时间单位为年,那就是2000年了。
;-------------------------------------------------------------------
; cylindrical cavity -- power law
;
;-------------------------------------------------------------------
;
title
Power-Law Creep Model --- Cylindrical Cavity
;
config creep
;
gen zone brick p0 1 0 0 p1 20 0 0 &
p3 0.98769 0.15643 0 p6 19.754 3.129 0 &
p2 1 0 -1 p4 20 0 -1 &
p5 0.98769 0.15643 -1 p7 19.754 3.129 -1 &
size 20 1 1 ratio 1.2 1 1
model power
;
; Properties and stresses in Pascal units (not MPa)
prop bulk=1e9 shear=3e8
prop a_1=1e-25 n_1=3
;
fix y z range y -0.01 0.01
apply nvel 0 plane dip 90 dd -9 orig 0 0 0 range y 0.1 100
apply dvel 0 plane dip 90 dd -9 orig 0 0 0 range y 0.1 100
;
apply nstr -100e6 range x 19 410
;
ini sxx -100e6 syy -100e6 szz -100e6
;
hist unbal
hist gp xdis id 1
hist gp xdis id 81
;
step 2000
;
save cyl_1.sav
;
hist reset
hist unbal
hist gp xvel id 1
hist gp xvel id 81
hist gp xdis id 1
hist gp xdis id 81
hist z sxx id 1
hist z syy id 1
hist z szz id 1
hist crtime
hist dt
;
ini xv 0 yv 0 zv 0
;
set creep on
set cr dt 1.0e-4
set cr mindt 1.0e-4
set cr maxdt 1.0
set cr dt auto on
;
solve age 2000
;
save cyl_2.sav
;
ca ex2_06.dat ; compare analytic solution......
;
save cyl_3.sav
;
table 1 name Flac3d
table 11 name analytical
table 12 name 'Flac3d radial stress'
table 2 name 'analytical radial stress'
table 13 name 'Flac3d hoop stress'
table 3 name 'analytical hoop stress'
;
plot tab 1 11 line
pause
plot tab 12 line black 2 13 line 3
;
save ex2_05.sav
ret