Flac学习心得

Flac&3dec学习心得

? ? ? ? ? ? ? ? ?

快捷键—显示位移等值线：pl cont disp 快捷键—显示塑性区：pl bl st 快捷键—显示网格：pl grid

显示某根锚杆的应力：pl sel cable stress range id xx 在显示等直线等时可以用modify中的outline显示单元的边界 平移模型的命令：ini x add 30 (x方向平移30m)

在三维计算中，想看某一剖面的图，可以用edit菜单中的cutting plane操作

任何操作如果不想对某个区域起作用，可以用not关键词，比如gr 5 not gr 6 not 改善锚杆作用的一种操作—把锚杆单元的link删掉，再重新见建一个link与周围的单元连接起来，如下所示：

sel del link range id 1

;sel del link range id 11 删掉link11 ;sel del link range id 12 ;sel del link range id 22

;sel link 2 11 target zone 在原来的link 处（此处的2是该处的node号）新建

link2 与单元连接

;sel link net rang id 11 ;sel link net rang id 12 ;sel link 13 22 target zone

;sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 11 any id 1 any id 12 any id 22

三种方式模拟预应力锚索：

方式1、通过删除－建立link模拟。

删除掉锚索端头原来CABLE和ZONE之间自动建立的连接，然后建立钢性连接。 示例命令如下：

sel 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 id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 & gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,10

sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 & gr_coh 10e5 gr_k 2e7 range cid 11,17 sel delete link range id 1 sel link id=100 1 target zone

sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 100 sel cable id=1 pretension 60e3 range cid 1,10

方式2、通过赋三段属性模拟。即将锚索的端头、自由段、锚固段赋不同的属性来模拟预应力锚索。示例命令如下： sel 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 & gr_coh 1 gr_k 1 gr_per 0.0785 range cid 2,10

sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 & gr_coh 10e5 gr_k 2e7 range cid 11,17

?

sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 & gr_coh 10e8 gr_k 2e10 range cid 1,1

sel cable id=1 pretension 60e3 range cid 1,10

方式3、借助别的结构单元（如liner单元）来模拟托盘。删除掉锚索端头的link，然后建立新的link，新link的target为liner上的node。示例命令如下：

sel 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 id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 & gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,10

sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 & gr_coh 10e5 gr_k 2e7 range cid 11,17 sel liner range y=-.1, .1 x=-1,1 z=-1,1

sel liner PROP iso=( 25e9, 0.15) thick=0.1 ; concrete sel liner PROP cs_nk=8e8 cs_sk=8e8 &

cs_ncut=0.0 cs_scoh=0.0 cs_scohres=0.0 cs_sfric=0.0 sel delete link range id 1

sel link id=100 1 target node tgt_num 18

sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 100 sel cable id=1 pretension 60e3 range cid 1,10

最后总结

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

2、预应力通过sel cable id=1 pretension 60e3 range cid 1,10方式施加 3、预应力加在自由 4、全锚预应力加全根。

? 在两个节点之间建立link的命令（特别适用于在锚杆和衬砌之间建立link）

sel set link node_tol 0.01 ;------0.01是误差范围，即相距0.01m之内才能建立

sel link xxx（id） 2 target=node tgt_num 47 ;--------xxx是建立的link编号，2是要建立link的节点，47是被连接的节点

sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid

? 用切片显示剖面上的等值线图等时，不要切在节点位置，而要避开，否则会出现显示上

的一些杂乱信息。

? attach face：在模型建好，调用后，使用attach face命令可以把那些不连续但节点坐标

都相同的单元的面粘贴在一起，使节点合并。但是对于节点坐标不同的面则无效。 ? 加法向面力的查看方法：在某个面上加了面力后，如果想查看，可以使用pl add fap命

令，如果加的比较小，那么需要在modify中将scale调大。注意，查看前需要先用step 0命令，使系统更新。同样，查看体力用pl add fob

? 建立interface的巧妙方法：如果模型已经建好，想增加interface的话，把模型的一部分

使用gen separate +group的编号，即可把该group的单元从整个模型中切割出来，然后使用interface xx（编号） wrap a b，其中a和b是要建立interface位置所在的两个group。 ? 用了separate切割后再用attach face粘贴，可以恢复切割前的状态，网格仍然连续。切

割后的模型，即使加了interface，仍然可以使用attach face将网格粘贴起来，这样interface虽然仍存在，但不起作用了。同样，虽然没有切割的网格是无法加interface的，但是切

割后再粘贴起来的网格，仍然可以在原位置加interface，但是不起作用，网格仍然为连续。

? attach face虽然能把网格粘在一起，但是编号不能合并，用gen merge 0.01 range xxx不

仅可以粘贴网格，而且可以合并节点编号，其中0.01是容差，可以自己设定。

3dec 3dec改内存的方式：通过dos启动3dec，并更改内存。具体操作—程序-附件-命令提示符-cd\\-cd\\program files\\itasca\\3dec400-输入“3dec.exe 512”。这样就把内存改为512M。

