接触材料表(Contact Material Table)的注意事项

1 引言

在小应变模式下，通过设置节理的本构关系(block contact jmodel assign)和属性(block contact property)【连续屈服节理模型(continuously yielding joint model)】即可进行计算，但在大应变模式下(model large-strain on), 必须为新产生的接触定义接触本构模型和属性【隧道开挖的BBM模拟---材料属性(空间变异性)】，否则不能进行计算，这个设置是通过接触材料表(Contact Material Table)来实现的(目前Itasca笔记共有contact material-table的条目52个)。

2 接触材料表

在大应变模式下，必须首先使用关键字defalut设置默认的节理模型和属性，如果不指定的话，默认的本构模型是mohr-coulomb。

block contact material-table default prop ...  stiff-norm 1e9 stiff-shear 1e9 fric 30

为了对新形成的接触设置不同的属性，需要使用block contact material-table命令，可以使用的关键字有：add, apply, default, list, modify, remove。add关键字用来设定一个条目(entry)，该条目由节理模型、节理属性和范围(Range)组成。

block contact material-table add jmodel mohr prop ...  stiff-norm 1e9 stiff-shear 1e9 fric 30 ...  range group-intersection 'slider' 'base1'

当使用add关键字时，range必不可少，否则会出现错误信息:"你需要指定一个过滤器来添加材料表的条目，或者使用block contact material-table default命令来修改默认槽位!(

You need to specify a filter to add an entry in the material table, or use the block contact material-table default command to modify default slot!)"。

下图所示的例题与【单个块体在边坡面上的稳定性模拟---滑动和倾倒破坏】讨论的问题类似。在一个由三种不同材料类型组成的边坡(block group)上，一个单个块体从边坡顶滑下。默认的摩擦角是30°，但是滑动块和边坡中间块base2之间接触的摩擦角是40°，由于边坡角是35°，所以块体在base1上滑动，但当它与sbase2接触后就会停止。

3 关于group-intersection

在上面改变接触的过程中，使用了range group-intersection关键字来改变slider与base2之间的摩擦属性，group-intersection只接受两个对象的相互作用，如果需要改变slider与base3之间的属性，必须另写一个contact material-table。此外，range group-intersection 'A' 'B' 只满足了语法上是否合理，并没有检查两个对象'A'和'B'是否真实存在或者合理地相互作用，如果这个设置不reasonable，那么相当于这条命令无效，按default的设置进行计算。