搜索
第一天看版主的初学者建议！初学者将FLAC/FLAC3D版所有帖子都浏览一遍；学有所得后，再浏览一遍，会发现又有新的收获。Flac3D中为岩土工程问题的求解开发了特有的本构模型，总共包含了10种材料模型：1.开挖模型]2.3个弹性模型（各向同性，横观各向同性和正交各向同性弹性模型）3.6个塑性模型（Drucker-Prager模型、Morh-Coulomb模型、应变硬化/软化模型、遍布节理模型、双线性应变硬化/软化遍布节理模型和修正的cam粘土模型）。因为我是学采矿专业的，我就先学习一下开挖模型吧！第二天学习建模FLAC3D建模的一般步骤如下：一、建模newgezobpoxyz&p1xyz&#Z/E7g/L(F)B%Ap2xyz&7F.m$^U-g5h'S7up3xyz&1i*d1K%C;c9y6R;asizexyzplotShowaddsurfacecoloraddaxescolor这样一个立体模型就建立出来了！二、定义材料参数momoh/null/...probu=~sh=~fric=~c=~ten=~其中bulk是指体积模量也就是K=E/3(1-2v),E是弹模v是泊松比，当v接近0.5时，不能盲目使用；shear是指剪切模量也就是G=E/2(1+v）；cohesion是指内聚力；tension是指抗拉强度。三、加载及边界setg00-10inde=~fixx/y/zrax/y/z~~其中要有一个面为自由面。四、求解7R(^:m9s%nC.fhn=10hunbhgpx/y/zdispxyzsetmecfo30so此为显式求解过程。五、结果分析监控变量：printhplothistory1history2位移等值线图：plotcr~addcondisp/...addaxescolorshow应力等值线图：Clearaddbconszz/...addaxes剖面上的等值线图：cr~setpdip=~dd=~o=x,y,zaddbobeaddbconszz/...Paddaxesshow有这些命令可以查看你想分析的FLAC图七、开挖proc=~ten=~monullrax=~,~y=~,~z=~,~setlarge/small1inxdisp=0ydisp=0zdisp=0系统中的所有节点位移全部清零。八、求解与后处理step~plotcrdispcontplotcopy~dispcontsettingsplotaddcondisppbeshadeonplotaddaxesplotshow开挖模型终于学习完了，经过学习我可以建立一个简单的开挖模型并计算、分析了。下一步的工作：我想建立一个曲面的立体模型来反映矿山地形地貌。研究中......经过学习，我可以建立一个分层的地质模型，怎样把地表设置成曲面呢？形成矿体上方的地形地貌!一个分层的地质模型：命令流如下：newgenzonebsize40101p0000p140000p202000p3001.21group泥岩1genzonebsize40102p0001.21p140001.21p202001.21p3007.07group灰岩1genzonebsize40103p0007.07p140007.07p202007.07p30013.47group粉砂岩1genzonebsize40101p00013.47p1400013.47p2020013.47p30014.96group灰岩2genzonebsize40104p00014.96p1400014.96p2020014.96p30031.4group粉砂岩2genzonebsize40101p00031.4p1400031.4p2020031.4p30032.46group泥岩2genzonebsize40103p00032.46p1400032.46p2020032.46p30039.73group粉砂岩3genzonebsize40106p00039.73p1400039.73p2020039.73p30051.93group细砂岩1genzonebsize40107p00051.93p1400051.93p2020051.93p30065group粉砂岩genzonebsize40103p00065p1400065p2020065p30067.93group10煤genzonebsize40104p00067.93p1400067.93p2020067.93p30071.97group泥岩genzonebsize40101p00071.97p1400071.97p2020071.97p30072.35group天然焦genzonebsize40102p00072.35p1400072.35p2020072.35p30075.15group岩浆岩genzonebsize40102p00075.15p1400075.15p2020075.15p30078.85group细砂岩2genzonebsize401025p00078.85p1400078.85p2020078.85p300132.6group泥岩4.genzonebsize40104p000132.6p14000132.6p20200132.6p300136.79group细砂岩3genzonebsize40103p000136.79p14000136.79p20200136.79p300140.24group泥岩5,genzonebsize40103p000140.24p14000140.24p20200140.24p300146.29group粉砂岩5genzonebsize40104p000146.29p14000146.29p20200146.29p300150group泥岩6plotShowaddsurfacegreenaddaxesred1FLAC建模方法1.1建模（1）设计计算模型的尺寸（2）规划计算网格数目和分布（3）安排工程对象（开挖、支护等）（4）给出材料的力学参数（5）确定边界条件（6）计算模拟1.2网格生成：Gridi,j例如：grid30,201.3网格规划：Genx1,y1x2,y2x3,y3x4,y4例如：Gen0,00,1010,2020,01.4分区规划网格。例如：GenxI1,yI1xI2,yI2xI3,yI3xI4,yI4i=1,10j=1,21（I区）GenxII1,yII1xII2,yII2xII3,yII3xII4,yII4i=10,20j=1,21(II区）1.5特殊形状的网格（1）圆形gencirclexc,ycrad（2）弧线genarcxc,ycxb,ybtheta（3）直线genlinex1,y1x2,y2（4）任意形状tab1x1,y1,x2,y2,¼,xn,yn,x1,y1gentab11.6赋给单元材料性质mode(弹性)propd1800e-6bu12.5sh5.77i=1,20j=1,10propd2400e-6bu1250sh577i=1,20j=11,20modm（弹塑性Mohr-Coulumb准则）propd1800e-6bu12.5sh5.77c0fri20ten0.015regi,j1.7赋给模型边界条件（1）固定边界（结点）Fixxi=1,j=1,21Fixyi=1,21j=1（2）施加边界力(结点)applyyf-10i=1,21j=21或applysyy-10i=1,21j=21或applyxf-5i=21,j=1,21或applysxx-5i=21,j=1,21（3）赋单元内应力（单元）inisxx-10i=1,20j=1,20inisyy-5var04i=1,21j=1,21.1.8计算Setgrav9.81SetlargeStep1000Savetest.sav1.9结果显示Plotgrid显示网格Plotbo显示边界Plotplas显示塑性区Plotsig1fi显示最大主应力s1Plotsig2fi显示最小主应力s2Plotsdiffi显示主应力差(s1-s2))Plotstr显示主应力矢量场Plotxdisfi显示X方向位移Plotydisfi显示Y方向位移Plotdisp显示位移矢量场1.10保存与调用结果命令Calltest.txt(或catest.dat)调用数据Savetest.sav保存结果New重新开始Resttest.sav调用结果Quit退出程序2.FLAC运用技巧2.1模型尺寸2.2模拟开挖Modnui=6,15j=5,12(或regioni,j)2.3模拟锚杆支护（端锚）structcablebegingridi1,j1endgridi2,j2segnprop1struprop1e2e5yield0.5a0.235e-3sbond0kbond0.0001struprop1sfri30peri0.2723den7.5e-32.4模拟锚杆支护（全长锚固）structcablebegingridi,jendx,ysegnprop2struprop1e2e5yield0.5a0.235e-3sbond0.42kbond5.37struprop1sfri30perimeter0.2723den7.5e-32.5模拟锚杆支护（预应力锚固）structcablebegingridi1,j1endi2,j2segnten10prop1structcablebeginnoden1endx,ysegnprop2struprop1e2e5yield0.5a0.235e-3sbond0.kbond0.1struprop1sfri0perimeter0.2723den7.5e-3struprop2e2e5yield0.5a0.235e-3sb0.42kb5.37struprop2sfri30perimeter0.2723den7.5e-32.6模拟断层Modnuj=37int1asfromi1,j1to,i1,j2bsfromi2,j1toi2,j2int1kn1000ks400fri25c0.01t0.00012.7其它技巧（1）网格优化（2）先弹性、后塑性（3）分步开挖（4）动力学问题模拟.（5）固流耦合问题模拟（6）合理解释结果（7）多种形式输出结果（8）报告格式第三、四天学习彭文斌的《FLAC3D实用教程》第5天运用FLAC做一个小项目第六天：分享一些FLAC3D命令的FAQ资料！1.如何调整视图的大小、角度？,答：综合使用x,y,z,m,Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键2.最先需要掌握的命令有哪些？答：需要掌握gen,ini,app,plo,solve等建模、初始条件、边界条件、后处理和求解的命令。3.怎样看模型的样子？答：ploblogro可以看到不同的group的颜色分布4.怎样看模型的边界情况？答：plogpfixredsk5.怎样看模型的体力分布？答：plofapredsk6.怎样看模型的云图？答：位移：plocondis(xdis,ydis,zdis)应力：ploconsz(sy,sx,sxy,syz,sxz)7.怎样看模型的矢量图？答：plodis(xdis,ydis,zdis8.怎样看模型有多少单元、节点？答：ploinfo9.怎样输出模型的后处理图？答：File/Printtype/Jpgfile，然后选择File/Print，将保存格式选择为jpg文件10.怎样调用一个文件？答：File/call或者call命令11.如何施加面力？!答：appnstress12.如何进行边界约束？答：fixxran（约束的是速度，在初始情况下约束等效于位移约束）13.如何知道每个单元的ID？答：用鼠标双击单元的表面，可以知道单元的ID和坐标14.如何进行切片？答：plosetplaneori(点坐标)norm(法向矢量)ploconszplane(显示z方向应力的切片）15.如何保存计算结果？答：save+文件名