zhuhaha 发表于 2006-12-21 11:28:31

我这儿有两种方法，不知道哪个好，还是都正确？
第一种：
mod mohr
prop bulk=6.67e9 shear=4e9 fric=45 coh=1e6 tens 2e6
inter 1 prop kn 1e8 ks 1e8 friction 15 cohesion 4e5
fix x range x -20.1 -19.9
fix x range x 19.9 20.1
fix z range z -20.1,-19.9
fix y range y -0.1,0.1
fix y range y 99.9,100.1
set grav 0,0,-10
ini density=2500
apply szz -2.5e6 range z 19.9 20.1
apply sxx 0.86e6 range x -20.1 -19.9
apply sxx -0.86e6 range x 19.9 20.1
;;;;;;;
set large
hist unbal
solve



第二种：

mod mohr
prop bulk=6.67e9 shear=4e9 fric=45 coh=1e6 tens 2e6
inter 1 prop kn 1e8 ks 1e8 friction 15 cohesion 4e5
fix x range x -20.1 -19.9
fix x range x 19.9 20.1
fix z range z -20.1,-19.9
fix y range y -0.1,0.1
fix y range y 99.9,100.1
set grav 0,0,-10
ini density=2500
ini szz -2.5e6 range z 19.9 20.1
ini sxx 1.25e6 range x -20.1 -19.9
ini sxx -1.25e6 range x 19.9 20.1
;;;;;;;
set large
hist unbal
solve

[ 本帖最后由 zhuhaha 于 2006-12-21 11:33 编辑 ]

zhuhaha 发表于 2006-12-21 11:32:27

手册上例子：
ini和apply不是重复了吗？
mod mohr
prop bulk=1e8 shear=3e8 fric=35 coh=1e3 tens=1e3
fix x range x -0.1,0.1
fix z range z -0.1,0.1
fix y range y 9.9,10.1
set grav 0,10,0
ini density=1000
ini syy=0.0 grad 0,-10000,0
ini sxx=0.0 grad 0,-5000,0
ini szz=0.0 grad 0,-5000,0
apply sxx=0.0 grad 0,-5000,0 range x 9.9,10.1
apply szz=0.0 grad 0,-5000,0 range z 9.9,10.1

[ 本帖最后由 zhuhaha 于 2006-12-21 11:34 编辑 ]

zhuhaha 发表于 2006-12-21 11:36:13

还是设成elas，参数加大？
mod elas
inter 1 prop kn 1e9 ks 1e9 friction 15 cohesion 4e5
fix x range x -20.1 -19.9
fix x range x 19.9 20.1
fix z range z -20.1,-19.9
fix y range y -0.1,0.1
fix y range y 99.9,100.1
set grav 0,0,-10
ini density=2500
apply szz -2.5e6 range z 19.9 20.1
apply sxx 1.25e6 range x -20.1 -19.9
apply sxx -1.25e6 range x 19.9 20.1

benjackxu 发表于 2006-12-21 11:38:20

方法都可，可以看LAKEWATER的PPT。岩土工程数值模拟不是追求精确度，而是力求反映规律。

hedan 发表于 2007-3-24 14:53:47

回复 #22 zhuhaha 的帖子

楼上把ini 和 apply命令混淆了，ini命令来对地层进行初始应力的施加，而APPLY是
施加应力边界条件

lzquan 发表于 2007-3-24 23:45:17

好贴阿 高手真多

zhangh2008 发表于 2007-6-13 01:42:15

FLAC3D沉降问题

我做了个桩承载力的模拟,主要是求沉降量,土体用mohr模型,但在自重下总得不到均匀沉降的位移云图,后来用elas模型初始平衡,平衡后再将模型改回mohr,但改回后如果elas中的参数与土参数一致则出现在加载荷时土体沉降不合理(在桩顶和桩中部出现异常,且沉降量过小),如果elas中的参数设的很大则在改回mohr模型时就出现不均匀沉降,但加力后沉降却合理.请高手指点

z1972 发表于 2007-6-13 21:00:55

GRAV用于平衡地应力场,若没有则初始状态将处于不平衡状态.

lyxwolf 发表于 2009-2-20 13:40:48

认真看了这个帖子，让我少了很多疑惑啊，好贴。

wangshuyang 发表于 2009-3-3 15:58:02

8# ruigrate
学习到了

mazhongyuan9049 发表于 2009-3-3 16:02:24

设定初始地应力场是十分必要的 没有地应力场进行计算会使模型在重力作用下发生过大的变形 从而不能与实际情况相符。

mazhongyuan9049 发表于 2009-3-3 16:07:06

关于初始地应力及其算法，有两个疑问一直没弄明白，请牛人帮忙解答一下：1. 为何一定要让模型达到初始平衡再进行计算？
2. 那种改大参数，让模型在自重下平衡得到初始应力场的做法中，要将哪些参数改大呢？bulshe...
real 发表于 2005-7-12 22:09 http://forum.simwe.com/images/common/back.gif

初始平衡的作用可以看做是地层在自然状态下经过很长时间的固结所成的状态
初始弹性计算可以将K和G取很大（10^9~10） 弹塑性计算将位移至零只留地应力场进行计算

arnoldray 发表于 2009-3-6 22:07:52

真是济济一堂的牛人啊

arnoldray 发表于 2009-3-6 22:07:20

真是济济一堂的牛人啊

whlxpp521 发表于 2009-3-31 20:35:41

m e 一般不用设置kn和ks吧

li_2001 发表于 2009-3-31 21:21:39

一般而言这样是为了计算比较快,因为这样设置以后,此时几乎已经平衡,了,计算步较少了,变形也就很少了嘛,当然最好再在以后计算时初始化位移为0,但初始化位移与否(在这种情况下),并不会明显影响计算结果的.当然还有另一种方法就是开始时设成弹性的(model elas)（或者还是mohr的但参数给的要大一些，目的是防止破坏），直接运行到平衡达来获取初始应力值，在真正计算时再改成mohr并将参数改为真实的数值，这种方法更适合于单元不太多，因为比较简单．但如果单元数目很多的话，那么此种方法速度会很慢的

li_2001 发表于 2009-3-31 21:22:25

apply是施加的外力ini是初始的内力

da2001881 发表于 2010-3-12 14:58:47

一般而言这样是为了计算比较快,因为这样设置以后,此时几乎已经平衡,了,计算步较少了,变形也就很少了嘛,当然最好再在以后计算时初始化位移为0,但初始化位移与否(在这种情况下),并不会明显影响计算结果的.当然还有另一 ...
li_2001 发表于 2009-3-31 21:21 http://forum.simwe.com/images/common/back.gif

这样做的话地应力不是被人为改变了吗？为什么要人为加大它的参数呢，很不合理

tianshen 发表于 2010-3-22 12:48:10

大家正好讨论地应力的问题，我也有个疑问，平衡地应力后，水平方向的地应力如何确定的？是由本构模型来确定的吗？

tianshen 发表于 2010-3-22 12:48:57

39# tianshen
我是说 内部的水平地应力是通过什么方式计算的
页: 1 [2] 3
查看完整版本: 手册上一个关于初始地应力的例子---我觉得有点奇怪哦