›› 2012, Vol. 33 ›› Issue (4): 1240-1246.

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

各向异性弹黏塑性模型在ABAQUS中的研发

钦亚洲1, 2,李 宁1, 2,许建聪1, 2   

  1. 1. 同济大学 岩土及地下工程教育部重点实验室,上海 200092;2. 同济大学 地下建筑与工程系,上海 200092
  • 收稿日期:2011-08-01 出版日期:2012-04-13 发布日期:2012-04-26
  • 作者简介:钦亚洲,男,1973年生,博士研究生,主要从事岩土本构理论及数值模拟方面的研究工作
  • 基金资助:

    国家自然科学基金资助(No. 40872179)。

Development of anisotropic elasto-viscoplastic model in ABAQUS software

QIN Ya-zhou1, 2,LI Ning1, 2,XU Jian-cong1, 2   

  1. 1. Key Laboratory of Geotechnical and Underground of Ministry of Education, Tongji University, Shanghai 200092, China; 2.Department of Geotechnical Engineering, Tongji University, Shanghai 200092, China
  • Received:2011-08-01 Online:2012-04-13 Published:2012-04-26

摘要: 通过将Perzyna过应力理论与临界状态理论相结合,并引入Wheeler旋转硬化法则,提出一个能描述土体初始各向异性及应力诱发各向异性的三维弹黏塑性本构模型。模型考虑流变发生的下限,在三维应力空间,模型存在形状相似的静屈服面及动态加载面。采用缩放形式的幂函数。本构模型数值算法采用回映算法,借助ABAQUS软件UMAT子程序接口实现。并通过对三轴不排水蠕变试验的模拟,确定合适的积分步长。此后,分别对三轴不排水蠕变试验及常应变率三轴不排水剪切试验进行了模拟。模拟中通过设置不同参数值,可将模型退化为各向同性模型,并对这两种模拟结果进行了比较。模拟结果表明:(1) 对于三轴不排水蠕变,在低剪应力水平下,各向同性模型和各向异性模型模拟的结果相差不大,而在高剪应力水平下,各向异性模型模拟结果更接近试验结果;(2) 对于常应变率加载试验的模拟,模型合理反映了土体不排水强度随着加载速率的增大而增大现象。

关键词: 各向异性, 弹黏塑性模型, 回映算法, ABAQUS软件, UMAT子程序

Abstract: This paper presents an improved three-dimensional anisotropic elasto-viscoplastic constitutive model which is based on overstress theory of Perzyna and critical state theory. In this model a rotational hardening law of Wheeler is also adopted to account for initial anisotropy and changes in anisotropy due to stress. The viscoplastic strain will not occur when the stress state is located within the static yield surface. A power-type scaling function is adopted for the viscoplastic strain-rate. The constitutive model is programmed in ABAQUS software by writing UMAT subroutine, which is integrated with Return Mapping Algorithm. Anisotropically consolidated undrained (CU) triaxial creep test for Sackville clay is simulated by proposed model and the suitable integration time step is determined. Then other CU triaxial creep tests and CU triaxial constant strain rate tests are simulated. By switching the anisotropic features off, the proposed model is degraded to the isotropic model. The simulation results show (1) The CU triaxial creep tests with rotational hardening law is more accurate than isotropic model in the high level of the shear stress; (2) The CU triaxial constant strain rate tests rationally reflect that the undrained strength of soil increase with loading rate.

Key words: Anisotropic, Elastic-viscoplastic, return mapping algorithm, ABAQUS software, UMAT subroutine

中图分类号: 

  • TU 443
