FLAC3D快速入门(手册翻译版).docx

快速入门(GETTINGSTARTED)版本：flac3d 3.0版FTD127 翻译：一米2023.06声明现在市面上关于FLAC3D软件的教材寥寥无几，在学习的过程中， 主要还是参考软件本身的使用手册，虽然读英文版手册有些吃力，但是它论述格外具体，我觉得是用户最好的教材。我在边看手册的时候边做了翻译，目前为止翻译完成了本局部的内容略去了局部内容和例子，还翻译了命令手册的前半局部内容，等翻译完成了，也会和网友共享，但是像本人这类英语水平一般的人做这样的翻译工作是比较辛苦的，我也不确定是否有毅力完成命令手册下半局部的内容。虽然这样的工作比较困难，但我觉得还是学到了不少东西，手册是最原始，最翔实的根底教材，看明白了手册，运用软件才会游刃有余。由于本人专业水平和英语力气的限制，存在问题是在所难免的， 有的地方甚至可能曲解了原意。考虑到时间因素，译文的措辞没有细 细斟酌，还请网友谅解。假设觉察译文中的错误，还请宽阔读者斧正。一米FLAC3D3.0 版本中文手册翻译：一米2 快速入门这一局部将向初次使用 flac3d 的用户介绍软件的根本使用方法。主要有以下内容：软件的安装与启动；用软件分析解决问题的步骤，在每一步的操作中，都有简洁例题来说明该步骤具体是如何操作的。假设你对软件比较生疏，但是现在很少用它来处理问题，那么这局部的内容尤其 2.7 节能很好的帮你回忆软件操作的要点。本局部 3.3 节全面具体的介绍了如何进展问题的求解。Flac3d 支持命令驱动和图形菜单驱动两种模式 *。在本手册中大局部的算例都承受了命令驱动模式。我们认为这种模式能给用户供给操作软件最清楚的思路。在 1.1 节中我们就已经提到了命令驱动模式使得 flac3d 在分析求解工程问题时成为了一个功能强大的“ 多面手” 。然而这种模式让用户，或者长时间未接触软件的老用户用起来有点不那么简洁。命令行必需用键盘输入，可以直接输入到软件的命令窗口，或者先保存为数据文件，再通过软件的相关命令进展读取。Flac3d 能识别超过 40 个主命令和 400 多个附属的关键词。本局部主要包括以下内容：1 在 2.1 节，手把手的教你们如何在自己的电脑上安装和启动 flac3d 软件。2 在 2.2 节，用一些简洁的教学案例帮组用户生疏一些常用的命令。3 在用户建立自己的模型并进展分析计算之前，有必要先了解 flac3d 的一些根本学问。在 2.3 节表达了 flac3d 的根本术语；在 2.4 节主要说明白有限差分网格的定义规章；而在 2.5 节阐述了输入命令的根本句法。4 在 2.6 节，阐述了 flac3d 的特点，比方创立、命名和使用对象，以便利用户进展问题的求解5 在 2.7 节，一步步的指导用户如何建模和分析问题，每一个步骤都分开论述，并供给简洁的例子帮助用户理解。6 2.8 节2.10 节分别论述了系统的符号商定、单位体系和精度限制7 2.11 节说明白软件中各种类型文件的创立和使用。8 2.12 节对图形菜单操作模式进展了简介。*：对于初级用户来说一般图形菜单驱动模式只进展图形输出或者文件操作。本章节的最终一局部将向用户呈现如何使用图形菜单驱动模式来操作软件。.2.2.1 安装启动程序2.1.1 系统要求安装运行 flac3d 需要的系统最低配置如下：处理器：时钟频率至少为 1GHZ，处理器的主频越高，那么 flac3d 的计算速度将越快。硬盘：安装软件至少需要 12MB 的硬盘空间。假设装载了在线的用户手册，那么还需要 16MB 的空间。留意默认状况下，安装软件时会自动装载用户手册。除此之外，还需要至少 100MB 的硬盘空间来存储分析计算时生成的各种文件。 内存启动软件至少需要 3MB 的内存。在建模过程中，软件所占用的内存，会 不断的发生变化见表 2.1WINDOW 操作系统还限定了软件建模时占用的内 存不能超过 2GB。显示器：推举 1024× 768 区分率，16 位彩色显示器。操作系统： FLAC3D 是 32 位操作系统的应用程序，所以基于 intel 技术的WINDOWS 98 及以上操作系统均支持软件的安装和使用。 输出设备：默认状况下，系统图形会输出到系统打印机上。也可以复制到剪贴板 上，或者保存为格式化的文件，这里所说的格式包括：加强型图元文件格式和位 图文件格式PCX/BMP/JPEG。用户可以使用 set plot 命令来指定输出的形式及 格式。2.1.2 软件的安装略2.1.3 组件软件的可执行文件为“ F3300.EXE” 。FLAC3D 是使用 VC+ 7.0 编写的。除了可执行程序外，还需要两套动态链接库DLL 文件，一套用来接入和存取各种各样的图形；另一套供给内置的各种本构模型。2.1.4 应用程序和图形处理设备在使用 FLAC3D 时，各种应用软件和图形处理设备会起到很大的关心作用。编辑器：任何以 ASCII 码为标准格式的文本编辑器都可以用来创立 FLAC3D 的数据文件。但是必需要留意一些“ 先进” 的文档编辑器如 WordPerfect, Word 等软件，这些编辑器会把格式说明信息编译成标准输出格式，这些说明信息并不能被 FLAC3D 识别，所以导入这类文档时会消灭错误。FLAC3D 输入的数据文件必需是标准 ASCII 码形式的文件。图形输出设备：FLAC3D 支持很多种类型的图形处理设备，默认状况下，生成的图形可以用“ Plot hardcopy” 命令来连接到系统默认的打印机以便输出。或者通过 FLAC3D 主窗口中 FILE 菜单栏下的 print-view 来设定“ Plot clipboard” 命令可以将显示的图形，存放到 WINDOWS 剪贴板上没有任何文件生成。该图形接着就可以以加强型图元文件格式被粘贴到其它兼容该格式的 WINDOWS 应用程序中去。“ Set plot metafile” 命令可以将图形以加强型图元格式存盘，以便作为计算的参考或日后插入到文档中去。通过命令：Set plot +关键词pcx, bitmap, bmp 或者 jpg可以存储为很多图像格式pcx,bmp,jpeg 等。输出的这些位图的区分率由命名行：Set plot <output type> size 来把握。固然也可以使用 Set plot avi 或者 Set plot dcx 以及 Set plot movie 命令将显示图形输出为视屏格式。无论是黑白的还是彩色的 postscript 打印机，都需要通过“ Set plot postscript” 命令来指定。打印图形将存储为文件，这样支持 postscript 格式的图形处理程序就可以读入并进展修改了。2.1.5 启动软件双击可执行文件“ F3300.EXE” 便启动了程序，接着会弹出一个 FLAC3D 的主窗口。在主窗口的最下面附带了一个命令窗口，我们可以把命令直接输入到命令窗口中来执行相关命令，命令窗口最初显示的提示符为：“ FLAC3D>” 。当软件启动后，它占用的系统内存是随着用户的操作而不断变化的比方说， 在建模过程中，系统所占用的内存会越来越多。我们可以在命令窗口中输入 print memory system 命令来查看现阶段程序已占用的内存及操作系统还可为软件供给的总内存。假设你在操作过程中觉察命令失效并不是命令错误，那么确定是系统可安排的内存太少了，软件所占用的内存过多。这个时候，最好退出并重启软件，以释放内存。表2.1 列出了一般建立摩尔库伦材料模型的单元数与软件占用的内存之间的大致对应关系。表 2.1FLAC3D 内存使用状况表单元数占用内存量MB20320303684031615033146035412.1.6 版本号说明略2.1.7 程序的初始化 刚翻开FLAC3D软件，它首先会在当前文件夹下查找“FLAC3D.INI”文件，如果没有找到，它就会到安装名目下查找。它的作用是存放用户设定的程序初始化模式的命令。以便每次翻开软件都载入用户的初始设置。假设“FLAC3D.INI”文件不存在，软件连续运行而不会提示出错信息，留意一点：一些存储在“FLAC3D.INI”里的命令，假设并不是设置初始化的命令，有可能导致错误的信息。2.1.8 运行FLAC3DFlac3d命令驱动模式包括两种方式：交互模式在命令窗口中输入命令行； 命令流模式将命令行保存在数据文件中，通过读入该文件执行相关命令。假设输入的命令存在错误，那么窗口中将会消灭错误提示。命令流文件一般通过文本编辑器创立和修改见2.14节，虽然命令流文件可以定义为任何文件名，但是最好设定其扩展名为“ .dat” ，以防止和flac3d其它类型的文件相混淆。要读入命令流文件可以使用以下命令：call file.dat其中，file.dat 指的是用户定义的命令流文件的文件名。一旦读入文件，你会觉察软件会将当前在文件中读入的命令行，显示在屏幕上。假设命令流文件保存在当前文件夹下*，那么在call命令后面只需输入完整的文件名即可，否则还应*笔者注：所谓的当前文件夹包括两种状况：1、没有读入任何数据时当前文件夹指的是软件应用程序所在的文家家。 2、假设已读入了数据，比方导入了模型信息文件“ .flac3d” 文件，这时当前文件夹指的就是用户之前读入文件所在的文件夹。在文件名前面加上文件的完整路径比方：c:我的文件夹 file.dat。除这种方法外，我们也可以菜单操作读入文件：依次点file-call按钮见2.12节 为便利起见，我们可以为应用程序创立快捷键，右键点击“ F3300.exe” 不放，并拖动到相应的创立快捷键的位置， 松手后会弹出一个对话框， 选中“ 在当前位置创立快捷方式” ，这样就生成了一个快捷方式。双击该快捷方式就可以启动软件。创立快捷方式的目的并不只在于便利翻开应用程序， 我们右键创立的快捷方式，选择“ 属性” ，接着在弹出的对话框中将“ 起始位置” 这个文本框中内容删除并点击左下角确实定按钮。这样当你双击该快捷方式启动应用程序时， 系统默认的“ 当前文件夹” 就是快捷方式所在的文件夹了。 我们可以将快捷方式和输入文件放在同一名目下，这样就便利了文件的输入。2.1.9 装载测试文件略FLAC3D 3.0 版本中文手册翻译：一米2.2 一个简洁的计算教程常用命令的使用这一局部主要是为那些刚接触FLAC3D，跃跃欲试的用户预备的。在这一局部，将通过一个简洁的例子来帮助用户学习一些求解问题的根本学问。例题的主要问题描述如下：在一块土体中一次性开挖一个 2m× 4m× 4m的沟渠，并对沟渠四周土体的变形作监测和分析。为了给用户供给便利，在安装名目中“TutorialBeginner”文件夹下的“ TUT.DAT” 数据文件里包含了本例题使用的全部命令。我们承受交互式的方法在FLAC3D中运行该算例。即：用键盘输入命令行， 在命令行输入完成后按回车键，然后直接查看输入命令后的结果首先让我们翻开FLAC3D软件，翻开后，便可以看到FLAC3D的主窗口，在主窗口下方，是命名窗口，命令窗口的初始提示符为“ FLAC3D>” 。 建立初始的有限差分网格的命令为GENERATE*：GEN ZONE BRICK SIZE 6 8 8这个命令会在软件中建立一个X方向有6格，Y方向有8格，Z方向有8格的三维长方体网格。在我们建立的这个模型中，Z方向为竖直方向。我们可以用PLOT 命令进入绘图模式，在该模式中可以观看所建立的网格。一旦输入PLOT命令，那么接下来所要输入的命令都必需为 PLOT子命令，绘图窗口默认定义为一个名为“ BASE/0” 的窗口，我们可以看到这时提示符变成了“ Plot Base/0>” 。为了显示彩色的网格体和模型的坐标系，就需要输入如下的命令行： 例2.2 创立和使用视窗create Trench创立了一个叫做“ TRENCH” 的视图窗口 add surface yellow外表显示黄色 add axes black显示黑色的坐标轴 show在屏幕上显示 创立的视窗自动的设置为当前活动视窗。按键盘上的x 、y、 z键小写状态可使得当前视图围着X Y Z轴旋转大*：该命令的具体信息请参见命令手册局部的第一章节。留意命令可以缩写，变量与其对应的参数值之间可以试用“ ” 符号来分隔见2.5节。.7.写状态时旋转方向相反使用m(M)键则可以对视图进展放大缩小,方向键则可以把握视图进展平移。固然反方向旋转也可以使用shift+把握键实现其原理就是大小写之分对视图的其它一些操作和把握，可以参考命令手册局部的表1.2。留意：我们也可以通过PLOT 命令直接指定放大的倍数和视角。每一个视图窗口的左边都会有该视图的具体信息：CENTER(中心);ROTATION旋转角度； DIST(距离)；MAG放大倍数等。我们现在通过交互模式建立我们想要的视图，而这些视图的数据信息其实还可以转变成数据文件，这样在其它视图中想要得到类似的视图效果，就可以通过命令输入一样的数据。实现过程见例题2.3 我们现在在例题2.2的根底上再建立一个视图窗口通过输入以下命令来创立一个视图。例2.3 指定视窗显示数据create Trench2 add surface yellow add axes blackset rotation x y z set center x y z set dist dset magnification mshow例题2.3中多输入了加了黑体字的局部(红色局部的数字等同于例题2.2中经用户调整后显示在屏幕左侧中对应的数据)，该局部等同于用键盘把握视图， 2.3和2.2的两个视图应当是一模一样的。 还是来看例子2.2：在Trench视窗中我们建立了一个6× 8× 8m坐标原点在0，0，0的网格。网格的尺寸，原点，以及网格的密度都是可以通过GEN+ 关键词的命令来更改的。这个问题将在稍后再作介绍。 现在我们要给整个模型空间定义本构模型和具体材料参数。在本例中我们定义为摩尔库伦弹塑性本构模型。反回到“ Flac3D>” 命令模式如在PLOT模式中，按回车键即可返回并输入：model mohr这个命令就定义了该网格体的本构模型为摩尔库伦模型。网格中的每块区域都可能需要定义成不同的本构模型和赐予不同的材料参数。但是假设MODEL 命令后没有特别指定命令作用的区域，则系统默认为作用的范围为整个模型。为了提高塑性解答的精度，每一个单元都会自动的被细分为两组四周体的子网格参见理论与背景局部的第一章。 想知道计算这个问题需要占用多少的内存，可以输入以下命令：print mem sav mem sys输入命令后软件便会供给一张在现阶段内存的统计列表。我们这个例子在现阶段或许需要占用0.6MB的内存来存储数据，大约需要4.2MB的内存来进展模型的计算。接着在命令窗口中输入材料属性赋值命令：prop bulk = 1e8 shear = 0.3e8 fric = 35 prop coh = 1e10 tens = 1e10以上的命令定义了材料的体积模量PA、剪切模量、内摩擦角、粘聚力和抗拉强度。用户会觉察粘聚力和抗拉强度值取得很大，这样取值的目的是为了模拟得到模型在自重作用下的初始应力状态 *。取值很大是防止了模型在初始加载分析时就到达塑性极限。一旦你拥有了很多的模拟阅历，你就会知道这样做的缘由是很明显的。在该样例中，土体仅受自重，在命令行中输入如下命令来加载重力：set grav 0, 0, -9.81重力加速度,单位：m/s2 ini dens = 1000密度，单位：Kg/m3 命令行中：SET命令设置了重力加速度，INI定义了初始网格的质量密度只有考虑重力的状况下才需要定义密度，在静力问题中，假设我们无视重力，那么就无需定义密度，这点与FLAC不同。接着我们用如下命令来定义模型的边界条件：fix x range x -0.1 0.1*：有很多方法都可以模拟到达模型的初始平衡状态，比方可以先给模型赋予弹性本构模型来计算初始应力状态，再重定义为摩尔库伦模型来分析计算模型的开挖、加载及各种条件的转变过程；还可以直接赐予摩尔库伦模型，然后承受solve elastic命令来分阶段求解。fix x range x 5.9 6.1fix y range y -0.1 0.1fix y range y 7.9 8.1fix z range z -0.1 0.1以上的命令固定了模型的五个面连杆支承，边界一旦被“ FIX” 固定后，在被固定的方向就不会发生位移和产生速度。FIX命令在例题中起到了如下作用：1、固定了边界面X=0,X=6上的全部节点X方向的位移，由于这两个边界面分别在range关键词所指的范围内；命令的前两行2、固定了边界面Y=0,Y=8上的全部节点Y方向的位移；第三、第四行的命令3、命令的最终一行固定了底部边界面Z=0Z方向的位移。 假设我们想在迭代计算过程中记录一些有用的变量或参数的变化，从而来判断分析是否已经到达平衡状态或者模型已经发生垮塌，就需要使用到HISTORYHIST命令。 在本算例中，我们接着在命令窗口输入如下命令：hist n= 5 hist unbalhist gp zdisp 4,4,8命令的第一行：指定了参数记录的步频： N(Nstep)=5指定了每迭代计算 5次记录一次相关的值假设不指定，则系统默认值为 10。命令的其次行指定记录最大不平衡力the maximum unbalanced force命令的第三行指定记录坐标值为4，4，8的节点Z方向的位移。在计算时记录最大不平衡力是一个很好的习惯，假设最大不平衡力接近一个很小的值，并且位移记录值不再发生变化，那就说明计算已经到达了平衡状态。现在已经预备好了求解模型初始状态的必要条件。由于FLAC3D计算的结果为显式的动态解，我们通过设定计算所需的时间步来把握计算的进程，模型的动能会渐渐衰减为零，这样就得到了我们所要的静态解。为了使单元体获得重力场， 我们需要设定计算的时间步来求解模型在自重的作用下初始平衡状态。软件中FLAC3D3.0 版本中文手册翻译：一米SOLVE命令就是按确定精度自动求解平衡状态的命令。了解了这些学问，我们接着回到例题中来，在命令窗口中输入如下命令：set mech force=50 solve这样，软件便开头了问题的求解，在求解过程中，会在程序窗口中显示计算的时间步数以及对应的最大不平衡力的值，当最大不平衡力小于所设定的上限值本例题中我们用SET命令设定了最大不平衡力的上限值为 50N,即：计算的最大不平衡力一旦小于此值，我们便人为模型到达了平衡状态。求解过程终止*。由于我们之前没有关闭绘图窗口，你会觉察显示在窗口中图像不断的被更。假设我们先关闭了绘图窗口，那么求解所需的时间会更短。在我们的模型中，计算停顿在351步*，对于一般主频为1GHZ的奔腾电脑， 完成这样的计算仅需要短短的几秒钟时间。我们可以通过查看最大不平衡力以及所记录节点Z方向的位移采样图来分析，模型是否到达了平衡状态。在命令窗口中输入：plot hist 1输入完命令后便可以在窗口中显示HIST 1(最大不平衡力)的采样记录图。按回车键，再输入如下命令：hist 2 可以查看所记录节点Z方向的位移采样图了。 最大不平衡力采样图图2.3中显示最大不平衡力接近于零，位移图图2.4中可以看出位移已经趋于一个固定值。这两个采样结果均说明系统已经到达了平衡状态。 上面两个采样图都是在FLAC3D默认视图窗口BASE/0中显示的。任何图形假设没有定义的视图窗口，都会直接在默认视图中加以显示，这并不会对已建立的自定义视图窗口产生影响。默认视图窗口就像一个“ 抓拍” 器，能马上*：计算精度的有关内容将在2.7.4节中具体争论。*：由于用户处理器的不同，实际计算步数和给出的计算步数会有微小差异。.10.FLAC3D 3.0 版本中文手册翻译：一米让你观看到你所需要的视图。图2.3最大不平衡力记录图图2.4节点4，4，8z方向位移记录图.11.FLAC3D 3.0 版本中文手册翻译：一米留意：输入的采样记录，系统按输入挨次从“ 1” 开头始终往下编号，这就是为什么输入hist 1命令显示的是最大不平衡力的采样记录图，输入hist 2命令可显示Z方向的位图采样图的缘由了。当采样参数较多，自己又不记得输入的挨次时，可以在“ FLAC3D>” 提示符状态下输入：print hist命令来显示全部采样参数对应的编号称为id号。接着我们通过下面几行PLOT命令来建立一个比较简洁的视图，以检查重力场是否已经施加到模型中：plot create GravVplot set plane dip=90 dd=0 origin=3,4,0 plot set rot 15 0 20plot set center 2.5 4.2 4.0 plot add bound behind plot add bcont szz plane plot add axesplot show留意以上命令均在“ FLAC3D>”提示符状态下输入，假设先输入plot命令， 进入“ plot” 相关的提示符状态，那么命令中全部的plot均可删除。执行以上命 令后，便会消灭一个名叫GravV的视窗，并且系统指定它为当前显示窗口。我们在命令中已经定义了当前视图的剖面：倾角为90度dip=90，倾向为0度dd=0， 通过点3，4，0。在命令中我们还定了视角plot set rot 15 0 20 ；plot set center2.5 4.2 4.0并且使用boundboundary命令来增加边界限框，线框只在平面的后面显示由关键字behind定义，注可简写为BE。plot add bcont szz plane 命令的作用是在当前平面中绘出垂直应力分布图，图形如图2.5所示。.12.图2.5 初始平衡状态z向应力剖面图现在，最好保存以上得出的初始状态，便于在将来任何时候返回到该状态， 重设定相关参数来做相应的争论。按回车键返回到“ FLAC3D>” 提示符状态， 输入下面的命令：save trench.sav 这样一个名叫trench的SAV类型文件就在当前文件夹下生成了。到此，我们在模型中已经建了两个视图窗口。我们可以用plot print命令来列出全部的视图名称，当前视图前会有“ =>” 这样的标志以示区分。也可以用plot print view命令来显示当前视窗的具体信息。在以上命令后加上视窗名或者其对应的id号，就可以显示相应视窗的具体信息，包括：全部的视图设置；记录了哪些参数。假设想从GravV视图返回到早先创立的Trench视窗，那么输入如下命令：plot current Trench。实现以上目的还有另一种方法，输入：plot show Trench。只要你知道你想设置为当前视图的视图名，你就可以通过以上两种方法实现。固然也可以承受菜单操作：当鼠标输入光标指示在命令窗口时，在主窗口中有一个 PLOT 菜单，在该菜单中排列了全部用户定义的视图，鼠标点到你想看到的视图，点选SHOW子选项，你就把该视图设置成了当前视图。我们再进展下一步工作，在土体中开挖沟渠，首先输入命令：prop coh=1e3 tens=1e3这个命令重定义了整个模型材料的粘聚力和抗拉强度为1000Pa，这样的取值其实已经能够保证在初始状态中不消灭错误了未开挖前，但是我们应当总是要通过一些计算步数的计算来检查在初始状态可能消灭的失稳垮塌。为了模拟开挖，我们只需把要开挖局部的网格体的本构模型设置成空模型 null即可， 命令如下：model null range x=2,4 y=2,6 z=5,10开挖局部即模型材料为NULL是通过限定x，y，z方向区域的方法来选取的。由于粘聚力比较小，而且沟渠壁没有支护，因此开挖后应当会发生土体垮塌的现象。我们想要模拟的是现实过程，那么通过上面的简洁分析，设定计算为大变形是合情合理的，设定的命令如下：set large我们需要看到的仅是开挖这个过程引起的位移变化。而不是从加载重力到开挖整个过程的位移变化，所以，系统中全部网格节点的位移应当全部清零 *命令如下：ini xdis=0 ydis=0 zdis=0我们有意取了较小的粘聚力值，以保证模型发生破坏。这样一来，我们就不能用SOVE最大不平衡力上限值推断到达平衡的条件这样的命令方式来求解问题了，由于模拟的计算将永久不收敛而不能到达平衡状态。我们可以通过限定每次计算的时间步数来把握模拟过程，从而了解在坍塌发生时的有关状况，这里才真正表达了显式求解方法的精华。计算过程永久不会渐渐趋于平衡，这也是有别于一般工程师所生疏的隐式求解线性方程组的过程。在FLAC3D中要求得这类解，就要用到STEP命令，在本例题中输入：step 2023，这样FLAC3D就会执行2023时间步的计算。对于主频为1GHZ的奔腾电脑，这个过程都花不了一分钟的时间然后我们要通过查看相关图形，来检查到现阶段为止的计算结果。现在我们先来查看位移的等值线图，输入如下命令：*：这个步骤不会影响到模型的计算，只是便利于用户查看需要的结果，因为在迭代计算过程中，没有用到位移。FLAC3D3.0 版本中文手册翻译：一米plotcreate DispContcopy GravV DispCont settingsadd cont disp plane behind shade on add axesshow在命令中我们拷贝了之前定义的GravV视图的设置为当前定义视图的设置具体来说就是平面的定位设置。在消灭的位移等值线图中图 2.6你会觉察，在开挖局部的四周一些网格开头发生变形。在图中你还可以看出因开挖，地面发生沉降的区域*。图2.6 开挖计算2023步后的位移等值线图 本例题临时就到此完毕，在接下来的内容中我们将介绍FLAC3D的其它一些学问，来解决本例题模型一些更简洁的问题，包括怎样对开挖的沟渠进展加固。我们建议初学者承受边看教学手册边操作的模式来学习接下来的学问。 假设你对FLAC3D的命令已经比较生疏了，那可以转到3.3章节来了解FLAC3D在求解问题时还需了解的其它细节信息。*：考虑到加快例题的计算时间，我们取的开挖土体边界很小，在实际状况中，开挖可能是在很大的一块地面上进展的。那么自然的，我们认为定义的边界会对计算结果产生影响。.15.2.3 根本术语FLAC3D使用的术语与一般的有限差分程序或者有限元程序是全都的， 为了使用户能够清楚的了解这些术语，现图文并茂加以说明，图2.7中标示了flac3d中的根本术语。图2.7根本术语说明图FLAC3D MODEL用户在FLAC3D软件里建立的一个物理模型， 用户通过输入命令的方式来定义模型所处的各种状态，已完成对真实问题的数值模拟过程。ZONE单元体，也叫ELEMENT,是分析现象如：应力应变关系时发生变 化的最小空间区域，空间单元可以有不同的外形如：六面体、楔体、锥体、四面体等，模型的网格是由这些最根本外形的单元体组成的。每一个单元都自动再细分为两套Sub-zone*，这个和用户建模无关，一般不需要去深入了解。GRIDPOINT节点单元体的角点，也叫node(或 nodal point)。任一个单元体都是由它上面的每个节点的三维坐标值而准确定位的。FINITE DIFFERENCE GRID有限差分网格，又叫MESH，是由各种外形的单元体组合而成的分析问题的区域，它也可以确定模型中变量的存储位置，全部的向量力、速度、位移等都存储在节点位置，全部标量和张量应力、材料参数都存储在单元体的中心位置。MODEL BOUNDARY即网格体外部边界和内部边界网格体内的孔洞 BOUNDARY CONDITION边界条件，模型边界的约束或者把握条件如位移的约束，在边界上施加的力，在渗流模型中的不透水边界，在热交换问题中的绝热边界等INITIAL CONDITIONS初始条件，在加载和扰动前全部变量的初始状态。CONSTITUTIVE MODEL本构模型，也叫MATERIAL MODEL，用理论公式表示模型材料的变形和受力特性。软件中附带有很多本构模型以用来模拟各种岩土材料的各种问题，每个单元都可以单独的定义本构模型和相关模型的材料参数。NULL ZONE空单元，用来定义计算时不需赋材料的网格区域，通常用来模拟开挖或者临时还未填筑的材料。SUB-GRID子网格，网格可以由子网格组成，子网格通常用来创立那些模型中不同外形的子区域比方：在建立大坝模型时，可独立建立坝体和地基的子网格，子网格是各自生成后再通过移动和绑定操作来组装成整个模型的。ATTACHED FACES绑定面，是指子网格被连接和绑定在一起的网格面， 两相互连接的子网格绑定面确定要共面和接触，但是不要求节点相互匹配，不同单元体划分密度的两个子网格也可以相互绑定在一起。INTERFACE接触面，是两个子网格间的一种连接面，通过这样的连接面连*笔者注：网上关于它的解析，subzone和flac3d的Mixed Discretization算法有关，比方说一个brick，计算时又自动划分为 2套overlay ，每套overlay中各有5个tetra四周体即subzone，据说这样可以更好地模拟材料的塑性变形。接的子网格在计算中可以被分开如产生相互滑动、开裂接触面通常用来模拟物质的不连续性，比方断层、实际接触面、两种不同材料间的分界面。 RANGE在FLAC3D中它用来定义所指向的区域，它限定了命令作用的区域。 GROUP组，组是FLAC3D中用来定义单元体集合的术语，组只能唯一命名， 不同的组要承受不同的组名，组用来限定命令有效性的空间域，命令后面加上了 组名， 就表示，该命令只对该组内的成员有效。ID NUMBER编号id号，模型中的各个元素都有一个数字编号，同类元素之间就是以这个编号加以区分和识别的。下面这些元素都具有ID标识号：接触面interface、节点gp、单元zone、体(volume)、参考点reference points、历史记录项history、表table、绘图项plot items以及构造单元包括： 梁、索、桩、壳、衬砌、土工栅格等可以用PRINT命令来得到元素的ID号。用户可以自行定义当前输入的接触面、参考点、历史记录项、表、绘图项的id号。每个构造单元中还有CID编号，CID号用来标识该构造体中各个元素，构造体中的节点，单元，以及之间联系部件等都有各自的CID号。STRUCTURAL ELEMENT构造单元， FLAC3D中有两种类型的构造单元， 两节点的直线段用来模拟梁、索、桩，三节点的三角形用来模拟壳、衬砌、土工栅格。构造单元通常用来模拟土或岩石与支承构造的相互作用。 这些单元是典型的非线性材料。每一种构造单元，都由以下三局部组成：节点、构造类型SELs、节点间的联系构造。STEP计算时间步，也叫time step或者cycle由于FLAC3D为显式码解法，求解一个问题时需要给定计算求解步数。一般简洁典型的问题求解只需要2023到4000 步，然而对于简洁的大型模型就可能需要上万步的求解才能到达稳定状态。假设是做动力分析，则必需在STEP命令后加上确定的求解步数。STATIC SOLUTION静态解，假设求解的相关变量的增量随着时间步的变化只发生微小的变化，而这个变化对于变量本身来说可以无视不计，那么我们就人为模型已经到达了相对平衡状态，计算的结果便是所求的静态解。在静态解模式中，要么模型到达了平衡状态，要么模型外荷载作用局部区域或者全部区域下到达了相对稳定的塑性流淌状态在发生滑坡等失稳状态下。静态解模式是FLAC3D 3.0 版本中文手册翻译：一米FLAC3D的默认求解模式*。静态解可以和地下不稳定渗流、热交换问题进展耦合。UNBALANCED FORCE不平衡力，也称为unbalanced或者out-of-balance force。是数值计算迭代过程中产生的系统内外力之差。在静态求解模式中用它来推断是否到达平衡状态或者模型刚开头发生塑性变形的状态假设模型网格的每个节点力都为0，那么模型应当到达了确定的平衡状态。在FLAC3D中用STEP 或SOLVE命令执行计算的过程中，系统会自动记录最大的节点力，并在屏幕上显示。最大的节点力也叫最大不平衡力，在数值分析中，最大不平衡力不行能为零，但是只要最大不平衡力与作用在体系上的外力相比小到可以无视不计时，我们便认为体系到达了平衡状态。假设在计算中，觉察最大不平衡力趋于一个非零值，那么就可以确定模型中发生了坍塌或者塑性变形。DYNAMIC SOLUTION动态解模式，用来求解动力学问题，一般适用于外力为高频率震荡且持续时间较短的状况比方：地震或者爆破。LARGE STRAIN/SMALL STRAIN大变形/小变形，在软件中，系统默认为小变形模式，即：系统中的节点坐标始终不发生变化， 即使计算所得的应变比较大与一般的单元尺寸相比。在大变形模式中，可能发生几何非线性的变化。2.4 有限差分网格有限差分网格掩盖了整个所要分析的物理区域，最小的网格可以只包含一个单元，但是大局部模型的网格都是由成千上万个单元组成的。有限元网格由全局的坐标系定义，全部的节点和单元体的重心都由x,y,z 坐标值定义。节点和单元也可通过ID号加以识别。ID号对应的节点全局坐标， 可以通过如下命令得到：print gp positionID号对应的单元重心全局坐标可以通过如下命令得到：*:在一些有限元的文献中存在这样的错误观念：有限元解法能够准确符合支配方程组，动态解法与有限元解法相比不行能到达真实的平衡状态。事实上，两种解法都只是近似的满足支配方程组，但是剩余误差可以被把握在尽可能小的范围内。在FLAC3D中这个误差被量化为节点上的不平衡力与作用在节点上全部