fby_engineer 发表于 2008-7-18 10:22:25

flac3d一维固结算例_计算结果供讨论

;一维固结例子
; 单级等速大面积竖向均布加载的单层均质地基压缩问题

;第一步平衡自重应力
new
config fluid
gen zone brick size 1 1 10 p0 0 0 -10 p1 1 0 -10 p2 0 1 -10 p3 0 0 0

model fl_iso
prop perm 1e-12 poro 0.3
ini fdensity 9.81e2
ini fmod 2e9 ftens -1e-3
set flow off
set mech on

;fix pp range z -10.1 -9.9

model elas
prop bulk 2512562.814shear 1152959.262 dens 1335.7
fix x range x -0.1 0.1
fix x range x 0.9 1.0
fix y range y -0.1 0.1
fix y range y 0.9 1.0
fix z range z -10.1 -9.9
ini szz 0 grad 0 0 15990.3 range z -10 0
ini syy 0 grad 0 0 12604.69 range z -10 0
ini sxx 0 grad 0 0 12604.69 range z -10 0
ini pp 0 grad 0 0 -9810 range z -10 0

;fix pp range z -0.1 0.1
;fix pp range z -10.1 -9.9
plot surface
set gravity 0 0 -9.81

solve
plot con pp
save gujie0.sav


;第二步流固耦合求解
set flow on
set mech on

model fl_iso
prop perm1e-12 poro 0.3
ini fdensity 981
ini fmod 2e9 ftens -1e-3
;app nstr -5e3 range z 0 x 0 1 y 0 1
;ini pp 1e4
fix pp 0 range z 0.1 -0.1



plot surface
;set mech substep 1000
;solve age 302400
plot con zdis


app nstr -1e4 range z 0 x 0 1 y 0 1
solve age 604800

app nstr -2e4 range z 0 x 0 1 y 0 1
solve age 1209600

app nstr -3e4 range z 0 x 0 1 y 0 1
solve age 1814400

app nstr -4e4 range z 0 x 0 1 y 0 1
solve age 2419200


app nstr -5e4 range z 0 x 0 1 y 0 1
solve age 3024000

app nstr -6e4 range z 0 x 0 1 y 0 1
solve age 3628800

app nstr -7e4 range z 0 x 0 1 y 0 1
solve age 4233600

app nstr -8e4 range z 0 x 0 1 y 0 1
solve age 4838400

app nstr -9e4 range z 0 x 0 1 y 0 1
solve age 5443200

app nstr -10e4 range z 0 x 0 1 y 0 1
solve age 6048000


solve age 6652800

solve age 7257600

solve age 7862400

solve age 8467200

      solve age 9072000

solve age 9504000

solve age 9676800

solve age 10886400

solve age 12096000

solve age 13305600

solve age 14515200

solve age 15724800

solve age 18144000

solve age 20563200

solve age 22982400

solve age 25401600

[ 本帖最后由 fby_engineer 于 2008-7-19 16:35 编辑 ]

fby_engineer 发表于 2008-7-18 10:24:23

计算的最终沉降与《岩土工程有限元原理与应用》的例子相同
但中间加载步，计算的沉降略微偏小，有兴趣的各位高手，可以帮忙分析下

fby_engineer 发表于 2008-7-19 16:38:05

时间/d        解析解/cm        本人的解/cm 解析解/本人的解
7        0.294        0.1905        1.543307087
14        0.835        0.53475        1.561477326
21        1.526        1.022        1.493150685
28        2.35        1.6428        1.430484539
35        3.284        2.3881        1.375151794
42        4.316        3.2495        1.328204339
49        5.439        4.2192        1.28910694
56        6.644        5.2899        1.255978374
63        7.926        6.4548        1.227923406
70        9.28        7.7075        1.204022056
77        10.405        8.867        1.173452126
84        11.352        9.9483        1.141099484
91        12.203        10.957        1.113717258
98        12.978        11.898        1.090771558
110        13.694        12.775        1.071937378
112        14.361        13.365        1.074523008
126        15.568        15.068        1.033182904
140        16.63        16.35        1.017125382
154        17.566        17.466        1.005725409
168        18.394        18.437        0.997667733
182        19.125        19.281        0.991909133
210        20.342        20.654        0.984893967
238        21.294        21.693        0.98160697
266        22.036        22.479        0.980292718
294        22.617        23.074        0.980194158

这是计算结果，开始误差很大，在50%，随着加载和时间的推移
到294天是能够比较接近，
请高人帮忙看一下，这个是什么原因啊？

xiaomag 发表于 2008-7-20 22:05:39

老大

老大，你算得真不错，能告诉我你的QQ吗，一起交流一下。

fby_engineer 发表于 2008-7-20 23:23:46

没问题啊，42123043
我还在搞这个问题
这两天又有些新体会，一起交流

070542 发表于 2008-7-21 13:14:14

任何软件做数值模拟开始结果都不好，这个是肯定的。尽量将拟合不好的那段时间缩短就行。

070542 发表于 2008-7-21 14:34:02

我以前做过一维固结的验证，发现还在，一维的情况还挺不错的。

fby_engineer 发表于 2008-7-21 15:44:33

好的，非常感谢，下下来看看

fby_engineer 发表于 2008-7-21 15:48:13

您做的结果非常不错，有没有源代码可以传上来，给大家分享分享
或者已经在里面了？
我只看到三张图

minhrac 发表于 2008-7-21 16:30:03

只有结果，还希望把命令流贴上来

flashchen 发表于 2009-12-3 18:02:41

很好，学习了

科技大虾 发表于 2010-5-22 23:48:10

外荷载随时间线性变化怎么加？

xjbincumt 发表于 2010-5-27 17:34:51

挺好的，学习了

1984 发表于 2010-5-28 15:25:27

好的，非常感谢，下下来看看

科技大虾 发表于 2010-5-29 23:49:10

楼主，利用你的命令，想知道个时刻的孔隙水压力变化应该怎么修改？
页: [1]
查看完整版本: flac3d一维固结算例_计算结果供讨论