岩土力学 ›› 2022, Vol. 43 ›› Issue (S2): 533-540.doi: 10.16285/j.rsm.2021.2089

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

基于罚函数偶应力理论的土体应变局部化研究

王冬勇1,陈曦2,王方宇1,彭丽云1,齐吉琳1   

  1. 1. 北京建筑大学 土木与交通工程学院,北京 102616;2. 北京交通大学 土木建筑工程学院,北京 100044
  • 收稿日期:2021-12-10 修回日期:2022-02-16 出版日期:2022-10-10 发布日期:2022-10-10
  • 通讯作者: 齐吉琳,男,1969年生,博士,研究员,主要从事寒区岩土工程方面的研究。E-mail: jilinqi@bucea.edu.cn E-mail:wangdongyong@bucea.edu.cn
  • 作者简介:王冬勇,男,1990年生,博士,博士后,主要从事计算岩土力学方面的研究。
  • 基金资助:
    国家自然科学基金(No.42002277,No.41972279,No.42172299);国家重点研发计划课题(No.2021YFF0306302)

Analysis of geotechnical strain localization based on penalty-based couple stress theory

WANG Dong-yong1, CHEN Xi2, WANG Fang-yu1, PENG Li-yun1, QI Ji-lin1   

  1. 1. School of Civil and Transportation Engineering, Beijing University of Civil Engineering and Architecture, Beijing 102616, China; 2. School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China
  • Received:2021-12-10 Revised:2022-02-16 Online:2022-10-10 Published:2022-10-10
  • Supported by:
    This work was supported by the National Natural Science Foundation of China (42002277,41972279,42172299) and the National Key R&D Program of China (2021YFF0306302).

摘要: 与 Cosserat 理论相比,偶应力理论在一定程度上可以降低数值框架的复杂度,已逐渐应用于岩土体应变局部化分析中。然而,一般的偶应力有限元法需要满足 C连续性,即单元内部和单元交界面上的应变都需要具有连续性。为了避免开发较为复杂的C型偶应力单元,在 Cosserat 连续体理论框架下,通过借助罚函数方法对 C连续性进行松弛来获得偶应力理论的逼近解,建立了基于罚函数的偶应力有限元方法 PCS-FEM。通过平面应变条件下的弹性圆孔应力集中问题对 PCS-FEM 方法的有效性进行了验证,并应用于土体应变局部化分析中。通过对Ottawa砂的平面应变试验进行数值模拟,发现 PCS-FEM 方法获得的应力−应变曲线及剪切带破坏形态与试验结果基本一致,且能够克服经典连续体理论病态的网格敏感性问题,保证应变局部化问题的适定性;通过对承受偏心荷载作用下的土坡应变局部化经典算例进行分析,发现 PCS-FEM 方法同样可以克服土坡应变软化阶段的网格敏感性问题,展现土体的渐进破坏过程。

关键词: Cosserat理论, 偶应力理论, 罚函数, 内部特征长度, 应变局部化

Abstract: Compared with the Cosserat theory, the complexity of the numerical framework can be reduced in the couple stress theory to some extent, and it has been gradually applied to the geotechnical strain localization analysis. However, the C1 continuity should be satisfied in the straight-forward couple stress finite element method, that is, the continuity of the strain inside the elements and on the interface of the elements should be satisfied. To avoid developing more complicated C1 couple stress elements, within the framework of Cosserat continuum theory, the approximate solution of the couple stress theory can be obtained by relaxing the C1 continuity with the aid of the penalty function method, and the penalty-based couple stress finite element method (named PCS-FEM) is developed. The effectiveness of the PCS-FEM is verified by the stress concentration problem of an elastic circular hole under plane strain condition, and it is further applied to the strain localization analysis of soil mass. Based on the numerical simulation of the plane strain experiment of Ottawa sand, it is found that the stress-strain curve and the failure form of shear band obtained by the PCS-FEM method are basically consistent with the test results, and the ill-conditioned mesh sensitivity problem of the classical continuum theory can be avoided by the PCS-FEM so that the well-posedness of the strain localization problem can be ensured. Based on the strain localization analysis of the soil slope subject to the eccentric load, it is found that the mesh sensitivity problem in the strain-softening stage of the soil slope can also be avoided, and the progressive failure process of the slope can be captured by the PCS-FEM.

Key words: Cosserat theory, couple stress theory, penalty function, internal characteristic length, strain localization

中图分类号: 

  • TU452
