Flac问题讲解备课讲稿

1、Good is good, but better carries it.精益求精，善益求善。Flac问题讲解-分差法岩土工程结构的数值解是建立在满足基本方程（平衡方程、几何方程、本构方程）和边界条件下推导的。由于基本方程和边界条件多以微分方程的形式出现，因此，将基本方程近假发改用差分方程（代数方程）表示，把求解微分方程的问题改换成求解代数方程的问题，这就是所谓的分差法。差分法由来已久，但差分法需要求解高阶代数方程组，只有在计算机的出现，才使该法得以实施和发展。FLAC3D（FastLagrangianAnalysisofContinua）由美国Itasca公司开发的。目前，FLAC有二维和三维

2、计算程序两个版本，二维计算程序V3.0以前的为DOS版本，V2.5版本仅仅能够使用计算机的基本内存64K），所以，程序求解的最大结点数仅限于2000个以内。1995年，FLAC2D已升级为V3.3的版本，其程序能够使用护展内存。因此，大大发护展了计算规模。FLAC3D是一个三维有限差分程序，目前已发展到V2。1版本。FLAC3D的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入各种命令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。因此，采用FLAC程序进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排列，形成可以完成一定计算任务的命令文件。F

3、LAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。调整三维网格中的多面体单元来拟合实际的结构。单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发变形和移动（大变形模式）。FLAC3D采用的显式拉格朗日算法和混合-离散分区技术能够非常准确发模拟材料的塑性破坏和流动。由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。FLAC3D采用ANSIC+语言编写的。FLAC3D有以下几个优点：1对模拟塑性破坏和塑性流动采用的是“混合离散法“。这种方法比有限元法中通常采用的“离散集成法“

4、更为准确、合理。2即使模拟的系统是静态的，仍采用了动态运动方程，这使得FLAC3D在模拟物理上的不稳定过程不存在数值上的障碍。3采用了一个“显式解“方案。因此，显式解方案对非线性的应力-应变关系的求解所花费的时间，几互与线性本构关系相同，而隐式求解方案将会花费较长的时间求解非线性问题。面且，它没有必要存储刚度矩阵，这就意味着；采用中等容量的内存可以求解多单元结构；模拟大变形问题几互并不比小变形问题多消耗更多的计算时间，因为没有任何刚度矩阵要被修改。当然，它也存在以下几个不足之处：1对于线性问题的求解，FLAC3D比有限元程序运行得要慢；因此，当进行大变形非线性问题或模拟实际可能出现不稳定问题时

5、，FLAC3D是最有效的工具。2用FLAC3D求解时间取决于最长的自然周期和最短的自然周期之比。但某些问题对模型是无效的。接上面的介绍：Flac3D中为岩土工程问题的求解开发了特有的本构模型，总共包含了10种材料模型：1.开挖模型null2.3个弹性模型（各向同性，横观各向同性和正交各向同性弹性模型）36个塑性模型（Drucker-Prager模型、Morh-Coulomb模型、应变硬化/软化模型、遍布节理模型、双线性应变硬化/软化遍布节理模型和修正的cam粘土模型）。Flac3D网格中的每个区域可以给以不同的材料模型，并且还允许指定材料参数的统计分布和变化梯度。还包含了节理单元，也称为界面单

6、元，能够模拟两种或多种材料界面不同材料性质的间断特性。节理允许发生滑动或分离，因此可以用来模拟岩体中的断层、节理或摩擦边界。FLAC3D中的网格生成器gen，通过匹配、连接由网格生成器生成局部网格，能够方便地生成所需要的三维结构网格。还可以自动产生交叉结构网格（比如说相交的巷道），三维网格由整体坐标系x,y,z系统所确定，不同于FLAC程序是由行列方式确定。这就提供了比较灵活的产生和定义三维空间参数。FLAC3D做计算分析的一般步骤：与大多数程序采用数据输入方式不同，FLAC采用的是命令驱动方式。命令字控制着程序的运行。在必要时，尤其是绘图，还可以启动FLAc用户交互式图形界面。为了建立FLA

