一种将线性粘弹微分型本构方程应用到ABAQUS的方法_潘晓明.pdfVIP

　第 31卷　第 5期 华 侨 大 学 学 报(自 然 科 学 版) Vo l.31　No.5　 　2010 年 9 月 Journal of H uaqiao Univer sity (Na tur al Science) Sep.2010　 　文章编号:　1000-5013(2010)05-0570-06 一种将线性粘弹微分型本构方程 应用到 ABAQUS的方法 潘晓明1 , 余俊2 , 杨钊1 , 孔娟1 (1.同济大学岩土及地下工程教育部重点实验室 , 上海 200092; 2.中南大学 土木建筑学院 , 湖南 长沙 410075) 摘要:　以 3 参量 Merchant模型 、Poyning-Thom son模型 , 以及 4 参量 Burger模型为例 , 推导 3 种常用微分型 本构方程的松弛剪切模量 P rony 级数形式 ,并确定其在 ABAQUS 中的输入参数 , 供 ABAQUS 调用分析.以 算例分析岩土流变性态 ,将 3 种模型的应变蠕变 、松弛数值解与一维情形下的解析解进行对比 , 结果吻合较 好.对土体采用 Merchant粘弹性流变模型的过江盾构隧道进行粘弹性数值分析 , 得到管片位移均趋于稳定 值 ,说明所提出的方法是正确和可靠的 , 可以应用于岩土工程数值计算. 关键词:　微分型本构;Laplace变换;P rony 级数;蠕变;松弛;ABAQUS 中图分类号:　TU 452 文献标识码:　A 对于岩土所呈现的各种不同的流变特性 ,一般采用 3种方法进行描述.(1)通过室内或现场试验和 实测[ 1-3] ,建立经验流变或蠕变公式;(2)建立积分型本构方程求解流变问题[ 4-6] ;(3)根据流变介质模 型 ,建立微分型本构方程 ,对岩土体流变进行分析.第 3种方法直观 、形象 ,且便于有限元编程计算 ,因 而得到广泛的应用 [ 7-9] .ABAQUS是目前国际上最为先进的大型通用有限元计算分析软件之一.它所 提供的粘弹本构模型是基于松弛剪切模量 Prony 级数形式的积分型本构方程 [ 10] ,这给传统岩土体线性 粘弹微分型本构方程的应用带来了不便.本文研究一种将线性粘弹性微分型本构方程转化成松弛剪切 模量的 Prony 级数形式的方法 ,以供 ABAQUS 调用分析. 1　线性粘弹微分型本构方程转换成 Prony级数形式 岩土流变模型是以反映蠕变 、应力松弛等为前提 ,利用弹性元件(胡克体)、粘性元件(牛顿体),通过 串并联关系形成的各种流变模型.文中讨论的 3参量 Merchant模型 、Poyning-Thomson模型 ,以及 4 参量 Burger模型 ,分别如图 1所示. 　　　　(a)Burger 模型　　　　　　　　　(b)Merchant模型　　　　　　(c)Poyning-Thomson 模型 图 1　流变模型 Fig.1　Rheological model 1.1　微分型本构方程的 Laplace变换[ 11-12] 线性粘弹性微分型本构关系可以表示为 　收稿日期:　2009-10-11 　通信作者:　潘晓明(1979-),男 ,博士研究生 , 主要从事隧道及地下结构方面的研究.E-mail:pxm155138@163.com. 　基金项目:　国家自然科学基金资助项目;中国博士后科研基金资助项目(20080440652) ∑ m k=0 pK dkσ dt k =∑ n k =0 qK dkε dt k , 　　n≥m. (1) 令 P=∑ m k=0 pK dkσ dtk ,Q=∑ n k=0 qK dkε dtk ,对式(1)进行 Laplace变换 ,可得 P(s) σ=Q(s) ε. (2) 为求松弛弹性模量 ,对恒定的应变ε=ε0 进行 Laplace 变换 ,则式(2)可转换为 Y(s)= Q(s) s· P(s). (3) 1.2　Burger模型 Burger 模型是由 Maxw ell体和 Kelvin体串联而成(图 1a),其一维本构关系为 σ+p1 ﹒σ+p2σ?=q1﹒ε+q2ε?. (4) 当以弹性模量表示时 ,有 p1 =ηM EM +ηM EK +ηK EM EMEK , 　　p2 = ηMηK EMEK , 　　　　q1 =ηM , 　　q2 =ηMηK EK . (5) 当以剪切模量表示三维形式时 ,则有 p1 =η′M GM +η′MGK +η′KGM GMGK , 　　p2 =η′Mη′K GMGK , 　　　　q1 =2η′M , 　　q2 =2 η′Mη′K GK , (6) GM = EM 2(1 +μM), 　　GK = EK 2(1+μK), η′M = ηM 2(1+μ′M), 　　η′K =