隧道及地下工程FLAC解析方法-常识、建模、常用命令流及其解释

1、隧道及地下工程FLAC解析方法 _常识、建模、常用命令流及其解释TOC o 1-1 h z u HYPERLINK l _Toc320528144 1.1 FLAC常识 PAGEREF _Toc320528144 h 2 HYPERLINK l _Toc320528145 1.2常用用命令流流 PAGEREF _Toc320528145 h 13 HYPERLINK l _Toc320528146 1.3建模模过程 PAGEREF _Toc320528146 h 21 HYPERLINK l _Toc320528147 2-1定义义一个FFISHH函数 PAGEREF _Toc32052814

2、7 h 24 HYPERLINK l _Toc320528148 2-2使用用一个变变量 PAGEREF _Toc320528148 h 24 HYPERLINK l _Toc320528149 2-3对变变量和函函数的理理解 PAGEREF _Toc320528149 h 24 HYPERLINK l _Toc320528150 2-4获取取变量的的历史记记录 PAGEREF _Toc320528150 h 24 HYPERLINK l _Toc320528151 2-5用FFISHH函数计计算体积积模量和和剪砌模模量 PAGEREF _Toc320528151 h 25 HYPERLINK

3、 l _Toc320528152 2-6 在在FLAAC输入入中使用用符号变变量 PAGEREF _Toc320528152 h 25 HYPERLINK l _Toc320528153 2-7 控控制循环环 PAGEREF _Toc320528153 h 26 HYPERLINK l _Toc320528154 2-8 拆拆分命令令行 PAGEREF _Toc320528154 h 26 HYPERLINK l _Toc320528155 2-9 变变量类型型 PAGEREF _Toc320528155 h 27 HYPERLINK l _Toc320528156 2-10 IF条条件语句句

4、 PAGEREF _Toc320528156 h 27 HYPERLINK l _Toc320528157 2-11 索单元元自动生生成 PAGEREF _Toc320528157 h 27 HYPERLINK l _Toc320528158 2-12圆圆形隧道道开挖模模拟计算算 PAGEREF _Toc320528158 h 28 HYPERLINK l _Toc3205281559 4-11数组 PAGEREF _Toc320528159 h 30 HYPERLINK l _Toc320528160 4-2函数数操作 PAGEREF _Toc320528160 h 31 HYPERLINK

5、 l _Toc320528161 4-3函数数删除与与重定义义 PAGEREF _Toc320528161 h 32 HYPERLINK l _Toc320528162 4-4字符符串 PAGEREF _Toc320528162 h 32 HYPERLINK l _Toc320528163 4-5马蹄蹄形隧道道网格 PAGEREF _Toc320528163 h 33 HYPERLINK l _Toc320528164 4-6复杂杂形状网网格生成成 PAGEREF _Toc320528164 h 33 HYPERLINK l _Toc320528165 4-7网格格连接 PAGEREF _To

6、c320528165 h 34 HYPERLINK l _Toc320528166 4-8立方方体洞穴穴网格生生成 PAGEREF _Toc320528166 h 35 HYPERLINK l _Toc320528167 4-9球体体洞穴网网格生成成 PAGEREF _Toc320528167 h 36 HYPERLINK l _Toc320528168 4-10应应力边界界 PAGEREF _Toc320528168 h 37 HYPERLINK l _Toc3205281669 4-111改变变应力边边界条件件 PAGEREF _Toc320528169 h 37 HYPERLINK l

7、_Toc320528170 4-12位位移边界界 PAGEREF _Toc320528170 h 39 HYPERLINK l _Toc320528171 4-13不不考虑重重力影响响的均匀匀应力 PAGEREF _Toc320528171 h 41 HYPERLINK l _Toc320528172 4-14考考虑应力力梯度的的均匀材材料 PAGEREF _Toc320528172 h 41 HYPERLINK l _Toc320528173 4-15考考虑应力力梯度的的非均匀匀材料 PAGEREF _Toc320528173 h 42 HYPERLINK l _Toc320528174 4

8、-16非非均匀网网格应力力初始化化 PAGEREF _Toc320528174 h 42 HYPERLINK l _Toc320528175 4-17不不规则自自由面应应力初始始化 PAGEREF _Toc320528175 h 43 HYPERLINK l _Toc320528176 4-18非非均网格格内部压压实 PAGEREF _Toc320528176 h 43 HYPERLINK l _Toc320528177 4-19模模型改变变后初始始应力变变化 PAGEREF _Toc320528177 h 44 HYPERLINK l _Toc320528178 4-20应应力与孔孔隙压力力