7、C计算模型，必须进行以下三个方面的工作：1.有限差分网格2.本构特性与材料性质3.边界条件与初始条件完成上述工作后，可以获得模型的初始平衡状态，也就是模拟开挖前的原岩应力状态。然后，进行工程开挖或改变边界条件来进行工程的响应分析，类似于FLAC的显式有限差分程序的问题求解。与传统的隐式求解程序不同，FLAC采用一种显式的时间步来求解代数方程。进行一系列计算步后达到问题的解。在FLAC中，达到问题所需的计算步能够通过程序或用户加以控制，但是，用户必须确定计算步是否已经达到问题的最终的解。最后，进行结果的分析与总结，这就是用户自己的事情了，呵呵，在手册中的图中有显示shell单元但我自己做的时候不

8、知用何命令才能在网格图中将的shell，beam等单元显示出来，高手指点！plotaddselgeom另外，还可在上面命令的后面加：cid=onnode=off.再问一下，shell是否只可以加在网格的表面？怎么在网格内部加不上？可以加上了，但是不知道厚度t如何表现。thickness啊gensurface有什么用？生成表面显示模型中用到setdyn=on,grav00-10,hist_rep=1是在动力分析中的问题，请问最后一个变量是什么含义，好像只能是整数，查了Fish说明也没有看到解释，请赐教！hist_repnHistoriesaresampledeveryntimesteps.The

9、defaultisn=10.ThesynonymHISTORYnstepmayalsobeused.每10个步时记录一次请教:(1)谁用过GENERATEsurfacexarc命令生成过圆弧面，能否贴出来让小弟学一下。我曾按照manual作过一些尝试但没有成功。(2)根据说明我曾生成过圆形隧道模型，但如何在隧道内壁上加正应力，请各位高手指点。加正应力应该是用nstress命令我没有做过曲面的例子，只用过规则的圆柱面请问各位大侠，groutstiffness是什么意思？1.gr_kgroutstiffnessperunitlength,kgF/L22.Ft=kg*ut3.Thegroutannu

10、lusisassumedtobehaveasanelastic-perfectlyplasticsolid.Asaresultofrelativesheardisplacement,ut,betweenthetendonsurfaceandtheboreholesurface,theshearforce,Ft,mobilizedperlengthofcableisrelatedtothegroutstiffness,kgi.e.,4.可见k_g是groutannulus与zone（soil）之间的剪切刚度。（单位剪切位移引起的应力）5.详细参考cable一节，有图形表示，很清楚。写了几个命令流

11、文件，但只要所分析问题的单元数一多，计算就很难收敛，总是在e-3或e-4之间跳来跳去，都迭代几万步了请教大侠，怎么加快收敛速度？一般一万五千单元左右要算多久？flac中的求解控制通常有1）setmechforce50solve2）solveratio1e-73）step30000等几种，不知这几种怎么选择？每种有什么特点？加快收敛速度，要让单元的尺寸相差不大，另外，ratio恒定，说明产生了塑性流动或者破坏了，看看约束条件对不对！材料模型对不对，材料参数是否太离谱？flac中的求解控制通常有1）setmechforce50solve2）solveratio1e-73）step30000第一种是

12、控制不平衡力，是给定具体的值，必须迭代到这个值才结束运算。这是绝对的值，个人认为只有针对具体问题才使用；第二种要重先设定，程序默认1e5，7也太小了，迭代的步数更多，更精确，一般情况下，个人认为5就够了；第三种是给定迭代步数，不至于30000吧，一般3000吧。但这种设定不一定能达到1e5，所以要在有把握的情况下使用。至于在循环中，solve与stepn到底有何差别，本人也不甚清楚。有时某一步可能需要7000步才能达到5次方，而其他的循环只需要1000步就达到了，这种情况下究竟该如何处理？请高手指点！solve是通过控制不平衡力或比例的大小决定计算步step，cyc等命令则是由作者控制计算迭代

13、步，在无法用step估计计算是否收敛情况下，多用solve，step自身的优势可以应用于逐步加载等方面。加大bulk，shear的量纲试一试，还有就是检查一下边界条件！;产生一个体，这个体里面包含隧道genzoneradbrick&p0(0,0,0)p1(10,0,0)p2(0,10,0)p3(0,0,10)&size3,5,5,7&ratio1,1,1,1.5&dim142fill;设置模型服从摩尔库仑屈服准则modmoh;定义bulk模量，剪切模量，摩擦角，粘聚力，和tenspropbulk1e8shear3e8fric35coh1e3tens1e3;施加约束fixxrangex-0.10

