flac3d讲座 - 图文

最后总结

1、 预应力和全锚均通过删除－建立LINK模拟。

4 ?; F4 f7 T7 b# {\7 \\2、 预应力通过sel cable id=1 pretension 60e3 range cid 1,10方式施加 3、 预应力加在自由端。4、 全锚预应力加全根

最后想提一点，预应力锚索的几个属性参数非常重要emod ytension xcarea gr_coh gr_k 模拟的时候务必搞清楚几个参数的真正含义！！许多在公开刊物（包括 岩石力学与工程学报）发表的论文，里面的参数都不恰当，有一篇模拟动载对锚杆影响的论文（不方便具体指出），给出的粘结强度竟然为100e5！这是什么概念，每米锚固剂内能承受的最大剪切力为1000吨力！我用这样的锚杆支护井下巷道的话，只怕再大的采动影响，巷道还是稳稳的不动。另外还有一点，另一个参数gr_fric，锚固剂的摩擦角，指定锚固剂内所能承受的最大剪力随着围压的增加而增加，个人觉得没有给出的必要。我查过许多文献，目前还没有查到锚固剂剪切力会随围压增加而增大的论证。也许是本人水平有限，请高手不吝指教。谢谢。

模型很简单，就是一块长方体，不设任何应力，只算附加应力。如下：

C) w5 }3 G8 S/ }* Y! K. x9 A2 ^) y& k [\

gen zone radtun p0 0,0,0 p1 25,0,0 p2 0,50, 0 p3 0,0,25 size 4 25 4 10 dim 4 4 4 4 ratio 1 1 1 1.1 fill

* R, p, C+ f* ]6 y, |3 cgen zone reflect normal 1 0 0 ori 0 0 0 mo mohr

% G2 U8 U' \\\

gen zone reflect normal 0 0 1 ori 0 0 0

pro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5 ini dens 2000

% h# [* ^( k7 l( y2 R/ fix x range x -25.1,-24.9 fix x range x 24.9 25.1

+ h3 l# {3 Z& w K+ j+

fix y range y 49.9 50.1 fix z range z -25.1 -24.9

/ e, X3 U( H- U& }. F6 O4 [

fix z range z 24.9 25.1

( C5 a! Z( t) d0 ?7 u3 K) ` o9 \\' t' n9 u7 [: Vsel cable id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10 sel cable id=1 beg 0,29,0 end 0,35,0 nseg 6

sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 &

9 L) d) A: R* `3 | gr_coh 1 gr_k 1 gr_per 0.0785 range cid 2,10

( P5 {* Z# N0 Z0 W7 s

sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 &

2 m/ N\9 O* e gr_coh 10e5 gr_k 2e7 range cid 11,17

) K0 p\ h

sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 &

! h0 Y0 n: p. D\ gr_coh 10e8 gr_k 2e10 range cid 1,1

5 S1 y0 B; j4 w+ T5 F6 w\& W3 A' ~% Q2 q H. p

\9 D( u# y- F% D$ ~sel cable id=1 pretension 60e3 range cid 1,10 step 2000

前面说的预应力加在自由段合适，现在有疑问：实际工况中，在对锚杆进行拉拔加载预应力时，锚固剂肯定是受力的（岩体通过锚固剂给锚杆剪力以平衡拉拔力），那么在模拟的时候，预应力是不是不能只加载到自由段，也要加载到锚固段？如果锚固段也加的话，肯定不是全锚固段都加，而是加靠近自由段的一部分，那么，应该加多长合适？？

其实如果只模拟巷道支护，不需要考虑这么细，加在自由段就行了。但如果要分析预应力锚杆锚固机理，上面的问题就必须要考虑了。

CAD TO ANSYS TO FLAC3D边坡稳定性分析全程揭密

最近看到不少朋友在问这个方面的问题，我现在把我过去几个月做的工程实例贡献出来，与大家共享。

（1）从CAD图获取等高线点的数据

5 V\

这个过程我是采用捕捉点获取数据的方式，不过是用软件实现的，否则手工操作十分费事。这样的软件在论坛上的ANSYS版块已经有了。有时候，CAD图也并非全有三维信息，不过只要知道二维信息，辅之以高程，就可以了，这个工作量不大。将得到的点的坐标导入到EXCEL中进行处理，转换为下面的格式

k, 1, 600 ,0, 0

\5 z, v

k, 2, 600 ,0, 700 k, 3, 0 ,0, 700 k, 4, 0 ,0, 0

: r' f3 e1 o! N/ d+ N9 _k, 5, 600 ,24, 0 k, 6, 600 ,24, 700 k, 7, 0 ,24, 700

- T+ E+ a( ]% v, nk, 8, 0 ,24, 0

9 `( |) E/ q: ?( G: D即关键点坐标，然后就可以由点生成线了。至于为什么要导入到EXCEL，因为得出的数据点格式未必是符合ANSYS的数据格式，采用她可以进行批量处理，减少工作量。处理完毕后，这些数据点格式可以直接拷贝到记事本中，就成为符合ANSYS输入要求的数据点格式了。

( {6 q- g/ w2 n0 z6 b（2）有了数据点，就进入到下一步，即由点在ANSYS中建立几何模型。

通过数据点生成几何模型的方法多种多样，因人而异，没有统一的方法。不过对于我来说，个人觉得岩土工程模型比较复杂，最好是采用点生成线、线生面，面生成体的办法。至于地质界面，我是采用先生成面，然后切割体的方式生成。附件是我的ANSYS命令流，里面配有说明，大家可以参考一下。附图是按上述步骤得到的结果图和网格模型图，说明一下，这些图和我现在贴出的命令流得到的结果有些差异，原因是我后来更改了模型，不过步骤完全一样；后面我提到的在ANSYS生成水面的方式也是采用面切割体的方式得到的。

请教：如何将table中存放的数据输出成txt文本文件？

在FLAC3D的后处理中，用fish计算一些数据存放在table中，这写数在FLAC命令窗中可以看见，但是怎样才能把这些数输出到文本文件中呢？用hist write nhis file filename 可以输出文件，但是我只是单纯的计算一些数，这些数据是fish中自定义的，没有办法用hist记录。希望高人可以帮忙解答一下，谢谢！ set log on print table n

8 |! c. f' \\0 G会生成一个log文件你打开后里面有table的数据不知道对不对 或者用open, write,close的文件操作写到纪事本里面

问：通过什么命令可以看到FLAC3D任意剖面的位移变形图？

plot set plane norm 1 0 0 origin 5 0 0plot con zdisp out on shade on plane

plot set plane ori (切面上任意一点坐标) norm （该切面法向向量） plot con zdis plane (behind)

* e$ i( O8 {4 d6 U: A( Z! v9 Y3 Y7 |+ v* N4 ]2 s

；behind 省略现实的是一个切片上的内容

edit/cutting plane/ori (输入切面上任意一分坐标，即或x或y或z) and norm （沿哪个方向切就把哪个设为1，即或x或y或z为1，然后其他设为0）