- 积分
- 0
- 注册时间
- 2008-12-26
- 仿真币
-
- 最后登录
- 1970-1-1
|
本帖最后由 marathon 于 2009-10-22 16:39 编辑
我采用了两种方法进行了接触面的建立,方法一接触面上有法向应力,(本人认为这种方法建立的接触面是正确的);方法二接触面上无法向应力,(本人认为这种方法建立的接触面是有问题的),现在请高手指点一下为什么这两种方法有差别。特别是希望高手指导一下陈育民博士的书中的“切割模型法”建立接触面,做一个比较详细的分析。重点说明一下“切割模型法”用到的关键命令:
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 |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|