|
利用Ansys对某基坑双排门架桩支护剖面计算 其实双排桩的计算,就从结构计算来说,是没什么新意了。只是帖子发了不少,土币没见长,所以也在这里凑合一下,赚点下载费。 我们公司承接了上海某船闸的基坑围护设计,船闸基坑分闸首和闸室,见下图。由于闸首基坑有一面是开口的(和闸室基坑相连),所以对面就没法做对撑,这一面采用了双排桩进行加强。我们是接人家转包给我们的,而且工程较大,人家帮我们出图章,其他部分都有启明星软件出的计算书,启明星已经有了双排桩模块,我们还没买。出章单位一定要看看我们的这个双排桩剖面计算成果,这样能有个底。所以我就用Ansys算了下。 前排桩为间距1000桩径800钻孔桩,后排桩为间距2000间距800钻孔桩,前后排桩上设850厚钢筋混凝土压顶板联系。前排桩腰间还有围囹梁。根据《建筑基坑支护技术规程》推荐排桩支护按m值的弹性地基梁法,我网上还找了下论文,对于前排桩能够形成一个很牢固的门板式的双排桩支护结构,大部分还是把桩间土和被动区土模拟成弹簧。所以模型是一个简单梁-弹簧结构,而且是线弹性静力分析,用Ansys可以方便的建模计算。 计算过程为,从支护体系中截取一个2m宽的单元,单元里包含两根前排桩,一根后排桩,该单元承受的土压力按《规程》计算公式算出,算出的土压力单位是kPa(见附件,这个也用程序算的)。把前后排桩每0.1m长度划分一个梁单元,土压力按宽2m,高0.1m的宽度转化为集中荷载作用在后排桩上,桩间土和被动区土等效成杆单元(即弹簧)。根据m值法,杆单元线刚度随深度线性增大,m值的计算根据《规程》附录C.3.2计算。压顶板也用梁单元模拟。 在Ansys中,桩、压顶板用beam3单元,每0.1m为一个单元。杆用link10单元并设置只受压不受拉,材料属性不变,线刚度增大采用截面积增大法。排桩只约束竖向位移不约束水平位移,土杆在远端约束竖向和水平位移。最后计算出的排桩弯矩和剪力如图。我后来又用类似命令模拟了单排桩的支护,和同济启明星进行了对比,结果弯矩和剪力只相差10%左右而最大位移却相差好几倍,还不知道什么原因。 命令流文件放附件老提示非法,只好考在下面了: /FILENAME,NEWANALYSIS /TITLE,NEWANALYSIS /UNITS,SI /PREP7 ET,1,BEAM3 ET,2,LINK10 !实常1后排钻孔桩,2是前排钻孔桩,3是压顶板,后面是土弹簧 R,1,0.502654816,0.02010619 R,2,1.005309632,0.04021238 R,3,1.70,0.08700104 KEYOPT,1,6,1 KEYOPT,2,3,1 !双排桩间土的杆面积 *DO,I,1,45 R,3+I,I*0.2651264 *ENDDO *DO,I,46,115 R,3+I,I*0.2448864 *ENDDO !被动区土的杆面积 *DO,I,1,4 R,118+I,I*3.0128 *ENDDO *DO,I,5,75 R,118+I,I*2.7828 *ENDDO !1是混凝土,2是钢材(本程序中无用),3是土杆的弹性模量,不是实际土弹模 !(m1=2.3,m2=5.6,m3=1.5,m4=6.1) MP,EX,1,30000000000 MP,PRXY,1,0.3 MP,EX,2,200000000000 MP,PRXY,3,0.3 MP,EX,3,1000000 MP,PRXY,3,0.3 !点个数=(挖深+嵌固深-放坡深)*10 *DO,I,1,116 K,I,0,-0.1*I,0 K,I+200,2.2,-0.1*I,0 *ENDDO *DO,I,41,116 K,I+400,12.2,-0.1*I,0 *ENDDO *DO,I,1,115 L,I,I+1 *ENDDO *DO,I,1,115 L,I+200,I+201 *ENDDO *DO,I,1,116 L,I,I+200 *ENDDO *DO,I,41,116 L,I+200,I+400 *ENDDO !划分单元前要重新选择所有部件 ALLSEL,ALL TYPE,1 MAT,1 REAL,1 ESIZE,0.1 LMESH,1,115,1 TYPE,1 MAT,1 REAL,2 ESIZE,0.1 LMESH,116,230,1 !压顶板 TYPE,1 MAT,1 REAL,3 ESIZE,0.1 LMESH,231 !划分排桩间土弹簧单元 *DO,I,1,115 TYPE,2 MAT,3 REAL,I+3 LESIZE,231+I,,,1 LMESH,231+I *ENDDO !划分被动区土弹簧单元 *DO,I,1,75 TYPE,2 MAT,3 REAL,I+74 LESIZE,347+I,,,1 LMESH,347+I *ENDDO !约束被动土弹簧远端 *DO,I,441,516 DK,I,UX,0 DK,I,UY,0 *ENDDO DK,116,UY,0 DK,316,UY,0 *DIM,EAJK,ARRAY,116 !由于不能从命令流读入土压力文件,所以这两句不输入,要靠界面操作读入土压力 *VREAD,EAJK(1),tuyali (F6.3) !以下继续输入命令 *DO,I,1,116 FK,200+I,FX,EAJK(I) *ENDDO FINISH /SOLU SOLVE !输出内力成果 /POST1 ETABLE,OUTI,SMISC,6 ETABLE,OUTJ,SMISC,12 ETABLE,REFL /TITLE,Axial force diagram PLLS,OUTI,OUTJ,1,0 /image,save,'Axial_force_%T%',jpg !命令结束 空间的图片,有时候链接不能直接显示,可以右键-属性把地址复制过来,用浏览器打开肯定能看到 以下是图片 1、基坑 2、截取一个支护单元 3、模型与土层参数 4、弯矩图 5、剪力图 [ 此帖被jedy54在2010-01-18 12:07重新编辑 ]
|