各位大神，小弟最近在模拟土钉墙支护，但是计算的结果显示坑底的土位移是向上的，不知道是什么问题，特请各位大神指导一下。命令
新命令.txt (6 K) 下载次数:30 [attachment=251839]和模型的cad图
模型坐标.rar (13 K) 下载次数:23 见附件。小弟恳请各位大神伸出援手
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
