- 积分
- 0
- 注册时间
- 2011-8-8
- 仿真币
-
- 最后登录
- 1970-1-1
|
本帖最后由 颠倒 于 2011-12-3 13:23 编辑
我是菜鸟,刚学flac3d,编了个程序但问题一大堆,自己解决不了,连哪个地方出问题了都不是很清楚,到这里求助前辈高手啊! 先把基本情况说说:这个程序模拟软土地基上箱涵的受力与软基变形的协调关系。假设涵洞和涵洞基础呈刚性,路堤范围内软基处理后土体性质与路堤以外部分有一定不同。模型中原点在涵洞底部中央,以下15m为软基,以上10.8m为涵洞和路堤部分。计算模型以路堤中线为分界截取一半进行研究。 - new
- config fluid
- generate zone radtunnel size 6 11 12 18 ratio 1 1 1 1.15 &
- p0=(0,0,0) p1=(27,0,0) p2=(0,29.2,0) p3=(0,0,10.8) p4=(27,29.2,0) p5=(0,23,10.8) p6=(27,0,10.8) p7=(27,23,10.8) &
- p8=(2.4,0,0) p9=(0,0,4.8) p10=(2.4,29.2,0) p11=(0,22,4.8) p12=(2.4,0,4.8) p13=(2.4,22,4.8) &
- group embankment
- generate zone brick size 6 11 12 &
- p0=(0,0,0) p1=(2.4,0,0) p2=(0,29.2,0) p3=(0,0,4.8) p4=(2.4,29.2,0) p5=(0,22,4.8) p6=(2.4,0,4.8) p7=(2.4,22,4.8) &
- group culvert
- generate zone brick size 10 11 6 ratio 1.15 1 1 &
- p0=(0,0,0), p1=(0,0,-15) p2=(0,29.2,0) p3=(2.4,0,0) &
- group foundation
- generate zone brick size 10 11 18 ratio 1.15 1 1.15 &
- p0=(2.4,0,0) p1=(2.4,0,-15) p2=(2.4,29.2,0) p3=(27,0,0) &
- group plastic_foundation
- generate zone brick size 10 6 6 ratio 1.15 1 1 &
- p0=(0,29.2,0) p1=(0,29.2,-15) p2=(0,42,0) p3=(2.4,29.2,0) &
- group ground1
- generate zone brick size 10 6 18 ratio 1.15 1 1.15 &
- p0=(2.4,29.2,0) p1=(2.4,29.2,-15) p2=(2.4,42,0) p3=(27,29.2,0) &
- group ground2
- group ground range group ground1 any group ground2 any
- generate zone reflect normal (1,0,0) origin =(0,0,0)
- model null range x -2 2 y 0 29.2 z 0.4 4.4 ; 这个命令的作用到后面怎么体现不出来了呀?算完初始应力,涵洞就变成实心的了。。。
- plot block group
- pause
- fix x y z range z -14.9 -15.1
- fix x range x -26.9 -27.1
- fix x range x 26.9 27.1
- fix y range y -0.1 0.1
- fix y range y 41.9 42.1
- model mohr
- model null range group embankment
- prop bulk 1.52e7 shear 1.05e7 coh 1e10 tens 1e10 range group culvert
- ini dens 2500 ran group culvert
- prop bulk 1.52e7 shear 1.05e7 coh 1e10 tens 1e10 range group foundation
- ini dens 2500 ran group foundation
- prop bulk 1.6e6 shear 1.2e6 coh 1e10 tension 1e10 range group ground
- ini dens 1800 ran group ground
- prop bulk 1.2e6 shear 9e5 coh 1e10 tension 1e10 range group plastic_foundation
- ini dens 1800 ran group plastic_foundation
- model fl_iso
- ini pp 0 grad 00 -10e3 range z 0 -15 ; 这样定义孔压是否得体?我的原点在涵洞底部中央,下面是软土。为什么我后面算出来初始应力的孔压图中只分了两层,体现不了渐变,而且凸出去的涵洞居然会存在孔压?很震惊。。。
- ini ftens -1e-3
- ini fdens 1000
- ini fmod 2e8
- prop perm 1e-11 poro 0.45 range group ground any group plastic_foundation any
- model fl_null range group foundation ; 这个命令说涵洞基础没有孔压,是不是应该把所有没孔压的地方都用这个命令给套一套?
- fix pp range x -26.9 -27.1 ; 我把四周都做了透水边界,合理不?
- fix pp range x 26.9 27.1
- fix pp range y -0.1 0.1
- fix pp range y 41.9 42.1
- set fluid off
- set gravity 0,0,-9.8
- plot con zdis sk
- hist id=1 unbal
- solve
- prop bulk 1.52e7 shear 1.05e7 fric 50 coh 7e7 range group culvert
- prop bulk 1.52e7 shear 1.05e7 fric 50 coh 7e7 range group foundation
- prop bulk 1.6e6 shear 1.2e6 fric 8 coh 1e4 range group ground
- prop bulk 1.2e6 shear 9e5 fric 8 coh 1e4 range group plastic_foundation ; 说实话我一直不明白为啥要赋值两变,而且项目还不一样(前一个没内摩擦角,后一个没抗拉强度)?书上说为了不屈服,但还是不咋理解,我参考的书上的。。。
- solve
- save 初始应力.sav
- plot con szz outline on
- plot con syy outline on
- ini xdis 0 ydis 0 zdis 0
- ini xvel 0 yvel 0 zvel 0
- set fluid on
- set mech force1.5e3
- set mech subs 10000 auto
- set fluid subs 30 ; 这三句话说的是主从进程法求解流固耦合问题吧?但具体赋值我真不知道该怎么弄。这里是参考的一位前辈大侠的,估计也不怎么适用。。。
- def day
- day=fltime/24/3600
- end
- day
- hist id=2 day
- model mohr ran z 0 1.2 group embankment ; 因为我考虑的不只是路堤的重量,还要考虑它对结构物受力的影响,所以就用的摩尔-库伦模型
- prop bulk 1.11e7 shear 8.3e6 coh 3e3 fric 36 ran z 0 1.2 group embankment
- ini dens 2100 ran z 0 1.2 group embankment
- hist id=3 gp zdis 4.2 0 0 ; 后面的很多个“hist”命令中好多id所表示的都是同一点的同一个应变,只是荷载不同所以大小不一样,不知道是否可以把它们归一?不知道这个命令是否占用很大内存,如果不占的话这样是否也行?
- hist id=4 gp zdis 6 0 0
- hist id=5 gp zdis 7.8 0 0
- hist id=6 gp xforce 2.4 0 0.4
- solve age 1.0368e6
- save fill-1.sav
- pause
- model mohr ran z 1.2 2.4 group embankment
- prop bulk 1.11e7 shear 8.3e6 coh 3e3 fric 36 ran z 1.2 2.4 group embankment
- ini dens 2100 ran 1.2 2.4 group embankment
- hist id=7 gp zdis 4.2 0 0
- hist id=8 gp zdis 6 0 0
- hist id=9 gp zdis 7.8 0 0
- hist id=10 gp xforce 2.4 0 0.4
- solve age 1.0368e6
- save fill-2.sav
- model mohr ran z 2.4 3.6 group embankment
- prop bulk 1.11e7 shear 8.3e6 coh 3e3 fric 36 ran z 2.4 3.6 group embankment
- ini dens 2100 ran z 2.4 3.6 group embankment
- hist id=11 gp zdis 4.2 0 0
- hist id=12 gp zdis 6 0 0
- hist id=13 gp zdis 7.8 0 0
- hist id=14 gp xforce 2.4 0 0.4
- hist id=15 gp xforce 2.4 0 2.4
- solve age 1.0368e6
- save fill-3.sav
- model mohr ran z 3.6 4.8 group embankment
- prop bulk 1.11e7 shear 8.3e6 coh 3e3 fric 36 ran z 3.6 4.8 group embankment
- ini dens 2100 ran z 3.6 4.8 group embankment
- hist id=16 gp zdis 4.2 0 0
- hist id=17 gp zdis 6 0 0
- hist id=18 gp zdis 7.8 0 0
- hist id=19 gp zdis 4.2 0 4.8
- hist id=20 gp zdis 6 0 4.8
- hist id=21 gp zdis 7.8 0 4.8
- hist id=22 gp xforce 2.4 0 0.4
- hist id=23 gp xforce 2.4 0 2.4
- hist id=24 gp xforce 2.4 0 4.4
- solve age 1.0368e6
- save fill-4.sav
- model mohr ran z 4.8 5.8 group embankment
- prop bulk 1.11e7 shear 8.3e6 coh 3e3 fric 36 ran z 4.8 5.8 group embankment
- ini dens 2100 ran z 4.8 5.8 group embankment
- hist id=25 gp zdis 4.2 0 0
- hist id=26 gp zdis 6 0 0
- hist id=27 gp zdis 7.8 0 0
- hist id=28 gp zdis 4.2 0 4.8
- hist id=29 gp zdis 6 0 4.8
- hist id=30 gp zdis 7.8 0 4.8
- hist id=31 gp xforce 2.4 0 0.4
- hist id=32 gp xforce 2.4 0 2.4
- hist id=33 gp xforce 2.4 0 4.4
- hist id=34 gp zforce 0 0 4.8
- hist id=35 gp zforce 2.0 0 4.8
- hist id=36 gp zforce 4.2 0 4.8
- hist id=37 gp zforce 6 0 4.8
- hist id=38 gp zforce 7.8 0 4.8
- solve age 8.64e5
- save fill-5.sav
复制代码
后面几层的加载编程就暂不列举了,反正也是差不多的循环。。。 但编程完毕后到目前为止,我就只算了个没有进行路堤加载时的初始应力图,也不知道有什么用,而且结果还是错的。。。接下来的第一层填土算了一天没有算完,估计也是错了,菜鸟呀这是。。。
问几个不懂的问题,有热心的前辈高手帮忙解答个,算是提携晚辈呀!!!
1. 为什么建模过程中我把涵洞中心用null命令掏空了,但在后来计算的初始应力图中仍然显示为填充状态?
2. 为什么程序中孔压定义在z=0到-15之间(ini pp 0 grad 0 0 -10e3 range z 0 -15),初始应力算完后涵洞上却分布有孔压?上面我没有设置水呀?
3. 这个命令(model fl_null range groupfoundation)是否可以让指定部分没有孔压了?这样的话是不是应该把所有不应该存在孔压的地方都用这个命令套一套?
4. 为什么我的初始应力计算出来以后孔压只分为两层,而且没有渐变的过程,显然不对,问题可能出在了哪呀?
5. 后面我记录采样的id有很多,但其中不少id所表示的都是同一个地点同一个变量,但是在不同荷载条件下大小可能不同。想问一下有没有这个必要,是不是能够只用一个id表示不同荷载下的同一个变量?如果用同一个id表示,后处理过程中取用是否方便?还有一点,记录采样是否会耗用很大的计算空间?如果不占什么地方的话是否就可以用多个id去表示了。。。。。。
6. 后面有个算流固耦合的命令不理解,直接copy别个的,也不知道可行不可行,我想自己算了一天也算不完第一层填土的原因就在这里:
“set mech force 1.5e3
set mech subs10000 auto
set fluid subs30”
书上介绍的好像是主从进程法,但看了半天不知所云,领悟力差呀,有前辈高手帮忙指点个,要是觉得痛心疾首了,想骂随便骂!
7. 还有好多错误我自己看不出来的,求前辈高手指正啊!!!感激不尽,感激不尽啊!!!!!
程序在这里:
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
评分
-
1
查看全部评分
-
|