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

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

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

;第一步  平衡自重应力
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.814  shear 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 perm  1e-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

楼主，利用你的命令，想知道个时刻的孔隙水压力变化应该怎么修改？

账号		自动登录	找回密码
密码			注册

账号		自动登录	找回密码
密码			立即注册

[流固耦合] flac3d一维固结算例_计算结果供讨论

相关帖子

评分

老大

本帖子中包含更多资源