为什么我的第二排锚杆拉力远小于第一排
本帖最后由 xuchao1122 于 2012-6-12 16:44 编辑new
gen zone brick p0 0 0 0 p1 10 0 0 p2 0 1.6 0 p3 0 0 30 size 10 2 30 ratio 0.9 1 1 group 1
gen zone brick p0 10 0 15 p1 10.6 0 15 p2 10 1.6 15 p3 10 0 30 size 1 2 15 group 2
gen zone brick p0 10 0 0 p1 10.6 0 0 p2 10 1.6 0 p3 10 0 15 size 1 2 15 group 4
gen zone brick p0 10.6 0 0 p1 30 0 0 p2 10.6 1.6 0 p3 10.6 0 30 size 20 2 30 ratio 1.1 1 1 & group 3
pl bl gr
fix x ran x -0.1 0.1
fix x ran x 29.9 30.1
fix z ran z -0.1 0.1
fix y
;初始应力场
model mohr
prop bulk 16.7e6 shear 7.69e6 coh 10e10 tension 1e10 ran z 0 30
ini dens 2000
set grav 0 0 -9.8
hist id=1 unbal
solve
;赋真实参数
prop bulk 16.7e6 shear 7.69e6 coh 0 fric 45 ran z 0 5.4
prop bulk 16.7e6 shear 7.69e6 coh 0 fric 40 ran z 5.4 21.4
prop bulk 16.7e6 shear 7.69e6 coh 0 fric 35 ran z 21.4 27
prop bulk 16.7e6 shear 7.69e6 coh 15e3 fric 20 ran z 27 28.1
prop bulk 16.7e6 shear 7.69e6 coh 8e3 fric 10 ran z 28.1 28.6
prop bulk 16.7e6 shear 7.69e6 coh 0 fric 8 ran z 28.6 30
hist id=2 gp xf 10.5 1 30
hist id=3 gp xf 10.5 1 28
hist id=4 gp xf 10.5 1 26
hist id=5 gp xf 10.5 1 24
hist id=6 gp xf 10.5 1 22
hist id=7 gp xf 10.5 1 20
hist id=8 gp xf 10.5 1 18
hist id=9 gp xf 10.5 1 16
hist id=10 gp xf 10.5 1 15
solve
save 1-1.sav
;桩
model elas range gr 2
prop bulk 66.7e9 shear 50e9 range gr 2
model null ran z 26.4 30 gr 1
hist id=11 gp xdis 10 1 30
hist id=12 gp xdis 10 1 28
hist id=13 gp xdis 10 1 26
hist id=14 gp xdis 10 1 24
hist id=15 gp xdis 10 1 22
hist id=16 gp xdis 10 1 20
hist id=17 gp xdis 10 1 18
hist id=18 gp xdis 10 1 16
solve
save 1-2.sav
;第一排锚索
pl ske
sel cable id=1 beg 10 0.8 26.7 end 15.909 0.8 25.658 nseg 6
sel cable id=1 beg 15.909 0.8 25.658 end 24.772 0.8 24.095 nseg 9
sel cable id=1 prop emod 200e9 ytension 1860e6 xcarea 0.5e-3 gr_coh 1 gr_k 1 &
gr_per 0.628 range cid 1 6
sel cable id=1 prop emod 200e9 ytension 1860e6 xcarea 0.5e-3 gr_coh 3e2 gr_k 2e7 &
gr_per 0.628 gr_fric 25 range cid 7 15
sel delete link range id 1
sel link id=100 1 target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 100
sel cable id=1 pretension 140e3 range cid 1 6
pl add sel geom scale 0.005
model null ran z 22.9 26.4 gr 1
solve
save 1-3.sav
pl sk
pl add sel cable force
;第二排锚索
sel cable id=2 beg 10 0.8 23.2 end 14.924 0.8 22.332 nseg 5
sel cable id=2 beg 14.924 0.8 22.332 end 24.772 0.8 20.595 nseg 10
sel cable id=2 prop emod 200e9 ytension 1860e6 xcarea 0.5e-3 gr_coh 1 gr_k 1 &
gr_per 0.628 range cid 1 5
sel cable id=2 prop emod 200e9 ytension 1860e6 xcarea 0.5e-3 gr_coh 3e2 gr_k 2e7 &
gr_per 0.628 gr_fric 25 range cid 6 15
sel delete link range id 101
sel link id=200 17 target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 200
sel cable id=2 pretension 175e3 range cid 1 5
model null ran z 19.3 22.9 gr 1
solve
save 1-4.sav
是跟我用实体单元模拟桩 有关吗
页:
[1]