FLAC3D网上学习资料整理3

上传人:qt****68 文档编号:40199149 上传时间:2018-05-24 格式:DOC 页数:5 大小:165.50KB
返回 下载 相关 举报
FLAC3D网上学习资料整理3_第1页
第1页 / 共5页
FLAC3D网上学习资料整理3_第2页
第2页 / 共5页
FLAC3D网上学习资料整理3_第3页
第3页 / 共5页
FLAC3D网上学习资料整理3_第4页
第4页 / 共5页
FLAC3D网上学习资料整理3_第5页
第5页 / 共5页
亲,该文档总共5页,全部预览完了,如果喜欢就下载吧!
资源描述

《FLAC3D网上学习资料整理3》由会员分享,可在线阅读,更多相关《FLAC3D网上学习资料整理3(5页珍藏版)》请在金锄头文库上搜索。

1、FLAC3D 资料1网上网上 FLAC3D 学习笔记学习笔记15.建立建立 interface 的几种方法的几种方法首先给出一个简单的例子 newgen zone brick p0 0 0 0 p1 111.4 0 0 p2 0 10 0 p3 0 0 45 size 11 1 5 group 1gen zone brick p0 0 0 45 p1 111.4 0 45 p2 0 10 45 p3 0 0 90 p4 111.4 10 45 interface 1 face rang plane norm 45,0,64.4 origin 47 5 90interface 1 face ra

2、nge plane norm 45,0,64.4 origin 47 5 90 dist 0.1ini z add -0.1 range group 4gen zone brick p0 111.4 0 0 p1 128 0 0 p2 111.4 10 0 p3 111.4 0 45 ini xvel 0 yvel 0 zvel 0;set largestep 1000;solve fos file cable.sav associated 修改上面蓝色字体部分第一种建立 interface 的方法(书上的移来移去式):gen zone brick p0 0 0 0 p1 111.4 0 0

3、p2 0 10 0 p3 0 0 45 size 11 1 5 group 1gen zone brick p0 0 0 45 p1 111.4 0 45 p2 0 10 45 p3 0 0 90 p4 111.4 10 45 interface 1 face rang plane norm 45,0,64.4 origin 47 5 90interface 1 face range plane norm 45,0,64.4 origin 47 5 90 dist=0.1ini z add -0.1 range group 4gen zone brick p0 111.4 0 0 p1 128

4、 0 0 p2 111.4 10 0 p3 111.4 0 45 注意在此这个值给的很大。Init dens 1000Set gravity 0 0 -10Fix x range x -0.1 0.1Fix x range x 5.9 6.1Fix y range y -0.1 0.1Fix y range y 7.9 8.1Fix z range z -0.1 0.1set mech force=50solve;- mode null -Prop coh 1e3 tens 1e3 ;改为正常值(在此例中我们故意给小值)Mode null range x 2 ,4 y 2 , 6 z 5, 1

5、0 Set largeIni xdis 0 ydis 0 zdis 0 ;清零,不影响结果,为画图方便。Setp 2000return方法 2:得到初始应力的一个不错办法是设为弹性模型,然后赋予地勘信息,so 以后就 ok 了,得到初始应力了,再把相应的位移清零,下面就可以正式进入所分析的内容了。如上例还可以先设为弹性 (MODEL elastic),之后在开挖之前再设为 Mohr-Coulomb model, 就行。19关于初始应力的问题关于初始应力的问题Example 3.21 Initial stress state with gravitational gradientfor exam

6、ple, a 20 m 20 m 20 m box ofhomogeneous material at a depth of 200 m underground, with fixed base and stress boundaries onthe other sides.gen zone brick size 10,10,10 p1 20,0,0 p2 0,20,0 p3 0,0,20model mohrprop bulk 5e9 shear 3e9 fric 35ini density 2500set gravity 0,0,-10fix x y z range z -0.1 0.1in

7、i szz = -5.0e6 grad 0,0,2.5e4ini sxx = -2.5e6 grad 0,0,1.25e4ini syy = -2.5e6 grad 0,0,1.25e4apply szz = -4.5e6 range z 19.9 20.1apply szz = -5.0e6 range z -0.1 0.1FLAC3D 资料3apply sxx = -2.5e6 grad 0,0,1.25e4 range x -0.1 0.1apply sxx = -2.5e6 grad 0,0,1.25e4 range x 19.9 20.1apply syy = -2.5e6 grad

8、 0,0,1.25e4 range y -0.1 0.1apply syy = -2.5e6 grad 0,0,1.25e4 range y 19.9 20.1它的意思是通过 ini szz 来赋予初始地应力,但是为何又有 set gravity 0,0,-10?这样不是重复了吗?ini 是模型的内力,set grav 是外力施加对模型的作用力,没有外力怎么有内力?一般而言,这样是为了计算比较快,因为这样设置以后,此时几乎已经平衡,了,计算步较少了,变形也就很少了嘛,当然最好再在以后计算时初始化位移为 0,但初始化位移与否(在这种情况下),并不会明显影响计算结果的.当然还有另一种方法就是开始时

