3DEC出现的诡异问题: 点不能正确地置于平面内

1. 引言

本文首先简要回顾了Itasca软件最新版本的新增功能，尤其强调了3DEC的新命令；其次记录了在3DEC中遇到的一个非常诡异的问题：点不能正确地置于平面内。

2. 9.00.165新增命令

(1) DFN断裂尺寸对数正态(log-normal)分布;【正态分布和对数正态分布的应用(Gauss and Log-normal Distribution)】

3. 为什么点不能置于平面内

这个问题的起源来自于使用block zone generate-new命令局部加密网格，本来想再现生成网格时出现的错误信息：Two nodes are too close to each other (same position or almost)，试验结果显示没有出现这个错误信息，但引发了一个非常奇怪的问题：点不能正确地置于平面(节理面)内。

在下面的图中，点4(pt=4)的坐标值取自平面内，这个值从Information窗口记录下来，应该是绝对正确的，但就是不能生成网格。

猜想这个问题可能与tolerance或gap有关，具体啥原因暂时还没有想出来。

4. 本试验命令

(1) block create brick 0 10

(2) block cut joint-set origin (0, 0, 0) dip 60 dip-direction 90

(3) [global p1 = vector (x, y, z)]

(4) data scalar add [p1] group "pt=1" 

(5) block zone size isolated-point [p1] edgelength 1

(6) block zone generate-new ...

    gradation-surface 0.2 ...

    gradation-volume 0.2 ...

    maximum-edge 10 ...

    minimum-edge 1

与generate-new相关的内容：

• 3DEC单元类型的选择(Zone Types)

• 自创与导入块体模型比较（create and import)

• 压缩试验模拟考虑的几个问题(本构模型和NMD算法)

• GeotechGPT训练产生的---带承压板的Voronoi块体单轴抗压强度试验