FLAC常用命令指南.docx

常用命令指南1 命令指南这一节包括flac使用的所有命令的详细信息。命令被描述为两个主要部分：首先，在节1.2中，这里有命令的概要汇总，它是以相关的模型函数来组织的。准备一个输入dat文件时建议先看一看这个命令汇总。其次，在节1.3中，按字母顺序对所有的命令进行详细的描述。节1.1中描述了一些共同的输入惯例和约定。在command and fish reference summary书中节1中同样提供了按字母顺序排例的命令汇总。1.1 共同的惯例和约定1.1.1 语法flac可以用“交互”的方式进行操作（例如，通过键盘输入命令）或“文件控制”方式（例如，数据保存和读取在一个数据文件中，而这个文件保存在磁盘或硬盘中）。无论哪种方式，解决一个问题的命令是相同的，数据输入的具体方法要看用户的偏爱。所有的命令都是单词导向的，它包含一个主命令词，根据需要后面跟着一个或多个关键词和数值。有些命令（例如plot）接受一个“开关”，它后面的关键字可以修改命令的动作。每一个命令都有如下的格式：command keyword value . 命令 关键词 值 所有的命令在输入行中逐字被输入。你可以注意到只有少数几个字母是黑体，这表示要这个程序认识这个命令，最少要输入黑体字母。同样地，小写字母显示的关键词也是逐字输入，只须要输入关键词的黑体字母长度即可。如果用户需要，可以将命令和关键词写完整。缺省情况下，单词是对大小写不敏感的-你可以使用大写字母，也可以使用小写字母。许多关键词后有一系列的数字（值），这是关键词所要求的。以黑体斜体表示的单词代表数值。当单词以i,j,m或n开始时，后面应接整型数；否则，后面应接浮点型数（或小数）。有时浮点型数中的小数点可以被省略，但是小数点不能出现在整型数中。命令、关键词和数值可以通过任何多个空格分开或者通过下列的分界符分开：() , =你将可以看到在一些输入参数后有一些附加的符号。它们是： 表示可选择的参数（尖括号不要输入） 表示可以给定这种参数的任意的数值。输入行中分号（；）后面的任何内容可以当作是注释，可以忽略。当运行于批处理模式（例如文件控制）时，在输入文件中这种注释是有用的，因为注释会在输出中出现。一个单独的输入行，包括注释，最多只能包含80个字符。一行的末尾出现&符号表示接下来的行是关键词或数值的继续。单个命令的最大长度（包括续行），是2000字符。在一个命令中最多可以有400个输入参数。1.1.2 交互式输入当flac在运算、从数据文件中读取或绘图时，flac可以支持用户提前输入的80个字符的缓冲。同样存在几个行编辑特征，它可以在交互地输入数据时使用。这些特征汇总在表1.1中。表1.1 交互输入编辑键键效果任何字符键在输入行中插入字符在输入行中向左移动光标在输入行中向右移动光标ctrl 光标跳至下一个输入参数的左边ctrl 光标跳至下一个输入参数的右边backspace删除光标左边的字符delete删除光标处的字符end移动光标到输入行的末尾esc取消行输入f3用上次键盘输入行替换输入行home移动光标到输入行的开始1.1.3 可选的关键词有些命令后面可以跟可选的关键词和数值，它们可以修改或限制命令执行的范围（网格点或区带）。特殊的命令和关键词汇总在表1.2中。表1.2 修改命令执行范围的关键词用法关键词用法命令i=i1,i2 j=j1,j2 (网格点/区带范围)apply，fix, free, generate, initial, interior, mark, model, print, property, unmarkfrom i=i1,i2 j=j1,j2 (边界路径范围)apply, attach, interface, initial, printvar vx, vy (递变)apply, initial, interior, propertymark (固定网格点)fix， free， initialregion i,j (区带范围)initial, model, print, propertyrdev s (随机偏差)property规则：1. 有两种方法来定义变量施加的范围：网格点/区带范围和边界路径范围。2. 当使用网格点/区带范围和边界路径范围时，它必须出现在输入行的最后。两种范围类型不能同时出现一个输入行中。3. 递变关键词后必须紧跟着数值，用它来指定递变方法。4. 递变的值，vx和vy，将赋值给网格点范围，既使变量是一个区带变量。5. mark关键词必须跟随指定网格点条件。6. region关键词必须跟随指定区带的条件。下面将对关键词用法作出更多的解释，并给出例子来说明。1.1.3.1 网格点/区带范围变量施加于上的网格点/区带范围通过关键词短语来定义：i=i1,i2 j=j1,j2这里，i1,j1是变量施加于上的边界范围的网格点/区带的起点，而i2,j2是网格点/区带的终点。范围对应于变量类型（例对于网格点变量是网格点范围，对于区带变量是区带范围）。范围必须在输入行的最后，但是i或j都可以列在前。如果i1=i2或j1=j2，只须给出一个数值即可。如果i1,i2或j1,j2就是网络的整个i或j范围，则在命令中可以省略关键词和相关的数值。下面这个例子说明了网格点/区带范围的使用。为了固定i=1到10 j=1的所有网格点在x方向的位移（或速度），下面的命令产生同样的效果。fix x i=1,10 j=1,1或 fix x j=1 i=1,10 或 f x i=1,10 j=11.1.3.2 边界路径范围如果边界是水平的或垂直的，这多发生在网格的外边界上，网格点/区带范围是一个简单的操作。如果边界是曲线或不规则的，则应用变量就会十分困难，例一个洞或挖空。因此，在两个特定的网格点之间提供一个面向边界路径的代码，边界路径范围的关键词短语为： from i1,j1 to i2,j2此外，在边界路径中还有3个可选参数long, short或both。这些参数混合在from i1,j1 to i2,j2范围描述中：long指定两点之间的最长路径，short指定路径的最短路径，both累加long和short路径。缺省路径是两点间的最短距离。注意：如果i1,j1和i2,j2是相同的，则整个封闭边界都将被当成路径范围。下面这个例子说明了边界路径范围的使用。沿着从点(6,3)到点(4.9)的边界施加一个机械压力，使用如下命令：apply pres=10 from 6,3 to 4,9此外，在两个网格点之间沿着最长路径施加一个不同的压力，使用：apply pres=30 long from 6,3 to 4,9注意：如果给定的点不在边界上，或如果两点之间没有路径，则flac会产生一个错误。如果在initial命令中使用边界路径范围，它将只被应用到网格点值。1.1.3.3 递变为了定义一个递变的应用，另一个可选关键词短语也可以用在一定的关键词后（表1.2）（例如边界上一个压力或应力递变）。这里，短语是： var vx vy它跟着指定递变的值。这个值指定在网格点/区带范围或指定边界路径上，将vx和vy值增加到终点上，vx和vy是指x轴的变化量和y轴的变化量。注意：递变区域的终点的值被赋值到更低或更高的网格点范围，甚至这个变化是一个0变化。1.1.3.4 标记网格点可选关键词mark让用户来限制指定条件到用mark命令来“标记”的网格点。例如，如果在指定的边界上指定x方向的速度，则应用如下命令：mark i=11 j=7,11mark i=1 j=7,11initial xv=-10 mark在这个例子中，mark命令定义了两个区域，而x速度值初始化中使用了mark关键词，当generate命令使用了arc, circle, line 或table关键词时，网格点被自动标记。1.1.3.5 区带范围表1.2中可选关键词region i,j让用户来限制命令的范围，这里必须有用一串标记的网格点封闭起来的封闭区带（例如，网格点用mark命令来标记）。给出标记区域中的任何一个区带i,j即可。例如，例子1.1中园形通道被挖空所使用的命令。例子1.1 定义一个园形通道作为一个区带范围grid 10,10m egen circ 5 5 3print markm null region 6,6命令gen circ 5 5 3 自动标记了园边界上的网格点。命令m null region 6,6将所有园形区域内的区带设为空物质。区带(6,6)是区域中的一个区带。用户可以体验定义一个区带宽的区域的难度。为了避免这种问题，总是在一系列的命令的最后才使用region命令来使行动产生的在一个元素宽的区域上。例如，例子1.2的数据文件产生一个元素宽的区域，在网格的最底部使用mohr-coulomb物质，在上半部使用elastic物质。例子1.2 一个元素宽的区域模型grid 10 10mod elasgen line 0 5 10 5gen 0,5 0,10 5,10 4,5 i=1,5 j=6,11gen 5,5 6,10 10,10 10,5 i=6,11 j=6,11mark i=5 j=6,11mark i=6 j=6,11mod mohr reg 3,3mod elas reg 6,8如果命令mod elas reg 6,8没有给出，model mohr命令的结果将会产生一个元素宽区域的“漏洞”。1.1.3.6 随机偏差当使用property命令时，对于可选关键词rdev s，紧跟着一个属性值, v, 使得属性值成为标准偏差s的一个随机偏差，平均值是v。在属性值之间没有空间的相关性。如果s是较大的，要确保属保属性值不是是负值。作为一个例子，下面的命令将给出40度的fric(摩擦角)存在一个5%的标准偏差： prop friction 40 rdev 2 1.1.3.7 一般注解捆绑的语言，fish语言，可以用来描述更复杂的变量和条件。例如，它可以提供非线性变化的属性或应力，随机属性变化，非均衡的边界条件（参见fish volume的第3节）。flac输入行中任何数字都可以被fish变量或函数名代替，当然事前得有一个define end语句定义该变量或函数名。任何输入行中的文件名都可以被名称代替，这里名称是一个fish的字符串变量或函数；名称中的字符顺序与文件名一致。更详细的内容请参见fish volume的第2.3节。1.1.4 在线帮助对于每一个flac命令，可以用help关键来得到该命令的所有可用的关键词列表。例如，为了了解print命令的所有关键词，输入 print help？可以用来代替help关键词。注意：有些通过help列出的关键词没有在手册中描述。1.2 按作用分类的命令下面是按照命令的功能而分类的命令输入顺序。通常地，命令可以以任何逻辑顺序给出；但是，有些命令必须在别一些命令之前出现。在这一节中可以区别它们。这里只给出用户最常用的最基本的命令。所有关键词的更详细的描述参见1.2.3节。1.2.1 特定的程序控制有一些命令可以让用户不用离开flac就可以启动一个新的分析或重新开始一个以前的模型模拟和从上一次分析结果接着计算。下面的命令提供程序控制。call 从用户准备的批输入数据文件中读到flac中，并执行这些命令。这就称为批处理模式。continue 在由于一个暂停或错误中止后，继续从批处理文件中读入。new 不用关闭flac而重新开始一个新问题的分析。pause在读一个批处理文件时暂停。quit停止flac的执行，返回操作系统控制状态。restore从已经存在的保存以前执行过的状态中恢复当时的状态。return让程序从批处理模式返回交互模式（或者在多个嵌套调用的情况下调用文件）。save保存当前的分析状态到一个文件中。set 设置基本的参数。stop 停止flac的执行，返回操作系统控制状态。最好给保存文件一个不同的扩展名（例如“.sav”），以区别输入文件（例如“.dat”），这样就可以避免在一个保存状态被restore或一个输入文件被call时产生混乱。1.2.2 指定特殊的计算模式flac在标准模式下执行机械（平面应力）计算。也可以用其它的计算模式，通过config命令来指定。如果config命令后面有下面的参数，则它必须在grid命令之前给定。下面的关键词可以给出一个或多个。axisymmetry轴对称的几何模型creep蠕动分析（只在蠕动模型中可用）dynamic完整的动力学分析（只在动力学模型中可用）extra捆绑语言（例fish）中使用额外的网格变量gwflow地下水流p_stress平面应力thermal热分析（只在热模型中可用）除了axisymmetry和planestress不能组合外，参数可以自由组合。1.2.3 问题几何学输入在计算网格时，必须通过grid命令给出网格的行和列的数目。下面的命令可以用来必变网格的形状：attach允许次网格沿着两个指定的边界结合generate用来改变网格的形状。网格产生变形有两个规则。首先，模型的区域（或次区域）通过指定一个坐标范围（或范围）可以再分区。其次，网格可以通过下面的关键词：arc, circle和line来定义一个曲面。或者可以用table命令和table关键词来产生用户所需的形状。initial 允许一个或一组网格点被移动。在模型开始计算以前，所有问题的几何形状命令都必须给出。这些包括，在后面的分析中要挖空的区域形状的产生。1.2.4 在模型中设定区域界线使用mark命令和在几个命令中使用mark关键词可以让用户定义一个区域来限制指定条件的作用（例如，定义区域中使用不同的物质属性）。已经存在的mark可以通过unmark命令取消。1.2.5 指定结构模型和属性一个结构模型是通过model命令与网格区域相关的。在theory and background的第2节中讨论了flac中可用的机械结构模型。下面的关键词可以用来指定适当的模型：机械模型：anisotropic横向等方弹性模型transversely isotropic elastic modelcam-clay修改的cam-clay塑性模型modified cam-clay plasticity modeldruckerdrucker-prager塑性模型drucker-prager plasticity modeldy双曲服塑性模型double-yield (cap) plasticity modelelastic等方弹性模型isotropic elastic modelmohr-coulmohr-coulomb塑性模型mohr-coulomb plasticity modelnull空模型，例没有物质null model -i.e., no materialss软应变塑性模型strain-softening plasticity modelsubiquitous双线性变硬/变软各向接合模型bilinear hardening/softening ubiquitous-joint modelubiquitous各向接合模型ubiquitous-joint model蠕动模型：cvisc挤压蠕动粘塑性模型burger-creep viscoplastic modelcwipp盐变形合并模型crushed-salt consolidation modelpower两组分动力模型two-component power lawpwippwipp蠕动粘塑性模型wipp-creep viscoplastic modelviscous传统粘性模型classical viscositywippwipp参考蠕动模型wipp reference creep model热模型：th_anisotropic各所异性热传导模型anisotropic thermal conductivityth_teneral依靠温度的热传导模型temperature-dependent thermal conductivityth_isotropic等方性热传导模型isotropic thermal conductivityth_null热传导的空区带null zone for heat conduction用户定义的模型：用户可以利用fish语言自己产生结构模型。（参见fish volume的第2.7节）。每个模型通过property命令指定属性。对于软应变、双曲服和双线性模型，属性依靠于通过table命令定义的塑性应力的累积。地下水的属性（密度、体积系数和张力强度）是通过water命令或initial命令（体积模量、张力界线）来指定的。渗透性和多孔性是通过property命令来指定的。1.2.6 指定初始条件问题初始化和模型解决条件是通过下面的命令指定的：initial 初始化一定的网格点和区带变量，例应力状态、孔隙压力、速度和温度。set通过选择一个或多个关键词，让用户来初始化问题和模型条件，例：creeptime设置蠕动打开的时间flow打开/关闭地下水流gravity指定重力implicit对地下水流或热模型打开/关闭隐式计算方法large/small选择大应变或小应变方式mech打开/关闭机械计算方法water table 为了定义压力分布指定一个井表面。（这个命令在应力计算时是有效的，但是当flac被配置成地下水流时是没效的。事实上，井表面是按流体的流动公式计算的。）水的table线是通过table命令来指定的。1.2.7 边界条件的应用flac中模型的边界条件是通过以下命令来描述的：apply施加机械的、地下水流和热条件到任何一个模型边界（外部或内部）。fix/free在选择的网格点上让速度、孔隙压力、饱和度或温度固定（例防止变化）或释放（例允许改变）。ieboundary指定flac网格一个无限弹性边界。在这种情况下网格边界必须是园形的（例如，通过使用generate circle命令）。interior施加机械的、地下水流和热条件到任何一个内部的网格点或区带。1.2.8 指定结构支撑通过structure命令可以指定四种类型的结构支撑。通过下列关键词表示元素类型是可用的：beam 指定一个beam（柱）元素cable指定一个cable（索）元素pile指定一个pile（堆）元素support指定一个support（一维的非线性弹性）元素每个元素类型的属性都是通过propetry关键词来定义的。theory and background中的第5节中有flac结构元素的描述。1.2.9 指定分界面或缝合面在flac网格中两个或更多部分之间定义分界面可以通过interface来实现。这些分界面是平面的，这里滑动和/或分离是允许的。分界面的属性也是通过这些命令来定义。分界面的具体描述参见theory and background的第4节。1.2.10 指定用户定义的变量或函数flac捆绑的语言fish或以调用用户定义的变量或函数。fish语句（参见fish volume 的第2节）可以出现在flac命令的define和end命令之间。opt命令可以优化fish函数。模型条件的变量也能通过table命令使用。1.2.11 在处理过程中监控模型条件在处理过程中可以用命令来监控模型的变化。这对确定是否达到平衡或故障状态是有用的，在计算过程中可以监控变量的变化，例如地下水的流动。history时间步处理过程中，记录变量的变化过程。绘制的结果图能帮助用户判断是否达到了平衡状态。track在流体计算瞬间（要求config gw），记录流动的轨迹。1.2.12 计算问题一旦flac模型的条件都定义好了，问题就可以通过一定的计算步得到解决。下面的命令允许flac模型自动计算或用户控制计算过程：cycle n执行n个时间步solve 自动计算到平衡状态。每次都与预先设置的限制条件进行比较，直至达到限制条件。限制条件可以通过solve或set命令来修改。在瞬时分析过程中例如地下水流动、热的传送和蠕动，和完整的动力学分析过程中，也可以通过关键词来让用户定义时间限制和时间步，step n 执行n个时间步。solve和step(或cycle)命令，最大不平衡压力率（或流体模型下的流体流动率，或热模型下的热变化率）不断地显示在屏幕上。在任何时候用户想中断计算，只需按键。当前的计算步完成后，flac会返回行控制状态；用户可以检查计算结果，保存状态或继续计算。关键词continue可以在适当的命令后给出，以继续执行命令。1.2.13 完成模型输出用户检查当前问题状态可以使用以下的几个命令：copy对上一个屏幕绘图产生一个硬拷贝label增加用户定义的标签到图上movie捕捉屏幕图并将它生成一个快速重放的动画文件（只有dos版本有此功能）plot对于不同的问题绘制图件，显示在屏幕上或硬拷贝设备上print输出问题条件和主网格变量sclin在等高线图上标记等高线set提供不同制图条件的几个控制title为后面的图件指定一个标题window改变显示窗口1.2.14 其它命令help提供可用命令的屏幕列表system允许访问dos命令或特性1.3 flac命令-详细描述按字母顺序对所有的flac命令进行详细的描述。（1）apply keyword range在某个范围中施加某一大小的变量值。apply命令用来施加机械的、地下水和热边界条件到模型网格的任何外部的或内部的边界上。（参见interior命令，它用来施加边界条件到内部网格点或区带）。用户必须指定施加的关键词类型（例xforce），数值（如果需要的话）和施加边界条件的范围。范围可以用两种方式表示：（参见1.1.3节）网格点范围： i=i1,i2 j=j1,j2或边界路径范围： from i1,j1 to i2,j2这里，在边界路径中可以使用3个可选关键词（long, short 或both），参见1.1.3.2节。缺省值是short。如果i1=i2和j1=j2，则表示范围是一个单个的网格点。如果i1=i2, j1=j2和指定了long，则范围是整个边界。apply命令可以使用两个可选的关键词短语：var vx,vy 短语可以在给定的范围内定义一个数值的线性变化。（参见1.1.3.3节的var的解释）。使用hist关键词可以施加一个历史的增加器到一个数值上。历史可以使用以下关键词中三种形式中的一种。n历史增加器是一个输入的历史（参见history命令），这里n是一个历史数字；它假设动力时间(dytime)作为时间轴（只在动力学参数中可用，参见optional features 第3节）。name 历史增加器是一个fish函数，这里name是函数名。table n 历史增加器应用于一个表（参见table命令），这里n表数，它假设动力时间(dytime)作为时间轴（只在动力学参数中可用，参见optional features 第3节）。在机械、地下水或热边界条件下，下面的关键词可以群组使用。例如，users guide的3.3节中给出的apply命令。在关键词后给出了施加边界条件的规则。机械边界条件：ff（自由free-field）、nacc(加速度，正常方向) v、nquiet(静止（粘性），正常方向)、nstress（应力，正常方向）v、nvelocity（速率，正常方向）v、pressure （机械压力）v、sacc (加速度，剪切方向) v、squiet(静止（粘性），剪切方向) 、sstress（应力，剪切方向) v、svelovity（速率，剪切方向） v、sxx（应力张量的xx组分) v、sxy（应力张量的xy组分) v、syy（应力张量的yy组分) v、xacc（加速度x组分，网格点) v、xforce（应力的x组分，网格点) v、xquiet（静止（粘性），x方向)、xtraction（总应力张量，x方向) v、xveocity（速率的x组分，网格点) v、yacc（加速度的y组分，网格点) v、yforce（压力的y组分，网格点) v、yquiet（静止（粘性），y方向) 、ytraction（总应力张量，y方向) v、yvelocity（速率的y组分，网格点) v。地下水边界条件：discharge(flux) v、pp(流体孔隙压力) v。热边界条件：convection v1 v2(v1对流发生时媒体的温度，是v2是热对流的传送系数w/m2)、flux v (w/m2)、radiation v1,v2 (v1是辐射发生时媒体的温度,v2是热福射的传送系数，黑色块体为5.668108w/m2k4)。取消模型边界条件： remove range施加模型边界条件的规则：1每个边界段只能有一个条件施加于每一个自由度。如果apply命令象前面一个命令一样给出了相同的关键词和相同的范围，则值（和/或变化）被更新。有两个例外：（1）与时间相关（例经过一个hist增加器）的压力（应力）被增加到同样范围的常压力（应力）上；（2）静止的边界增加到压力和应力边界。2对于每一段，边界条件必须用术语x-y轴或s-n轴给出。例如，即使边界条件是垂直的，关键词xvel将消去已经存在的svel和nvel条件。一个xvel条件可能施加于y方向条件上（例yvel, ytractuion或yforce）-而不能是一个常规或剪切边界条件。表1.3（略）显示了机械关键词的兼容性。3用应力矢量组分sxx,sxy和syy术语指定的应力边界可施加到两个自由度（独立于边界方向）。应力组分将代替已存存在的边界条件。唯一的例外是，静止的边界可以增加到sxx, sxy和syy边界。4下面的边界条件关键词可应用于边缘，至少需要两个网格点范围：sxx, sxy, syysstr, nstr, pressurextraction, ytractionxquiet, yquiet, squiet, nquiet5下面的边界条件关键词可应用于单个的边界网格点：xvel, yvel, svel, nvelxacc, yacc, sacc, naccxforce, yforce6如果一个“点类型”条件（参见规则5）应用于一个相邻网格点的一个“边缘类型”条件上（参见规则4），点类型条件将被更新，但是边缘类型条件将不受影响。7当静止边界应用到一个已经施加了速率或加速度的网格点上，或一个fix的网格点上，压力的反作用将自动计算，一直到达到静态平衡状态。fix条件被取消。（输入plot apply可以查看压力的反作用）。压力反作用只能被remove 关键词取消（例，后来施加的负载将不取消压力反作用）。如果在静止边界应用前，取消了施加的速率，则不会再有压力反作用。8命令fix和free可以应用于施加于压力或应力的网格点。但是，如果速率或加速度已经施加于直交方向上，则这些命令不能用来固定某一方向的网格点。如果已经施加了速率或加速度，再指定fix网格点，则已经存在于直交方向的fix条件将被取消。9边界路径范围或网格点范围不能通过施加于边界上attach过的网格点，如果施加一个压力、速率或加速度边界的话。attach过的网格点不能被apply命令识别。命令将分别施加于attach的两边的网格点，只有一边可以施加压力、速度或加速度。10当指定倾斜的边界的静止边界条件时，总是同时施加nquiet和squiet。11apply ff命令取消i=1和i=igp（这里igp是i方向的网格点总数）处的网格点已经存在的边界条件，指定自由边界条件。固定的网格点被取消，除了底部网格点（1，1）和（igp，1）的y固定条件。所有的区带的第一个和最后一个区带数据都被拷贝到自由区域。通过网格的第一个和最后下个列的当前值，自由区域的负载、施加的速率和静止边界都是自动更新的。横向边界必须是垂直的，位于i=1和i=igp列的网格点。底部区带j=1 (i=1和i=igp)不能是空的。当apply ff命令给出时，模型将处于一个静态平衡状态。自由区域只是针对动力学分析，在自由区域施加之前，要使用set dyn on命令。在改变其它边界条件到动力状态时之前，必须给出apply ff命令。自由区域条件可以通过以下的命令取消：apply remove i =1apply remove i =igp当命令发出后，左边界和右边界将变成自由的。可以施加静态边界直到模型达到平衡。分界面和附属线不能移到自由区域网格，因此，interface或attach条件不能扩展到自由区域边界。与分界面有相同属性的区带层能重现分界面的效果。（2）attach aside from i1, j1 to i2, j2 bside from i3, j3 to i4, j4attach命令让用户将aside和bside指定的两个边界配属成一对网格点。这个命令对产生线性发射的网格是有效的。详细参见theory and background的第3节。（3）call filename调用一个已经编辑好的数据文件来处理一个模型。call命令中可以包括一个输入文件名。在这个文件中可以输入任何指令，因此flac可以运行在无人的状态。文件可以调用另一个文件，嵌套没有限制。但是，要防止相互递归的调用（例，文件“abc”调用“def”，而“def”调用“abc”），这样结果是不正确的，计算机会出错。call命令执行后，文件中的return命令会返回用户控制状态。按键可以中止从文件中输入，返回用户控制状态。如果call命令执行后，文件中有plot命令，即使用户没有输入return命令，flac也会在显示图形后，返回交互方式。但是，用户可以在plot命令中使用hold关键词来取消绘图。如果movie打开了，则屏幕绘图将会保存到动画文件中。数据文件中插入pause命令可以让用户去检查中间结果（参见pause命令）。文件名应该给出全路径文件名，例: call c:temptest1.dat（4）config keyword 这个命令让用户指定需要分配每个区带或网格点的额外内存计算模式。参数可以是：axisymmetric geometry, groundwater flow, plane stress, heat transfer, fully dynamicanalysis, creep analysis和fish语言的特殊网格变量。如果需要给出关键词，则config命令必须在grid命令之前给出。关键词用法如下：axisymmetry 轴对称creep 蠕动物质分析dynamic 全动力学分析extra n fish语言使用的额外的网格变量n gwflow 地下水流动p_stress 平面应力(只对model elastic, model anisotropic and model ss)thermal 热分析只允许一个config命令，尽管它可以有多个参数（例,config ax th，同时具有轴对称和热参数）。如果总是需要某几个参数，最好将config命令写在flac.ini文件中，这样flac就会自动进行正确的设置。总是使用config extra来适应执行fish函数时额外网格变量是一个好的习惯。有些参数不能同时使用，下面的混合使用是不允许的：axisymmetry and p_stress此外，在一定的模型中有些参数也不能运行。在mohr,ubi或dy模型中不允许有p_stress。在平面应力使用model ss来执行mohr-coulomb塑性分析。如果在model anisotropic 或model ubiquitous模型中使用axisymmetry，则各向异性角度必须等于0。对于纯动力学模型，一般按惯例设为config extra 2即可。对于动力学+热模型，一般按惯例设为config extra 4 thermal即可。（5）continue这个命令让用户继续读一个数据文件。有两种情况可以中止读入数据文件，一种是pause命令（或return命令），一种是命令行出错。如果命令行出错，用户将使用交互输入编辑键来修改输入行来纠正错误。continue命令可以让flac继续读入数据文件中的下一行。注意：这个修改的输入行是暂时的，不会写到源文件中去。（6）copy 执行这个命令将产生一个硬拷贝。如果指定了文件名，则plot命令的输出将直接拷贝到该文件中。如果不指定文件名，plot数据将会送至当前设置的设备中（参见set out和set plot命令）。交互式操作下的缺省设备是lpt1（打印机）。（7）cycle ncycle是step的同义词，这里n是执行计算的步数。如果在执行过程中按键，则当前步完成后会返回用户控制状态。如果循环是通过fish函数驱动的，按一次，将中止计算步，但是函数的下一条命令将执行。按两次将退出回到命令状态（一次是计算，一次是fish循环）。可选的关键词continue表示在用户按了键中止后继续计算。剩下的步数将被执行。如果在一个数据文件中执行cycle命令，则当循环完成后会继续读入剩下的数据文件。（8）define function-name enddefine和end命令用来定义一个fish语言函数。define和end命令之间的语句被编译和以简单的方式保存。语句处理过程中会报告编译错误。这些“源”语句不被flac保存，因此，fish函数可以象数据文件一样准备好，如果发现错误可以进行纠正和修改。fish可以用来产生一个新的变量来打印或绘图、在flac执行过程中控制条件、产生指定分布的属性、以某些特定方式分析flac输出，或产生用户定义的结构模型。fish volume的第2节详细描述了fish的操作和使用。fish volume的第3节包含了fish函数的通用库。（9）fix keyword . 使用这个命令，可以阻止所选网格点上的速度、孔隙度、饱和度或温度产生变化。如果需要固定位移，则相应的速率要设置为0（在开始时0速率是缺省值）。使用apply命令来提供一个刚性的、移动的边界条件。注意，如果给出了mark关键词，则只有网格点范围内被标记的网格点才被固定。可以使用下面的关键词：pp固定孔隙压力（饱和度也必须固定，fix saturation）saturation固定饱和度。饱和度可以被固定为任何存在的值。这对应于水源（sat=1）或者吸水的设备（sat0）。temperature 固定温度。如果给出了value值，则温度固定为这个值。如果没有指定，则value值为0，除非通过initial temp命令来指定一个固定的网格点。x固定x速率y固定y速率网格点的约束可以通过free命令解除。fix命令可以施加到已经指定了压力或应力的网格点（参见apply命令）。但是，如果速率或加速度已经施加在直角方向，则不能使用fix命令在一个方向上固定网格点。如果施加的速率或加速度被指定到固定了的网格点，则在直交方向上已经存在的固定条件将会取消。（10）free keyword 这个命令用来取消通过fix命令设置的约束。如果给出了mark关键词，则只有网格点范围中的被标记的网格点被释放。可以使用的关键词与fix命令的一样。（11）generate orkeyword value generate命令是用在grid命令之后用来设定模型的几何形状的。有两种形式可以使用：网格的外形被坐标区域再分区，网格的外形变成曲线。坐标可以产生整个网格或次区域。如果省略了坐标值，则通过grid命令指定了一个缺省的区带坐标，在坐标之间使用一个单位长度。如果给出了generate的第一种形式，则四边形区域（决定于行和列的范围）被再分区。范围中最小行和列数的坐标值为（x1,y1）；按顺时针方向，剩下的几个角的坐标为（x2,y2），（x3,y3），（x4,y4）。任何一个坐标对可以用same来代替，在这种情况下，特定的角将保持它当前的坐标。四边形必须按顺时针方向设定，内部的角度必须小于180度。四边形区域内的网格点是均匀分布的，除非使用了ratio关键词。使用ratio关键词后，区带的大小可以按照ri和rj参数递增或递减。例如，ratio 1.1, 1.1将使得每个连续的区带在i或j方向都是前一个区带大小的1.1倍。此外，网格还可以被曲线变形。这种形式的generate命令可以指定网格的形状为园形、弧形或线状。这些命令可以用来定义网格的表面形状或内部形状，用它来定义挖空或构造，例层状堆积或矿体。为了这些函数能被执行，首先要指定model命令。在定义形状时可以使用下面的关键词：arc xc,yc xb,yb theta 网格被映射成一段园弧，园心为（xc,yc），起始的边界点为（xh,yh），theta是逆时针方向的角度。circle xc,yc rad 网格被映射成一个园，半径为rad，园心为（xc,yc）。line x1,y1 x2,y2 网格被画成一条线，两个端点分别为（x1,y1）和（x2,y2）。table n 网格点被移动以适应指定的坐标对（参见table命令）。注意，这个命令必须使用一个封闭的路径，第一点和最后一点必须是相同的。注意，这种方式形成的网格点可以被“标记”（参见mark和unmark命令）。如果标记的网格点形成了一个封闭区域，则在initial、model、print和property命令中可以使用region关键词。例如，园形的区域可以用命令model null reg=i,j来挖空，这里i,j是园形内的任一个区带。使用generate命令时有几个注意事项。有些几何形状是不能用四边形元素构建的（例，有直角边的长、薄的三角形区域）。此外，用户应该注意，一旦网格点被标记，它们不能被另一个generate命令移动。如果用户希望移动以前标记过的网格点，首先必须使用unmark命令。最后，可以通过离散化上面形成的网格形状而园滑网格。可以使用下面的关键词：adjust 网格进行自动调整以得到较园滑的网格。标记的网格点和边界点不能移动。连续的generate adjust命令可以使网格更加园滑。在users guide的3.2节中可以发现网格产生的例子。（12）grid icol jrow指定计算网格的总行、列数。其中icol为列数，jrow为行数。能产生的区带总数（icol,jrow）决定于计算机的内存大小和问题的类型（机械、流体力学、热力学），也就是选择的结构模型。当运行mohr-coulomb模型的机械问题时，8m的计算机内存可以产生大约30000个区带。如果config命令带有参数的话，必须在grid命令之前给出。（13）helphelp命令提供一个可用的flac命令的屏幕列表。help也可以象关键词一样在命令中使用。例如，print help将列出print命令可用的所有关键词。“？”可以用来代替help命令或关键词。（14）history (i1,j1)网格点或区带中的网格变量的历史每隔n时间步就会保存下来。历史的顺序是自己从1开始累加。也可以指定一个数值nh。使用print history或history list命令可以在任何时候显示历史请求摘要。nstep对所有的历史是相同的，只能给出一次，并且要在计算步开始之前。n的缺省值是10。变量对时间步或其它的历史可以被绘图（参见plot history命令）。注意，每一个history命令只能给出一个区带（i,j），尽管它可以在同一行中给出几个关键词。注意不要混淆网格点和区带关键词。历史值保存在硬盘中直接存取的文件中。当停止flac时，这个文件被删除，因此，如果用户希望保存历史，必须产生一个save文件或历史必须写到一个已格式化的磁盘文件中（参见history write n命令）。网格点或区带的有效历史关键词是：angle（x轴正向反时针主应力角度）、esxx（有效的xx应力）