FLAC3D建模方法探讨

FLAC3D建模方法探讨

摘要：针对FLAC3D前处理方面的功能不足，本文列举了目前FLAC3D建模的常用方法以及其建模思路，供需要者参考和选择，具有一定的实用价值。

关键词：FLAC3D；建模方法

1.引言

随着计算机技术的发展，数值模拟方法已广泛应用于岩土及地下工程的研究和设计中，FLAC3D是三维快速拉格朗日法(Fast Lagrangian Analysis of Continua-3D)的缩写，它是一种基于三维显式有限差分法的数值分析方法，可以模拟岩土或其他材料的三维力学行为[1]。这种分析方法非常适合于模拟大变形问题，尤其在材料的弹塑性分析、大变形分析以及模拟施工过程等领域有其独到的优点，但FLAC3D非完全可视化的建模方式和网格划分方法过于繁琐，限制了其通用性[2]。

2.建模方法探讨

随着FLAC3D的应用逐步普遍，出现了FLAC3D的一些建模的方法，在此列举出来，以供需要者参考和选择。

2.1 FLAC3D中INTERFACE建模方法

INTERFACE单元在数字模拟中应用相当广泛，很多数值模拟软件都包含有这种单元。它可以模拟结构的接缝，不同材料、不同块体的接触面，以及岩石的节理面等。

在FLAC3D中使用INTERFACE建模，目前来说基本上有三种方法，分离法、移动组合法和依次建模法。

分离法的思路是先把块体单元分组建好，然后将界面上两部分的接触网格分离为两部分，最后设置INTERFACE单元。

移动组合法的实现步骤是先在不同的位置分别建立两组块体，然后在一个块体的表面建立INTERFACE，最后通过移动使两组块体在INTERFACE处组合起来。

依次建模法是建立一组块体，然后在其相应表面建立INTERFACE，再建立与之接触的另外一组块体。对于复杂的模型，交替建立块体和INTERFACE[3]。

2.2 基于Midas／Gts的FLAC3D的建模方法

flac3d常用命令

1、最先需要掌握的命令有哪些？ 答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。 2、怎样输出模型的后处理图？ 答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件。 3、怎样调用一个文件？ 答：File/call或者call命令 4、如何施加面力？ 答：app nstress 5、如何调整视图的大小、角度？ 答：综合使用x, y, z, m, Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键。 6、如何进行边界约束？ 答：fix x ran （约束的是速度，在初始情况下约束等效于位移约束）。 7、如何知道每个单元的ID？ 答：用鼠标双击单元的表面，可以知道单元的ID和坐标。 8、如何进行切片？ 答：plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片) 9、如何保存计算结果？ 答：save +文件名 10、如何调用已保存的结果？ 答：rest +文件名；或者File / Restor 11、如何暂停计算？ 答：Esc 12、如何在程序中进行暂停，并可恢复计算？ 答：在命令中加入pause命令，用continue进行继续。 在我们分步求解中想得到某一个过程中的结果，不用等到全求完，还可以在分布求解错误的时候就进行改正，而不是等到结果出来。 13、如何跳过某个计算步？ 答：在计算中按空格键跳过本次计算，自动进入下一步 14、Fish是什么东西？Fish是否一定要学？

答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能。Fish可以不用学，需要的时候查Mannual获得需要的变量就可以了。 15、FLAC3D允许的命令文件格式有哪些？ 答：无所谓，只要是文本文件，什么后缀都可以。 16、如何调用一些可选模块？ 答：config dyn (fluid, creep, cppudm) 17、如何在圆柱体四周如何施加约束条件？ 可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2

flac3d建模方法

