flac3d常用命令
1、最先需要掌握的命令有哪些?答:需要掌握 ge n, in i, ap p, pl o, so lve 等建模、初始条件、边界条件、后处理和求解的命令。
2、怎样输出模型的后处理图?答:Fil e/Pri nt ty pe/Jp g fil e,然后选择 Fil e/Pri nt,将保存格式选择为 jpe文件。矢量图是什么意思
3、怎样调用一个文件?答:F ile/c all 或者 cal l 命令
4、如何施加面力?答:app nstr ess
5、如何调整视图的大小、角度?答:综合使用 x,y, z, m, S hift键,配合使用
Ctr l+R,C trl+Z等快捷键。 6、如何进行边界约束?答:fixx ran(约束的是速度,在初始情况下约束等效于位移约束)。 7、如何知道每个单元的ID?答:用鼠标双击单元的表面,可以知道单元的ID 和坐标。 8、如何进行切片?答:
plo s et pl ane o ri (点坐标) n orm (法向矢量) plocon s z pla ne (显示 z 方向应力的切片) 9、如何保存计算结果?答:sav e +文件名 10、如何调用已保存的结果?
答:r est +文件名;或者 Fil e / R estor 11、如何暂停计算?答:E sc 12、如何在程序中进行暂停,并可恢复计算?答:在命令中加入p ause命令,用conti nue 进行继续。在我们分步求解中想得到某一个过程中的结果,不用等到全求完,还可以在分布求解错误的时候就进行改正,而不是等到结果出来。13、如何跳过某个计算步?答:在计算中按空格键跳过本次计算,自动进入下一步14、Fis h 是什么东西?Fi sh 是否一定要学?答:是FLAC3D 的内置语言,可以用来进行参数化模
型、完成命令本身不能进行的功能。 Fish可以不用学,需要的时候查Mannu al 获得需要的变量就可以了。 15、F LAC3D允许的命令文件格式有哪些?答:无所谓,只要是文本文件,什么后缀都可以。 16、如何调用一些可选模块?答:
confi g dyn (flu id, c reep, cppu dm) 17、如何在圆柱体四周如何施加约束条件?可以
用f ix ... ran cyli nderend1end2radiu s r1cylin der e nd1 e nd2 r adius r2 n ot,其中 r219、关于 dd、d ip 的理解: di p-dir ectio n ang le, m easur ed in thegloba l xy-plane
cloc kwise from theposit ive y-axis这个是从 y 轴到 x 顺时针方向所转过的角度; dipdip a ngle, meas uredin th e neg ative z-di recti on fr om th e glo bal x y-pla ne 这个是 xy平面与z轴负方向的夹角;dip 表示对称参照面与 xy平面的夹角,对称参照面与x y 平面的夹角在x y 平面的投影是一条射线;d d 指的是从 y 轴正方向按顺时针
(cl ockwi se)方向转向所要确定面的法线方向在xy 平面上的投影的夹角,然后再确定dip,指的是从xy 平面转向所要dip 确定的平面的角。 20、rang e 的命令:
①取圆球面的范围,就是同一个中心,不同半径之间夹的那一部分,用在给球面施加力,约束时用该命令,命令:r angeannul us ce nterxyz(中心坐标)r adius r1 r2(两半径),想要表达的是球面时,r1r2 取值相差小点即可,如radiu s 2.9 3.1这时范围就是半径为3 的面。②圆柱的环面施加力或约束时也用到该命令。 ran ge
cy linde r end1 xyz end2 xyzradiu s r(e nd 是表示柱体断面的中心点坐标) ra nge p laneabove表示的是 z 大于 0 的上部结构r angeplane blow表示的是 z 小于 0 的下部结
构r angedd adip b orig in xy z 表示的是走向角度 a,倾向角度 b的面 dd表示,当面与 xo y 面相交时的垂线与 y 轴正向的夹角,顺为正,说白话,就是面的指向方向面是朝那个方向的意思。d ip 面与 xoy面的夹角,与 z 轴负向为正,origi n 表示面经过 xy z 点;range norm al xy z 法线方向为 xy z 的面range not没试用过(范围外的东西,想想用到圆上取范围非常好)。 21、如何调用 FISH文件?新手调用FISH文件出错的原因:是:文件名为fishc all.f is.tx t(因为在记事本里这个 tx t 的后缀是默
认的)。修改方法是:将fishc all.f is.tx t 里面的内容全部复制到新的记事本里,再按保存,保存时候将保存类型改为所有文件(*,*),这时你便可在文件名栏内输入
fish call.fis,类型也就正确了。 22、如何理解 Tens ion-p shea r-p t ensio n-n s hear-n 表示的含义?在计算的循环里面,每个循环中,每个 zo ne 都依据 fai lurecrite rion处于不同的状态,s hear和 ten sion分别表示因受剪和受拉而处于塑性状态;n表示
no w, p表示 pr eviou s,即分别表示在现在和以前的循环当中处于塑性状态。 23、绘制相关曲线命令的一些解释。1) 在p lot h ist m vs n的形式里,m 代表 y 轴,n 代表x 轴(不管 m,n的正负); 2)"-"表示对其值作"mirro r"镜像。比如上面的 hist 1 gp
ydis p 的值是 0~100,那么 vs -1 就变成-100~0.以此类推; 3)为什么要这么做呢?主要是为了符合岩土工程的习惯需要,例如桩载荷试验曲线,沉降曲线等都不是画在第一象限内。24、如何查看剖面上任一点的位移值?p rintgp di sp ra nge x() y( ) z( )或p rintgp di sp ra nge i d plo t gpdisprange x()y( )z( )或 plot gp d isp r angeid
25、问 a式与 b式有何区别? fix x yz ran ge x-0.10.1 ..........................a fix x ra nge x -0.1 0.1................................b a表示在x=0 这个平面是固定的;b 表示在 x=0这个平面只 x 方向不能发生运动。26、ini tial与 app ly 有何区别?i nitia l 初始化命令,如初始化计算体的应力状态等; ap ply 边界条件限制命令,如施加边界的力、位移等约束等;i nitia l 的应力状态会随计算过程的发生而发生改变,一般体力需要初始化,而 appl y 施加的边界条件不会发生变化。 27、如何查看各个时段不平衡力的具体
数值?采用 his来记录计算,包括位移应力等命令 hi s unb al;h is gp(zone) zdi s ran ge (0 0 0)或者i d=?;导出数据命令 his writ e n v s m b egin时步 en d 时步+ file file name.his n表示纪录的 idm 表示时步;要导出不平衡力的具体数值 hisunbal step 100000 or solv e his writ e 1 v s ste p beg in 1end 1000 f ile 123.hi s 使用上述命令就可以查看各个时步下的不平衡力的具体数值。28、如何显示变形轮廓线的命令? plo skemagf10 其中 10 为放大系数29、怎样查看模型的矢量图?pl gr id ve l (di sp) r ed pl o blo gro可以看到不同的 gr oup 的颜分布(看模型的样子);plo g pfixred s k 看模型的边界情况 plofap r ed sk看模型的体力分布plo c on di s (xd is, y dis,zdis)看模型的云图:位移 plocon s z (sy, sx, sxy, syz, sxz)看模型的云图:应力 plodis (xdis, ydis, zdi s) 看模型的矢量图 ploinfo看模型有多少单元、节点 plo t gri d 可以查看网格,p lot g rid n um 可以查看节点号。 plo t fix查看边界条件 pl ot gr id ma rk 查看预先留设的开挖边界pl sy y (sx x xdi s ydi s sig1..)fill;查看模型的云图。注意,带 fi ll 关键字就是查看云图,不带 fil l 关键字就是查看等直线图。plo a dd sk e 剖面的沉降云图上增加网格plo a dd di s pla ne 剖面上的矢量plo a dd ax e 剖面上坐标值30、怎样查看 INT ERFAC E?答:pl gr id if ace i d red prin t ifa ce id ; id是 in terfa ce 的编号 31、FLAC怎么调整视图?答: ctrl+z 或shift+z ,然后点击要放大部位的左上角和右下角就可以放大某个
部位,如果要缩小,双击屏幕。32、判断模型中网格点是一个还是两个的方法。显示模型网格点的 ID,如果有两个ID 号,说明是两个。节理面上的网格点应该有两个 ID 号。 33、改变模型应力边界时FLAC
与 3DE C 的不同之处。在 3DEC 中,当改变应力边界时,模型表面的应力直接被更换,FLAC中改变应力边界
时,模型表面的应力为两者应力相加。 34、f lac 中影响运行时间的因素(1)单元材料,结构面之间的刚度对比很大;(2)划分单元尺寸对比很大,尽量避免长细
比大于 5:1的细长单元。 35、FLAC中有多少中材料模型。(1)零模型:用于开挖回填的区域(2)各向同性弹性模型(3)正交各向异性弹性模型(4)横观各向同性弹性模型(5)德鲁克-普拉格塑性模型(6)莫尔-库仑塑性模型,一般岩土体的力学行为(7)节理化塑性模型:用于研究材料破坏后行为(如边破稳定,屈服煤矿开采)(8)应变硬化/软化模型-库仑塑性模型:薄板层状材料破坏后力学行为(9)双线性应变硬化/软化节理化塑性模型临近层状土层开挖
(10)双屈服塑性模型(11)修正的剑桥粘土模型。 36、对不平衡力的理解。 flac 中用最大不平衡力和典型内力大的比值R表示模型的不平衡力的相对大小,该值为百分数,且从不会减小到零,因此,模型不会达到绝对稳定的平横状态,可以根据要求的精度设定R 值,模型默认的R 值是1e-5。注意 R值小仅表示所有网格点的力都平衡,但模型还可能发生塑性流动,通过plotvel 命令查看网格点速
度,一个低振幅随即速度场是没有塑性流动的平衡状态的可靠的指示器。 37、简单命令解释: plo t add cont disp outon 时时显示云图plotblock grou p gpn um on显示节点编号 pl ot in ter r ed sk e 看接触
面的效果. delrange grou p 1 n ot 不删除 gro up 信息delrange grou p 2 删除 gro up 信息 iniz add数值大小 rang e gro up na me 顺z 方向移动建立接触面 in terfa ce 1facerange x(范围) y(范围) z(范围)或者 int erfac e
1 f ace r angeplane orig(x yz 坐标)norm(法线方向)water tabl e cle ar 清除水位线p l set back whit e 设置背景 set plot bitm ap si ze (800,510) 设置图片大小plotset c aptio n
siz e 37设置图例数字大小(defau lt=35范围10-50) plot setcente r 100 0 40 fixx y z rang e x -0.1 0.1 表示在 x=0这个平面是固定的fixed boun daryfix x rang e x -0.1 0.1 表示在 x=0这个平面只 x 方向不能发生运动r oller boun dary。 38、在固定面上施加
位移时,fix 与 appl y 不能同时使用,应该用 in i 和f ix, 才能施加上。 39、你可以用键盘上的上、下、左、右肩在窗口中移动模型。40、在进行三轴压缩时候,如何施加围压?答:使用 appl y sxx =-1e5 szz=-1e5 rang e cyl end1 0,0,0 end2 0,2,0 ra dius1,其
中 X 和Z 轴是在我们眼睛正对的面。42、在进行三轴压缩时候,如何对轴向力进行位移控制?答:使用 ini yvel 1e-7 rang e y -.1 .1 iniyvel1e-7range y 1.9 2.1即可,数值可以变化。 43、在进行边坡计算的时候,为什么用 sol, fos命令运行?答:应该输入s ol fo s slo pe3df os.sa v ass ociat ed(不要也行);再用 re store命令调用即可。44、在FLAC3D 计算中,如何记录程序中的值,如果不想时时记录该怎么弄?
答:用h ist w titen1 vs n2 f ile.t xt,想间隔记录可用 set hist_repn 45、想记录接触面上的相关力学参数该怎么办?答:his t int er instre ss /n dis/s stres s/sdi sp 46、接触面设置参数包括哪几个?答:coh dial fric kn k s ten. 47、如何输出单元格的数据?答;e xpgri d .fl acsd48、如何做切片?答:plo t set plan e ori x yz nor m x y z pl ot co n zdplane这两行命令,用来显示指定面上的位移变化云图(把第二行改改,还可以显示应力云图,等其他); pl ot co n syy y 向应力云图p lot c on sz z out lineon z向应力云图加显示网格 plot conszz e ffect ive o n z 向有效应力plotcon s zz in ter 1e4 云图显示增量(后面的值越小,图上颜条带就越细,) plo t con szzmin -10e4显示最小值,图上小于这个数的值都用同一表示 pl ot co n szz max10e4显示最大值,图上大于这个数的值都用同一表示 pl ot fi p out sideon 显示力的箭头现在看边界,如果是不动面,fix x rang e…如果边界带速率,用或位移用 in i… se t gra v
0 0,-9.81 ini dens 1000设置重力就一定要设置密度,不考虑重力就不要设计密度。两者共存的,设计程序的时候一定要记的。 se t mec h for ce 50 solv e 用S ET 命令设定了最大不平衡力的上限值为50N,即:计算的最大不平衡力一旦小于此值,我们便人为模型达到了平衡状态。求解过程终止。设置这个目的是在一定精度要求下,加快计算(值越大,算的量越少,花的时间也越少)。 49、FLAC结果数据的导出p l set back whit e 设置背景 set plot bitm ap si ze (800,510) 设置图片大小plotset c aptio n siz e 37设置图例数字大小(defau lt=35范围10-50) main win s ize 1 1
po sitio n 0 0.8 ;comma nd wi ndowplotset w indow size 0.80.8 p ositi on 0.1 0.0 ; gr aphic s win dow p
lot s et ce nter100 0 40 根据情况更改。
50、一段 fi sh 命令的解释d ef au x1 定义函数 lo op n(1,npt ab1)开始一个循环,从1
循环到变量 npt ab1 r r=xta ble(1,n) 从 tabl e 1 的 x 列中读第 n个数付给变量 rr
zz=yt able(1,n)从 tab le 1的 y 列中读第n个数付给变量 zz loop k(1,nprof) 开始第二个循环,从 1 循环到变量nprof case_of k下面是一个选择结构,利用变量 k 判断alfa=0.5*(pi+.1)*fl oat(k-2)/f loat(nprof-2) x table(n+10,k)=r r*cos(alfa)
yta ble(n+10,k)=rr*sin(a lfa)case1 当k=1 的时候执行下面语句 xt able(n+10,k)=rr
tabl e n+10 的x列第k个值等于变量 rr ytab le(n+10,k)=-40.0 tab le n+10 的y 列第k 个值等于变量-40 cas e 2 当 k=2的时候执行下面语句xtabl e(n+10,k)=rr ta ble n+10 的x 列第 k 个值等于变量rr yt able(n+10,k)=0.0 tab le n+10 的y 列第k 个值等于变量0 end_case结束选择结构 end_loop结束第二个循环e nd_lo op 结束第一个循环 end函数定义结束 aux1 执行什么定义的函数 51、有关命令的解释 1)apply nvel 0.1plane dip60 dd 270range left_boun appl y 施加边界条件,p lane法向速度为 0.1,plane面的位置:倾角60(在x oy 平面从-Z 轴方向逆时针旋转),dd 走向 270度,从+y 方向绕z 轴顺时针旋转;2)ini szz-5e4grad0 0 -1e4 i ni 初始边界条件,szz=-
5e4+x*0+y*0+-z*-1e4,注:-5e4 是r ange范围起点坐标的S ZZ 值,x,y,z所施加应力网格点的坐标。52、FLA C 程式使用前准备步骤步骤1:依比例画出所欲分析的资料在纸上画出地点的位置、地层资料、并简单标示距离及深度资料。步骤 2:换算输入资料成同一单位需谨慎检查输入资料的单位,如因单位不同而造成过大或过小的值,将会造成 FL AC 无法计算,而产生 ERR OR 信息。步骤3:应用公式简略计算应用公式或依据经验,简略算出 FLAC输出资料的范围,以作为
De bug 及输出资料分析时验证。步骤4:建立x, y坐标与n ode i, j 之间的关系在图上距离及深度的关系,建立x, y坐标系统,再由 x, y 坐标系统,转换与网格关系,为了便于以后输出资料的分析,故应确实掌握网格的位置及其相对应的 x, y 坐标。建议在敏感区域使用较密的网格,其它地方则使用较疏的网格,刚开始执行程序时,不宜使用网格太大的网格数目,因尽量使网格总数少于 1000,以节省时间。
53、F LAC 3D 程序的编写步骤: (1)建立 x,y 坐标与网格的关系,建议使用 Gen指示:G en x1,y1 x2,y2,x3,y3 x4,y4 i=i0,i1 j=j0,j1 (2)设定材料性质:p rop (3)设定外力:Set Grav, App ly Pr essur e, in i sxx, Syy (4)设定边界条件:fix, free (5)求起始的应力平衡:so lve (6)储存:Save54、建模的过程中,或同一个文件中有重复应用的命令,就可以用宏命令了,宏命令: macr o 名称?命令的内容? 例:macro ship ing ? plot conzdipoutli ne on? 用s hipin g 代替后面的引号里的东西。在建模中只要有相同
的都可以用s hipin g 代替,缩减了书写命令的时间,整体上也能让人看的更明白。注意,你命名的符号
最好跟 fa lc 自带的命令不要重复,处理方法是,名称后带数字。如s hipin g1;如果你想修改宏对象,改变命令的内容,在名称上加单引号就是替换的意思。例 macr o ?sh iping1? ?新的命令?。 55、很多时候我们要对单一的对象进行操作,例如进行模拟运算后,只想单独的显示一部分的内容,单一的土层,锚杆,等就要用到 gro up。g roup给对象命名,例:gen z one b rick------ siz e ---- gro up na me na me 随便取,在(不要跟 flac的已有的命令重合就行)这里的意思就是对这个体进行了命名,接下来想把这一块模型挖掉加: mode l nul l ran ge
gr oup n ame g roup的用法:①group直接加到模型的后面相当于给整个模型命名;
②gr oup n ame r angex xyy zz 给一个范围的体命名。 56、上面说过建模一般只建一半然后用refle ct 命令(镜像)。命令:refle ct no rmalxyz o rigin xyz(根据法线和过一点建的面)最常用至于实际使用过程中有人问道:norm al -1 0 0与 nor mal 1 0 0的区别,我试过两者的效果是一样的,没区别,虽然方向不同,但表达的是同一个面。refle ct dd a di p b o rigin xyz这个就是用里建立斜面镜像的。有一点注意了,镜像的时候如果是模型镜像的话,相当于在建立一个模型,故,在 re flect的前面要加建模的关键字,ge n zon e ref lect…… plo t set rota tion20 030 ;设置观测角度的plo t set cent er au to ;设置观测中心点 plo t set magn ifica tion1.5 ;设置放大倍数的
pl ot ad d sur faceyello w plo t add axes plot show;显示这一整部分就是我开头说的要自己整理出来的
句段了,等以后用到的时候直接插入,作用就是显示模型的。任何模型的显示都可以用他来表示。很好很强大。如果还要显示其他功能就把相应的命令在插入上面那部分就行了。这个就是我常用来查看模型的命令,大家也可以自己整理一份。57、ad d 与s et 的区别 add原来有的加新的用add
s et 原来没有的,设置用 se t plo t cre ate G ravVplotset p lanedip=90 dd=0 ori gin=3,4,0plotset r ot 15 0 20 plot setcente r 2.5 4.24.0 p lot a dd bo und b ehind plot addbcont szzplane plot addaxesplotshow。我们想显示别的图像,又不想把这个弄好的图像丢掉,则用p lot c reate name建立一个图像,在建立新的图像时,p lot c reate name1 (新的名)显示以前的图像,那么输入如下命令:plot curr ent n ame。另一种方
法,输入:pl ot sh ow na me 58、不同的区域定义不同的本构模型 mode l moh r ran ge z=0,5 p rop C layey Grave l ran ge z=0,5 在定义参数后面加个r ange范围,为了避免在输入命令时总要通过 ra nge+全局坐标域的方法来限定该命令执行的区域,我们可以先通过 ra nge 命令来给不同的区域命名,然后就可以用r ange+区域名的方法来限定范围了。这样就避免了每次输入命令要限定在该区域执行时,都要输入一遍坐标范围的麻烦。具体如何给区域命名,以及日后如何用区域名来限定该区域,可以看下面这个例子:range name tren ch x=-1,1y=0,4 z=-2,2 mo del n ull r angetrenc h 59、右手法则:建立模型时,各点的排放时有规律的。方法如下:拿出右手伸直了 4 指合并指向就是 p0到 p1的方向,掌心指向就是 p0到 p2的方向,拇指与 4指垂直的方向就是p0 到p3 的方向。很多人说网格后的云图有带齿什么的,
改进的方法就是重复用下面的命令Gen A DJust用此命令来微调网格, 使之离散平滑。该命令可以连续使用以增进平滑效果的。 60、如果想看模型在计算过程中的变化情况S et mo vie a vi st ep 1file1.avi;定义动画,动画记录的频率为 1 时步一频率(计算量大的时候就没必要1时步的记录,可以10,100 时步的记录,减少电脑的使
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。
发表评论