14、.1fixzrangez-0.10.1fixyrangey9.910.1;施加重力setgrav0100;设置密度inidens1000;设置初始应力，这个有什么用呢？inisyy0.0grad0-100000inisxx0.0grad0-50000iniszz0.0grad0-50000applysxx0.0grad0-50000rangex9.910.1applyszz0.0grad0-50000rangez9.910.1;挖坑的命令，使用fish语言defdigHolecommandrangenametrenchx01y04z02;挖了modelnulrangetrenchsetlarg

15、ehistgpxdisp100endcommandend;表面蓝色显示plotsurfblue;停顿pause;开挖digHole;设置大变形setlarge;记录位移变化histgpxdisp1,0,0;.andthencalculatetheresponse;设置计算步为1000step1000;我不知道上述解释是否明白，是否对头？cjiao2000wrote:;产生一个体，这个体里面包含隧道genzoneradbrick&p0(0,0,0)p1(10,0,0)p2(0,10,0)p3(0,0,10)&size3,5,5,7&ratio1,1,1,1.5&dim142fill;设置模型服从

16、摩尔库仑屈服准则modmoh;定义bulk模量，剪切模量，摩擦角，粘聚力，和tenspropbulk1e8shear3e8fric35coh1e3tens1e3;施加约束fixxrangex-0.10.1fixzrangez-0.10.1fixyrangey9.910.1;施加重力setgrav0100;设置密度inidens1000;设置初始应力，这个有什么用呢？inisyy0.0grad0-100000inisxx0.0grad0-50000iniszz0.0grad0-50000别的软件很少可以加地应力的，怎么说没用呢？这可是Flac的专长呀，象ansys也没的办法的。applysxx0

17、.0grad0-50000rangex9.910.1applyszz0.0grad0-50000rangez9.910.1;挖坑的命令，使用fish语言defdigHolecommandrangenametrenchx01y04z02;挖了modelnulrangetrenchsetlargehistgpxdisp100endcommandend;表面蓝色显示plotsurfblue;休息pause;开挖digHole;设置大变形setlarge;这个好像是施加位移我认为是记录坐标是（1，0，0，）点的位移历史记录吧，不是施加位移吧。权当抛砖引玉了，各位大侠不要见笑histgpxdisp1,0

18、,0;.andthencalculatetheresponse;设置计算步为1000step1000;我不知道上述解释是否明白，是否对头？首先确定位移是z向，还是x向histgpxdisp1,0,0有时候会出现位移都为0的情况,怎么回事?明明那是有位移的啊?是啊，还有向上得哦！;设置初始应力，这个有什么用呢？inisyy0.0grad0-100000inisxx0.0grad0-50000iniszz0.0grad0-50000别的软件很少可以加地应力的，怎么说没用呢？这可是Flac的专长呀，象ansys也没的办法的。设置初始应力做什么哦？一般情况下围岩地应力包括自重应力、构造应力、采掘应力。

19、你所描速的情况应是在挖掘深度不大的情况下，只考虑土体的自重应力对开挖的影响。在古典压力理论中，认为地层的水平压力的侧压系数为u/1-u，这里取0.5。apply描述的模型的边界条件。histgpxdisp1,0,0是用于记录点（1，0，0）在开挖过程中的位移情况，模拟的目的就是要研究目标位置的位移或者应力变化情况。setlarge设置大变形的话,应该是按大变形理论吧!最近做了算例，发现将cable单元的pretension设置为0或者100没有什么区别，大家看一下是怎么回事啊？;generatethemodletitleVerticalAnchoredGravityWall;setmechar

20、atio1e-5;setlarge;Zonegenzoneradtunnelp01.5,0,-20p11.5,20,-20p20,0,-20p31.5,0,0&p40,20,-20p50,0,0p61.5,20,0p70,20,0&p81.56-20p91.50-10.5p1006-20p1100-10.5p121.56-10.5p1306-10.5&size12,3,19,20ratio1,1,1,1dim9.559.55fill;Material;Forbalancemodelmohrpropdensity=1.80e3bulk=12.5e6shear=6.5e6coh10.0e3tens