利用FLAC3D 进行数值分析的第一步便是如何将物理系统转化为由实体单元和结构单元所组合的网格模型（Modeling ），该模型与分析对象的几何外形特征相一致。目前，FLAC3D 网格模型的建立方法可分为两种，即直接法及间接法，直接法是按照分析对象的几何形状利用FLAC3D 内置的网格生成器建模，网格和几何模型同时生成，该方法较适用于简单几何外形的物理系统；与之不同，间接法则适用于复杂的、单元数目较多的物理系统，该方法建立网格模型时，像一般计算机绘图软件一样，通过点、线、面、体，先建立对象的几何外形，再进行实体模型的分网（Meshing ），以完成网格模型的建立，FLAC3D 自身不具备间接法建模功能，读者可借助第三方软件与FLAC3D 的接入轻松实现。本章主要介绍FLAC3D 的网格建模方法，包括利用网格生成器建立简单网格、利用第三方软件进行模型导入以及复杂模型的方法。 本章要点： z FLAC3D 网格单元的基本类型 z 网格的连接 z FLAC3D 网格的数据格式 z 常用有限元模型与FLAC3D 的接入 z 复杂模型的建立 5.1 简单网格的建立 5.1.1 基本网格的形状 FLAC3D 内置网格生成器中的基本形状网格有13种，通过匹配、连接这些基本形状网格单元，能够生成一些较为复杂的三维结构网格。网格单元的基本类型和特征如表5-1所示，基本可以归为四大类，即六面块体网格、退化网格、放射网格和交叉网格。 5 FLAC3D 建模方法

表5-1 FLAC3D 基本形状网格的基本特征

5.1.2 单元网格的生成 生成块体网格（Brick ）的命令格式如下： generate zone brick p0 x0 y0 z0 p1 x1 y1 z1 …… p7 x7 y7 z7 size n1 n2 n3 ratio r1 r2 r3 或者 generate zone brick p0 x0 y0 z0 p1 add x1 y1 z1 …… p7 add x7 y7 z7 size n1 n2 n3 ratio r1 r2 r3 在该命令中，generate 为“生成网格”之意，可以缩写为gen ，zone 表示该命令文件生成的是实体单元，brick 关键词表明建立的网格采用的是brick 基本形状，p0，p1……p7是块体单元的8个控制点，其后跟这些点的三维坐标值（xn, yn, zn ），含义是由8个点可确定一个六面体网格。不过，p0～p7各点的定义需遵从“右手法则”，不能随意颠倒顺序。如果采用全局坐标系，三维坐标值应为建模空间内的全局三维坐标值；若采用局部坐标系，则除p0点采用全局三维坐标值外，其他点的坐标值都必须取其相对于点p0的三维坐标值，且在点编号后加关键词add （见本节第2行命令）。size 为定义坐标轴（x ，y ，z ）方向网格单元数目的关键词，其后跟划分的单元数目（n1，n2，n3）；ratio 为定义相邻单元尺寸大小比率的关键词，其后跟坐标轴方向相邻网格单元的比率（r1，r2，r3）。 如果生成的是长方体网格，前述命令可以简化为： generate zone brick p0 x0 y0 z0 p1 x1 y1 z1 p2 x2 y2 z2 p3 x3 y3 z3 size n1 n2 n3 ratio r1 r2 r3 或者 generate zone brick p0 x0 y0 z0 p1 add x1 y1 z1 p2 add x2 y2 z2 p3 add x3 y3 z3 & size n1 n2 n3 ratio r1 r2 r3 即只需采用4个控制点即可确定该长方体。 此外，当网格的几何形状为立方体时，上述命令文件可以用下列命令替代，进一步简化，关键词edge 后跟的evalue 是立方体的边长。 generate zone brick p0 x0 y0 z0 edge evalue size n1 n2 n3 ratio r1 r2 r3

FLAC3D学习笔记(自己总结版)

