flac使用基础知识.doc

I. SXX是指X方向的正应力，而SZZ是指Z方向的正应力2. gp_head结点指针循环，zone_head单元指针循环3. grad线性梯度丿V力的关系4. apply施加边界条件，initial施加初始条件5. dim就是dimension,尺寸般指内部尺寸，比女U「adcyl内部的隧道的尺寸6. norm是表示法向量,dist是interface的厚度,norm是表示法向量与X、Y、Z交角的余旋7. 检测某点的最大主应力和最小主应力:hist zone smax(smin) id ...8. apply sxx 1.0 hist x_stress 就是把 x_stress 的丿力史记录"i成一个力施加给 xstress,hist x_stress 前面的1表示1倍9•各点变形量用文件形式输出set log onset logfile gp-disp.txtset log off10. 显示報性区plo bl sta she-n当前处于剪切破坏plo b! sta she・p当前处于弹性，以前处于剪切破坏plo bl sta ten-n当前处于抗拉破坏plo bl sta ten-p沟前处丁•弹性，以前处丁•抗拉破坏这跟flac3d的运算原理有关，它实际上是一个平衡计算扩散的求解过稈。

与有限元的求解不 同：有-限元的计算是先组成总体的刚度矩阵，也就是模型有任何一个扰动，模型计算都要进 行整体的应力平衡，这样很费内存，也是所有隐式计算稈序都使用的方法，这不太符合实际 岩体或土体的应力传播实际而fiac3d软件是采用显式计算方法进行的编程，不用形成总体 刚度矩阵，节省内存用量模型中的应力、位移传播、平衡过程比较符合丁程实际以前处 于犁性状态实际上是计算过程中(模型中的应力、位移传播、平衡过程中)局部平衡过程中 出现的艰性状态在不断扩大的计算求解中可能该部位又一次调整为了弹性状态，也就是现 在处于弹性状态，不过展示舉性区时也要算上该区域！II. 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 in elude interface可以对interface参数进行折减13.ini s⑻eO指的是将模型中的槊性区清0,也就是强制使计算模型中不出现卑•性状态。

14. 一般的岩土问题川set small （默认），特殊情况下的土体用set large,进行模型的塑性计 算等大变形计算时设置SET largeo也就是允许模型发生诸如弟分网格的扭曲等大的变形计 算15. 显示施加在模型上的力和边界条件pl fappl gpfix16. config zextra 10对于毎个zone多分配10个单位的内存这样可以保存一些白定义的和单元相关的数据 并通过fish进行访问如zextra_：可以按照zone的指针索引你H己定义的数据grextra :可以按照grid point的指针索引你自己定义的数据17. 应力分布plot cont smax确定应力分布方向principal stress tensor主应力的方向还是得根据弹郵!形力学计算得出18. set mech step stepjim force 50solve step 1000 force 50the 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", stepjim: This may be the variable defined by the "Fish".If the max step has reached the value "step_lim\the programme will stop without con sideratio n of the status of the system.solve step 1000 force: max. step = 1000, min. unbalanee force=50.19.ini(或 apply) zvel -0.5e-5 range x -.1 3.1 z 9.9 10.1用apply和ini的区别：边界条件和初始条件，边界条件不随时步变化，初始条件会变apply的话，那些位置的zvel就不会变了。

2().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,限定命令 发生作用的范I韦Iplane 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米的范围之内同时又在()

