flac3d常见问题总结

1、-. z.flac3d学习有一段时间，过一段时间没学究忘了很多，所以学习最好不要放下来。对自己总结，也是学习的一种。flac3d命令很多，用到的不多，就则几个，根底的就不多说了进入正题：个人认为如果学习flac3d最好把flac的命令进展总结一下。做模型就好比拼盘：要地应力，就把地应力的那段命令摘入，要出图，参加plot的相关命令，要出动画，加个动画命令，所以总结好命令后，以后做模型，只要把命令组一条条的对应放入，一个程序就出来了，后期调试一下，改相关参数就可以了。接下来，我把自己弄好的整理一下，按照我上面的思路进展建模什么人建模都会有自己的喜好，下面说一下我自己建模第一步就是建立一个全新的t

2、*t文件不管命令大小都是如此，这样的好处是，模型可以修改方便快捷。而且还可以自我查错建立新的t*t后，第一行命令就是new，加new的目的是为了防止其他命令的影响相当于刷新再执行其他命令前刷新一下。不用想，接下来就是建模了，建立模型，不同人有不同的方法，有的人喜欢建完后在对模型开挖最终得到自己想要的模型。model null 等；有的人喜欢用各种带缺角的组合成最终的模型。具体建模看各自的喜好，我觉得，不同点是建模的网格划分有区别，用带缺角的组合模型网格划分更好最好把常用的集中模型记下来组成各种句子，以后用的时候直接插入例如：gen zone brick/radbrick size *,*,*,

3、* ratio *,*,*,*建模的过程中，或同一个文件中有重复应用的命令，就可以用宏命令了，宏命令macro 名称 命令的容例：macro shiping plot con zdip outline on用shiping 代替后面的引号里的东西。在建模中只要有一样的都可以用shiping代替，缩减了书写命令的时间，整体上也能让人看的更明白注意点是，你命名的符号最好跟falc自带的命令不要重复，我的处理方法是，名称后带数字如shiping1如果你想修改宏对象，改变命令的容，在名称上加单引号就是替换的意思例macro shiping1 新的命令这些很多课本里都有讲，不在表达很多时候我们要对单一的

4、对象进展操作，例如进展模拟运算后，只想单独的显示一局部的容，单一的土层，锚杆，等就要用到groupgroup给对象命名，例：在gen zone brick - size - group namename随便取，不要跟flac的已有的命令重合就行这里的意思就是对这个体进展了命名接下来想把这一块模型挖掉加model null range group namegroup的用法：1。group直接加到模型的后面相当于给整个模型命名2.group name range * *yy zz给一个围的体命名继续学习，昨天说到了建模，我是搞巷道的，重点讲巷道。巷道建模，用到brick，radbrick，radt

5、unnel，radcylinder这4个模型比拟多，一般模型都是建一半，让后用镜像命令复制另一半，建模也是大家有各自不同的习惯，总的有两种，一种用点法，就是吧模型的所有关键点都输入进去，这种方法可以用来建立不规则的模型，弊端就是工作量大，不当在设点的时候考虑多，输入的时候也很麻烦。第二种。按照模型只要取三个点确定*yz三周的方法，有点输入简单，弊端，只能建立规则的模型，像要建立，路堤模型，边坡模型，就只能是按照数点的方法建模range的命令：书上有，但是我用自己的话说一下他的用途1：取圆球面的围，就是同一个中心，不同半径之间夹的那一局部，用在给球面施加力，约束时用该命令，命令：range an

6、nulus center *yz中心坐标radius r1 r2两半径想要表达的是球面时，r1r2取值相差小点即可，如radius 2.9 3.1 这时围就是半径为3的面。2。圆柱的环面施加力或约束时也用到该命令。命令range cylinder end1 *yz end2 *yz radius r(end是表示柱体断面的中心点坐标)range plane above 表示的是z大于0的上部构造range plane blow 表示的是z小于0的下部构造range dd a dip b origin *yz 表示的是走向角度a，倾向角度b的面dd表示，当面与*oy面相交时的线的垂线与y轴正向的

7、夹角，顺为正，说白话，就是面的指向方向面是朝那个方向的意思。dip面与*oy面的夹角，与z轴负向为正，origin表示面经过*yz点range normal *yz 法线方向为*yz的面range not 没试用过围外的东西，想想用到圆上取围非常好现在回到建模上面说过建模一般只建一半然后用reflect命令镜像镜像首先得有个面吧，就目前我用的面都是*oy面yoz面。估计大家也很少用到斜面镜像吧面的命令上面range里说过了这里就不讲了命令：reflect normal *yz origin *yz 根据法线和过一点建的面最常用至于实际使用过程中有人问道：normal -1 0 0 与norma

