地基承载力与基础宽度之间的关系
本帖最后由 xingchen866 于 2011-11-26 13:20 编辑最近在做地基承载力的模拟,为什么我算出来的结果,基础宽度越宽承载力越低啊。我是分级加载的。每一级载荷作用下,监测位移,得出P-S曲线 肯定是错的具体要看你的程序 看了之后才能知道 n
res grav-bal-pla.sav
;************
;归零
;************
ini xdis=0 ydis=0 zdis=0
ini xvel=0 yvel=0 zvel=0
ini state 0
m e range group conctrete
prop bu=b2 sh=k2 range group conctrete
ini dens dens2 range group conctrete
so
pl con szz ou on
;************
;监测与计算
;************
call app_load1.dat
---------------------------------------------------------以上是主计算文件,app_load1.dat文件是以下程序
table 1 name load_settlement
set mech ratio 1.0e-5
pl add table 1 both
range name concrete group conctrete
def add_load
ss=b*0.5
njj=2+2*b ;基础底面节点数
mm = 0
force=(-30.0e3*ss)/njj
ytable(1,1) = 0
xtable(1,1) = 0
;****************
loop n (1,1000)
loop nn (1,3)
mm = mm+1
app_load=mm*(force)
command
app zforce app_load range z h1hd h1hu x bzz byy
solve step 10000
endcommand
;****************
mzdisp1=0
mxdisp1=0
mdisp1 =0
mmm =0 ;检查底面点是否正确
pgt = gp_head
loop while pgt # null
if inrange('concrete',pgt) = 1
gpz = gp_zpos(pgt)
if gpz < h1u
mmm=mmm+1
mzdisp1=mzdisp1+gp_zdisp(pgt)
mxdisp1=mxdisp1+gp_xdisp(pgt)
mdisp1=sqrt(mxdisp1^2+mzdisp1^2)
endif
endif
mzdisp=mzdisp1/njj
mxdisp=mxdisp1/njj
mdisp =mdisp1/njj
pgt = gp_next(pgt)
endloop
;****************
ytable(1,mm+1)=-(1.0*app_load*njj/ss)/1000
xtable(1,mm+1)=-mzdisp
if mech_ratio > 1.0e-5
exit
endif
endloop
endloop
end
add_load
sav result.sav
我现在有个问题,我这种加载方式对么?问题出在哪里,谢谢
本帖最后由 xingchen866 于 2011-12-1 16:15 编辑
miaorenfeng2008 发表于 2011-11-30 18:45 http://forum.simwe.com/static/image/common/back.gif
看了之后才能知道
将基础建在地基里面,不知道用何种方式施加边界条件来模拟刚性基础呢?求教~谢谢,能不能留下QQ方便交流一下?我的是710171305
页:
[1]