[分享]：FLac 讲义

freestart

看到有人在讨论Flac，我把收集到的一些Flac资料供大家参考：

[B]FLAC　讲义 [/B]

一、什么是FLAC

1.1 FLAC之字义

　　F(Fast)L(Lagrangian)A(Analysis of)C(Continua). Lagrangian相对于Eulerian为每一时阶(timestep)之位 移在Lagrangian之公式中，需对网格之座标予以更新， 而Eulerian之公式则不予更新。

1. 2 FLAC之运算流程

1.3 FLAC 基本单元

1.4 分析模式大小与RAM之关系

1.5 单位

1.6 正负号方向

(1)应力－正号代表张力，负号代表压力

(2)剪应力－详见下图，图中所示剪应力为正号

(3)应变－正的应变表示伸长，负的应变代表压缩

(4)剪应变－剪应变的正负号与剪应力相同

(5)孔隙压力－孔隙压力永远为正

(6)重力－正号的重力物质往下拉，负号的重力将物质 往上提。



二、FLAC内建之组合律

FLAC内建之组合律有：

1． 空洞模式(null model)

使用于土壤被卸载或开挖

2． 弹性模式

3． 塑性模式，

包括

a. Drucker －Prager model

b. Mohr－Coulomb model

c. ubiquitous－joint model

d. strain－hardening/softening model

e. bilinear strain－hardening/softening model

f. double－yield model

g modified cam－clay model

此外，另有选购(option)模式，包括：

1. 动力模式(Dynamic Option)

2. 热力模式(Thermal Option)

3. 潜变模式 (Creep Option)

　　使用者另可使用FISH语言去建构独特的组合律以符 合所需。

三、FLAC－以命令为输入语法

请查阅相关手册



四、FLAC程序之使用步骤

4.1 FLAC程序使用前准备步骤

步骤1：依比例画出所欲分析之资料

于纸上画出地点之位置、地层资料、并简标示距离及深度资料。

步骤2：换算输入资料成同一单位

