摘要
近年来,由于我国高层建筑的发展,桩筏基础被广泛采用。但由于桩筏基础是一个十分复杂的系统,许多方面的问题尚未解决。因此,发展桩土相互作用的实用、有效和可靠的理论和数值分析方法,进一步认识桩筏与地基土的共同作用机制,充分考虑筏板直接对地基土的荷载传递等因素,从而改进桩筏基础的传统设计方法,仍是目前迫切需要解决的课题。
ANSYS软件在应用于岩土工程领域时,对岩土等颗粒性材料,其程序内部仅提供了近似的Drucker-Prager材料模型,是一种理想的弹塑性模型,不能够很好地反映岩土材料的非线性本构关系,尚无国内岩土工程领域中广泛采用的反映岩土材料的非线性本构关系的Duncan-Chang模型。这一缺憾,无疑使ANSYS软件在土工有限元分析问题中的应用范围和适用性受到了局限。为弥补这一不足,本文首先基于ANSYS的用户可编程特性(UPFs )的二次开发功能,成功导入了Duncan-Chang非线性本构模型,并且对二次开发的成果进行了一系列的验证和比较。计算结果表明,模型算例的求解效率和精度令人满意,有望成为强大的解决大型复杂土工问题的非线性数值分析工具。
与此同时,由于桩土相互作用问题的复杂性,要获得解析解是十分困难的,所以利用非线性有限元法对该问题进行数值分析是非常必要的。本文建立了桩土相互作用的三维有限元-接触面单元的非线性耦合数值模型来模拟桩与土的几何、力学特性,基于ANSYS二次开发上编制了桩土相互作用的三维非线性有限元分析程序PSIFEM。利用该程序作为本课题研究的基本工具,对竖向荷载作用下的单桩、群桩与土相互作用机理进行了探讨,在计算中,主要考虑不同桩数、桩径、桩间距、桩长径比及承台刚度等条件下对桩筏基础体系的受力和变形影响,提出了桩筏基础设计的改进建议。
最后,结合本文所编的ANSYS二次开发程序,将本文的计算理论在天津的实际工程中进行了应用,取得了较好的效果。
In recent years, the piled raft foundation has been applied widely because of the development of high buildings in our country. But for the reason that the piled raft foundation is a very complicated system, there are many aspects of problems that have not been resolved. Therefore, it is an urgent task to develop effective and reliable theories and numerical analysis methods for the pile-soil interaction, understand the mechanism of the interaction between the piled raft foundation and the soil deeply, make clear about the load-transfer factors of the raft-soil and improve the traditional designing method of the piled raft foundation.
In ANSYS, only the Drucker-Prager model is provided for soil simulation. As considering soil as a perfect elastoplastic material, this model can only give an approximate simulation. With a nonlinear constitutive law, the Duncan-Chang model can make a good representation of soil, but it has not been implemented into ANSYS yet, which may limit ANSYS to more geotechnical research fields. Therefore, one purpose of this work is to implement the Duncan-Chang model into ANSYS based on the User Programmable Features (UPFs) of ANSYS. The compatibility between this implemented model and the platform of ANSYS has been verified by some numerical analysis cases, and the efficiency and accuracy of this program have been manifested by validation case studies. It could be expected that this program may be served as a powerful tool for complex geotechnical problems solving.
Meanwhile, the problem of the pile-soil interaction is so complex that it is impossible to gain an analytical result. Thus, it is necessary to employ the nonlinear finite element method to analyze this problem. In this paper, a three-dimensional nonlinear numeric model of the coupled finite element-contact interface element is established to simulate the geometrical and mechanical characteristic of piles and the soil. A program named PSIFEM, which is used to analyze the pile-soil interaction, is compiled with the secondary development of ANSYS. Based on this program, the interaction mechanism between soil and single pile or group piles under vertical load was studied. The influence of different parameters on the stress and deformation of piled raft foundation was investigated, which include the number of piles, pile diameter, pile spacing, ratio of pile-length to pile-diameter and the stiffness of cap. Several pieces of improvement advice about the optimization design of piled raft foundations are suggested.
Finally, the developed program has been applied to practical engineering problems in Tianjin, and good performance has been achieved.
引文
[1] Clancy, P. & Randolph, M. F.,A contribution to optimum design of piled rafts, Geotechnique, 1998,48(3): 301~317
[2] Patka,J.& Naborczyk, J., Bearing capacity of piles with the spot footing , Proc. of nineth ICSMFE,1977,1:683~686
[3] Akinmusuru,J.O.,Interaction of piles and cap in piled footings,Proc. of ASCE,J. of Geotech.Engng.Div.,1980,106(11):1263~1268
[4] 周福田,低桩基础的垂直承载力,岩土工程学报,1984, 6(4): 86~94
[5] 杨克己等,基础-桩-土共同作用的特性与承载力研究,岩土工程学报,1988,No.1
[6] Koirumi,Y.& Iro,K.,Field tests with regard to pile driving and bearing capacity of piled foundations,Soils and Foundation, JSSMFE,1967, 7(3)
[7] Seed, H.B.& Reese, L.C.,The action of soft clay along friction piles, Proc.ASCE , 1955,Vol.81
[8] Poulos,H.G. & Davis,E.H., Pile foundation analysis and design, New York: Wiley, 1980
[9] Ta, L.D.& Small, J.C., Analysis of pole raft systems in layered soils, Int. J. Numer. Analytic . Meth . Geometh ., 1996,20
[10] 费勤发,马海龙,分层总和法在桩非线性分析中的应用-变形总和法,上海力学,1997
[11] Cooke, R.W., The settlement of friction pile foundations, Proc. Conf. Tall Buildings , Kuala Lumpur, 1974,7~9
[12] Randolph,M.F. and Wroth,C.P., Analysis of deformation of vertically loaded piles , J. Geotech .Eng. Div , ASCE,1978, 104(12):1465~1488
[13] Cooke,R.W.,et al., Jacked piles in London clay: a study of load transfer and settlement under working condition, Geotechnique , 1979, 29(2):113~147
[14] 宰金珉,杨嵘昌,桩周土非线性分析的广义剪切位移法,南京建筑工程学院学报,1993,No.1
[15] Vesic,A.S.,史佩栋译,桩基设计原理,杭州:浙江省建筑科学研究所,1982
[16] Ellison,R.D.、Dappolonia,E,G.R., Load-deformation mechanism for bored piles, Proc.J.Soil Mechanics and Foundation Div.ASCE,1971,97(4):661~678
[17] 陈雨孙,周红,纯摩擦桩荷载-沉降曲线的数值计算方法,岩土工程学报,1987,9(2):661~678
[18] Ottaviani,M.,Three dimensional finite element analysis of vertically loaded pile groups , Geotechnique , 1975, 25(2): 159~174
[19] 倪新华,筏基-桩基-土体共同作用的数值分析[博士学位论文] ,上海:同济大学,1990
[20] Butterfield R.& Banerjee P.K., The problem of pile-group-cap interaction, Geotechnique ,1971, 21(2):135~142
[21] Banerjee,P.k. & Davis ,T.G., Analysis of pile groups embedded in Gibson soil , Proc, 9th .Int. Conf. Soil Mech . Fdn Engng , Tokyo 1, 1977,381~386
[22] Poulos, H.G., Pile behavior-theory and application , Geotechnique, 1989, 39(3): 365~415
[23] Kuwabara,F., An elastic anaylsis for piled raft foundations in a homogeneous soil, Soils and Foundations,1989,28(1),82~92
[24] Poulos, H.G., Piled raft in swelling or consolidated soils, Proc. J. Geot. Eng. Div. ASCE,1993,119(2):374~380
[25] 张问清,赵锡宏,宰金珉,任意力系作用下层状弹性半空间的有限层分析方法,岩土工程学报,1981, 3(2):27~42
[26] 张崇文等,有限层有限元混合法研究桩土相互作用,天津大学学报,1995,26(6)
[27] 王玉杰,严宗达,有限棱柱法在桩土相互作用分析中的应用,天津大学学报,1997,30(2)
[28] 水利部西北水利科学研究所试验中心,心墙料坝基砂静力特性试验报告,2002 年 1 月
[29] 冯卫星,常绍东等,北京细砂土邓肯-张模型参数试验研究,岩土力学与工程学报,1999,18(3):327~330
[30] Cooke, R.W., Piled raft foundation on stiff claysa-a contribution to design philosophy , Geotechnique , 1986, 36(2): 169~203
[31] Hooper, J.A., et. al, Observations on the behaviour of a piled-raft foundation on Londaon clay,Proc. of Institute of Civil Engineers,55,Dec. 1973,855~877
[32] Chow, Y. K.. Analysis of vertically loaded pile groups, Int. J. Numer. Analytic . Meth . Geometh ., 1986, 10(1): 59~72
[33] Chow . Y. K., Three-dimensional analysis of pile groups , J. Geotechnical Eng. 1987, 113(6): 637~651
[34] Mandolini, A.&Viggiani, C., Settlement of piled foundations,Geotechnique ,1997, 47(4): 791~816
[35] Kuwabara, F., Settlement behavior of non–linear soil around single piles subjected to vertical loads, Soils and Foundations , 1991,31(1): 39~46
[36] 雷晓燕,三维接触问题新模型研究,土木工程学报,1996, 29(3): 24~33
[37] 徐斌,王大通,高大钊, 群桩沉降验算中接触单元模型应用的若干问题,同济大学学报,1998, 26(2):149~152
[38] 陈仁朋,梁国钱,考虑桩土相互滑移的单桩和群桩的非线性分析,浙江大学学报(工学版),2002,11: 668~673
[39] 殷宗泽,朱泓,许国华, 土与结构材料接触面的变形及其数学模拟,岩土工程学报,1994,16(3): 14~21
[40] 于丙子,陈连礼,叶伟泉,三维空间的节理单元,岩土工程学报,1983,15(3): 1~11
[41] 高文华,杨林德,叶为民,板—土三维等厚接触单元模型及其应用,同济大学学报,1999,27(3): 297~300
[42] 雷晓燕,Swoboda G.,杜庆华,接触摩擦单元的理论及其应用,岩土工程学报,1994,16(3): 23~31
[43] 简政,刘用暖,桩-土三维等厚度接触单元模型及其工作性状,水利学报,2002,11: 113~118
[44] 高俊合,于海学,赵维炳,土与混凝土接触面特性的大型单剪试验研究及数值模拟, 土木工程学报,2000, 33(4): 42~46
[45] 张冬霖,卢廷浩,一种土与结构接触面的建立及其应用,岩土工程学报,1998, 21(3): 222~226
[46] Goodman , R.E. ,Taylor , R.L. & Breekle ,T.L., A model for the mechanics of joints and interfaces, J. Soil Meth . Found . ASCE,94,1968
[47] Duncan J. & Chang C.Y., Nonlinear analysis of stress and strain in soils, J.SMFD, ASCE,1970,1629~1653
[48] Clough, G.W. & Hansen, L.A., Clay anisotreopy and braced wall behavior, J. Geotech.Engneering.Div.1981,107(7):893~914
[49] Desai C.S.,Zaman M.M., Thin Layer Element for Interfaces and Joints,Int Journ.For Num.&Analy Meth in Geomech,1984,8(1):19~43
[50] 朱百里,沈珠江等,计算土力学,上海:上海科学技术出版社,第一版,1990
[51] 刘金砺,袁振隆,章国平,钻孔群桩工作机理及承载力研究,科学建筑研究报告,1984,No.2-1
[52] 刘金砺,黄强,竖向荷载下群桩变形性状及沉降计算,岩土工程学报,1995,17(6):1~13
[53] 宰金珉,宰金璋,高层建筑基础分析与设计,北京:中国建筑工业出版社,第一版,1993
[54] 宰金珉,杨嵘昌,群桩与土和承台的非线性共同作用分析,南京建筑工程学院学报,1993,No.2: 1~14
[55] 郑俊杰,彭小荣,桩土共同作用设计理论研究,岩土力学,2003,24(2):241~245
[56] 郑颖人,龚晓南,岩土塑性力学基础,北京:中国建筑工业出版社,1989
[57] 龚晓南,土塑性力学,杭州:浙江大学出版社,1991
[58] 王瑁成,邵敏,有限单元法基本原理和数值方法,北京:清华大学出版社,第二版,1997
[59] 陈环,叶国良,关于桩土相互作用问题,港口工程,1990,No.5: 1~6
[60] 张学言,岩土塑性力学,北京:人民交通出版社,1993
[61] 钱家欢,殷宗泽,土工原理与计算,北京:中国水利水电出版社,第二版,1996
[62] 罗定安,工程结构数值分析方法与程序设计,天津:天津大学出版社,1995
[63] 华南理工大学等四院校编,地基与基础,北京:中国建筑工业出版社,1981
[64] 朱伯芳,有限单元法原理与应用,北京:中国水利水电出版社,第二版,1998
[65] 龙驭球,有限元法概论,北京:人民教育出版社,1978
[66] 卢世深,林亚超编译,桩基础的计算和分析,北京:人民交通出版社,1987
[67] 胡汉兵,余祖明,杨升威,竖向荷载下群桩承载特性的弹塑性分析,长江科学院院报,2001,18(1): 41~44
[68] 高大钊,赵春风,徐斌,桩基础的设计方法与施工技术,北京:机械工业出版社,1999
[69] 王晓东,魏道垛,宰金珉,群桩—土—承台结构共同作用的数值分析,岩土工程学报,1996, 18(4): 27~33
[70] 张保良,赵锡宏等,上部结构-筏-桩-地基共同作用分析,建筑结构学报,1997, 18(2): 72~78
[71] 杨敏,艾智勇,桩土相互作用理论研究与按沉降控制设计桩基础的工程实践,第六届全国岩土力学数值分析与解析方法讨论会,中国,广州,1998,43~59
[72] 刘前曦,侯学渊,筏-桩-土共同作用分析方法,同济大学学报,1996,24(6): 625~630
[73] 李宁,韩洹,单桩复合地基加固机理数值试验研究,岩土力学,1999,20(4): 42~48
[74] 杨敏,王树娟,桩筏基础相互作用下土中应力场的变化规律,岩土工程学报,1999,21(1): 26~30
[75] 左名麒,胡人礼,毛洪渊,桩基础工程设计施工监测,北京:中国铁道出版社,1995
[76] 桩筏基础在上下部共同作用下的极限分析,严平,龚晓南,土木工程学报,2000,33(2): 87~95
[77] 赖伟,熊峰,群桩-土-承台结构共同作用的数值分析,四川大学学报(工学科学版),2001,33(4): 17~20
[78] 周翠英,谭详韶等,动静荷载综合作用下桩土共同作用的机理研究,中山大学学报(自然科学版),2001,40(4): 97~101
[79] 杨进,彭苏萍,群桩条件下桩土相互作用实验研究,岩土力学,2004,25(2): 312~315
[80] 张社荣,李光华,竖向荷载下桩土相互作用的有限元分析,中国农村水利水电,2005,12: 51~54
[81] 袁凡凡,栾茂田,筏板-桩-土相互作用的三维弹塑性有限元分析,岩石力学与工程学报,2005,24(18): 3332~3336
[82] 张建辉,邓安福,桩筏基础的新型分析方法,土木工程学报,2002,35(4): 102~108
[83] 林天健,熊厚金,桩基础设计指南,北京:中国建筑工业出版社,1999
[84] 黄强,桩基工程若干热点技术问题,北京:中国建材工业出版社,1996
[85] 刘利民,舒翔等,桩基工程的理论进展与工程实践,北京:中国建材工业出版社,2002
[86] П.Г.库兹明,清华大学工程地质及基础工程教研组译,土力学讲义,北京:高等教育出版社,1959
[87] 陈希哲,土力学地基基础,北京:清华大学出版社,第三版,1997
[88]孙钧,汪炳鉴,地下结构有限元法解析,上海:同济大学出版社,第二版,1988
[89]陈仲颐, 叶书麟,基础工程学,北京:中国建筑工业出版社,1990
[90] 杨敏,减少沉降桩基础的设计理论和工程应用,岩土工程青年专家学术论坛文集,北京:中国建筑工业出版社,1998,129~136
[91] 黄绍铭,裴捷,贾宗元等,软土中桩基沉降计算,中国土木工程学会第四届土力学及基础工程学术会议论文选集,北京:中国建筑工业出版社,1986,237~243
[92] 黄绍铭,王迪民等,减少沉降量桩基的设计与初步实践,中国土木工程学会第六届土力学及基础工程学术会议论文集,北京:中国建筑工业出版社,1991,405~410
[93] 刘金砺,迟铃泉,桩土变形计算模型和变刚度调平设计,岩土工程学报,2000,22(3):151~157
[94] 管自立,疏桩基础设计实例分析与探讨(一),建筑结构,1993(10):26~31
[95] 管自立,疏桩基础设计实例分析与探讨(二),建筑结构,1993(11):42~46
[96] 赵锡宏,董建国等,桩基减少桩数与沉降问题的研究,土木工程学报,2000,33(6):71~74
[97] 董建国,赵锡宏,高层建筑地基基础-共同作用理论与实践,上海:同济大学出版社,1997,292~314
[98] Poulos,H.G.&Davis,E.H.编,铁天石等译,桩基础的分析和设计,北京:中国建筑科学研究院科技资料交流部,1991
[99] 刘金砺等,桩基工程技术,北京:中国建材工业出版社,1996
[100] 刘金砺,桩基础设计与计算,北京:中国建筑工业出版社,1990
[101] 暴雨,桩土相互作用的三维非线性有限元分析[硕士学位论文],天津:天津大学,2002
[102] 陈明中,群桩沉降计算理论及桩筏基础优化设计研究 [博士学位论文] ,杭州:浙江大学,2000
[103] 刘利民,桩基础内力分析方法的研究[博士学位论文] ,上海:同济大学,1997
[104] 《桩基工程手册》编写委员会,桩基工程手册,北京:中国建筑工业出版社,1995
[105] 史佩栋,实用桩基工程手册,北京:中国建筑工业出版社,1999
[106] 中华人民共和国建设部,国家质量监督检验检疫总局,建筑地基基础设计规范(GB50007-2002),北京:中国建筑工业出版社,2002
[107] 中华人民共和国建设部,建筑地基基础设计规范(GBJ7-89),北京:中国建筑工业出版社,1989
[108] 中国建筑科学研究院,建筑桩基技术规范(JGJ94-94),北京:中国建筑工业出版社,1995
[109] 唐章宏,薛赛男,冯峰等,Visual Fortran 程序设计,北京:人民邮电出版社,2000
[110] 邓巍巍,王越男,Visual Fortran 编程指南,北京:人民邮电出版社,2000
[111] 王国强,实用工程数值技术及其在 ANSYS 上的实践,西安:西北工业大学出版社,2000
[112] 陈精一,蔡国忠,电脑辅助工程分析:ANSYS 使用指南,北京:中国铁道出版社,2001
[113] 嘉木工作室,ANSYS 有限元分析教程,北京:机械工业出版社,2002
[114] 刘涛,杨凤鹏,精通 ANSYS,北京:清华大学出版社,2002
[115] 博嘉科技,有限元分析软件:ANSYS 融汇与贯通,北京:中国水利水电出版社,2002
[116] 李权,ANSYS 在土木工程中的应用,北京:人民邮电出版社,2005
[117] ANSYS Inc., ANSYS theory reference, 1994
[118] ANSYS Inc., Guide to ANSYS user programmable features, 1994
[119] ANSYS Inc., Guide to Interface with ANSYS, 1994
[120] ANSYS Inc., The APDL Programmer’s Guide, 1994
[121] ANSYS Inc., ANSYS 技术报告,北京,1998
[122] ANSYS Inc., ANSYS APDL 使用指南,北京,2000
[123] ANSYS Inc., ANSYS 结构非线性分析培训资料,北京,2000
[124] ANSYS Inc., ANSYS 单元手册摘要,北京,2000