21、50e10&friction=10.0dilation=0.0rangez-20.0-18.0propdensity=2.00e3bulk=20.0e6shear=12.0e6coh15.0e3tens50e10&friction=25.0dilation=5.0rangez-18.0-10.0propdensity=2.00e3bulk=44.4e6shear=33.3e6coh0.0e3tens50e10&friction=30.0dilation=10.0rangez-10.00.0;Boundaryconditonfixxrangex-0.10.1fixxrangex0.41.5fix

22、yrangey-0.1.1fixyrangey19.920.1fixzrangez-0.10.1;Setgravitionalaccsetgrav0,0,10;Initialstressiniszz-4.0e5grad00-2.0e4inisxx-2.0e5grad00-1.00e4inisyy-2.0e5grad00-1.00e4;InitialState,Solveforbalancesolveinistate0inixdisp0.0ydisp0.0zdisp0.0;Reinstallthesoilmodelmohrpropdensity=1.80e3bulk=12.5e6shear=6.

23、5e6coh10.0e3tens58e3&friction=10.0dilation=0.0rangez-20.0-18.0propdensity=2.00e3bulk=20.0e6shear=12.0e6coh15.0e3tens32.6e3&friction=25.0dilation=5.0rangez-18.0-10.0propdensity=2.00e3bulk=44.4e6shear=33.3e6coh0.0e3tens5e3&friction=30.0dilation=10.0rangez-10.00.0;Setthesurchargeapplyszz10e3rangey620z-

24、20.1-19.9solve;Begintoexcavate;Step1:1.5moffmodnullrangey06.0z-18.5-20.0sellinerrangey5.96.1z-18.5-20.0sellinerPROPiso=(25e9,0.15)thick=0.1;concretesellinerPROPcs_nk=8e8cs_sk=8e8&cs_ncut=4e6cs_scoh=4e4cs_scohres=2e6cs_sfric=20.0solvesavestep1.sav;Step2:setthesoilnail,and1.5moffselcableid=1begin=(0.7

25、5,6.0,-18.5)end=(0.75,15.5,-17.5)nseg=12pretention50e3selcablepropxcarea=5.1e-3emod=200e9ytens=210e6&gr_k=6e6gr_coh=40e3selcablepropgr_per=0.314gr_fric=25.0modnullrangey06.0z-17.0-18.5sellinerrangey5.96.1z-17.0-18.5sellinerPROPiso=(25e9,0.15)thick=0.1;concretesellinerPROPcs_nk=8e8cs_sk=8e8&cs_ncut=4

26、e6cs_scoh=4e4cs_scohres=2e6cs_sfric=20.0solvesavestep2.sav;Step3:setthesoilnail,and1.5moffselcableid=2begin=(0.75,6.0,-17.0)end=(0.75,15.5,-16.0)nseg=12pretention50e3selcablepropxcarea=5.1e-3emod=200e9ytens=210e6&gr_k=6e6gr_coh=40e3selcablepropgr_per=0.314gr_fric=25.0modnullrangey06.0z-15.5-17.0sell

27、inerrangey5.96.1z-15.5-17.0sellinerPROPiso=(25e9,0.15)thick=0.1;concretesellinerPROPcs_nk=8e8cs_sk=8e8&cs_ncut=4e6cs_scoh=4e4cs_scohres=2e6cs_sfric=20.0solvesavestep3.sav;Step4:setthesoilnail,and1.5moffselcableid=3begin=(0.75,6.0,-15.5)end=(0.75,15.5,-14.5)nseg=12pretention50e3selcablepropxcarea=5.1

28、e-3emod=200e9ytens=210e6&gr_k=6e6gr_coh=40e3selcablepropgr_per=0.314gr_fric=25.0modnullrangey06.0z-14.0-15.5sellinerrangey5.96.1z-14.0-15.5sellinerPROPiso=(25e9,0.15)thick=0.1;concretesellinerPROPcs_nk=8e8cs_sk=8e8&cs_ncut=4e6cs_scoh=4e4cs_scohres=2e6cs_sfric=20.0solvesavestep4.sav;Step5:setthesoiln