将现有地层资料，如 Density,Bulk modulus,Young`s modulus, tension,cohesion, friction Angle 等资料，换算成同一单位。　

附注 ：需谨慎检查输入资料之单位，如因单位 不同而造成过大或过小的值，将会造成 FLAC无法计算，而产生ERROR讯息。　 　　　　

步骤3：应用公式简略计算

应用公式或依据经验，简略算出FLAC 输 出资料的范围，以做为 Debug 及输出资 料分析时验证。

步骤4：建立x,y座标与node i,j之关系

于图上距离及深度之关系，建立x,y座标系统，再由x,y座标系统，转换与网格间系，为了便于以后输出资料的分析，故应确实掌握网格之位置及其相对应的x,y 座标。

建议在敏感区域使用较密之网格，其它地 区则使用较疏之网格，刚开始跑程序时， 不宜使用网格太大的网格数目，因尽量使 网格总数少于1000，以节省时间。

4.2 FLAC输入程序编写步骤

　　FLAC程序编写顺序需依实际工程进行之逻 辑步骤建立。如欲分析开挖题目，则需先求出未挖前之应力分布，再以此应力分布求开挖后之土壤变形及是否会崩毁，如因开挖后造成崩毁， 则FLAC程序将无法继续，屏幕将出现ERROR 讯息，如

Bad Geometry Zone ，---------'-----------．

FLAC程序前几行之顺序为

Config　________
Grid　　________
Model 　________　　　　
以后各行编写之顺序，则建议以下之步骤：

1. 求起始之应力平衡

(1)建立x,y座标与网格之关系，建议使用Gen指示：

Gen x1,y1 x2,y2 ,x3,y3 x4,y4 i=i0,i1 j=j0,j1 详细指令参见使用手册，FLAC程序可自动产生x,y 座标与网格之关系，但由于产生之网格座标不易控 制，将对其它之工作产生负面影响，故依使用前步骤 4所建立之关系，将网格依其疏密程度需要之不同， 实际控制网格之座标。

(2)设定材料性质：prop

(3)设定外力：Set Grav, Apply Pressure, ini sxx, Syy

(4)设定边界条件：fix, free

(5)求起始之应力平衡：solve

(6)储存：Save



２求工程之影响

　　求出区域内之应力分布情况后，再依工程之流程及步骤阶段执行各工程进行过程之影响，建议使用以下之步 骤：

(1)叫出起初之应力平衡：re_____ .sav

(2)设定新的材料性质：model,prop

(3)设定新的支撑性质：struct

(4)设定新的外力

(5)设定边界条件

(6)求工程时之应力平衡

(7)储存



五、分析结果之印出及绘制

　　FLAC V.3.4 在分析成果之绘制上，较先前之各版本 有一明显及方便的设计，因为V.3.4（CONSOLE）版 本。

5.1 分析成果绘图

a.直接绘图（不存图档）

flac：set plot windows

flac：plot pen GR YD(欲画出格网及y向变位)

说明:输入上述两行指令后，连结之打印机会直接印出图形， 约占半页A4之纸张，如欲印出全页，则须设定打印机横向 打印。

b.绘图(存图档)

flac：set out YD.EMF(设以下要画之内容档名为.emf)

flac：set plot emf color(设彩色印制)

flac：plot pen GR YD(画出之内容为格网及Y向变位)

说明：YD.EMF可用Word软件叫出并绘图

5.2 印出分析内容或成果数据

flac：set log YD.TXT(YD.TXT为所要打印内容之档名)

flac：print yd i=1,10 j=1,5(设要印出y向位移量)

flac：set log off

说明：输入上述三个指令后，可用Word，Nodepad等软件叫出YD.TXT并打印。



六、分析范例－边坡稳定分析

本范例共分析三种情况即

(1)粒性土壤 C=0 ；s13.sav

(2)C∮土壤C≠0 ∮＝0；s14.sav

(3)考虑水位线 ；s15.sav

其分析网格之建立如下二图所示



--------------------------------------------------------------------------------

01 title

02 SLOPE UNDER GRAVITATIONAL LOAD

03 grid 20,10

04 ;Mohr-Coulomb model

05 m m

06 ;soil properties-note large cohesion to force initial elastic

07 ;behavior for determining initial stress sate. This will prevent

08 ;slope failure when initializing the gravity stresses

09 prop s=.3e8 b=1e8 d=1500 fri=20 coh=1e10 ten=1e10

10 ;warp grid to form a slope:

11 gen 0,0 0,3 20,3 20,0 j 1 4

12 gen same 9,10 20,10 same i 6 21 j 4 11

13 mark i=1,6 j=4

14 mark i=6, j=4,11

15 model null region 1,10

16 ;displacement boundary conditions

17 fix x i=1

18 fix x i=21

19 fix x y j=1

20 ;apply gravity

21 set grav=9.81

22 ;displacement history of slope

23 his ydis i=10 j=10

24 ;solve for initial gravity stresses

25 slove

26 ;save initial stae

27 save sll.sav

28 ;reset displacement components to zero

29 ini xdis=0 ydis=0

30 ;set cohesion to 0

31 prop coh=0

32 ;use large strain logic

33 set large

34 step 200

35 plot hold bo dis xvel min-2.25e-4 max 0 int 2.5e-5 zero

36 save s12.sav

37 step 800

38 plot hold bo dis xvel min –1e-3 max 0 int 2e-4 zero

39 save s13.sav

40 rest sll.sav

41 ini xdis=0 ydis=0

42 prop coh=le4 tens 0.0

43 set large

44 solve

45 save sl4.sav

46 ini xdis=0.0 ydis=0.0

47 ;install phreatic surface in slope

48 water table 1 den 1000

49 table 1 (0,5) (6.11,5) (20,9)

50 def wet den

51 　loop i (1,izones)

52　　loop j (1,jzones)

53　　　if mode(i,j)>1　then

54　　　　xa=(x(i,j)+x(i+1,j)+(i+1,j+1)+x(i,j+1)

55　　　　xc=0.25*xa

56　　　　ya=(y(i,j)+y(i+1,j)+y(i+1,j+1)+y(i,j+1)

57　　　　 yc=0.25*ya

58　　　　 if yc<table (1, xc) then

59　　　　　density(i,j)=1800

60　　　　end if

61　　　end if

62 　　end loop

63　 end loop

64 end

65 wet den

66 apply press 2e4 var 0 –2e4 from 1,4 to 6,6

67 plot den blo water apply

68 step 6000

69 sclin 1 19 0 19 10

70 plot hold bou vel pp

71 save s15.sav

72 return

程序说明

Line 1 － Line 2

本分析之抬头

Line 3

建立20 ╳ 10之格网

Line 4

以；为开始之指令为说明内容，FLAC不读；以后之指令

Line 5

分析模式为模耳－库伦，即Model Mohr

Line 6 － Line 8

说明：底下之土壤系数中之C值及张力强度值故意放大， 以免土壤在初始状况时即产生破坏

Line 9

土壤性质 properties

ｓ Shear Modulus G=0.3 ╳ 10^8 N/m^2

ｂ Bulk Modulus B=1 ╳ 10^8 N/m^2

ｄ density ｄ=1500 kg/m^3

fri friction angle ∮=20^。

Coh Cohesion C=1 ╳ 10^8 N/m^2

Ten Tension Strength T=1 ╳ 10^8 N/m^2

G=E/2(1+v) 　k=E/(3(1-2v)

Line 10 － Line 12

格网座标化

Line 13 － Line 14

以mark指令将格网区域化为两部份

Line 15

挖除由mark区分后包含Region(1,10)的那部份将边坡 〝制造〞出来

Line 16 － Line 19

边界条件之设定 i=1 及 21 之边界只容许上下位移，最 底部之边界则固定

Line 20 － Line 21

加入重力于整个网格系统

Line 22 － Line 23

选择格点(10,10)检核其y向之位移历线

Line 24 － Line 25

开始计算求取答案

Line 26 － Line 27

储存初始平衡之条件。唯一之力为重力。此时之土体内 应力为边坡既有之情况

Line 28 － Line 29

将x向及y向之位移归零。因边坡在未被碰触之前其位 移应为零，但是其应力为力入重力后之自然现象，不应 归零

Line 30 － Line 31

探讨CASE 1，纯粹性土壤之情况。将土壤之property 中 之Cohesion 设为零，其余参数不变

Line 32 － Line 33

使用〝大应变〞之模式，则每一个step其格网座标自 动更新

Line 34

求解。除了使用slove之外，可令step=n，

Line 35

绘出位移的等高线

Line 37

设求解之step=800，以便观察800个step后之情况。 因本CASE中土壤之C=0 ∮=20，很明显此边坡会破坏， step=800已够发现此现象。

Line 40

再叫出 sll.sav 即初始平衡况以便进行CASE 2 之分析

Line 41

再设X及Y向位移为零

Line 42

变更Cohesion为1 ╳ 10^4 N/m^2

Line 47 － Line 49

设定水位线为table 1，水的密度是1000 kg/cm^3

Line 50 － Line 64

使用FLAC去定义function wet___ den，此自设之 function乃在定义水位线下之单位重为 1800/cm^3 。 注意：有几个 loop 就要有几个 end loop，有几个if 就要有几个end if，结尾必定要以end结束。

Line 65

前面之 define wet___ den 是在定义FLAC function之内 容。而此行才是要命令FLAC在此时执行此function .必需先定义出function之内容才能令FLAC执行。

Line 66

加外力于边界上，如在边界内部则使用interio指令。 此行即沿着i=1,4在j=6之边界加上2e4 Newton渐消 减为0之外应力。

Line 69

定出画图之扫瞄线scan line，此行定出一条扫瞄线 连接(19,0)(19,10)两点，只要图中之等高线与此扫瞄 线相交，即会标出，方便看图时之参考。

Line 72

跳出FLAC执行档

lzfat

天！这么复杂？你会，干脆教我们得了，看的眼花花。

freestart

Flac 使用指南之简体中文版

FLAC输入命令

     FLAC 的输入和一般的数值模拟的程序不一样, 它可以用交互的方式从键盘输入各个命令, 也可以写成命令文件, 类似于批处理, 由文件来驱动。
     FLAC 命令大小写一样。所有的命令可以附带若干个关键词和有关的数值。在下面的命令解释中, 只有大写的字母起作用, 小写的字母写不写、写多少个都没有关系。i,j,m 和 n 开始的变量要求整型数, 否则要求实型数。 实型数的小数点可以忽略, 但是整型数不能带小数点。数值间可以用空格隔开, 空格的数目不限, 也可以用下面的分隔符隔开:
                       ( ) , / =
     < > 表示可选的参数, 输入时括号不用输入;
     ... 表示可以有任意个参数。
     由 * 号开始到行末为注释, FLAC 在执行时不理会。

     下面的 FLAC 命令按字母排列。

Apply 关键词 = 数值 <关键词 = 数值 ...> <范围>
          可以有下面的关键词:
          Pressure 压力
          XForce X-方向的力
          YForce Y-方向的力

ATtach 该命令可以将一条线上的结点和另一条线上的结点互相接合在一起, 用以形成复杂的网格形状。

Call 文件名
          写成的命令文件可以用 Call 命令来调用, 命令文件的最后一行必须是RETURN, 以返回到交互方式。命令文件中不能有 CALL 命令本身。

Config 关键词
          FLAC 用以解平面应变问题, 但经过配置命令也可以用于解平面应力问题或轴对称问题。 需要时应在形成网格之前发。关键词有:
          P_STR 平面应力问题
          AX 轴对称问题

CYC n
          该命令同 STEP, 为执行 n 个时步的循环运算。

Fix X <Mark> <范围>
          Y
          X Y
          用此命令可以使 <范围> 内结点的 x- 或 y- 方向的速度保持不变。<范围> 的格式可以是 I = i1,i2, J = j1,j2; i 和 j 何者先输入没有关系。如果要求位移固定, 则必须将速度初始化为零 (开始时速度的 缺省值为零)。FIX 和 INI XV, YV 可以联合使用以提供一个刚性移动的边界条件。如果有 MARK 的关键词, 则只有在此范围内被标记的结点将被固定。

FRee X <Mark> <范围>
          Y
          X Y
          该命令与 FIX 命令相反, 用于放松对结点的约束。

GEn x1,y1 x2,y2 x3,y3 x4,y4 <Ratio ri rj> <range>
         <Same> <Same> <Same> <Same>
          用于在全域或局域中产生网格。(x1,y1)的点放在左下侧, (x2,y2)、
          (x3,y3) 和 (x4,y4) 为顺时针排列。如果点和前一个 GEn 命令中的相同则用 SAME 来代替。如果网格不是均匀排列, 则用关键词 RATIO, 在 I 和 J 方向的比例系数分别为 ri 和 rj.

GEn Circle xc,yc rad
          Arc xc,yc xb,yb theta
          Line x1,y1 x2,y2
          该命令产生园、弧或直线。
          ARC 弧的圆心为 (xc,yc), 起始点为 (xb,yb), 逆时针的角为
                    theta 度
          CIRCLE 园的中心为 (xc,yc), 半径为 rad
          LINE 直线的两端为 (x1,y1) 和 (x2,y2).
          用 ARC 等关键词所形成的边界系被标记 (参见命令 MARK 和 UNMARK)过, 标记结点所围住的区域在以后用命令 INITAL, MODEL 和 PROPERTY 时可以用关键词 REGION 来说明。

GEn ADJust
          用此命令来微调网格, 使之离散平滑。该命令可以连续使用以增进平滑效果的。

Grid icol jrow
          此命令用以产生数目为 icol 列和 jrow 行的网格。对于 640k 内存的FLAC 版本, 如用摩尔-库仑的本构模型约可分两千个单元。

Help 帮助命令, 可在屏幕上显示命令表。

His <Nstep=n> <关键词 . . . I=i1 J=j1>
          每隔 NSTEP 时步, 记录一下关键词所示项目在 i1 和 ji 结点或单元的数值。 NSTEP 的缺省值为 10. 用户应记住所要求 His (历史) 的顺序, 因以后在打印或绘图时要用。历史值在 FLAC 停止运行时删除, 因 此如果要保留的话, 要用到命令 HIS WRITE N (见下面)。
          可以用的关键词有:
          Ang 单元 i,j 内的最小主应力和 X-轴所形成的夹角
          PP 单元 i,j 内的孔隙压力
          SIG1 单元 i,j 的最大主应力
          SIG2 单元 i,j 的最小主应力
          SXX 单元 i,j 的 xx-应力
          SYY 单元 i,j 的 yy-应力
          SXY 单元 i,j 的 xy-应力
          X 结点 i,j 的 x-坐标
          Y 结点 i,j 的 y-坐标
          XDis 结点 i,j 的 x-位移
          YDis 结点 i,j 的 y-位移
          Unbal 最大不平衡力
          XVel 结点 i,j 的 x-速度
          YVel 结点 i,j 的 y-速度
          XXA 单元 i,j 三角形 a 的 xx-应力
          XYA 单元 i,j 三角形 a 的 xy-应力
          YYA 单元 i,j 三角形 a 的 yy-应力
          XXB 单元 i,j 三角形 b 的 xx-应力
          XYB 单元 i,j 三角形 b 的 xy-应力
          YYB 单元 i,j 三角形 b 的 yy-应力
          XXC 单元 i,j 三角形 c 的 xx-应力
          XYC 单元 i,j 三角形 c 的 xy-应力
          YYC 单元 i,j 三角形 c 的 yy-应力
          XXD 单元 i,j 三角形 d 的 xx-应力
          XYD 单元 i,j 三角形 d 的 xy-应力
          YYD 单元 i,j 三角形 d 的 yy-应力
His 关键词
          关键词可以是:
          Dump nhis 将第 nhis 的历史写屏;
          Write nhis 将第 nhis 的历史写在文件 FLAC.HIS 上, 该文件可在 FLAC 结束后打印出来。连续执行 HIS WRITE 命令可将结果顺序写在 FLAC.HIS 上, 但是首次执行此命令会将以前盘上所存的同名文件冲掉;
          Reset 所有的历史都清除掉。

Initial 关键词 = 值 <. . . > <范围>
          某些结点值可以给初值:
               PP 空隙压力
               X x-坐标
               Y y-坐标
               SXX xx-应力
               SYY yy-应力
               SXY xy-应力
               XDisp x-位移
               YDisp y-位移
               XVel X-速度
               YVel y-速度
          还有几个可选用的关键词以协助实现 INITIAL 命令:
          Mark 只有标记的结点值将被初始化
          Region i j
                    为标记结点所围住的区域将被初始化。i,j 为标记区域中的任一个单元
          Var xv, yv
                    在一定范围内的参数的变差值, xv 和 yv 为 x-方向和 y-方向的变差值。
          范围可以有形式 I=i1,j1 J=j1,j2, 其中 i 和 j 的先后次序无关。

INTerface
          n 关键词 <i1,j1> <i2,j2>
          n 关键词 = 值
          INTERFACE 为交界面命令。网格的一部分可以和网格的另一部分通过交界面而相互作用。界面的性质由其刚度, 粘结力和摩擦力来表征。用于本命令的关键词有:
          Aside i1,j1 <i2,j2 i3,j3 ...>
          Bdide i1,j1 <i2,j2 i3,j3 ...>
          Cohesion 粘结力值
          Friction 摩擦力值
          Glued
          KN 法向刚度值
          KS 切向刚度值
          TBond 抗拉强度
          Unglued
          胶结命令 GLUE 将使交界面的上下两部分胶结在一起, 使之不能滑动和分离, 命令 UNGLUE 则为命令 GLUE 之逆。

MArk <范围>
          在给定范围内的结点将被标记。计算的进行与结点的标记与否无关, 但标记的结点可以限定一个区域作为命令 INITIAL, PROP 和 MODEL 的作用域。命令 GEN 将自动对结点进行标记。<范围>的形式可以是 I= i1, i2, J=j1,j2, i 和 j 的顺序无关。

Model 关键词 <REGION=i,j> <范围>
          该命令对于给定的区域或范围赋予有关的本构模型。关键词有:
          ANisotropic 横观各向同性弹性模型
          Elastic 弹性, 各同性本构模型
          Mohr-Coulomb 摩尔-库伦塑性模型
          Mull 零模型, 用于开挖掉的单元
          SS 应变软化
          Ubiquitous 彻体节理模型

NEW 该命令可以不用退出 FLAC 重新开始一个新问题。

Plot 关键词 <开关 <= v>...> <关键词 ...>
          用缺省值。可以在一行上连写几个关键词以在一个图上画出几个变量。在绘图前先要用 MODEL 命令赋予网格以一定的本构模型。 关键词及其意义为:BEam画出结构单元的几何形状
          Boundary 画出网格的内外边界
          Cable 画出锚束的几何形状
          Disp 位移矢量
          E_p 塑性应变的轮廓线 (只限于应变软化材料)
          Grid 画出
                    画出纪录在第 nhis 个历史上的变量值
          PP 孔隙压力的等值线
          RF 以矢量形式画出的固定结点的反力
          STAte 画出单元中心当前的屈服状态
          STress 主应力矢量
          SXX xx-应力的等值线 (全应力)
          SYY yy-应力的等值线 (全应力)
          SXY xy-应力的等值线
          Velocity 以箭头画出的速度矢量
          WAter 水位线
          XDisp x-位移的等值线
          YDisp y-位移的等值线
          XVel x-速度
          YVel y-速度
          所谓开关其本身也是关键词, 用于设置图形的的某些特征, 有:
          color 如红为 Red, 绿为 GREen, 黄为 Yellow 等。
          Interval = c
                    将等值线的间距置为 c. 在 PLOT 命令前发扫描线命令SCLIN 可以画出等值线的值。
          Max = v 在画矢量时, 置箭头的最大长度为 v.
          Noh 不画标题
          Zero 不画零值线
          如果要放大或缩小图形, 可以在 PLOT 命令前用 WINDOW 命令来设置窗口的大小。

Print 关键词 <关键词>... <Region=i,j> <范围>
          和绘图命令 PLOT 一样, 结点变量的值只有在给定的材料模型和性质后才可以打印出来。关键词有:
          Apply 打印出所施加力或压力的大小和范围
          Fix 打印出固定 x 或 y 的结点
          Limits 对 SOLVE 命令的限制值
          MArk 打印出标记结点
          MEm 占用的内存
          STruct 打印出结构单元上有关结点的力, 力矩和位移
          Interface 打印出交界面的数据, 包括结点力和单位法向矢量
          主要的网格关键词
          Bulk 体积模量
          COhesion 粘结力
          Den 质量密度
          Dilation 剪胀角
          E_p 全塑性应变 (应变软化模型)
          Friction 摩擦系数
          SHear 剪切模量
          X X-坐标
          Y Y-坐标
          PP 孔隙压力
          JFric 节理摩擦力
          JCoh 节理粘结力
          JAngle 节理角度
          NUYx y-x 泊松比
          NUZx z-x 泊松比
          SIG1 最大主应力
          SIG2 最小主应力
          SXX XX-应力
          SYY YY-应力
          SXY XY-应力
          Theta 最小主应力与 x-轴的夹角
          XDis X-位移
          YDis Y-位移
          XMod X-模量
          YMod Y-模量
          XVel X-速度
          YVel Y-速度
          ASXX XX-应力 (三角形 A)
          BSXX XX-应力 (三角形  
          CSXX XX-应力 (三角形 C)
          DSXX XX-应力 (三角形 D)
          ASYY YY-应力 (三角形 A)
          BSYY YY-应力 (三角形  
          CSYY YY-应力 (三角形 C)
          DSYY YY-应力 (三角形 D)
          ASXY XY-应力 (三角形 A)
          BSXY XY-应力 (三角形  
          CSXY XY-应力 (三角形 C)
          DSXY XY-应力 (三角形 D)
          State 塑性状态
                    0 弹性
                    1 正在屈服中
                    2 曾经屈服, 现为弹性
                    3 已超过单轴抗张力
                    4 屈服并超过单轴抗张力
                    5 已经超过抗张力
                    6 彻体节理正在屈服
                    7 彻体节理过去屈服过, 现为弹性状态
          Tables 打印所存贮的表格
          XReaction X-反力
          YReaction Y-反力

PROp 关键词 =值 < . . . > <范围>
          本命令为 MODEL 命令赋予 材料的性质。下面的关键词对 PROPERTY 命令加以补充:
          Region i,j
                    所有为标记结点所连续包围的区域将赋予该性质。i,j 单元可为标记区域中的任一单元
          Var xv yv
                    性质可以有变差, xv 和 yv 分别为 x-方向和 y- 方向的变差值
          <范围> 可为 I=i1,i2, J=j1,j2, i 和 j 孰者在前无妨, 但 <范围>必需在输入行的最后。
          各种模型所需输入的性质为:
          弹性
               (1) 剪切模量
               (2) 体积模量
               (3) 密度
          摩尔-库伦
               (1) 剪切模量
               (2) 体积模量
               (3) 密度
               (4) 摩擦角
               (5) 粘结力
               (6) 剪胀角 (任选)
          横观各向同性
               (1) 剪切模量
               (2) x-模量
               (3) y-模量
               (4) 密度
               (5) NUYx
               (6) NUYz
          彻体节理
               (1) 剪切模量
               (2) 体积模量
               (3) 密度
               (4) 粘结力 (整体材料)
               (5) 摩擦力 (整体材料)
               (6) 节理粘结力
               (7) 节理摩擦力
               (8) 节理角
               (9) 整体材料的剪胀角 (任选)
          应变-硬化/软化
               (1) 剪切模量
               (2) 体积模量
               (3) 密度
               (4) 初始摩擦角 (任选)
               (5) 初始粘结力 (任选)
               (6) 摩擦角和塑性应变的关系表的表号
               (7) 粘结力和塑性应变的关系表的表号
               (8) 剪胀角和塑性应变的关系表的表号
          如果表号给 0, 则取关键词 COHESION, DILATION 或 FRICTION 所给出的值。

Quit 该命令同 STOP, FLAC 终止运行。

REstore 文件名
          把以前用 SAVE 命令所存的文件读入内存恢复现场。

RETurn 该命令应是输入数据文件中的最后一个命令。

SAve 文件名
          将内存中问题的现场存入文件中, 如果已经有同名的文件, 则该文件将被复盖。

SClin n x1,y1 x2,y2
          该命令在屏幕上作一个扫描线, 以切割等值线, 交点以 A 到 Z 的字母来表示。在一张图上最多可以有 5 条扫描线。扫描线的参数为:
               n = 扫描线的编号 (必需为 1,2,3,4 或 5)
           x1,y1 = 扫描线的起点坐标
           x2,y2 = 扫描线的终点坐标

[B]Flac的使用最好参考随机PDF电子文档[B]
附件为Flac2d的介绍，E文

zuicai

哇哇哇，好棒，就需要这个资料，谢！

cdut2003

呵呵，好打消积极性哦。有程序、有说明（即所谓的楼上同志的Flac 使用指南之简体中文版 ），一分没有给，楼上的兄弟却给了20个币。

唉，斑竹同志应一视同仁嘛。
算了，flac3d不用上传了。

http://www.xdcad.net/forum/showt ... 45&pagenumber=2

lzwge

OH！My ghost！那么复杂，而且上传这东东（f2d403.rar）又是英文版，看来又得加强学习英语了！

cad517

最初由 cdut2003 发布
呵呵，好打消积极性哦。有程序、有说明（即所谓的楼上同志的Flac 使用指南之简体中文版 ），一分没有给，楼上的兄弟却给了20个币。

cdut2003，已经给你添加：）
好的东西一定要嘉奖！（你的帖子漏看了）

如果你觉得自己的帖子可以加币或者加分，请给我发悄悄话。
请看下贴
http://www.xdcad.net/forum/showthread.php?threadid=50829