flac学习笔记Word文档格式.docx

1、选择面里面的一个网格，右键，edit block pro或者按空格键-zone multiplier-左边三角形-group命名（绿色表示分组成功）。分组可以是平面分组，也可以纵向分组。show group 显示分组。none组包括没有被分组的网格。像上一次才能同步（delete）。分组撤销-object properties-group减号可以对进深分组，重点将某一区域加密做重点研究。变化率：object properties-ratio或者与分组一样，使网格有变化，距离明显。坐标输入点，移动点，删除点，建立其他线类型：object properties-line type案例：画50*10

2、0直角三角形，在三边画点并将上面两点连成直线。将两个区域做成面。划分各线网格数量及长度。进入进深按钮设置。回到网格页面，分组A、B。1.2 内置命令建模.表示与下面衔接（一个命令不完整才可以用）。敲前敲后要有空格英文字符和中文字符，注意加空格5.0没有 plot show ，而是plot zone设置透明：attribute-transparent坐标：list-Axis，坐标可以拖动并且放大缩小ctrl+R 恢复原来开始位置x左右 Y前后 Z上下 z=-0.01 0.01 （二者有空格）把0包括其中就可以了szz s应力 zz表示z轴方向可上可下plot hist 1；list可以查看很多选

3、项出图右键1.3 其他要点彭文斌第二章（完成一个简单分析计算样例）generatege zonezo brickbadd增加条目在键盘上分别按x y z小写字母键，可旋转当前视图中的X Y Z轴，相应的大写字母键则相反字母m键或M键则放大或缩小当前视图modelmo mohrmohpropertyPro bulkbu shearshfrictionfric cohesionctensionten（该值大，可以防止在初试加载就达到塑性极限）initialin（初始化网格相关值） densityde网格质量密度historyh采样或记录迭代时模型中变量的值，一次只能设置一个变量nstep迭代次数，

4、默认10unbalancenub最大不平衡力gridpointgp 采样网格节点的相关点zdisplacementzdis为z轴方向的位移mechanicalmec设置静态力学分析的有关参数forcefo 最大不平衡力solvesocontourcon 在当前视图中显示等值线图 disp位移量等值线图bcontourbcon绘制指定区域的等值线图planeP设置一个剖平面，剖面参数由后面的关键字确定boundarybo在视图中增加面的边界线框behindbe当前视图剖平面后面二 flac基础知识2.1操作和术语操作代码鼠标右键复制代码右键拖动模型旋转，点空白处选择其他命令轮滑放大缩小左键选择术

5、语Finite different grid 有限差分网格，也叫Mesh，是横跨一个或多个单元体分析物理的区域，意义存储很多参数Model boundary 模型边界。网格内的孔也是，叫做内部边界Constitution model 本构模型，可分配到单元体Null zone 空单元，即无用Sub-grid 子单元，可以组成网格。如水坝子网格放在地基子网格上Attached faces 绑定面，与子网格连接和绑定。绑定面一定要共面和接触，但每个面上的节点并不要求匹配，不同单元体密度的子网格可以绑定。Interface 分界面，是网格面，由子网格组成。在计算过程中可以分开（滑动、张开）2.2 约

6、定和默认命令句法必选参数、可选参数、任意个参数命令、关键字、数字可以用空格隔开或使用符号：（）、=分号（；）在前在后都表示不执行&为行连续。每一行不能超过80字符，不要用中文字符。每行连续总的加起来不能超过1024个字符Null可以简单施加到开挖的组名上Model null range group Tunnel对象命名采用大小写混用的形式，防止和命令及关键字冲突。方向约定坐标方向，应力和弹性力学方向一致。拉正压负，应变、孔隙压力也是单位约定不建议以mm为单位。非固定单位制常见：kg/N/m/s负号表示和正轴相反默认最大不平衡力ratio=1e-5.激活窗口才能用callsolve前可以对计算进

7、行规定：step 50000查看history:attribute-chart云图（contour过渡性的）和方块（zone以方块为单位，没有月牙）三 实体建模技术3.1 初级建模技术n1长度，r1网格数size。p1,p2顺序按照右手螺旋size 20 2 6 表示x向20个网格，以此类推。d1内部p1=（22 -10 0）也可以写成p1 22 -10 0ge zo生成基本形状网格常用关键词：add：用于以p0为原点的局部坐标系建模dimension:定义内部区域的尺寸edge：定义网格边长fill：定义网格内部填充区域默认基本网格建立：只是输入一个点，指定一个基本网格，默认长度10Ge z

8、o pyramid p0 0 12 局部坐标系用addge zo b p0 12 12 0 .p1 add 8 0 0 .p2 add 0 8 0 .p3 add 0 0 10 .size 8 8 10最后按照绝对坐标系：p0 12 12 0；p1 20 12 0 /p2 12 20 0 /p3 12 12 10利用edge生成网格：ge zo b p0 0 12 0 .p1 0 12 0 .edge 6 生成复杂的基本网格newge zo radb p0 12 -10 0 .p1 22 -10 0 .p2 12 0 0 .p3 12 -10 10 .d 3 6 4 .内部的长方体长宽高siz

