FLAC 讲义
一、什么是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
1.4 分析模式大小与RAM之关系
1.5 单位
1.6 正负号方向
(1) 应力－正号代表张力，负号代表压力
(2) 剪应力－详见下图，图中所示剪应力为正号 (3) 应变－正的应变表示伸长，负的应变代表压缩 (4) 剪应变－剪应变的正负号与剪应力相同 (5) 孔隙压力－孔隙压力永远为正
(6) 重力－正号的重力物质往下拉，负号的重力将物质 往上提。
2
二、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)
3
使用者另可使用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程式编写顺序需依实际工程进行之逻 辑步骤建立。如欲分析开挖题目，则需先求出未挖前之应力分布，再以此应力分布求开挖后之土壤变形及是否会崩
4
毁，如因开挖后造成崩毁， 则FLAC程式将无法继续，萤幕将出现ERROR 讯息，如
Bad Geometry Zone ，---------'-----------．
FLAC程式前几行之顺序为
o Config ________ o Grid ________
o 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
FLAC 讲义(含地下水的边坡例子)