8、l 1 0 0的区别，我试过两者的效果是一样的，没区别，虽然方向不同，但表达的是同一个面reflect dd a dip b origin *yz 这个就是用里建立斜面镜像的，各个符号的意思上面讲过，不重复了有一点注意了，镜像的时候如果是模型镜像的话，相当于在建立一个模型，故，在reflect的前面要加建模的关键字，gen zone reflect。先讲点后期出图的，然后总结一下plot set rotation 20 0 30 设置观测角度的plot set center auto 设置观测中心点plot set magnification 1.5 设置放大倍数的plot show显示这一整

9、局部就是我开头说的要自己整理出来的句段了，等以后用到的时候直接插入，作用就是显示模型的。任何模型的显示都可以用他来表示。很好很强大。如果还要显示其他功能就把相应的命令在插入上面那局部就行了plot set rotation 20 0 30 plot set center autoplot set magnification 1.5 plot add surface yellowplot add a*esplot show这个就是我常用来查看模型的命令，大家也可以自己整理一分，总结，做个实例，给大家看看隧道，半圆直墙隧道分两步一，建模，二，显示建模 建一半然后镜像对吧尺寸，自己在草纸上画好，点坐

10、标都弄清楚了把以前的句子直接套下来gen zone radcylinder size 4 25 6 14网格数量，自己看着办ratio1 1 1 1.2&p0 0 0 0 p1 50 0 0 p2 0 100 0 p3 0 0 50 dimension 2 2 2这是建上半局部gen zone rad tunnel size 4 25 6 14 ratio 1 1 1 1.2 &p0 0 0 0 p1 0 0 -50 p2 0 100 0 p3 50 0 0 dimension 1.8 2 2建立下半部局部建一半后，用镜像命令也讲过了直接插入，gen zone reflect normal -

11、1 0 0 origin 0 0 0最后一步就是显示给大家看那一打包显示直接插入plot set rotation 20 0 30 plot set center autoplot set magnification 1.5 plot add surface yellowplot add a*esplot show。整体的建模加显示，的命令就完了合一块 建立一个新的t*t 开头加一个newnewgen zone radcylinder size 4 25 6 14 ratio1 1 1 1.2&p0 0 0 0 p1 50 0 0 p2 0 100 0 p3 0 0 50 dimension

12、2 2 2gen zone rad tunnel size 4 25 6 14 ratio 1 1 1 1.2 &p0 0 0 0 p1 0 0 -50 p2 0 100 0 p3 50 0 0 dimension 1.8 2 2gen zone reflect normal -1 0 0 origin 0 0 0plot set rotation 20 0 30 plot set center autoplot set magnification 1.5 plot add surface yellowplot add a*esplot show整体就这样了，大家可以粘贴一下看看。plot命令

13、大集合：显示切片的命令plot set plane ori 0 0 0 norm 0 1 0plot con zd plane这两行命令，用来显示指定面上的位移变化云图。吧第二行改改，还可以显示应力云图，等其他plot con syy y向应力云图plot con szz outline on z向应力云图加显示网格plot con szz effective on z向有效应力plot con szz inter 1e4 云图显示增量后面的值越小，图上颜色条带就越细，plot con szz min -10e4 显示最小值，图上小于这个数的值都用同一色表示plot con szz ma* 1

14、0e4 显示最大值，图上大于这个数的值都用同一色表示plot fip outside on 显示力的箭头建模完事后就要给模型赋值，赋参数了，关键词，modelmodel mor赋摩尔库伦模型，想要什么模型改相应的词pro den 2500密度bulk2e8体积模量，shear 2e8 剪切模量co 2e6聚力fric45摩擦角ten1e6 抗拉强度。range.(我们要赋值的围)给模型赋值就是hi这么简单，难在参数的选取。现在看边界，如果是不动面，用fi* * range。如果边界带速率，或位移用ini。set grav 0 0，-9.81ini dens 1000设置重力就一定要设置密度不考

