请教:如何模拟压力分散型锚索

liqimin85

我最近在模拟压力分散型锚索，可是结果不理想，请各位高手指教下！谢谢
gen zone radcyl p0 0 0 0 p1 0.4 0 0 p2 0 5 0 p3 0 0 0.4 &
                dim 0.065 0.065 0.065 0.065 size 5 10 4 5 fill
gen zon reflect dip 90 dd 90 orig 0 0 0
gen zon reflect dip 0 orig -0.4 0 0
mo mohr
pro bulk 2.6e9 shear 1.2e9 fric 19 coh 52.4e3 ten 5e9
fix x range x -0.39 -0.41
fix x range x 0.39 0.41
fix y range y 4.99 5.01
fix z range z -0.39 -0.41
fix z range z 0.39 0.41
set grav 0 0 -10
ini density 1900
;solve
ini xdisp 0 ydisp 0 zdisp 0
model null range cyl end1 0 0 0 end2 0 4.5 0 rad 0.065
;围岩
interface 1 face range cylinder end1 (0,0,0) end2 (0,4.5,0) radius .067 &
            cylinder end1 (0,0,0) end2 (0,4.5,0) radius .063 not
interface 1 prop kn 1.21e10 ks 1.21e10 fric 30 coh 30000 bslip on
gen zone cylinder p0 0 6 0 p1 0.065 6 0 p2 0 10.5 0 p3 0 6 0.065 size 5 9 4
；接触面
gen zon reflect dip 90 dd 90 orig 0 6 0 range y 6 10.5
gen zon reflect dip 0 orig -0.065 6 0 range y 6 10.5
sel liner id 1 elemtype dkt_cst range y=10.49 10.51
delete range y 6 10.5
sel node ini ypos add -6 range y 6 10.5
sel liner PROP iso=( 25e9, 0.15) thick=0.1
sel liner PROP cs_nk=8e8 cs_sk=8e8 &
               cs_ncut=0.0 cs_scoh=0.0 cs_scohres=0.0 cs_sfric=0.0
；承压板
sel cable id 1 begin 0.04,0,0  end 0.04,4.5,0 nseg 5
sel cable id 2 begin -0.04,0,0  end -0.04,4.5,0 nseg 5
sel cable id 1 pro emod 1.95e11 xcarea 1.632e-3 ytens 186e4 &
              gr_per 0.162024 gr_k 0 gr_c 0
sel cable id 2 pro emod 1.95e11 xcarea 1.632e-3 ytens 186e4 &
              gr_per 0.162024 gr_k 0 gr_c 0
sel cablesel nodes 6 83
sel cablesel nodes 31 88
；锚索
gen zone cylinder p0 0 6 0 p1 0.065 6 0 p2 0 10.5 0 p3 0 6 0.065 size 5 9 4
gen zon reflect dip 90 dd 90 orig 0 6 0 range y 6 10.5
gen zon reflect dip 0 orig -0.065 6 0 range y 6 10.5
group grout1 range y 6 10.5
ini y add -6 range group grout1
group rock range group grout1 not
；注浆体
mo mohr
pro bulk 2.6e9 shear 1.2e9 fric 19 coh 52.4e3 ten 5e9 range group rock
pro bulk 1.44e10 shear 1.2e10 fric 32.5 coh 2.0e6 ten 6e9 range group grout1
interface 1 prop kn 1.21e10 ks 1.21e10 fric 20 coh 30000 bslip on
sel cable id=1 pretension 15e3 range cid 149,145
sel cable id=2 pretension 15e3 range cid 150,154                                          
sel delete link range y 0.01 -0.01
sel link net range y 0.01 -0.01
；垫片  
solve

liqimin85

自己顶一下

xunxunyuanyuan

楼主，我也在做这方面的模拟，能不能交流一下
QQ 385637092

xunxunyuanyuan

我按你的命令流算了一下，确实很不理想啊，比我的还不理想·····