flac学习吐血经验.doc

第一天 看版主的初学者建议！ 初学者将 FLAC/FLAC3D 版所有帖子都浏览一遍；学有所得后，再浏览一遍，会发现又有新的收获 Flac3D 中为岩土工程问题的求解开发了特有的本构模型，总共包含了 10 种材料模型： 1. 开挖模型 ] 2. 3 个弹性模型（各向同性，横观各向同性和正交各向同性弹性模型） 3. 6 个塑性模型（Drucker-Prager 模型、Morh-Coulomb 模型、应变硬化/ 软化模型、遍布节理模型、双线性应变硬化/软化遍布节理模型和修正的 cam 粘土模型） 因为我是学采矿专业的，我就先学习一下开挖模型吧！ 第二天 学习建模 FLAC3D 建模的一般步骤如下： 一、建模 new ge zo b po x y z & p1 x y z &# Z/ E7 g/ L( F) B% A p2 x y z &7 F. m$ ^ U- g5 h' S7 u p3 x y z &1 i* d1 K% C; c9 y6 R; a size x y z plotShowadd surface color add axes color 这样一个立体模型就建立出来了！ 2、定义材料参数mo moh/null/...pro bu=~ sh=~ fric=~ c=~ ten=~ 其中 bulk 是指体积模量也就是 K=E/3(1-2v),E 是弹模 v 是泊松比，当 v 接近 0.5 时，不能盲目使用；shear 是指剪切模量也就是 G=E/2(1+v） ； cohesion 是指内聚力； tension 是指抗拉强度。

