还剩20页未读,继续阅读
本资源只提供10页预览,全部文档请下载后查看!喜欢就下载吧,查找使用更方便
文本内容:
分差法岩土工程结构的数值解是建立在满足基本方程(平衡方程、几何方程、本构方程)和边界条件下推导的由于基本方程和边界条件多以微分方程的形式出现,因此,将基本方程近假发改用差分方程(代数方程)表示,把求解微分方程的问题改换成求解代数方程的问题,这就是所谓的分差法差分法由来已久,但差分法需要求解高阶代数方程组,只有在计算机的出现,才使该法得以实施和发展FLAC3D(FastLagrangian____ysisofContinua)由美国Itasca公司__的目前,FLAC有二维和三维计算程序两个版本,二维计算程序V
3.0以前的为DOS版本,V
2.5版本仅仅能够使用计算机的基本内存64K),所以,程序求解的最大结点数仅限于2000个以内1995年,FLAC2D已升级为V
3.3的版本,其程序能够使用护展内存因此,大大发护展了计算规模FLAC3D是一个三维有限差分程序,目前已发展到V21版本FLAC3D的输入和一般的数值分析程序不同,它可以用交互的方式,从键盘输入各种命令,也可以写成命令(集)文件,类似于批处理,由文件来驱动因此,采用FLAC程序进行计算,必须了解各种命令关键词的功能,然后,按照计算顺序,将命令按先后,依次排列,形成可以完成一定计算任务的命令文件FLAC3D是二维的有限差分程序FLAC2D的护展,能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析调整三维网格中的多面体单元来拟合实际的结构单元材料可采用线性或非线性本构模型,在外力作用下,当材料发生屈服流动后,网格能够相应发变形和__(大变形模式)FLAC3D采用的显式拉格朗日算法和混合-离散分区技术能够非常准确发模拟材料的塑性破坏和流动由于无须形成刚度矩阵,因此,基于较小内存空间就能够求解大范围的三维问题FLAC3D采用ANSIC++语言编写的FLAC3D有以下几个优点1对模拟塑性破坏和塑性流动采用的是“混合离散法“这种方法比有限元法中通常采用的“离散集成法“更为准确、合理2即使模拟的系统是静态的,仍采用了动态运动方程,这使得FLAC3D在模拟物理上的不稳定过程不存在数值上的障碍3采用了一个“显式解“方案因此,显式解方案对非线性的应力-应变关系的求解所花费的时间,几互与线性本构关系相同,而隐式求解方案将会花费较长的时间求解非线性问题面且,它没有必要存储刚度矩阵,这就意味着;采用中等容量的内存可以求解多单元结构;模拟大变形问题几互并不比小变形问题多消耗更多的计算时间,因为没有任何刚度矩阵要被修改当然,它也存在以下几个不足之处1对于线性问题的求解,FLAC3D比有限元程序运行得要慢;因此,当进行大变形非线性问题或模拟实际可能出现不稳定问题时,FLAC3D是最有效的工具2用FLAC3D求解时间取决于最长的自然周期和最短的自然周期之比但某些问题对模型是无效的接上面的介绍Flac3D中为岩土工程问题的求解__了特有的本构模型,总共包含了10种材料模型
1.开挖模型null
2.3个弹性模型(各向同性,横观各向同性和正交各向同性弹性模型)36个塑性模型(Drucker-Prager模型、Morh-Coulomb模型、应变硬化/软化模型、遍布节理模型、双线性应变硬化/软化遍布节理模型和修正的cam粘土模型)Flac3D网格中的每个区域可以给以不同的材料模型,并且还允许指定材料参数的统计分布和变化梯度还包含了节理单元,也称为界面单元,能够模拟两种或多种材料界面不同材料性质的间断特性节理允许发生滑动或分离,因此可以用来模拟岩体中的断层、节理或摩擦边界FLAC3D中的网格生成器gen,通过匹配、连接由网格生成器生成局部网格,能够方便地生成所需要的三维结构网格还可以自动产生交叉结构网格(比如说相交的巷道),三维网格由整体坐标系xyz系统所确定,不同于FLAC程序是由行列方式确定这就提供了比较灵活的产生和定义三维空间参数FLAC3D做计算分析的一般步骤与大多数程序采用数据输入方式不同,FLAC采用的是命令驱动方式命令字控制着程序的运行在必要时,尤其是绘图,还可以启动FLAc用户交互式图形界面为了建立FLAC计算模型,必须进行以下三个方面的工作
1.有限差分网格
2.本构特性与材料性质
3.边界条件与初始条件完成上述工作后,可以获得模型的初始平衡状态,也就是模拟开挖前的原岩应力状态然后,进行工程开挖或改变边界条件来进行工程的响应分析,类似于FLAC的显式有限差分程序的问题求解与传统的隐式求解程序不同,FLAC采用一种显式的时间步来求解代数方程进行一系列计算步后达到问题的解在FLAC中,达到问题所需的计算步能够通过程序或用户加以控制,但是,用户必须确定计算步是否已经达到问题的最终的解最后,进行结果的分析与总结,这就是用户自己的事情了,呵呵,在手册中的图中有显示shell单元但我自己做的时候不知用何命令才能在网格图中将的shell,beam等单元显示出来,高手指点!!plotaddselgeom另外,还可在上面命令的后面加cid=onnode=off.....再问一下,shell是否只可以加在网格的表面?怎么在网格内部加不上?可以加上了,但是不知道厚度t如何表现thickness啊gensu_____有什么用?生成表面显示模型中用到setdyn=ongr__00-10hist_rep=1是在动力分析中的问题,请问最后一个变量是什么含义,好像只能是整数,查了Fish说明也没有看到解释,请赐教!hist_repnHistoriesaresampledeveryntimesteps.Thedefaultisn=
10.ThesynonymHISTORYnstep__yalsobeused.每10个步时记录一次请教:1谁用过GENERATEsu_____xarc命令生成过圆弧面,能否贴出来让小弟学一下我曾按照__nual作过一些尝试但没有成功2根据说明我曾生成过圆形__模型,但如何在__内壁上加正应力,请各位高手指点加正应力应该是用nstress命令___做过曲面的例子,只用过规则的圆柱面请问各位大侠,groutstiffness是什么意思?
1.gr_kgroutstiffnessperunitlengthkg[F/L2]
2.Ft=kg*ut
3.Thegroutannulusisassumedtobeh__easanelastic-perfectlyplasticsolid.Asaresultofrelativesheardispla__mentutbetweenthetendonsu_____andtheboreholesu_____theshearfor__Ftmobilizedperlengthofcableisrelatedtothegroutstiffnesskg—i.e.
4.可见k_g是groutannulus与zone(soil)之间的剪切刚度(单位剪切位移引起的应力)
5.详细参考cable一节,有图形表示,很清楚写了几个命令流文件,但只要所分析问题的单元数一多,计算就很难收敛,总是在e-3或e-4之间跳来跳去,都迭代几万步了请教大侠,怎么加快收敛速度?一般一万五千单元左右要算多久?flac中的求解控制通常有1)setmechfor__50solve2)solveratio1e-73)step30000等几种,不知这几种怎么选择?每种有什么特点??加快收敛速度,要让单元的尺寸相差不大,另外,ratio恒定,说明产生了塑性流动或者破坏了,看看约束条件对不对!!材料模型对不对,材料参数是否太离谱??flac中的求解控制通常有1)setmechfor__50solve2)solveratio1e-73)step30000第一种是控制不平衡力,是给定具体的值,必须迭代到这个值才结束运算这是绝对的值,个人认为只有针对具体问题才使用;第二种要重先设定,程序默认1e-5,-7也太小了,迭代的步数更多,更精确,一般情况下,个人认为-5就够了;第三种是给定迭代步数,不至于30000吧,一般3000吧但这种设定不一定能达到1e-5,所以要在有把握的情况下使用至于在循环中,solve与stepn到底有何差别,本人也不甚清楚有时某一步可能需要7000步才能达到-5次方,而其他的循环只需要1000步就达到了,这种情况下究竟该如何处理?请高手指点!solve是通过控制不平衡力或比例的大小决定计算步step,cyc等命令则是由__控制计算迭代步,在无法用step估计计算是否收敛情况下,多用solve,step自身的优势可以应用于逐步加载等方面加大bulk,shear的量纲试一试,还有就是检查一下边界条件!;产生一个体,这个体里面包含__genzoneradbrickp0000p11000p20100p30010size3557ratio
1111.5dim142fill;设置模型服从摩尔库仑屈服准则modmoh;定义bulk模量,剪切模量,摩擦角,粘聚力,和tenspropbulk1e8shear3e8fric35coh1e3tens1e3;施加约束fixxrangex-
0.
10.1fixzrangez-
0.
10.1fixyrangey
9.
910.1;施加重力setgr__0100;设置密度inidens1000;设置初始应力,这个有什么用呢????inisyy
0.0grad0-_____0inisxx
0.0grad0-50000iniszz
0.0grad0-50000applysxx
0.0grad0-50000rangex
9.
910.1applyszz
0.0grad0-50000rangez
9.
910.1;挖坑的命令,使用fish语言defdigHolecom__ndrangenametrenchx01y04z02;挖了modelnulrangetrenchsetlargehistgpxdisp100endcom__ndend;表面蓝色显示plotsurfblue;停顿pause;开挖digHole;设置大变形setlarge;记录位移变化histgpxdisp100;...andthencalculatetheresponse;设置计算步为1000step1000;我不知道上述解释是否明白,是否对头??cjiao2000wrote:;产生一个体,这个体里面包含__genzoneradbrickp0000p11000p20100p30010size3557ratio
1111.5dim142fill;设置模型服从摩尔库仑屈服准则modmoh;定义bulk模量,剪切模量,摩擦角,粘聚力,和tenspropbulk1e8shear3e8fric35coh1e3tens1e3;施加约束fixxrangex-
0.
10.1fixzrangez-
0.
10.1fixyrangey
9.
910.1;施加重力setgr__0100;设置密度inidens1000;设置初始应力,这个有什么用呢????inisyy
0.0grad0-_____0inisxx
0.0grad0-50000iniszz
0.0grad0-50000别的软件很少可以加地应力的,怎么说没用呢?这可是Flac的专长呀,象ansys也没的办法的applysxx
0.0grad0-50000rangex
9.
910.1applyszz
0.0grad0-50000rangez
9.
910.1;挖坑的命令,使用fish语言defdigHolecom__ndrangenametrenchx01y04z02;挖了modelnulrangetrenchsetlargehistgpxdisp100endcom__ndend;表面蓝色显示plotsurfblue;休息pause;开挖digHole;设置大变形setlarge;这个好像是施加位移我认为是记录坐标是(1,0,0,)点的位移历史记录吧,不是施加位移吧权当抛砖引玉了,各位大侠不要见笑histgpxdisp100;...andthencalculatetheresponse;设置计算步为1000step1000;我不知道上述解释是否明白,是否对头??首先确定位移是z向,还是x向histgpxdisp100有时候会出现位移都为0的情况怎么回事明明那是有位移的啊是啊,还有向上得哦!!!;设置初始应力,这个有什么用呢????inisyy
0.0grad0-_____0inisxx
0.0grad0-50000iniszz
0.0grad0-50000别的软件很少可以加地应力的,怎么说没用呢?这可是Flac的专长呀,象ansys也没的办法的设置初始应力做什么哦???一般情况下围岩地应力包括自重应力、构造应力、采掘应力你所描速的情况应是在挖掘深度不大的情况下,只考虑土体的自重应力对开挖的影响在古典压力理论中,认为地层的水平压力的侧压系数为u/1-u,这里取
0.5apply描述的模型的边界条件histgpxdisp100是用于记录点(1,0,0)在开挖过程中的位移情况,模拟的目的就是要研究目标位置的位移或者应力变化情况setlarge设置大变形的话应该是按大变形理论吧!最近做了算例,发现将cable单元的pretension设置为0或者100没有什么区别,大家看一下是怎么回事啊?;generatethemodletitleVerticalAnchoredGr__ityWall;setmecharatio1e-5;setlarge;Zonegenzoneradtunnelp
01.50-20p
11.520-20p200-20p
31.500p4020-20p5000p
61.5200p70200p
81.56-20p
91.50-
10.5p1006-20p1100-
10.5p
121.56-
10.5p1306-
10.5size1231920ratio1111dim
9.
559.55fill;__terial;Forbalan__modelmohrpropdensity=
1.80e3bulk=
12.5e6shear=
6.5e6coh
10.0e3tens50e10friction=
10.0dilation=
0.0rangez-
20.0-
18.0propdensity=
2.00e3bulk=
20.0e6shear=
12.0e6coh
15.0e3tens50e10friction=
25.0dilation=
5.0rangez-
18.0-
10.0propdensity=
2.00e3bulk=
44.4e6shear=
33.3e6coh
0.0e3tens50e10friction=
30.0dilation=
10.0rangez-
10.
00.0;Boundaryconditonfixxrangex-
0.
10.1fixxrangex
0.
41.5fixyrangey-
0.
1.1fixyrangey
19.
920.1fixzrangez-
0.
10.1;Setgr__itionalaccsetgr__0010;Initialstressiniszz-
4.0e5grad00-
2.0e4inisxx-
2.0e5grad00-
1.00e4inisyy-
2.0e5grad00-
1.00e4;InitialStateSolveforbalan__solveinistate0inixdisp
0.0ydisp
0.0zdisp
0.0;Reinstallthesoilmodelmohrpropdensity=
1.80e3bulk=
12.5e6shear=
6.5e6coh
10.0e3tens58e3friction=
10.0dilation=
0.0rangez-
20.0-
18.0propdensity=
2.00e3bulk=
20.0e6shear=
12.0e6coh
15.0e3tens
32.6e3friction=
25.0dilation=
5.0rangez-
18.0-
10.0propdensity=
2.00e3bulk=
44.4e6shear=
33.3e6coh
0.0e3tens5e3friction=
30.0dilation=
10.0rangez-
10.
00.0;Setthesurchargeapplyszz10e3rangey620z-
20.1-
19.9solve;Begintoexc__ate;Step1:
1.5moffmodnullrangey
06.0z-
18.5-
20.0sellinerrangey
5.
96.1z-
18.5-
20.0sellinerPROPiso=25e
90.15thick=
0.1;concretesellinerPROPcs_nk=8e8cs_sk=8e8cs_ncut=4e6cs_scoh=4e4cs_scohres=2e6cs_sfric=
20.0solves__estep
1.s__;Step2:setthesoilnailand
1.5moffselcableid=1begin=
0.
756.0-
18.5end=
0.
7515.5-
17.5nseg=12pretention50e3selcablepropxcarea=
5.1e-3emod=200e9ytens=210e6gr_k=6e6gr_coh=40e3selcablepropgr_per=
0.___gr_fric=
25.0modnullrangey
06.0z-
17.0-
18.5sellinerrangey
5.
96.1z-
17.0-
18.5sellinerPROPiso=25e
90.15thick=
0.1;concretesellinerPROPcs_nk=8e8cs_sk=8e8cs_ncut=4e6cs_scoh=4e4cs_scohres=2e6cs_sfric=
20.0solves__estep
2.s__;Step3:setthesoilnailand
1.5moffselcableid=2begin=
0.
756.0-
17.0end=
0.
7515.5-
16.0nseg=12pretention50e3selcablepropxcarea=
5.1e-3emod=200e9ytens=210e6gr_k=6e6gr_coh=40e3selcablepropgr_per=
0.___gr_fric=
25.0modnullrangey
06.0z-
15.5-
17.0sellinerrangey
5.
96.1z-
15.5-
17.0sellinerPROPiso=25e
90.15thick=
0.1;concretesellinerPROPcs_nk=8e8cs_sk=8e8cs_ncut=4e6cs_scoh=4e4cs_scohres=2e6cs_sfric=
20.0solves__estep
3.s__;Step4:setthesoilnailand
1.5moffselcableid=3begin=
0.
756.0-
15.5end=
0.
7515.5-
14.5nseg=12pretention50e3selcablepropxcarea=
5.1e-3emod=200e9ytens=210e6gr_k=6e6gr_coh=40e3selcablepropgr_per=
0.___gr_fric=
25.0modnullrangey
06.0z-
14.0-
15.5sellinerrangey
5.
96.1z-
14.0-
15.5sellinerPROPiso=25e
90.15thick=
0.1;concretesellinerPROPcs_nk=8e8cs_sk=8e8cs_ncut=4e6cs_scoh=4e4cs_scohres=2e6cs_sfric=
20.0solves__estep
4.s__;Step5:setthesoilnailand
1.5moffselcableid=3begin=
0.
756.0-
14.0end=
0.
7515.5-
13.0nseg=12pretention50e3selcablepropxcarea=
5.1e-3emod=200e9ytens=210e6gr_k=6e6gr_coh=40e3selcablepropgr_per=
0.___gr_fric=
25.0modnullrangey
06.0z-
12.5-
14.0sellinerrangey
5.
96.1z-
12.5-
14.0sellinerPROPiso=25e
90.15thick=
0.1;concretesellinerPROPcs_nk=8e8cs_sk=8e8cs_ncut=4e6cs_scoh=4e4cs_scohres=2e6cs_sfric=
20.0solves__estep
5.s__;Step6:setthesoilnailand
1.5moffselcableid=3begin=
0.
756.0-
12.5end=
0.
7515.5-
11.5nseg=12pretention50e3selcablepropxcarea=
5.1e-3emod=200e9ytens=210e6gr_k=6e6gr_coh=40e3selcablepropgr_per=
0.___gr_fric=
25.0modnullrangey
06.0z-
11.5-
12.5sellinerrangey
5.
96.1z-
11.5-
12.5sellinerPROPiso=25e
90.15thick=
0.1;concretesellinerPROPcs_nk=8e8cs_sk=8e8cs_ncut=4e6cs_scoh=4e4cs_scohres=2e6cs_sfric=
20.0solves__estep
6.s__;Step7:setthesoilnailand
1.5moffselcableid=3begin=
0.
756.0-
11.0end=
0.
7515.5-
10.0nseg=12pretention50e3selcablepropxcarea=
5.1e-3emod=200e9ytens=210e6gr_k=6e6gr_coh=40e3selcablepropgr_per=
0.___gr_fric=
25.0modnullrangey
06.0z-
10.5-
11.5sellinerrangey
5.
96.1z-
10.5-
11.5sellinerPROPiso=25e
90.15thick=
0.1;concretesellinerPROPcs_nk=8e8cs_sk=8e8cs_ncut=4e6cs_scoh=4e4cs_scohres=2e6cs_sfric=
20.0solves__estep
7.s__flac怎么作一个剖面!就是我建了一个模型,___一个切面的视图plsetplaneori000nor100plcondispplane在所建立的模型中有两个命令;就是printgpposition与printzonegp不知道这两个命令中的gp是否有__?按找手册中的解释,printgpposition是用来输出gp的位置坐标的;printzonegp是用来输出与zone相关的gp的id值的这样我就有一个上面的疑问,两个命令中的gp都有id值,是不是不同命令中的gp的id值相同的话,他们的坐标就可以看作是一样的啊?是一个gp有且只有一个id谢谢dyna__x兄这下我就明白了,如果要画等值线的话,先输出gp的坐标来,将其另存为一个文件,再输出zone的应力值,将与zone对应的gp的id与前面输出的gp的id值对应起来,找到gp的坐标用ex__l打开上面的文件进行编辑,最后用surfer画等值线这样的办法太笨,不知到有没有更好的办法啊?欢迎前来讨论你的思路很对!要简化操作,如果你对fish熟悉,可以直接按surfer的文件格式输出如果会用tecplot,可以参考前面我那篇【原创】用tecplot对flac3d计算结果进行后处理surfer我不熟,原则应该是一样的看了你的哪个原创,发现里面有这样几个语句不太理解,IO_READ=0IO_WRITE=1IO_FISH=0IO_ASCII=1其中的0,1具体指的是什么啊?我发现一般zone中显示应力的gp的id往往不是全部模型的gp,那么如果想要用fish语句输出这些gp的坐标,该如和判断呢?不知道在fish语句中如何将zone中的gp的坐标给出可以将gp的坐标和应力值存到一个数据文件,然后用surfer调用数据文件即可,不须在做额外的编辑比如这个命令printzones输出了所有区的应力,可是要如果把它保存到指定写字板呢setlogonsetlogfile*.logprintzones依次执行以上命令便可以实现你的目的其实简单的理解就是dd映射平面法向与y轴正向夹角(顺时针)dip映射平面法向与z轴负向夹角dip是倾角,就是plane和水平面的夹角dd是dipdirection,倾角的方向,注意不是走向就是你用来看倾角的那根线,在水平面上的投影,与y轴正向的交角从y轴正向开始算那位知道flac3D与ansys中单元不等比例剖分中网格长度与比例系数之间的计算公式s=[a1-q^n-1]/1-qa第一个网格长度q比例系数不知版主问的是不是这个公式等比公式在做边坡稳定性分析建模的时候一定会遇到非自立边坡此时生成模型的时候有几种方法.大家用什么方法
1.用不同的zonebrick拼装起来;
2.用planeabove该命令切出来的斜坡呈锯齿状很难看除非网格很细否则与实际偏离较大
3.使用fish手册上有一个例子通过改变各个grid的坐标来实现__的球形开挖面依此类推一定可以修改grid的坐标实现斜坡.用fish或用zonebrick拼装生成非90度的基坑边坡,能给个例子吗?下面就是我编的,请大家帮我看看,到底哪里有问题,怎么老是有问题(我观察到第一步开挖后,图上显示只挖掉了一排网格,我要挖的可是两排在第二次开挖前,shellSEL竟然只布置上去2个!)我改了好多次都是类似结果,不知到底我哪弄错了!genzoneradbrickp16000p20600p30050size20201032dim202010rat
1111.05fillmodelmohrpropbulk
1.25e7shear
5.77e6fric20coh
1.5e4rangez
0.
050.0;边界条件fixxrangex-
0.
010.01fixyrangey-
0.
010.01fixxyzrangez
49.
9950.01fixxyrangex
59.
9960.01fixxyrangey
59.
9960.01;初始条件inidensity
1.9e3rangez
0.
0050.00setgr__ity0010iniszz0grad
001.9e4rangez
0.
0050.00inisxx0grad
007.6e3rangez
0.
0050.00inisyy0grad
007.6e3rangez
0.
0050.00setlargehistid1unbalplothist1setmechratio=1e-4solves__ej_ini.s__inixdis=0ydis=0zdis=0grouplod1rangex2030y030z-
0.
010.01grouplod2rangex030y2030z-
0.
010.01grouplodrangegrouplod1anygrouplod2anyapplyszz10e3rangegrouplodmodelnullrangex020y020z
0.
02.0step2000s__ej_eca
1.s__selshellid=1rangex
19.
920.1y020z
0.
02.0selshellid=2rangex020y
19.
920.1z
0.
02.0selnodefixxyrzrrangex-
0.
010.01selnodefixyxrzrrangey-
0.
010.01selshellid=1propiso=
26.6e
90.25thick=
0.1selshellid=2propiso=
26.6e
90.25thick=
0.1modelnullrangex020y020z
2.
04.0step2000s__ej_eca
2.s__selshellid=3rangex
19.
920.1y020z
2.
04.0selshellid=4rangex020y
19.
920.1z
2.
04.0selnodefixxyrzrrangex-
0.
010.01selnodefixyxrzrrangey-
0.
010.01selshellid=3propiso=
26.6e
90.25thick=
0.1selshellid=4propiso=
26.6e
90.25thick=
0.1modelnullrangex020y020z
4.
06.0step2000s__ej_eca
3.s__selshellid=5rangex
19.
920.1y020z
4.
06.0selshellid=6rangex020y
19.
920.1z
4.
06.0selnodefixxyrzrrangex-
0.
010.01selnodefixyxrzrrangey-
0.
010.01selshellid=5propiso=
26.6e
90.25thick=
0.1selshellid=6propiso=
26.6e
90.25thick=
0.1modelnullrangex020y020z
6.
08.0step2000s__ej_eca
4.s__selshellid=7rangex
19.
920.1y020z
6.
08.0selshellid=8rangex020y
19.
920.1z
6.
08.0selnodefixxyrzrrangex-
0.
010.01selnodefixyxrzrrangey-
0.
010.01selshellid=7propiso=
26.6e
90.25thick=
0.1selshellid=8propiso=
26.6e
90.25thick=
0.1modelnullrangex020y020z
8.
010.0step2000s__ej_eca
5.s__selshellid=9rangex
19.
920.1y020z
8.
010.0selshellid=10rangex020y
19.
920.1z
8.
010.0selnodefixxyrzrrangex-
0.
010.01selnodefixyxrzrrangey-
0.
010.01selshellid=9propiso=
26.6e
90.25thick=
0.1selshellid=10propiso=
26.6e
90.25thick=
0.1step2000s__ej_finl.s__完整的命令流见附件genzoneradbrickp16000p20600p30050size20201032dim202010rat
1111.05fillmodelmohrpropbulk
1.25e7shear
5.77e6fric20coh
1.5e4rangez
0.
050.0;边界条件fixxrangex-
0.
010.01fixyrangey-
0.
010.01fixxyzrangez
49.
9950.01fixxyrangex
59.
9960.01fixxyrangey
59.
9960.01;初始条件inidensity
1.9e3rangez
0.
0050.00setgr__ity0010iniszz0grad
001.9e4rangez
0.
0050.00inisxx0grad
007.6e3rangez
0.
0050.00inisyy0grad
007.6e3rangez
0.
0050.00setlargehistid1unbalplothist1setmechratio=1e-4solves__ej_ini.s__inixdis=0ydis=0zdis=0grouplod1rangex2030y030z-
0.
010.01grouplod2rangex030y2030z-
0.
010.01grouplodrangegrouplod1anygrouplod2anyapplyszz10e3rangegrouplodmodelnullrangex020y020z
0.
02.0step2000s__ej_eca
1.s____crotudingemod2e11xcarea
3.81e-4ytens2e8gr_coh
2.5e4gr_k
6.33e6gr_per
0.346gr_fric25;设置cabledefpla___cable1xx1=1loopwhilexx120com__ndselcablebeg20xx
11.5end
25.908847xx
12.5418__nseg6selcablebegxx
1201.5endxx
125.
9088472.5418__nseg6end_com__ndxx1=xx1+
1.5end_loopendpla___cable1selcableproptudingselshellid=1rangex
19.
920.1y020z
0.
02.0selshellid=2rangex020y
19.
920.1z
0.
02.0selnodefixxyrzrrangex-
0.
010.01selnodefixyxrzrrangey-
0.
010.01selshellid=1propiso=
26.6e
90.25thick=
0.1selshellid=2propiso=
26.6e
90.25thick=
0.1modelnullrangex020y020z
2.
04.0step2000s__ej_eca
2.s__defpla___cable2xx2=
0.5loopwhilexx220com__ndselcablebeg20xx
23.0end
25.908847xx
24.0418__nseg6selcablebegxx
2203.0endxx
225.
9088474.0418__nseg6end_com__ndxx2=xx2+
1.5end_loopendpla___cable2selcableproptudingselshellid=3rangex
19.
920.1y020z
2.
04.0selshellid=4rangex020y
19.
920.1z
2.
04.0selnodefixxyrzrrangex-
0.
010.01selnodefixyxrzrrangey-
0.
010.01selshellid=3propiso=
26.6e
90.25thick=
0.1selshellid=4propiso=
26.6e
90.25thick=
0.1modelnullrangex020y020z
4.
06.0step2000s__ej_eca
3.s__defpla___cable3xx3=1loopwhilexx320com__ndselcablebeg20xx
34.5end
25.908847xx
35.5418__nseg6selcablebegxx
3204.5endxx
325.
9088475.5418__nseg6end_com__ndxx3=xx3+
1.5end_loopendpla___cable3selcableproptudingselshellid=5rangex
19.
920.1y020z
4.
06.0selshellid=6rangex020y
19.
920.1z
4.
06.0selnodefixxyrzrrangex-
0.
010.01selnodefixyxrzrrangey-
0.
010.01selshellid=5propiso=
26.6e
90.25thick=
0.1selshellid=6propiso=
26.6e
90.25thick=
0.1modelnullrangex020y020z
6.
08.0step2000s__ej_eca
4.s__defpla___cable4xx4=
0.5loopwhilexx420com__ndselcablebeg20xx
46.0end
25.908847xx
47.0418__nseg6selcablebegxx
4206.0endxx
425.
9088477.0418__nseg6end_com__ndxx4=xx4+
1.5end_loopendpla___cable4selcableproptudingdefpla___cable5xx5=1loopwhilexx520com__ndselcablebeg20xx
57.5end
25.908847xx
58.5418__nseg6selcablebegxx
5207.5endxx
525.
9088478.5418__nseg6end_com__ndxx5=xx5+
1.5end_loopendpla___cable5selcableproptudingselshellid=7rangex
19.
920.1y020z
6.
08.0selshellid=8rangex020y
19.
920.1z
6.
08.0selnodefixxyrzrrangex-
0.
010.01selnodefixyxrzrrangey-
0.
010.01selshellid=7propiso=
26.6e
90.25thick=
0.1selshellid=8propiso=
26.6e
90.25thick=
0.1modelnullrangex020y020z
8.
010.0step2000s__ej_eca
5.s__defpla___cable6xx6=1loopwhilexx620com__ndselcablebeg20xx
69.0end
25.908847xx
610.0418__nseg6selcablebegxx
6209.0endxx
625.
90884710.0418__nseg6end_com__ndxx6=xx6+
1.5end_loopendpla___cable6selcableproptudingselshellid=9rangex
19.
920.1y020z
8.
010.0selshellid=10rangex020y
19.
920.1z
8.
010.0selnodefixxyrzrrangex-
0.
010.01selnodefixyxrzrrangey-
0.
010.01selshellid=9propiso=
26.6e
90.25thick=
0.1selshellid=10propiso=
26.6e
90.25thick=
0.1step2000;水平位移变化plotcreatedispla__plotsetcapsize20plotsetrot2000195plotaddcontxdispoutlonplotaddselcablefor__plotaddaxesblackplotcreatetu_dispplotsetcapsize20plotsetrot1800180plotaddsketchplotadddispredplotaddselgeombluebluenode=offs__ej_finl.s__tension=cohesion/tanfriction好像在flac的theory里面有,就是莫尔园的外包线和x轴的交点横坐标mohr模型有shearfailure和tensionfailure,如果把tension设置很小,有可能拉应力很小就达到塑性了,所以最好还是按前面那个关系设回eplier 你的程序___运行过,但根据我个人的推断,或许是你的模型在达到初始平衡后已经破坏,故如果继续进行后续计算,就可能产生你所述的情况 所以,如将初始粘聚力值增大,使之在初始平衡时不会产生塑性破坏,在进行后续计算之前,在将粘聚力恢复为你所选择的参数,应该可以避免你的情况 你可以试试,也可以参考一下__nual中相关的说明 如果有什么结果,还望在此告知你的边界条件设置的合理吗?;边界条件fixxrangex-
0.
010.01fixyrangey-
0.
010.01fixxyzrangez
49.
9950.01fixxyrangex
59.
9960.01fixxyrangey
59.
9960.01;初始条件inidensity
1.9e3rangez
0.
0050.00setgr__ity0010iniszz0grad
001.9e4rangez
0.
0050.00inisxx0grad
007.6e3rangez
0.
0050.00inisyy0grad
007.6e3rangez
0.
0050.00你所建立的模型应该是真实工程中的一半吧,你这个模型不复杂,你可以将整个模型建立起来,进行运算,如果单就半个模型来说,我想应该在x-z平面加上应力边界条件,当初建议__用set__all计算只是为了让其发现单元开挖一直有问题这个具体问题实际应用中,大变形是flac的特色,而且如果分析对象是大变形问题,建议还是应用setlarge至于单元开挖有误,说明__程序肯定有问题最近我正仔细研究flac计算原理,并阅读连续介质力学计算方法,收获很大,觉得flac的__还需要很多工作要做我正在整理,大约不长时间将会展开专门的详细论述对的,基坑开挖是不可以用set__all的,必须用setlarge,因为基坑在变形破坏时一定是单元与原来的网格发生了比较大的位移,如果用小应变模式,那就是你限制了单元不让其破坏,这是强制性的,模型重力往上,别人很不习惯问题1K0值为何取042由于初始自平衡导致网格已经下沉,所以实际开挖只有一层3shell设置不上都与网格变形有关4建议用set__all测试是说这句setmechratio=1e-4这样设能少花算的时间!set__all不知斑竹能否告知和setlarge的结果有多大差别!我看网上的别人的文件,别人在用modelnull时,range好象比实际都大一些,是不是我实际开挖x方向0-20mz方向0-2m应该rangex-
0.
120.1z-
0.
12.1qinjianshe斑竹我去掉了setlarge计算,开挖的数目就对了,shell布置的数目也对了!qinjianshe的作法无可非议我觉得还有一个方法,因为从你们的讨论似乎可以看出,你的模型在达到初始平衡状态时已产生破坏,故在你以后的计算中,就产生了上述的问题如果确是如此,可将coh设的大些,比如1e10,在模型达到初始平衡之后,在将其改为
1.5e4,然后进行后续计算,或许能够解决你的问题,而不改变setlarge是不是可以得出这样的结论我计算出问题和我所设置的土性参数有关!和我划分网格的方式,以及我所用的命令流格式无关问题1eplier,,,你的文件里prop时没有tens这个值,这个可以省掉吗?问题2深基坑模拟时用大变形模式和小变形模式对结果的影响到底是多大?可以肯定,大部分的程序问题出在材料参数有关,少量程序出问题在网格和命令格式上最近学习flac3d,遇到一个困难,就是边坡有不同的岩层分布,同一高层有不同的岩层,不知道如何正确赋材料参数,请高手指教一二感激之至建模的时候先定义出不同的group然后直接对group赋材料参数用Geo-slope/sig__模拟基坑开挖支护时可以设立时间步骤每一步开挖多少单元把这些单元设为相应的步骤即可不知Flac3d有没有类似的功能还是要多次使用Free命令请指教.感谢taizhong朋友的回复.用modelnullrangegroupXXX可以模拟开挖我还见到一个例题使用free命令的.不知有哪位大牛比较过这两种方法下面是使用free命令模拟开挖的命令流:newtitleSoilnailingusingcablesgenzonebricksize11111modelmohrpropbulk5e9shear1e9coh4e4fric30inidens2000fixxzrangez-.
1.1fixxrangex-.
1.1fixxrangex
10.
911.1fixy;iniszz=-
2.2e5grad002e4inisxx=-
1.32e5grad
001.2e4inisyy=-
0.88e5grad
000.8e4setgr__00-10;histunbalhistgpxdisp0011;solve;freexrangex-.
1.1;selcableid=1begin=
0.
00.
53.5end=
80.
53.5nseg=8selcableid=2begin=
0.
00.
56.5end=
80.
56.5nseg=8selcableid=3begin=
0.
00.
59.5end=
80.
59.5nseg=8selcablepropxcarea=
8.5e-3emod=200e9ytens=1e10gr_k=7e6gr_coh=1e2selcablepropgr_per=
0.___gr_fric=25;plotcreatesnailplotsetcapsize25plotaddsketchplotadddisp__x=2e-2plotaddselcablefor__yellowredplotaddselcablegroutslipplotaddselgeomnode=offblackblackplotshow;solves__ecable-sn.s__我在求解是打算用最大不平衡力来确定计算结束的时间,不知那位高手能说说不平衡力的确定方法谢谢了!!!!!!显示有限差分法与隐示有限差分法的区别是什么?显示有限差分法与隐示有限差分法的定义分别是什么?显式方程式右边全是已知量,隐式反之ntitle壳支撑genzonebrickp13000p20200p30010size15105rat111grouptuplotblgrplotshowmodelelasrangegrouptuprops=119047619b=
2083333.3dens=2200rangegrouptugroupexcarangex
15.
930.1y
9.
920.1z
5.
910.1fixxrangex-
0.
10.1fixxrangex
29.
930.1fixzrangez-
0.
10.1fixzrangez
9.
910.1fixxyzrangey-
0.
10.1setlarge;设置大变形setgr__0-
9.80step1500s__tu.s__;这是不加shell单元,去掉分号后为加上壳单元,怎么加不加没啥区别呢?;selshellid=1turangex
15.
930.1y
2.
920.1z
5.
96.1;selshellid=2turangex
15.
916.1y
2.
920.1z
5.
910.1;selshellpropisotropic
10.5e
90.25thickness
0.3density2500;selnodefixzrangeid=1z
9.
910.1;selnodefixxrangeid=2x
29.
930.1;plotaddselgeomblueyellowplotshowinitxdisp=0ydisp=0zdisp=0modelnullrangegroupexcastep1500s__exca.s__最大不平衡力除以施加的初始力(主意,应力×作用__)可以得到一个比值,计算过程中,当这个比值达到1e-5时,就可以认为状态stable了,这个值是可以事先设定的反推回去,就可以得到所要的最大不平衡力但比值一定不能设定太大。