FLAC-3D讲义.docVIP

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.4 分析模式大小与RAM之关系 1.5 单位 1.6 正负号方向 (1)应力－正号代表张力，负号代表压力 (2)剪应力－详见下图，图中所示剪应力为正号 (3)应变－正的应变表示伸长，负的应变代表压缩 (4)剪应变－剪应变的正负号与剪应力相同 (5)孔隙压力－孔隙压力永远为正 (6)重力－正号的重力物质往下拉，负号的重力将物质 往上提。 ? 二、FLAC内建之组合律 FLAC内建之组合律有： 空洞模式(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)设