1、FLAC软件教学,第一讲 FLAC简介,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 的先后次序无关。,Aside i1,j1 Bdide i1,j1 Cohesion 粘结力值Friction 摩擦力值Glued
2、 KN 法向刚度值KS 切向刚度值TBond 抗拉强度,Unglued 胶结命令 GLUE 将使交界面的上下两部分胶结在一起, 使之不能滑动和分离,命令 UNGLUE 则为命令 GLUE 之逆。 MArk 在给定范围内的结点将被标记。计算的进行与结点的标记与否无关,但标记的结点可以限定一个区域作为命令 INITIAL,PROP 和 MODEL 的作用域。命令 GEN 将自动对结点进行标记。的形式可以是 I=i1,i2, J=j1,j2, i 和 j 的顺序无关。 Model 关键词 该命令对于给定的区域或范围赋予有关的本构模型。关键词有:,ANisotropic 横观各向同性弹性模型 Elas
3、tic 弹性, 各同性本构模型 Mohr-Coulomb 摩尔-库伦塑性模型 Null 零模型, 用于开挖掉的单元 SS 应变软化 Ubiquitous 彻体节理模型 NEW 该命令可以不用退出 FLAC 重新开始一个新问题。,Plot 关键词. 用缺省值。可以在一行上连写几个关键词以在一个图上画出几个变量。在绘图前先要用 MODEL 命令赋予网格以一定的本构模型。关键词及其意义为:BEam画出结构单元的几何形状 Boundary 画出网格的内外边界 Cable 画出锚束的几何形状,Disp 位移矢量 E_p 塑性应变的轮廓线 (只限于应变软化材料) Grid 画出 画出纪录在第 nhis 个
4、历史上的变量值 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. 在
5、 PLOT 命令前发扫描线命令SCLIN 可以画出等值线的值。 Max = v 在画矢量时, 置箭头的最大长度为 v. Noh 不画标题Zero 不画零值线,如果要放大或缩小图形, 可以在 PLOT 命令前用 WINDOW 命令来设置窗口的大小。 Print关键词. 和绘图命令 PLOT 一样, 结点变量的值只有在给定的材料模型和性质后才可以打印出来。关键词有:,Apply打印出所施加力或压力的大小和范围 Fix 打印出固定 x 或 y 的结点 Limits 对 SOLVE 命令的限制值 MArk 打印出标记结点 MEm 占用的内存 STruct 打印出结构单元上有关结点的力, 力矩和位移,I
6、nterface 打印出交界面的数据, 包括结点力和单位法向矢量 主要的网格关键词 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
7、Y-位移 XMod X-模量 YMod Y-模量 XVel X-速度 YVel Y-速度 ASXX XX-应力 (三角形 A) BSXX XX-应力 (三角形 B),CSXX XX-应力 (三角形 C) DSXX XX-应力 (三角形 D) ASYY YY-应力 (三角形 A) BSYY YY-应力 (三角形 B) CSYY YY-应力 (三角形 C) DSYY YY-应力 (三角形 D) ASXY XY-应力 (三角形 A) BSXY XY-应力 (三角形 B),CSXY XY-应力 (三角形 C) DSXY XY-应力 (三角形 D) State 塑性状态0 弹性1 正在屈服中2 曾经屈服
8、, 现为弹性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 孰者在前
9、无妨, 但 必需在输入行的最后。 各种模型所需输入的性质为:,弹性 (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)
10、初始粘结力 (6) 摩擦角和塑性应变的关系表的表号 (7) 粘结力和塑性应变的关系表的表号 (8) 剪胀角和塑性应变的关系表的表号 如果表号给 0, 则取关键词 COHESION, DILATION 或 FRICTION 所给出的值。,Quit 该命令同 STOP, FLAC 终止运行。 REstore 文件名 把以前用 SAVE 命令所存的文件读入内存恢复现场。 RETurn 该命令应是输入数据文件中的最后一个命令。 SAve 文件名将内存中问题的现场存入文件中, 如果已经有同名的文件, 则该文件将被复盖。,SClin n x1,y1 x2,y2 该命令在屏幕上作一个扫描线, 以切割等值线,
11、 交点以 A 到 Z 的字母来表示。在一张图上最多可以有 5 条扫描线。扫描线的参数为: n=扫描线的编号(必需为 1,2,3,4 或 5) x1,y1 = 扫描线的起点坐标 x2,y2 = 扫描线的终点坐标,SEt 该命令设置整体条件, 其关键词有: Aspect a 调整屏幕的高宽比, 以使输出不变形Col nn 为打印输出的最大列数。用户必需在运行 FLAC 前设置好打印机的列数。,EGA 设置 640x350 的 EGA 图形模式 Force f 设置 SOLVE 命令时的不平衡力的限值。 Gravity g 重力加速度 g 及其与负 y-轴的夹角Large 大变形 (坐标更新),Lo
12、g ON 在当前盘上打开名为 FLAC.LOG 的文件, 以纪录 FLAC 的运行进程。如果 FLAC.LOG 文件已经存在, 则将被复盖。 OFF 关闭 FLAC.LOG 文件。如果在稍后再 SET LOG ON, 则屏幕显示将继续纪录在 FLAC.LOG 文件上。,Small 小变形 (坐标不更新) STep n SOLVE 命令进行计算循环的限制数目。 Time t SOLVE 命令进行计算的限制时间的分钟数。,SOLve 用本命令进行计算循环, 关键词为:Step 时步数目Time 运行时间 (分钟数)Force 不平衡力 缺省值为: S = 500 时步, T = 5 分钟, F =
13、 100,最长的时限为 1440 分钟 (24 小时), 如果要求运算超过一天, 则可以用几个SOLVE 的命令。FLAC 在运行中如要中止, 可以按退出键 Esc 退出循环, 将控制转给用户, 进行交互操作。 STEp n执行 n 个时步的循环运算。,STop FLAC 终止运行。 STRuct 关键词 . . . STRUCT 命令用于定义结构单元的几何形状, 性质等条件。STRUCT 命令要求说明结构单元的类型和形状以及它与 FLAC 网格的连接情况。定义结构单元的关键词有:,(1) 结构类型关键词梁 一个梁单元系在关键词 BEAM 后用梁的几何形状和性质的关键词及其值来表征。锚束 一个
14、锚束单元系在关键词 CABLE 后以其几何形状和性质的关键词及其值来表征。,(2) 结构单元的几何形状及其结点联结 关键词 BEGIN 表征梁或锚束的起始点 关键词 END 表征梁或锚束的终了点 结构单元的端点可以在 BEGIN 和 END 关键词后直接用下面三个之一的关键词来表征:,Grid i,j 结构单元的起点或终点为结点 i,j Node n FLAC 将自动按照构件单元输入的顺序为结构单元的结点编号, 此处n 为所编的第 n 号结点。x,y 直接给端点以 x 和 y 坐标。,(3) 结构单元的性质类型关键词Prop n 结构单元的类型要赋予性质数n,用 struct prop=n 再
15、加上关键词及其值来说明。 (4) 结构单元分段的关键词Seg k将结构单元等分成 k 个分段。,(5) 锚束预加张力的关键词Tension t STRuct Node=n 关键词 每个结构结点用下列的关键词可以有不同的边界或几何条件。 FFIx 可以将 n 结点的 x- 和/或 y- 位移或转动固定。,FRee 此命令为命令 FIx 之逆。Pi 使结点 n 为铰接, 即放松力矩。Load fx, fy, m可对结点 n 施加 x- 和/或 y- 方向的力或力矩。x, y n 结点的 x-坐标和 y-坐标。,STRuct Prop=n 关键词 = 值 性质的关键词有:E = 值I = 值 (梁的
16、惯性矩)Area = 断面值 KBond = 灰浆的粘结刚度 SBond = 灰浆的粘结强度 Yield = 锚束的屈服强度 (力),SYstem 不退出 FLAC 暂时进入 DOS 系统, 可以执行的DOS 命令有: DIR , COPY, REN, DEL 和 TYPE. 用回车键可以回到 FLAC. Table n x1 y1 建成一个点的坐标的表, 以供 FLAC 引用。n 是表号。,TItle 在标题命令 TITLE 的下一行输入标题的名称。 Unmark 去掉在给定范围内结点的标记。 WAter 引入水面下的孔隙压力, 孔隙压力为正。命令 WATER TABLE =0 将撤去所有结
17、点的孔隙压力。,WIndow 该命令按照 xlo (x-低坐标), xhi (x-高坐标) 和 ylo (y-低坐标), yhi (y-高坐标) 建立一个窗口, 否则为自动建立 (AUTO).,第二讲 FLAC实例,FLAC 基本单元,分析模式大小与RAM之关系,Flac3.3（解密）计算规模（不限） 内存64MB 弹性,各同性本构模型模型 360万个单元 摩尔-库伦塑性模型 302万个单元,Flac3.3（解密）计算规模 内存128MB 弹性,各同性本构模型模型 722万个单元 摩尔-库伦塑性模型 562万个单元 内存180MB 弹性,各同性本构模型模型 1000万个单元 摩尔-库伦塑性模型
18、 810万个单元,计算机计算规模的确定,flacw_ sp 180 Gr 1000,1000 M e Pr mem,Flac3.3（解密）使用方法,首次运行时，点击flac33.bat，按提示输入：960510，就ok了,改批处理文件,flac33.bat : echo off setback date 5-10-96 copy flaccfg.org flac.cfg flac 64 setoday,Gr 900,900 M m Pr mem,单 位,正负号方向,(1)应力正号代表张力，负号代表压力 (2)剪应力详见下图，图中所示剪应力为正号,(3)应变正应变表示伸长，负应变代表压缩 (4)
19、剪应变剪应变的正负号与剪应力相同 (5)孔隙压力孔隙压力永远为正 (6)重力正号的重力物质往下拉，负号的重力将物质 往上提。,FLAC 模型,1空洞模式(null model) 使用于土壤被移除或开挖 2 弹性模式 3 塑性模式，包括 a. Drucker Prager model b. MohrCoulomb model c. ubiquitousjoint model d. strainhardening/softening model e. bilinear strainhardening/softening model f. doubleyield model g modified c
20、amclay model,此外，另有选购(option)模式，包括： 1. 动力模式(Dynamic Option) 2. 热力模式(Thermal Option) 3. 潜变模式 (Creep Option) 使用者另可使用FISH语言去建构独特需要的模型。,FLAC使用前准备步骤,步骤1：依比例画出所欲分析之资料 于纸上画出地点之位置、地层资料、并简标示距离及深度资料。 步骤2：换算输入资料成同一单位 将现有地层资料，如 Density,Bulk modulus,Youngs modulus, tension,cohesion, friction Angle 等资料，换算成同一单位。 附注
21、 ：需谨慎检查输入资料之单位，如因单位 不同而造成过大或过小的值，将会造成 FLAC无法计算，而产生ERROR讯息。,步骤3：应用公式简略计算 应用公式或依据经验，简略算出FLAC 输 出资料的范围。 步骤4：建立x,y座标与node i,j之关系 于图上距离及深度之关系，建立x,y座标系统，再由x,y座标系统，转换与网格间系，为了便于以后输出资料的分析，故应确实掌握网格之位置及其相对应的x,y 座标。 建议在敏感区域使用较密之网格，其它地 区则使用较疏之网格。,FLAC 输入数据文件编写步骤,FLAC 只认文本文件,可以用记事本、EDIT、PE2等编辑，但不能用word文件,FLAC数据文件
22、编写顺序需依实际工程进行之逻 辑步骤建立。 如欲分析开挖题目，则需先求出未挖前之应力分布， 再以此应力分布求开挖后之土壤变形及是否会崩毁， 如因开挖后造成崩毁， 则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 座标与网格之关系，将网格
23、依其疏密程度需要之不同， 实际控制网格之座标。,(2)设定材料性质：prop (3)设定外力：Set Gravity, Apply Pressure, ini sxx, Syy (4)设定边界条件：fix, free (5)求起始之应力平衡：solve (6)储存：Save,求工程之影响 求出区域内之应力分布情况后，再依工程之流程及步骤阶段执行各工程进行过程之影响，建议使用以下之步 骤： (1) 调出起初之应力平衡状态：re_ .sav (2)设定新的材料性质：model,prop (3)设定新的支撑性质：struct (4)设定新的外力 (5)设定边界条件 (6)求工程时之应力平衡 (7)储
24、存,分析结果输出及图形绘制,分析成果绘图 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软体叫出并绘图 印出分析内容或成果数据 flac：s
25、et log YD.TXT(YD.TXT为所要列印内容之档名) flac：print yd i=1,10 j=1,5(设要印出y向位移量) flac：set log off 说明：输入上述三个指令后，可用Word，Nodepad等软体调出YD.TXT并打印。,例 一,;call ex1.txt grid 6,5 m m plot grid number plot grid gnum,;call ex2.txt ;Function to calculate K and G from E and nu. ;Remember to set e_mod and p_ratio !def DeriveK
26、Gs_mod=e_mod/(2.0*(1.0+p_ratio)b_mod=e_mod/(3.0*(1.0-2.0*p_ratio) end,;call ex3.txt set e_mod=5e8 p_ratio=0.25;call ex4.txt DeriveKG;call ex5.txt property dens=1000 bulk=b_mod shear=s_mod,;call ex6.txt;sxx=0.5*syy def k0_setloop i (1,izones) loop j (1,jzones)sxx(i,j)=0.5*syy(x,y)end_loopend_loop End
27、,例 二,; A Simple Analysis- Failure of a Sandy Slope;call s1.txt grid 40 20 model mohr gen -15,-16 -15,20 60,20 60,-16;call s2.txt table 1 -15,0 0,0 34.641,20 60,20 gen table 1,;call s3.txt model null reg 1,15;call s4.txt title Slope Stability Analysis plot grid,;call s5.txt fix x i=1 fix x y j=1 fix
28、x i=41 prop bulk=2e8 she=1e8 den=2500 prop fric=45 coh=1e6 ten=1e6 set grav=9.81 hist unbal hist xd i=25 j=19 hist yd i=25 j=19 solve,;call s6.txt prop coh=0 ten=0 solve save save sl_ini.sav,;call s7.txt ;Equation 1:FoS=tan(actual friction angle)/tan(friction angle at failure) def ssolvefric_act=fri
29、ction(1,1)loop k (1,int_fric)fname=sl_+string(k)+.savfric_fac=fric_inc*(float(k)-1.0)fric_cal=fric_ref-fric_facFoS =tan(fric_act*degrad)/tan(fric_cal*degrad),commandini xdis=0 ydis=0ini xvel=0 yvel=0print fric_calprop fric=fric_calsolve step 4000end_command xtable(1,k)=fric_calytable(1,k)=xvel(25,19
30、)xtable(2,k)=fric_calytable(2,k)=yvel(25,19),commandsave fnameend_command ;stop series of runs if unbalance force 100if unbal 100 thencommandprint fric_cal FoSquitend_commandend_ifend_loop end,;call s8.txt set fric_ref=31 fric_inc=0.2 int_fric=12 ssolve;call s9.txt rest sl_9.sav print table 1 print
31、table 2 window 10,50 -10 40 plot boun yvel zero,例 题,例3 ;Modeling collapse of an unsupported trench ; create model new grid 5,5,; assign material model and properties model mohr prop bulk=1e8 shear=.3e8 fric=35 prop dens=1000 coh=1e10 ten=1e10 ; print and plot grid print x y,plot grid,;assign boundar
32、y conditions and apply gravity fix y j=1 fix x i=1 fix x i=6 set grav=9.81 Plot gr fix,; monitor history of movement his nstep=5 his ydis i=2 j=6 ; solve for the initial equilibrium state solve ; plot and print the results plot his 1,plot his 1,title a simple trench excavation example sclin 1 (1,0)
33、(1,5) plot hold syy yel bou gre plot hold sxx yel bou gre print sxx syy,save trench.sav ; save the initial equilibrium state ; change the material strength and excavate the trench prop coh=0 set large init xdis=0 ydis=0,model null i=3 j=3,5,; calculate the progressive failure of the trench ; and plo
34、t the results pause step 100 plot hold plastic boundary,plot hold xv yell int=5e-6 dis red max=1e-2 bou green,step 400 plot hold grid,plot hold plas bou,plot hold syy zero int=2500 disp max=0.2 mage bou green,now we make a movie of the progressive failure ; start from the initial equilibrium state ;
35、 and excavate the trench as before pause rest trench.sav,prop coh=0 set large Ini xdis=0 ydis=0 model null i=3 j=3,5 ; use a FISH function to capture the movie plots frames automatically,def showif n_step = 0 thenstep_int = 25elsestep_int = n_stependifloop k(1,snap_shot)command,step step_intplot xd
36、z int=.03 max .25 min -.25 fill dis red max=0.5 bou iwhendcommandendloop end ; make 30 plots at a 30 step interval for the movie,set snap_shot = 30 n_step = 30 movie on pause show movie off ; view the movie ; press the space bar to pause the movie,; press the space bar again to continue the movie ;
37、press the Enter key to exit the movie pause movie view ; now install horizontal struts to support the excavation,; start from the initial equilibrium state ; and excavate the trench as before rest trench.sav prop coh=0 set large init xdis=0 ydis=0 hist reset,his xdis i=3 j=4 model null i=3 j=3,5 ; i
38、nstall three struts across the trench struct beam begin grid 3,4 end grid 4,4 seg=2 struct beam begin grid 3,5 end grid 4,5 seg=2,struct beam begin grid 3,6 end grid 4,6 seg=2 struct prop 1 e=200e9 i 2.3e-5 area 4.8e-3 step 400 ; plot the results: the excavation is now stable,plot str disp blue 25 s
39、truct beam,plot hold hist 1,分析范例边坡稳定分析,分析网格之建立如下二图所示,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 init
40、ializing 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
41、 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
42、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 4
43、5 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
44、,j+1)+y(i,j+1) 57 yc=0.25*ya 58 if yctable (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,Lin
45、e 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 108 N/m2 Bulk Modulus B=1 108 N/m2 density =1500 kg/m3 fri friction angle =20。 Coh Cohesion C=1 108 N/
46、m2 Ten Tension Strength T=1 108 N/m2 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 开
47、始计算求取答案 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
48、 设求解之step=800，以便观察800个step后之情况。 因本CASE中土壤之C=0 =20，很明显此边坡会破坏， step=800已够发现此现象。 Line 40 再叫出 sll.sav 即初始平衡况以便进行CASE 2 之分析 Line 41 再设X及Y向位移为零 Line 42 变更Cohesion为1 104 N/m2 Line 47 Line 49 设定水位线为table 1，水的密度是1000 kg/cm3,Line 50 Line 64 使用FLAC去定义function wet_ den，此自设之 function乃在定义水位线下之单位重为 1800/cm3 。 注意：有几个 loop 就要有几个 end loop，有几个if 就要有几个end if，结尾必定要以end结束。 Line 65 前面之 define wet_ den 是在定义FLAC function之内 容。而此行才是要命令FLAC在此时执行此function .必需先定义出function之内容才能令FLAC执行。,