9、e 3 6 4 7 外环绕单元个数，按照表格取。3 6 4是指内部的长方体，和外面没有缺口的面的网格数。其他的有缺口的大的面并保证不和内部长方体的网格数一致则按照7取网格。3.0版本用的是plot show，而不是plot zone3.2 中级建模技术镜像命令：ge zo reflect noramal（镜像平面单位法向量，如果其定义，则dd和dip不需要） （-1 0 0） origin（原点）=（12 -10 0）（p0点）：一个坐标一个点ge zo radc p0 0 0 0 .p1 50 0 0 .p2 0 100 0 .p3 0 0 50 .d 2 2 2 .ge zo radt p

10、0 0 0 0 .p1 0 0 -50 .p3 50 0 0 .d 1.8 2 2 .其中，1.8为前面一端的高度，后面的2是指后面一段的高度。中间的2是整个内部长方体进深多少。ge zo ref nor （-1 0 0） o = （0 0 0） （-1 0 0）是镜像平面的法向量，o = （0 0 0）是镜像平面过的点。原点可以是赋值的，不一定是全零ge zo reflect dd 0 dip 0dip 表示镜像平面的倾角，以x-y平面为基准，z轴方向为负dd指镜像平面的倾向，以正y轴方向为基准（0度），顺时针方向为正ratio 是指变化率，个数及方向按照网格库中原始形状及相关参数确定。 拉

11、伸用进深。orientation 改变方向edge 定义网格的边长，如果p1 p2 p3没给出，则由edge与p0确定。attach连接具有不同单元体大小的网格，相邻单元体大小相差需要一定范围，最好是成整数merge 合并（融合），小于设定的容差separate 分离接触面切片的方式：list-cutting planesActant：八分圆；成45视角的位置Plane：面Wedge：楔子；楔形物attribute 改变方向及其他参数出图的时候可能会有透明面，应该选择关闭，如需加网格之类的选择list的选项即可3.3 如何分组分组方法：内置命令，Extrusion，外接建模分组ge zo b

12、p0 0 0 0 p1 20 0 0 p2 0 100 p3 0 0 30 size 10 5 20 group AAgroup zone BB range z 5 30 x 5 15 y 0 10group TT range group BB not（要先命名一个BB分组，即取其否定反向）del zone range group BB删除分组内的内容斜面分组group zone GG range plane above（顶面） normal -2 0 -1.origin 15 0 0 plane normal -2 0 -1.origin 5 0 0 below（底面）default未分组p

13、rint information 查看节点数，网格数等分组命名时大小写都是一样的import 导入网格group zone 10 range y -100 -360 group 7 给分组7的某个区域在分组并命名10。四 本构模型与材料参数4.1 本构模型MODEL命令格式：MODEL overlay n range 或MODEL 应变硬化/软化塑性模型 strainsoftening4.2 材料参数定义材料参数命令PROPERTY 关键字 值 range . . .附加关键字Gauss_dev sGradient gx gy gz Uniform_dev s材料变形参数体积模量切变模量当泊松

14、比接近0.5时，K值显得很大，此时应该选取K的经验值。对于横向同性材料，需要5个弹性参数；对于各向异性需要9个弹性参数。后破坏参数的选取剪切膨胀：库伦摩尔；多节理；应变软化多节理剪切硬化/软化：应变软化；多节理模型体积硬化/软化：修正剑桥抗拉软化：多节理赋值方式：全局，范围，分组，fish赋值全局：model mech mohr attributes-histories可以选择在坐标中显示的量。范围：pro density 2000 range z 0 2输入命令后有刷新系统。pro young 2.7e10 ran group 177 ini 输入之后，后面可以变化。property反之fi

15、sh赋值：地下水default未分组五 边界条件格式APPLY：施加力学，流体，热边界APPLY keyword value 或APPLY remove FIX：选定节点速率不变，或使孔隙压力和温度不发生变化FIX FREE：对节点释放约束FREE 注意：1、 对节点、面和单元体施加的条件是相互独立的2、 除了在新值之前有add关键字或新值之后有history关键字之外，对于相同的关键字，新值替换旧值，如对一个节点，有一个恒定力（没用history）和一个步进变化力（用history）是可以叠加的。3、 一个条件只能用remove才能移去4、 施加速度总与节点相关5、 对于节点，指定速度时必须

16、基于相同的坐标；要么全局坐标，要么局部坐标6、 局部坐标不能施加到任一方向已固定的任何节点，全局坐标速度也不能施加在某一方向已固定的节点的相同方向上7、 对于面，指定应力时必须基于相同坐标应力边界szz/syy/szz/sxy关键字指定应力张量的某个分量施加面边界条件：nstress施加法向应力到一个面；dstress、sstress施加切应力到面，dstress为局部面倾向，sstress为局部面走向。全局和局部坐标应力不能施加在相同的面上。使用xforce、yforce、zforce关键字可以对网格施加单个力，此时不计边界面的面积，指定力给给定节点。apply nstress -1e6 r

