桩基沉降计算的理论研究及ANSYS数值仿真分析
|
摘要
目前工程界桩基沉降计算的基本方法为“地基规范”及“桩基规范”，而规范的计算方法理论基础是浅层基础沉降理论，可以说，浅层基础沉降计算的缺点同样带给了桩基沉降计算，从一些文献及本文的分析可以发现，目前桩基沉降计算主要存在如下缺限：1) 浅层基础沉降计算的基本理论是各向同性均质弹性理论，该理论忽略了基础的刚度影响；2) 地基深度确定的不合理；3) 分层土压缩模量E_s的取值问题；4) 现有文献对Geddes解的应用并不合理，无法解决桩侧摩擦阻力呈非线性分布的情况单桩沉降问题。
本文主要针对上述问题，运用弹性理论并结合传统分层总和法对桩基沉降计算进行分析探讨，借助重要数学软件Matlab及有限元计算软件ANSYS对桩基沉降问题进行数值分析和数值模拟分析，得到较多有益于工程实际的结论。
主要工作及结论如下：
1)分析了现行国家规范法计算地基沉降与弹性理论解在理论上的差异以及与工程实测沉降产生偏差的根本原因，阐明了基础宽度、基础埋深、基础刚度、土的泊松比等因素对地基沉降的影响，并证实了基于Mindlin积分解的徐志英解比基于Boussinesq解的积分解合理得多。编制了新的地基变形深度计算程序。
2)当基础尺寸较小、基础刚度较大时，柔性基础与刚性基础的沉降差异较大。
3)通过线积分完整地推导出了桩侧摩阻均匀分布、桩侧摩阻随深度呈线性增加分布以及桩侧摩阻随深度呈任意形状分布的纯摩擦单桩竖向应力及竖向位移的计算公式。
4)在分层总和法基础上提出了一种新的群桩沉降计算模式，并从多方面予以论证，且编制了群桩计算程序。
5)借助通用有限元软件ANSYS模拟分析了摩擦桩、端承桩及摩擦端承桩的沉降性状，ANSYS计算的地基沉降及桩基沉降结果与经典弹性理论计算的结果非常接近，但与弹塑性地基有一定的差距。
At present the essential method for the pile foundation settlement calculation are the ground standard and the pile foundation standard in the engineering field, but the base theory of standard computational method is the settlement theory of shallow layer foundation, we can say, the shortcoming of shallow layer foundation settlement computation has taken to the pile foundation settlement computation similarly, from some literature and this article's analysis we can find, at present the pile foundation settlement calculation mainly existence as follows lacks limits: 1) the elementary settlement calculation theory of shallow layer foundation is the isotropic homogenous theory of elasticity, this theory has neglected the influence of foundation rigidity; 2) the determination of the ground calculation depth is unreasonable; 3) the modulus of compressibility of stratified soil has some problem; 4) the existing literature is unreasonable to the application of Geddes' solution and unable to solve the settlement calculation problem of the single pile when the skin friction is on non-linear distribution.
This article mainly aims at the above question, analyzes and discuss the settlement calculation of the pile foundation by using the theory of elasticity and the traditional layerwise summation method, and carries on the numerical analysis and the value simulation analysis with the aid of the important mathematics software Matlab and the finite element computation software ANSYS, then obtains many conclusion which is beneficial in the actual project.
Main work and result as follows:
1) Through comparison between layerwise summation method and classical elasticity displacement solution, analyze the difference of calculating foundation settlement between the method of Code for Design of Building Foundation and the theory of elasticity ,so does between calculated and observed settlement Expounds the effect of foundation width embedded depth of foundation and Poisson's ratio of soil on foundation settlement, and proves that the Xu Zhi-ying's solution for calculating foundation settlement based on Mindlin's equation is more perfect than the integral solution based on Boussinesq's equation. And a new
本文主要针对上述问题，运用弹性理论并结合传统分层总和法对桩基沉降计算进行分析探讨，借助重要数学软件Matlab及有限元计算软件ANSYS对桩基沉降问题进行数值分析和数值模拟分析，得到较多有益于工程实际的结论。
主要工作及结论如下：
1)分析了现行国家规范法计算地基沉降与弹性理论解在理论上的差异以及与工程实测沉降产生偏差的根本原因，阐明了基础宽度、基础埋深、基础刚度、土的泊松比等因素对地基沉降的影响，并证实了基于Mindlin积分解的徐志英解比基于Boussinesq解的积分解合理得多。编制了新的地基变形深度计算程序。
2)当基础尺寸较小、基础刚度较大时，柔性基础与刚性基础的沉降差异较大。
3)通过线积分完整地推导出了桩侧摩阻均匀分布、桩侧摩阻随深度呈线性增加分布以及桩侧摩阻随深度呈任意形状分布的纯摩擦单桩竖向应力及竖向位移的计算公式。
4)在分层总和法基础上提出了一种新的群桩沉降计算模式，并从多方面予以论证，且编制了群桩计算程序。
5)借助通用有限元软件ANSYS模拟分析了摩擦桩、端承桩及摩擦端承桩的沉降性状，ANSYS计算的地基沉降及桩基沉降结果与经典弹性理论计算的结果非常接近，但与弹塑性地基有一定的差距。
At present the essential method for the pile foundation settlement calculation are the ground standard and the pile foundation standard in the engineering field, but the base theory of standard computational method is the settlement theory of shallow layer foundation, we can say, the shortcoming of shallow layer foundation settlement computation has taken to the pile foundation settlement computation similarly, from some literature and this article's analysis we can find, at present the pile foundation settlement calculation mainly existence as follows lacks limits: 1) the elementary settlement calculation theory of shallow layer foundation is the isotropic homogenous theory of elasticity, this theory has neglected the influence of foundation rigidity; 2) the determination of the ground calculation depth is unreasonable; 3) the modulus of compressibility of stratified soil has some problem; 4) the existing literature is unreasonable to the application of Geddes' solution and unable to solve the settlement calculation problem of the single pile when the skin friction is on non-linear distribution.
This article mainly aims at the above question, analyzes and discuss the settlement calculation of the pile foundation by using the theory of elasticity and the traditional layerwise summation method, and carries on the numerical analysis and the value simulation analysis with the aid of the important mathematics software Matlab and the finite element computation software ANSYS, then obtains many conclusion which is beneficial in the actual project.
Main work and result as follows:
1) Through comparison between layerwise summation method and classical elasticity displacement solution, analyze the difference of calculating foundation settlement between the method of Code for Design of Building Foundation and the theory of elasticity ,so does between calculated and observed settlement Expounds the effect of foundation width embedded depth of foundation and Poisson's ratio of soil on foundation settlement, and proves that the Xu Zhi-ying's solution for calculating foundation settlement based on Mindlin's equation is more perfect than the integral solution based on Boussinesq's equation. And a new
引文
[1] 桩基工程手册编写委员会．桩基工程手册[M]．北京：中国建筑工业出版社，1993
[2] 廖向东，桩基沉降理论的应用研究：[硕士学位论文][D]．上海：同济大学，2001
[3] 陈祥福，沉降计算理论及工程实例[M]．北京：科学出版社，2005
[4] [澳]H．G Poulos；E．H．Davis著，孙幼兰译．岩土力学弹性解[M]．徐州：中国矿业大学出版社，1990
[5] GB 50007-2002 建筑地基基础设计规范[S]．北京：中国建筑工业出版社，2002
[6] JGJ 94-94 建筑桩基技术规范[S]．北京：中国建筑工业出版社，1994
[7] 徐志英．以明德林公式为根据的地基中垂直应力的计算公式[J]，土木工程学报，1957，(4)：485～497
[8] 袁聚云，赵锡宏．竖向均布荷载作用在地基内部时的土中应力公式[J]，上海力学，1995，第16卷第3期：213～222
[9] 王士杰，张梅，张吉占．工程荷载下的明氏应力实用理论解(I)[J]，四川建筑科学研究，2001，第27卷第2期：36～37
[10] 王士杰，张梅，张吉占．工程荷载下的明氏应力实用理论解(Ⅱ)[J]，四川建筑科学研究，2001，第27卷第4期：60～61
[11] 张子明。求解地基位移的影响函数法[J]，河海大学学报，1990，18(2)：25～31
[12] 陈甦．均质地基中的桩基础沉降计算[J]，岩土力学，1996，第17卷第1期：61～70
[13] K．L．Johnson蓍，徐秉业等译．接触力学[M]．北京：高等教育出版社，1992
[14] 陈希哲．土力学地基基础(第四版)[M]．北京：清华大学出版社，2004
[15] 王铁行，赵树德．计算地基沉降分层总和法的缺陷的分析与改进[J]，西安建筑科技大学学报，1996，第28卷第2期：179～182．
[16] Joseph E. Bowles. Foundations analysis and design (the third edition) [M]. NEW YORK: McGraw-Hill Book Company. 1982.
[17] Harry G. Poulos, Edward H. Davis. Elastic solutions for soil and rock mechanics [M]. NEW YORK: John Wiley & Sons, Inc. 1974.
[18] 应永法，泮威风．圆形板基础下地基强度和变形非线性性状研究[J]，岩土力学，2002，第23卷第2期：250～257
[19] 黄绍铭，裴捷，贾宗元等．软土中桩基沉降估算[J]，第四届全国土力学及基础工程学术会议选集，1986，237～243
[20] 杨敏，王树娟，王伯钧等．使用Geddes应力系数公式求解单桩沉降[J]，同济大学学报，1997，第25卷第4期：379～385
[21] 肖拥军，刘新华，杨荣丰．砂土地基桩基础承载性状试验研究[J]，湖南科技大学学报(自然科学版)，2004，第19卷第2期：60～63
[22] 高强，杨龙才，王炳龙．软粘土钻孔灌注桩承载特性的试验研究[J]，岩土工程技术，2004，第18卷第1期：35～39
[23] 钟闻华，石名磊，刘松玉．超长桩荷载传递性状研究[J]，岩土力学，2005，第26卷第2期：316～318
[24] 王东红，谢星，张炜等．黄土地区超长钻孔灌注桩荷载传递性状试验研究[J]，工程地质学报，1004-9665／2005／13(01)-0117-07
[25] 钟龙辉．西安地区高层建筑桩基沉降计算[J]，勘察科学技术，2001，第3期：9～12
[26] 单玉川、郑学斌．论桩基的沉降量[J]，建筑结构，2001，第31卷第10期：26～Z7
[27] 李云翼、陈国栋、甘德福等．桩基最终沉降量的Gan-Chen模式计算方法[J]，工程勘察，2002，第6期：4～6
[28] 陈仁朋．软弱地基中桩筏基础工作性状和分析设计方法研究：[博士学位论文][D]．杭州：浙江大学，2001
[29] 陈国兴，金永彬，宰金珉．高层建筑桩基础沉降分析与优化设计方法研究[J]，南京建筑工程学院学报，2000，第1期：1～9
[30] 陈仁朋，凌道盛，陈云敏．群桩基础沉降计算中的几个问题[J]，土木工程学报，2003，第36卷第10期：89～94
[31] 魏东．桩筏基础沉降研究及应用：[硕士学位论文][D]．杭州：浙江大学，2003
[32] 熊大阅，施云华．上海地区双层地基对超高层建筑桩基沉降影响与探析[J]，岩土工程技术，1996，第4期：3～8
[33] 吴京生．高层建筑桩基沉降计算的探讨[J]，武汉城市建设学院学报，1998，第15卷第2期：52～56
[34] 董炳炎，甘德福．论上海软土地基⑦⑧⑨层压缩模量Es的定值问题[J]，水文地质工程地质，2003，第4期：59～62
[35] 张建勋．桩基沉降的非线性性态研究：[博士学位论文][D]．上海：同济大学，2001
[36] 宰金珉，宰金璋．高层建筑基础分析与设计[M]．北京：中国建筑工业出版社，1993
[37] 林天健等．桩基础设计指南[M]．北京：中国建筑工业出版社，1999
[38] 袁建新，朱国甫．桩的大应变法波动分析优化算法程序[J]，岩土工程学报，1990，(6)
[39] 张学言，闫澍旺．岩土塑性力学基础[M]．天津：天津大学出版社，2004
[40] Goodman R F, Tailor R L, Brekke T L. A model for the mechanics ofjointed rock[J]. Journ Soil Mech & Found Div ASCE, 1968,94(3):637～660
[41] C.S.Desai, Zaman M M. Thin layer element for interface and joints[J]. Int Journ for Num & Analy Mech in Geomech, 1984,8(1): 19～43
[42] C.S.Desai et al. Modeling for cyclic normal and shear behavior of interface[J]. ASCE,1958,vol.34, No.7:398～1217
[43] C.S.Desai et al. Cyclic testing and modeling of interface[J]. ASCE,1985,vol.31,No.6:793～815
[44] E.D.Drumm and C.S.Desai, Determination of parameters for a model of the cyclic behavior of interfaces[J]. Earth.Eng.&S.Dyn., 1986,vol. 14: 1～18
[45] 殷宗泽，朱泓，许国华．土与结构材料接触面的变形及其数学模拟[J]，岩土工程学报，1994，16(3)：14～22
[46] 姚纬明，李同春，任旭华．带软弱结构面岩体的弹塑性有限元分析[J]，河海大学学报，1999，27(3)：34～38
[47] 嘉木工作室．ANSYS 5.7有限实例分析教程[M]．北京：机械工业出版社，2002
[48] 龚曙光，谢桂兰．ANSYS操作命令与参数化编程[M]．北京：机械工业出版社，2005
[49] 美国ANSYS 公司北京办事处．ANSYS非线性分析指南(上、下)．1998
[50] Benson, David J. and Hallquist, John O. A Single Surface Contact Algorithm for the Post-Buckling Analysis of Shell Structures. Computer Methods in Applied Mechanics and Engineering, Vol. 78, No.2. 1990
[51] Cescotto, S. and Zhu, Y.Y. Large Strain Dynamic Analysis Using Solid and Contact Finite Elements Based on a Mixed Formulation—Application to Metalforming. Journal of Metals Processing Technology, Vol. 45,1994: 657～663.
[52] Cescotto, S. and Charilier, R. Frictional Contact Finite Elements Based on Mixed Variational Principles. International Journal for Numerical Method in Engineering, Vol. 36, 1992: 1681～1701.
[53] Giannakopoulos. The Return Mapping Method for the Integration of Friction Constitutive Relations. Computers and Structures, Vol. 32, 1989: 157-167.
[54] Laursen, T.A. and Simo, J.C. Algorithmic Symmetrization of Coulomb Frictional Problems Using Augmented Lagrangians. Computers Methods in Applied Mechanics and Engineering, Vol. 108, No. 1 & 2, 1993: 133-146.
[55] 张季容、朱向荣．简明建筑基础计算与设计手册[M]．北京：中国建筑工业出版社，1997
[56] 林柏、黄建鑫、郑雪芳等．实体深基础法计算桩基沉降经验系数之探讨[J]，建筑施工，2004，第26卷第2期：114～116
[57] 段文峰，廖雄华，金菊顺等．桩—土界面的数值模拟与单桩Q～S曲线的数值分析[J]，哈尔滨建筑大学学报，2001，第34卷第5期：34～38
[58] 薛定宇，陈阳泉．高等应用数学问题的Matlab求解[M]．北京：清华大学出版社，2004
[59] 王沫然．MATLAB6.0与科学计算[M]．北京：电子工业出版社，2001
[60] 郝文化等．ANSYS土木工程应用实例[M]．北京：中国水利水电出版社，2005
[61] 王勖成，邵敏．有限单元法基本原理和数值方法(第二版)[M]．北京：清华大学出版社，1995
[62] 赵均海，王敏强，魏雪英．高等有限元[M]．武汉：武汉理工大学出版社，2004
[63] 曾祥勇．竖向荷载作用下浅基础的层状地基强度和变形有限元分析：[硕士学位论文][D]．重庆：重庆大学，2001
[64] 李婉．南宁地区软岩嵌岩桩的ANSYS计算机模拟：[硕士学位论文][D]．南宁：广西大学，2001
[65] 李爱民．三维弹性接触问题的数值模拟技术：[硕士学位论文][D]．西安：西北工业大学，2004
[2] 廖向东，桩基沉降理论的应用研究：[硕士学位论文][D]．上海：同济大学，2001
[3] 陈祥福，沉降计算理论及工程实例[M]．北京：科学出版社，2005
[4] [澳]H．G Poulos；E．H．Davis著，孙幼兰译．岩土力学弹性解[M]．徐州：中国矿业大学出版社，1990
[5] GB 50007-2002 建筑地基基础设计规范[S]．北京：中国建筑工业出版社，2002
[6] JGJ 94-94 建筑桩基技术规范[S]．北京：中国建筑工业出版社，1994
[7] 徐志英．以明德林公式为根据的地基中垂直应力的计算公式[J]，土木工程学报，1957，(4)：485～497
[8] 袁聚云，赵锡宏．竖向均布荷载作用在地基内部时的土中应力公式[J]，上海力学，1995，第16卷第3期：213～222
[9] 王士杰，张梅，张吉占．工程荷载下的明氏应力实用理论解(I)[J]，四川建筑科学研究，2001，第27卷第2期：36～37
[10] 王士杰，张梅，张吉占．工程荷载下的明氏应力实用理论解(Ⅱ)[J]，四川建筑科学研究，2001，第27卷第4期：60～61
[11] 张子明。求解地基位移的影响函数法[J]，河海大学学报，1990，18(2)：25～31
[12] 陈甦．均质地基中的桩基础沉降计算[J]，岩土力学，1996，第17卷第1期：61～70
[13] K．L．Johnson蓍，徐秉业等译．接触力学[M]．北京：高等教育出版社，1992
[14] 陈希哲．土力学地基基础(第四版)[M]．北京：清华大学出版社，2004
[15] 王铁行，赵树德．计算地基沉降分层总和法的缺陷的分析与改进[J]，西安建筑科技大学学报，1996，第28卷第2期：179～182．
[16] Joseph E. Bowles. Foundations analysis and design (the third edition) [M]. NEW YORK: McGraw-Hill Book Company. 1982.
[17] Harry G. Poulos, Edward H. Davis. Elastic solutions for soil and rock mechanics [M]. NEW YORK: John Wiley & Sons, Inc. 1974.
[18] 应永法，泮威风．圆形板基础下地基强度和变形非线性性状研究[J]，岩土力学，2002，第23卷第2期：250～257
[19] 黄绍铭，裴捷，贾宗元等．软土中桩基沉降估算[J]，第四届全国土力学及基础工程学术会议选集，1986，237～243
[20] 杨敏，王树娟，王伯钧等．使用Geddes应力系数公式求解单桩沉降[J]，同济大学学报，1997，第25卷第4期：379～385
[21] 肖拥军，刘新华，杨荣丰．砂土地基桩基础承载性状试验研究[J]，湖南科技大学学报(自然科学版)，2004，第19卷第2期：60～63
[22] 高强，杨龙才，王炳龙．软粘土钻孔灌注桩承载特性的试验研究[J]，岩土工程技术，2004，第18卷第1期：35～39
[23] 钟闻华，石名磊，刘松玉．超长桩荷载传递性状研究[J]，岩土力学，2005，第26卷第2期：316～318
[24] 王东红，谢星，张炜等．黄土地区超长钻孔灌注桩荷载传递性状试验研究[J]，工程地质学报，1004-9665／2005／13(01)-0117-07
[25] 钟龙辉．西安地区高层建筑桩基沉降计算[J]，勘察科学技术，2001，第3期：9～12
[26] 单玉川、郑学斌．论桩基的沉降量[J]，建筑结构，2001，第31卷第10期：26～Z7
[27] 李云翼、陈国栋、甘德福等．桩基最终沉降量的Gan-Chen模式计算方法[J]，工程勘察，2002，第6期：4～6
[28] 陈仁朋．软弱地基中桩筏基础工作性状和分析设计方法研究：[博士学位论文][D]．杭州：浙江大学，2001
[29] 陈国兴，金永彬，宰金珉．高层建筑桩基础沉降分析与优化设计方法研究[J]，南京建筑工程学院学报，2000，第1期：1～9
[30] 陈仁朋，凌道盛，陈云敏．群桩基础沉降计算中的几个问题[J]，土木工程学报，2003，第36卷第10期：89～94
[31] 魏东．桩筏基础沉降研究及应用：[硕士学位论文][D]．杭州：浙江大学，2003
[32] 熊大阅，施云华．上海地区双层地基对超高层建筑桩基沉降影响与探析[J]，岩土工程技术，1996，第4期：3～8
[33] 吴京生．高层建筑桩基沉降计算的探讨[J]，武汉城市建设学院学报，1998，第15卷第2期：52～56
[34] 董炳炎，甘德福．论上海软土地基⑦⑧⑨层压缩模量Es的定值问题[J]，水文地质工程地质，2003，第4期：59～62
[35] 张建勋．桩基沉降的非线性性态研究：[博士学位论文][D]．上海：同济大学，2001
[36] 宰金珉，宰金璋．高层建筑基础分析与设计[M]．北京：中国建筑工业出版社，1993
[37] 林天健等．桩基础设计指南[M]．北京：中国建筑工业出版社，1999
[38] 袁建新，朱国甫．桩的大应变法波动分析优化算法程序[J]，岩土工程学报，1990，(6)
[39] 张学言，闫澍旺．岩土塑性力学基础[M]．天津：天津大学出版社，2004
[40] Goodman R F, Tailor R L, Brekke T L. A model for the mechanics ofjointed rock[J]. Journ Soil Mech & Found Div ASCE, 1968,94(3):637～660
[41] C.S.Desai, Zaman M M. Thin layer element for interface and joints[J]. Int Journ for Num & Analy Mech in Geomech, 1984,8(1): 19～43
[42] C.S.Desai et al. Modeling for cyclic normal and shear behavior of interface[J]. ASCE,1958,vol.34, No.7:398～1217
[43] C.S.Desai et al. Cyclic testing and modeling of interface[J]. ASCE,1985,vol.31,No.6:793～815
[44] E.D.Drumm and C.S.Desai, Determination of parameters for a model of the cyclic behavior of interfaces[J]. Earth.Eng.&S.Dyn., 1986,vol. 14: 1～18
[45] 殷宗泽，朱泓，许国华．土与结构材料接触面的变形及其数学模拟[J]，岩土工程学报，1994，16(3)：14～22
[46] 姚纬明，李同春，任旭华．带软弱结构面岩体的弹塑性有限元分析[J]，河海大学学报，1999，27(3)：34～38
[47] 嘉木工作室．ANSYS 5.7有限实例分析教程[M]．北京：机械工业出版社，2002
[48] 龚曙光，谢桂兰．ANSYS操作命令与参数化编程[M]．北京：机械工业出版社，2005
[49] 美国ANSYS 公司北京办事处．ANSYS非线性分析指南(上、下)．1998
[50] Benson, David J. and Hallquist, John O. A Single Surface Contact Algorithm for the Post-Buckling Analysis of Shell Structures. Computer Methods in Applied Mechanics and Engineering, Vol. 78, No.2. 1990
[51] Cescotto, S. and Zhu, Y.Y. Large Strain Dynamic Analysis Using Solid and Contact Finite Elements Based on a Mixed Formulation—Application to Metalforming. Journal of Metals Processing Technology, Vol. 45,1994: 657～663.
[52] Cescotto, S. and Charilier, R. Frictional Contact Finite Elements Based on Mixed Variational Principles. International Journal for Numerical Method in Engineering, Vol. 36, 1992: 1681～1701.
[53] Giannakopoulos. The Return Mapping Method for the Integration of Friction Constitutive Relations. Computers and Structures, Vol. 32, 1989: 157-167.
[54] Laursen, T.A. and Simo, J.C. Algorithmic Symmetrization of Coulomb Frictional Problems Using Augmented Lagrangians. Computers Methods in Applied Mechanics and Engineering, Vol. 108, No. 1 & 2, 1993: 133-146.
[55] 张季容、朱向荣．简明建筑基础计算与设计手册[M]．北京：中国建筑工业出版社，1997
[56] 林柏、黄建鑫、郑雪芳等．实体深基础法计算桩基沉降经验系数之探讨[J]，建筑施工，2004，第26卷第2期：114～116
[57] 段文峰，廖雄华，金菊顺等．桩—土界面的数值模拟与单桩Q～S曲线的数值分析[J]，哈尔滨建筑大学学报，2001，第34卷第5期：34～38
[58] 薛定宇，陈阳泉．高等应用数学问题的Matlab求解[M]．北京：清华大学出版社，2004
[59] 王沫然．MATLAB6.0与科学计算[M]．北京：电子工业出版社，2001
[60] 郝文化等．ANSYS土木工程应用实例[M]．北京：中国水利水电出版社，2005
[61] 王勖成，邵敏．有限单元法基本原理和数值方法(第二版)[M]．北京：清华大学出版社，1995
[62] 赵均海，王敏强，魏雪英．高等有限元[M]．武汉：武汉理工大学出版社，2004
[63] 曾祥勇．竖向荷载作用下浅基础的层状地基强度和变形有限元分析：[硕士学位论文][D]．重庆：重庆大学，2001
[64] 李婉．南宁地区软岩嵌岩桩的ANSYS计算机模拟：[硕士学位论文][D]．南宁：广西大学，2001
[65] 李爱民．三维弹性接触问题的数值模拟技术：[硕士学位论文][D]．西安：西北工业大学，2004