dip表示对称参照面与xy平面的夹角,对称参照面与xy平面的夹角在xy平面的投彫 是一条射线，dd表示y轴正向顺时针到那条射线的夹角,dip认为是平血从xoy平面向z 的负方向旋转角度，dd是平血的法向在xoy上的投影与y轴正向的夹角，注意是沿顺时针 旋转25. table相当于数已经赋值的数组26. 设置边长 gen zone brick(wedge...) size ... pO pl edge 10027. 显尔 shell 单元 plot add sei geom cid=on node=off.....厚度 t J IJ thickness 表现28. gen surface生成表面，显示模型中用到29. 显示单元数pr info30. 开挖和开挖后的冋填 model null range .....、model mohr range ..…31. n是通过节点面的外法向s是平行于这两个面的交线方向d是这个面上梯度最大的方向32. 显示模型中的节点号和单元号：plot modi i gpnum onplot modi i node on id onplot sei geomplot modi i id on —zone33. 输出实体单元一点的应力pr zone stress range 单元号也可以，限定x y z也可以34. 得到纟fl合力 plot sei recover就是弯矩，薄膜力得等売类求解顺序：先得到薄膜力，再由此得到売类应力35.ini state 0初始化舉性变形36.group right range group left not右边的组中不包含左边的组37.ini z mul」中mul是个乘子，这个命令可以用来移动网格，此句表示就是说移动10%看 原来的z坐标情况决定向那个方向移动。

38. pile表示结构类的桩单元39. plot surf yellow sei geom velo white将几何体的速度矢量显示为H色40. print sei pos得到节点的坐标41. fix x range x -0.1 0.1是将坐标在-0.1>x>0.1范|节I内的节点在x方向进行约束42. 设置内存(2D) "C:\Program Files\Itasca\flac500\exe\flacv_dp.exe" giic 30()43. 双击鼠标左键，可获得单元的id和positiono44. 如果不晓得一个命令后面跟什么关键字，可以在命令后加？可以获得一点help!比如 键入plot ?然后冋车，看会有什么提示 45. W1+Z可以选择El己想变大的区域46. F3,提取前命令47. ctrl+G：有彩色变黑白，有色曲线变成佶区别的点画线ctrl+R：变I川原来的状态ctrl+c： brings up the camera dialogctrl+1: cuttingplane4&育接按delete,图形就缩小,M是图形放大49. 不平衡力是由于数值计算处理中产生的系统内、外力Z差，通常用不平衡应力比来表征 迭代计算的精度，类似与ansys等求解,FLAC求解时，通过计算每一步的不平衡力，再与前面 的状态进行叠加，肓到求得的不平衡力小于规定的值.假如经过8000步的替代,达到我们容许 的精度，我们认为800()步Z后就达到了平衡状态.那么，这Z前没有达到平衡状态(如在100步 时)的应力或其他量值是否可用?手册中一个例子是在一个给定的速度下叠代到容许的精度， 得到了平均剪应力和平均剪应变的1111线.那么我理解在达到平衡状态Z前,「111线的应力应变 值都是不可信的”.显然这种理解不正确.敬请备位大侠答疑解惑！1. 在FLAC3D中，毎一个结点都是由8个单元I韦I绕，它们都对结点贡献结点力。

在平衡状 态时，这些力的代数和几乎等于零(即作用在这个结点的一个面上的力几乎与另一个面上的 力平衡)如果这个不平衡力达到了一个非零常数，则表明该模型正出现破坏和槊性流动 在计算时间步中，祭个网格结点的最大不平衡力被实时计算，并做为记录加以保存和图形显 示2. 常用最大不平衡力和节点内力的比值R来控制运算的阶段.3. wwiem说在(100步)时候，这些数字我认为这些过程值,在分析变化的过程中存在的，且可做 分析用的，你说的应力应变曲线是最终的结果，还是过稈曲线了.他是逐渐变化的过程，就象岩石压缩曲线样了，那你说这个曲线是没有可信度么贝能说是在R 小与某一个值的时候，表示节点力平衡或者处于没有加速度的稳定槊性流.这要结合速度的 变化了.如果速度his趋近于零,说明处于绝对平衡;稳定与某一个值，表明达到一个稳定状态.每一个 结点都是由8个单元围绕，它们都对结点贡献结点力不平衡力就是这些点做的作用力的合 力了.sub zone对gridpiont的内力。