9、的初始始化 PAGEREF _Toc320528178 h 44 HYPERLINK l _Toc320528179 4-21加加载顺序序 PAGEREF _Toc320528179 h 45 HYPERLINK l _Toc320528180 6-1 VV级围岩岩施工过过程模拟拟 PAGEREF _Toc320528180 h 47 HYPERLINK l _Toc320528181 6.2 IIV级围围岩施工工过程模模拟 PAGEREF _Toc320528181 h 50 HYPERLINK l _Toc320528182 6.3 IIII级级围岩施施工过程程模拟 PAGEREF _To

10、c320528182 h 52 HYPERLINK l _Toc320528183 第七章 命命令流按按照顺序序进行 PAGEREF _Toc320528183 h 54 HYPERLINK l _Toc320528184 ;右隧道中中地层网网格 PAGEREF _Toc320528184 h 54 HYPERLINK l _Toc320528185 ;右隧道上上侧土体体网格 PAGEREF _Toc320528185 h 54 HYPERLINK l _Toc320528186 ;右隧道下下侧土体体网格 PAGEREF _Toc320528186 h 54 HYPERLINK l _Toc3

11、20528187 第七章 命命令流按按照顺序序进行 PAGEREF _Toc320528187 h 571.1 FFLACC常识1. FLLAC33D是有有限元程程序吗？答：不是！是有限限差分法法。2. 最先先需要掌掌握的命命令有哪哪些？答：需要掌掌握geen, inii, aapp, pllo, sollve等等建模、初初始条件件、边界界条件、后后处理和和求解的的命令。3. 怎样样看模型型的样子子？答：ploo bllo ggro可可以看到到不同的的grooup的的颜色分分布4. 怎样样看模型型的边界界情况？答：ploo gppfixx reed5. 怎样样看模型型的体力力分布？答：ploo

12、 faap rred6. 怎样样看模型型的云图图？答：位移：ploo coon ddis (xddis, yddis, zddis) 应应力：pplo conn szz (ssy, sx, sxxy, syzz, ssxz)7. 怎样样看模型型的矢量量图？答：ploo diis (xdiis, ydiis, zdiis)8. 怎样样看模型型有多少少单元、节节点？答：prii innfo9. 怎样样输出模模型的后后处理图图？答：Fille/PPrinnt ttypee/Jppg ffilee，然后后选择FFilee/Prrintt，将保保存格式式选择为为jpee文件10. 怎怎样调用用一个文文件

13、？答：Fille/ccalll或者ccalll命令10. 如如何施加加面力？答：appp nsstreess11. 如如何调整整视图的的大小、角角度？答：综合使使用x, y, z, m, Shhiftt键，配配合使用用Ctrrl+RR，Cttrl+Z等快快捷键12. 如如何进行行边界约约束？答：fixx x rann （约约束的是是速度，在在初始情情况下约约束等效效于位移移约束）13. 如如何知道道每个单单元的IID？答：用鼠标标双击单单元的表表面，可可以知道道单元的的ID和和坐标14. 如如何进行行切片？答：ploo seet pplanne oori (点坐坐标) norrm (法向矢矢量

14、) pplo conn szz pllanee (显显示z方方向应力力的切片片) 15. 如如何保存存计算结结果？答：savve +文件名名.16. 如如何调用用已保存存的结果果？答：resst +文件名名；或者者Fille / Reestoore17. 如如何暂停停计算？答：Escc18. 如如何在程程序中进进行暂停停，并可可恢复计计算？答：在命令令中加入入pauuse命命令，用用conntinnue进进行继续续19. 如如何跳过过某个计计算步？答：在计算算中按空空格键跳跳过本次次计算，自自动进入入下一步步20. FFishh是什么么东西？答：是FLLAC33D的内内置语言言，可以以用来进进

15、行参数数化模型型、完成成命令本本身不能能进行的的功能21. FFishh是否一一定要学学？答：可以不不用，需需要的时时候查MMannnuall获得需需要的变变量就可可以了22. FFLACC3D允允许的命命令文件件格式有有哪些？答：无所谓谓，只要要是文本本文件，什什么后缀缀都可以以23. 如如何调用用一些可可选模块块？答：connfigg dyyn (fluuid, crreepp, ccppuudm)24 .如如何在圆圆柱体四四周如何何施加约约束条件件？答：可以用用fixx rran cyllindder endd1 eend22 raadiuus rr1 ccyliindeer eend1

