求助--路基分层填筑过程沉降变形突变？
求助--路基分层填筑过程中考虑固结--加载时沉降变形突变，不仅仅在地表，地下10mi处也是在突变，不知何原因，请求各位帮忙给看看附：命令流（FLAC3D）new
title preloading_by_hanwj
config fluid
set fluid biot on
;; ------------ 建立网格------------
gen zone bri p0 0 0 0 p1 20 0 0 p2 0 2 0 p3 0 0 10 size 20 2 10group soil-5-clay
gen zone bri p0 20 0 0 p1 80 0 0 p2 20 2 0 p3 20 0 10 size 60 2 10group soil-5-clay
gen zone bri p0 0 0 10 p1 20 0 10 p2 0 2 10 p3 0 0 20 size 20 2 10group soil-4-clay
gen zone bri p0 20 0 10 p1 80 0 10 p2 20 2 10 p3 20 0 20 size 60 2 10group soil-4-clay
gen zone bri p0 0 0 20 p1 20 0 20 p2 0 2 20 p3 0 0 24 size 20 2 4group soil-3-clay
gen zone bri p0 20 0 20 p1 80 0 20 p2 20 2 20 p3 20 0 24 size 60 2 4group soil-3-clay
gen zone bri p0 0 0 24 p1 20 0 24 p2 0 2 24 p3 0 0 36 size 20 2 12group soil-2-clay
gen zone bri p0 20 0 24 p1 80 0 24 p2 20 2 24 p3 20 0 36 size 60 2 12group soil-2-clay
gen zone bri p0 0 0 36 p1 20 0 36 p2 0 2 36 p3 0 0 40 size 20 2 4group soil-1-clay
gen zone bri p0 20 0 36 p1 80 0 36 p2 20 2 36 p3 20 0 40 size 60 2 4group soil-1-clay
;; ------------建立路堤
gen zone bri p0 0 0 40 p1 20 0 40 p2 0 2 40 p3 0 0 43.72 p4 20 2 40 p5 0 2 43.72 &
p6 15.75 0 43.72 p7 15.75 2 43.72 size 20 2 6group soil-fill;;六层,每层0.62m
; ------------去掉路堤并将地基赋mc模型
model mohr
model null range z 39.9 43.73
;; ------------地基土赋参数;虚假值;干密度
prop bulk 5.16667E6 shear 1.10714E6 cohesion 18e10 friction 12.6 dilation 0.0 tension 1e10 range group soil-1-clay
ini density 1275.8 range group soil-1-clay
prop bulk 3.66667E6 shear 7.857E5 cohesion 14e10 friction 6.8 dilation 0.0 tension 1e10 range group soil-2-clay
ini density 1197.9 range group soil-2-clay
prop bulk 1.78337E7 shear 3.82142E6 cohesion 11e10 friction 31.0 dilation 0.0 tension 1e10 range group soil-3-clay
ini density 1662.4 range group soil-3-clay
prop bulk 9.5002E6 shear 2.03571E6 cohesion 24e10 friction 18.0 dilation 0.0 tension 1e10 range group soil-4-clay
ini density 1594.3 range group soil-4-clay
prop bulk 1.45E7 shear 3.10714E6 cohesion 14e10 friction 25.6 dilation 0.0 tension 1e10 range group soil-5-clay
ini density 1657.7 range group soil-5-clay
;;; ------------比奥参数设置
model fl_iso
model fl_null range group soil-fill
ini ftens -5e5
ini fdens 1000
ini biot_mod 0 ;;防止体积应变对孔压影响
prop biot_c 1
prop por=0.5342 perm=8.364E-12 range group soil-1-clay
prop por=0.5721 perm=1.377E-12 range group soil-2-clay
prop por=0.3876 perm=1.5496E-12 range group soil-3-clay
prop por=0.3857 perm=1.7238E-12 range group soil-4-clay
prop por=0.3523 perm=2.0604E-12 range group soil-5-clay
; ------------固定边界
fix x y z range z -0.1 0.1
fix x range x -0.1 0.1
fix x range x 79.9 80.1
fix y range y -0.1 0.01
fix y range y1.99 2.1
; ------------
;plot bl gr
; ------------
set grav 0 0 -10.0
water density 1000
water table origin 0 0 40 norm 0 0 1
; ------------; ------------
ini pp 400e3 grad 0 0 -10e3 range z 0 40
; ------------
set fluid off
solve
;------------地基土赋参数-真实值------------
prop bulk 5.16667E6 shear 1.10714E6 cohesion 18e3 friction 12.6 range group soil-1-clay
prop bulk 3.66667E6 shear 7.857E5 cohesion 14e3 friction 6.8 range group soil-2-clay
prop bulk 1.78337E7 shear 3.82142E6 cohesion 11e3 friction 31.0 range group soil-3-clay
prop bulk 9.5002E6 shear 2.03571E6 cohesion 24e3 friction 18.0 range group soil-4-clay
prop bulk 1.45E7 shear 3.10714E6 cohesion 14e3 friction 25.6 range group soil-5-clay
solve
;------------位移速度置0------------
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
ini state 0
;plot con szz outline on
;
fix pp 0 range z 39.9 40.1
;fix pp range x 79.9 80.1
;------------开启渗流并设置参数------------
set fluid on
ini biot_mod 1e8 range z -0.1 40.01
;ini biot_mod 6.63e6 rang group soil-1-clay ;modify biot_mod
;ini biot_mod 8.32e6 rang group soil-2-clay
;ini biot_mod 1.44e7 rang group soil-3-clay
;ini biot_mod 7.66e6 rang group soil-4-clay
;ini biot_mod 1.17e7 rang group soil-5-clay
;------------ 收敛条件及流体力学子步设置------------
;set mech force 1.5e3
set mech ratio 1e-4
set mech subs 10000 auto ;slave
set fluid subs 100
;------------ 记录------------
def day
day=fltime/24/3600
end
day
his id=1 unbal
his id=2 day
his id=3 fltime
his id=4 gp zd 1 1 38
his id=5 gp xd 1 1 38
his id=6 gp zd 1 1 30
his id=7 gp xd 1 1 30
his id=8 gp pp 1 1 38
his id=9 gp pp 1 1 30
his id=10 gp pp 1 1 10
his id=11 gp pp 1 1 0
plot bl state
plo add hist 4 6 vs 2
plo add hist 8 vs 2
save D:\FLAC4D-analysis\Consolidation-zizhong.sav
;------------------------------------
;第一层填筑
model mohr range z 39.9 40.63
prop bulk 8.33333E6 shear 1.78571E6 cohesion 50e4 friction 12.6range z 39.9 40.63
ini density 2000.0 range z 39.9 40.63
; --- solve to 15 days ---
sol age 1.296e6
save D:\FLAC4D-analysis\Consolidation-fill-1.sav
;
;第二层填筑
model mohr range z 40.61 41.25
prop bulk=8.33333E6 shear=1.78571E6 cohesion 50e4 friction 12.6range z 40.61 41.25
ini density 2000.0 range z 40.61 41.25
; --- solve to 30 days ---
sol age 2.592e6
save D:\FLAC4D-analysis\Consolidation-fill-2.sav
;
;第三层填筑
model mohr range z 41.23 41.87
prop bulk=8.33333E6 shear=1.78571E6 cohesion 50e4 friction 12.6 range z 41.23 41.87
ini density 2000.0 range z 41.23 41.87
; --- solve to 32 days ---
sol age 2764800
save D:\FLAC4D-analysis\Consolidation-fill-3.sav
;
;第四层填筑
model mohr range z 41.85 42.49
prop bulk=8.33333E6 shear=1.78571E6 cohesion 50e4 friction 12.6 range z 41.85 42.49
ini density 2000.0 range z 41.85 42.49
; --- solve to 34 days ---
sol age 2937600
save D:\FLAC4D-analysis\Consolidation-fill-4.sav
;
;第五层填筑
model mohr range z 42.47 43.11
prop bulk=8.33333E6 shear=1.78571E6 cohesion 50e4 friction 12.6 range z 42.47 43.11
ini density 2000.0 range z 42.47 43.11
; --- solve to 36 days ---
sol age 3110400
save D:\FLAC4D-analysis\Consolidation-fill-5.sav
;
;第六层填筑
model mohr range z 43.10 43.73
prop bulk=8.33333E6 shear=1.78571E6 cohesion 50e4 friction 12.6 range z 43.10 43.73
ini density 2000.0 range z 43.10 43.73
; --- solve to 180 days ---
sol age 15552000
save D:\FLAC4D-analysis\Consolidation-fill-6.sav
;
; --- solve to 180 days ---
sol age 9.5e7
save D:\FLAC4D-analysis\Consolidation-fill-7.sav
自己顶起来
没人关注吗 请问，您解决这个问题了吗 还是无人问津啦
页:
[1]