[1] 王刚, 韦林邑, 魏星, 张建民, . 压实黏土三轴压缩变形过程中的渗透性变化规律[J]. 岩土力学, 2020, 41(1): 32-38.
[2] 邓华锋, 支永艳, 段玲玲, 潘 登, 李建林. 水−岩作用下砂岩力学特性及微细观结构损伤演化[J]. 岩土力学, 2019, 40(9): 3447-3456.
[3] 王 杰, 贡京伟, 赵泽印. 单轴压类岩石试件应变局部化位置、方向及 预警应用研究[J]. 岩土力学, 2018, 39(S2): 186-194.
[4] 王洪波,张庆松,刘人太,李术才,张乐文,郑 卓,张连震. 基于压水试验的地层渗流场反分析[J]. , 2018, 39(3): 985-992.
[5] 高军程 ,郭 莹 ,贾金青 ,涂兵雄,. 基于数字图像测量系统的饱和细砂渐进破坏特性研究[J]. , 2016, 37(5): 1343-1350.
[6] 余 村 ,楚锡华 ,唐洪祥 ,徐远杰 . 基于Cosserat连续体的颗粒破碎影响研究[J]. , 2013, 34(S1): 67-72.
[7] 陈 龙,楚锡华,张明龙,徐远杰. 基于CLoE与Gudehus-Bauer亚塑性模型的颗粒材料应变局部化对比研究[J]. , 2013, 34(11): 3306-3314.
[8] 郑利涛 ,胡志强 ,唐洪祥 . 强间断分析方法在土工结构物渐进破坏过程中的应用[J]. , 2012, 33(9): 2771-2780.
[9] 姚 池 ,姜清辉 ,叶祖洋 ,周创兵 . 裂隙网络无压渗流分析的初流量法[J]. , 2012, 33(6): 1896-1903.
[10] 马 刚 ,常晓林 ,周 伟 ,周创兵 . 基于Cosserat理论的重力坝深层抗滑稳定分析[J]. , 2012, 33(5): 1505-1512.
[11] 贾善坡 ,伍国军 ,陈卫忠. 基于粒子群算法与混合罚函数法的有限元优化反演模型及应用[J]. , 2011, 32(S2): 598-603.
[12] 王小平,封金财. 基于非局部化方法的边坡稳定性分析[J]. , 2011, 32(S1): 247-252.
[13] 陶 帅,王学滨,潘一山,王 玮. 基于摩尔-库仑模型的非线性本构模型的开发及其在应变局部化中的应用[J]. , 2011, 32(S1): 403-0407.
[14] 司 鹄,谢延明,杨春和. 磨料水射流作用下岩石损伤场的数值模拟[J]. , 2011, 32(3): 935-940.
[15] 吕玺琳,黄茂松,钱建固. 真三维状态下土体强度及变形分叉影响[J]. , 2011, 32(1): 21-26.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 姚仰平,侯 伟. 土的基本力学特性及其弹塑性描述[J]. , 2009, 30(10): 2881 -2902 .
[2] 徐金明,羌培,张鹏飞. 粉质黏土图像的纹理特征分析[J]. , 2009, 30(10): 2903 -2907 .
[3] 向天兵,冯夏庭,陈炳瑞,江 权,张传庆. 三向应力状态下单结构面岩石试样破坏机制与真三轴试验研究[J]. , 2009, 30(10): 2908 -2916 .
[4] 石玉玲,门玉明,彭建兵,黄强兵,刘洪佳. 地裂缝对不同结构形式桥梁桥面的破坏试验研究[J]. , 2009, 30(10): 2917 -2922 .
[5] 夏栋舟,何益斌,刘建华. 土-结构动力相互作用体系阻尼及地震反应分析[J]. , 2009, 30(10): 2923 -2928 .
[6] 徐速超,冯夏庭,陈炳瑞. 矽卡岩单轴循环加卸载试验及声发射特性研究[J]. , 2009, 30(10): 2929 -2934 .
[7] 张力霆,齐清兰,魏静,霍倩,周国斌. 淤填黏土固结过程中孔隙比的变化规律[J]. , 2009, 30(10): 2935 -2939 .
[8] 张其一. 复合加载模式下地基失效机制研究[J]. , 2009, 30(10): 2940 -2944 .
[9] 易 俊,姜永东,鲜学福,罗 云,张 瑜. 声场促进煤层气渗流的应力-温度-渗流压力场的流固动态耦合模型[J]. , 2009, 30(10): 2945 -2949 .
[10] 陶干强,杨仕教,任凤玉. 崩落矿岩散粒体流动性能试验研究[J]. , 2009, 30(10): 2950 -2954 .