FLAC若干问题的解答

把结果中的信息输出到指定文件 分为两种：

1、把单元信息输出到指定文件 set log on

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMO5R(EO4S$s-Sa y rMP#U1B

set log on finame.dat;可以输入文件路径，否则按当前目录处理 print zone

j4}$`H0_wSimW|*O仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentWw3Nj { set log off 通过反复使用该命令，可以把不同信息输出到不同的文件。2、把节点信息输出到指定文件 set log on

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent\^M2k9S%W\SimWe仿真论坛)k_$KW

SimWe仿真论坛8gS+|)~- set log on finame.dat;可以输入文件路径，否则按当前目录处理 print gp set log off

请教一下

print gp

\q0lW里的keyword该如何填。比如想输出某点的zdis 谢谢 print gp disp print gp position

log 文件中有，然后处理，这些论坛里都有的

关于flac3d—内置fish语言精讲 FISH语言是FLAC3D程序的内置编程语言

|Simwe.com|仿真设计有限元虚拟仪器d7AO'm6y SimWe仿真论坛U/t6vm4\\5B因为FLAC3D的最佳操作方式是命令流文件方式

这一点与ANSYS很相似，而FISH就相当于ANSYS的APDL语言。它包括循环、判断等结构。 如果你还用过其它高级语言，那么从形式上讲你也可以把它理解为 子函数。 FISH语言的引入极大的方便了用户进行复杂的程序建模 它不但可以嵌入命令流文件里工作

SimWe仿真论坛` Fg|Simwe.com|仿真设计有限元虚拟仪器

而且还可以引用FLAC3D本身的任何命令 所以说它实现了对FLAC3D的完全控制 引入FISH语言的作法值得其它通用软件效仿。 以上仅是个人的使用体会 也供初学者入门之用。

-Ow-k|Simwe.com|仿真设计有限元虚拟仪器&V仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent)M0z7S-a W0T7]Y Q

请版主与管理员指教

