锚索抗滑桩支档边坡 [复制链接]

我遇到的问题是


我的命令流是
new
gen zone brick p0 0 0 0 p1 18 0 0 p2 0 5 0 p3 0 0 8 size 17 10 8  ratio 1.03 1 1 group 1  
gen zone brick p0 18 0 0 p1 24 0 0 p2 18 5 0 p3 18 0 8 size 9 10 8  group  2
gen zone brick p0 0 0 8 p1 18 0 8 p2 0 5 8 p3 0 0 18 p4 18 5 8 p5 0 5 18 p6 8 0 18 p7 8 5 18 size 17 10 17 ratio 1.03 1 1  group 3
group 0 range group 1 any group 2 any group 3 any
;====================================
;设置边界条件
fix x y z range z -0.1 0.1
fix x range x 23.9 24.1
fix x range x -0.1 0.1
fix y range y -0.1 0.1
fix y range y 4.9 5.1
;====================================
;初始地应力的生成
model elas
prop density 2000 bulk 3e9 shear 1e9
set gravity 0 0 -10.0
solve
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
save 初始应力.sav
;====================================
;实际参数求解
model mohr
prop density 2000 bulk 1.0e8 shear 3.0e7 coh 11500.0 tens 1.0e7  fric 20 dila 30
;plot fos cont ssi outline on vel red
solve
save 支护前.sav

; =======================================================
sel cable id 1 beg 16 2.5 13 end 12 2.5 11 nseg 4
sel cable id 1 beg 12 2.5 11 end 6 2.5 8 nseg 6
sel cable id 1 prop emod 200e9 xcarea 3.8e-3 ytens 1.414e5 gr_k 1 gr_coh 1 gr_fric 1 gr_per 0.314 range cid 1 4
sel cable id 1 prop emod 200e9 xcarea 3.8e-3 ytens 1.414e5 gr_k 7e6 gr_coh 1e2 gr_fric 25 gr_per 0.314 range cid 5 10
; =======================================================
sel pile id=1 begin=(16 2.5 13) end=(16 2.5 4) nseg=9
sel pile prop Emod=8.0e10 Nu=0.30 XCArea=0.7854 XCJ=9.82e-2 XCIy=4.91e-2  XCIz=4.91e-2 Per=3.14 CS_sK=1.3e11 CS_sCoh=0.0 CS_sFric=10.0
&
CS_nK=1.3e11 CS_nCoh=0.0 CS_nFric= 0.0 CS_nGap=off
sel delete link range id 1
sel link id 1 1 target node tgt_num 12
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 1  
sel cable id 1 pretension 10e3 range cid 1 4 ;在自由段施加预应力
; =======================================================
sel node fix y z range id 12
sel node init xvel=0.707e-8 yvel=0.707e-8 range id 12 ;
plot con sxx ou on
;plot con ssi ou on vel red ;剪切应变增量及速度矢量图
;plot con xd magf 5
;plot set mag 2.8
plot add sel cable stress sk   ;锚杆应力分布图
plot add sel pile coupling stress normal  ;桩法向受力图
hist id=2 sel node xdisp id=1
hist id=3 sel cable force cid 1
hist id=4 sel pile  force  fy cid 2
hist id=5 sel pile moment mz(16 2.5 6)
solve
我的问题是：为什么出现could not find pilesel？桩怎么找不到啊？可是我用plo sel geo 看的时候是可以看到建立好的桩。

问题是hist id=4 sel pile  force  fy cid 2上，你plot一下CID号就知道了，cable的cid是从1到10，pile的cid是11到19，你的cid号指定错了，pile没有cid 2。

谢谢@wqq8889   ，我知道了，这个问题已经解决了。可是还有个问题，就是我改过之后，加上锚索和抗滑桩后为什么会和没加之前的效果一样？位移和应力几乎没什么变化，这是什么原因？

修改之后的命令流为：
new
gen zone brick p0 0 0 0 p1 18 0 0 p2 0 5 0 p3 0 0 8 size 17 10 8   group 1    
gen zone brick p0 18 0 0 p1 24 0 0 p2 18 5 0 p3 18 0 8 size 9 10 8  group  2
gen zone brick p0 0 0 8 p1 18 0 8 p2 0 5 8 p3 0 0 18 p4 18 5 8 p5 0 5 18 p6 8 0 18 p7 8 5 18 size 17 10 17  group 3
group 0 range group 1 any group 2 any group 3 any
plo blo gro
;====================================
;设置边界条件
fix x y z range z -0.1 0.1
fix x range x 23.9 24.1
fix x range x -0.1 0.1
fix y range y -0.1 0.1
fix y range y 4.9 5.1
;====================================
;初始地应力的生成
model elas
prop density 2000 bulk 3e9 shear 1e9  range group 0
set gravity 0 0 -10.0
solve  
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zdisp 0
model mohr
prop density 2000 bulk 1.0e8 shear 3.0e7 coh 12100 tens 1.0e6  fric 20 dila 20  range group 0
plot cont ssi outline on
solve  
save 初始应力.sav
;====================================
sel pile id=1 begin=(13 2.5 13) end=(13 2.5 4) nseg=12
sel pile prop  Emod=8.0e10 Nu=0.30 XCArea=12.56 XCJ=0 XCIy=0  XCIz=0 Per=3.14 CS_sK=1.3e11 CS_sCoh=10e3 CS_sFric=10.0 &
CS_nK=1.3e11 CS_nCoh=10e3 CS_nFric= 20 CS_nGap=off
plo con ssi
solve
; =======================================================
sel cable id 1 beg 13 2.5 13  end 4 2.5 9 nseg 8
sel cable id 1 prop emod 200e9 xcarea 5.54e-3 ytens 15.34e5 gr_k 0 gr_coh 0 gr_fric 0  range cid 1 3
sel cable id 1 prop emod 200e9 xcarea 5.54e-3 ytens 30e5 gr_k 0.56e9  gr_coh 1.5e10 gr_fric 25 gr_per 0.264 range cid 4 8
sel delete link range id 1
sel link id 1 1 target node tgt_num 14
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 3  
sel cable id 1 pretension 10e3 range cid 1 3 ;在自由段施加预应力
solve

上面那个sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 3中的id 3 应为id1

谢谢楼主，好东西。

这样抗滑桩还是没起作用啊