›› 2012, Vol. 33 ›› Issue (2): 590-596.

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

针对砂土应变软化强非线性问题的动态松弛有限元法研究

彭芳乐1, 2,李福林1, 3,白晓宇1, 2   

  1. 1.同济大学 地下建筑与工程系,上海 200092;2.同济大学 岩土及地下工程教育部重点实验室,上海 200092; 3. 中国矿业大学 力学与建筑工程学院,江苏 徐州 221116
  • 收稿日期:2010-07-16 出版日期:2012-02-10 发布日期:2012-02-14
  • 作者简介:彭芳乐,男,1965年生,博士,教授,博士生导师,主要从事岩土力学、地下空间与工程等方面的研究工作。
  • 基金资助:

    国家自然科学基金(No. 40972176);教育部长江学者和创新团队发展计划资助(No. IRT1029);上海市重点学科建设项目(No. B308)。

A dynamic relaxation - finite element method for strong nonlinearity caused by post-peak strain softening of sands

PENG Fang-le1, 2, LI Fu-lin1, 3, BAI Xiao-yu1, 2   

  1. 1. Department of Geotechnical Engineering, Tongji University, Shanghai 200092, China; 2. Key Laboratory of Geotechnical and Underground Engineering of Ministry of Education, Tongji University, Shanghai 200092, China; 3. School of Mechanics and Civil Engineering, China University of Mining & Technology, Xuzhou, Jiangsu 221116, China
  • Received:2010-07-16 Online:2012-02-10 Published:2012-02-14

摘要: 对于密实砂土峰值后应变软化特性引起的材料强非线性问题,利用常规有限元所采用的隐式算法很难求解,而动态松弛法在求解这类非线性问题方面却具有独特的优势。针对砂土应变软化强非线性问题,将动态松弛法与有限单元法相结合,提出了一种新的动态松弛有限元法。该方法根据动态松弛法的显式特性,由中心差分法导出了动态松弛有限元法的基本控制方程,并实现了对应力-应变空间中整个平衡路径的追踪。将动态松弛法嵌入到非线性有限元程序中,即可对砂土材料应变软化引起的强非线性问题进行有限元数值计算。该有限元程序中,应力更新采用回归映射算法。最后通过对砂土平面应变压缩试验进行有限元模拟,验证了动态松弛有限元法在求解材料强非线性问题方面的优越性。

关键词: 砂土, 应变软化, 有限单元法, 强非线性, 动态松弛, 显式算法

Abstract: In the presence of the strong material nonlinearity caused by post-peak strain softening of dense sandy soils, the solution of traditional finite element method (FEM) with implicit algorithm often becomes intractable. The dynamic relaxation (DR) method has a great reputation in solving highly nonlinear equations. A new dynamic relaxation-finite element method (DR-FEM) for strong nonlinearity caused by post-peak strain softening of sands is proposed, which takes into account the advantages of DR and FEM together. According to the explicit nature of DR and the central difference technique, the general governing equations of DR-FEM has been derived. The way of tracing the whole equilibrium curve in the stress-strain space has also been presented; and then the DR method is implemented into a general nonlinear finite element codes. The return mapping algorithm is used for stress updating, which is a first-order approximated Euler backward integration. Therefore, the strong material nonlinearity caused by post-peak strain softening of sandy soils can be simulated by the proposed DR-FEM combined with the corresponding material model. The DR-FEM is validated by simulating the result of physical plane strain compression test performed on sands. It is shown that the DR method has a superiority to solve the material nonlinearity.

Key words: sands, strain softening, finite element method, strong nonlinearity, dynamic relaxation, explicit algorithm

中图分类号: 

  • O 242.21
