建立interface的几种方法的讨论
首先给出一个简单的例子 ，大家讨论讨论哪种方法最好：new
gen zone brick p0 0 0 0 p1 111.4 0 0 p2 0 10 0 p3 0 0 45 size 11 1 5 group 1
gen zone brick p0 0 0 45 p1 111.4 0 45 p2 0 10 45 p3 0 0 90 p4 111.4 10 45 &
p5 0 10 90 p6 47 0 90p7 47 10 90 size 11 1 9 group 3
gen zone wedge p0 66.4 0 90.1 p1 47 0 90.1 p2 66.4 10 90.1 p3 111.4 0 45.1 &
p4 47 10 90.1 p5 111.4 10 45.1 size 5 1 9 group 4
;interface 1 face rang plane norm 45,0,64.4 origin 47 5 90
interface 1 face range plane norm45,0,64.4 origin 47 5 90 dist 0.1
ini z add -0.1 range group 4
gen zone brick p0 111.4 0 0 p1 128 0 0 p2 111.4 10 0 p3 111.4 0 45 &
size 2 1 5 group 2
model mohr
pro bulk=1e8 she=3e7 dens 2e3 coh 12e3 &
friction 20 dilation 0. tension 1.0e10range group 1
pro bulk=1e8 she=3e7 dens 2e3 coh 12e3 &
friction 20 dilation 0. tension 1.0e10range group 2
int 1 prop ks 2e12 kn 2e12 fric 36
pro bulk=1e8 she=3e7 dens 2e3 coh 12e3 &
friction 20 dilation 0. tension 1.0e10range group 3
pro bulk=1e8 she=3e7 dens 2e3 coh 12e3 &
friction 20 dilation 0. tension 1.0e10range group 4
fix z range z -0.1 0.1
fix x range x -0.1 0.1
fix x range x 127.9 128.1
fix y
set gravity 0 0 -10
step 1000
ini state 0
ini xdis 0 ydis 0 zdis 0
;ini xvel 0 yvel 0 zvel 0
;set large
step 1000
;solve fosfile cable.sav associated
Re:建立interface的几种方法的讨论
修改上面蓝色字体部分第一种建立interface的方法（书上的移来移去式）：
gen zone brick p0 0 0 0 p1 111.4 0 0 p2 0 10 0 p3 0 0 45 size 11 1 5 group 1
gen zone brick p0 0 0 45 p1 111.4 0 45 p2 0 10 45 p3 0 0 90 p4 111.4 10 45 &
p5 0 10 90 p6 47 0 90p7 47 10 90 size 11 1 9 group 3
gen zone wedge p0 66.4 0 90.1 p1 47 0 90.1 p2 66.4 10 90.1 p3 111.4 0 45.1 &
p4 47 10 90.1 p5 111.4 10 45.1 size 5 1 9 group 4
;interface 1 face rang plane norm 45,0,64.4 origin 47 5 90
interface 1 face range plane norm45,0,64.4 origin 47 5 90 dist=0.1
ini z add -0.1 range group 4
gen zone brick p0 111.4 0 0 p1 128 0 0 p2 111.4 10 0 p3 111.4 0 45 &
size 2 1 5 group 2
通过计算的到xdis 的云图
Re:建立interface的几种方法的讨论
第二种方法(lakewater的方法)：gen zone brick p0 0 0 0 p1 111.4 0 0 p2 0 10 0 p3 0 0 45 size 11 1 5 group 1
gen zone brick p0 0 0 45 p1 111.4 0 45 p2 0 10 45 p3 0 0 90 p4 111.4 10 45 &
p5 0 10 90 p6 47 0 90p7 47 10 90 size 11 1 9 group 3
;interface 1 face rang plane norm 45,0,64.4 origin 47 5 90
gen zone wedge p0 66.4 0 90 p1 47 0 90 p2 66.4 10 90 p3 111.4 0 45 &
p4 47 10 90 p5 111.4 10 45 size 5 1 9 group 4
gen zone brick p0 111.4 0 0 p1 128 0 0 p2 111.4 10 0 p3 111.4 0 45 &
size 2 1 5 group 2
save model.sav
del ran group 4 any group 2 any
interface 1 face rang plane norm 45,0,64.4 origin 47 5 90
save joint1.sav
rest model.sav
delete ran group 1 any group 3 any
expgrid 1.flac3d
rest joint1.sav
impgrid 1.flac3d
save model1.sav
得到的xdis的云图
Re:建立interface的几种方法的讨论
第三种方法：先建立interface 下面的
建立interface
再建立interface 上面的
gen zone brick p0 0 0 0 p1 111.4 0 0 p2 0 10 0 p3 0 0 45 size 11 1 5 group 1
gen zone brick p0 0 0 45 p1 111.4 0 45 p2 0 10 45 p3 0 0 90 p4 111.4 10 45 &
p5 0 10 90 p6 47 0 90p7 47 10 90 size 11 1 9 group 3
interface 1 face rang plane norm 45,0,64.4 origin 47 5 90
gen zone wedge p0 66.4 0 90 p1 47 0 90 p2 66.4 10 90 p3 111.4 0 45 &
p4 47 10 90 p5 111.4 10 45 size 5 1 9 group 4
gen zone brick p0 111.4 0 0 p1 128 0 0 p2 111.4 10 0 p3 111.4 0 45 &
size 2 1 5 group 2
得到的xdis 的云图
Re:建立interface的几种方法的讨论
首先，你的比较研究的思路是值得学习的！但是，在你所说的三种比较程序里，接触面的设置都是错误的。
第一种得到的结果是：
Re:建立interface的几种方法的讨论
第二种和第三种得到的接触面设置也不正确，只有一个面上：Re:建立interface的几种方法的讨论
正确的设置你可以再研究一下。另外，需要提醒的是，这种直接通道坡脚的接触面可能会有些问题。因为，在坡趾的这一点接触情况没有明确的界定，可能是我对interface corner这个细节还不是很了解。
很抱歉，我也没能给出正确解答，大家都需要再努力！
Re:建立interface的几种方法的讨论
按道理，各种接触面的设置方法只要正确，得出的结果应该是一样的。如果结果不同，那么肯定这种设置方法本身，或者命令编写有问题。接触面的设置最终的结果就是产生接触单元，简单的说就是在同一个结点的位置有两个点，这样才可能相对滑动、拉裂等等。由于FLAC3D本身的接触面设置比较麻烦，所以我才建议通过expgrid的方法建立，原理上是一样的。
这个问题还有继续讨论的必要和意义。
Re:建立interface的几种方法的讨论
我的和你的结果不一样哦第三种方法
Re:建立interface的几种方法的讨论
第二种方法Re:建立interface的几种方法的讨论
plo inter yel sk这个命令？
Re:建立interface的几种方法的讨论
还有第一种方法 贴错了啊要加 dist=0.1 就可以了
Re:建立interface的几种方法的讨论
pl int yel sket继续看看为什么第一种方法要加 dist 0.1？？ 。。
还有就是加入节理以后，这个例子就不能计算
fos 了
？？？
大家试试
Re:建立interface的几种方法的讨论
你的第2种，把group2和group4删掉，这样造成的新的模型group2和group4有重复节点，这样也是不正确的。Re:建立interface的几种方法的讨论
lakewater wrote:你的第2种，把group2和group4删掉，这样造成的新的模型group2和group4有重复节点，这样也是不正确的。
呵呵。。第二种好像是你的方法啊
要怎么改 ？？
期待中。。。
好像就我们两个在这里讨论哦 。。。
Re:建立interface的几种方法的讨论
toss wrote:呵呵。。第二种好像是你的方法啊
要怎么改 ？？
期待中。。。
好像就我们两个在这里讨论哦 。。。
你可以看看以前我发的那个expgrid和impgrid命令的讨论的帖子。
我的方法是在一个实体的外表面全部都设置了接触面，就是可以将两个体完全分开。而你的这个接触问题，在坡脚的位置很难处理，或者说你处理得不当。应该只用把斜坡上的部分删掉，然后导入就可以了。
Re:建立interface的几种方法的讨论
是这样的把interface右边的 group 2 group 4 删掉
然后加interface
这样是不是也把 他们给隔开了？
然后通过 导入节点group 1 group 3
得到整个实体 、、、
Re:建立interface的几种方法的讨论
toss wrote:是这样的
把interface右边的 group 2 group 4 删掉
然后加interface
这样是不是也把 他们给隔开了？
然后通过 导入节点group 1 group 3
得到整个实体 、、、
这样你人为的把group1和group2之间的交界面上多出了节点，下图是两个group的节点图，是分开显示的，可以看到两个共同面上界面标号不同，这样是不对的。计算结果也同样是错误的。
Re:建立interface的几种方法的讨论
好像这样建模的结果都是会有这样的啊可以试试
n
gen zone brick p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 group 1
gen zone brick p0 10 0 0 p1 15 0 0 p2 10 10 0 p3 10 0 10 group 2
这样得到他们的交界面上结点的编号也不一样的
不知道 我有没有弄懂你的意思？？
Re:建立interface的几种方法的讨论
toss wrote:好像这样建模的结果都是会有这样的啊
可以试试
n
gen zone brick p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 group 1
gen zone brick p0 10 0 0 p1 15 0 0 p2 10 10 0 p3 10 0 10 group 2
这样得到他们的交界面上结点的编号也不一样的
不知道 我有没有弄懂你的意思？？
这样建，在x=10的面上节点就是1个了。节点编号是一致的。