请教关于基坑土钉墙支护分析的问题

ziguayida

各位大神，小弟最近在模拟土钉墙支护，但是计算的结果显示坑底的土位移是向上的，不知道是什么问题，特请各位大神指导一下。命令和模型的cad图见附件。小弟恳请各位大神伸出援手
gen zo b p0=(-3.5 0 0) p1=(0 0 0) p2=(-3.5 0.108 0) p3=(-3.5 3.14 7.85)  p4=(0 0.108 0) p5=(-3.5 3.248 7.85) p6=(0 3.14 7.85) p7=(0 3.248 7.85) size 7 1 16 group penshehunningtu
sel shell id=1 elemtype=dkt crossdiag range group penshehunningtu
sel group shotceng range z=0 7.85
sel cable id=1 begin=(-0.89,0.1,0.25) end=(-0.89,6.909,-1.575) nseg=10
sel cable id=2 begin=(-2.610,1.140,1.55) end=(-2.610,7.429,-0.275) nseg=10
sel cable id=3 begin=(-0.89,1.14,2.85) end=(-0.89,8.922,0.765) nseg=10
sel cable id=4 begin=(-2.61,1.7,4.25) end=(-2.61,9.482,2.165) nseg=10
sel cable id=5 begin=(-0.89,2.26,5.65) end=(-0.89,10.42,3.565) nseg=10
sel cable id=6 begin=(-2.610,2.7,6.75) end=(-2.61,10.482,4.665) nseg=10
sel group nail_1 range id=1
sel group nail_2 range id=2
sel group nail_3 range id=3
sel group nail_4 range id=4
sel group nail_5 range id=5
sel group nail_6 range id=6
gen zo b p0=(-3.5 -15 -15) p1=(0 -15 -15) p2=(-3.5 0 -15) p3=(-3.5 -15 0) size 7 30 30 group xiawoceng1
gen zo b p0=(-3.5 0 -15) p1=(0 0 -15) p2=(-3.5 0.108 -15) p3=(-3.5 0 0)  size 7 1 30 group xiawoceng2
gen zo b p0=(-3.5 0.108 -15) p1=(0 0.108 -15) p2=(-3.5 20 -15) p3=(-3.5 0.108 0) size 7 40 30 group xiawoceng3
gen zo b p0=(-3.5 -15 0) p1=(0 -15 0) p2=(-3.5 0 0) p3=(-3.5 -15 7.85)  p4=(0 0 0) p5=(-3.5 3.14 7.85) p6=(0 -15 7.85) p7=(0 3.14 7.85) size 7 30 16 group watu
gen zo b p0=(-3.5 0.108 0) p1=(0 0.108 0) p2=(-3.5 20 0) p3=(-3.5 3.248 7.85)  p4=(0 20 0) p5=(-3.5 20 7.85) p6=(0 3.248 7.85) p7=(0 20 7.85) size 7 40 16 group bianpo
save model.sav
new
restore model.sav
set echo on
set plot jpg
set plot jpg si 1600 1000
plot set ba wh
plot set cap right
plot set cap si 25
plot set rot 30 0 30