[1] 徐杰, 周建, 罗凌晖, 余良贵, . 高岭-蒙脱混合黏土渗透各向异性模型研究[J]. 岩土力学, 2020, 41(2): 469-476.
[2] 范运辉, 朱其志, 倪涛, 张坤, 张振南, . 基于弹性张量离散化的脆延转变本构模型研究[J]. 岩土力学, 2019, 40(S1): 181-188.
[3] 彭守建, 岳雨晴, 刘义鑫, 许江, . 不同成因结构面各向异性特征及其剪切力学特性[J]. 岩土力学, 2019, 40(9): 3291-3299.
[4] 尹晓萌, 晏鄂川, 王鲁男, 陈 利, . 各向异性片岩的微观组构信息定量 提取与断面形貌特征分析[J]. 岩土力学, 2019, 40(7): 2617-2627.
[5] 加瑞, 雷华阳, . 有明黏土各向异性固结特性的试验研究[J]. 岩土力学, 2019, 40(6): 2231-2238.
[6] 张坤勇, 臧振君, 李 威, 文德宝, Charkley Frederick Nai, . 土体三维卸荷弹塑性模型及其试验验证[J]. 岩土力学, 2019, 40(4): 1313-1323.
[7] 张玉伟, 翁效林, 宋战平, 谢永利, . 考虑黄土结构性和各向异性的修正剑桥模型[J]. 岩土力学, 2019, 40(3): 1030-1038.
[8] 柯志强, 王环玲, 徐卫亚, 林志南, 吉 华, . 含横向节理的柱状节理岩体力学特性试验研究[J]. 岩土力学, 2019, 40(2): 660-667.
[9] 周建, 蔡露, 罗凌晖, 应宏伟, . 各向异性软土基坑抗隆起稳定极限平衡分析[J]. 岩土力学, 2019, 40(12): 4848-4856.
[10] 田雨, 姚仰平, 路德春, 杜修力, . 基于修正应力法的横观各向同性摩尔-库仑 准则及被动土压力公式[J]. 岩土力学, 2019, 40(10): 3945-3950.
[11] 唐洪祥, 韦文成. 耦合强度各向异性与应变软化的边坡稳定 有限元分析[J]. 岩土力学, 2019, 40(10): 4092-4100.
[12] 周 辉, 程广坦, 朱 勇, 张春生, 卢景景, 张传庆, 章颖辉, . 基于3D雕刻技术的岩体结构面剪切各向异性研究[J]. 岩土力学, 2019, 40(1): 118-126.
[13] 李深圳, 沙 鹏, 伍法权, 伍 劼. 层状结构岩体变形的各向异性特征分析[J]. 岩土力学, 2018, 39(S2): 366-373.
[14] 张坤勇,李 威,Charkley Nai Frederick,陈 恕,. 小主应力方向加载条件下的掺砾黏土真三轴试验[J]. , 2018, 39(9): 3270-3276.
[15] 杨建民,霍王文,. 渗透性水平向各向异性导致椭圆形地面沉降漏斗的一个性质[J]. , 2018, 39(8): 2960-2976.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 徐金明,羌培,张鹏飞. 粉质黏土图像的纹理特征分析[J]. , 2009, 30(10): 2903 -2907 .
[2] 董 诚,郑颖人,陈新颖,唐晓松. 深基坑土钉和预应力锚杆复合支护方式的探讨[J]. , 2009, 30(12): 3793 -3796 .
[3] 梁桂兰,徐卫亚,谈小龙. 基于熵权的可拓理论在岩体质量评价中的应用[J]. , 2010, 31(2): 535 -540 .
[4] 李荣涛. 一种高温下混凝土化学塑性-损伤耦合本构模型[J]. , 2010, 31(5): 1585 -1591 .
[5] 马文涛. 基于灰色最小二乘支持向量机的边坡位移预测[J]. , 2010, 31(5): 1670 -1674 .
[6] 于琳琳,徐学燕,邱明国,闫自利,李鹏飞. 冻融作用对饱和粉质黏土抗剪性能的影响[J]. , 2010, 31(8): 2448 -2452 .
[7] 王协群,张有祥,邹维列,熊海帆. 降雨入渗条件下非饱和路堤变形与边坡的稳定数值模拟[J]. , 2010, 31(11): 3640 -3644 .
[8] 王 伟,刘必灯,周正华,王玉石,赵纪生. 刚度和阻尼频率相关的等效线性化方法[J]. , 2010, 31(12): 3928 -3933 .
[9] 王海波,徐 明,宋二祥. 基于硬化土模型的小应变本构模型研究[J]. , 2011, 32(1): 39 -43 .
[10] 曹光栩,宋二祥,徐 明. 山区机场高填方地基工后沉降变形简化算法[J]. , 2011, 32(S1): 1 -5 .