29、ail,and1.5moffselcableid=3begin=(0.75,6.0,-14.0)end=(0.75,15.5,-13.0)nseg=12pretention50e3selcablepropxcarea=5.1e-3emod=200e9ytens=210e6&gr_k=6e6gr_coh=40e3selcablepropgr_per=0.314gr_fric=25.0modnullrangey06.0z-12.5-14.0sellinerrangey5.96.1z-12.5-14.0sellinerPROPiso=(25e9,0.15)thick=0.1;concretesell

30、inerPROPcs_nk=8e8cs_sk=8e8&cs_ncut=4e6cs_scoh=4e4cs_scohres=2e6cs_sfric=20.0solvesavestep5.sav;Step6:setthesoilnail,and1.5moffselcableid=3begin=(0.75,6.0,-12.5)end=(0.75,15.5,-11.5)nseg=12pretention50e3selcablepropxcarea=5.1e-3emod=200e9ytens=210e6&gr_k=6e6gr_coh=40e3selcablepropgr_per=0.314gr_fric=

31、25.0modnullrangey06.0z-11.5-12.5sellinerrangey5.96.1z-11.5-12.5sellinerPROPiso=(25e9,0.15)thick=0.1;concretesellinerPROPcs_nk=8e8cs_sk=8e8&cs_ncut=4e6cs_scoh=4e4cs_scohres=2e6cs_sfric=20.0solvesavestep6.sav;Step7:setthesoilnail,and1.5moffselcableid=3begin=(0.75,6.0,-11.0)end=(0.75,15.5,-10.0)nseg=12

32、pretention50e3selcablepropxcarea=5.1e-3emod=200e9ytens=210e6&gr_k=6e6gr_coh=40e3selcablepropgr_per=0.314gr_fric=25.0modnullrangey06.0z-10.5-11.5sellinerrangey5.96.1z-10.5-11.5sellinerPROPiso=(25e9,0.15)thick=0.1;concretesellinerPROPcs_nk=8e8cs_sk=8e8&cs_ncut=4e6cs_scoh=4e4cs_scohres=2e6cs_sfric=20.0

33、solvesavestep7.savflac怎么作一个剖面！就是我建了一个模型，怎么看一个切面的视图。plsetplaneori000nor100plcondispplane在所建立的模型中有两个命令；就是printgpposition与printzonegp。不知道这两个命令中的gp是否有联系？按找手册中的解释，printgpposition是用来输出gp的位置坐标的；printzonegp是用来输出与zone相关的gp的id值的。这样我就有一个上面的疑问，两个命令中的gp都有id值，是不是不同命令中的gp的id值相同的话，他们的坐标就可以看作是一样的啊？是一个gp有且只有一个id谢谢dyn

34、amax兄。这下我就明白了，如果要画等值线的话，先输出gp的坐标来，将其另存为一个文件，再输出zone的应力值，将与zone对应的gp的id与前面输出的gp的id值对应起来，找到gp的坐标。用excel打开上面的文件进行编辑，最后用surfer画等值线。这样的办法太笨，不知到有没有更好的办法啊？欢迎前来讨论。你的思路很对！要简化操作，如果你对fish熟悉，可以直接按surfer的文件格式输出如果会用tecplot，可以参考前面我那篇【原创】用tecplot对flac3d计算结果进行后处理surfer我不熟，原则应该是一样的看了你的哪个原创，发现里面有这样几个语句不太理解，IO_READ=0IO

35、_WRITE=1IO_FISH=0IO_ASCII=1其中的0，1具体指的是什么啊？我发现一般zone中显示应力的gp的id往往不是全部模型的gp，那么如果想要用fish语句输出这些gp的坐标，该如和判断呢？不知道在fish语句中如何将zone中的gp的坐标给出。可以将gp的坐标和应力值存到一个数据文件，然后用surfer调用数据文件即可，不须在做额外的编辑比如这个命令：printzones输出了所有区的应力，可是要如果把它保存到指定写字板呢setlogonsetlogfile*.logprintzones依次执行以上命令便可以实现你的目的。其实简单的理解就是：dd：映射平面法向与y轴正向夹角

