flac学习笔记教学内容.docx
《flac学习笔记教学内容.docx》由会员分享，可在线阅读，更多相关《flac学习笔记教学内容.docx（17页珍藏版）》请在冰豆网上搜索。
flac学习笔记教学内容
一初次建模
1.1extrusion建模
新建命名字母开头
菜单栏：
import导入
叉号：
删除文件（撤销——编辑ctrl+z）
画点先划线，点在线上成点。
空中画点（右键输坐标如0,0）。
删除线——delete。
画水平竖直线：
按住ctrl。
右键设置显示菜单和ctrl时定位角度
objectproperties显示、修改信息。
选择线不改大小，改点坐标
控制三个按钮，控制参数、视角、information
建立面（只能是三线或者是四线，否则会报错），默认网格数3，橘黄色（线比其他地方粗或密），增加网格先选择某一方向的线再用objectproperties-zones增加数目.面内划线按ctrl
autozone按钮:
改变单元大小，即网格长度、进深，网格总数。
进深：
ZY方向，改变Y方向长度，zone表示线的长度。
extrude（串联关系）:
由上一部命令激活，但autozone改变大小时并不同步，此时再按该按钮就会同步（delete）,改变进深之后和模型（plot）一致。
改了就要同步。
分组：
选择面里面的一个网格，右键，editblockpro或者按空格键--zonemultiplier--左边三角形-group命名（绿色表示分组成功）。
分组可以是平面分组，也可以纵向分组。
showgroup显示分组。
none组包括没有被分组的网格。
像上一次才能同步（delete）。
分组撤销-objectproperties-group减号
可以对进深分组，重点将某一区域加密做重点研究。
变化率：
objectproperties-ratio或者与分组一样，使网格有变化，距离明显。
坐标输入点，移动点，删除点，建立其他线类型：
objectproperties-linetype
案例：
画50*100直角三角形，在三边画点并将上面两点连成直线。
将两个区域做成面。
划分各线网格数量及长度。
进入进深按钮设置。
回到网格页面，分组A、B。
1.2内置命令建模
...表示与下面衔接（一个命令不完整才可以用）。
敲前敲后要有空格
英文字符和中文字符，注意加空格
5.0没有plotshow，而是plotzone
设置透明：
attribute-transparent
坐标：
list-Axis，坐标可以拖动并且放大缩小
ctrl+R恢复原来开始位置
x——左右Y——前后Z——上下z=-0.010.01（二者有空格）把0包括其中就可以了
szzs应力zz表示z轴方向可上可下
plothist1；list可以查看很多选项
出图右键
1.3其他要点
彭文斌第二章（完成一个简单分析计算样例）
generate——gezone——zobrick——b
add——增加条目
在键盘上分别按xyz小写字母键，可旋转当前视图中的XYZ轴，相应的大写字母键则相反
字母m键或M键则放大或缩小当前视图
model——momohr——moh
property——Probulk——bushear——sh
friction——friccohesion——c
tension——ten（该值大，可以防止在初试加载就达到塑性极限）
initial——in（初始化网格相关值）density——de网格质量密度
history——h采样或记录迭代时模型中变量的值，一次只能设置一个变量
nstep迭代次数，默认10
unbalance——nub最大不平衡力
gridpoint——gp采样网格节点的相关点
zdisplacement——zdis为z轴方向的位移
mechanical——mec设置静态力学分析的有关参数
force——fo最大不平衡力
solve——so
contour——con在当前视图中显示等值线图disp位移量等值线图
bcontour——bcon绘制指定区域的等值线图
plane——P设置一个剖平面，剖面参数由后面的关键字确定
boundary——bo在视图中增加面的边界线框
behind——be当前视图剖平面后面
二flac基础知识
2.1操作和术语
操作
代码鼠标右键复制代码
右键拖动模型旋转，点空白处选择其他命令
轮滑放大缩小
左键选择
术语
Finitedifferentgrid有限差分网格，也叫Mesh，是横跨一个或多个单元体分析物理的区域，意义存储很多参数
Modelboundary模型边界。
网格内的孔也是，叫做内部边界
Constitutionmodel本构模型，可分配到单元体
Nullzone空单元，即无用
Sub-grid子单元，可以组成网格。
如水坝子网格放在地基子网格上
Attachedfaces绑定面，与子网格连接和绑定。
绑定面一定要共面和接触，但每个面上的节点并不要求匹配，不同单元体密度的子网格可以绑定。
Interface分界面，是网格面，由子网格组成。
在计算过程中可以分开（滑动、张开）
2.2约定和默认
命令句法
<>必选参数、[]可选参数、…任意个参数
命令、关键字、数字可以用空格隔开或使用符号：
（）、=
分号（；）在前在后都表示不执行
&为行连续。
每一行不能超过80字符，不要用中文字符。
每行连续总的加起来不能超过1024个字符
Null可以简单施加到开挖的组名上
ModelnullrangegroupTunnel
对象命名采用大小写混用的形式，防止和命令及关键字冲突。
方向约定
坐标方向，应力和弹性力学方向一致。
拉正压负，应变、孔隙压力也是
单位约定
不建议以mm为单位。
非固定单位制
常见：
kg/N/m/s
负号表示和正轴相反
默认最大不平衡力ratio=1e-5.激活窗口才能用call
solve前可以对计算进行规定：
step50000
查看history:
attribute-chart
云图（contour过渡性的）和方块（zone以方块为单位，没有月牙）
三实体建模技术
3.1初级建模技术
n1长度，r1网格数size。
p1,p2顺序按照右手螺旋
size2026表示x向20个网格，以此类推。
d1内部
p1=（22-100）也可以写成p122-100
gezo生成基本形状网格常用关键词：
add：
用于以p0为原点的局部坐标系建模
dimension:
定义内部区域的尺寸
edge：
定义网格边长
fill：
定义网格内部填充区域
默认基本网格建立：
只是输入一个点，指定一个基本网格，默认长度10
Gezopyramidp0012
局部坐标系用add
gezobp012120...
p1add800...
p2add080...
p3add0010...
size8810
最后按照绝对坐标系：
p012120；p120120/p212200/p3121210\
利用edge生成网格：
gezobp00120...
p10120...
edge6…
size8810
生成复杂的基本网格
new
gezoradbp012-100...
p122-100...
p21200...
p312-1010...
d364...内部的长方体长宽高
size3647外环绕单元个数，按照表格取。
364是指内部的长方体，和外面没有缺口的面的网格数。
其他的有缺口的大的面并保证不和内部长方体的网格数一致则按照7取网格。
3.0版本用的是plotshow，而不是plotzone
3.2中级建模技术
镜像命令：
gezoreflectnoramal（镜像平面单位法向量，如果其定义，则dd和dip不需要）（-100）origin（原点）=（12-100）（p0点）：
一个坐标一个点
new
gezoradcp0000...
p15000...
p201000...
p30050...
d222...
gezoradtp0000...
p100-50...
p201000...
p35000...
d1.822...其中，1.8为前面一端的高度，后面的2是指后面一段的高度。
中间的2是整个内部长方体进深多少。
gezorefnor（-100）o=（000）（-100）是镜像平面的法向量，o=（000）是镜像平面过的点。
原点可以是赋值的，不一定是全零
gezoreflectdd0dip0
dip表示镜像平面的倾角，以x-y平面为基准，z轴方向为负
dd指镜像平面的倾向，以正y轴方向为基准（0度），顺时针方向为正
ratio是指变化率，个数及方向按照网格库中原始形状及相关参数确定。
拉伸用进深。
orientation改变方向
edge定义网格的边长，如果p1p2p3没给出，则由edge与p0确定。
attach连接具有不同单元体大小的网格，相邻单元体大小相差需要一定范围，最好是成整数
merge合并（融合），小于设定的容差
separate分离接触面
切片的方式：
list-cuttingplanes
Actant：
八分圆；成45°视角的位置
Plane：
面
Wedge：
楔子；楔形物
attribute改变方向及其他参数
出图的时候可能会有透明面，应该选择关闭，如需加网格之类的选择list的选项即可
3.3如何分组
分组方法：
内置命令，Extrusion，外接建模分组
gezobp0000p12000p20100p30030size10520groupAA
groupzoneBBrangez530x515y010
groupTTrangegroupBBnot（要先命名一个BB分组，即取其否定反向）
delzonerangegroupBB删除分组内的内容
斜面分组groupzoneGGrangeplaneabove（顶面）normal-20-1...
origin1500planenormal-20-1...
origin500below（底面）
default未分组
printinformation查看节点数，网格数等
分组命名时大小写都是一样的
import导入网格
groupzone10rangey-100-360group7给分组7的某个区域在分组并命名10。
四本构模型与材料参数
4.1本构模型
MODEL命令格式：
MODEL<关键字>[overlayn][range…]
或MODEL[overlayn][range…]
应变硬化/软化塑性模型strainsoftening
4.2材料参数
定义材料参数命令
PROPERTY关键字值<关键字...><...>
附加关键字
Gauss_devs
Gradientgxgygz
Uniform_devs
材料变形参数
体积模量
切变模量
当泊松比接近0.5时，K值显得很大，此时应该选取K的经验值。
对于横向同性材料，需要5个弹性参数；对于各向异性需要9个弹性参数。
后破坏参数的选取
剪切膨胀：
库伦摩尔；多节理；应变软化多节理
剪切硬化/软化：
应变软化；多节理模型
体积硬化/软化：
修正剑桥
抗拉软化：
应变软化；多节理
赋值方式：
全局，范围，分组，fish赋值
全局：
modelmechmohr
attributes-histories可以选择在坐标中显示的量。
范围：
prodensity2000rangez02
输入命令后有刷新系统。
分组：
proyoung2.7e10rangroup177
ini输入之后，后面可以变化。
property反之
fish赋值：
地下水default未分组
printinformation查看节点数，网格数等
分组命名时大小写都是一样的
import导入网格
groupzone10rangey-100-360group7给分组7的某个区域在分组并命名10。
五边界条件
格式
APPLY：
施加力学，流体，热边界
APPLYkeywordvalue
或APPLYremove
FIX：
选定节点速率不变，或使孔隙压力和温度不发生变化
FIX…
FREE：
对节点释放约束
FREE…
注意：
1、对节点、面和单元体施加的条件是相互独立的
2、除了在新值之前有add关键字或新值之后有history关键字之外，对于相同的关键字，新值替换旧值，如对一个节点，有一个恒定力（没用history）和一个步进变化力（用history）是可以叠加的。
3、一个条件只能用remove才能移去
4、施加速度总与节点相关
5、对于节点，指定速度时必须基于相同的坐标；要么全局坐标，要么局部坐标
6、局部坐标不能施加到任一方向已固定的任何节点，全局坐标速度也不能施加在某一方向已固定的节点的相同方向上
7、对于面，指定应力时必须基于相同坐标
应力边界
szz/syy/szz/sxy关键字指定应力张量的某个分量
施加面边界条件：
nstress——施加法向应力到一个面；dstress、sstress——施加切应力到面，dstress为局部面倾向，sstress为局部面走向。
全局和局部坐标应力不能施加在相同的面上。
使用xforce、yforce、zforce关键字可以对网格施加单个力，此时不计边界面的面积，指定力给给定节点。
applynstress-1e6rangeplanedip60dd270origin0.1,0,0below
施加渐变应力
applysxx-1e6gradient0,1e5,0rangez-100,0；在y方向渐变
改变边界应力
再一次用apply命令，有时需要使用applyremove移除边界条件
警告和建议
对无位移约束的边界施加应力，物体可能像真世界一样反应。
承受了边界应力的材料被开挖后也出现类似现象，但更隐蔽，物体在重力作用下初始化是平衡的，但移除材料后降低了重力，则物体向上运动。
此时应该用固定底部边界来解决。
位移边界
全局坐标系——xvel、yvel、zvel
局部坐标系——dvel、svel、nvel
局部坐标由节点的法向量定义，节点的法线方向是交汇于节点的所有面的平均法向矢量，对某一点所施加的速率坐标系统相同。
在局部坐标，使用plane任意指定节点法线方向。
resolution——setstheresolutionincomparingtwobracketing（stableandunstable）factors.（Thedefaultis.005×themeanvalue.）
在apply之后要加上step才能显示处力，位移、边界。
ManualP153
genzonebricksize4,4,4p32,0,3.464
modelmechelas
propbulk1e8shear.3e8
rangenameleft_bounplanedip60dd270origin0.1,0,0above
rangenameright_bounplanedip60dd270origin3.9,0,0below
saveini
;usingthelocalaxesplane
applynvel0.1planedip60dd270rangenrangeleft_boun
applynvel0.1planedip120dd90rangenrangeright_boun
step1
savelocal
;notusingthelocalaxesplane
resini
applynvel0.1rangenrangeleft_boun
applynvel0.1rangenrangeright_boun
step1
genzonebricksize444
modelmechelastic
propbulk1e8shear.3e8
applynstress-1e6rangey3.94.1
;--thefollowingconditionsarenotpossible:
;applynvel0.0planedip90dd180rangey-.1.1
;applynvel0.0planedip90dd270rangex-.1.1
;--thefollowingconditionsarepossible:
applysvel0.0planedip0dd0rangey-.1.1
applydvel0.0planedip0dd0rangex-.1.1
运行第二个APPLYnvelocity时出现错误信息，因为x=0,y=0模型边界上的节点只能有一个定义的局部坐标。
为避免这个错误，局部坐标系统只能有相同的dip和dd。
相同的节点不允许用不同的局部坐标。
任何节点只要用FIX命令固定，则任何方向的局部坐标速率也不能施加，同样，全局坐标在相同方向不能施加（如，APPLYxvel就不能在FIXx之后使用）。
既可以使用FIX，也可以使用INI来对节点的x、y、z方向施加速率。
在计算过程中，无论什么力施加到FIXedgridpoints上，thevelocityatthestartofsteppingisretained。
如果在之后的计算过程中需要零位移边界，可以在FIXed节点之前将合适的节点位移设为零。
为了模拟试样之上压板的刚性运动，有如下的程序片段：
fixxzrangey-0.10.1
applyyvel=1e-5rangey-0.10.1
or
fixxyzrangey-0.10.1
iniyvel=1e-5rangey-0.10.1
or
applysvel=0planedip=0rangey-0.10.1
applydvel=0planedip=0rangey-0.10.1
applynvel=1e-5planedip=0rangey-0.10.1
APPLY...interior可以施加速率到内部节点，FIX和INITIAL可以在任何节点、边界、内部operate。
节点移动时，可以简单地用FISH获得xforce、yforce、zforce反作用力的大小。
如果要求速率不均匀，可以用gradient关键字，但对于复杂的速率，或过程中进行改变，则需要编写fish函数。
通常，与受载体相比十分坚硬的（硬20），处理呈刚性边界；反之（1/20）为应力边界.
人工边界
对称平面——FIXx或INITIALxvelocity=0.0。
APPLYnvelocity=0
剖断面——
六初始条件
格式：
2003年，全年商品消费价格总水平比上年上升1%。
消费品市场销售平稳增长。
全年完成社会消费品零售总额2220.64亿元，比上年增长9.1%。
INITIAL<关键字>[关键字]<值>[gradgxgygz][range]
300-400元1632%均匀应力——无重力，manualP159
（二）大学生对DIY手工艺品消费态度分析inisxx=-5e6syy=-1e7szz=-5e6
将
、
、
设置成压缩应力，贯穿整个网格。
但其不能保证应力平衡。
第一，给定应力可能违反分配到网格非线性模型的屈服标准，step后可能塑性流动，使其重新分布。
此时可以用PLOTzonecolorbystate）检查；第二，边界指定应力不等于给定的初始应力，step后就开始移动，可以用PLOTVelocity检查。
在调查中我们注意到大多数同学都比较注重工艺品的价格，点面氛围及服务。
genzonebricksize666
modelmechelas
inisxx=-5e6syy=-1e7szz=-2e7
因为是连锁店，老板的“野心”是开到便利店那样随处可见。
所以办了积分卡，方便女孩子到任何一家“漂亮女生”购物，以求便宜再便宜。
applysxx=-5e6rangex-0.10.1
“碧芝”的成功归于他的唯一，这独一无二的物品就吸引了各种女性的眼光。
applysxx=-5e6rangex5.96.1
applysyy=-1e7rangey-0.10.1
applysyy=-1e7rangey5.96.1
applyszz=-2e7rangez-0.10.1
applyszz=-2e7rangez5.96.1
但这些困难并非能够否定我们创业项目的可行性。
盖茨是由一个普通退学学生变成了世界首富，李嘉诚是由一个穷人变成了华人富豪第一人，他们的成功表述一个简单的道理：
如果你有能力，你可以从身无分文变成超级富豪；如果你无能，你也可以从超级富豪变成穷光蛋。
如果边界用FIX固定，而不是应力控制，则初始应力将自动平衡，不需要APPLY。
渐变应力
为此，装潢美观，亮丽，富有个性化的店面环境，能引起消费者的注意，从而刺激顾客的消费欲望。
这些问题在今后经营中我们将慎重考虑的。
均值材料
除了“漂亮女生”形成的价格，优惠等条件的威胁外，还有“碧芝”的物品的新颖性，创意的独特性等，我们必须充分预见到。
七求解
八后处理
升级会员 