›› 2014, Vol. 35 ›› Issue (6): 1795-1800.

• 数值分析 • 上一篇    下一篇

基于同伦分析法的非饱和土层固结问题的级数解

李纪伟1, 2,汪华斌1, 2,张 玲1, 2   

  1. 1. 华中科技大学 土木工程与力学学院,武汉 430074;2. 华中科技大学 岩土与地下工程研究所,武汉 430074
  • 收稿日期:2013-04-17 出版日期:2014-06-10 发布日期:2014-06-20
  • 作者简介:李纪伟,男,1984年生,博士研究生,主要从事岩土工程数值模拟、计算力学与应用数学相关研究。
  • 基金资助:

    国家自然科学基金(41372296);国家科技部“十二五”支撑计划项目(No. 2012BAK10B02-2);广东省交通厅科技项目(No. 2011-02-028)

Series solution to one-dimensional consolidation in unsaturated soils based on the homology analysis method

LI Ji-wei1, 2,WANG Hua-bin1, 2,ZHANG Ling1, 2   

  1. 1. College of Civil Engineering and Mechanics, Huazhong University of Science and Technology, Wuhan 430074, China; 2. Institute of Geotechnical and Underground Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
  • Received:2013-04-17 Online:2014-06-10 Published:2014-06-20

摘要: 由于非饱和土的渗透系数是基质吸力的函数,使得控制方程带有强非线性的特征,进而使得控制方程的解析求解变得十分困难。同伦分析法对级数基函数和辅助线性算子的选择具有更大的自由性、灵活性,且收敛性的控制和调节更加容易实现,求解强非线性微分方程时在选择线性算子以及辅助参数上具有明显的优势。因此,针对非饱和土固结方程的非线性特征,对于处于地表浅层的非饱和土层,假设孔隙气压力为大气压力,在Richard经验公式与非饱和土一维固结理论的基础上,推导了非饱和一维固结无量纲控制方程;应用同伦分析法,通过选取适当的初始猜测解与辅助参数,将该非线性方程转换为线性的微分方程组并求解得到固结问题的级数解。此外,以压实高岭土为研究对象,在收集相关试验参数基础之上,将由同伦分析法求得的固结问题的近似解析解与有限差分法数值结果相对比,分析结果验证了解析解的正确性。

关键词: 非饱和土, 一维固结, 级数解, 同伦分析法

Abstract: In a one-dimensional consolidation equation, the permeability coefficient is the function of matric suction in an unsaturated soil. The consolidation equation in one dimension is strongly nonlinear due to the presence of the permeability coefficient. As an analytic solution to nonlinear problems, a homolopy analysis method (HAM) is efficient in the selection of series basis functions and auxiliary linear operators, and easily make the solution convergence. In order to analytically solve the equation, the HAM was introduced in the present study. During the solution, the pore air pressure was assumed as the atmospheric pressure. The equations with two unknown variables were then reduced to dimensionless pore water pressure as only one basic unknown variable subjected to constant pore air pressure. Firstly, a governing equation in a dimensionless form was derived from the basic one-dimensional consolidation theory after the integration with the Richard empirical formula. The method was then used for a mapping technique to transfer the original nonlinear differential equations to a number of linear differential equations. These differential equations are not dependent on any small parameters, which is convenient to control the convergence region. After this transferring, a series solution to the equations was then obtained by the HAM after selection of auxiliary linear operator parameters. Finally, comparisons were carried out between the analytical solutions and the finite difference method in case of compacted kaolin. It can be found that the series solutions indicate that the pore water pressure increases firstly, and then decreases with the depth after the consolidation of the compacted kaolin. The results indicate that the analytical solution in the present study is reasonable.

Key words: unsaturated soil, one-dimensional consolidation, series solution, homology analysis method (HAM)

中图分类号: 

  • TU 431