www.simwe.comO6f~[7G

谢谢！

Fish函数增加了以下几种新的特性：

?,k-Vwww.simwe.comr-a 1.增加了fish变量来获取结点、单元和界面变量 2.fish提供了获取结构单元变量的途径

3.休单元和面单元性质目前可以通过单元变量名z_prop(i_z string)和界面单元名i_prop(i_z string)分别加以识别

\j!Q4E9^Z%i

4.fish函数可以获取单元应变和应力速率，还提供了全应变增量张量和应变速率张量

1

5.提供了fish绘图子程序函数够生用户定义的图形内容 6.fish函数已经增加了从文件读、写数据的能力

www.simwe.com8~`-Um*}

FLAC3D是一个强大的软件，但是不得不承认的，它的界面没有ansys好用，商业运行不是很强的说也来凑个热闹，呵呵注意点：

1.fish函数可以嵌套使用;

2.以save命令保存模型时，fish函数和变量也同时保存；

SimWe仿真论坛$T M$No(VSimWe仿真论坛,?-n?N |5G仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM'T2Z*Rv,f/LA,s

3.fish函数不支持缩写，这与flac3d命令不同，另外所有的fish函数或变量不区分大小写，程序同意转化为大写进行编译，当然也可以通过执行set case-sensitivity on来区分大小写；

航仿真分析有限元模拟计算力学航空4.变量或函数名不能以一个数字开头也不能是下列字符：

. , * / ^ = > < # ( ) [ ] @ ; “ '

5.如果用命令set safe on 指定了编译安全模式，则用户调用fish函数时，函数名前必须加@；

SimWe6.如果变量不曾赋值，则系统默认为零（整形），如果赋值，其类型由值的类型决定； 7.fish函数的调用方法：

.可以出项在其他fish函数的单独行中； .可以出现在其他fish函数的表达式中； .出现在flac3d的命令行中； .作为命令set,print,hist的参数。 8.fish函数不支持递归调用;

9.fish函数可以重新定义，其中变量初始化为原先值； 10.数组名不得和已定义fish函数同名（和参变量不同）；

11.在flac3d命令提示下，set命令可以用已存在的参数对另一参数进行赋值，但数组类不支持 如果用命令set safe on 指定了编译安全模式，则用户调用fish函数时，函数名前必须加@； 今天看了手册，从我的理解来看，应该是对于函数名与FLAC3D中的命令相同时，可以在函数名前加@以示区别。比如initial是FLAC命令，而@initial则为用户自定义的函数名，这样在计算的时候不易出错。但前提是必须set safe on。是这样的吗？

什么情况下设置接触面单元？

比如不同土层之间好像不用，土钉锚杆和土之间是否需要呢？桩和土之间呢？

设置接触面单元的优缺点是什么？谢谢高手答复

如果考虑接触界面对问题的影响或接触界面的变化特征时就的设置接触单元。

如桩与土之间时，由于桩和土的刚度相差较大，所以它们之间的接触单元定义后可以减少由刚度差别较大产生的误差影响（如位移）。当然在定义接触单元后降低了刚度的影响也会提高软件的运行速度。 个人看法供参考

接触单元一般设置在刚度相差较大得界面间，如桩土，地连墙与土，由于刚度相差太大，位移不连续，有限元中尤其如此；设置接触可以反映接触面之间得剪切滑移或侵入行为，如节理面，断层，我想；

add，dim的含义

gen zone radcyl p0 15 0 -15.75 p1 add 15 0 0 p2 add 0 6 0 p3 add 0 0 15.75 &

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent-d7H8rSimWe仿真论坛!qPD$X`/eC3n:J+e9y

|Simwe.com|仿真设计有限元虚拟仪器C'~\~1L:XE

SimWe仿真论坛.hO\\-o [

2

size 3 4 8 5 dim 2.25 2.25 2.25 2.25 rat 1 1 1 1.1 fill

|Simwe.com|仿真设计有限元虚拟仪器V'R(x&s%R

其中add是什么，dim又是做什么用的，已经有了size定义尺寸了 add 是在p0点的坐标上对应增加的坐标大小，dim是圆的半径。

g4hj)x0[3T2C仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM7q 建议好好看一下帮助。

gen zone brick p0 10,0,0 p1 add 10,0,0 p2 add 0,10,0 p3 add 0,0,10 size 10,10,10 这条命令相当于在p0坐标的基础上加变量值，上面的命令相当于

:O7SL.Cl;Igen zone brick p0 10,0,0 p1 20,0,0 p2 10,10,0 p3 10,0,10 size 10,10,10

集思广益之五：关于Flac中几个“力”的理解

学习Flac已有一两个月，但对Flac中几个力的理解还是不到位，在手册也没找到相关阐述，特在此向各位讨教一二

。

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM9e t_{\r+H

第一个力，不平衡力即unbalence force，对这个平衡力，接触也相当多，也知道在计算的时候不平衡力越小越好，但是就是不明白这个unbalence force的实际意思或计算表达式，我的理解是：当在模型上施加了边界条件后，模型内初始应力不能够使之处于平衡状态，不平衡力可能就此产生

，但具体的算法就更不知道了

，与之相关还有一个不平衡力比例，即求解的ratio，

还有就是

这个ratio是怎么得到的，节点的最大不平衡力与所以不平衡力的比值？gp_xfunbal(p_gp)，字面意思当为 节点处x向的不平衡力，这个力 怎么来理解？

第二个力——gp_xfapp(p_gp)，这个力意义又是什么呢？我的理解是作用在节点上的x方向的力

结合具体例子：假如是一个基坑模型，我要得到作用在支护结构上的土压力（假设为x方向），那么因该使用上述的哪个 力合适？我的理解是 gp_xfapp(p_gp)，同样的道理，假如有一根基础桩，桩 上有结构物，现在我想得到桩上某一横截面所受的荷载，应该 使用同样的力 进行操作么？

为更好的学习Flac，理解这几个力 还是挺有帮助的，希望大家不吝赐教，踊跃发言，谢谢unbalence force的实际意思或计算表达式

可以参看 theory background部分, 当对系统施加外力时,

所施加的力要逐步 传给 系统的每个节点,在传递过程中, 就存在不平衡力.随着计算的进行,不平衡力逐渐消散,

gp_xfunbal(p_gp)，是 x方向的 不平衡力,其大小应该是逐渐减小的,

gp_xfapp(p_gp)，是 节点上的外力,包括 体积力,这个力应该是 逐渐变大,或者是不变的...

集思广益之四：关于apply

本来准备对某一节点施加一位移，结果发现好像没有这个命令，只能用施加速度来间接施加位移，用apply zvel range……，当达到我想要的位移时，我想将速度归零，使用apply remove，这是时就出现问题了，使用弹性材料时，在速度归零后，开始产生的位移迅速回弹，可以理解，当我采用摩尔库仑材料时，发现似乎apply remove并未产生作用，监测点的位移还是以同样的速度

3

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM;iW[\]6c

发展，命令流如下，大家可以试试看。

SimWe仿真论坛8Fz'Xo;C(m/Inew

gen zone brick size 3 3 3 model mohr

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM+C`'~a'Y

;}wV o(S eprop bulk 8e6 shear 5e6 fric 15 coh 1e4 fix x range x -0.1 0.1 any x 2.9 3.1 any fix y range y -0.1 0.1 any y 2.9 3.1 anyfix z range z -0.1 0.1 hist unbal

] O&e仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluentw2kU/]!SimWe仿真论坛!R1M_'K3B

apply zvel 0.001 range z 2.9 3.1

|Simwe.com|仿真设计有限元虚拟仪器[&P*h;mhist gp zdisp 2 2 3 set largeplo his 2step 100

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMi.Y\\ step 100

Z/[v6[$K仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentH\

apply remove

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent0QM[7_ 不知道这段命令流问题出在哪里，特向各位请教。可能是我用的Flac版本的问题（盗版的嘛），所以我向给位验证一下，不知道给位怎么理解

还有一个问题就是，当我想让位移就固定在我取消速度的瞬间，而不想让它回弹，这时应该怎么做呢？用fix？谢谢 new

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMX!]-P\r-i$V$ja'a

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMk;fx-e,f0E'XF \\gen zone brick size 3 3 3model m

prop bulk 8e6 shear 5e6 fric 25 coh 5e4 dens 2000fix x range x -0.1 0.1 any x 2.9 3.1 anyfix z range z -0.1 0.1

apply zvel 0.001 range z 2.9 3.1 hist unbalset large step 100

www.simwe.com\SimW |Simwe.com|仿真设计有限元虚拟仪器P'_Hmi(C(i;^SimWe仿真论坛*jb3U'yr0p

fix y range y -0.1 0.1 any y 2.9 3.1 any

www.simwe.comC)V|c!r.V8R Z

hist gp zdisp 2 2 3

|Simwe.com|仿真设计有限元虚拟仪器[+b-a

plo his 2

apply remove zvel range z 2.9 3.1 ini xvel 0 yvel 0 zvel 0step 100 pl cont zvelpause key pl cont zdis

以上对于mc模型可以，对于弹性模型不行，仍存在回弹现象，这个倒是挺有意思的，值得深究。。。

2lR*W|Simwe.com|仿真设计有限元虚拟仪器't+Egf#S-i仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM y$y\B.X6{6x

Q-?-^

4

继续。。。

还有 apply remove zvel range z 2.9 3.1 和 apply remove 的效果是不一样的。。为什么？？？

多谢版主给出提示，按手册上说 apply remove就是将先前所有施加的条件全部移除：If ‘ APPLY remove’ isspecified,all applied conditions are removed.

关于回弹，我想弹性材料是可以理解的，就像对弹性材料加载卸载，变形又恢复了，这就郁闷了，如果我想通过位移加载，那就实现不了啊，按照这种理解，mc材料在将速度归零后，也应该有回弹的，下面我来验证一下

，不知道情况如何？

SimWe仿真论坛/]8Rwww.simwe.com~-

还有哪位大侠有位移加载的方法，不妨一起讨论讨论 唉，结果让我失望，那么我理解也就错了

，mc材料根本没有回弹，为什么两种材料性质会

如此不同呢，难道对于mc材料来说，没有回弹？

进一步的实验发现，apply remove zvel根本不起作用，也就是说只要ini zvel 0就行了，那么apply remove意义何在？

|Simwe.com|仿真设计有限元虚拟仪器3n6Ln)k$G

|Simwe.com|仿真设计有限元虚拟仪器7if+Ks9U:r8W关于mc 材料的回弹可以参考下面一段命令流new

SimWe仿真论坛6n&v4a-H

gen zone brick size 3 3 3 model mo

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluente4UL:e6K7lWm

prop bulk 8e6 shear 5e6 fric 15 coh 1e4 dens 1500 fix x range x -0.1 0.1 any x 2.9 3.1 any fix y range y -0.1 0.1 any y 2.9 3.1 any fix z range z -0.1 0.1set grav 0 0 -10hist gp zdisp 2 2 3 ;set large plo hist 2solve

F\\(CS-x|Simwe.com|仿真设计有限元虚拟仪器0t8@仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent9w-_e:] s

;apply zvel 0.000001 range z 2.9 3.1 hist unbal

|Simwe.com|

;step 100 plo his 2

SimWe仿真论坛;Ys+m/d5g6R8k

|Simwe.com|Simwe.com|仿真设计有限元虚拟仪器;apply remove zvel range z 2.9 3.1set grav 0 0 -5;ini xvel 0 yvel 0 zvel 0 step 100

我通过减小z向加速度来实现卸载，回弹现象还是很明显的，图见下帖，

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMw)NJ8m

，

那么到这里当初的那个矛盾还没有解决，即对弹性材料和mc材料所表现出的不同的反应还在思考中

接上帖，加卸载位移曲线

求助：流变模型与屈服准则是怎么对应的呢?

5

FLAC里静力分析的部分,先把本构方程列出来,然后给出了不同的屈服准则,这样才有了M-C,D-P等模型今天看了些流变的内容,只看到了本构方程,难道在计算这个的时候不考虑屈服准则吗?

仿学航空航天,ANSYS,M 仿真分析有限元模拟计算力今天搞清楚这个问题啦，有些脸红，呵呵，一些东西没有搞清楚．

屈服准则一般是塑性的本构关系中才出现的问题，所以，只有塑性模型才给出了屈服准则，flac里Ｍ-C准则就说的是,M-C plastic model, 同样,在流变模型里,只有出现了塑性元件的时候才会考虑屈服准则,所以,有个模型是关于D-P准则的,便指明了这个准则.各位看看下面一段命令流，发现有何问题啊gen zon brick size 10 10 10gro 1 range x 0 5save 1.savexpgrid 2 new

|Simwe.com仿真设计有限元虚拟仪器$f.wF0H |Simwe.com仿真设计有限元虚拟仪器6w\8FO8w(_仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent+JC3M x!Z V\www.simwe.c

gro 2 range gro 1 not

SimWe仿真论坛3^n R0Z仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMU0n+IS!A(f[#q,f

model null range gro 1

res 1.sav

|Simwe.com仿真设计有限元虚拟仪器3R2|model null range gro 2impgrid 2

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM v7B5U9kJ8b3]#g.Q

inter 1 face range x 4.9 5.1

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent z.zD4m!tfix x range x -0.1 0.1 any x 9.9 10.1 any fix y range y -.01 0.1 any y 9.9 10.1 any fix z range z -0.1 0.1 model mo

prop bulk 1.078e7 shear 4.14e6 coh 1.3e4 fric 12.5 dens 1740 set grav 0 0 -10

inter 1 prop kn 1.078e8 shear 4.14e7 coh 1.3e4 fric 12.5 plo inter ns out on

'D%I\step 100

www.simwe.comF0N2^L;H+_3@\\-S$w+{

仿真分析为何结果显示接触面的法向应力 为0呢？model null range gro 2 inter 1 face range x 4.9 5.1;impgrid 2

SimW |Simwe.com仿真设计有限元虚拟仪器E0j-@ 用这种方法删除 group 2 实际上网格都在，只是网格属性消失了， 要分别建两个模型，然后移来移去

建立接触面首先保证这个面是单独的，没有其他节点连接

model null 仅仅是将单元的材料属性设置为null，并不删除节点和单元 del range 删除节点，单元和材料属性严重同意楼上的说法

$|/`

怎么显示出最大、最小、竖直、水平等应力的方向？用箭头表示的，可以显示在模型上或云图上。

仿真最大。最小用 principal stress tensor，在plot item中。

www.simwe.co3至于竖直和水平的，就是竖直和水平的啊 plot stensor

www.simwe.comh I%| x1^#x

显示应力张量，good lucky!

6

FLAC里静力分析的部分,先把本构方程列出来,然后给出了不同的屈服准则,这样才有了M-C,D-P等模型今天看了些流变的内容,只看到了本构方程,难道在计算这个的时候不考虑屈服准则吗?

仿学航空航天,ANSYS,M 仿真分析有限元模拟计算力今天搞清楚这个问题啦，有些脸红，呵呵，一些东西没有搞清楚．

屈服准则一般是塑性的本构关系中才出现的问题，所以，只有塑性模型才给出了屈服准则，flac里Ｍ-C准则就说的是,M-C plastic model, 同样,在流变模型里,只有出现了塑性元件的时候才会考虑屈服准则,所以,有个模型是关于D-P准则的,便指明了这个准则.各位看看下面一段命令流，发现有何问题啊gen zon brick size 10 10 10gro 1 range x 0 5save 1.savexpgrid 2 new

|Simwe.com仿真设计有限元虚拟仪器$f.wF0H |Simwe.com仿真设计有限元虚拟仪器6w\8FO8w(_仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent+JC3M x!Z V\www.simwe.c

gro 2 range gro 1 not

SimWe仿真论坛3^n R0Z仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMU0n+IS!A(f[#q,f

model null range gro 1

res 1.sav

|Simwe.com仿真设计有限元虚拟仪器3R2|model null range gro 2impgrid 2

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM v7B5U9kJ8b3]#g.Q

inter 1 face range x 4.9 5.1

仿真分析有限元模拟计算力学航空航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent z.zD4m!tfix x range x -0.1 0.1 any x 9.9 10.1 any fix y range y -.01 0.1 any y 9.9 10.1 any fix z range z -0.1 0.1 model mo

prop bulk 1.078e7 shear 4.14e6 coh 1.3e4 fric 12.5 dens 1740 set grav 0 0 -10

inter 1 prop kn 1.078e8 shear 4.14e7 coh 1.3e4 fric 12.5 plo inter ns out on

'D%I\step 100

www.simwe.comF0N2^L;H+_3@\\-S$w+{

仿真分析为何结果显示接触面的法向应力 为0呢？model null range gro 2 inter 1 face range x 4.9 5.1;impgrid 2

SimW |Simwe.com仿真设计有限元虚拟仪器E0j-@ 用这种方法删除 group 2 实际上网格都在，只是网格属性消失了， 要分别建两个模型，然后移来移去

建立接触面首先保证这个面是单独的，没有其他节点连接

model null 仅仅是将单元的材料属性设置为null，并不删除节点和单元 del range 删除节点，单元和材料属性严重同意楼上的说法

$|/`

怎么显示出最大、最小、竖直、水平等应力的方向？用箭头表示的，可以显示在模型上或云图上。

仿真最大。最小用 principal stress tensor，在plot item中。

www.simwe.co3至于竖直和水平的，就是竖直和水平的啊 plot stensor

www.simwe.comh I%| x1^#x

显示应力张量，good lucky!

6

本文来源：https://www.bwwdw.com/article/j2rh.html