36、（顺时针）dip：映射平面法向与z轴负向夹角dip是倾角，就是plane和水平面的夹角dd是dipdirection，倾角的方向，注意不是走向就是你用来看倾角的那根线，在水平面上的投影，与y轴正向的交角从y轴正向开始算那位知道flac3D与ansys中单元不等比例剖分中网格长度与比例系数之间的计算公式。s=a(1-qn-1)/1-qa第一个网格长度q比例系数不知版主问的是不是这个公式等比公式在做边坡稳定性分析建模的时候,一定会遇到非自立边坡,此时生成模型的时候有几种方法.大家用什么方法?1.用不同的zonebrick拼装起来;2.用planeabove(该命令切出来的斜坡呈锯齿状,很难看,除非

37、网格很细,否则与实际偏离较大)3.使用fish,手册上有一个例子通过改变各个grid的坐标来实现隧道的球形开挖面,依此类推,一定可以修改grid的坐标实现斜坡.用fish或用zonebrick拼装生成非90度的基坑边坡，能给个例子吗？下面就是我编的，请大家帮我看看，到底哪里有问题，怎么老是有问题。（我观察到第一步开挖后，图上显示只挖掉了一排网格，我要挖的可是两排。在第二次开挖前，shellSEL竟然只布置上去2个！）我改了好多次都是类似结果，不知到底我哪弄错了!genzoneradbrickp16000p20600p30050size20201032&dim202010rat1111.05fi

38、llmodelmohrpropbulk1.25e7shear5.77e6fric20coh1.5e4rangez0.050.0;边界条件fixxrangex-0.010.01fixyrangey-0.010.01fixxyzrangez49.9950.01fixxyrangex59.9960.01fixxyrangey59.9960.01;初始条件inidensity1.9e3rangez0.0050.00setgravity0010iniszz0grad001.9e4rangez0.0050.00inisxx0grad007.6e3rangez0.0050.00inisyy0grad007.

39、6e3rangez0.0050.00setlargehistid1unbalplothist1setmechratio=1e-4solvesavej_ini.savinixdis=0ydis=0zdis=0grouplod1rangex2030y030z-0.010.01grouplod2rangex030y2030z-0.010.01grouplodrangegrouplod1anygrouplod2anyapplyszz10e3rangegrouplodmodelnullrangex020y020z0.02.0step2000savej_eca1.savselshellid=1rangex

40、19.920.1y020z0.02.0selshellid=2rangex020y19.920.1z0.02.0selnodefixxyrzrrangex-0.010.01selnodefixyxrzrrangey-0.010.01selshellid=1propiso=(26.6e9,0.25)thick=0.1selshellid=2propiso=(26.6e9,0.25)thick=0.1modelnullrangex020y020z2.04.0step2000savej_eca2.savselshellid=3rangex19.920.1y020z2.04.0selshellid=4

41、rangex020y19.920.1z2.04.0selnodefixxyrzrrangex-0.010.01selnodefixyxrzrrangey-0.010.01selshellid=3propiso=(26.6e9,0.25)thick=0.1selshellid=4propiso=(26.6e9,0.25)thick=0.1modelnullrangex020y020z4.06.0step2000savej_eca3.savselshellid=5rangex19.920.1y020z4.06.0selshellid=6rangex020y19.920.1z4.06.0selnod

42、efixxyrzrrangex-0.010.01selnodefixyxrzrrangey-0.010.01selshellid=5propiso=(26.6e9,0.25)thick=0.1selshellid=6propiso=(26.6e9,0.25)thick=0.1modelnullrangex020y020z6.08.0step2000savej_eca4.savselshellid=7rangex19.920.1y020z6.08.0selshellid=8rangex020y19.920.1z6.08.0selnodefixxyrzrrangex-0.010.01selnode

43、fixyxrzrrangey-0.010.01selshellid=7propiso=(26.6e9,0.25)thick=0.1selshellid=8propiso=(26.6e9,0.25)thick=0.1modelnullrangex020y020z8.010.0step2000savej_eca5.savselshellid=9rangex19.920.1y020z8.010.0selshellid=10rangex020y19.920.1z8.010.0selnodefixxyrzrrangex-0.010.01selnodefixyxrzrrangey-0.010.01sels

44、hellid=9propiso=(26.6e9,0.25)thick=0.1selshellid=10propiso=(26.6e9,0.25)thick=0.1step2000savej_finl.sav完整的命令流见附件genzoneradbrickp16000p20600p30050size20201032&dim202010rat1111.05fillmodelmohrpropbulk1.25e7shear5.77e6fric20coh1.5e4rangez0.050.0;边界条件fixxrangex-0.010.01fixyrangey-0.010.01fixxyzrangez49.

