flac3d常见问题总结.docx
《flac3d常见问题总结.docx》由会员分享,可在线阅读,更多相关《flac3d常见问题总结.docx(13页珍藏版)》请在冰豆网上搜索。
flac3d常见问题总结
flac3d学习有一段时间,过一段时间没学究忘了很多,所以学习最好不要放下来。
对自己总结,也就是学习的一种。
flac3d命令很多,用到的不多,就那么几个,
基础的就不多说了进入正题:
个人认为如果学习flac3d最好把flac的命令进行总结一下。
做模型就好比拼盘:
要地应力,就把地应力的那段命令摘入,要出图,加入plot的相关命令,要出动画,加个动画命令,所以总结好命令后,以后做模型,只要把命令组一条条的对应放入,一个程序就出来了,后期调试一下,改相关参数就可以了。
接下来,我把自己弄好的整理一下,按照我上面的思路进行建模
什么人建模都会有自己的喜好,下面说一下我自己
建模第一步就就是建立一个全新的txt文件(不论命令大小都就是如此),这样的好处就是,模型可以修改
方便快捷。
而且还可以自我查错
建立新的txt后,第一行命令就就是new,加new的目的就是为了防止其她命令的影响(相当于刷新)再执行其她命令前刷新一下。
不用想,接下来就就是建模了,建立模型,不同人有不同的方法,有的人喜欢建完后在对模型开挖最终得到自己想要的模型。
modelnull等;有的人喜欢用各种带缺角的组合成最终的模型。
具体建模瞧各自的喜好,我觉得,不同点就是建模的网格划分有区别,用带缺角的组合模型网格划分更好
最好把常用的集中模型记下来组成各种句子,以后用的时候直接插入
例如:
genzonebrick/radbricksize*,*,*,*ratio*,*,*,*
建模的过程中,或同一个文件中有重复应用的命令,就可以用宏命令了,
宏命令
macro名称‘命令的内容’
例:
macroshiping’plotconzdipoutlineon’
用shiping代替后面的引号里的东西。
在建模中只要有相同的都可以用shiping代替,缩减了书写命令的时间,整体上也能让人瞧的更明白
注意点就是,您命名的符号最好跟falc自带的命令不要重复,我的处理方法就是,名称后带数字
如shiping1
如果您想修改宏对象,改变命令的内容,
在名称上加单引号就就是替换的意思
例
macro’shiping1‘‘新的命令’
这些很多课本里都有讲,不在叙述
很多时候我们要对单一的对象进行操作,例如进行模拟运算后,只想单独的显示一部分的内容,单一的土层,锚杆,等就要用到group
group给对象命名,
例:
在genzonebrick------size----groupname(name随便取,不要跟flac的已有的命令重合就行)这里的意思就就是对这个体进行了命名
接下来想把这一块模型挖掉加
modelnullrangegroupname
group的用法:
1。
group直接加到模型的后面相当于给整个模型命名
2、groupnamerangexxyyzz给一个范围的体命名
继续学习,昨天说到了建模,我就是搞巷道的,重点讲巷道。
巷道建模,用到brick,radbrick,
radtunnel,radcylinder这4个模型比较多,一般模型都就是建一半,让后用镜像命令复制另一半,
建模也就是大家有各自不同的习惯,总的有两种,一种用点法,就就是吧模型的所有关键点都输入进去,这种方法可以用来建立不规则的模型,弊端就就是工作量大,不当在设点的时候考虑多,输入的时候也很麻烦。
第二种。
按照模型只要取三个点确定xyz三周的方法,有点输入简单,弊端,只能建立规则的模型,像要建立,路堤模型,边坡模型,就只能就是按照数点的方法建模
range的命令:
书上有,但就是我用自己的话说一下她的用途
1:
取圆球面的范围,就就是同一个中心,不同半径之间夹的那一部分,用在给球面施加力,约束时用该命令,命令:
rangeannuluscenterxyz(中心坐标)radiusr1r2(两半径)
想要表达的就是球面时,r1r2取值相差小点即可,如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)显示
这一整部分就就是我开头说的要自己整理出来的句段了,等以后用到的时候直接插入,作用就就是显示模型的。
任何模型的显示都可以用她来表示。
很好很强大。
如果还要显示其她功能就把相应的命令在插入上面那部分就行了
plotsetrotation20030
plotsetcenterauto
plotsetmagnification1、5
plotaddsurfaceyellow
plotaddaxes
plotshow
这个就就是我常用来查瞧模型的命令,大家也可以自己整理一分,
总结,做个实例,给大家瞧瞧
隧道,半圆直墙隧道
分两步
一,建模,
二,显示
建模建一半然后镜像对吧
尺寸,自己在草纸上画好,点坐标都弄清楚了
把以前的句子直接套下来
genzoneradcylindersize425614(网格数量,自己瞧着办)ratio1111、2&
p0000p15000p201000p30050dimension222
这就是建上半部分
genzoneradtunnelsize425614ratio1111、2&
p0000p100-50p201000p35000dimension1、822
建立下半部部分
建一半后,用镜像
命令也讲过了直接插入,
genzonereflectnormal-100origin000
最后一步就就是显示给大家瞧
那一打包显示直接插入
plotsetrotation20030
plotsetcenterauto
plotsetmagnification1、5
plotaddsurfaceyellow
plotaddaxes
plotshow
。
整体的建模加显示,的命令就完了
合一块建立一个新的txt开头加一个new
new
genzoneradcylindersize425614ratio1111、2&
p0000p15000p201000p30050dimension222
genzoneradtunnelsize425614ratio1111、2&
p0000p100-50p201000p35000dimension1、822
genzonereflectnormal-100origin000
plotsetrotation20030
plotsetcenterauto
plotsetmagnification1、5
plotaddsurfaceyellow
plotaddaxes
plotshow
整体就这样了,大家可以粘贴一下瞧瞧。
plot命令大集合:
显示切片的命令
plotsetplaneori000norm010
plotconzdplane
这两行命令,用来显示指定面上的位移变化云图。
(吧第二行改改,还可以显示应力云图,等其她)
plotconsyyy向应力云图
plotconszzoutlineonz向应力云图加显示网格
plotconszzeffectiveonz向有效应力
plotconszzinter1e4云图显示增量(后面的值越小,图上颜色条带就越细,)
plotconszzmin-10e4显示最小值,图上小于这个数的值都用同一色表示
plotconszzmax10e4显示最大值,图上大于这个数的值都用同一色表示
plotfipoutsideon显示力的箭头
建模完事后就要给模型赋值,赋参数了,关键词,model
modelmor(赋摩尔库伦模型,想要什么模型改相应的词)
proden2500(密度)bulk2e8体积模量,shear2e8(剪切模量)co2e6(内聚力)fric45(摩擦角)ten1e6抗拉强度。
range、、、、、、、、、、、(我们要赋值的范围)
给模型赋值就就是hi这么简单,难在参数的选取。
现在瞧边界,如果就是不动面,用fixxrange。
。
。
。
。
如果边界带速率,或位移用ini。
。
。
setgrav00,-9、81
inidens1000
设置重力就一定要设置密度
不考虑重力就不要设计密度。
两者共存的,设计程序的时候一定要记的
setmechforce50
solve
用SET命令设定了最大不平衡力的上限值为50N,即:
计算的最大不平衡力一旦小于此值,我们便人为模型达到了平衡状态。
求解过程终止。
设置这个目的就是在一定精度要求下,加快计算(值越大,算的量越少,花的时间也越少)
说到这里,add与set的区别
add原来有的加新的用add
set原来没有的,设置用set
plotcreateGravV
plotsetplanedip=90dd=0origin=3,4,0
plotsetrot15020
plotsetcenter2、54、24、0
plotaddboundbehind
plotaddbcontszzplane
plotaddaxes
plotshow。
如这个命令,显示图像的时候命令有很多,
那么我们想显示别的图像,又不想吧这个弄好的图像丢掉
则用plotcreatename建立一个图像,在建立新的图像时
plotcreatename1(新的名)
显示以前的图像
那么输入如下命令:
plotcurrentname。
另一种方法,输入:
plotshowname
这里有一点最好开始的时候将网格剖分的相对稀疏一些(1000~1500个单元
左右),先求解一下,瞧瞧结果,然后再根据大致的结果有针对性的增加模型的
网格数,以得到更精确的解答。
更快的检验一下自己的模型就是否建错了
不同的区域定义不同的本构模型
modelmohrrangez=0,5
propClayeyGravelrangez=0,5
在定义参数后面加个range范围
为了避免在输入命令时总要通过range+全局坐标域的方法来限定该命令执
行的区域,我们可以先通过range命令来给不同的区域命名,然后就可以用range
+区域名的方法来限定范围了。
这样就避免了每次输入命令要限定在该区域执行
时,都要输入一遍坐标范围的麻烦。
具体如何给区域命名,以及日后如何用区域
名来限定该区域,可以瞧下面这个例子:
rangenametrenchx=-1,1y=0,4z=-2,2
modelnullrangetrench
今天做方案,没时间,就讲几点得了,1:
右手法则,建立模型时,的点的排放时有规律的,
方法如下,拿出右手
伸直了
4指合并指向就就是p0到p1的方向
掌心指向就就是p0到p2的方向
拇指与4指垂直的方向就就是p0到p3的方向
很多人说网格后的云图有带齿什么的,改进的方法就就是重复用下面的命令
GenADJust
用此命令来微调网格,使之离散平滑。
该命令可以连续使用以增进平滑效果的。
ploblogro可以瞧到不同的group的颜色分布(例不同的土层定义组后,就可以用此命令显示不同的色了)
如何查瞧剖面上任一点的位移值?
printgpdisprangex()y()z()或
printgpdisprangeid
plotgpdisprangex()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-L
2)"-"表示对其值作"mirror",如果大家使用AUTOCAD的话就明白"mirror"的含义了,比如上面的hist1gpydisp的值就是0~100,那么vs-1就变成-100~0、以此类推、
3)为什么要这么做呢?
主要就是为了符合岩土工程的习惯需要,例如桩载荷试验曲线,沉降曲线等都不就是画在第一象限内、
最近瞧有很多对后处理中的收敛有很多不同的想法,
自己也总结一下
plotskdisscale0、07显示位移矢量箭头的后面的scale默认的就是0、05
3越大箭头越大
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、设定初始地应力场就是十分必要的没有地应力场进行计算会使模型在重力作用下发生过大的变形从而不能与实际情况相符
2。
apply就是施加的外力ini就是初始的内力。
3、有人问如果改参数这样做的话地应力不就是被人为改变了不?
为什么要人为加大它的参数呢,很不合理?
回答:
不能改变bulk与shear模量,因为这两个参数就是用来计算应力的。
而con与fric就是用来判断破坏状态的。
如果将bulk与shear模量改变了,那您在重力条件下计算的应力就不真实了。
咱们改大参数的目的就是让模型在不破坏的条件下得到正确的应力,所以只有将con与fric变大才能达到目的。
当然如果在弹性模型下计算初始应力就没有这种问题了。
因为弹性模型就是没有破坏准侧的当然也没有描述达到破坏状态的参数了。
4、大家正好讨论地应力的问题,我也有个疑问,平衡地应力后,水平方向的地应力如何确定的?
就是由本构模型来确定的不?
这个问题现在我也没答案也就是我的问题
5、genzonebricksize10,10,10p120,0,0p20,20,0p30,0,20
modelmohr
propbulk5e9shear3e9fric35
inidensity2500
setgravity0,0,-10
fixxyzrangez-0、10、1
iniszz=-5、0e6grad0,0,2、5e4
inisxx=-2、5e6grad0,0,1、25e4
inisyy=-2、5e6grad0,0,1、25e
applyszz=-4、5e6rangez19、920、1
applyszz=-5、0e6rangez-0、10、17
applysxx=-2、5e6grad0,0,1、25e4rangex-0、10、1
applysxx=-2、5e6grad0,0,1、25e4rangex19、920、1
applysyy=-2、5e6grad0,0,1、25e4rangey-0、10、1
applysyy=-2、5e6grad0,0,1、25e4rangey19、920、1
U#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,-25,20
modelelas
propbulk5e9shear3e9
inidensity1600rangey-10,0
inidensity2000rangey-15,-10
inidensity2200rangey-25,-15
setgravity0,-10
fixxrangex-、1、1
fixxrangex19、920、1
fixzrangez-、1、1
fixzrangez19、920、1
fixyrangey-25、1-24、9
solve
如果这样就是不就是就完成了初应力的计算了?
当然在这个基础上进一步计算需初始位移为0:
如果加入下一段代码来计算初应力:
inisyy0、0grad0,1、6e4,0rangey-10,0
inisyy4e4grad0,2、0e4,0rangey-15,-10
inisyy7e4grad0,2、2e4,0rangey-15,-25
以上添加了力边界,在这种情况下,得到的初应力就是不就是与前一个方法得到的相同呢?
另外,这样做就是不就是就可以不初始(置零)位移呢?
我总感觉这个initial命令象流体软件fluent里面求解前设置初始条件,它的最重要的作用就是帮助解收敛的快一些,
a)不知道我的理解对不对,斑竹及各位大哥帮帮忙吧,我这个可怜的新手,谢谢!
回答:
(1)我试着计算了这两种情况,对于ini应力的作用,我得出以下结论:
A、对于差分法,为了加快收敛速度,在开始求解之前,最好能设置初始条件(如fluent