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

第3l卷第5期 华侨大学学报(自然科学版) VOL31 No．5 of Science) sep．2010 2010年9月 JournalHuaqiaoUniversity(Natural 文章编号：1000-5013{2010)05—0570-06 一种将线性粘弹微分型本构方程 应用到ABAQUS的方法 潘晓明1，余俊2，杨钊1，孔娟1 (1．同济大学岩土及地下工程教育部重点实验室，上海200092； 2．中南大学土木建筑学院，湖南长沙410075) 本构方程的松弛剪切模量Prony级数形式，并确定其在ABAQUS中的输入参数，供ABAQUS调用分析．以 算例分析岩土流变性态，将3种模型的应变蠕变、松弛数值解与一维情形下的解析解进行对比，结果吻合较 好．对土体采用Merchant粘弹性流变模型的过江盾构隧道进行粘弹性数值分析，得到管片位移均趋于稳定 值，说明所提出的方法是正确和可靠的，可以应用于岩土工程数值计算． 关键词：微分型本构，Laplace变换；Prony级数；蠕变；松弛；ABAQUS 452 中图分类号：TU 文献标识码：A 对于岩土所呈现的各种不同的流变特性，一般采用3种方法进行描述．(1)通过室内或现场试验和 实测[1喝]，建立经验流变或蠕变公式；(2)建立积分型本构方程求解流变问题[4剞；(3)根据流变介质模 型，建立微分型本构方程，对岩土体流变进行分析．第3种方法直观、形象，且便于有限元编程计算，因 而得到广泛的应用[7引．ABAQUS是目前国际上最为先进的大型通用有限元计算分析软件之一．它所 提供的粘弹本构模型是基于松弛剪切模量Prony级数形式的积分型本构方程[10l，这给传统岩土体线性 粘弹微分型本构方程的应用带来了不便．本文研究一种将线性粘弹性微分型本构方程转化成松弛剪切 模量的Prony级数形式的方法，以供ABAQUS调用分析． 1 线性粘弹微分型本构方程转换成Prony级数形式 岩土流变模型是以反映蠕变、应力松弛等为前提，利用弹性元件(胡克体)、粘性元件(牛顿体)，通过 参量Burger模型，分别如图1所示． (a)Burger模型 (b)Merchant模型 (c)Poyning—Thomson模型 图1流变模型 model Fig．1Rheological 1．1微分型本构方程的Laplace变换‘11删 线性粘弹性微分型本构关系可以表示为 收稿日期：2009-10-11 通信作者： 基金项目：国家自然科学基金资助项目；中国博士后科研基金资助项目(20080440652) 万方数据 第5期 潘晓明，等：一种将线性粘弹微分型本构方程应用到ABAQUS的方法 571 薹抓器=客瓠参， 竹≥仇． (1) 令 Pll gK謇，对式(1)进行Laplace变换，可得 。∑㈣ 九矿dkdr，Q=奎k=o P(s)孑=Q(s)三． (2) 为求松弛弹性模量，对恒定的应变8：EO进行Laplace变抉，则式(2)司转抉为 乳j=蔫． (3) 1．2 Burger模型 Burger模型是由Ma