FLAC2D模拟代码解释.docx
《FLAC2D模拟代码解释.docx》由会员分享，可在线阅读，更多相关《FLAC2D模拟代码解释.docx（10页珍藏版）》请在冰豆网上搜索。
FLAC2D模拟代码解释
1实例
new
title新建命名
grid6,5;定义网格数目
gen-15,-16-15,2060,2060,-16;四角点坐标
gen00010100101000i=1,201j=1,11;把大的范围浓缩到小范围内
propbulk=1e9s=5e8d=2500f=30c=5e5;定义物质属性
propbulk=1e8shear=.3e8fric=35;体积模量、剪切模量、内摩擦角
propdens=1000coh=1e10ten=1e10;密度、粘聚力、抗拉强度
interface1asidefrom1,10to20,10bsidefrom1,11to20,11;定义节理，将节理放大，随后跟节理的属性
mnj=11;interfaces,把j=11断开
genline50.021.553.521.5;画线
mssj3335
propbulk=10e9s=4.6e9d=1400f=30c=5e5i=12,20j=33,35;赋予这个范围属性
propftab1ctab2j3335
tab10270.01251.023
tab205e50.013e51.02e5;从mss开始是一组固定语言，应变软化模型先不学
mnregioni=1020j=1020;开挖
mnregion10234;开挖，跟前面的genline组合使用，开挖一个块就会挖整个
mnregion10230;不定义genline会开挖一层
inisyy=-13e8;初始设置，y方向为压应力
fixy=1;固定y=1这条边界，不变形
setg=9.81;定义重力加速度
setaspecta;调整屏幕的高宽比，使输出不变形
setlarge;大变形，坐标更新
setsmall;小变形，坐标不更新
setlogonoff;在当前盘上打开flac.log的文件，记录FLAC的运行进程
settime
setstep;设置时间与步数的限制
appsyy=-12.23e6j=129;使y=129边界（即上边界）受向下的力
histnstep=10;每十步记录一次
HISUNBAL;设置不平衡力
HISsyyi=100j=68;记录此点的三向应力变化，默认编号为1
HISsxxi=100j=68;
HISszzi=100j=68;
str;支护，后跟属性或支护材料
step;计算步数
solve;计算
solvestep;时步数目，缺省500时步
solvetime;运行分钟数，缺省5分钟
solveforce;不平衡力，缺省100
pause;暂停
continue;继续
save;保存
mm;modelmohr摩尔模型
seg;分节点
node1tonode;托盘的节点，分两段就有三个节点
plboustrubeamnodered;即可查看托盘节点号
plotgridnumber;标注坐标在格子中间
plotgridgnum;标注在点上，查看时用的，简化为plgnum或gnu
plotgridnum;可以查看节点号。
plotfix;查看边界条件
plotgridmark;查看预先留设的开挖边界
plgrifred;查看interface及编号
plsyy（sxxsdisydissigl）fill;查看模型的云图，带fill为查看云图，不带fill为查看等值线图
plgrpl;查看塑性区图
plgridvel（disp）red;查看模型的矢量图
plgridifaceidred;查看INTERFACE，id为其编号
plogrognumon;查看点的坐标
plotkeywordmagfacvalue;显示计算后模型的变形状态，value取大于0的
setplotdxf256;回车
copy1.dxf;保存图片
title;定义文件名为“斜坡稳定性分析”
SlopeStabilityAnalysis
plotgrid即plg;查看
SClinnx1,y1x2,y2;屏幕上做扫描线，编号为1到5，起始点坐标
plothistn;查看位移量，hist的编号
2FLAC2D学习心得
new
title
ABCDE（命名文件）
gridab（网格数a×b）
1，对于大规模建模时，如何利用GEN命令生成局域性的板块
如果其中用到到SAME时，是如何完成建模的?
1）当用GRIDIJ生成网络后,程序默认原点是网格左下角.
2）GEN命令,对网格坐标进行重分配时,如果新的分区坐标与GRID生成网格中坐标重合时,用SAME代替
2，如何建立起模型的破坏？
赋值时如何体现出来
3，如何设置不同层次岩土的参数?
4,GEN命令将使网格沿x,y向以1.25传于前一个网跨的比率扩展至1000个单位。
注意，如果比例介于0～1之间，网格的尺寸将随坐标值的增加而减少。
例如，发布以下命令：
new
grid10，10
me
gen–100，0–100，1000，1000,0rat0.80，1.25
plotgrid
你将会看见沿x,y向逐渐变小的网格。
5,移动网点
可以用Initial命令，把一个点或若干点从当前位置移到一个新位置。
以下各命令建立了一个网格，并用Initial命令来使它变形：
grid5,5
modelelastic
gen0,00,1010,1010,0（沿x,y向规定网点坐标为0～10）
Inix＝－12I＝1j＝6（把左上角水平移动－2个单元）
Inix＝12i＝6（把各右边界网点各移2个单元，注意，因为没有给j的范围，取为整个范围）
注意，可以用Initial命令，把任何网点移至任何位置。
当然，各单元不能重迭。
如果重迭了将给出一个为“BAD.GEOMETRY”的警告信息，直到校正了网格结构中的错误为止，flac不会继速执行的。
为了得到合理的解，应把各域的高宽比的实际范围，保持在约1:
10或更小。
6,应力荷载
grid10,10
me
props=5.7e9b=11.1e9d=2000
fixxyi=1
fixxyj=1
fixxyi=11
fixxyj=11
intsxx=－30e6syy=－20e6（场地应力）
hisnstep=1
hisxdisi=5j=5（回采工作与壁处的x位移）
solve
;newexcavate
Mnulli=3.4j=3.6（开挖非1）
plotgrid
step200（步入平衡）
plohis1（校核平衡）
plotgridstr（绘应力，网格）
modmohri=3.4j=3.6（用sangfill特征）
modnulli=3.4j=7（开挖非2）
inixd=0yd=0（初始化位移，以致于看出开挖引起的变化）
step200（再次步入接近平衡）
plohis1
plotgridstr（绘应力及网格）
plotcoh（绘不同的材料类型）
3问答
.FLAC3D是有限元程序吗？
答：
不是！
是有限差分法.
仿2.最先需要掌握的命令有哪些？
www
.答：
需要掌握gen,ini,app,plo,solve等建模、初始条件、边界条件、后处理和求解的命令。
3.怎样看模型的样子?
~5i
答：
ploblogro可以看到不同的group的颜色分布
4.怎样看模型的边界情况？
答：
plogpfixredsk
5.怎样看模型的体力分布?
c
答：
plofapredsk
6.怎样看模型的云图？
答：
位移：
plocondis（xdis,ydis,zdis）
应力：
ploconsz（sy,sx,sxy,syz,sxz）
7.怎样看模型的矢量图？
答：
plodis（xdis,ydis,zdis）
8.怎样看模型有多少单元、节点？
答：
ploinfo仿
9.怎样输出模型的后处理图？
答：
File/Printtype/Jpgfile，然后选择File/Print，将保存格式选择为jpe文件
10.怎样调用一个文件？
答：
File/call或者call命令
11.如何施加面力？
答：
appnstress,
12.如何调整视图的大小、角度？
答：
综合使用x,y,z,m,Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键
13.如何进行边界约束？
答：
fixxran（约束的是速度，在初始情况下约束等效于位移约束）
14.如何知道每个单元的ID？
仿真
答：
用鼠标双击单元的表面，可以知道单元的ID和坐标仿真
15.如何进行切片？
答：
plosetplaneori（点坐标）norm（法向矢量）
ploconszplane（显示z方向应力的切片）
16.如何保存计算结果？
答：
save+文件名
17.如何调用已保存的结果？
答：
rest+文件名；或者File/Restor仿k,Z4\5rD+D
18.如何暂停计算？
答：
Esc
19.如何在程序中进行暂停，并可恢复计算？
答：
在命令中加入pause命令，用continue进行继续仿真仿
20.如何跳过某个计算步?
答：
在计算中按空格键跳过本次计算，自动进入下一步
21.Fish是什么东西？
答：
是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能
22.Fish是否一定要学？
答：
可以不用，需要的时候查Mannual获得需要的变量就可以了仿
23.FLAC3D允许的命令文件格式有哪些？
答：
无所谓，只要是文本文件，什么后缀都可以
24.如何调用一些可选模块？
答：
configdyn（fluid,creep,cppudm）
25.如何在圆柱体四周如何施加约束条件？
答：
可以用fix...rancylinderend1end2radiusr1cylinderend1end2radiusr2not，其中r226.如何能把一个PLOT的图像数据导出来以便用其他软件绘图？
答：
用setlogon命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了
27.关于dd、dip的理解：
dip-directionangle,measuredintheglobalxy-planeclockwisefromthepositivey-axis这个是从y轴到x顺时针方向所转过的角度；dipdipangle,measuredinthenegativez-directionfrom theglobalxy-plane这个是xy平面与z轴负方向的夹角。
dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线，dd表示y轴正向顺时针到那条射线的夹角。
首先应该按照dd的方向大体确定这个面的朝向，dd指的是从y轴正方向按顺时针（clockwise）方向转向所要确定面的法线方向在xy平面上的投影的夹角，然后再确定dip，dip指的是从xy平面转向所要确定的平面的角.
28. 如何调用FISH文件？
新手调用FISH文件出错的原因：
是：
文件名为fishcall.fis.txt（因为在记事本里这个txt的后缀是默认的）。
修改方法是：
将fishcall.fis.txt里面的内容全部复制到新的记事本里，再按保存，保存时候将保存类型改为所有文件（*,*）,这时你便可在文件名栏内输入fishcall.fis，类型也就正确了.
29. 如何理解Tension-pshear-ptension-nshear-n表示的含义？
:
在计算的循环里面,每个循环中,每个zone都依据failurecriterion处于不同的状态,shear和tension分别表示因受剪和受拉而处于塑性状态.n表示now,p表示previous,即分别表示在现在和以前的循环当中处于塑性状态.
30.绘制相关曲线命令的一些解释
1）在plothistmvsn的形式里,m代表y轴,n代表x轴（不管m,n的正负）;
2）"-"表示对其值作"mirror",如果大家使用AUTOCAD的话就明白"mirror"的含义了,比如上面的hist1gpydisp的值是0~100,那么vs-1就变成-100~0.以此类推.
3）为什么要这么做呢?
主要是为了符合岩土工程的习惯需要,例如桩载荷试验曲线,沉降曲线等都不是画在第一象限内.
31.如何查看剖面上任一点的位移值？
printgpdisprangex（）y（）z（）或
printgpdisprangeid
plotgpdisprangex（）y（）z（）或
plotgpdisprangeid
32.fixxyzrangex-0.10.1..........................aS
fixxrangex-0.10.1...............................b
请问a式与b式有何区别？
a表示在x=0这个平面是固定的
b表示在x=0这个平面只x方向不能发生运动
33.initial与apply有何区别？
initial初始化命令，如初始化计算体的应力状态等；
apply边界条件限制命令，如施加边界的力、位移等约束等。
initial的应力状态会随计算过程的发生而发生改变，一般体力需要初始化，而apply施加的边界条件不会发生变化。
34.如何查看各个时段不平衡力的具体数值？
采用his来记录计算，包括位移应力等命令hisunbal
hisgp（zone）zdisrange（000）或者id=？
导出数据命令hiswritenvsmbegin时步 end 时步filefilename.his
n表示纪录的idm表示时步
要导出不平衡力的具体数值
hisunbal
step100000orsolve
hiswrite1vsstepbegin1end1000file123.his
使用上述命令就可以查看各个时步下的不平衡力的具体数值
35.如何显示变形轮廓线的命令？
ploskemagf10 其中10为放大系数
4补充：
FLAC2d命令的FAQ：
1、怎样开始入门？
仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM1giE:
R.LUp+K:
U
答：
初步学习首先看说明书User'sgurid的Section2Gettingstarted和Section3ProblemslovingwithFLAC
2、哪些命令需要掌握？
R-E
C8h@ch
答：
gridgenpropinterfaceinifixsethiststepsolvepausecontsav等，在学习的过程中遇到这些命令和不会的命令应该查看说明书的CommandReference，里面讲的很清楚，特别需要注意的是一定要仔细看一下CommandReference开头的OptionalRangePhrases，这是讲范围的，因为很多命令的施加都需要加到模型的一部分，这里就告诉你怎样来表示这一部分。
这一原则同样适合与其他5个ITASCA软件。
还有一点就是六个软件（FLACFALC3DUDEC3DECPFCPFC3D）的range大不相同，要注意区分。
3、FLAC和FLAC3D建模有什么区别？
SimWe仿真论坛tR7Jx.V|0|
答：
FLAC和FLAC3D基本的建模思想大致是一样的，不同的是FLAC是首先用GRID命令把整体的网格画出来，然后用MODELLNULL命令把网格分成许多部分，再用GEN命令把这些部分一块一块搭接起来。
然后再用GENLINE，GENARC等命令画出要开挖的边界，在FLAC中产生网格的命令只有一个GRID，而且只能用一次。
而FLAC3D是直接用GEN命令一块一块搭接，GEN可以用很多次。
4、怎样查看模型？
答：
plotgrid可以查看网格，plotgridnum可以查看节点号。
仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent]&ldUun
5、怎样查看边界条件？
仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent[3dJaK
答：
plotfix~2e6t6_qa:
s
6、怎样查看预先留设的开挖边界？
答：
plotgridmark。
7、怎么查看模型的云图？
答：
plsyy（sxxxdisydissig1..）fill，注意，带fill关键字就是查看云图，不带fill关键字就是查看等直线图。
8、怎样查看模型的矢量图？
答：
plgridvel（disp）red
9、怎样查看INTERFACE？
答：
plgridifaceidred
printifaceid;id是interface的编号
10、FLAC怎么调整视图？
tP3d,Y:
@N-Z
答：
ctrl+z或shift+z,然后点击要放大部位的左上角和右下角就可以放大某个部位，如果要缩小，双击屏幕。
升级会员 