45、9950.01fixxyrangex59.9960.01fixxyrangey59.9960.01;初始条件inidensity1.9e3rangez0.0050.00setgravity0010iniszz0grad001.9e4rangez0.0050.00inisxx0grad007.6e3rangez0.0050.00inisyy0grad007.6e3rangez0.0050.00setlargehistid1unbalplothist1setmechratio=1e-4solvesavej_ini.savinixdis=0ydis=0zdis=0grouplod1rangex203

46、0y030z-0.010.01grouplod2rangex030y2030z-0.010.01grouplodrangegrouplod1anygrouplod2anyapplyszz10e3rangegrouplodmodelnullrangex020y020z0.02.0step2000savej_eca1.savmacrotudingemod2e11xcarea3.81e-4ytens2e8gr_coh2.5e4gr_k6.33e6gr_per0.346gr_fric25;设置cabledefplace_cable1xx1=1loopwhilexx120commandselcableb

47、eg20 xx11.5end25.908847xx12.541889nseg6selcablebegxx1201.5endxx125.9088472.541889nseg6end_commandxx1=xx1+1.5end_loopendplace_cable1selcableproptudingselshellid=1rangex19.920.1y020z0.02.0selshellid=2rangex020y19.920.1z0.02.0selnodefixxyrzrrangex-0.010.01selnodefixyxrzrrangey-0.010.01selshellid=1propi

48、so=(26.6e9,0.25)thick=0.1selshellid=2propiso=(26.6e9,0.25)thick=0.1modelnullrangex020y020z2.04.0step2000savej_eca2.savdefplace_cable2xx2=0.5loopwhilexx220commandselcablebeg20 xx23.0end25.908847xx24.041889nseg6selcablebegxx2203.0endxx225.9088474.041889nseg6end_commandxx2=xx2+1.5end_loopendplace_cable

49、2selcableproptudingselshellid=3rangex19.920.1y020z2.04.0selshellid=4rangex020y19.920.1z2.04.0selnodefixxyrzrrangex-0.010.01selnodefixyxrzrrangey-0.010.01selshellid=3propiso=(26.6e9,0.25)thick=0.1selshellid=4propiso=(26.6e9,0.25)thick=0.1modelnullrangex020y020z4.06.0step2000savej_eca3.savdefplace_cab

50、le3xx3=1loopwhilexx320commandselcablebeg20 xx34.5end25.908847xx35.541889nseg6selcablebegxx3204.5endxx325.9088475.541889nseg6end_commandxx3=xx3+1.5end_loopendplace_cable3selcableproptudingselshellid=5rangex19.920.1y020z4.06.0selshellid=6rangex020y19.920.1z4.06.0selnodefixxyrzrrangex-0.010.01selnodefi

51、xyxrzrrangey-0.010.01selshellid=5propiso=(26.6e9,0.25)thick=0.1selshellid=6propiso=(26.6e9,0.25)thick=0.1modelnullrangex020y020z6.08.0step2000savej_eca4.savdefplace_cable4xx4=0.5loopwhilexx420commandselcablebeg20 xx46.0end25.908847xx47.041889nseg6selcablebegxx4206.0endxx425.9088477.041889nseg6end_co

52、mmandxx4=xx4+1.5end_loopendplace_cable4selcableproptudingdefplace_cable5xx5=1loopwhilexx520commandselcablebeg20 xx57.5end25.908847xx58.541889nseg6selcablebegxx5207.5endxx525.9088478.541889nseg6end_commandxx5=xx5+1.5end_loopendplace_cable5selcableproptudingselshellid=7rangex19.920.1y020z6.08.0selshel

53、lid=8rangex020y19.920.1z6.08.0selnodefixxyrzrrangex-0.010.01selnodefixyxrzrrangey-0.010.01selshellid=7propiso=(26.6e9,0.25)thick=0.1selshellid=8propiso=(26.6e9,0.25)thick=0.1modelnullrangex020y020z8.010.0step2000savej_eca5.savdefplace_cable6xx6=1loopwhilexx620commandselcablebeg20 xx69.0end25.908847x

