FLAC3D学习疑点详解Word文件下载.docx

1、初始网格gen zone radcyl p0 15 0 -15.75 p1 add 15 0 0 p2 add 0 6 0 p3 add 0 0 15.75 &size 3 4 8 5 dim 2.25 2.25 2.25 2.25 rat 1 1 1 1.1 fillgen zone radtun size 3 4 3 5 dim 0.75 2.25 0.75 2.25 rat 1 1 1 1.1 &p0 15 0 -15.75 p1 add 0 0 -14.25 p2 add 0 6 0 p3 add 15 0 0 fillgen zone reflect dip 90 dd 270 or

2、i 15 0 -15.75 摩尔库仑模型model mohrmacro niyan bulk 3.2825e3 shear 1.5150e3 fric 20 coh 0.07macro meiceng bulk 2.1735e3 shear 8.8917e2 fric 19 coh 0.05prop niyan range z =-13.5 0prop meiceng range z =-15.5 -13.5prop niyan range z = -30 -15.5为开挖定义整体部分range name in z -15.75 0 cyl end1 15 0 -15.75 end2 15 6

3、 -15.75 rad 2.25 range name hangdao x 12.75 17.25 y 0 6 z -16.50 -15.75进行加载fix x y z range z -30.001 -29.999app sxx=-20 range x -0.001 0.001app sxx=-20 range x 29.999 30.001app syy=-20 range y -0.001 0.001app syy=-20 range y 5.999 6.001app szz=-20 range z -0.001 0.001 solve（1）model null range in mod

4、el null range hangdaorange name in z -15.75 0 cyl end1 15 0 -15.75 end2 15 6 -15.75 rad 2.25range name hangdaoliangce1 x 12.749 12.751 y 0 6 z -16.50 -15.75range name hangdaoliangce2 x 17.249 17.251 y 0 6 z -16.50 -15.75range name hangdao x 12.75 17.25 y 0 6 z -16.501 -16.499apply nstress -0.05 plan

5、e dip 90 dd 0 range in 给巷道底部加压力app nstress -0.05 plane dip 90 dd 0 range hangdao 给巷道直墙部分（左侧）加压力apply nstress -0.05 plane dip 90 dd 0 range hangdaoliangce2 给巷道直墙部分（右侧）加压力 apply nstress -0.05 plane dip 90 dd 0 range hangdaoliangce1;对点进行监测hist unbalhist gp xdisp 12.75 0 -15.75:监测左帮hist gp zdisp 15 0 -1

6、3.5 ;监测顶板hist gp zdisp 15 0 -16.5 ;监测底板solve问：第一个solve平衡了,是第二个solve不能平衡，这是为什么？如果我把支护力加大，也不能平衡我遇到的不收敛的原因，有两个：1）检查你的模型是不是没有建好，相互之间是不是没有连接好；2）仔细检查你的参数 ，是不是已经发生了破坏。力施加有问题，隧道底部隆起破坏，所以无法计算我这个模型就是该模型就是要出现巷道底臌现象，如果按照你的意思，力应该怎样施加才对？ 如果无法收敛，应该用step命令，关键问题是要运行多少步。3.问：g 10 10m edef distortloop i （2,jzones）loop

7、j （2,jzones）x（i,j） =x（i,j）+mag*（urand-0.5）y（i,j）=y（i,j）+mag*（urand-0.5）end_loopendprop d 1000 sh 1e8 bu 2e8其中：x（i,j） =x（i,j）+mag*（urand-0.5）中的urand是否是随机生成数的命令？urand ：random number drawn from uniform distribution between0.0 and 1.0；grand ：random number drawn from normal distribution: mean = 0.0; stan

8、dard deviation = 1.04. 问：怎么能把一个PLOT的图像数据导出来以便用其他软件绘图？用set log on 命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。5. 问：用命令建立模型后，如何显示点的坐标？使用plo blo gro gpnum on命令6.interface问题我加了一个层面，初始平衡就出现了这么大的位移？是怎么原因？No1gen zone brick p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 8 &size 8 8 6 rat 1.02 1.02 1.02gen zone brick

