基础理论与实验研究

基于上下负荷面的弹黏塑性本构模型及应力积分算法实现

展开
  • 1.同济大学 地下建筑与工程系,上海 200092;2.同济大学 土木工程防灾国家重点实验室,上海 200092; 3.湘潭大学 岩土力学与工程安全湖南省重点实验室,湖南 湘潭 411105;4.华东勘测设计研究院交通与市政工程院,浙江 杭州 310014
王智超,男,1980年生,博士后,副教授,主要从事岩土流变学的研究

收稿日期: 2014-10-13

  网络出版日期: 2018-06-09

基金资助

国家杰出青年科学基金(No. 51025932);国家自然科学基金(No. 51308485);湖南省自然科学基金资助项目(No. 12JJ4006)。

An elasto-viscoplastic constitutive model and its stress integration algorithm based on super-subloading yield surface

Expand
  • 1. Department of Geotechnical Engineering, Tongji University, Shanghai 200092, China; 2. State Key Laboratory of Disaster Reduction in Civil Engineering, Tongji University, Shanghai 200092, China; 3. Hunan Key Laboratory of Geomechanics and Engineering Safety, Xiangtan University, Xiangtan, Hunan 411105, China; 4. Transportation and Municipal Engineering Institute, HydroChina Huadong Engineering Corporation, Hangzhou, Zhejiang 310014, China

Received date: 2014-10-13

  Online published: 2018-06-09

Supported by

This work was supported by the National Science Foundation for distinguished Young Scholars(51025932), the National Natural Science Foundation of China(51308485) and the Natural Science Foundation of Hunan Province of China(12JJ4006).

摘要

为了描述天然软土的时间相依以及结构性特征,提出了一种能考虑土体超固结和结构性的实用弹黏塑性本构模型。它以Asaoka和Hashiguchi的上下负荷面作为某一应变速率下的参考屈服面,按照相对过应力的基本思路,新引入了两个能通过不同应变速率三轴压缩试验测定的率敏性参数 和 ,建立了以当前应力、黏塑性应变以及黏塑性应变速率为状态变量的动屈服准则函数,并给出了基于Newton-Raphson迭代的应力积分算法,且成功地将其嵌入到大型有限元软件ABAQUS中。最终通过数值算例来验证模型的正确性以及应力积分算法的可靠性。结果表明:该模型能同时描述土体的率敏性、蠕变以及结构性特征,模型参数物理意义明确、易懂可测,预测结果与试验数据吻合良好,可用于复杂边值问题的有限元计算。

本文引用格式

王智超 ,蒋明镜 ,陈双林 ,蔡 军, . 基于上下负荷面的弹黏塑性本构模型及应力积分算法实现[J]. 岩土力学, 2016 , 37(2) : 357 -366 . DOI: 10.16285/j.rsm.2016.02.007

Abstract

To describe the time-dependent, over-consolidated and structural characteristics of natural soft clay, a simple elasto-viscoplastic constitutive model including the influence of structural behavior of clay is proposed. In the new model, the Asaoka's superloading and Hashiguchi's subloading yield surface is used as a reference yield surface at a specified strain rate. A relative overstress relation is used to obtain the dynamic loading surface at any strain rate based on the reference yield surface. A dynamic equilibrium condition, using the current stress, viscoplastic strain and viscoplastic strain rate as state variables, is adopted as a convergence criterion. Two material parameters and about the rate sensitivity are added into the new model, which can be determined by triaxial compression tests with different strain rates. The new model has been implemented into ABAQUS by a stress integration algorithm using Newton-Raphson iteration. The validity of the model and the reliability of the stress integration algorithm are verified by numerical calculation. The numerical results show that the model can describe the time-dependent behavior of clay such as rate sensitivity, structural characteristics and creep. The material parameters are easy to understand and measure with clear physical meaning. Furthermore, the prediction results of the model are in good agreement with experimental data. The model can be used in a finite element calculation under the complex boundary value problem.
文章导航

/