FLAC3D接触问题

marathon

在FLAC3D建模过程中设置了接触，为了计算初始应力场，除了用initial命令付值外， 我想用solve 或者solve elastic来求解初始应力场，请问接触面的参数应该如何取值才对初始应力场无影响？谢谢

tg0215

flac3d的interface也可以进行应力初始化。。。

marathon

tg0215，以不能举个例子呀？谢谢！

li_2001

举个例子 大家看看

marathon

;---------------------------------------------------
;
;  Vertically and horizontally loaded pile in clay
;
;---------------------------------------------------
gen zone radcyl p0  (0,0,0)  p1  (8,0,0)  p2  (0,0,-5)  p3  (0,8,0)  &
                p4  (8,0,-5) p5  (0,8,-5) p6  (8,8,0)   p7  (8,8,-5) &
                p8  (.3,0,0)  p9  (0,.3,0)  p10 (.3,0,-5)  p11 (0,.3,-5) &
                size 3 10 6 15  ratio 1 1 1 1.15
gen zone radcyl p0  (0,0,-5)  p1  (8,0,-5)  p2  (0,0,-8)  p3  (0,8,-5)  &
                p4  (8,0,-8)  p5  (0,8,-8)  p6  (8,8,-5)  p7  (8,8,-8)  &
                p8  (.3,0,-5) p9  (0,.3,-5) p10 (.3,0,-8) p11 (0,.3,-8) &
                size 3 6 6 15  ratio 1 1 1 1.15 fill
gen zone reflect dd 270 dip 90
group clay
;
interface 1 face range cylinder end1 (0,0,0)  end2 (0,0,-5.1) radius .31 &
                      cylinder end1 (0,0,0)  end2 (0,0,-5.1) radius .29 not
interface 2 face range cylinder end1 (0,0,-4.9)  end2 (0,0,-5.1) radius .31
;
gen zone cyl p0 (0,0,6) p1 (.3,0,6) p2 (0,0,1)  p3 (0,.3,6)  &
             p4 (.3,0,1) p5 (0,.3,1) &
             size 3 10 6
gen zone cyl p0 (0,0,6.1) p1 (.3,0,6.1) p2 (0,0,6)  p3 (0,.3,6.1)  &
             p4 (.3,0,6) p5 (0,.3,6) &
             size 3 1 6
gen zone reflect dd 270 dip 90 range z 1 6.1
group pile range z 1 6.1
ini z add -6.0 range group pile
save pile_geom.sav

;pause

;
model mohr                                          range group clay
prop bulk 8.333e7  shear 3.846e7  coh 30000 fric 0  range group clay
model elas                        range group pile
prop bulk 8.333e7  shear 3.846e7  range group pile


interface 1 prop kn 0  ks 0  fric 0 coh 30000
interface 2 prop kn 0  ks 0  fric 0 coh 30000
;
ini dens 1230 range group clay
ini dens 1230 range group pile
model null range z -0.1 0.15
;
fix z range z -8.1 -7.9
fix x range x -8.1 -7.9
fix x range x 7.9 8.1
fix y range y -.1 .1
fix y range y 7.9 8.1
set grav 0 0 -10
solve

比如这个例子中，如何设置接触面参数来计算初始应力场，谢谢

tianshen

我觉得 接触面的 参数 会影响初始地应力的快慢，个人感觉最好是设置大一点的接触面参数，比如摩擦角或者粘聚力等 使其不影响土体间的应力传递，

feigle

学习下  ^_^

月三更

:)谢谢楼主分享

dead5845065

好有用呀！谢谢啦！

xzh2012

最近也在学习FLAC，谢谢各位楼主！[

wx86_68xw

marathon 发表于 2009-10-20 16:51
;---------------------------------------------------
;
;  Vertically and horizontally loaded pile in ...

请问interface 2 face range cylinder end1 (0,0,-4.9)  end2 (0,0,-5.1) radius .31
的目的是什么？
以及gen zone cyl p0 (0,0,6.1) p1 (.3,0,6.1) p2 (0,0,6)  p3 (0,.3,6.1)  &
             p4 (.3,0,6) p5 (0,.3,6) & size 3 1 6
的目的是什么？
一个初学者，请教一下高手，呵呵