1、Flacsel 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
2、=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 10sel cable id=1 beg 0,29,0 end 0,35,0 nseg 6sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 concret
3、e sel liner PROP cs_nk=8e8 cs_sk=8e8 -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：在模型
4、建好，调用后，使用 attach face 命令可以把那些不连续但节点坐标都相同的单元的面粘贴在一起，使节点合并。但是对于节点坐标不同的面则无效。 加法向面力的查看方法：在某个面上加了面力后，如果想查看，可以使用 pl add fap 命令，如果加的比较小，那么需要在 modify 中将 scale 调大。注意，查看前需要先用step 0 命令，使系统更新。同样，查看体力用 pl add fob 建立 interface 的巧妙方法：如果模型已经建好，想增加 interface 的话，把模型的一部分使用 gen separate +group 的编号，即可把该 group 的单元从整个模型中切
5、割出来，然后使用 interface xx（编号） wrap a b，其中 a 和 b 是要建立 interface 位置所在的两个group。 用了 separate 切割后再用 attach face 粘贴，可以恢复切割前的状态，网格仍然连续。切割后的模型，即使加了 interface，仍然可以使用 attach face 将网格粘贴起来，这样interface 虽然仍存在，但不起作用了。同样，虽然没有切割的网格是无法加 interface的，但是切割后再粘贴起来的网格，仍然可以在原位置加 interface，但是不起作用，网格仍然为连续。 attach face 虽然能把网格粘在一起，但
6、是编号不能合并，用 gen merge 0.01 range xxx 不仅可以粘贴网格，而且可以合并节点编号，其中 0.01 是容差，可以自己设定。3dec 3dec 改内存的方式：通过 dos 启动 3dec，并更改内存。具体操作程序-附件-命令提示符-cd-cdprogram filesitasca3dec400-输入“3dec.exe 512”。这样就把内存改为 512M。破解版的 3dec 里不要使用 save 命令，否则将导致内存不够用，而且保存的.sav 文件也无法 restore。BOUNDARY 命令如果在同一边界第二次使用，则第二次的值加在初值上，而不是替换。对于可变形块体，
7、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。推荐值为典型块体的 1jset 30 0 50 0 0 0
8、4 0 range x 30 40 y 92 97 在（30，92） （40，97）范围内生成节理，节理倾角+30，迹长 50，贯通，间距 4jset -60 0 50 0 0 0 4 0 range x 30 40 y 92 97 在（30，92） （40，97）范围内生成节理，节理倾角-60，迹长 50，贯通，间距 4change mat 1 定义模型本构关系为线弹性change jmat 1 jcons 2 定义节理本构关系为完全弹性模型，set jmatdf 1 jcondf 2 定义计算开始后的新接触关系 prop m 1 d 0.0025 jkn 5000 jks 5000 jfr
9、ic 30 各种参数 damp auto 粘性阻尼damp .23 225 stiff ？set ovtol .1 接触叠加范围pl bl vel iw max 5 pl bl vel hold3dec 中关于结构面参数的确定前面看到有人说 3dec 结构面参数比较难以确定，我觉得还是可以确定，并可以和有限元的结果相媲美的，我曾经计算过大型工程的应力结果，其和有限元结果基本一致。一般而言，真实结构的变形由两方面组成，一是块体本身的变形，二是结构面之间的相互嵌入，也正是因为后者，很多人认为离散元结果人为性太大。其实不是这样的，如果是人为切割的分割面（真实结构中下不存在的） ，为了避免结构面的嵌入
10、带来的变形干扰正常的计算结果，一般只要保证 E/skn 约 20 倍左右即可，此时的结构面互相嵌入对于整个系统变形的影响几乎可以忽略不计。如果是真实存在的结构面，比如是天然状态下山体的节理裂隙，期间可能充填各种杂质，比如夹泥等等，此时需要使用较为真实的 kn 参数，一般取为 E/T，ks 可取 kn/(1-v） ，也可以取与 kn 一样，这个关系一般不大。上面的说明文字中，E 弹模，大家都知道s 接触面面积T 裂隙厚度设置出图精度： set nx=100 ny=100在力学分析角度上，离散元对三大定律的满足上与有限元方法不同。从平衡方程上看，离散元采用牛顿第二定律来控制，按围绕各刚性块体形心的
11、力平衡和力矩平衡来满足。从变形协调方程上看，各刚性块体间不再位移连续，而是允许大变形和断裂分开，可以模拟岩体不连续结构面的滑移与开裂。其位移联系为刚性块体间各接触点的相对速度关系。从材料本构关系上看，离散元法避开了复杂的本构关系推导，采用在刚性块体间设置不同种类弹簧和阻尼（法向刚度和阻尼、切向刚度和阻尼）来反映材料的应力位移关系。同时，在求解方法上，离散元法采用动态的松弛方法求解，从而避开了大规模刚度矩阵的形成，不须求解大型联立方程组，能够高效的求解非线性大变形问题。图形放大的快捷键是“M” ，缩小的快捷键是“U” 。单纯 jset 命令默认为 jset d 0 dd 0 org 0 0 03
12、dec 中建立的模型默认为 region 0delete 删掉的部分不能再回填， excavate 删掉的可以在用 fill 命令回填，remove 删掉的不可以回填，但是可以显示出来。显示塑性区：pl plasticity 显示位移：pl disp 显示速率：pl velocity 显示应力 pl syy 显示主应力矢量图：pl principal 显示主应力等值线：pl sintsmaxsmin 显示单元：pl zone3dec 命令流中每一行长度不宜超过 80 个字符，否则会出错，但注释长度可以任意如果因为块体形状或大小影响，单元划分不过去，有两种处理办法：一是在 gen edge xx
13、 后加关键词 alt ，如 gen edge 12.0 alt。或者把不能划分的小块体直接删除，如 del block 133395施加衬砌或喷层的命令流：config linerstruct liner radial_gen 0 0 -.9 0 0 .9 seg 2 8 prop 1 range .43struct prop 1 fric 60 coh .5e6 kn 1e9 ks 1e9struct prop 1 nu .15 thick .20 e 15e9 tens .3e6其中 range .43 是给定能够包住洞子临空面范围的半径。如果并不是对全断面施加liner，则还要给定起止平面上的各一个点。
