饱和砂土局部变形带模拟的有限元数值实现

›› 2009, Vol. 30 ›› Issue (6): 1837-1842.

饱和砂土局部变形带模拟的有限元数值实现

黄林冲^{1, 2}，徐志胜¹，孙瑞义³

1. 中南大学土木建筑学院，长沙 410075；2. 西北大学土木与环境工程系，美国伊利诺伊埃文斯顿 60208； 3. 广梅汕铁路有限责任公司，广州 510600

收稿日期:2008-08-01 出版日期:2009-06-10 发布日期:2011-03-09
作者简介:黄林冲，男，1981年生，博士研究生，主要从事岩土工程与地下结构工程的数值计算与安全性研究。
基金资助:
中南大学优秀博士论文创新选题项目(No. 77206)；铁道部科技司重点课题(No. 2006G007-C)；美国国家科学基金项目（No. CMMI-0726908）。

Finite element numerical implementation for simulation of deformation band in saturated sand

HUANG Lin-chong ^{1, 2}，XU Zhi-sheng¹，SUN Rui-yi³

1. School of Civil and Architectural Engineering, Central South University. Changsha 410075, China; 2. Department of Civil and Environmental Engineering, Northwestern University, Evanston, IL 60208, USA; 3. Guangzhou-Meizhou-Shantou Railway Co., Ltd., Guangzhou 510600, China

Received:2008-08-01 Online:2009-06-10 Published:2011-03-09

摘要/Abstract

摘要：

基于有限变形理论，推导了Newton-Raphson 迭代算法在k+1步增量表达的矩阵形式，实现了饱和砂土变形局部化的有限元数值计算，得到了饱和砂土发生局部化变形的准则。基于Galerkin 方法，得到了位移场和应力场的空间离散化矩阵方程；由土体局部变形带的连续性条件，引入第1切线算子，推导出了砂土等颗粒状媒介发生局部化变形的必要条件。基于此核心算法，编制了有限元计算程序，模拟了饱和砂土在不排水条件下平面压缩过程中剪切带的形成与发展；通过比较分析，研究了有限元网格粗细对于土体局部变形带的影响，结果表明，网格粗细的病态依赖只是微小的，它只与变形条带的宽度有关，对于土体所表现出来的其他力学特性没有影响。

关键词: 砂土, 变形带, 局部化, 有限变形理论, 网格

Abstract:

The matrix form of Newton-Raphson incremental solution at the k+1 iteration is updated based on the finite deformation theory, leading the finite element implementation to the simulation of the deformation band in the fluid-saturated sand, and the criteria for localized deformation under drained and undrained conditions are derived and utilized to detect instabilities. The spatial discretization of the displacements and the pressure fields are furnished by the classical Galerkin method; according to the continuity of tractions across the deformation band, the necessary condition for localization is obtained by introducing the first tangent operator. Based on the algebraic manipulation above, a finite element numerical computational program is coded by FORTRAN; and then numerical simulations on saturated sand under undrained condition is performed to study the onset and development of the shear band. Also, the mesh sensitivity is studied; and the results show that the pathological mesh dependence is mild, which is related to the apparent width of the shear band, while the different meshes display the same overall mechanical behavior.

Key words: sand, deformation band, localization, finite deformation theory, mesh

中图分类号:

TB 115

黄林冲，徐志胜，孙瑞义. 饱和砂土局部变形带模拟的有限元数值实现[J]. , 2009, 30(6): 1837-1842.

HUANG Lin-chong，XU Zhi-sheng，SUN Rui-yi. Finite element numerical implementation for simulation of deformation band in saturated sand[J]. , 2009, 30(6): 1837-1842.

参考文献

相关文章 15

[1]	王刚, 韦林邑, 魏星, 张建民, . 压实黏土三轴压缩变形过程中的渗透性变化规律[J]. 岩土力学, 2020, 41(1): 32-38.
[2]	覃玉兰, 邹新军, 曹雄. 均质砂土中水平简谐荷载与扭矩联合受荷单桩内力、位移分析[J]. 岩土力学, 2020, 41(1): 147-156.
[3]	孔亮, 刘文卓, 袁庆盟, 董彤, . 常剪应力路径下含气砂土的三轴试验[J]. 岩土力学, 2019, 40(9): 3319-3326.
[4]	邓华锋, 支永艳, 段玲玲, 潘登, 李建林. 水−岩作用下砂岩力学特性及微细观结构损伤演化[J]. 岩土力学, 2019, 40(9): 3447-3456.
[5]	孙逸飞, 陈成, . 无状态变量的状态依赖剪胀方程及其本构模型[J]. 岩土力学, 2019, 40(5): 1813-1822.
[6]	刘登学, 张友良, 丁秀丽, 黄书岭, 裴启涛, . 数值流形法中基于适合分析T样条的局部网格加密算法[J]. 岩土力学, 2019, 40(4): 1584-1595.
[7]	魏星, 张昭, 王刚, 张建民, . 饱和砂土液化后大变形机制的离散元细观分析[J]. 岩土力学, 2019, 40(4): 1596-1602.
[8]	庄海洋, 付继赛, 陈苏, 陈国兴, 王雪剑, . 微倾斜场地中地铁地下结构周围地基液化与变形特性振动台模型试验研究[J]. 岩土力学, 2019, 40(4): 1263-1272.
[9]	张勋, 黄茂松, 胡志平, . 砂土中单桩水平循环累积变形特性模型试验[J]. 岩土力学, 2019, 40(3): 933-941.
[10]	陆勇, 周国庆, 杨冬英, 宋家庆, . 砂土剪胀软化、剪缩硬化统一本构的显式计算[J]. 岩土力学, 2019, 40(3): 978-986.
[11]	张成功, 尹振宇, 吴则祥, 金银富, . 颗粒形状对粒状材料圆柱塌落影响的三维离散元模拟 [J]. 岩土力学, 2019, 40(3): 1197-1203.
[12]	董建勋, 刘海笑, 李洲. 适用于砂土循环加载分析的边界面塑性模型[J]. 岩土力学, 2019, 40(2): 684-692.
[13]	周翠英, 赵珊珊, 杨旭, 刘镇, . 生态酯类材料砂土改良及工程护坡应用[J]. 岩土力学, 2019, 40(12): 4828-4837.
[14]	王东伟, 陆武萍, 唐朝生, 赵红崴, 李胜杰, 林銮, 冷挺, . 砂土微观结构样品制备技术及量化方法研究[J]. 岩土力学, 2019, 40(12): 4783-4792.
[15]	沙飞, 李术才, 林春金, 刘人太, 张庆松, 杨磊, 李召峰. 砂土介质注浆渗透扩散试验与加固机制研究[J]. 岩土力学, 2019, 40(11): 4259-4269.

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed