flac--学习.docVIP

FLAC原理实例与应用指南(附光盘一张)FLAC——Fast Lagrangian Analysis of Continua FLAC建立在拉格朗日算法基础上，采用有限差分显式算法来获得模型全部运动方程(包括内变量)的时间步长解，从而可以追踪材料的渐进破坏和垮落，这对研究岩土工程设计是非常重要的。 FLAC适用模拟计算岩土材料力学行为，特别适合模拟大变形和扭曲，包括材料的高度非线性(应变硬化/软化)、不可逆剪切破坏和压密、粘弹(蠕变)、孔隙介质的应力—渗流耦合、热—力耦合以及动力学问题等。 2. FLAC软件应用 2.1 建模 （1）设计计算模型的尺寸 （2）规划计算网格数目和分布 （3）安排工程对象（开挖、支护等） （4）给出材料的力学参数 （5）确定边界条件 2.2 网格生成 Grid i,j 例如：grid 30,20 2.3 网格规划：Gen x1,y1 x2,y2 x3,y3 x4,y4 例如：Gen 0,0 0,10 10,20 20,0 2.4 赋给单元材料性质 mod e (弹性) prop d 1800e-6 bu 12.5 sh 5.77 i=1,20 j=1,10 prop d 2400e-6 bu 1250 sh 577 i=1,20 j=11,20 mod m （弹塑性Mohr-Coulumb准则） prop d 1800e-6 bu 12.5 sh 5.77 c 0 fri 20 ten 0.015 reg i,j 2.5 赋给模型边界条件 （1）固定边界 （结点） Fix x i=1, j=1,21 Fix y i=1,21 j=1 （2）施加边界力 (结点) apply yf -10 i=1,21 j=21 或 apply syy -10 i=1,21 j=21 apply xf -5 i=21, j=1,21 或 apply sxx -5 i=21, j=1,21 2.6 赋给模型初始条件 赋单元内应力 （单元） ini sxx -10 i=1,20 j=1,20 ini syy -5 var 0 4 i=1,21 j=1,21 2.7 计算 Set grav 9.81 Set large Step 1000 Save test.sav 2.8 结果显示 Plot grid 显示网格 Plot bo 显示边界 Plot plas 显示塑性区 Plot sig1 fi 显示最大主应力?1 Plot sig2 fi 显示最小主应力?2 Plot sdif fi 显示主应力差(?1- ?2) Plot str 显示主应力矢量场 Plot xdis fi 显示X方向位移 Plot ydis fi 显示Y方向位移 Plot disp 显示位移矢量场 Plot e_p--塑性应变的轮廓线（只限于应变软化材料）； ff—以矢量形式画出固定结点的反力；pp—孔隙压力的等值线；wate—水位线； FLAC 使用步骤 FLAC 使用步骤 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，以节省时间。 2.FLAC 3D 程序的编写步骤： 1 Config　________ 2 Grid　　________ 3 Model 　________ 4求起始的应力平衡 (1)建立x, y坐标与网格的关系，建议使用Gen指示： Gen x1,y1 x2,y2 ,x3,y3 x4,y4 i=i0,i1 j=j0,j1 详细指令参见使用手册，FLAC程式可自动产生x, y 坐标与网格的关系，但由于产生的网格