破解版的3dec里不要使用save命令，否则将导致内存不够用，而且保存的.sav文件也无法restore。

BOUNDARY命令如果在同一边界第二次使用，则第二次的值加在初值上，而不是替换。 对于可变形块体，GENERATE quad命令一般用于建立内部或重点区域的网格，而GENERATE edge则用于建立外部不重要的区域网格。

速率边界（Velocity boundary）必须在应力边界（stress boundary）之后添加，否则速率边界将被取消。

不同prop对应的颜色是不同的。 str prop 1 ...... 对应的是青色。 str prop 2 ...... 对应的是绿色。 str prop 3 ...... 对应的是红色。 高手的程序解释

config cell 配置相互接触关系，允许块体与另一个块体接触后弹起 round 0.001 定义角点圆滑情况，默认0.5。推荐值为典型块体的1％ jset 30 0 50 0 0 0 4 0 range x 30 40 y 92 97 在（30，92）（40，97）范围内生成节理，节理倾角+30，迹长50，贯通，间距4

jset -60 0 50 0 0 0 4 0 range x 30 40 y 92 97 在（30，92）（40，97）范围内生成节理，节理倾角-60，迹长50，贯通，间距4

change mat 1 定义模型本构关系为线弹性

change jmat 1 jcons 2 定义节理本构关系为完全弹性模型， set jmatdf 1 jcondf 2 定义计算开始后的新接触关系 prop m 1 d 0.0025 jkn 5000 jks 5000 jfric 30 各种参数 damp auto 粘性阻尼 damp .23 225 stiff ？？ set ovtol .1 接触叠加范围 pl bl vel iw max 5 pl bl vel hold

3dec中关于结构面参数的确定

前面看到有人说3dec结构面参数比较难以确定，我觉得还是可以确定，并可以和有限元的结果相媲美的，我曾经计算过大型工程的应力结果，其和有限元结果基本一致。

一般而言，真实结构的变形由两方面组成，一是块体本身的变形，二是结构面之间的相互嵌入，也正是因为后者，很多人认为离散元结果人为性太大。其实不是这样的，如果是人为切割的分割面（真实结构中下不存在的），为了避免结构面的嵌入带来的变形干扰正常的计算结果，一般只要保证E/skn约20倍左右即可，此时的结构面互相嵌入对于整个系统变形的影响几乎可以忽略不计。

如果是真实存在的结构面，比如是天然状态下山体的节理裂隙，期间可能充填各种杂质，比如夹泥等等，此时需要使用较为真实的kn参数，一般取为E/T，ks可取kn/(1-v），也可以取与kn一样，这个关系一般不大。

上面的说明文字中，E弹模，大家都知道 s 接触面面积 T 裂隙厚度

设置出图精度： set nx=100 ny=100

在力学分析角度上，离散元对三大定律的满足上与有限元方法不同。从平衡方程上看，离散元采用牛顿第二定律来控制，按围绕各刚性块体形心的力平衡和力矩平衡来满足。从变形协调方程上看，各刚性块体间不再位移连续，而是允许大变形和断裂分开，可以模拟岩体不连续结构面的滑移与开裂。其位移联系为刚性块体间各接触点的相对速度关系。从材料本构关系上看，离散元法避开了复杂的本构关系推导，采用在刚性块体间设置不同种类弹簧和阻尼（法向刚度和阻尼、切向刚度和阻尼）来反映材料的应力位移关系。同时，在求解方法上，离散元法采用动态的松弛方法求解，从而避开了大规模刚度矩阵的形成，不须求解大型联立方程组，能够高效的求解非线性大变形问题。 图形放大的快捷键是“M”，缩小的快捷键是“U”。 单纯jset命令默认为jset d 0 dd 0 org 0 0 0 3dec中建立的模型默认为region 0

delete删掉的部分不能再回填，excavate删掉的可以在用fill命令回填，remove删掉的不可以回填，但是可以显示出来。

显示塑性区：pl plasticity 显示位移：pl disp 显示速率：pl velocity 显示应力 pl syy 显示主应力矢量图：pl principal 显示主应力等值线：pl sint\\smax\\smin 显示单元：pl zone

3dec命令流中每一行长度不宜超过80个字符，否则会出错，但注释长度可以任意 如果因为块体形状或大小影响，单元划分不过去，有两种处理办法：一是在gen edge xx 后加关键词alt ，如gen edge 12.0 alt。或者把不能划分的小块体直接删除，如del block 133395

施加衬砌或喷层的命令流： config liner

struct liner radial_gen 0 0 -.9 0 0 .9 seg 2 8 prop 1 range .43 struct prop 1 fric 60 coh .5e6 kn 1e9 ks 1e9 struct prop 1 nu .15 thick .20 e 15e9 tens .3e6

其中range .43 是给定能够包住洞子临空面范围的半径。如果并不是对全断面施加liner，则还要给定起止平面上的各一个点。

本文来源：https://www.bwwdw.com/article/s1u2.html