初始地应力场生成程序求解 [复制链接]

;网格建立
gen zone brick p0 3 0 -9 p1 16.5 0 -9 p2 3 1 -9 p3 3 0 0 size 27 2 18 rat 1 1 1
gen zone brick p0 16.5 0 -9 p1 25 0 -9 p2 16.5 1 -9 p3 16.5 0 0 size 10 2 18 rat 1 1 1
gen zone brick p0 0 0 -14 p1 3 0 -14 p2 0 1 -14 p3 0 0 -9 size 6 2 10 rat 1 1 1
gen zone brick p0 3 0 -14 p1 16.5 0 -14 p2 3 1 -14 p3 3 0 -9 size 27 2 10 rat 1 1 1
gen zone brick p0 16.5 0 -14 p1 25 0 -14 p2 16.5 1 -14 p3 16.5 0 -9 size 10 2 10 rat 1 1 1
gen zone reflect ori 0 0 0 norm 1 0 0
;生成初始应力场
model mohr
prop bulk 500e6 sh 330e6 coh 2e3 fri 25
fix x y z ran z -14.1 -13.9
fix x ran x 24.9 25.1
fix x ran x -25.1 -24.9
fix y ran y 0.9 1.1
fix y ran y -1.1 -0.9
set grav 0 0 -10
ini dens 2000
set large
hist unbal
set mech ratio 1e-4
solve elas
最后出现Illegal geometry in zone 4223，是因为模型已经破坏还是网格建立有问题？

变形太大，网格就有问题了。可以使用Plot查看，肯定变形异常啦！

变形异常？？

你再重新显示一下block group，就会看到网格已经变形了。