9、 p0 0 0 9 p1 10 0 9 p2 0 10 9 p3 0 0 11 &size 8 8 2 rat 1.02 1.02 1.02gen zone reflect dip 90 dd 90 origin（0,0,6） range z 0 11group layer1 range z 9 11interface 1 face range plane norm （0,0,1） origin （0,0,8） dist 0.01ini z add -1.0 range group layer1No2gen zone brick p0 0 0 11 p1 10 0 11 p2 0 10 11

10、p3 0 0 13 &gen zone reflect dip 90 dd 90 origin（0,0,1） range z 11 13group layer2 range z 11 13interface 2 face range plane norm （0,0,1） origin （0,0,10） dist 0.01ini z add -1.0 range group layer2No3gen zone brick p0 0 0 13 p1 10 0 13 p2 0 10 13 p3 0 0 14 &size 8 8 1 rat 1.02 1.02 1.02gen zone reflect

11、 dip 90 dd 90 origin（0,0,1） range z 13 14group layer3 range z 13 14interface 3 face range plane norm （0,0,1） origin （0,0,12） dist 0.01ini z add -1.0 range group layer3No4gen zone brick p0 0 0 14 p1 10 0 14 p2 0 10 14 p3 0 0 15 &gen zone reflect dip 90 dd 90 origin（0,0,1） range z 14 15group layer4 ra

12、nge z 14 15interface 4 face range plane norm （0,0,1） origin （0,0,13） dist 0.01ini z add -1.0 range group layer4No5gen zone brick p0 0 0 15 p1 10 0 15 p2 0 10 15 p3 0 0 17 &gen zone reflect dip 90 dd 90 origin（0,0,15） range z 15 17group layer5 range z 15 17interface 5 face range plane norm （0,0,1） or

13、igin （0,0,14） dist 0.01ini z add -1.0 range group layer5No6gen zone brick p0 0 0 17 p1 10 0 17 p2 0 10 17 p3 0 0 18 &gen zone reflect dip 90 dd 90 origin（0,0,1） range z 17 18group layer6 range z 17 18interface 6 face range plane norm （0,0,1） origin （0,0,16） dist 0.01ini z add -1.0 range group layer6

14、No7gen zone brick p0 0 0 18 p1 10 0 18 p2 0 10 18 p3 0 0 20 &gen zone reflect dip 90 dd 90 origin（0,0,1） range z 18 20group layer7 range z 18 20interface 7 face range plane norm （0,0,1） origin （0,0,17） dist 0.01ini z add -1.0 range group layer7No8gen zone brick p0 0 0 20 p1 10 0 20 p2 0 10 20 p3 0 0

15、 24 &size 8 8 3 rat 1.02 1.02 1.02gen zone reflect dip 90 dd 90 origin（0,0,17） range z 20 24group layer8 range z 20 24interface 8 face range plane norm （0,0,1） origin （0,0,19） dist 0.01ini z add -1.0 range group layer8No9gen zone brick p0 0 0 24 p1 10 0 24 p2 0 10 24 p3 0 0 29 &size 8 8 4 rat 1.02 1

16、.02 1.02gen zone reflect dip 90 dd 90 origin（0,0,0） range z 24 29group layer9 range z 24 29interface 9 face range plane norm （0,0,1） origin （0,0,23） dist 0.01ini z add -1.0 range group layer9-边界条件-fix x range x -20.1 -19.9fix x range x 19.9 20.1fix y range y -0.1 0.1fix y range y 9.9 10.1fix z range

17、 z -0.1 0.1-初始条件-apply szz -26e6 range z 27.9 28.1initial szz -1.08e6 grad 0 0 2.6e4 range z 0 28set gravity 0 0 -10initial density 2.6e3initial sxx -45.8e6 grad 0 0 2.6e4 range z 0 28initial syy -45.8e6 grad 0 0 2.6e4 range z 0 28-材料属性-No1-pro bulk 8e9 shear 7e9 fric 25 coh 7.5e6 ten 2.6e6 range z