[1] 程昊, 唐辉明, 吴琼, 雷国平. 一种考虑水力滞回效应的非饱和土弹塑性扩展 剑桥本构模型显式算法有限元实现[J]. 岩土力学, 2020, 41(2): 676-686.
[2] 程涛, 晏克勤, 胡仁杰, 郑俊杰, 张欢, 陈合龙, 江志杰, 刘强, . 非饱和土拟二维平面应变固结问题的解析计算方法[J]. 岩土力学, 2020, 41(2): 453-460.
[3] 邓子千, 陈嘉帅, 王建伟, 刘小文, . 基于SFG模型的统一屈服面本构模型与试验研究[J]. 岩土力学, 2020, 41(2): 527-534.
[4] 李潇旋, 李涛, 彭丽云, . 控制吸力循环荷载下非饱和黏性土 的弹塑性双面模型[J]. 岩土力学, 2020, 41(2): 552-560.
[5] 刘丽, 吴羊, 陈立宏, 刘建坤, . 基于数值模拟的湿润锋前进法测量精度分析[J]. 岩土力学, 2019, 40(S1): 341-349.
[6] 周凤玺, 柳鸿博, . 非饱和土中Rayleigh波的传播特性分析[J]. 岩土力学, 2019, 40(8): 3218-3226.
[7] 詹良通, 胡英涛, 刘小川, 陈捷, 王瀚霖, 朱斌, 陈云敏. 非饱和黄土地基降雨入渗离心模型试验 及多物理量联合监测[J]. 岩土力学, 2019, 40(7): 2478-2486.
[8] 周凤玺, 高国耀, . 非饱和土中热−湿−盐耦合作用的稳态分析[J]. 岩土力学, 2019, 40(6): 2050-2058.
[9] 陶高梁, 吴小康, 甘世朝, 肖衡林, 马 强, 罗晨晨, . 不同初始孔隙比下非饱和黏土渗透性 试验研究及模型预测[J]. 岩土力学, 2019, 40(5): 1761-1770.
[10] 童立红, 王 珏, 郭生根, 朱怀龙, 徐长节, . 变荷载下连续排水边界黏弹性地基 一维固结性状分析[J]. 岩土力学, 2019, 40(5): 1862-1868.
[11] 方瑾瑾, 冯以鑫, 赵伟龙, 王立平, 余永强, . 真三轴条件下原状黄土的非线性本构模型[J]. 岩土力学, 2019, 40(2): 517-528.
[12] 蒙宇涵, 陈征, 冯健雪, 李红坡, 梅国雄, . 初始孔压非均布下一维地基砂垫层优化研究[J]. 岩土力学, 2019, 40(12): 4793-4800.
[13] 杨明辉, 陈贺, 陈可. 基于分形理论的SWCC边界曲线 滞后效应模型研究[J]. 岩土力学, 2019, 40(10): 3805-3812.
[14] 陈正汉, 郭 楠、. 非饱和土与特殊土力学及工程应用研究的新进展[J]. 岩土力学, 2019, 40(1): 1-54.
[15] 段晓梦,曾立峰, . 非饱和土的承载结构与岩土广义结构性[J]. , 2018, 39(9): 3103-3112.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 刘 晓,唐辉明,刘 瑜. 基于集对分析和模糊马尔可夫链的滑坡变形预测新方法研究[J]. , 2009, 30(11): 3399 -3405 .
[2] 胡大伟,周 辉,谢守益,张 凯,邵建富,冯夏庭. 大理岩破坏阶段Biot系数研究[J]. , 2009, 30(12): 3727 -3732 .
[3] 师旭超,韩 阳. 卸荷作用下软黏土回弹吸水试验研究[J]. , 2010, 31(3): 732 -736 .
[4] 朱建明,彭新坡,姚仰平,徐金海. SMP准则在计算煤柱极限强度中的应用[J]. , 2010, 31(9): 2987 -2990 .
[5] 原喜忠,李 宁,赵秀云,杨银涛. 东北多年冻土地区地基承载力对气候变化敏感性分析[J]. , 2010, 31(10): 3265 -3272 .
[6] 白 冰,李小春,石 露,唐礼忠. 弹塑性应力-应变曲线的斜率恒等式及其验证和应用[J]. , 2010, 31(12): 3789 -3792 .
[7] 唐利民. 地基沉降预测模型的正则化算法[J]. , 2010, 31(12): 3945 -3948 .
[8] 李占海,朱万成,冯夏庭,李绍军,周 辉,陈炳瑞. 侧压力系数对马蹄形隧道损伤破坏的影响研究[J]. , 2010, 31(S2): 434 -441 .
[9] 蔡辉腾,危福泉,蔡宗文. 重庆主城区粉质黏土动力特性研究[J]. , 2009, 30(S2): 224 -228 .
[10] 宋 玲 ,刘奉银 ,李 宁 . 旋压入土式静力触探机制研究[J]. , 2011, 32(S1): 787 -0792 .