南水模型在ABAQUS中实现及在工程中应用.docVIP

南水模型在ABAQUS中实现及在工程中应用 　　摘要:ABAQUS作为一款强大的通用有限元分析程序,具有很强的非线性计算功能和前、后处理能力。标准ABAQUS程序中包含一些材料的通用本构模型,但土体具有弹塑性及剪胀性等特性,因而有必要在ABAQUS中开发出适应土体的本构模型,南水模型则可以很好地反映土体剪胀等特性。为拓展ABAQUS在岩土工程领域的计算能力,利用二次开发工具UMAT和数据接口将南水模型连接到标准ABAQUS程序中。计算表明,这样不仅可以充分利用ABAQUS程序方便、快捷的前、后处理和强大的非线性求解平台,而且还可以完成更复杂、更切合实际土体应力应变有限元的数值计算。 　　关键词:有限元;南水模型; ABAQUS;二次开发 　　中图分类号:O241.8文献标识码:A文章编号:1672-1683(2010)01-0052-04 　　 　　Development of NHRI Constitutive Model in ABAQUS and Application in Engineering 　　SI Hai-bao1,HUA Xi-ting2 　　(1.Nanjing Hydraulic Research Institute,Nanjing 210024,China;2.Huainan Normal University ,Huainan 232001,China) 　　Abstract: The software of ABAQUS, which is one of the most powerful nonlinear finite element analysis implemented in the world including many different kinds of material models. However, the release of ABAQUS currently absents NHRI constitutive model, which has been widely used to simulate the stress-strain behavior of soil in geo-engineering. The secondary development tools-UMAT and data interface have been used to add NHRI model to standard ABAQUS procedure. All numerical results exhibit that it will be possible to accomplish some complex stress-strain analysis for soil engineering problems by fully making use of the powerful solving platform. 　　Key words: finite element;NHRI model ;ABAQUS;secondary development 　　 　　0 引言 　　作为一款功能强大的有限元分析软件,ABAQUS占有及其重要的位置[1]。它融结构力学、传热力学、流体力学、声学、电学以及热耦合、流固耦合、热电耦合、声电耦合于一体。可以分析复杂的固体力学、结构力学、土力学系统,并能够驾驭非常庞大复杂的工程和模拟高度非线性问题,因而得到了广泛的应用[2]。 　　由于土体具有非线性、非弹性、塑性体应变和剪胀性、塑性剪应变以及硬化与软化等复杂特性, 而ABAQUS仅包括一些通用本构模型,缺少反映土体这些特征的本构模型。在国内,岩土工程领域广泛采用的一些土体本构,如Duncan-Chang本构模型[3](1980),南京水利科学研究院沈珠江的南水模型[4](1986),香港理工大学Li-Dafalias的状态相关砂土本构模型[5](2000),都很好地反映相关土体这些特征,这使得ABAQUS在实际工程应用中受到了极大的限制,仅徐远杰等(2004)在ABAQUS基础上,完成了Duncan-Chang本构模型接口开发[6]。笔者在ABAQUS开发平台的基础上完成了土的本构模型库的建立,包括Duncan-Chang模型、南水模型、状态相关砂土本构模型。 　　计算结果表明,不仅计算效率和精度满足工程需要,而且极大降低了开发、维护工作量和时间周期,是解决复杂土工数值计算的有效途径。 　　1 南水模型本构关系 　　南水模型是把屈服面看作弹性区域的边界由椭圆函数和幂函数组成[4,7],见图1。 　　模型的屈服面方程为: 　　f??1=p2+r2q2(1)