18、0 8 -No2-pro bulk 7.8e9 shear 3.5e9 fric 30 coh 6.5e6 ten 2.4e6 range z 8 10-No3-pro bulk 7e9 shear 3.0e9 fric 22 coh 5.5e6 ten 2.2e6 range z 10 12-No4-pro bulk 6e9 shear 2.3e9 fric 26 coh 5.5e6 ten 2.0e6 range z 12 13-No5-pro bulk 6.2e9 shear 2.0e9 fric 28 coh 5.2e6 ten 1.2e6 range z 13 14-No6-pro

19、bulk 5.8e9 shear 2.5e9 fric 35 coh 6.0e6 ten 1.3e6 range z 14 16-No7-pro bulk 9e9 shear 3.2e9 fric 30 coh 10.3e6 ten 1.2e6 range z 16 17-No8-pro bulk 8.1e9 shear 3.1e9 fric 26 coh 8.0e6 ten 1.5e6 range z 17 19-No9-pro bulk 7.6e9 shear 2.6e9 fric 23 coh 8.5e6 ten 1.8e6 range z 19 23-No10-pro bulk 7.5

20、e9 shear 2.1e9 fric 26 coh 8.0e6 ten 1.9e6 range z 23 28int 1 prop ks 8e10 kn 8e10 fric 30 dil 6 tension 2e7 bslip=onint 2 prop ks 8e10 kn 8e10 fric 30 dil 6 tension 2e7 bslip=onint 3 prop ks 8e10 kn 8e10 fric 35 dil 8 tension 2e7 bslip=onint 4 prop ks 9e11 kn 9e11 fric 32 dil 5 tension 2e7 bslip=on

21、int 5 prop ks 9e11 kn 9e11 fric 32 dil 5 tension 2e7 bslip=onint 6 prop ks 8e11 kn 8e11 fric 33 dil 5 tension 2e7 bslip=onint 7 prop ks 8e10 kn 8e10 fric 33 dil 5 tension 2e7 bslip=onint 8 prop ks 8e10 kn 8e10 fric 33 dil 5 tension 2e7 bslip=onint 9 prop ks 8e10 kn 8e10 fric 33 dil 5 tension 2e7 bsl

22、ip=onstep 5000 这是给的模型加的初始地应力，难道这样加地应力，模型就自动平衡了吗？你用ini平衡，为什么还要计算5000步？我的INI是用来施加初始地应力的，难道用INI那样加应力边界条件就不用step 或solve算平衡了吗？ 模型就能自动平衡吗？把上面的sxx和syy的剃度2.6e4改为9e3试一下，xx和yy方向的梯度不可能与zz方向的梯度一样的。7.问：prop mat 和zone定义的材料特性有什么区别？前者是对整个模型（如果整个模型只分成一个region的话），后者是对具体的zone的。8.关于ff边界设置的问题自由场地边界（FF边界）在动力分析中经常采用，它的原理是

23、采用粘滞阻尼器来模拟静止边界，将自由场节点的不平衡力加到主体的网格边界上。注意：模型底部的动力边界条件应当在app ff之前，进行app ff以后底部的边界条件就自动转为自由场地边界。现举出一简单例子,比如静力模型的尺寸为10*10*10,有以下几种边界施加方式1）free x y z rang z -.1 .1apply squiet dquiet nquiet rang z -.1 .1apply nstress value hist wave rang z -.1 .1apply ff2）3）可以验证,以上的施加方式得到的结果是完全一样的.所谓的自由场边界最重要的特征就是复制特征,现对以上各例子做简要分析.1）第二三两句就是在10*10的底面范围内分别设置静止边界和动力条件,当加了ff后,ff的网格的条件即从临近的maingrid的zone进行复制.所以ff的底面的位移边界和动力和maingrid完全一样.2）第二句直接加ff边界,注意这时ff底面是free的,而且动力为0,但是接下来的3.4两句仅仅给定的是z=-.1至.1的范围,导致ff的底面也包括在内了,所以作用方式和命令流1完全相同.3）同2的分析,同样可以得到结论但是对于以上三个命令流中的两句apply nst