16、1 ennd2 raddiuss r22 noot，其其中r2rr1，其其实就是是选择rrangge的时时候选两两个圆柱柱的差，即即得到边边界。25.如何何能把一一个PLLOT的的图像数数据导出出来以便便用其他他软件绘绘图？答：用seet llog on 命令，把把数据导导出来，转转到exxcell里处理理一下，然然后用ssurffer或或者什么么作图软软件绘制就行了。26.关于于dd、ddip的的理解：答： 首先先你要确确定你定定义面的的法线方方向，ddip 是倾角角的意思思（xooy面和和定义面面的夹角角），法法线向上上dipp是锐角角，法线线向下ddip 时钝角角。ddd是y轴轴顺时针针

17、转到法法线在xxoy面面上的投投影线时时，所转转过的角角度。27.如何何调用FFISHH文件？答：新手调调用FIISH文文件出错错的原因因：是：文件名名为 ffishhcalll.ffis.txtt(因为为在记事事本里这这个txxt的后后缀是默认的的）。修修改方法法是：将将 fiishccalll.fiis.ttxt里里面的内内容全部部复制到到新的记记事本里里，再按按保存，保存时候候将保存存类型改改为所有有文件(*,*),这这时你便便可在文文件名栏栏内输入入 fiishccalll.fiis，类类型也就就正确了。28.如何何理解TTenssionn-p sheear-p ttenssionn-

18、n sheear-n表示示的含义义？答：在计算算的循环环里面,每个循循环中,每个zzonee都依据据faiilurre ccritteriion处处于不同同的状态态,shhearr 和tensiion分分别表示示因受剪剪和受拉拉而处于于塑性状状态.nn表示nnow,p表示示preevioous,即分别别表示在在现在和和以前的循环当中中处于塑塑性状态态.29.绘制制相关曲曲线命令令的一些些解释在plott hiist m vvs nn的形式式里,mm代表yy轴,nn代表xx轴(不不管m,n的正正负);-表示示对其值值作mmirrror,如果果大家使使用AUUTOCCAD的的话就明明白mmirrr

19、or的含义义了30.如何何查看剖剖面上任任一点的的位移值值？答：priint gp dissp rrangge xx() y( ) zz( )或printt gpp diisp rannge id plot gp dissp rrangge xx() y( ) zz( )或plot gp dissp rrangge iid31. aa: ffix x yy z rannge x -0.11 0.1 bb: ffix x rrangge xx -00.1 0.11 请问问a式与与b式有有何区别别？答：a 表表示 在在x=00 这个个平面 是固定定的；bb 表示示在x=0 这这个平面面 只xx方向

20、不不能发生生运动32.innitiial 与 aapplly 有有何区别别？答：初始化化命令，如如初始化化计算体体的应力力状态等等；applyy边界条条件限制制命令，如如施加边边界的力力、位移移等约束束等。initiial的的应力状状态会随随计算过过程的发发生而发发生改变变，一般般体力需需要初始始化，而而appply施施加的边边界条件件不会发生变变化。33.如何何查看各各个时段段不平衡衡力的具具体数值值？答：采用hhis来来记录计计算，包包括位移移应力等等 命令令hiss unnball his ggp(zzonee) zzdiss raangee (00 0 0) 或者iid=？导出数据命命

21、令hiis wwritte nn vss m beggin 时步eend 时步 +fiile fillenaame.hiss要导出不平平衡力的的具体数数值his uunbaalstep 10000000 orr soolveehis wwritte 11 vss sttep beggin 1 eend 10000 ffilee 1223.hhis使用上述命命令就可可以查看看各个时时步下的的不平衡衡力的具具体数值值34.如何何显示变变形轮廓廓线的命命令？ plo sske maggf 110 其中中10为为放大系系数35、怎样样开始入入门？答：初步学学习首先先看说明明书 UUserrs gurr

22、id的的Secctioon2 Getttinng sstarrtedd和Seectiion 3 PProbblemmsloviing witth FFLACC36、哪些些命令需需要掌握握？答：griid ggen proop iinteerfaace inii fiix sset hisst sstepp soolvee paausee coont savv 等，在在学习的的过程中遇到这些些命令和和不会的的命令应应该查看看说明书书的Coommaand Reffereencee，里面面讲的很很清楚，特特别需要要注意的是一定定要仔细细看一下下Commmannd RRefeerennce开开头的OO

23、ptiionaal RRangge PPhraasess，这是是讲范围围的，因因为很多命令的的施加都都需要加加到模型型的一部部分，这这里就告告诉你怎怎样来表表示这一一部分。这这一原则则同样适适合与其他5个IITASSCA软软件。还还有一点点就是六六个软件件（FLLAC FALLC3DD UDDEC 3DEEC PPFC PFCC3D）的的rannge大大不相37、FLLAC和和FLAAC3DD建模有有什么区区别？答：FLAAC和FFLACC3D基基本的建建模思想想大致是是一样的的，不同同的是FFLACC是首先先用GRRID命命令把整整体的网网格画出出来，然后用用MODDELLL NUULL命命

24、令把网网格分成成许多部部分，再再用GEEN命令令把这些些部分一一块一块块搭接起起来。然然后再用GENN LIINE，GGEN ARCC等命令令画出要要开挖的的边界，在在FLAAC中产产生网格格的命令令只有一一个GRRID，而而且只能用一次。而而FLAAC3DD是直接接用GEEN 命命令一块块一块搭搭接，GGEN可可以用很很多次。38、怎样样查看模模型？答：ploot ggridd 可以以查看网网格，pplott grrid numm 可以以查看节节点号。39、怎样样查看边边界条件件？答：ploot ffix40、怎样样查看预预先留设设的开挖挖边界？答：ploot ggridd maark。41

25、、怎么么查看模模型的云云图？答：pl syyy (ssxx xdiis yydiss siig1.) filll ，注注意，带带filll关键键字就是是查看云云图，不不带fiill 关键字字就是查看等等直线图图。42、怎样样查看模模型的矢矢量图？答：pl griid vvel (diisp) reed43、怎样样查看IINTEERFAACE？答：pl griid iifacce iid rred printt iffacee idd ; id 是innterrfacce的编编号44、FLLAC怎怎么调整整视图？答：ctrrl+zz 或sshifft+zz ,然然后点击击要放大大部位的的左上角角

26、和右下下角就可可以放大大某个部部位，如如果要缩缩小，双击屏幕。45、判断断模型中中网格点点是一个个还是两两个的方方法 显示示模型网网格点的的ID，如如果有两两个IDD号，说说明是两两个。节节理面上上的网格格点应该该有两个个ID号号。46、改变变模型应应力边界界时FLLAC与与3DEEC的不不同之处处在3DECC中，当当改变应应力边界界时，模模型表面面的应力力直接被被更换，FFLACC中改变变应力边边界时，模模型表面的应力力为两者者应力相相加。47、fllac中中影响运运行时间间的因素素（1）单元元材料，结结构面之之间的刚刚度对比比很大；（2）划分分单元尺尺寸对比比很大，尽尽量避免免长细比比大于

27、55：1的的细长单单元。48 FLLAC中中有多少少中材料料模型内置11中中本构模模型%_4 W/NN,n R5qq（1）零模模型 用用于开挖挖回填的的区域（2）各向向同性弹弹性模型型（3）正交交各向异异性弹性性模型（4）横观观各向同同性弹性性模型（5）德鲁鲁克普普拉格塑塑性模型型（6）莫尔尔库仑仑塑性模模型 一一般岩土土体的力力学行为为（7）节理理化塑性性模型 用于研研究材料料破坏后后行为（如如边破稳稳定，屈屈服煤矿矿开采）（8）应变变硬化/软化模模型库库仑塑性性模型 薄薄板层状状材料破破坏后力力学行为为（9）双线线性应变变硬化/软化节节理化塑塑性模型型 临近近层状土土层开挖挖（10）双双屈

28、服塑塑性模型型（11）修修正的剑剑桥粘土土模型49、对不不平衡力力的理解解flac中中用最大大不平衡衡力和典典型内力力大的比比值R表表示模型型的不平平衡力的的相对大大小，该该值为百百分数，且且从不会减小到到零，因因此，模模型不会会达到绝绝对稳定定的平横横状态，可可以根据据要求的的精度设设定R值值，模型型默认的R值是11e-55。注意R值小小仅表示示所有网网格点的的力都平平衡，但但模型还还可能发发生塑性性流动，通通过 pplott veel 命命令查看看网格点速度度，一个个低振幅幅随即速速度场是是没有塑塑性流动动的平衡衡状态的的可靠的的指示器器。50、简单单命令解解释：plot addd coo

29、nt dissp oout on 时时显显示云图图plot bloock grooup gpnnum on 显示节节点编号号plot intter redd skke 看接接触面的的效果.del rrangge ggrouup 11 noot 不不删除ggrouup信息息 del rrangge ggrouup 22 删除除grooup信信息ini zz addd 数数值大小小 raangee grroupp naame 顺z方方向移动动建立接触面面 innterrfacce 11 faace rannge x(范范围）yy(范围围） zz(范围围） 或或者 iinteerfaace 1 ff

30、acee raangee pllanee orrig(x yy z 坐标） norrm（法法线方向向）waterr taablee cllearr 清除除水位线线pl seet bbackk whhitee ;设设置背景景set pplott biitmaap ssizee (8800,5100) ;设置图图片大小小 plot sett caaptiion sizze 337 ; 设置置图例数数字大小小 (ddefaaultt=355 范围围10-50)plot sett ceenteer 1100 0 440 fix xx y z rrangge xx -00.1 0.11 表示示在x=0

31、这这个平面面是固定定的 ffixeed bbounndarryfix xx raangee x -0.1 00.1表表示在xx=0 这个平平面只xx方向不不能发生生运动rrolller bouundaary51、fiish命命令的解解释：def aaux11 定义义函数loop n(11,npptabb1) 开始一一个循环环，从11循环到到变量nnptaab1rr=xttablle(11,n) 从ttablle 11 的xx列中读读第n个个数付给给变量rrrzz=yttablle(11,n) 从ttablle 11 的yy列中读读第n个个数付给给变量zzzloop k(11,npproff)

32、开开始第二二个循环环，从11循环到到变量nnproofcase_of k 下下面是一一个选择择结构，利利用变量量k判断断i F aalfaa=0.5*(pi+.1)*flloatt(k-2)/flooat(nprrof-2)xtablle(nn+100,k)=rrr*coos(aalfaa)ytablle(nn+100,k)=rrr*siin(aalfaa)case 1 当当k11的时候候执行下下面语句句xtablle(nn+100,k)=rrr taablee n+10的的x列第第k个值值等于变变量rrrytablle(nn+100,k)=-440.00 taablee n+10的的y列第第

33、k个值值等于变变量-440case 2 当当k22的时候候执行下下面语句句xtablle(nn+100,k)=rrr taablee n+10的的x列第第k个值值等于变变量rrrytablle(nn+100,k)=0.0 ttablle nn+100的y列列第k个个值等于于变量00end_ccasee 结束束选择结结构end_lloopp 结束第第二个循循环end_lloopp 结束束第一个个循环end 函函数定义义结束aux1 执行什什么定义义的函数数52. 在在固定面面上施加加位移时时，fiix与aapplly不能能同时使使用，应应该用iini和和fixx,才能能施加上上。53.菜单单 p

34、llot iteems/addd/sttruccturre eelemmentts/ggeommetrry/结结构单元元的形式式，把对对话框里里的maag ffac的的值改成成不同的的数，就就可以看看到位移移放大后后的结构构单元。54.你可可以用键键盘上的的上、下下、左、右右肩在窗窗口中移移动模型型。55.如何何能把一一个PLLOT的的图像数数据导出出来以便便用其他他软件绘绘图？答：用seet llog on 命令，把把数据导导出来，转转到exxcell里处理理一下，然然后用ssurffer或或者什么么作图软软件绘制制就行了了。56 在进进行三轴轴压缩时时候，如如何施加加围压？答：使用aapp

35、lly ssxx =-11e5 szzz=-11e5 rannge cyll ennd1 0,00,0 endd2 00,2,0 rradiius 1，其其中X和和Z轴是是在我们们眼睛正正对的面面。57 在进进行三轴轴压缩时时候，如如何对轴轴向力进进行位移移控制？答：使用 inii yvvel 1e-7 rrangge yy -.1 .1 iini yveel 11e-77 raangee y 1.99 2.1 即即可，数数值可以以变化58 在进进行边坡坡计算的的时候，为为什么用用soll, ffos 命令运运行不行行？答：应该输输入sool ffos sloope33dfoos.ssav a

36、sssociiateed(不不要也行行） 再再用reestoore 命令调调用即可可。59 在FFLACC3D计计算中，如如何记录录程序中中的值，如如果不想想时时记记录该怎怎么弄？ 答：用 hhistt wttitee n11 vss n22 fiile.txtt， 想想间隔记记录可用用sett hiist_repp n。60 想记记录接触触面上的的相关力力学参数数该怎么么办？答：hisst iinteer ii nsstreess /nddis/ssttresss/ssdissp61 接触触面设置置参数包包括哪几几个？答：cohh diial friic kkn kks tten.62 如何

37、何输出单单元格的的数据？答；exppgriid .flaac3dd63 如何何在剖面面的沉降降云图上上增加网网格、剖剖面上的的矢量及及坐标值值？ DD Q Z9ss1 e4YY U%q 答：ploo addd sske plo aadd diss pllaneeplo aadd axee64 如何何做切片片？ 答答：pllot sett pllanee orri xx y z nnormm x y zz65 举一一个隧道道支护开开挖的例例子？gen zzonee raadbrrickk p00 (00,0,0) p1 (100,0,0) p2 (0,10,0) p3 (0,0,110) &si

38、ze 3 55 5 7 rrat 1 11 1 1.55 diim 11 4 2 ffilllmodell moohrprop bullk 11e8 sheear .3ee8 ffricc 355prop cohh 1ee3 ttenss 1ee3fix xx raangee x -0.1 00.1fix zz raangee z -0.1 00.1fix yy raangee y 9.99 100.1set ggravv 0 10 0ini ddenss 10000 ini ssyy=0.00 grrad 0,-100000,0ini ssxx=0.00 grrad 0,-50000,00i

39、ni sszz=0.00 grrad 0,-50000,00applyy sxxx=00.0 graad 00,-550000,0 rannge x 99.9,10.1applyy szzz=00.0 graad 00,-550000,0 rannge z 99.9,10.1set llarggehist unbbalsolveesave tr_eq.savv; monnitoor ttrennch exccavaatioonini xxdiss=0 ydiis=00 zddis=0hist gp xdiisp 1,00,0hist gp zdiisp 0,00,2 ; exccavaatio

40、on sstepps rangee naame excc1 xx 0 1 yy 0.0 00.8 z 00 2 rangee naame excc2 xx 0 1 yy 0.8 11.6 z 00 2rangee naame excc3 xx 0 1 yy 1.6 22.4 z 00 2rangee naame excc4 xx 0 1 yy 2.4 33.2 z 00 2rangee naame excc5 xx 0 1 yy 3.2 44.0 z 00 2; exccavaatioon sstepp 1modell nuull rannge excc1step 100000save exc

41、c1.ssav; exccavaatioon sstepp 2modell nuull rannge excc2step 10000save excc2.ssav; exccavaatioon sstepp 3modell nuull rannge excc3step 10000save excc3.ssav ; exccavaatioon sstepp 4modell nuull rannge excc4step 10000save excc4.ssav; exccavaatioon sstepp 5modell nuull rannge excc5step 10000 save excc5

42、.ssav; insstalll ccablle ssuppporttrest excc2.ssavsel ccablle bbeg 1.00,0.4,11.5 endd 5.0,00.4,1.55 nsseg 4sel ccablle bbeg 1.00,0.4,00.5 endd 5.0,00.4,0.55 nsseg 4sel ccablle bbeg 1.00,1.2,11.5 endd 5.0,11.2,1.55 nsseg 4sel ccablle bbeg 1.00,1.2,00.5 endd 5.0,11.2,0.55 nsseg 4sel ccablle ppropp emm

43、od 2e99 yttenssionn 1ee8 xxcarrea 1.00 grr_cooh 11e100 grr_k 2e99 grr_peer 11.0modell nuull rannge excc3step 10000save cabb3.ssavsel ccablle bbeg 1.00,2.0,11.5 endd 5.0,22.0,1.55 nsseg 4sel ccablle bbeg 1.00,2.0,00.5 endd 5.0,22.0,0.55 nsseg 4sel ccablle ppropp emmod 2e99 yttenssionn 1ee8 xxcarrea 1

44、.00 grr_cooh 11e100 grr_k 2e99 grr_peer 11.0modell nuull rannge excc4;step 10000save cabb4.ssavsel ccablle bbeg 1.00,2.8,11.5 endd 5.0,22.8,1.55 nsseg 4sel ccablle bbeg 1.00,2.8,00.5 endd 5.0,22.8,0.55 nsseg 4sel ccablle ppropp emmod 2e99 yttenssionn 1ee8 xxcarrea 1.00 grr_cooh 11e100 grr_k 2e99 grr

45、_peer 11.0modell nuull rannge excc5step 10000save cabb5.ssavplot creeatee TrrencchSuuppoort plot sett root 663.889, 3344.355, 33.4plot sett ceent 2.8893, 3.4466, 33.8999plot sett siize 4.1136plot addd coont dissp oout onplot addd seel ccablle fforcceplot addd axxes blaack66 举一一个三轴轴压缩的的例子？gen zzonee

46、cyyl pp0 00 0 0 pp1 11 0 0 pp2 00 2 0 pp3 00 0 1 ssizee 4 5 44 gen zzonee reefleect norrm 11,0,0gen zzonee reefleect norrm 00,0,19modell moohr prop bullk 11.199e100 shhearr 1.1e110prop cohh 2.72ee5 ffricc 444 teen 22e5fix xx y z rrangge yy -.1 .1fix xx y z rrangge yy 1.9 22.1applyy sxxx =-1ee5 sszz=

47、-1ee5 rrangge ccyl endd1 00,0,0 eend22 0,2,00 raadiuus 11ini yyvell 0 graad 00 -11e-44 0 rannge y 00 2sol1.FLAAC程式式使用前前准备步步骤步骤1：依依比例画画出所欲欲分析的的资料在纸上画出出地点的的位置、地地层资料料、并简简单标示示距离及及深度资资料。步骤2：换换算输入入资料成成同一单单位将现有地层层资料，如如 Deensiity, Buulk modduluus, Youungs mmoduuluss, ttenssionn, ccoheesioon, friictiion Angg

48、le 等资料料，换算算成同一一单 位位。附注 ：需需谨慎检检查输入入资料的的单位，如如因单位位不同而而造成过过大或过过小的值值，将会会造成 FLAAC无法法计算，而而产生EERROOR信息息。 步骤3：应应用公式式简略计计算应用公式或或依据经经验，简简略算出出FLAAC输出出资料的的范围，以以作为DDebuug及输输出资料料分析时时验证。步骤4：建建立x, y坐坐标与nnodee i, j之之间的关关系在图上距离离及深度度的关系系，建立立x, y坐标标系统，再再由x, y坐坐标系统统，转换换与网格格间系，为为了便于于以后输输出资料料的分析析，故应应确实 掌握网网格的位位置及其其相对应应的x,

49、y 坐标。建议在敏感感区域使使用较密密的网格格，其它它地方则则使用较较疏的网网格，刚刚开始执执行程序序时，不不宜使用用网格太太大的网网格数目目，因尽尽量使网网 格总总数少于于10000，以以节省时时间。2.FLAAC 33D 程程序的编编写步骤骤：1 Connfigg _2 Griid _3 Moddel _4 求起始始的应力力平衡_ (1)建建立x, y坐坐标与网网格的关关系，建建议使用用Genn指示：Gen xx1,yy1 xx2,yy2 ,x3,y3 x4,y4 i=ii0,ii1 jj=j00,j11 详细细指令参参见使用用手册，FFLACC程式可可自动产产生x, y 坐标与与网格的的

50、关系，但但由 于于产生的的网格坐坐标不易易控制，将将对其它它的工作作产生负负面影响响，故依依使用前前步骤44所建立立的关系系，将网网格依其其疏密程程度需要要的不同同，实 际控制制网格的的坐标。 (2)设设定材料料性质：proop (3)设设定外力力：Seet GGravv, AApplly PPresssurre, inii sxxx, Syyy (4)设设定边界界条件：fixx, ffreee (5)求求起始的的应力平平衡：ssolvve (6)储储存：SSavee5 求工程程的影响响_求出区域内内的应力力分布情情况后，再再依工程程的流程程及步骤骤阶段执执行各工工程进行行过程的的影响，建建议

51、使用用以下的的步骤： (1)调调出起初初的应力力平衡：re_ .savv (2)设设定新的的材料性性质：mmodeel,ppropp (3)设设定新的的支撑性性质：sstruuct (4)设设定新的的外力 (5)设设定边界界条件 (6)求求工程时时的应力力平衡(7)储存存1.2常用用命令流流;模型镜像像gen zzonee raadcyylinnderr siize 25 1 225 225gen zzonee reefleect norrmall -11 0 0 ooriggin x yy z(面上一一点）;沿X轴轴镜像,通过对对称平面面法线向向量确定定对称面面gen zzonee reef

52、leect norrmall 0 0 -1 ;沿z轴轴镜像;绘图控制制pl coontoour szzz ouutliine on;在模型中中显示位位移-应应变曲线线hist gp ydiisp 0,00,0hist zonne ssyy 0,11,0hist zonne ssyy 1,11,0pl hiis -2 -3 vvs 11 ;在在ploot hhistt m vs n的形形式里,m代表表y轴,n代表表x轴(不管mm,n的的正负); -表表示对其其值作mirrrorr;对模型进进行压缩缩实验的的方法;即在模型型两侧施施加相反反方向的的速度ini yyvell 1ee-7 rannge

53、 y -.1 .1ini yyvell -11e-77 raangee y 1.99 2.1;修改模型型的坐标标值ini xx addd -1000 y addd -1100 z aadd -1000;显示云图图的同时时也显示示模型网网格轮廓廓plot addd coont dissp ooutllinee onn ;ggraddiennt更精精确;输入角度度、弧度度方法 pi=，900为900.0*deggradddef sset_vallsptA = 255.0 * ssin(pi/2);ptAA=255.0ptB = 255.0 * ccos( 600.0*deggradd );ptBB

54、=122.5ptC=ppi;pptC=3.1141559266endset_vvalssprintt pttA pptB ptCC;施加结构构单元方方法sel sshelll iid=55 raangee cyylinnderr ennd1=(0.0, 0.00,0.0) &end2=(0.0,225.00,0.0) raddiuss=244.5 nottplot addd seel ggeomm bllackk bllackk ciid oon sscalle=00.033sel nnodee innit zpoos aadd -255.0;如何显示示某一平平面plot creeatee n

55、aame_plaaneplot sett pllanee ooriggin 3 44 0 norrmall 1 0 00plot addd coont dissp pplanne bbehiind shaade onplot addd seel ggeomm bllackk plot addd axxes reddplot currrennt nnamee_pllaneeplot shoow;如何施加加法向应应力并显显示所施施加的力力applyy nsstreess -3.0E77 raangee cyylinnderr ennd1 0 -1 00 ennd20010 raddiuss 399

56、.9 nott;柱体体外表面面施加法法向力applyy nsstreess -5.0E66 raangee cyylinnderr ennd1 0 -1 00 ennd20010 raddiuss 2.1;柱柱体内表表面施加加法向力力plot griid ggreeen ffap blaack;显示施施加的力力;获取某网网格点的的ID，坐坐标值def aaaget_iid=ggp_iid(ggp_nnearr(3,0,00)endaaprintt geet_iidprintt gpp poosittionn raangee idd 3;保存命令令流文件件set llog onset llogf

57、filee huua.llog;文件存存放的路路径与调调用的文文本路径径相同set llog offf;fishh函数的的说明与与解释def ffindd_addd ;定义义fissh函数数finnd_aadd head = nnulll ;给hheadd赋值p_gp = ggp_hheadd ;第一个个网格结结点的指指针赋给给p_ggploop_whiile p_ggp # nuull ;当当p_ggp值不不为nuull时时作循环环x_poss = gp_xpoos(pp_gpp) ;将指指针为pp_gpp的结点点的x坐坐标值赋赋给x_possif x_poss = 10.0 tthenn

58、;如果xx_poos = 100.0 则(执执行)new = geet_mmem(2) ;从从主内存存空间里里得到22个fiish变变量对象象并返回回第一个个对象的的开始地地址mem(nnew) = heaad ;将hheadd类型和和数值置置于地址址为neew的ffishh变量mem(nnew+1)=p_ggp ;第二二个变量量存储符符合条件件的zoone地地址head = nnew ;将nnew值值赋给hheaddendiffp_gp = ggp_nnextt(p_gp) ;将结点点指针为为p_ggp的下下一个结结点的指指针赋给给p_ggpendlooop;结束循循环end ;结束束fis

59、sh函数数;实际上这这个fiish函函数为满满足条件件（x坐坐标为110的）的的zonne的地地址开辟辟一定的的地址空空间，各各地址之之间存在在一定的的联系;找到第一一个符合合条件的的zonne地址址后，用用neww = gett_meem(22) 从从主内存存空间里里得到22个fiish变变量对象象并返回回第一个个对象的的开始地地址，并并用;下下面的mmem(neww) = heead ，将将第一次次开辟的的两个变变量的第第一个变变量存储储地址hheadd（注意意第一次次heaadnnulll），第第二个变变量存储储第一;个符合合条件的的zonne地址址，并将将第一个个变量的的地址赋赋予he

60、ead（hheaddneew），第第一次循循环结束束；下次次循环，同同样开辟辟两个变变量对象象，第;一个变变量对象象记录上上次循环环开辟的的第一个个变量的的地址，第第二个记记录第二二个符合合条件的的zonne地址址，其余余循环依依次类推推，这样样子就建建立了;一个符符合条件件的zoone地地址链条条，方便便以后使使用。；FLACC后处理理；输出JPPG图片片set pplott jppgset pplott quualiity 1000plot harrd ffilee XXX.jppg;输出TXXT文件件hist wriite X ffilee xxx.txxt ；(xx是历史史记录IID）