基于ABAQUS的颗粒材料亚塑性模型的数值实现及探讨  

Numerical implementations and discussion on hypoplastic model for granular materials based on ABAQUS

作　　者：楚锡华[1] 徐远杰[1] 张明龙[1] 余村[1] 

机构地区：[1]武汉大学土木建筑工程学院,武汉430072

出　　处：《岩土力学》2009年第S1期215-219,224,共6页Rock and Soil Mechanics

基　　金：国家自然科学基金(No.10802060)

摘　　要：颗粒材料的亚塑性模型以Jaumann应力率张量及变形率张量描述本构关系,基于Cauchy应力的Jaumann速率及变形率给出了有限元的切线刚度矩阵,由此指出在严格意义上基于ABAQUS二次开发的亚塑性模型数值实现必需借助UEL接口。为了简化程序开发,文中建议了与有限元相结合的两种近似切线模量矩阵,即基于矩阵广义逆一致性切向模量矩阵与近似切向模量矩阵,从而形成了通过UMAT接口实现亚塑性模型的数值方案。由此降低了程序开发的难度,同时也可借助ABAQUS的后处理功能,提高了工作效率。数值算例表明了所开发程序的正确性以及所建议方案的可行性。The constitutive relation of hypoplastic model for granular materials is represented by Jaumann stress rate and deformation rate,the tangential stiffness matrix for finite element simulation of hypoplasticity is given based on Gudehus-Bauer model with Jaumann rate of Cauchy stress and deformation rate,and it is illustrated that the numerical implementations of hypoplastic model based on ABAQUS must resort to UEL interface proper.For simplifying the of programming development,two methods for tangential modulus matrix are suggested combing with finite element method,viz.tangential modulus matrix based on general invert of matrix and approximate tangential modulus matrix,so numerical methods for implementations of hypoplastic model by UMAT interface are formed.By the approach suggested;the complexity of programming is reduced;and the post processing of ABAQUS can be utilized,so the work efficiency is improved.Numerical examples show the validity of methods suggested and fensibility of programming code.

关 键 词：颗粒材料 亚塑性模型 切向刚度阵 切向模量矩阵 ABAQUS 

分 类 号：TU43[建筑科学—岩土工程]