请高手指点一下接触面

marathon

我采用了两种方法进行了接触面的建立，方法一接触面上有法向应力，（本人认为这种方法建立的接触面是正确的）；方法二接触面上无法向应力，（本人认为这种方法建立的接触面是有问题的），现在请高手指点一下为什么这两种方法有差别。特别是希望高手指导一下陈育民博士的书中的“切割模型法”建立接触面，做一个比较详细的分析。重点说明一下“切割模型法”用到的关键命令：
gen separate group 1
int 1 wrap group 1 group 2

方法一：（彭文斌《FLAC3D实用教程》P278例题14-2
new
gen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,2

gen zone brick size 8 8 4 p0 0,0,3 p1 4,0,3 p2 0,4,3 p3 0,0,5

inter 1 face range z 2.9,3.1
inter 1 prop kn 300e9 ks 300e9 tens 1e10 SBRATIO = 1

ini z add -1.0 range z 2.9,5.1
model elas
prop bulk 8e9 shear 5e9
fix z range z -.1 .1
fix x range x -.1 .1
fix x range x 3.9 4.1
fix y range y -.1 .1
fix y range y 3.9 4.1
apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2
hist unbal
;pause
solve
方法二：（参照“切割模型法”）
new
gen zone brick size 4 4 4 p0 0,0,0 p1 4,0,0 p2 0,4,0 p3 0,0,2 &
group g1
gen zone brick size 8 8 4 p0 0,0,2 p1 4,0,2 p2 0,4,2 p3 0,0,4 &
group g2
gen separate g2
inter 1 face range plane ori 0,0,2 normal 0,0,1 above z=2

inter 1 prop kn 300e9 ks 300e9 tens 1e10 SBRATIO = 1

;ini z add -1.0 range z 2.9,5.1
model elas
prop bulk 8e9 shear 5e9
fix z range z -.1 .1
fix x range x -.1 .1
fix x range x 3.9 4.1
fix y range y -.1 .1
fix y range y 3.9 4.1
apply szz -1e6 range z 3.9 4.1 x 0,2 y 0,2
hist unbal
;pause
solve

marathon

为了方便各位高手指导，现上传建模命令流文件，免得复制上面的命令流时出错乱码。谢谢

franco_chao

我用的是第一种方法，既移来移去。很方便。也很好用，第二种没用过。

cloud328

陈博士提出的几种方法都没有问题的。

marathon

4# cloud328

我也知道陈博士提出的几种方法没有问题，但是上面的例子，你可不可以帮我解释一下原因呀，谢谢

abblcn

移来移去法，不过切割法更加的方便撒

huhaixin2009

切割法方便！