flac3d常见问题总结

flac3d学习有一段时间,过一段时间没学究忘了很多,所以学习最好不要放下来。对自己总结,也就是学习的一种。flac3d命令很多,用到的不多,就那么几个,基础的就不多说了进入正题:个人认为如果学习flac3d最好把flac的命令进行总结一下。做模型就好比拼盘:要地应力，就把地应力的那段命令摘入,要出图,加入plot的相关命令,要出动画,加个动画命令,所以总结好命令后,以后做模型,只要把命令组一条条的对应放入，一个程序就出来了，后期调试一下，改相关参数就可以了。接下来，我把自己弄好的整理一下，按照我上面的思路进行建模什么人建模都会有自己的喜好，下面说一下我自己建模第一步就就是建立一个全新的txt文件（不论命令大小都就是如此），这样的好处就是，模型可以修改方便快捷。而且还可以自我查错建立新的txt后，第一行命令就就是new,加new的目的就是为了防止其她命令的影响（相当于刷新）再执行其她命令前刷新一下。不用想，接下来就就是建模了，建立模型,不同人有不同的方法，有的人喜欢建完后在对模型开挖最终得到自己想要的模型。modelnull等;有的人喜欢用各种带缺角的组合成最终的模型。具体建模瞧各自的喜好，我觉得，不同点就是建模的网格划分有区别,用带缺角的组合模型网格划分更好最好把常用的集中模型记下来组成各种句子，以后用的时候直接插入例如：genzonebrick/radbricksize*,*,*,*ratio*,*,*,*建模的过程中，或同一个文件中有重复应用的命令,就可以用宏命令了，宏命令macro名称‘命令的内容’例：macroshiping'plotconzdipoutlineon'用shiping代替后面的引号里的东西。在建模中只要有相同的都可以用shiping代替,缩减了书写命令的时间,整体上也能让人瞧的更明白注意点就是,您命名的符号最好跟falc自带的命令不要重复，我的处理方法就是，名称后带数字如shiping1如果您想修改宏对象，改变命令的内容，在名称上加单引号就就是替换的意思例macro'shiping1‘‘新的命令’这些很多课本里都有讲，不在叙述很多时候我们要对单一的对象进行操作，例如进行模拟运算后，只想单独的显示一部分的内容，单一的土层，锚杆，等就要用到groupgroup给对象命名，例:在genzonebricksize----groupname（name随便取，不要跟flac的已有的命令重合就行）这里的意思就就是对这个体进行了命名接下来想把这一块模型挖掉加modelnullrangegroupnamegroup的用法:1。group直接加到模型的后面相当于给整个模型命名2、groupnamerangexxyyzz给一个范围的体命名继续学习，昨天说到了建模，我就是搞巷道的，重点讲巷道。巷道建模，用到brick,radbrick,radtunnel,radcylinder这4个模型比较多，一般模型都就是建一半，让后用镜像命令复制另一半，建模也就是大家有各自不同的习惯，总的有两种，一种用点法，就就是吧模型的所有关键点都输入进去，这种方法可以用来建立不规则的模型,弊端就就是工作量大，不当在设点的时候考虑多，输入的时候也很麻烦。第二种。按照模型只要取三个点确定xyz三周的方法,有点输入简单，弊端，只能建立规则的模型，像要建立,路堤模型，边坡模型，就只能就是按照数点的方法建模range的命令:书上有，但就是我用自己的话说一下她的用途1:取圆球面的范围，就就是同一个中心，不同半径之间夹的那一部分，用在给球面施加力，约束时用该命令，命令：rangeannuluscenterxyz（中心坐标）radiusr1r2（两半径）想要表达的就是球面时，门r2取值相差小点即可，如radius2、93、1这时范围就就是半径为3的面。2。圆柱的环面施加力或约束时也用到该命令。命令rangecylinderend1xyzend2xyzradiusr（end就是表示柱体断面的中心点坐标）rangeplaneabove表示的就是z大于0的上部结构rangeplaneblow表示的就是z小于0的下部结构rangeddadipboriginxyz表示的就是走向角度a,倾向角度b的面dd表示,当面与xoy面相交时的线的垂线与y轴正向的夹角，顺为正，说白话,就就是面的指向方向面就是朝那个方向的意思。dip面与xoy面的夹角，与z轴负向为正，origin表示面经过xyz点rangenormalxyz法线方向为xyz的面rangenot没试用过（范围外的东西，想想用到圆上取范围非常好）现在回到建模上面说过建模一般只建一半然后用reflect命令（镜像）镜像首先得有个面吧，就目前我用的面都就是xoy面yoz面。估计大家也很少用到斜面镜像吧面的命令上面range里说过了这里就不讲了命令:reflectnormalxyzoriginxyz（根据法线与过一点建的面）最常用至于实际使用过程中有人问道:normal-100与normal100的区别，我试过两者的效果就是一样的,没区别，虽然方向不同，但表达的就是同一个面reflectddadipboriginxyz这个就就是用里建立斜面镜像的，各个符号的意思上面讲过，不重复了有一点注意了，镜像的时候如果就是模型镜像的话，相当于在建立一个模型，故，在reflect的前面要加建模的关键字,genzonereflect。。。。。。先讲点后期出图的，然后总结一下（plotsetrotation20030设置观测角度的plotsetcenterauto设置观测中心点plotsetmagnification1、5设置放大倍数的plotshow）显示这一整部分就就是我开头说的要自己整理出来的句段了，等以后用到的时候直接插入，作用就就是显示模型的。任何模型的显示都可以用她来表示。很好很强大。如果还要显示其她功能就把相应的命令在插入上面那部分就行了plotsetrotation20030plotsetcenterautoplotsetmagnification1、5plotaddsurfaceyellowplotaddaxesplotshow这个就就是我常用来查瞧模型的命令,大家也可以自己整理一分,总结，做个实例，给大家瞧瞧隧道,半圆直墙隧道分两步一,建模，二，显示建模建一半然后镜像对吧尺寸，自己在草纸上画好，点坐标都弄清楚了把以前的句子直接套下来genzoneradcylindersize425614（网格数量，自己瞧着办）ratio1111、2&p0000p15000p201000p30050dimension222这就是建上半部分genzoneradtunnelsize425614ratio1111、2&p0000p100-50p201000p35000dimension1、822建立下半部部分建一半后,用镜像命令也讲过了直接插入，genzonereflectnormal-100origin000最后一步就就是显示给大家瞧那一打包显示直接插入plotsetrotation20030plotsetcenterautoplotsetmagnification1、5plotaddsurfaceyellowplotaddaxesplotshow。整体的建模加显示，的命令就完了合一块建立一个新的txt开头加一个newnewgenzoneradcylindersize425614ratio1111、2&TOC\o"1-5"\h\zp0000p15000p201000p30050dimension222genzoneradtunnelsize425614ratio1111、2&p0000p100-50p201000p35000dimension1、822genzonereflectnormal-100origin000plotsetrotation20030plotsetcenterautoplotsetmagnification1、5plotaddsurfaceyellowplotaddaxesplotshow整体就这样了，大家可以粘贴一下瞧瞧。plot命令大集合：显示切片的命令plotsetplaneori000norm010plotconzdplane这两行命令，用来显示指定面上的位移变化云图。（吧第二行改改,还可以显示应力云图，等其她）plotconsyyy向应力云图plotconszzoutlineonz向应力云图加显示网格plotconszzeffectiveonz向有效应力plotconszzinter1e4云图显示增量（后面的值越小,图上颜色条带就越细,）plotconszzmin-10e4显示最小值，图上小于这个数的值都用同一色表示plotconszzmax10e4显示最大值，图上大于这个数的值都用同一色表示plotfipoutsideon显示力的箭头建模完事后就要给模型赋值，赋参数了,关键词,modelmodelmor（赋摩尔库伦模型,想要什么模型改相应的词）proden2500（密度）bulk2e8体积模量,shear2e8（剪切模量）co2e6（内聚力）fric45（摩擦角）ten1e6抗拉强度。range、、、、、、、、、、（我们要赋值的范围）给模型赋值就就是hi这么简单，难在参数的选取。现在瞧边界，如果就是不动面，用fixxrange。。。。。如果边界带速率，或位移用ini。。。setgrav00,-9、81inidens1000设置重力就一定要设置密度不考虑重力就不要设计密度。两者共存的，设计程序的时候一定要记的setmechforce50solve用SET命令设定了最大不平衡力的上限值为50N,即:计算的最大不平衡力一旦小于此值，我们便人为模型达到了平衡状态。求解过程终止。设置这个目的就是在一定精度要求下，加快计算（值越大，算的量越少，花的时间也越少）说到这里,add与set的区别add原来有的加新的用addset原来没有的,设置用setplotcreateGravVplotsetplanedip=90dd=0origin=3,4,0plotsetrot15020plotsetcenter2、54、24、0plotaddboundbehindplotaddbcontszzplaneplotaddaxesplotshow。如这个命令，显示图像的时候命令有很多，那么我们想显示别的图像，又不想吧这个弄好的图像丢掉则用plotcreatename建立一个图像，在建立新的图像时plotcreatename1（新的名）显示以前的图像那么输入如下命令:plotcurrentname。另一种方法，输入:plotshowname这里有一点最好开始的时候将网格剖分的相对稀疏一些（1000〜1500个单元左右）,先求解一下,瞧瞧结果,然后再根据大致的结果有针对性的增加模型的网格数,以得到更精确的解答。更快的检验一下自己的模型就是否建错了不同的区域定义不同的本构模型modelmohrrangez=0,5propClayeyGravelrangez=0,5在定义参数后面加个range范围为了避免在输入命令时总要通过range+全局坐标域的方法来限定该命令执行的区域，我们可以先通过range命令来给不同的区域命名,然后就可以用range+区域名的方法来限定范围了。这样就避免了每次输入命令要限定在该区域执行时，都要输入一遍坐标范围的麻烦。具体如何给区域命名,以及日后如何用区域名来限定该区域，可以瞧下面这个例子：rangenametrenchx=-1,1y=0,4z=-2,2modelnullrangetrench今天做方案,没时间，就讲几点得了,1:右手法则，建立模型时，的点的排放时有规律的，方法如下，拿出右手伸直了4指合并指向就就是pO到p1的方向掌心指向就就是pO到p2的方向拇指与4指垂直的方向就就是p0到p3的方向很多人说网格后的云图有带齿什么的,改进的方法就就是重复用下面的命令GenADJust用此命令来微调网格，使之离散平滑。该命令可以连续使用以增进平滑效果的。ploblogro可以瞧到不同的group的颜色分布(例不同的土层定义组后,就可以用此命令显示不同的色了)如何查瞧剖面上任一点的位移值？printgpdisprangex()y()z()或printgpdisprangeidplotgpdisprangex()y()z()或plotgpdisprangeid很多人都就是这么做的，我也就是学习，initial与apply有何区别？initial初始化命令，如初始化计算体的应力状态等;apply边界条件限制命令，如施加边界的力、位移等约束等。initial的应力状态会随计算过程的发生而发生改变，一般体力需要初始化，而apply施加的边界条件不会发生变化如果想瞧模型在计算过程中的变化情况Setmovieavistep1淀义动画，动画记录的频率为1时步一频率(计算量大的时候就没必要1时步的记录,可以10,100时步的记录，减少电脑的使用)Moviestart;动画开始histn=5:N（Nstep）=5指定了每迭代计算5次记录一次相关的值（如果不指定，则系统默认值为10）histunbal命令的第二行指定记录最大不平衡力（themaximumunbalancedforce）histgpzdisp4,4,8记录坐标值为（4,4,8）的节点Z方向的位移（这个用的多，查瞧收敛量时就就是这个命令，上面两个个人人认为用处不大，也很少用过）如何在程序中进行暂停，并可恢复计算？答:在命令中加入pause命令，用continue进行继续个人认为这个命令非常好用在我们分步求解中想得到某一个过程中的结果，不用等到全求完，还可以在分布求解错误的时候就进行改正。而不就是等到结果出来关于dd、dip的理解：不同人有不同的理解，我的理解就就是dip面与水平面的夹角dd就是面与水平面的交线的垂线与y轴的夹角1）在plothistmvsn的形式里,m代表y轴,n代表x轴（不管m,n的正负）;5L0D#_1~1C-L2）"-"表示对其值作"mirror",如果大家使用AUTOCAD的话就明白"mirror"的含义了，比如上面的hist1gpydisp的值就是0~100,那么vs-1就变成-100~0、以此类推、3）为什么要这么做呢？主要就是为了符合岩土工程的习惯需要，例如桩载荷试验曲线,沉降曲线等都不就是画在第一象限内、最近瞧有很多对后处理中的收敛有很多不同的想法，自己也总结一下plotskdisscale0、07显示位移矢量箭头的后面的scale默认的就是0、053越大箭头越大magf也就是放大变形的意思，但就是不能放在skdis的后面，计算机会出错回到主题，收敛数据,我们想得到一条收敛量随时间变化，或者不同应力直间关系的图线方法1用setlogon然后把数据导入excel或者其她的处理软件进行处理2、用histid=1（不写id的话,电脑会自己按123456的顺序进行排序）gp（点，如果就是体的话用zone）zdisxyz（xyz点处z向位移）这里应力应变之类的都就是用zone，位移之类的都就是用gp最后相要出相应的结果用命令plothist1;如果相要的不同的id之间的关系（比如说竖向应力szz与剪应力sxz的关系）命令plothist2v3（2,3就是对应应力开始设的id号）flac3d程序里面自动设置的步长就是10，也就就是说每隔10步记录一次,这样如果模型小得出的结果就数据少图线不光滑，如果模型大，得到的数据多计算文件又太大，浪费空间。设置步长的命令sethist_repn（n就就是步长）上面说的plot多就是flac里显示的如果想把结果存到别处去用命令histwrite2v3（其实2相当于x轴,3相当于y轴）设置影片格式上面见过，温故而知新setmovieavistep1（存储步长，可以任意设，就就是每隔几步存储的意思）moviestart分步开挖就就是分布计算每多一步（多添一块土，或多开挖，多掘进等）就计算一下,加一个solve即可如果想保存每一步的结果与图，那就在每一步的solve后面添加保存命令设立相应的文件，以后调用即可1、设定初始地应力场就是十分必要的没有地应力场进行计算会使模型在重力作用下发生过大的变形从而不能与实际情况相符2oapply就是施加的外力ini就是初始的内力。3、有人问如果改参数这样做的话地应力不就是被人为改变了不？为什么要人为加大它的参数呢，很不合理？回答:不能改变bulk与shear模量，因为这两个参数就是用来计算应力的。而con与fric就是用来判断破坏状态的。如果将bulk与shear模量改变了，那您在重力条件下计算的应力就不真实了。咱们改大参数的目的就是让模型在不破坏的条件下得到正确的应力，所以只有将con与fric变大才能达到目的。当然如果在弹性模型下计算初始应力就没有这种问题了。因为弹性模型就是没有破坏准侧的当然也没有描述达到破坏状态的参数了。4、大家正好讨论地应力的问题，我也有个疑问,平衡地应力后，水平方向的地应力如何确定的？就是由本构模型来确定的不？这个问题现在我也没答案也就是我的问题5、genzonebricksize10,10,10pl20,0,0p20,20,0p30,0,20modelmohrpropbulk5e9shear3e9fric35inidensity2500setgravity0,0,-10fixxyzrangez-0、10、1iniszz=-5、0e6grad0,0,2、5e4inisxx=-2、5e6grad0,0,1、25e4inisyy=-2、5e6grad0,0,1、25eapplyszz=-4、5e6rangez19、920、1applyszz=-5、0e6rangez-0、10、17applysxx=-2、5e6grad0,0,1、25e4rangex-0、10、1applysxx=-2、5e6grad0,0,1、25e4rangex19、920、1applysyy=-2、5e6grad0,0,1、25e4rangey-0、10、1applysyy=-2、5e6grad0,0,1、25e4rangey19、920、1U#P!y-P$_+D它的意思就是通过iniszz来赋予初始地应力，但就是为何又有setgravity0,0,-10这样不就是重复了不？ini就是的确就是模型的内力,这个内力通过setgrav，然后就像一些帖子上介绍的那样，让程序计算不就行了不?干嘛还有手动在设定一次呢？我感觉它这样的做法有点重复了、回答：(1)它这样做其实只就是减少程序的计算步骤、5_这样设置以后整个模型就已经就是平衡了(2)设置重力的原因不就是为了平衡减少程序计算步骤,而就是因为如果未初始化应力的话，在setgrav的时候,模型会在自重应力下发生应力变形。所以，初始化应力的原因就是防止模型在自重应力下发生变形(我们认为地层自重变形已经发生，当然欠固结等现象就是另外一回事)。3)其实就让模型发生变形，然后把位移清零也可以的6。ini命令来对地层进行初始应力的施加，而APPLY就是~7u(施加应力边界条件7。用到ini初始平衡的时候要与重力计算的相同才有用8。genzonebricksize101010&p00,-25,0p120,-25,0p20,0,0p30,-