土体非线性双曲线本构模型的二次开发.docxVIP

土体非线性双曲线本构模型的二次开发 1 abaqus的发展 abaqus由美国hks开发。它被认为是世界上最具功能的有限软件之一。分析各种固体力学和结构力学系统，尤其是能够处理非常复杂和高非线性问题。近些年来,ABAQUS优秀的分析能力和模拟复杂系统的可靠性使得ABAQUS被各国的工业和研究机构广泛采用。自从1997年进入中国以来,越来越多的国内企业和高校也采用ABAQUS作为产品研发和科学研究的工具,内容涉及机械、材料、土木、电子等专业。 为方便用户开发自己需要或感兴趣的模型,ABAQUS提供了方便灵活的二次开发平台,包括若干用户子程序(user subroutines)以及在编程时可以调用的实用程序(utility routines)。Duncan?Chang模型是一种能够很好地反映土体非线性变形特征的本构模型,其参数具有比较明确的物理意义,且可由常规三轴试验得出,在实际工程中得到了广泛应用。对Duncan?Chang模型进行二次开发,将为ABAQUS在岩土工程领域中的应用提供更大的空间。 2 dunpan:开放模型 康纳(Kondner)指出,粘土和砂土的非线性应力应变关系可以用双曲线方程近似描述,即 Duncan和Chang(1970)改进了这一观点,他们采用Mohr?Coulomb破坏准则,建立了一个土体非线性双曲线本构模型,即著名的Duncan?Chang模型,其表达式为 式中:Ei为初始切线弹性模量,由试验确定;Rf为破坏比,表示破坏时的主应力差(σ1-σ3)f与主应力差渐进值(σ1-σ3)ult的比值;ε1为轴向应变。 Duncan?Chang较常用的两种模型分别是E?B模型和E?v模型,两者的主要区别在于选择了不同的三轴试验体变曲线函数形式。E?v模型假定三轴试验中ε3?ε1关系曲线为双曲线,而E?B模型则假定三轴试验中εv?ε1关系曲线为双曲线。在E?v模型中,切线泊松比vt随围压σ3和应力水平(σ1-σ3)/(σ1-σ3)f而变,但体积模量只与σ3有关,在参数整理上比较简单、方便,并有足够的精确度。E?B模型更能符合增量弹性理论,即体积变化只与σ3有关,剪应力不引起体积变化。基于上述原因,本文对E?B模型进行了二次开发。 Duncan?Chang非线性弹性模型虽然反映了土体变形的主要规律,但同时也存在不足:模型不能反映固结压力增加与降低的差别,也没有反映加荷和卸荷对v的影响;同时,没有反映中主应力对E,v和强度指标的影响,不能反映剪胀性;只考虑硬化,不能反映软化等。所以Duncan?Chang模型适用于荷载较小,即破坏之前的情况。 3 模型建立和求解 ABAQUS作为一个大型通用有限元软件,采用的是弹塑性理论的应力符号约定,但Duncan?Chang模型遵从土力学的应力符号约定,两者有所不同。 据ABAQUS二次开发的规定,开发者需要利用UMAT用户子程序来定义材料模型,主要包括: (1) 文头部分。用户借助文头部分来说明子程序的体名、参数和变量,其一般格式为:SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD, (2) 调用外部参数。对Duncan?ChangE?B模型而言,需要输入的参数共有9个:K,n,Rf,c,《,Pa,Kb,m,Aur。其中K,n,Kb,m为试验数据;Rf为破坏比;Aur为卸荷比。 (3) 计算当前应力水平。首先调用实用程序(utility routines)SPRING读取当前应力分量,然后求得各主应力,并对其进行大小、符号的调整。如果是拉应力,还需要做拉应力的修正。 (4) 计算本构关系中各参数,其中包括加卸载准则的选择。 (5) 求解Jacobian矩阵。 (6) 更新应力分量。 ABAQUS是基于几何非线性、材料非线性和接触非线性的通用有限元分析平台,按照ABAQUS 非线性增量加载技术和平衡迭代求解算法,如果ABAQUS在运行过程中发现UAMT模块的存在,则软件将在每个载荷增量步和迭代步中调用这个模块。需要注意的是,在每个增量步的首次迭代中,每个单元积分点都需要调用UMAT 2次,以便形成正确的初始刚度矩阵。ABAQUS的脚本语言Python接口非常友好,是一种功能强大、具有良好开放性的面向对象程序设计语言,移植性也较强。但基于Fortran语言应用的广泛性,本文采用Fortran语言来编写UMAT模块。 4 模型材料参数对模拟结果的影响 模型二次开发完成后,还需采取正确的方法对程序进行验证。Duncan?Chang模型参数是基于常规三轴压缩试验得出的,所以本构的二次开发是否成功,能否应用于工程计算,可以通过它能否模拟常规三轴压缩试验来确定。 试验参数参考了由冯卫星等人编写的《北京细砂土邓肯?张模型参数试验研究》,土样属于