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

中国力学学会学术大会2009 （CCTAM2009 ） 基于ABAQUS 的颗粒材料亚塑性模型的数 值实现及探讨1 楚锡华, 徐远杰,张明龙,余村 （武汉大学土木建筑工程学院，武汉，430072 ） 摘要：颗粒材料的亚塑性模型以Jaumann 应力率张量及变形率张量描述本构关系，本文基于Cauchy 应力的 Jaumann 速率及变形率给出了有限元的切线刚度矩阵，由此指出在严格意义上基于ABAQUS 二次开发的亚 塑性模型数值实现必需借助UEL 接口。为了简化程序开发，本文建议了与有限元相结合的两种近似切线模 量矩阵，即基于矩阵广义逆一致性切向模量矩阵与近似切向模量矩阵，从而形成了通过UMAT 接口实现亚 塑性模型的数值方案。由此降低了程序开发的难度，同时也可借助ABAQUS 的后处理功能，提高了工作效 率。数值算例表明了所开发程序的正确性以及所建议方案的可行性。 关键词：颗粒材料; 亚塑性模型; 切向刚度阵; 切向模量矩阵，ABAQUS 中图分类号：O347;TU 443 文献标识码：A 引 言 颗粒材料是有大量离散固体颗粒组成。其宏观力学行为与自身的微观结构（包括颗粒的 形状及排列）及颗粒的物理性质密切相关，且具有强烈的非线性。考虑到颗粒材料的离散特 性，基于离散颗粒模型的颗粒材料数值模拟得到了快速发展[1-7]，但庞大的计算量限制了其 在实际工程中的应用[8]，为此发展研究基于连续模型的颗粒材料本-构理论具有非常重要的 意义。基于连续模型为更精确的模拟颗粒材料的力学行为，通常需要在经典模型中引入能够 1如Cosserat模型为代表的微极理论[9]，或者引入大量唯象参 表征其微观结构的内尺度变量， 数，如扰动状态概念模型[10]。上述模型在考虑不可恢复的塑性变形时需引入屈服函数以及 相应的塑性势函数，而屈服函数、塑性势函数、以及硬化函数中参数的可变性（有些参数实 验测定存在极大的困难）在某种程度上已成为这些模型的缺点[11]。基于速率型张量函数建 立的亚塑性模型则属于一套新的本-构框架，它扬弃了经典弹塑性理论中应变的弹塑性分解， 因此不需引入屈服函数及塑性势函数的概念，但在该理论中经典弹塑性的基本概念都能找到 相应解释，此外该理论具有严密的数学证明及简洁表达式，且其模型参数具有明确的测定程 序[12] ，因此该模型自提出已在颗粒材料力学行为模拟的各个方面得到了一定的应用 [13-17]。有关亚塑性模型的详细介绍可参看[17]。本文研究亚塑性模型的数值实现及其在 ABAQUS二次开发时的相应处理技术，重点为针对亚塑性模型的率型本构方程如何近似获 得一致性切向模量矩阵。对复杂本构关系，切向模量矩阵的理论推导存在一定的困难，甚至 难以获得理论解[18]，因此如何通过数值的方式来获取切向模量矩阵对数值计算非常重要。 Fellin[18] 以亚塑性模型为例结合ABAQUS 的UMAT接口程序的信息传递特点给出了一种高 效获取近似切向模量矩阵的数值方案，但其存在一定的局限性（具体见下文分析）。本文从 另一个角度给出了获取近似切向模量矩阵的方法，即首先从切向刚度阵的推导出发，从理论 角度给出与有限元相结合的切向模量矩阵，具体工作为针对Gudehus-Bauer模型给出了基于 Cauchy应力的Jaumann速率及变形率张量给出了有限元的切线刚度矩阵。由于切向刚度矩阵 的复杂性（具体可参看下文），指出了在ABAQUS 的二次开发时不能直接应用UMAT接口， 而必须应用UEL接口开发相应的单元。进一步考虑到UEL相对UMAT 的复杂程度，建议了与 有限元相结合的两种近似切线模量矩阵，即基于矩阵广义逆一致性切向模量矩阵与近似切向 1 基金项目：国家自然科学基金 中国力学学会学术大会2009 （CCTAM2009 ） 模量矩阵，从而形成了通过UMAT接口实现亚塑性模型的数值方案。由此降低了程序开发的 难度，同时也可借助ABAQUS 的后处理功能，提高了工作效率。数值算例表明了所建议方 案的可行性以及所开发程序的正确性。 1 亚塑性本构模型 亚塑性本-构模