用Tecplot对Flac3D计算结果进行后处理
FLAC3D知识介绍
FLAC 3D基础知识介绍一、概述FLAC（Fast Lagrangian Analysis of Continua）由美国Itasca公司开发的。
目前，FLAC有二维和三维计算程序两个版本，二维计算程序V3.0以前的为DOS版本，V2.5版本仅仅能够使用计算机的基本内存64K），所以，程序求解的最大结点数仅限于2000个以内。
1995年，FLAC2D已升级为V3.3的版本，其程序能够使用护展内存。
因此，大大发护展了计算规模。
FLAC3D是一个三维有限差分程序，目前已发展到V3.0版本。
FLAC3D的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入各种命令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。
因此，采用FLAC程序进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排列，形成可以完成一定计算任务的命令文件。
FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。
调整三维网格中的多面体单元来拟合实际的结构。
单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发生变形和移动（大变形模式）。
FLAC3D采用的显式拉格朗日算法和混合-离散分区技术，能够非常准确的模拟材料的塑性破坏和流动。
由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。
三维快速拉格朗日法是一种基于三维显式有限差分法的数值分析方法，它可以模拟岩土或其他材料的三维力学行为。
三维快速拉格朗日分析将计算区域划分为若干四面体单元，每个单元在给定的边界条件下遵循指定的线性或非线性本构关系，如果单元应力使得材料屈服或产生塑性流动，则单元网格可以随着材料的变形而变形，这就是所谓的拉格朗日算法，这种算法非常适合于模拟大变形问题。
三维快速拉格朗日分析采用了显式有限差分格式来求解场的控制微分方程，并应用了混合单元离散模型，可以准确地模拟材料的屈服、塑性流动、软化直至大变形，尤其在材料的弹塑性分析、大变形分析以及模拟施工过程等领域有其独到的优点。
FLAC3d后处理步骤
1、首先rest123.sav，然后再callFLAC3d2Tecplop生成Tec10.dat2、打开Tecplot程序，执行【File】/【Load Data File】命令，打开数据文件对话框，选择刚刚生成的tec10.dat文件，并执行【打开】命令，在出现的Select Initial Plot 对话框中，单击OK按钮。
并在Tecplot左上角显示的3维笛卡尔坐标系下选择【3D Cartasion 】3、在Tecplot中选中【Data】→【Extract】→【Slice from Plane】，打开切面对话框。
在对话框中我们可以选中x=6/y=5/z=7等切片，并单击Extract.4、在Tecplot中选中【File】→【White Data File】，打开【White Data FileOptions】对话框，在对话框中选中ASCII、Point。
并在下面的Zone（s）选项中选中你所要的2 SLc=3，在Variable(s)中同时选中1 x（m）,2 y(m)及你所要的13 szz（pa）等，并单击OK按钮。
保存的文件自己最好独自命名（如tec-90.dat格式）5、以txt格式打开刚才命名的tec-90.dat，把里面的三列数复制到txt里面后再导入Excel里面，同时在Excel里面把szz的值变成相反的值保存在D列，（txt里面注意不是三列的都不要，只复制三列数据，把前面的和后面的都去掉）。
6、打开surfer，【网格】→【数据】，选中刚才生成的123.xls(即刚才生成的Excel数据，选中相应的A、B和D列),并单击打开。
对于【网格化数据】选项单则单击【确定】，则会生成与123.xls相对应的123.grd（即与刚才生成的Excel 文件名相同）。
7、在【地图】→【表面图】单击，在【打开网格对话框】中选中生成的123.grd文件，并单击打开。
8、右击图片选中【属性】并打击打开，可以对网格、颜色进行设置（颜色谱选安装文件夹中的rainbow.CLR即可,具体是在打开的【属性对话框】→【常规】→【材质颜色】→【上部】单击上部中的颜色，在颜色谱中选择【载入】，在打开里面选中安装文件夹里面的rainbow.CLR即可，在多次确认）9、此外也可以在surfer里面对其他参数进行设置。
深部大断面硐室围岩变形及控制技术
深部大断面硐室围岩变形及控制技术苏晓建【摘要】采用数值模拟方法分析了赵固二矿－800 m泵房围岩变形原因及控制技术。
认为深部大断面硐室围岩应力集中区距离硐室中心较远，但底鼓较为明显，支护方案设计时应重点关注。
基于此提出了主动支护与被动支护联合支护方案，首先锚网喷支护，其次双底拱刚性支架联合支护，第三围岩注浆加固。
工程应用表明，锚注＋刚性支架支护能有效控制围岩变形，满足泵房保持长期稳定的要求。
%Deformation mechanism and control technology of surrounding rocks of -800 m water pump house in Zhaogu No. 2 Coal Mine were analyzed by using numerical simulation. It was thought that the distance from stress concentration area to the center of chamber was rela-tively far,so the obvious floor heave should be the key point of supporting scheme. Based on the condition,combined supporting scheme of initiative and passive support was proposed which in-cluded bolt-mesh-spurting supporting,rigid support with double bottom arch and grouting rein-forcement. The result showed bolting-grouting and rigid support could effectively control the surrounding rock deformation and meet the requirement of chamber's long-term stability.【期刊名称】《中国煤炭》【年(卷),期】2017(043)001【总页数】4页(P75-78)【关键词】深部硐室;大断面;围岩控制;联合支护;数值模拟【作者】苏晓建【作者单位】中国矿业大学北京资源与安全工程学院，北京市海淀区，100083; 河南能源化工集团焦煤公司赵固二矿，河南省辉县市，453633【正文语种】中文【中图分类】TD353赵固二矿Ⅰ盘区－800 m泵房布置在二1煤层下方20 m左右的砂质泥岩中,埋深约900 m。
06-边坡变形-flac3d前后处理
• ini sxx=-5.0e6
• fix x range x -0.1 0.1
• fix x range x 5.9 6.1 • pause
• pl con sxx • solve • pl con sxx
• 注意是约束x方向
第三十八页，共84页。
第三十二页，共84页。
赋予初始应力
• gen zone brick size 6 6 6 • model elas • prop bulk 5e9 shear 3e9
• ini sxx=-5e6
• pause
• solve
只用ini，所有单元将有sxx应力， 但solve变为0，因为没有约束， 力引起位移最后平衡.
。 体）和Tetrahedral（四面体）单元
第二十四页，共84页。
第二十五页，共84页。
在ansys选择菜 单:
Input from
选择APDL命令流 文件(在仿真论坛 可以下载此程序)
执行后在电脑的D:\ 下生成ansys单元节 点信息的两个文件
01_node.dat
02_ele.dat
第二十六页，共84页。
• apply sxx=-5e6 range x 5.9 6.1 • pause
• pl con sxx • solve
• pl con sxx
• 注意是约束x方向---------变形参数并不影响应力场，它只影响位移场。如果考虑重力的时候， 变形参数影响泊松比，从而影响侧压力系数，才会影响水平的应力场分布。
平衡力大于容许值时，才会进行cyc。设置不同的变形参数并不影响位移结果，因为系统平衡了。 如果人为设置执行计算流程（如cyc 10000），由于位移边界条件约束，边界上位移为0，内 部则会慢慢随时间产生位移，但应力场是不会发生变化的。
Tecplot画等值线的操作步骤详细介绍00
FLAC3D计算结果—Tecplot画位移等值线，比较详细的操作步骤介绍！！！主要内容：详细介绍在FLAC3D中计算完毕后，把计算结果导入Tecplot中，绘制位移等值线的方法。
第一部分：Tecplot画位移等值线的具体操作步骤。
1. 运行FLAC3D程序，先调用FLAC3D计算结果的保存文件xxxx.sav，命令为：restore xxxx.sav然后，调用dynamax 编制的flac3d到tecplot的转换程序Flac3d2T ecplot.dat，命令为：call Flac3d2Tecplot.dat则在flac3d计算结果文件的当前目录下，生成一个名为tec10.dat文件，这个文件就是将flac3d计算结果转换成符合tecplot 10的数据格式后的文件，可供tecplot 10调用。
2. 在tecplot 10中，调用第一步生成的tec10.dat文件，将FLAC3D计算结果导入到tecplot 中，具体操作为：运行tecplot 10程序，点击File主菜单，在下拉菜单中选择Load Data File(s)选项，在弹出的Load Data File对话框中，选中上一步生成的tec10.dat文件，点击“打开”，之后一路点击默认选项，得到下图，3.），显示模型内部的一个指定剖面。
（相后，图形显示只选择on plane）4. 点击Iinteractive slicing tool*弹出如下对话框，显示切片工具的详细设置，先勾选Show Slices选项（显示切片）；再选择切片的方向（X-Planes和Y-Planes等分别表示法线方向平行于X轴和Y轴等方向的切片）；然后指定切片的位置，在Position后面输入数字，该数字为沿着所选切片的法线方向上的切片位置坐标。
切片设置好后，点击Close，关闭设置对话框。
5. 将切片截面上的数据读出，形成除整体zone外的另一格zone，操作如下，之后弹出如下对话框，点击Extract，弹出如下信息对话框，显示切片上数据提取成功。
tecplot生成等高线的步骤
1、先用FLAC把结果计算出来，可以保存几张图片（zdisp,xdisp,ydisp等），再接着Call“Flac3d2Tecplot位移应力.dat”这个文件，则生成了“tec10.dat”这个文件。
2、打开Tecplot10.3这个程序，点File>>Load data files选择“tec10.dat”文件。
在左边位置的zone下的“Contour”勾上、“Boundary”勾上,zone effects下的“Lighting”勾上,其余的不选。
3、“Plot”>>contour出现对话框，点“Show Lables”，可将等高线的数字显示出来，同时可设置字体的颜色和格式“Color,Font”
3、主菜单的“Data”>>Extract>>Slice from plane，出现对话框，在对话框中“Slice Plane”下的“constant”选择例如“X=0.5”就是X=0.5这个平面。
点“Extract”
4、点在左边的“Zone style”出现一个框，选择第一个“Tecplot v10”，点“zone show”选择“Deactive”；只让需要的面击活，其余的都再点后面的Deactive ，选择需要的那个面，在“Contour Type”选择“Line”，再点后面的“Line color”选择线的颜色。
5、点“Plot”>>Axis出现对话框，把“Show X-Axis”前的勾取消（X、Y、Z 三个都要取消）
6、点“View” >>“Rotate”出现对话框，点“Preset views”可选择视图的方向，点“Fit to full size”则可将图形填充在整个区域。
flac3d后处理
;Sk是草图sketch的缩写，作用使显示模型的外围网格线~
; ··变形矢量图··
;pl sk dis
;变形矢量图可以形象的表征变形的发生方向和相对大小
;pl sk dis scale 0.07
;scale关键字用来调整矢量线头的的大小
;pl block state
;显示塑性区域
;可在pl bl state 添加关键字
;pl bl state shear ; 获得剪切屈服的单元
;pl bl sta tension-p; 获得过去拉伸屈服的单元
;pl bl state now tension-n; 获得当前处于塑性状态的单元
;pl bl state past tension-p; 获得过去处于塑性状态的单元
new
gen zon b size 3 3 3
mod mohr
prop sh 1e6 coh 10e3 fric 15 bu 3e6
fix z ra z -.1 .1
fix x ra x -.1 .1
fix x ra x 2.9 3.1
fix y ra y -.1 .1
;
;save 6-2.sav
;施加额外的边界条件，
;为了更形象的描述模型变形后的结果，常常使用变形后的网格图命令如下
;··变形云图··
;pl con zd ou on magf 20
;magf 为放大倍数，不能缩写为mag，倍数自定，pl bcon不可以跟随 magf
;PL SK MAGF 20
;pl con szz max -10e3
;max 表示显示的应力数值最大不超过-10e3，默认为0
“用tecplot对flac3d计算结果进行后处理”相关讨论帖整理(第二部分)
第二部分：对dynamax大牛的原创主题帖“【原创】用tecplot对flac3d计算结果进行后处理”及相关讨论帖的整理和总结。
1. dynamax大牛的原创主题帖“【原创】用tecplot对flac3d计算结果进行后处理”的链接地址为：SimWe仿真论坛» D01：FLAC3D/FLAC2D » 【原创】用tecplot对flac3d计算结果进行后处理说明：该flac3d-tecplot转换程序，由dynamax大人利用fish编写，总共有三个版本，不断改进。
最初只能导入位移数据，画出位移等值线，并且不能剔除挖空单元（null），只能导入全部网格；第二版不再导出null zone，并且可以指定显示范围；第三版则可进一步同时导出位移和应力结果。
dynamax发表于2004-7-11 17:45[前后处理] 【原创】用tecplot对flac3d计算结果进行后处理fish写的因为flac3d中不能标注等值线图将位移计算结果导入tecplot中进行处理;可以绘出漂亮的等值线图熟悉tecplot的可以试一下不然得到转换后的数据文件还是不知道如何处理给出了简单的测试实例283356-tecplot.rar (45.85 KB) （第一版转换程序）dynamax切片功能，从菜单中选择，操作为：data = >extract = >slice from planedynamax上图显示成2D Cartesian可能更好要显示数值的话，plot-> Contour-> More-> Labels-> Show Labels坐标轴可以关掉plot-> axis-> show %-axisdynamax我给的例子只能导出全部网格的信息(针对第一版转换程序而言)如果要只导出开挖后剩下的部分因为节点编号的关系要另外做处理才行否则没有办法做出像下面这样的图我看代码有点不对，稍加修改应该没有问题，主要是把null部分的网格信息排除我们努力看看行不行。
FLAC3D基础知识介绍解析
FLAC 3D基础知识介绍一、概述FLAC（Fast Lagrangian Analysis of Continua）由美国Itasca公司开发的。
目前，FLAC有二维和三维计算程序两个版本，二维计算程序V3.0以前的为DOS版本，V2.5版本仅仅能够使用计算机的基本内存64K），所以，程序求解的最大结点数仅限于2000个以内。
1995年，FLAC2D已升级为V3.3的版本，其程序能够使用护展内存。
因此，大大发护展了计算规模。
FLAC3D是一个三维有限差分程序，目前已发展到V3.0版本。
FLAC3D的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入各种命令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。
因此，采用FLAC程序进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排列，形成可以完成一定计算任务的命令文件。
FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。
调整三维网格中的多面体单元来拟合实际的结构。
单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发生变形和移动（大变形模式）。
FLAC3D采用的显式拉格朗日算法和混合-离散分区技术，能够非常准确的模拟材料的塑性破坏和流动。
由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。
三维快速拉格朗日法是一种基于三维显式有限差分法的数值分析方法，它可以模拟岩土或其他材料的三维力学行为。
三维快速拉格朗日分析将计算区域划分为若干四面体单元，每个单元在给定的边界条件下遵循指定的线性或非线性本构关系，如果单元应力使得材料屈服或产生塑性流动，则单元网格可以随着材料的变形而变形，这就是所谓的拉格朗日算法，这种算法非常适合于模拟大变形问题。
三维快速拉格朗日分析采用了显式有限差分格式来求解场的控制微分方程，并应用了混合单元离散模型，可以准确地模拟材料的屈服、塑性流动、软化直至大变形，尤其在材料的弹塑性分析、大变形分析以及模拟施工过程等领域有其独到的优点。
裂隙岩体应力—损伤—渗流耦合理论、试验及工程应用研究
裂隙岩体应力—损伤—渗流耦合理论、试验及工程应用研究裂隙岩体内部富含各种缺陷,包括微裂纹、孔隙以及节理裂隙等宏观非连续面,这些存在的缺陷不但大大地改变了岩体的力学性质,而且也严重影响着岩体的渗透特性,因而裂隙岩体具有复杂的力学特性和渗透特性。
应力-损伤-渗流耦合是指裂隙中的渗流水压力加剧岩体裂隙的起裂、扩展、贯通,导致岩体中应力场和损伤场的改变；而岩体应力的改变和岩体裂隙的损伤扩展,又导致裂隙岩体渗透特性变化,进而改变渗流场的分布。
本文重点研究了裂隙岩体的应力-损伤-渗流耦合过程,揭示了岩体裂隙在应力-损伤-渗流耦合状态下裂隙的起裂、扩展、贯通规律,探求了裂隙岩体在应力-损伤-渗流耦合状态下由于其裂隙的扩展贯通而导致破坏的机理,进而揭示了含水裂隙岩体失稳破坏的原因。
该耦合理论在预防采矿工程中的矿井突水、水利水电工程中的大坝失稳破坏以及水下隧道工程突(涌)水等方面有着广泛的应用前景。
本文结合国家重点基础研究发展973计划项目(2007CB209402：矿井突水的动力学特征及控制因素)和山东省黄金集团相关科研项目的共同资助下开展理论和应用研究,主要开展的研究工作如下：(1)理论分析研究：将断裂力学和损伤力学相结合,对渗透水压作用下裂隙岩体的损伤断裂力学性质进行研究；探讨压剪、拉剪应力状态下及渗透水压力作用下裂隙岩体中张开型裂纹闭合、起裂、分支裂纹扩展、相邻裂纹间岩桥断裂破坏及贯穿等的规律；给出渗透水压力-压(拉)剪应力共同作用下岩石裂隙变化发展规律。
(2)力学试验研究：设计制作含预置裂纹圆柱体岩石试样(φ50mm×100mm)的模具,可根据需要制作单、双预置裂纹的类岩石试件；在原有2000KN伺服试验机设备的基础上,设计开发三轴室、水压加载系统和围压加载系统,改进后的试验装置可完成完整类岩石试件应力应变-渗透率试验和渗透水压-应力耦合作用下含预置裂隙岩石试件破坏过程的试验；进行含预置裂隙类岩石试件的常规单轴试验并使用声发射仪对试验过程进行监测,并用离散元软件PFC对单轴试验过程进行数值模拟；运用RFPA-flow模拟分析在渗透水压力作用下,含预置裂隙岩石的破坏过程,并与试验结果进行了对比,通过应力场和渗流场的图形显示观察到岩石试件中裂纹的萌生、扩展、贯通直至整个结构破裂的全过程。
FLAC3D基础知识
FLAC3D基础知识FLAC 3D基础知识介绍一、概述FLAC（Fast Lagrangian Analysis of Continua）由美国Itasca 公司开发的。
目前，FLAC有二维和三维计算程序两个版本，二维计算程序V3.0以前的为DOS版本，V2.5版本仅仅能够使用计算机的基本内存64K），所以，程序求解的最大结点数仅限于2000个以内。
1995年，FLAC2D已升级为V3.3的版本，其程序能够使用护展内存。
因此，大大发护展了计算规模。
FLAC3D是一个三维有限差分程序，目前已发展到V3.0版本。
FLAC3D的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入各种命令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。
因此，采用FLAC程序进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排列，形成可以完成一定计算任务的命令文件。
FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。
调整三维网格中的多面体单元来拟合实际的结构。
单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发生变形和移动（大变形模式）。
FLAC3 D采用的显式拉格朗日算法和混合-离散分区技术，能够非常准确的模拟材料的塑性破坏和流动。
由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。
三维快速拉格朗日法是一种基于三维显式有限差分法的数值分析方法，它可以模拟岩土或其他材料的三维力学行为。
三维快速拉格朗日分析将计算区域划分为若干四面体单元，每个单元在给定的边界条件下遵循指定的线性或非线性本构关系，如果单元应力使得材料屈服或产生塑性流动，则单元网格可以随着材料的变形而变形，这就是所谓的拉格朗日算法，这种算法非常适合于模拟大变形问题。
三维快速拉格朗日分析采用了显式有限差分格式来求解场的控制微分方程，并应用了混合单元离散模型，可以准确地模拟材料的屈服、塑性流动、软化直至大变形，尤其在材料的弹塑性分析、大变形分析以及模拟施工过程等领域有其独到的优点。
flac-3d的常用后处理命令
flac-3d的常用后处理命令flac-3d的常用后处理命令•;包括各种应力，位移的云图，塑性区显示和图像的存储，矢量图显示，历史变量显示及其txt文件存储，;应力和位移结果输出及其log文件的存储，切片云图显示，求解过程的动画设置和存储newgen zon bri size 3 3 3model mohrprop bu 3e6 sh 1e6 coh 10e3 fric 15fix z ran z -0.1 0.1fix x ran x -0.1 0.1fix x ran x 2.9 3.1fix y ran y -0.1 0.1fix y ran y 2.9 3.1ini dens 2000hist unbaset grav 10solve elasticsave 6-1.savrest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2hist id=2 gp zdis 1 1 3hist id=3 gp zdis 1 1 2hist id=4 gp xdis 1 1 3hist id=6 zone szz 1 1 3hist id=7 zone szz 1.5 1.5 2.5hist id=8 zone sxz 1.5 1.5 2.5set hist_rep 1 ;历史变量存储间隔设置solvesave 6-2.savplot hist 2plot hist 7 v 8 ;变量7和8对比显示plot set rot 20 0 30plot set mag 0.8set plot jpgset plot quality 20plot hard file 003.jpghist write 7 v 8 file 7v8hist.txt ;变量写如数据文件plot surplot block group ;显示块体分组情况plot con szz ou on magf 20plot con zdisp ou on magf 20plot add axeplot set rot 20 0 30 ;输出前设置视图角度plot set mag 0.8set plot jpg ;jpg格式的图形输出和设置set plot quality 20plot hard file 001.jpgplot block group id on ;显示单元的id的网格模型plot con smaxplot con smin ;显示最大、最小、中间主应力云图plot con smidplot con szz ef on ou on ;输出有效应力plot add axeplot set rot 20 0 30plot set mag 0.8setplot jpg set plot quality 20plot hard file 002.jpgplot bcon prop bu ;块云图应用之检查模型参数的赋值情况plot sk magf 20 ;显示模型外围网格线plot sk dis scale 0.07 ;矢量图显示plot add axeplot block state ;显示块体的塑性区分布plot block state nowplot block state pastplot block state shear-nset log on ;通过日志文件输出单元和节点结果数据set logfile yingli&weiyi.logprint zone stressprint gp dishist write 2 3 4 5 6 vs 7 file 1.datset log offplot set plane ori 0 1.5 0 norm 0 1 0 ;设置切面切面中一点坐标ori和切面的法向norm确定plot con zd planeplot add ske ;增加网格线plot add dis plane ;增加矢量plot add axe ;增加坐标系plo set pla ori 1.5 1.5 1.5 norm 1 1 1 ;剖面图显示状态下可以直接更改剖面位置rest 6-1.sav ;动画的设置ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2plot set rot 20 0 30plot con szz ou on magf 10plot add hist 1set movie avi step 1 file szz&unforce.avimoviestartsolvemoviefinishsave 6-3.sav。
FLAC3D基础知识介绍
FLAC 3D基础知识介绍一、概述FLAC（Fast Lagrangian Analysis of Continua）由美国Itasca公司开发的。
目前，FLAC有二维和三维计算程序两个版本，二维计算程序V3.0以前的为DOS版本，V2.5版本仅仅能够使用计算机的基本内存64K），所以，程序求解的最大结点数仅限于2000个以内。
1995年，FLAC2D已升级为V3.3的版本，其程序能够使用护展内存。
因此，大大发护展了计算规模。
FLAC3D是一个三维有限差分程序，目前已发展到V3.0版本。
FLAC3D的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入各种命令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。
因此，采用FLAC程序进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排列，形成可以完成一定计算任务的命令文件。
FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。
调整三维网格中的多面体单元来拟合实际的结构。
单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发生变形和移动（大变形模式）。
FLAC3D采用的显式拉格朗日算法和混合-离散分区技术，能够非常准确的模拟材料的塑性破坏和流动。
由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。
三维快速拉格朗日法是一种基于三维显式有限差分法的数值分析方法，它可以模拟岩土或其他材料的三维力学行为。
三维快速拉格朗日分析将计算区域划分为若干四面体单元，每个单元在给定的边界条件下遵循指定的线性或非线性本构关系，如果单元应力使得材料屈服或产生塑性流动，则单元网格可以随着材料的变形而变形，这就是所谓的拉格朗日算法，这种算法非常适合于模拟大变形问题。
三维快速拉格朗日分析采用了显式有限差分格式来求解场的控制微分方程，并应用了混合单元离散模型，可以准确地模拟材料的屈服、塑性流动、软化直至大变形，尤其在材料的弹塑性分析、大变形分析以及模拟施工过程等领域有其独到的优点。
FLAC3D后处理命令集
;0.)实体建模的显示plot set rotation 0 0 45plot block group ;显示块体plot add axes redplot add surface yellowplot show;(1.形成切面的操作);以0为中心取切面plot set plane ori 0 0 0 norm 0 1 0 ;输出szz方向的云图plot con szz plane;添加网格线plot add ske;在切面上添加位移矢量plot add dis plane;在切面上增加坐标系plot add axe;在切面上绘制竖向位移云图plot con zd plane;(2.结果输出命令:print(输出关键词)，gp（节点），zone(单元)，history （监测历史）。
);单元应力结果输出print zone stress;节点变形的信息print gp dis;存储节点信息(可以用记事本打开.log)set log onset 文件名.logprint zone stressprint gp disset log off;(3.变形云图）;生产变形矢量图plot sk dis;输出变形网格plot sk magf 20;(变形倍数);变形后的网格与应力云图结合起来plot con szz ou on magef 20;（4.形成一半的剖面图plot set plane ori 0 0 0 norm 0 1 0plot add surf plane behind yell;在剖面后显示位移云图plot add cont disp plane behind shade on ;;显示单元的id的网格模型plot block group id on;显示最大、最小、中间主应力云图plot con smaxplot con sminplot con smid。
FLAC3D常用命令总结
FLAC3D常用命令总结1、模型建好如何检验所建模型是否有问题注意：在建模的第一步，我们不能建立太大的模型。
我们可以建立非常小的模型，以确保模型的准确性。
这样，我们就可以尽快开始后面的步骤。
在后面的计算和分析步骤快速完成且没有重大问题后，我们可以替换精确的模型并重新计算，以获得更精确的结果。
2.如何检查应力边界、位移边界和加载速率是否准确plotaddgpfix：看位移边界Ploaddfap：检查应力边界和加载方向是否准确。
Plotaddvel：检查速率和加载速率是否准确。
3.如何检查计算结果ploconszzoutlineon：看z向应力云图，并显示边界网格ploadddisp：看位移Ploconzdispoutlineonmagf5：看看z方向的位移云，变形扩大了5倍。
Ploaddbo:添加边界ploblostate:查看是否屈服ploaddstensor：查看大主应力方向ploconsmin：小主应力ploconsmax：大主应力Ploconssi（SurristaReCuMe）：查看剪切应变增量云图：查看单元4的小主应力。
如何制作剖面并查看剖面结果在编辑“edit“菜单中设置切面位置，在上述3命令后加上plane（pl）即可。
如：ploconsminpl：查看切面小主应力5、后处理将模型导入Tecplot制作剖面、添加等值线等，并进一步处理结果以获得更漂亮的结果。
其他重要命令总结：1.如何在圆柱体周围应用约束？可以用fix...rancylinderend1end2radiusr1cylinderend1end2radiusr2not，其中r22、如何能把一个plot的图像数据导出来以便用其他软件绘图？答：用setlogon命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。
3、绘制相关曲线命令的一些解释。
1)在plothistmvsn的形式里,m代表y轴,n代表x轴(不管m,n的正负)；2)\表示对其值作\镜像。
Tecplot画等值线的操作步骤详细介绍(第一部分)
FLAC3D计算结果—Tecplot画位移等值线，比较详细的操作步骤介绍！！！本文主要内容：(1)详细介绍在FLAC3D中计算完毕后，把计算结果导入Tecplot中，绘制位移等值线的方法。
(2)对dynamax大牛的原创主题帖“【原创】用tecplot对flac3d计算结果进行后处理”及相关讨论帖的整理和总结。
第一部分：Tecplot画位移等值线的具体操作步骤。
1. 运行FLAC3D程序，先调用FLAC3D计算结果的保存文件xxxx.sav，命令为：restore xxxx.sav然后，调用dynamax 编制的flac3d到tecplot的转换程序Flac3d2Tecplot.dat，命令为：call Flac3d2Tecplot.dat则在flac3d计算结果文件的当前目录下，生成一个名为tec10.dat文件，这个文件就是将flac3d计算结果转换成符合tecplot 10的数据格式后的文件，可供tecplot 10调用。
2. 在tecplot 10中，调用第一步生成的tec10.dat文件，将FLAC3D计算结果导入到tecplot 中，具体操作为：运行tecplot 10程序，点击File主菜单，在下拉菜单中选择Load Data File(s)选项，在弹出的Load Data File对话框中，选中上一步生成的tec10.dat文件，点击“打开”，之后一路点击默认选项，得到下图，3.），显示模型内部的一个指定剖面。
（相后，图形显示只选择on plane）4. 点击Iinteractive slicing tool之后，再点击下面的，弹出如下对话框，显示切片工具的详细设置，先勾选Show Slices选项（显示切片）；再选择切片的方向（X-Planes和Y-Planes等分别表示法线方向平行于X轴和Y轴等方向的切片）；然后指定切片的位置，在Position后面输入数字，该数字为沿着所选切片的法线方向上的切片位置坐标。
基于FLAC3D的近距离极薄煤层开采数值模拟研究
基于FLAC3D的近距离极薄煤层开采数值模拟研究杨松;杨邦友【摘要】通过FLAC3D数值模拟软件模拟研究了近距离倾斜薄煤层开采引起的工作面上覆岩层运移形态规律。
得出:工作面回采过程中,工作面顶板受开采卸荷作用的影响,顶板出现垂直向下的卸载膨胀变形,引起顶板的离层、冒落。
底板出现垂直向上的卸荷膨胀变形,煤层顶板变形大于底板;随着工作面的不断推进,煤层回采后,工作面前后方出现应力集中现象,随着开采进行,工作面前后方支承应力是移动型支承应力;煤层回采后,地表变形逐渐增大,且最大沉降与煤层倾角有关,约在煤层中部偏采区深部一侧。
【期刊名称】《矿山工程》【年(卷),期】2019(007)001【总页数】7页(P69-75)【关键词】数值模拟;应力集中;地表沉降【作者】杨松;杨邦友【作者单位】[1]川煤集团芙蓉公司四川省叙永煤矿,四川泸州;[2]重庆大学,资源及环境科学学院,重庆;【正文语种】中文【中图分类】TD31.引言在地下开采过程中，采空区应运而生。
开采过程破坏了开采区域的原岩应力状态。
随着开采推进，采场围岩因为应力改变发生变形破坏。
采场变形破坏引起地表沉陷，进而可能会对地表建筑产生不利影响。
空区治理在矿山的安全生产中具有重大的意义。
基于拉格朗日算法的FLAC3D程序特别适合于求解岩土力学工程中非线性的大变形问题，这点是非常适合采矿工程的设计研究工作的。
张海波等[1]以石人沟铁矿为工程背景，运用FLAC3D对采空区的稳定性进行了数值模拟分析，并得到空区围岩的应力应变分布规律。
黄庆享[2]根据不同条件的浅埋煤层工作面矿压实测结果，得出了浅埋矿区顶板破断规律；武崇福等[3]通过观测采空区破坏的形式及采空区破坏的特征，为治理采空区提供了理论依据。
宋广朋等[4]分析倾斜煤层回采巷道矿压显现规律研究。
庞德等[5]通过理论分析以及FLAC3D数值模拟的方法，总结了谢桥矿煤层开采引起的覆岩和地表沉陷移动规律，为预测地表沉陷提供依据。
用Tecplot对Flac3D计算结果进行后处理
⽤Tecplot对Flac3D计算结果进⾏后处理
⽤Tecplot对Flac3D计算结果进⾏后处理
第⼀步：⽤FLAC3D打开已经完成的*.sav⽂件；
第⼆步：然后【File】/【Call】运⾏⽂件，运⾏结束时会在命令⾏⾥显⽰写⼊成功，在命令⽂件所在的⽂件夹⾥会发现数据⽂件（可以根据⾃⼰的需要更改⽂件名）。
第三步：打开Tecplot程序，执⾏【File】/【Load Data Files(s)】命令打开数据⽂件。
以下是画位移等值线图的步骤：
1.三维图形选择3D Cartesian，⼆维选择2D Cartesian；
2.在边栏左上⽅Zone Layers栏⽬中,去掉Mesh前的√，选中Contour前的√。
然后按下
more按钮，在Var选择⾃⼰需要显⽰的量：
3.选择切⽚⼯具，并打开下⾯的Tool Details，可以设置切⽚的⽅向和位置等；
4.将切⽚截⾯上的数据读出，形成整体zone外的另⼀个zone；
5.点击zone style，编辑zone；
6.点击surfaces，发现有两个zone，在zoneshow中激活Slice：z=0，关掉tecplot v10；
7.点击contour，contour type选取lines；
8.然后再等值线上显⽰标注；
9.然后可以进⾏去掉颜⾊、⽹格等编辑。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
用Tecplot对Flac3D计算结果进行后处理
第一步：用FLAC3D打开已经完成的*.sav文件；
第二步：然后【File】/【Call】运行Flac3d2Tecplot.dat文件，运行结束时会在命令行里显示写入成功，在命令文件所在的文件夹里会发现tec10.dat数据文件（可以根据自己的需要更改文件名）。
第三步：打开Tecplot程序，执行【File】/【Load Data Files(s)】命令打开tec10.dat数据文件。
以下是画位移等值线图的步骤：
1.三维图形选择3D Cartesian，二维选择2D Cartesian；
2.在边栏左上方Zone Layers栏目中,去掉Mesh前的√，选中Contour前的√。
然后按下
more按钮，在Var选择自己需要显示的量：
3.选择切片工具，并打开下面的Tool Details，可以设置切片的方向和位置等；
4.将切片截面上的数据读出，形成整体zone外的另一个zone；
5.点击zone style，编辑zone；
6.点击surfaces，发现有两个zone，在zoneshow中激活Slice：z=0，关掉tecplot v10；
7.点击contour，contour type选取lines；
8.然后再等值线上显示标注；
9.然后可以进行去掉颜色、网格等编辑。