Flac使用基础知识
1.sxx是指x方向的正应力，而szz是指z方向的正应力
2.gp_head 结点指针循环，zone_head单元指针循环
4.apply施加边界条件，initial 施加初始条件。
5.dim就是dimension，尺寸。一般指内部尺寸，比如radcyl内部的隧道的尺寸。
<是表示法向量, dist是interface的厚度,norm是表示法向量与X、Y、Z 交角的余旋
programme用法7.检测某点的最大主应力和最小主应力:hist zone sm ax(smin) id …
8.apply sxx 1.0 hist x_stress就是把x_stress的历史记录当成一个力施加给xstress，hist x_stress前面的１表示１倍
9.各点变形量用文件形式输出
set log on
set
set log off
10.显示塑性区
plo bl sta she-n 当前处于剪切破坏
plo bl sta she-p 当前处于弹性，以前处于剪切破坏
plo bl sta ten-n 当前处于抗拉破坏
plo bl sta ten-p 当前处于弹性，以前处于抗拉破坏
这跟flac3d的运算原理有关，它实际上是一个平衡计算扩散的求解过程。与有限元的求解不同：有限元的计算是先组成总体的刚度矩阵，也就是模型有任何一个扰动，模型计算都要进行整体的应力平衡，这样很费内存，也是所有隐式计算程序都使用的方法，这不太符合实际岩体或土体的应力传播实际。而flac3d软件是采用显式计算方法进行的编程，不用形成总体刚度矩阵，节省内存用量。模型中的应力、位移传播、平衡过程比较符合工程实际。以前处于塑性状态实际上是计算过程中（模型中的应力、位移传播、平衡过程中）局部平衡过程中出现的塑性状态。在不断扩大的计算求解中可能该部位又一次调整为了弹性状态，也就是现在处于弹性状态，不过展示塑性区时也要算上该区域！
11.id是指在整个结构中的编号，而cid是指在某一类比如说cable中的编号。拿cable 中的一个单元来说，它既有自己在整个结构中的cd，又有自己在cable 中的cid。一个结构单元类型对应一个ID，该类型单元的具体的每个单元又对应一个子ID，就是CID.打个比方，把人分为两类，男人和女人，男人的ID为1，女人的ID为2，在男人中有10个人，他们的CID就分别为1,2到10，在女人中又有5人，她们的CID分别为1,2,3到5.
12.plot fos显示安全系数
13.求解安全系数并保存到文件名为slope3dfos的文件中
solve fos file slope3dfos.sav associated ,FOS只能用M-M模型
solve fos
默认情况下只折减土体参数
solve fos include interface
可以对interface参数进行折减
13.ini state 0指的是将模型中的塑性区清0，也就是强制使计算模型中不出现塑性状态。
14.一般的岩土问题用set small（默认），特殊情况下的土体用set large，进行模型的塑性计算等大变形计算时设置SET large。也就是允许模型发生诸如差分网格的扭曲等大的变形计算。
15.显示施加在模型上的力和边界条件
pl fap
pl gpfix
对于每个zone多分配10个单位的内存这样可以保存一些自定义的和单元相关的数据
并通过fish进行访问如zextra_：可以按照zone的指针索引你自己定义的数据。grextra——：可以按照grid point 的指针索引你自己定义的数据。
17.应力分布
plot cont smax
确定应力分布方向
principal stress tensor主应力的方向还是得根据弹塑形力学计算得出
18. set mech step step_lim force 50
solve step 1000 force 50
the min. unbalanced force you set.If the min. unbalace force small than this value ,the programme will stop.As a result the programme think the system has reached the status of "balance".
step_lim: This may be the variable defined by the "Fish".If the max step has reached the value "step_lim",the programme will stop without consideration of the status of the system.
solve step 1000 force: max. step = 1000, min. unbalance force=50.  19.ini(或apply) zvel -0.5e-5 range x -.1 3.1 z 9.9 10.1
用apply和ini的区别：
边界条件和初始条件，边界条件不随时步变化，初始条件会变
apply的话，那些位置的zvel就不会变了。
20.mag fac 是变形系数
21.interf 1 face ran plane ori (4,0,0) nor (-5,0,2) dist 0.01 z (0,5) y (1,6)ran表示range，限定命令发生作用的范围
plane ori (4,0,0) nor (-5,0,2) dist 0.01 z (0,5) y (1,6)
表示在通过点(4,0,0)，法向(-5,0,2) 的平面，沿法向 0.01米的范围之内
同时又在 0<z<5，1<y <6范围之内的面
22.solve age=25的意思是进行25个时间单位的动力计算，一般单位是s。
23.把图形保存下来
set plot bitmap（还可以是其它的图形格式，flac里有四种选择）
plot hard 'base'（自己定义的图形显示窗口的名字，默认为base） file 1.bmp
24.plot hist 3 vs -1 显示历史曲线，其中历史3为横坐标，1为纵坐标，其中历史1的数值取负号（即原来为正值的变成负值）
25.dip 就是一个倾斜平面内的倾斜线和负z轴的夹角（锐角）；
dd 就是一个倾斜平面和水平面（x-y面）相交的水平线和正y轴的夹角；dip
是倾角，dd是倾向。dip表示对称参照面与xy平面的夹角，对称参照面与xy
平面的夹角在xy平面的投影是一条射线， dd表示y轴正向顺时针到那条射线的夹角,dip认为是平面从xoy平面向z的负方向旋转角度，dd是平面的法向在xoy上的投影与y轴正向的夹角，注意是沿顺时针旋转。
25.table相当于数已经赋值的数组
26.设置边长gen zone ) size ... p0 p1 edge 100
27.显示shell单元plot add sel geom cid=on node=off ..... 厚度t用thickness表现
< surface生成表面，显示模型中用到
29.显示单元数pr info
30.开挖和开挖后的回填model null range ..... 、model mohr range .....
31.n是通过节点面的外法向
s是平行于这两个面的交线方向d是这个面上梯度最大的方向
32.显示模型中的节点号和单元号:
plot modi i gpnum on
plot modi i node on id on
plot sel geom
plot modi i id on --zone
33.输出实体单元一点的应力
pr zone stress range …… 单元号也可以,
限定x y z 也可以
34.得到组合力 plot sel recover
就是弯矩，薄膜力得等
壳类求解顺序：先得到薄膜力，再由此得到壳类应力
35.ini state 0初始化塑性变形
37.ini z mul .1中mul是个乘子，这个命令可以用来移动网格，此句表示就是说移动10%看原来的z坐标情况决定向那个方向移动。
38.pile表示结构类的桩单元。
39.plot surf yellow sel geom  velo white将几何体的速度矢量显示为白
40.print sel pos得到节点的坐标
41.fix x range x -0.1 0.1 是将坐标在-0.1>x>0.1范围内的节点在x方向进行约束
42.设置内存（2D）"C:\Program Files\Itasca\flac500\exe\" giic 300
43.双击鼠标左键，可获得单元的 id 和 position。
44.如果不晓得一个命令后面跟什么关键字，可以在命令后加？可以获得一点help！比如键入plot ？然后回车，看会有什么提示..........
46.F3，提取前命令
ctrl+R：变回原来的状态
ctrl+c：brings up the camera dialog
ctrl+l: cuttingplane
48.直接按delete，图形就缩小,M是图形放大
49.不平衡力是由于数值计算处理中产生的系统内、外力之差，通常用不平衡应力比来表征迭代计算的精度，类似与ansys等求解,FLAC求解时，通过计算每一步的不平衡力，再与前面的状态进行叠加，直到求得的不平衡力小于规定的值. 假如经过8000步的叠代，达到我们容许的精度，我们认为8000步之后就达
到了平衡状态.那么，这之前没有达到平衡状态(如在100步时)的应力或其他量值是否可用?手册中一个例子是在一个给定的速度下叠代到容许的精度，得到了平均剪应力和平均剪应变的曲线.那么我理解在达到平衡状态之前, "曲线的应力应变值都是不可信的". 显然这种理解不正确.敬请各位大侠答疑解惑!
1. 在FLAC3D中，每一个结点都是由8个单元围绕，它们都对结点贡献结点力。在平衡状态时，这些力的代数和几乎等于零（即作用在这个结点的一个面上的力几乎与另一个面上的力平衡）。如果这个不平衡力达到了一个非零常数，则表明该模型正出现破坏和塑性流动。在计算时间步中，整个网格结点的最大不平衡力被实时计算，并做为记录加以保存和图形显示。
2. 常用最大不平衡力和节点内力的比值R来控制运算的阶段.
3. wwiem 说在(100步)时候，这些数字我认为这些过程值，在分析变化的过程中存在的，且可做分析用的，你说的应力应变曲线是最终的结果，还是过程曲线了. 他是逐渐变化的过程，就象岩石压缩曲线样了，那你说这个曲线是没有可信度么.只能说是在R小与某一个值的时候 ，表示节点力平衡或者处于没有加速度的稳定塑性流.这要结合速度的变化了.
如果速度his趋近于零，说明处于绝对平衡；稳定与某一个值，表明达到一个稳定状态. 每一个结点都是由8个单元围绕，它们都对结点贡献结点力。不平衡力就是这些点做的作用力的合力了.
sub zone对gridpiont的内力贡献称之为等效节点力，等效节点力与外力在该节点贡献的集中力之差称为不平衡力
50.kn =EEr/[s (Er-E)]，
ks =GGr/[s (Gr-G)]

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。