54、x610.041889nseg6selcablebegxx6209.0endxx625.90884710.041889nseg6end_commandxx6=xx6+1.5end_loopendplace_cable6selcableproptudingselshellid=9rangex19.920.1y020z8.010.0selshellid=10rangex020y19.920.1z8.010.0selnodefixxyrzrrangex-0.010.01selnodefixyxrzrrangey-0.010.01selshellid=9propiso=(26.6e9,0.25)thi

55、ck=0.1selshellid=10propiso=(26.6e9,0.25)thick=0.1step2000;水平位移变化plotcreatedisplaceplotsetcapsize20plotsetrot2000195plotaddcontxdispoutlonplotaddselcableforceplotaddaxesblackplotcreatetu_dispplotsetcapsize20plotsetrot1800180plotaddsketchplotadddispredplotaddselgeombluebluenode=offsavej_finl.savtensio

56、n=cohesion/tan(friction)好像在flac的theory里面有，就是莫尔园的外包线和x轴的交点横坐标。mohr模型有shearfailure和tensionfailure，如果把tension设置很小，有可能拉应力很小就达到塑性了，所以最好还是按前面那个关系设。回eplier：你的程序我没有运行过，但根据我个人的推断，或许是你的模型在达到初始平衡后已经破坏，故如果继续进行后续计算，就可能产生你所述的情况。所以，如将初始粘聚力值增大，使之在初始平衡时不会产生塑性破坏，在进行后续计算之前，在将粘聚力恢复为你所选择的参数，应该可以避免你的情况。你可以试试，也可以参考一下manua

57、l中相关的说明。如果有什么结果，还望在此告知。你的边界条件设置的合理吗？;边界条件fixxrangex-0.010.01fixyrangey-0.010.01fixxyzrangez49.9950.01fixxyrangex59.9960.01fixxyrangey59.9960.01;初始条件inidensity1.9e3rangez0.0050.00setgravity0010iniszz0grad001.9e4rangez0.0050.00inisxx0grad007.6e3rangez0.0050.00inisyy0grad007.6e3rangez0.0050.00你所建立的模型应该

58、是真实工程中的一半吧，你这个模型不复杂，你可以将整个模型建立起来，进行运算，如果单就半个模型来说，我想应该在x-z平面加上应力边界条件，当初建议作者用setsmall计算只是为了让其发现单元开挖一直有问题这个具体问题。实际应用中，大变形是flac的特色，而且如果分析对象是大变形问题，建议还是应用setlarge。至于单元开挖有误，说明作者程序肯定有问题。最近我正仔细研究flac计算原理，并阅读连续介质力学计算方法，收获很大，觉得flac的开发还需要很多工作要做。我正在整理，大约不长时间将会展开专门的详细论述。对的，基坑开挖是不可以用setsmall的，必须用setlarge，因为基坑在变形破坏

59、时一定是单元与原来的网格发生了比较大的位移，如果用小应变模式，那就是你限制了单元不让其破坏，这是强制性的，模型重力往上，别人很不习惯。问题：1。K0值为何取0。42。由于初始自平衡导致网格已经下沉，所以实际开挖只有一层3。shell设置不上都与网格变形有关4。建议用setsmall测试是说这句：setmechratio=1e-4这样设能少花算的时间！setsmall,不知斑竹能否告知和setlarge的结果有多大差别！我看网上的别人的文件，别人在用modelnull时，range好象比实际都大一些，是不是我实际开挖x方向0-20m,z方向0-2m,应该rangex-0.120.1z-0.12.1qinjianshe斑竹：我去掉了setlarge计算，开挖的数目就对了，shell布置的数目也对了！qinjianshe的作法无可非议我觉得还有一个方法，因为从你们的讨论似乎可以看出，你的模型在达到初始平衡状态时已产生破坏，故在你以后的计算中，就产生了上述的问题。如果确是如此，可将coh设的大些，比如1e10，在模型达到初始平衡之后，在将其改为1.5e4，然后进行后续计算，或许能够解决你的问题，而不改变setlarge。是不是可以得出这样的结论：我计算出问题和我所设置的土性参数有关！和我划分网格的方式，以及我所用的命令流格式无关。问题1：eplier，你的文件里prop时没有ten