[1] 覃玉兰, 邹新军, 曹雄. 均质砂土中水平简谐荷载与扭矩联合 受荷单桩内力、位移分析[J]. 岩土力学, 2020, 41(1): 147-156.
[2] 王伟, 陈国庆, 郑水全, 张广泽, 王栋, . 考虑张拉-剪切渐进破坏的边坡矢量和法研究[J]. 岩土力学, 2019, 40(S1): 468-476.
[3] 孔亮, 刘文卓, 袁庆盟, 董彤, . 常剪应力路径下含气砂土的三轴试验[J]. 岩土力学, 2019, 40(9): 3319-3326.
[4] 金俊超, 佘成学, 尚朋阳. 基于应变软化指标的岩石非线性蠕变模型[J]. 岩土力学, 2019, 40(6): 2239-2246.
[5] 孙逸飞, 陈 成, . 无状态变量的状态依赖剪胀方程及其本构模型[J]. 岩土力学, 2019, 40(5): 1813-1822.
[6] 徐 鹏, 蒋关鲁, 雷 涛, 刘 琪, 王智猛, 刘 勇, . 考虑填土强度的加筋土挡墙动位移计算[J]. 岩土力学, 2019, 40(5): 1841-1846.
[7] 王凤云, 钱德玲, . 基于统一强度理论深埋圆形隧道围岩的剪胀分析[J]. 岩土力学, 2019, 40(5): 1966-1976.
[8] 魏 星, 张 昭, 王 刚, 张建民, . 饱和砂土液化后大变形机制的离散元细观分析[J]. 岩土力学, 2019, 40(4): 1596-1602.
[9] 庄海洋, 付继赛, 陈 苏, 陈国兴, 王雪剑, . 微倾斜场地中地铁地下结构周围地基液化与变形特性振动台模型试验研究[J]. 岩土力学, 2019, 40(4): 1263-1272.
[10] 王 腾, 吴 瑞. 黏土中海底管线竖向贯入阻力研究[J]. 岩土力学, 2019, 40(3): 871-878.
[11] 张 勋, 黄茂松, 胡志平, . 砂土中单桩水平循环累积变形特性模型试验[J]. 岩土力学, 2019, 40(3): 933-941.
[12] 陆 勇, 周国庆, 杨冬英, 宋家庆, . 砂土剪胀软化、剪缩硬化统一本构的显式计算[J]. 岩土力学, 2019, 40(3): 978-986.
[13] 张成功, 尹振宇, 吴则祥, 金银富, . 颗粒形状对粒状材料圆柱塌落影响的 三维离散元模拟 [J]. 岩土力学, 2019, 40(3): 1197-1203.
[14] 董建勋, 刘海笑, 李 洲. 适用于砂土循环加载分析的边界面塑性模型[J]. 岩土力学, 2019, 40(2): 684-692.
[15] 周翠英, 赵珊珊, 杨旭, 刘镇, . 生态酯类材料砂土改良及工程护坡应用[J]. 岩土力学, 2019, 40(12): 4828-4837.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 梁桂兰,徐卫亚,谈小龙. 基于熵权的可拓理论在岩体质量评价中的应用[J]. , 2010, 31(2): 535 -540 .
[2] 李荣涛. 一种高温下混凝土化学塑性-损伤耦合本构模型[J]. , 2010, 31(5): 1585 -1591 .
[3] 马文涛. 基于灰色最小二乘支持向量机的边坡位移预测[J]. , 2010, 31(5): 1670 -1674 .
[4] 于琳琳,徐学燕,邱明国,闫自利,李鹏飞. 冻融作用对饱和粉质黏土抗剪性能的影响[J]. , 2010, 31(8): 2448 -2452 .
[5] 王 伟,刘必灯,周正华,王玉石,赵纪生. 刚度和阻尼频率相关的等效线性化方法[J]. , 2010, 31(12): 3928 -3933 .
[6] 王海波,徐 明,宋二祥. 基于硬化土模型的小应变本构模型研究[J]. , 2011, 32(1): 39 -43 .
[7] 曹光栩,宋二祥,徐 明. 山区机场高填方地基工后沉降变形简化算法[J]. , 2011, 32(S1): 1 -5 .
[8] 刘华丽 ,朱大勇 ,钱七虎 ,李宏伟. 边坡三维端部效应分析[J]. , 2011, 32(6): 1905 -1909 .
[9] 刘年平 ,王宏图 ,袁志刚 ,刘竟成. 砂土液化预测的Fisher判别模型及应用[J]. , 2012, 33(2): 554 -557 .
[10] 王卫东 ,李永辉 ,吴江斌 . 超长灌注桩桩-土界面剪切模型及其有限元模拟[J]. , 2012, 33(12): 3818 -3824 .