15、虑重力就不要设计密度。两者共存的，设计程序的时候一定要记的set mech force 50solve用SET命令设定了最大不平衡力的上限值为50N,即：计算的最大不平衡力一旦小于此值，我们便人为模型到达了平衡状态。求解过程终止。设置这个目的是在一定精度要求下，加快计算值越大，算的量越少，花的时间也越少说到这里，add与set的区别add原来有的加新的用addset原来没有的，设置用setplot create GravVplot set plane dip=90 dd=0 origin=3,4,0plot set rot 15 0 20plot set center 2.5 4.2 4.0p

16、lot add bound behindplot add bcont szz planeplot add a*esplot show 。如这个命令，显示图像的时候命令有很多，则我们想显示别的图像，又不想吧这个弄好的图像丢掉则用plot create name建立一个图像，在建立新的图像时plot create name1新的名显示以前的图像则输入如下命令：plot current name。另一种方法，输入：plot show name这里有一点最好开场的时候将网格剖分的相对稀疏一些10001500个单元左右，先求解一下，看看结果，然后再根据大致的结果有针对性的增加模型的网格数，以得到更准确的

17、解答。更快的检验一下自己 的模型是否建错了不同的区域定义不同的本构模型model mohr range z=0,5prop ClayeyGravel range z=0,5在定义参数后面加个range围为了防止在输入命令时总要通过range全局坐标域的方法来限定该命令执行的区域，我们可以先通过range命令来给不同的区域命名，然后就可以用range区域名的方法来限定围了。 这样就防止了每次输入命令要限定在该区域执行时，都要输入一遍坐标围的麻烦。具体如何给区域命名，以及日后如何用区域名来限定该区域，可以看下面这个例子：range name trench *=-1,1 y=0,4 z=-2,2mo

18、del null range trench今天做方案，没时间，就讲几点得了，1：右手法则，建立模型时，的点的排放时有规律的，方法如下，拿出右手伸直了4指合并指向就是p0到p1的方向掌心指向就是p0到p2的方向拇指与4指垂直的方向就是p0到p3的方向很多人说网格后的云图有带齿什么的，改良的方法就是重复用下面的命令Gen ADJust 用此命令来微调网格, 使之离散平滑。该命令可以连续使用以增进平滑效果的。plo blo gro可以看到不同的group的颜色分布例不同的土层定义组后，就可以用此命令显示不同的色了如何查看剖面上任一点的位移值？print gp disp range *() y( )

19、z( )或print gp disp range idplot gp disp range *() y( ) z( )或plot gp disp range id很多人都是这么做的，我也是学习，initial 与 apply 有何区别？initial初始化命令，如初始化计算体的应力状态等；apply边界条件限制命令，如施加边界的力、位移等约束等。initial的应力状态会随计算过程的发生而发生改变，一般体力需要初始化，而apply施加的边界条件不会发生变化如果想看模型在计算过程中的变化情况Set movie avi step 1 file 1.avi；定义动画，动画记录的频率为1时步一频率计算

20、量大的时候就没必要1时步的记录，可以10,100时步的记录，减少电脑的使用Movie start；动画开场hist n= 5 ：N(Nstep)=5指定了每迭代计算5次记录一次相关的值如果不指定，则系统默认值为10hist unbal 命令的第二行指定记录最大不平衡力the ma*imum unbalanced forcehist gp zdisp 4,4,8 记录坐标值为4，4，8的节点Z方向的位移这个用的多，查看收敛量时就是这个命令，上面两个个人人认为用处不大，也很少用过如何在程序中进展暂停，并可恢复计算？答：在命令中参加pause命令，用continue进展继续个人认为这个命令非常好用在

21、我们分步求解中想得到*一个过程中的结果，不用等到全求完，还可以在分布求解错误的时候就进展改正。而不是等到结果出来关于dd、dip的理解：不同人有不同的理解，我的理解就是dip面与水平面的夹角dd是面与水平面的交线的垂线与y轴的夹角1) 在plot hist m vs n的形式里,m代表y轴,n代表*轴(不管m,n的正负);5 L0 D# _1 1 C- L2) -表示对其值作mirror,如果大家使用AUTOCAD的话就明白mirror的含义了,比方上面的 hist 1 gp ydisp的值是0100,则vs -1就变成-1000.以此类推.3) 为什么要这么做呢主要是为了符合岩土工程的习惯需

22、要,例如桩载荷试验曲线,沉降曲线等都不是画在第一象限.最近看有很多对后处理中的收敛有很多不同的想法，自己也总结一下plot sk dis scale 0.07显示位移矢量箭头的后面的scale默认的是0.053越大箭头越大magf也是放大变形的意思，但是不能放在sk dis的后面，计算时机出错回到主题，收敛数据，我们想得到一条收敛量随时间变化，或者不同应力直间关系的图线方法1用set log on然后把数据导入e*cel或者其他的处理软件进展处理2.用hist id=1不写id的话，电脑会自己按123456的顺序进展排序gp点，如果是体的话用zonezdis *yz*yz点处z向位移这里应力应

23、变之类的都是用zone，位移之类的都是用gp最后相要出相应的结果用命令plot hist 1；如果相要的不同的id之间的关系比方说竖向应力szz与剪应力s*z的关系命令plot hist 2 v 3 2,3是对应应力开场设的id号flac3d程序里面自动设置的步长是10，也就是说每隔10步记录一次，这样如果模型小得出的结果就数据少图线不光滑，如果模型大，得到的数据多计算文件又太大，浪费空间。设置步长的命令set hist_rep n(n就是步长)上面说的plot多是flac里显示的如果想把结果存到别处去用命令hist write 2 v 3 file namehist.t*t(其实2相当于*轴

24、，3相当于y轴)设置影片格式上面见过，温故而知新set movie avi step 1 存储步长，可以任意设，就是每隔几步存储的意思file name.avimovie start分步开挖就是分布计算每多一步多添一块土，或多开挖，多掘进等就计算一下，加一个solve即可如果想保存每一步的结果和图，那就在每一步的solve后面添加保存命令设立相应的文件，以后调用即可1.设定初始地应力场是十分必要的 没有地应力场进展计算会使模型在重力作用下发生过大的变形 从而不能与实际情况相符2。apply是施加的外力 ini是初始的力。3. 有人问如果改参数这样做的话地应力不是被人为改变了吗？为什么要人为加大

25、它的参数呢，很不合理？答复：不能改变bulk和shear模量，因为这两个参数是用来计算应力的。而con和fric是用来判断破坏状态的。如果将bulk和shear模量改变了，那你在重力条件下计算的应力就不真实了。咱们改大参数的目的是让模型在不破坏的条件下得到正确的应力，所以只有将con和fric变大才能到达目的。 当然如果在弹性模型下计算初始应力就没有这种问题了。因为弹性模型是没有破坏准侧的当然也没有描述到达破坏状态的参数了。4. 大家正好讨论地应力的问题，我也有个疑问，平衡地应力后，水平方向的地应力如何确定的？是由本构模型来确定的吗？这个问题现在我也没答案也是我的问题5. gen zone b

26、rick size 10,10,10 p1 20,0,0 p2 0,20,0 p3 0,0,20model mohrprop bulk 5e9 shear 3e9 fric 35ini density 2500 set gravity 0,0,-10fi* * y z range z -0.1 0.1ini szz = -5.0e6 grad 0,0,2.5e4ini s* = -2.5e6 grad 0,0,1.25e4ini syy = -2.5e6 grad 0,0,1.25eapply szz = -4.5e6 range z 19.9 20.1apply szz = -5.0e6 r

27、ange z -0.1 0.17apply s* = -2.5e6 grad 0,0,1.25e4 range * -0.1 0.1apply s* = -2.5e6 grad 0,0,1.25e4 range * 19.9 20.1apply syy = -2.5e6 grad 0,0,1.25e4 range y -0.1 0.1apply syy = -2.5e6 grad 0,0,1.25e4 range y 19.9 20.1 U# P! y- P$ _+ D它的意思是通过ini szz来赋予初始地应力,但是为何又有set gravity 0,0,-10这样不是重复了吗 ini是确实

28、是模型的力，这个力通过set grav ,然后就像一些帖子上介绍的那样,让程序计算不就行了吗干嘛还有手动在设定一次呢我感觉它这样的做法有点重复了. 答复：1它这样做其实只是减少程序的计算步骤.5 _这样设置以后 整个模型就已经是平衡了2设置重力的原因不是为了平衡减少程序计算步骤，而是因为如果未初始化应力的话，在set grav的时候，模型会在自重应力下发生应力变形。所以，初始化应力的原因是防止模型在自重应力下发生变形我们认为地层自重变形已经发生，当然欠固结等现象是另外一回事。3其实就让模型发生变形，然后把位移清零也可以的 6。ini命令来对地层进展初始应力的施加，而APPLY是7 u(施加应力边界条件7. 用到ini初始平衡的时候要和重力计算的一样才有用8。gen zone brick size 10 10 10 &p0 0,-25,0 p1 20,-25,0 p2 0,0,0 p3 0,-25,