9、设成弹性的(model elas)(或者还是 mohr的但参数给的要大一些,目的是防止破坏) ,直接运行到平衡达来获取初始应力值,在真正计算时再改成 mohr 并将参数改为真实的数值,这种方法更适合于单元不太多,因为比较简单但如果单元数目很多的话,那么此种方法速度会很慢的关于初始地应力及其算法,有两个疑问:1)为何一定要让模型达到初始平衡再进行计算?2.)那种改大参数,让模型在自重下平衡得到初始应力场的做法中,要将哪些参数改大呢?bul she coh ?还是改其他哪几个?回答:1)这是模拟多数现场施工前的实际情况,并非总是如此,例如土钉墙施工,如果开挖后立即加土钉,就无须先平衡再加土钉再SO

10、LVE.2) 取决于具体模型,核心思想是用该模型在此阶段模拟弹性模型.对于 MOHR-COULOMB 模型, 这样做(将COHESION,TENSION 设置很大的值)与 SOLVE ELASTIC 是等价的.但对于别的模型,SOLVE ELASTIC 并不适用,最好用改变参数的方法来加快计算速度和防止屈服.另外一点: 在FLAC/FLAC3D 中, SET 命令作用于整个网格(赋值全局变量),INI 则初始化指定范围内网格点/域的偏移量.初始应力不仅仅是重力引起的。它还包括土压力,水浮力,上方建筑物的力,这些信息的勘测一般都不容易获得。通常可以把模型材料设为弹性,再把两个模量设到很大 sol

11、ve 后再把模量和材料还原可得。具体可以看以前关于初始应力的帖子。不只是这样的,set grav 不是直接将力仅仅加在底层上,而是作用在上面,如果没有初始应力的话,地层就会有相应的变形直到地层变形和产生的力与之平衡。所以应当通过 initial 来初始化应力,而这个力与 density 相匹配。我是做地铁开挖的如果没有初始化的话地层在开挖前就有了沉降。这与实际不符合。设置重力的原因不是为了平衡减少程序计算步骤,而是因为如果未初始化应力的话,在 set grav 的时候,模型会在自重应力下发生应力变形。所以,初始化应力的原因是防止模型在自重应力下发生变形(我们认为地层自重变形已经发生,当然欠固结

12、等现象是另外一回事) 。其实就是让模型发生变形,然后把位移清零也可以的。有种方法,不少人也在用,就是为了得到初始应力值,先将材料参数放大到很大,获得初始值以后再恢复到正常参数。20. Tension-p shear-p tension-n shear-n 表示的含义表示的含义在计算的循环里面,每个循环中,每个 zone 都依据 failure criterion 处于不同的状态,shear 和 tension 分别表示因受剪和受拉而处于塑性状态.n 表示 now,p 表示 previous,即分别表示在现在和以前的循环当中处于塑性状态.21. 有一个例子有一个例子:gen zone cyl p

13、0 0 0 0 p1 1 0 0 p2 0 2 0 p3 0 0 1 size 4 5 4gen zone reflect norm 1,0,0gen zone reflect norm 0,0,1model mohrprop bulk 1.19e10 shear 1.1e10prop coh 2.72e5 fric 44 ten 2e5fix x y z range y -.1 .1fix x y z range y 1.9 2.1ini yvel 1e-7 range y -.1 .1ini yvel -1e-7 range y 1.9 2.1hist gp ydisp 0,0,0hist

14、 zone syy 0,1,0hist zone syy 1,1,0hist write 1 file 原点位移.dathist write 2 file 中部单元应力.dathist write 3 file 边缘单元应力.dathist write 1 table 1 hist write 2 table 2hist write 3 table 3hist write 2 3 vs 1 file 单压.datstep 3000结果在“原点位移.dat”文件中只有:“Step 1 Y-Displace” 而没有保存下数值,怎么回事?我怎么才能获得 table 中的数据?修改:gen zone

15、 cyl p0 0 0 0 p1 1 0 0 p2 0 2 0 p3 0 0 1 size 4 5 4gen zone reflect norm 1,0,0gen zone reflect norm 0,0,1model mohrprop bulk 1.19e10 shear 1.1e10prop coh 2.72e5 fric 44 ten 2e5fix x y z range y -.1 .1fix x y z range y 1.9 2.1ini yvel 1e-7 range y -.1 .1ini yvel -1e-7 range y 1.9 2.1hist gp ydisp 0,0

16、,0hist zone syy 0,1,0hist zone syy 1,1,0FLAC3D 资料4step 3000hist write 1 file 原点位移.dathist write 2 file 中部单元应力.dathist write 3 file 边缘单元应力.dathist write 1 table 1 hist write 2 table 2hist write 3 table 3hist write 2 3 vs 1 file 单压.datsave a print tabel n will solve your second qusetion22.这个这个 plot 是什么意思?是什么意思? hist n 1hist gp ydisp 0,0,0hist ax_strhist gp xdisp 1,1,0step 1500plot hist 2 vs -1 ;axial stress vs axial disp.plot hist 3 vs -1 ; circum

展开阅读全文
相关资源
正为您匹配相似的精品文档
相关搜索

最新文档


当前位置:首页 > 高等教育 > 教育学

电脑版 |金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号