17、ange plane dip 60 dd 270 origin 0.1,0,0 below施加渐变应力apply sxx -1e6 gradient 0,1e5,0 range z -100,0 ；在y方向渐变改变边界应力再一次用apply命令，有时需要使用apply remove 移除边界条件警告和建议对无位移约束的边界施加应力，物体可能像真世界一样反应。承受了边界应力的材料被开挖后也出现类似现象，但更隐蔽，物体在重力作用下初始化是平衡的，但移除材料后降低了重力，则物体向上运动。此时应该用固定底部边界来解决。位移边界全局坐标系xvel、yvel、zvel局部坐标系dvel、svel、nvel

18、局部坐标由节点的法向量定义，节点的法线方向是交汇于节点的所有面的平均法向矢量，对某一点所施加的速率坐标系统相同。在局部坐标，使用plane任意指定节点法线方向。resolutionsets the resolution in comparing two bracketing （stable and unstable） factors. （The default is .005 the mean value.）在apply之后要加上step才能显示处力，位移、边界。Manual P153gen zone brick size 4,4,4 p3 2,0,3.464model mech elaspr

19、op bulk 1e8 shear .3e8range name left_boun plane dip 60 dd 270 origin 0.1,0,0 aboverange name right_boun plane dip 60 dd 270 origin 3.9,0,0 belowsave ini; using the local axes planeapply nvel 0.1 plane dip 60 dd 270 range nrange left_bounapply nvel 0.1 plane dip 120 dd 90 range nrange right_bounstep

20、 1save local not using the local axes planeres iniapply nvel 0.1 range nrange left_bounapply nvel 0.1 range nrange right_boungen zone brick size 4 4 4model mech elasticapply nstress -1e6 range y 3.9 4.1 - the following conditions are not possible:apply nvel 0.0 plane dip 90 dd 180 range y -.1 .1appl

21、y nvel 0.0 plane dip 90 dd 270 range x -.1 .1 - the following conditions are possible:apply svel 0.0 plane dip 0 dd 0 range y -.1 .1apply dvel 0.0 plane dip 0 dd 0 range x -.1 .1运行第二个APPLY nvelocity时出现错误信息，因为x = 0, y = 0模型边界上的节点只能有一个定义的局部坐标。为避免这个错误，局部坐标系统只能有相同的dip和dd。相同的节点不允许用不同的局部坐标。任何节点只要用FIX命令固定，

22、则任何方向的局部坐标速率也不能施加，同样，全局坐标在相同方向不能施加（如，APPLY xvel就不能在FIX x之后使用）。既可以使用FIX，也可以使用INI来对节点的x、y、z方向施加速率。在计算过程中，无论什么力施加到 FIXed gridpoints上，the velocity at the start of stepping is retained。如果在之后的计算过程中需要零位移边界，可以在FIXed节点之前将合适的节点位移设为零。为了模拟试样之上压板的刚性运动，有如下的程序片段：fix x z range y -0.1 0.1apply yvel = 1e-5 range y -0

23、.1 0.1orfix x y z range y -0.1 0.1ini yvel = 1e-5 range y -0.1 0.1apply svel = 0 plane dip = 0 range y -0.1 0.1apply dvel = 0 plane dip = 0 range y -0.1 0.1apply nvel = 1e-5 plane dip = 0 range y -0.1 0.1APPLY. . . interior可以施加速率到内部节点，FIX 和 INITIAL可以在任何节点、边界、内部operate。节点移动时，可以简单地用FISH获得xforce、yforce

24、、zforce反作用力的大小。如果要求速率不均匀，可以用gradient关键字，但对于复杂的速率，或过程中进行改变，则需要编写fish函数。通常，与受载体相比十分坚硬的（硬20），处理呈刚性边界；反之（1/20）为应力边界.人工边界对称平面FIX x或INITIAL xvelocity = 0.0。APPLY nvelocity = 0剖断面六 初始条件格式：INITIAL 关键字 grad gx gy gz range均匀应力无重力，manual P159ini sxx=-5e6 syy=-1e7 szz=-5e6将、设置成压缩应力，贯穿整个网格。但其不能保证应力平衡。第一，给定应力可能违反

25、分配到网格非线性模型的屈服标准，step后可能塑性流动，使其重新分布。此时可以用PLOT zone colorby state）检查；第二，边界指定应力不等于给定的初始应力，step后就开始移动，可以用PLOT Velocity检查。gen zone brick size 6 6 6ini sxx=-5e6 syy=-1e7 szz=-2e7apply sxx=-5e6 range x -0.1 0.1apply sxx=-5e6 range x 5.9 6.1apply syy=-1e7 range y -0.1 0.1apply syy=-1e7 range y 5.9 6.1apply szz=-2e7 range z -0.1 0.1apply szz=-2e7 range z 5.9 6.1如果边界用FIX固定，而不是应力控制，则初始应力将自动平衡，不需要APPLY。渐变应力均值材料七 求解八 后处理