【求助】为何初始应力场不能平衡？
我的算例在生成处理应力场的时候总是在最大不平衡力10-4的地方开始不能继续下去，无法达到10-5。不知道为什么我的模型首先将粘聚力和抗拉强度设的很大生成土体的初始应力，但是不能平衡，我试着变换四个参数的值，但是总是的不到平衡,
; gen geometry
;--------------------------------------------------------------------------
;gen radbrick
;zuo
gen zone brick p0 16 0 0 p1 80 0 0 p2 16 0 -18 p3 16 9 0 p4 80 0 -80 p5 16 9 -18 &
p6 80 45 0 p7 80 45 -80 size 20 18 9 rat 1.1 1 1
gen zone brick p0 16 9 0 p1 80 45 0 p2 16 9 -18 p3 16 10 0 p4 80 45 -80 p5 16 10 -18 &
p6 80 50 0 p7 80 50 -80 size 20 18 2 rat 1.1 1 1
;you
gen zone brick p0 0 10 0 p1 15 10 0 p2 0 10 -18 p3 0 50 0 p4 15 10 -18 p5 0 50 -80 &
p6 75 50 0 p7 75 50 -80 size 15 18 20 rat 1 1 1.1
gen zone brick p0 15 10 0 p1 16 10 0 p2 15 10 -18 p3 75 50 0 p4 16 10 -18 p5 75 50 -80 &
p6 80 50 0 p7 80 50 -80 size 2 18 20 rat 1 1 1.1
;xia
gen zone brick p0 0 0 -18 p1 15 0 -18 p2 0 0 -80 p3 0 9 -18 p4 75 0 -80 p5 0 45 -80 &
p6 15 9 -18 p7 75 45 -80 size 15 20 9 rat 1 1.1 1
gen zone brick p0 15 0 -18 p1 16 0 -18 p2 75 0 -80 p3 15 9 -18 p4 80 0 -80 p5 75 45 -80 &
p6 16 9 -18 p7 80 45 -80 size 2 20 9 rat 1 1.1 1
gen zone brick p0 0 9 -18 p1 15 9 -18 p2 0 45 -80 p3 0 10 -18 p4 75 45 -80 p5 0 50 -80 &
p6 15 10 -18 p7 75 50 -80 size 15 20 2 rat 1 1.1 1
gen zone brick p0 15 9 -18 p1 16 9 -18 p2 75 45 -80 p3 15 10 -18 p4 80 45 -80 p5 75 50 -80 &
p6 16 10 -18 p7 80 50 -80 size 2 20 2 rat 1 1.1 1
;top
gen zone brick p0 0 0 30 p1 15 0 30 p2 0 0 12 p3 0 9 30size 15 18 9
gen zone brick p0 15 0 30 p1 16 0 30 p2 15 0 12 p3 15 9 30 size 2 18 9
gen zone brick p0 0 9 30 p1 15 9 30 p2 0 9 12 p3 0 10 30 size 15 18 2
gen zone brick p0 15 9 30 p1 16 9 30 p2 15 9 12 p3 15 10 30 size 2 18 2
group soil
;def interface
interface 1 face range x 16 16 y 0 10 z -18 0
interface 2 face range x 0 16 y 10 10 z -18 0
ini z add -30 range x 0 16 y 0 10 z 12 30
;wall and cuttion
group wal1 range x 15.5 16 y 0 10 z -18 0
group wal2 range x 0 15.5 y 9.5 10 z -18 0
group wall rangegroup wal1 any group wal2 any
;
mod mohr range group soil
pro bulk=6e6shear=2e6f=20 coh=1e9ten=1e9range group soil
mod mohr range group wall
pro bulk=6e6shear=2e6f=20 coh=1e9ten=1e9range group wall
;
interface 1 prokn 1e8ks 1e8 fric 15 coh 5e3
interface 2 prokn 1e8ks 1e8 fric 15 coh 5e3
;
ini dens 1600 range group soil
ini dens 1600 range group wall
ini szz 0.0 grad 0 0 1.6e4 range z -80 0
;--------------boun----------------
fix x range x -0.1 0.1
fix x range x 79.9 80.1
fix y range y -0.1 0.1
fix y range y 49.9 50.1
fix x y z range z -80.1-79.9
;
set g 0 0 -10
set small
;
hist unbal
solve
save inisoil.sav
;---------ini stress-wall
mod mohr range group soil
pro bulk=6e6shear=2e6f=20coh=1e4 range group soil
mod elas range group wall
pro bulk=13.9e9shear=10.4e9 range group wall
ini dens 2400 range group wall
hist unbal
solve
save iniwall.sav
回复： 【求助】为何初始应力场不能平衡？
这是几何模型回复： 【求助】为何初始应力场不能平衡？
我的建议：1，看看boundary是否设置正确；
2，去掉interface看看能够收敛，如果收敛，说明interface的参数有问题；
3，既然你把c，fi，tension设置10e9，应该不会产生塑性流，但还是建议检查一下是否产生塑性流；
4，网格（特别是地下墙）的网格的b/h》》5，是否可以优化一下。
回复：starsmoon 【求助】为何初始应力场不能平衡？
谢谢starsmoon ：1边界条件不应该又大的问题，
2.interface的参数设置我还是不太懂，具体应该怎么样设？遵守怎样的原则？
3.塑性流动如何查询？我没有做过。
4.网格b/h>>5?????????,在哪里？
另外，我做了一下尝试，在我另外一个帖子里提到了，我吧命令流中的
ini z add -30 range x 0 16 y 0 10 z 12 30
改成直接生成网格，就可以在在第一步中收敛，
这个我不明白为什么？
但在第二步中，就是我吧混凝土墙加进去后，在到了2.5e-5后就不能在边小了，这又意味着什么呢？
请指教
页:
[1]