set plot jpg
set plot jpg si 2880 1800
plot set mag=1.1
plot bl gr axes
set outp soil_group.jpg
pl ha
pl ske axes
pl add sel geom node off sh 0
set outp model_geom.jpg
pl ha
plot set cent 0 20 0 mag=3.2
plot sel group  alias '结构单元分组' axes blue
set outp sel_group.jpg
pl ha
sel group shot1 range group shotceng z=(6.35 7.85)
sel group shot2 range group shotceng z=(4.85 6.35)
sel group shot3 range group shotceng z=(3.35 4.85)
sel group shot4 range group shotceng z=(1.85 3.35)
sel group shot5 range group shotceng z=(0.85 1.85)
sel group shot6 range group shotceng z=(0 0.85)
;group wazou range group7
group watu1 range group watu z=(6.35 7.85)
group watu2 range group watu z=(4.85 6.35)
group watu3 range group watu z=(3.35 4.85)
group watu4 range group watu z=(1.85 3.35)
group watu5 range group watu z=(0.85 1.85)
group watu6 range group watu z=(0 0.85)
range name s1 z=(-15 1.55)
range name s2 z=(1.55 3.35)
range name s3 z=(3.35 6.35)
range name s4 z=(6.35 7.85)
;=================initial state===========
model elas
prop y=40e6 p=0.35 c=8e3 f=10.0 t=1e3 density=1.78e3 range s1
prop y=60e6 p=0.35 c=12e3 f=14.0 t=1e3 density=1.83e3 range s2
prop y=60e6 p=0.35 c=21e3 f=21.0 t=1e3 density=1.91e3 range s3
prop y=60e6 p=0.30 c=5e3 f=32.0 t=1e3 density=1.99e3 range s4
;=========================================
fix x range x=(-3.6 -3.4) any x=(-0.1 0.1)
fix y range y=(-15.1 -14.9) any y=(19.9 20.1)
fix z range z=(-15.1 -14.9)
set gravity=(0 0 -10)
ini szz 0.0 gradient=(0 0 18e3)
ini sxx 0.0 gradient=(0 0 9e3)
ini syy 0.0 gradient=(0 0 9e3)
;===================施加地面15kpa超载====
apply szz -15e3 range y=(3.14,20) z=(7.75,7.95)
;=========================================
cyc 1
sel node fix lsys
;=============重先建立link================
sel node ini ypos add 0.0
;==============结构边界===================
sel node fix x yr zr range sel shell
;==============结构单元赋值==============
sel shell prop iso=(1e-10,0.2) density=0.0
sel cable prop emod=1e-10 density=0.0 gr_coh=1e20 gr_fri=14 gr_k=14e9 gr_per=0.314 slide=on slide_tol=1e-6 ycomp=1e20 ytens=1e20
;========================================
set mech damp combined
sel set damp combined
hist unbal
pl hist 1
step 3000
save initial.sav
;====================开挖================
new
res initial.sav
ini xdisp=0 ydisp=0 zdisp=0
ini xvel=0 yvel=0 zvel=0
sel node ini xdisp=0 ydisp=0 zdisp=0
sel node ini xvel=0 yvel=0 zvel=0
ini state=0
;====================土钉轴力============
;===================stage1===============
model null range group watu1
step 2000
sel cable prop emod=210e9 density=7.8e3 gr_coh=3.14e4 gr_fri=10 gr_k=24e9 gr_per=0.314 slide=on slide_tol=1e-6 ycomp=16.5e6 ytens=16.5e6
sel shell prop iso=(10e9 0.2) density=2.5e3 thick=0.1 range group shot1
save stage1.sav
;===================stage2===============
model null range group watu2
step 2000
sel cable prop emod=210e9 density=7.8e3 gr_coh=3.14e4 gr_fri=10 gr_k=24e9 gr_per=0.314 slide=on slide_tol=1e-6 ycomp=16.5e6 ytens=16.5e6
sel shell prop iso=(10e9 0.2) density=2.5e3 thick=0.1 range group shot2
save stage2.sav
;===================stage3===============
model null range group watu3
step 2000
sel cable prop emod=210e9 density=7.8e3 gr_coh=3.14e4 gr_fri=10 gr_k=24e9 gr_per=0.314 slide=on slide_tol=1e-6 ycomp=16.5e6 ytens=16.5e6
sel shell prop iso=(10e9 0.2) density=2.5e3 thick=0.1 range group shot3
save stage3.sav
;===================stage4===============
model null range group watu4
step 2000
sel cable prop emod=210e9 density=7.8e3 gr_coh=3.14e4 gr_fri=10 gr_k=24e9 gr_per=0.314 slide=on slide_tol=1e-6 ycomp=16.5e6 ytens=16.5e6
sel shell prop iso=(10e9 0.2) density=2.5e3 thick=0.1 range group shot4
save stage4.sav
;===================stage5===============
model null range group watu5
step 2000
sel cable prop emod=210e9 density=7.8e3 gr_coh=3.14e4 gr_fri=10 gr_k=24e9 gr_per=0.314 slide=on slide_tol=1e-6 ycomp=16.5e6 ytens=16.5e6
sel shell prop iso=(10e9 0.2) density=2.5e3 thick=0.1 range group shot5
save stage5.sav
;===================stage6===============
model null range group watu6
step 2000
sel cable prop emod=210e9 density=7.8e3 gr_coh=3.14e4 gr_fri=10 gr_k=24e9 gr_per=0.314 slide=on slide_tol=1e-6 ycomp=16.5e6 ytens=16.5e6
sel shell prop iso=(10e9 0.2) density=2.5e3 thick=0.1 range group shot6
save stage6.sav

ziguayida

c沉降图是这样的

ziguayida

y向位移是这样的

ziguayida

切应变增量是这样的

ziguayida

算出来的结果看上去觉得土体是分离开的，请大神指导一下

ziguayida

跪求各位大神的指导

ziguayida

在线等指导

ziguayida

请懂的大神指导一二

ziguayida

求指导啊！！

ziguayida

没人施加援手啊，郁闷

gongzhenyu

基坑开挖之后，上面的压力不存在了，基坑底部当然会隆起了
这个很正常，据我所知，在FLAC里面计算，这个是有这么个问题的，不用担心
这个也是FLAC程序的一个弱点。

ziguayida

gongzhenyu 发表于 2012-12-26 11:37
基坑开挖之后，上面的压力不存在了，基坑底部当然会隆起了
这个很正常，据我所知，在FLAC里面计算，这个是 ...

但是挖完基坑后，ssi和yzdisp真的很不正常

ziguayida

ziguayida 发表于 2012-12-24 09:26
没人施加援手啊，郁闷

大神再给点指点呗

ziguayida

gongzhenyu 发表于 2012-12-26 11:37
基坑开挖之后，上面的压力不存在了，基坑底部当然会隆起了
这个很正常，据我所知，在FLAC里面计算，这个是 ...

大神再给点指点呗

gongzhenyu

ziguayida 发表于 2012-12-26 11:47
大神再给点指点呗

很多论文里面的计算，都有这个问题，这个问题其实很正常，也很好理解
未开挖之前，是个平衡状态，开挖之后，上面的荷载释放了，当然会隆起了
这个是FLAC程序的一个计算弱点，不能克服，只要你会解释就行
去下载点相关的论文看看吧。

ziguayida

gongzhenyu 发表于 2012-12-26 11:49
很多论文里面的计算，都有这个问题，这个问题其实很正常，也很好理解
未开挖之前，是个平衡状态，开挖之 ...

大神，我明白你说的意思，这个问题也好理解，我说的是，我的计算结果不正常，12楼是ssi，ydisp，zdisp。我觉得程序有些东西没设置好，您能帮忙看看吗？

ziguayida

谁能告诉一下model mohr
prop y=40e6 p=0.35 c=5e3 f=10.5 t=1e3 density=1.82e3 range s1中y和p的是什么意思，为什么mohr模型中没有体积模量和切变模量？