一、建立模型笔记 1键盘上的键的功能:X,Y,Z分别控制所建立的模型围绕X,Y,Z三个轴做逆时针旋转.如果打开大小写锁定键,分边按XYZ时建立的模型会做顺时针旋转.+号的功能是加大模型旋转的角度，-号的功能是减小模型旋转的角度。 2键盘上的上下左右四个键可以实现图形的向上向下,向左和向右的移动. 3当在一个区域建立多个模型的时候,有时候为了镜像单个模块儿,需要预先编组,编组的方法如下: gen zon tunint p0 0 0 0 p1 0 10 0 p2 5 5 0 p3 0 0 -5 p4 5 10 0 & p5 5 5 -5 p6 0 10 -5 p7 5 10 -5 dim sq2 sq2 1 sq2 1 1 1 group(编组) z1(组名为:z1) & fill group tunnel 4镜像整个模型的方法:首先判断即将镜像的模型在原模型的哪个方向(即三个轴的方向),然后在轴所在的方向上去两个点.其中一个是原点(origin(xyz)). 5镜像部分模型的方法: gen zone reflect normal (x,y,z) origin (x,y,z) range group z1就可以了. 6特别要注意的一点:交叉巷道的巷道充填和巷道的网格是两个组,所以在镜像的时候不要忘了给充填部分镜像. 7 group radcyl range group cylint not(编组的时候不包括cylint) 7 快捷键总结 Ctrl+Z 局部缩放的功能。Ctrl+R的功能是恢复到初始状态。 8模型建立的方法: 方法1:利用默认值生成网格, 各个默认值:(1)p0的默认值是(0,0,0),网格的每边的单元格数size默认值是10,网格每边的长度默认值为size的值. 方法2:利用4个点的坐标来生成矩形网格.p0~p3 size 方法3:利用edge来确定边长生成矩形网格. 方法4:利用参数ratio来确定单元体几何边画边生成矩形网格. 方法5:利用参数add(相对坐标)来生成矩形网格. 方法6:利用八个点的左边来生成矩形网格. 9、用户自定义模型的方法 我来贡献一点自己的成果： FLAC3D的二次开发环境提供了开放的用户接口，在软件安装文件中包含了软件自带所有本构模型的源代码，且给出了Mohr-Coulomb模型和应变软化模型的编译示例，因此可以方便地进行本构模型的修改与开发。为了方便起见，下面的说明以建立UserModel模型为例。 (1) 在模型头文件(usermodel.h)中进行新的本构模型派生类的声明，修改模型的ID(为避免与已有模型冲突，一般要求大于100)、名称和版本，修改派生类的私有成员，主要包括模型的基本参数及程序执行过程中主要的中间变量。( ? p6 u' J5 Q3 y( a (2) 在程序C++文件(usermodel.cpp)中修改模型结构(UserModel::UserModel(bool bRegister): Constit -utive Model)的定义，这是一个空函数，主要功能是给(1)中定义的所有私有成员赋初值，一般均赋值为0.0。(3) 修改const char **UserModel: roperties()函数，该函数包含了给定模型的参数名称字符串，在FLAC3D的计算命令中需要用到这些字符串进行模型参数赋值。 (4) const char **UserModel::States()函数是单元在计算过程中的状态指示器，可以按照需要进行修改指示器的内容。 (5) 按照派生类中定义的模型参数变量修改double UserModel::GetProperty()和void UserModel:: SetProperty()函数，这两个函数共同完成模型参数的赋值功能。' U; e G' W" Q# R4 q/ @9 G" h (6) const char * UserModel::Initialize()函数在执行CYCLE命令或大应变模式下对于每个模型单元(zone)调用一次，主要执行参数和状态指示器的初始化，并对派生类声明中定义的私有变量进行赋值。值得注意的是，Initialize()函数调用时没有定义应变分量，但可以调用应力分量，但不能对应力进行修改。9 n# e8 |' c- B/ q, B B5 m

FLAC3D建模(复杂边坡)问题

FLAC3D建模（复杂边坡）问题 总在论坛里面看帖子，学习到了很多建模的技巧，觉得基本上是通过ansys与surfer导网格到FLAC3D进行建模，这其中就存在一些问题，如果模型比较复杂的话，可能会严重失真，例如边坡（我是岩土工程专业的，比较关注边坡），如果边坡内部存在多条断层，且存在崩积物、坡积物叠加分布的现象，我们假设存在基岩上覆盖10层崩坡积物，那用ansys与surfer 导入的网格是否能够真实反映边坡的真实情况呢？（我没有具体用ansys与surfer做过，所以并不知道是否可行！），最近做的一个边坡，强风化基岩上存在15层叠加分布的崩坡积物，且边坡内部存在3条断层，考虑了很久，觉得用ansys与surfer做的话可能不合适（如果有朋友能做的话，请告诉我，先谢谢了！），最后考虑用AutoCAD建模，然后导入FLAC3D的方法，模型建立起来后，觉得效果还不错，下面就说下我的方法，有什么不恰当之处，还请各位批评指正。 1、首先在AutoCAD中建模，因为我觉得我们的剖面图基本上是AutoCAD绘制的，尽量就CAD建模的话可以省下不少事（个人观点），建立3维模型需要尽量多的剖面，因此，在边坡表面地形变化处、坡体内部不同成因岩土体处都需要做剖面； 2、剖面做好后，进行剖面分层，就是在剖面上尽量细致的划分出不同成因性质的岩土体层，例如崩积物、坡积物，之后利用AutoCAD的3dmesh功能进行网格划分； 3、基本上只需要brick与wedge就可以了，自己编写一个VB的程序，可以读出3dmesh网格节点并输出成*.dat文件； 例如：Z方向为单位长度，读取3dmesh（brick）程序 Dim AcadApp As AutoCAD.AcadApplication Set AcadApp = GetObject(, "AutoCAD.Application") Dim Mydocument As AcadDocument Set Mydocument = AcadApp.ActiveDocument Dim Myentity As AcadPolygonMesh Dim Mysel As AcadSelectionSet Dim fil_type(0) As Integer Dim fil_data(0) As Variant Dim Mycoordinates As Variant

FLAC3D和3DEC的区别(itasca论坛)

FLAC3D和3DEC的区别（itasca论坛） 在“三维可视化”专区中一位坛友讲到了3DEC建模和参数取值困难的问题，且被业界所“诟病”。楼主是3DEC的忠实用户，也用过FLAC3D和PFC，就3DEC自身的这些“问题”，略谈一二。 在建模环节上，3DEC是典型的“欺生”，即生手上来的时候可能觉得无从下手，而掌握以后觉得游刃有余。有兴趣的坛友可以访问Itasca网站中咨询研究专题下的一些应用实例，看看那里的3DEC模型，所有这些模型，除其中一个以外，建模时间都在5天以内，一般为3天，不助任何第三方软件。大家可以比较一下，这些模型用其他软件建模时需要花多长时间。 与其他ITASCA软件如FLAC3D一致地，3DEC采用命令流建模，可以嵌入FISH。这对初学者而言有些困难，但如果掌握了就很灵活： 相比较菜单操作而言，命令流显然要困难一些，因此难学。但一旦掌握，修改模型只需要修改命令流，无需进行重复的菜单操作，熟手因此多喜欢命令流的方式； FISH本身就不是为初学者准备的，但掌握了FISH以后，任何有规律性对象的建模就容易得多。比如多机组的地下厂房洞室群，用FISH建立起来了其中一个单元以后，修改参数即可很快获得其他单元。建模过程中FISH的应用可以大大提高效率，是命令流、特别是菜单操作所难以比拟的。 因此，3DEC是典型的“欺生”，学习3DEC和训一匹“烈马”有得

一比，如果你不能驯服它，你只能是望而却步。 学习3DEC还需要“洗脑”，即固有的理论和思维模式可能会有所障碍。 不知道哪位坛友在基础理论学习过程中系统学习过非连续力学方法、或者说学习过离散元方法。以楼主的理解，国内一些关于离散元的文献中，只要谈深一点，就不乏误解和错误。其中的原因是多方面的，一是没有跟上国际潮流，没有真正弄懂就开始写文章。二是既往的教育背景基本都是连续理论和有限元方法，当用这种思维定式来应用离散元程序如3DEC时，可能会受到一些制约。 3DEC处理的对象包括两大部分，即块体和接触（结构面），相比较而言，有限元乃至FLAC3D针对的对象只是其中的块体，即便存在结构面，但在程序结构中的处理方式也与3DEC有着本质的差别。 如果说单元网格是数值模型的基本单元，那么，3DEC中包括两种网格体系，即接触网格和块体网格。这要求用户在脑子里建立这种概念，即3DEC程序结构中有两个系列的网格，而不是传统的一个系列。当然，从应用角度，如果不理会这一点也无所谓，但如果是达到写文章和理解程序的深度，则需要理解这种差别及其带来的不同。 在3DEC中，接触被作为块体的边界，因此，两个块体相接触时，无需节点之间有任何的对应关系。在FLAC专区中有一幅帖子问“这样的网格行不行”，如果放在3DEC中，这不是问题。这一差别显

FLAC3D50模型及输入参数说明

1.1模型参数代码 可参考manual中各个章节的command命令及说明,注意单位。用prop 赋值。

16 ttable 塑性拉应变-抗拉强度的表号 下列参数可以显示、绘图与通过fish访问 1 es_plastic 塑性切应变 2 et_plastic 塑性拉应变 3 ff_count 检测切应变反向的数 4 ff_cvd 体应变,εvd 经典粘弹性模型的材料参数(Classical Viscoelastic (Maxwell Substance) –MODEL mechanical viscous) 1 bulk 弹性体积模量,K 2 shear 弹性剪切模量,G 3 viscosity 动力粘度,η 粘弹性模型的材料参数(Burgers Model –MODEL mechanical burgers) 1 bulk 弹性体积模量,K 2 kshear Kelvin弹性剪切模量,G K 3 kviscosity Kelvin动力粘度,ηK 4 mkshear Maxwell切边模量,G M 5 mviscosity Maxwell动力粘度,ηM 二分幂律模型的材料参数(Power Law –MODEL mechanical power) 1 a_1 常数,A1 2 a_2 常数,A2 3 bulk 弹性体积模量,K 4 n_1 指数,n1 5 n_2 指数,n2 6 rs_1 参考应力,σ1ref 7 rs_2 参考应力,σ2ref 8 shear 弹性剪切模量,G 蠕变模型材料参数(WIPP Model –MODEL mechanical wipp) 1 act_energy 活化能,Q 2 a_wipp 常数,A 3 b_wipp 常数,B 4 bulk 弹性体积模量,K 5 d_wipp 常数,D 6 e_dot_star 临界稳定状态蠕变率, 7 gas_c 气体常数,R 8 n_wipp 指数,n 9 shear 弹性剪切模量,G 10 temp 温度,T 下列参数可以显示、绘图与通过fish访问

FLAC3D学习笔记(自己总结版)

FLAC3D学习笔记(自己总结版) 一、建立模型笔记1键盘上的键的功能:X,Y,Z分别控制所建立的模型围绕X,Y,Z三个轴做逆时针旋转、如果打开大小写锁定键,分边按XYZ时建立的模型会做顺时针旋转、+号的功能是加大模型旋转的角度，-号的功能是减小模型旋转的角度。2键盘上的上下左右四个键可以实现图形的向上向下,向左和向右的移动、3当在一个区域建立多个模型的时候,有时候为了镜像单个模块儿,需要预先编组,编组的方法如下:gen zon tunint p0 0 0 0 p1 010 0 p255 0 p3 0 05 p6 0105 dim sq2 sq21 sq2111 group(编组) z1(组名为:z1) &fill group tunnel4镜像整个模型的方法:首先判断即将镜像的模型在原模型的哪个方向(即三个轴的方向),然后在轴所在的方向上去两个点、其中一个是原点(origin(xyz))、5镜像部分模型的方法: gen zone reflect normal (x,y,z) origin (x,y,z) range group z1就可以了、6特别要注意的一点:交叉巷道的巷道充填和巷道的网格是两个组,所以在镜像的时候不要忘了给充填部分镜像、7 group radcyl range group cylint not (编组的时候不包括cylint)7 快捷键总结Ctrl+Z 局部缩放的功能。 Ctrl+R的功能是恢复到初始状态。8模型建立的方法:方法1:利用

默认值生成网格,各个默认值:(1)p0的默认值是(0,0,0),网格的每边的单元格数size默认值是10,网格每边的长度默认值为size的值、方法2:利用4个点的坐标来生成矩形网格、p0~p3 size方法3:利用edge来确定边长生成矩形网格、方法4:利用参数ratio来确定单元体几何边画边生成矩形网格、方法5:利用参数add(相对坐标)来生成矩形网格、方法6:利用八个点的左边来生成矩形网格、9、用户自定义模型的方法我来贡献一点自己的成果：; U" j/ ?7 L+ a0 uFLAC3D的二次开发环境提供了开放的用户接口，在软件安装文件中包含了软件自带所有本构模型的源代码，且给出了Mohr-Coulomb模型和应变软化模型的编译示例，因此可以方便地进行本构模型的修改与开发。为了方便起见，下面的说明以建立UserModel模型为例。c % b1 h、 \ K6 Q: l(1) 在模型头文件(usermodel、h)中进行新的本构模型派生类的声明，修改模型的ID(为避免与已有模型冲突，一般要求大于100)、名称和版本，修改派生类的私有成员，主要包括模型的基本参数及程序执行过程中主要的中间变量。( ? p6 u J5 Q3 y( a(2) 在程序C++文件(usermodel、cpp)中修改模型结构(UserModel::UserModel(bool bRegister): Constit L, t% V、D(4)

多心圆隧道用FLAC3D直接建模的方法

多心圆隧道直接用FLAC3D建模 丁其乐2013/6/18 多心圆隧道模型使用FLAC3D建模的难点在于FLAC3D中并没有以不规则曲线为边界的原始3D基本网格，隧道的断面图见图一 O4 图一 隧道内轮廓线是由多个圆弧做成的，所以内部区域模型只能通过3D基本网格拼接而成，这时我想到了使用cylinder（圆柱体形网格）拼接，但是这个圆弧的圆心并不是一个点，怎么办？这时我想到了通过fish来调整圆心位置，使得各圆弧的圆心都调整到O1位置，同时又要保证各圆弧的圆弧线不变，且内部网格点的相对位置不变。 下面我以一个简单的例子来说明我调整的过程： O1 O2 M1 M2 A B 图二 在图二中圆弧O1AB 的圆心为O1，圆心坐标为（X1,Z1）,现在通过坐标调整把圆弧O1AB 所占区域调整到O2AB（如上图所示），O2的坐标为（X2，Z2），同时O1A上面一点M1，调整之后在O2A上为M2。A和M1的坐标分别为A(Ax,Az)和M1(M1x,M1z)。 设：

|11| |1| O M k O A = （1） '|1| 1|1| M A k k O A = =- （2） M2的坐标为M2x,M2y 则： '21(21)M x M x k O x O x -=- （3） '21(21)M z M z k O z O z -=- （4） 由（3）、（4）就以求出来M2的坐标。 下面用一个简单的例子说明实现代码。 new title 调整圆弧区域到指定位置 ;定义参数 def para ;O1的平面坐标 O1x=-2.9 O1z=0 ;O2的平面坐标 O2x=0 O2z=0 r1=5.9 ；O1的半径 a1=0*degrad ;A 点与X 轴的夹角 a2=30*degrad ;B 点与X 轴的夹角 ;A B 点的坐标 Ax=r1*cos(a1)+O1x Az=r1*sin(a1)+O1z Bx=r1*cos(a2)+O1x Bz=r1*sin(a2)+O1z end para ;生产模型 gen zo cyl p0 (O1x,0,O1z) p1 (Ax,0,Az) p2(O1x,1,O1z) p3(Bx,0,Bz) size 4 1 6 group '调整' gen zo cyl p0 (O1x,2,O1z) p1 (Ax,2,Az) p2(O1x,3,O1z) p3(Bx,2,Bz) size 4 1 6 group '原始';为了对比而建 ;调整节点位置 def ding pg=gp_head loop while gp_group(pg,1) = '调整' xx=gp_xpos(pg)

Flac3D复杂三维模型建立_jx

基于Gocad-Surfer-犀牛-Ansys的Flac3d建模步骤 （金星整理） 准备软件： Autocad→Gocad→Surfer→Rhinoceros（犀牛）→Ansys→Ansys to Flac3D→Flac3D 建模思路： 以专业地质建模软件Gocad为基础，借助Ansys的强大建模能力，为Flac3d的建立复杂的三维地质模型，整个过程的文件目录尽量不要出现中文。 建模步骤： 1 AutoCAD 目的：将CAD文件的dwg格式转换成dxf格式文件。 (1) 清除无关图层 只留下含有等高线数据的图层。（建议将等高线另外复制出去到一个新的CAD文件中进行操作）。 (2) 赋予高程 一般地形图的等高线上都赋有高程信息，保证每条等高线都有。 (3) 保存留于GoCAD调用 另存文件为XXX.dxf文件（2007、2004、2000版本都可以），作为GoCAD软件的处理对象。 2 GoCAD 目的：将CAD中等高线信息转成txt格式文件。 (1) 打开GoCAD软件 打开软件，New project（自己命名），点保存，再确认后，弹出窗口选择Select all，再次确认。

(2) 导入文件 操作过程：File→Import Objects→Cultural Data→DXF，导入之后右下角有进度条，直到100％，可能会需要一点时间。 (3) 显示等高线图 操作过程：Objects→选中所需等高线图层。 (4) 保存文件 操作过程：File→Export Objects→Surface→Custom ASCII→Output File（选到文件夹，并键入名字）→增加X、Y、Z→OK，导出XXX.dat文件，直接不用处理，在Surfer中调用。

FLAC3D命令流(整理版)

1、怎样查看模型？ 答：plot grid 可以查看网格，plot grid num 可以查看节点号。 2、请问在圆柱体四周如何施加约束条件？ 答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2

FLAC3D双线隧道建模

new gen zone radcyl p0 9 0 0 p1 add 6 0 0 p2 add 0 200 0 p3 add 0 0 6 & size 3 20 6 3 dim 3 3 3 3 rat 1 1 1 1 group inner fill group out gen zone reflect norm 0 0 -1 origin 9 0 0 gen zone reflect norm -1 0 0 origin 9 0 0 gen zone reflect norm -1 0 0 origin 0 0 0 gen zone brick p0 -15 0 6 p1 add 30 0 0 p2 add 0 200 0 p3 add 0 0 12 & size 15 20 6 group j1 gen zone brick p0 -3 0 -6 p1 add 6 0 0 p2 add 0 200 0 p3 add 0 0 12 & size 3 20 6 group j2 gen zone brick p0 -15 0 -43 p1 add 30 0 0 p2 add 0 200 0 p3 add 0 0 37 & size 15 20 18 group j3 gen zone brick p0 15 0 -6 p1 add 37 0 0 p2 add 0 200 0 p3 add 0 0 24 & size 18 20 12 group j4 gen zone brick p0 15 0 -43 p1 add 37 0 0 p2 add 0 200 0 p3 add 0 0 37 & size 18 20 18 gen zone brick p0 -39 0 -6 p1 add 24 0 0 p2 add 0 200 0 p3 add 0 0 24 & size 12 20 12 gen zone brick p0 -39 0 -43 p1 add 24 0 0 p2 add 0 200 0 p3 add 0 0 37 & size 12 20 18 plot show add surface yellow add axes black quit save tun1.sav model mohr pro bulk=0.6e9 shear=0.36e9 friction=35 pro cohesion=1e5 tension=1e5 set gravity 0 0 -9.81 initial density=1000 apply szz -1.4e6 range z 17.9 18.1 fix x range x -39.1 38.9 fix x range x 51.9 52.1 fix y range y -0.1 0.1 fix y range y 199.9 200.1 fix z range z -43.1 -42.9 fix z range z 17.9 18.1 ini sxx -1.4e6 syy -1.4e6 szz -1.4e6 history unb h gp xdis 6 0 -6 h gp zdis 0 0 6 h gp xdis 6 25 -6