3、加载及边界 set g 0 0 -10 in de=~fix x/y/z ra x/y/z ~ ~ 其中要有一个面为自由面四、求解 7 R( ^: m9 s% n C. fh n=10h unb h gp x/y/zdisp x y z set mec fo 30so 此为显式求解过程 五、结果分析监控变量：print hplot history 1history 2位移等值线图： plot cr ~add con disp/... add axes color show 应力等值线图： Clearadd bcon szz/... add axes剖面上的等值线图： cr ~ set p dip=~ dd=~ o=x,y,z add bo beadd bcon szz/... Padd axesshow 有这些命令可以查看你想分析的 FLAC 图 七、开挖pro c=~ ten=~ mo null ra x=~,~ y=~,~ z=~,~set large/small1in xdisp=0 ydisp=0 zdisp=0 系统中的所有节点位移全部清零8、求解与后处理step~ plot cr dispcont plot copy ~ dispcont settingsplot add con disp p be shade on plot add axes plot show开挖模型终于学习完了，经过学习我可以建立一个简单的开挖模型并计算、分析了。

下一步的工作：我想建立一个曲面的立体模型来反映矿山地形地貌研究中...... 经过学习，我可以建立一个分层的地质模型，怎样把地表设置成曲面呢？形成矿体上方的地形地貌! 一个分层的地质模型： 命令流如下：new gen zone b size 40 10 1 p0 0 0 0 p1 400 0 0 p2 0 200 0 p3 0 0 1.21 group 泥岩 1 gen zone b size 40 10 2 p0 0 0 1.21 p1 400 0 1.21 p2 0 200 1.21 p3 0 0 7.07 group 灰岩 1 gen zone b size 40 10 3 p0 0 0 7.07 p1 400 0 7.07 p2 0 200 7.07 p3 0 0 13.47 group 粉砂岩 1 gen zone b size 40 10 1 p0 0 0 13.47 p1 400 0 13.47 p2 0 200 13.47 p3 0 0 14.96 group 灰岩 2 gen zone b size 40 10 4 p0 0 0 14.96 p1 400 0 14.96 p2 0 200 14.96 p3 0 0 31.4 group 粉砂岩 2gen zone b size 40 10 1 p0 0 0 31.4 p1 400 0 31.4 p2 0 200 31.4 p3 0 0 32.46 group 泥岩 2 gen zone b size 40 10 3 p0 0 0 32.46 p1 400 0 32.46 p2 0 200 32.46 p3 0 0 39.73 group 粉砂岩 3 gen zone b size 40 10 6 p0 0 0 39.73 p1 400 0 39.73 p2 0 200 39.73 p3 0 0 51.93 group 细砂岩 1 gen zone b size 40 10 7 p0 0 0 51.93 p1 400 0 51.93 p2 0 200 51.93 p3 0 0 65 group 粉砂岩gen zone b size 40 10 3 p0 0 0 65 p1 400 0 65 p2 0 200 65 p3 0 0 67.93 group 10 煤gen zone b size 40 10 4 p0 0 0 67.93 p1 400 0 67.93 p2 0 200 67.93 p3 0 0 71.97 group 泥岩gen zone b size 40 10 1 p0 0 0 71.97 p1 400 0 71.97 p2 0 200 71.97 p3 0 0 72.35 group 天然焦 gen zone b size 40 10 2 p0 0 0 72.35 p1 400 0 72.35 p2 0 200 72.35 p3 0 0 75.15 group 岩浆岩gen zone b size 40 10 2 p0 0 0 75.15 p1 400 0 75.15 p2 0 200 75.15 p3 0 0 78.85 group 细砂岩 2 gen zone b size 40 10 25 p0 0 0 78.85 p1 400 0 78.85 p2 0 200 78.85 p3 0 0 132.6 group 泥岩 4. gen zone b size 40 10 4 p0 0 0 132.6 p1 400 0 132.6 p2 0 200 132.6 p3 0 0 136.79 group 细砂岩 3 gen zone b size 40 10 3 p0 0 0 136.79 p1 400 0 136.79 p2 0 200 136.79 p3 0 0 140.24 group 泥岩 5, gen zone b size 40 10 3 p0 0 0 140.24 p1 400 0 140.24 p2 0 200 140.24 p3 0 0 146.29 group 粉砂岩 5 gen zone b size 40 10 4 p0 0 0 146.29 p1 400 0 146.29 p2 0 200 146.29 p3 0 0 150 group 泥岩 6 plotShowadd surface greenadd axes red1 FLAC 建模方法1.1 建模（1）设计计算模型的尺寸（2）规划计算网格数目和分布（3）安排工程对象（开挖、支护等） （4）给出材料的力学参数 （5）确定边界条件（6）计算模拟 1.2 网格生成：Grid i,j 例如：grid 30,20 1.3 网格规划： Gen x1,y1 x2,y2 x3,y3 x4,y4 例如：Gen 0,0 0,10 10,20 20,0 1.4 分区规划网格。

例如：Gen xI1,yI1 xI2,yI2 xI3,yI3 xI4,yI4 i=1,10 j=1,21 （I 区）Gen xII1,yII1 xII2,yII2 xII3,yII3 xII4,yII4 i=10,20 j=1,21 (II 区） 1.5 特殊形状的网格（1）圆形 gen circle xc,yc rad（2）弧线 gen arc xc,yc xb,yb theta（3）直线 gen line x1,y1 x2,y2 （4）任意形状 tab 1 x1,y1, x2,y2, ¼,xn,yn, x1,y1 gen tab 11.6 赋给单元材料性质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,20mod m （弹塑性 Mohr-Coulumb 准则）prop d 1800e-6 bu 12.5 sh 5.77 c 0 fri 20 ten 0.015 reg i,j 1.7 赋给模型边界条件（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 （3）赋单元内应力 （单元） ini sxx -10 i=1,20 j=1,20ini syy -5 var 0 4 i=1,21 j=1,21. 1.8 计算 Set grav 9.81 Set largeStep 1000 Save test.sav1.9 结果显示 Plot grid 显示网格Plot bo 显示边界 Plot plas 显示塑性区Plot sig1 fi 显示最大主应力 s1Plot sig2 fi 显示最小主应力 s2Plot sdif fi 显示主应力差(s1- s2)) Plot str 显示主应力矢量场Plot xdis fi 显示 X 方向位移Plot ydis fi 显示 Y 方向位移 Plot disp 显示位移矢量场1.10 保存与调用结果命令Call test.txt (或 ca test.dat) 调用数据 Save test.sav 保存结果 New 重新开始Rest test.sav 调用结果 Quit 退出程序 2. FLAC 运用技巧 2.1 模型尺寸 2.2 模拟开挖 Mod nu i=6,15 j=5,12 (或 region i,j) 2.3 模拟锚杆支护（端锚）struct cable begin grid i1,j1 end grid i2,j2 seg n prop 1 stru prop 1 e 2e5 yield 0.5 a 0.235e-3 sbond 0 kbond 0.0001 stru prop 1 sfri 30 peri 0.2723 den 7.5e-。