双排桩支护出现的弯矩等问题,恳请各位给予指导
本人做一个双排桩的基坑支护,现在在结果分析中遇到如下几个问题:1)在表示沿深度方向的弯矩时候应该是Mx还是Mz(局部和全局的问题)
2)在我的想法中,应该是每排桩的那几根桩的受力和所受到的弯矩应该是同样分布的或者是对称分布的,这个Mx,Mz图离我的想法有出入.
3)我这里面的桩单元和梁单元没有设置link,以前曾经试过link,但是不能够计算,应该怎么设置那?非常希望各位能够给我以提示,为感! 程序如下
new
;输出到文本
set log on
set logfile liangplot.log
plot sk
;建立土体模型
gen zone brick size 30 10 30
group clay
mod mohr range group clay
pro bulk 8e7 shear 3.8e7 coh 30000 fric 20 range group clay
ini dens 1230 range group clay
;设置边界条件和初始应力
fix x range x -.1 .1
fix x range x 29.9 30.1
fix y range y -.1 .1
fix y range y 9.9 10.1
fix z range z -0.1 0.1
ini szz -369000 grad 0 0 12300 range z 0 30
ini sxx -158130 grad 0 0 5271range z 0 30
ini syy -158130 grad 0 0 5271range z 0 30
set gravity 0 0 -10.0
solve
save chushi.sav
;设置结构单元
sel pile id 1 begin 10 0 30 end 10 0 13 nseg 17
sel pile id 1 begin 11.2 0 30 end 11.2 0 13 nseg 17
sel pile id 1 begin 10 2.5 30 end 10 2.5 13 nseg 17
sel pile id 1 begin 11.2 2.5 30 end 11.2 2.5 13 nseg 17
sel pile id 1 begin 10 5 30 end 10 5 13 nseg 17
sel pile id 1 begin 11.2 5 30 end 11.2 5 13 nseg 17
sel pile id 1 begin 10 7.5 30 end 10 7.5 13 nseg 17
sel pile id 1 begin 11.2 7.5 30 end 11.2 7.5 13 nseg 17
sel pile id 1 begin 10 10 30 end 10 10 13 nseg 17
sel pile id 1 begin 11.2 10 30 end 11.2 10 13 nseg 17
sel nod id 501 10.6 0 30
sel nod id 502 10.6 2.5 30
sel nod id 503 10.6 5 30
sel nod id 504 10.6 7.5 30
sel nod id 505 10.6 10 30
sel nod id 601 10 1.25 30
sel nod id 602 10 3.75 30
sel nod id 603 10 6.25 30
sel nod id 604 10 8.75 30
sel nod id 701 11.2 1.25 30
sel nod id 702 11.2 3.75 30
sel nod id 703 11.2 6.25 30
sel nod id 704 11.2 8.75 30
def set_beam
loop ii (1,5)
iid1=1+36*(ii-1);节点号
iid2=500+ii
iid3=19+36*(ii-1)
command
selbeamselid 2nodes iid1iid2 ;刚性梁
selbeamselid 2nodes iid2iid3
endcommand
endloop
end
set_beam
def qianliang
loop ii (1,4)
iid1=1+36*(ii-1) ;节点号
iid2=600+ii
iid3=37+36*(ii-1);节点号
command
selbeamselid 3nodes iid1 iid2 ;前梁
selbeamselid 3nodes iid2 iid3
endcommand
endloop
end
qianliang
def houliang
loop ii (1,4)
iid1=19+36*(ii-1);节点号
iid2=700+ii
iid3=55+36*(ii-1);节点号
command
selbeamselid 3nodes iid1 iid2 ;后梁
selbeamselid 3nodes iid2 iid3
endcommand
endloop
end
houliang
sel pile id = 1 prop emod 3.0e10 nu 0.16 xcarea 0.5024 dens 2700 &
xciy 0.001256xciz 0.001256 xcj 0.002512&
per 2.512 &
cs_sk0.11e9 cs_scoh 40000 cs_sfric 0 &
cs_nk0.67e9 cs_ncoh 20000 cs_nfric 0.0 cs_ngap off
sel beam id = 2 prop emod 3.0e10 nu 0.16 xcarea 0.5024 dens 2700 &
xciy 0.001256xciz 0.001256 xcj 0.002512
sel beam id = 3 prop emod 3.0e10 nu 0.16 xcarea 0.5024 dens 2700 &
xciy 0.001256xciz 0.001256 xcj 0.002512
;是否还需要设置link???
;pr gp pos
;pr zone gp
;pr zone str
ini xdis 0 ydis 0 zdis 0
solve
save jiazhuang.sav
;开挖模拟
ini xdis 0 ydis 0 zdis 0
model null range x 0 10 y 0 10 z 20 30
hist unbal
hist write file zuida.dat
solve
print hist
save liang.sav
set log off 做了修改之后发现弯矩比较对称了,但是还是不怎么合理,
另还想请教下:桩的正负弯矩是否是和桩的弯曲方向有关,即桩的水平位移.因为弯矩平时都不强调正负性, 你桩的弹簧刚度取得对吗?有什么根据吗?
我最近也做一个桩锚支护的例子?
有机会QQ聊一下? 不好意思,刚度取得不是很科学,只是综合了别人的数据,自己再调试了下,我觉得现在最大的问题就是link的问题,桩和beam的连接,这个问题一直没有解决,望请各位能够给予帮助下,谢谢,7附:本人的qq:304975361 设置了连梁和桩的link之后,整个Mx,Fx,Fz都为零了,急盼各位给予解答,毕业在即,还望各位帮忙! 桩沿深度的弯矩看到某实例中用的是My,通过修改弯矩的值能够比较好的吻合,现在反而是位移的值不太合理了,差了几倍... 看来楼主也是个高手。
页:
[1]