1 各向同性
各向同性线弹性材料的弹性矩阵为:

式中拉梅常数的表达式为:

因此在编写各向同性材料的umat时,需要两个材料参数,在这里我们使用杨氏模量E和泊松比v。
2 横观各向同性
横观各向同性线弹性材料的弹性矩阵为:

并有关系式:


可见其弹性矩阵需要5个独立的参数,为下列5个工程常数:

下标a代表轴向,下标t代表横向。
3 正交各向异性
正交各向异性线弹性材料的弹性矩阵为:

并有关系式:

因此对于正交各向异性材料,其弹性矩阵需要9个工程常数来确定:

4 程序
使用Fortran90编写umat程序。由于Abaqus默认的umat子程序为Fortran77,因此为了使用f90程序,使用命令:
abaqus make library=xxx.f90
该命令可以生成相应的后缀为obj的文件,之后使用该文件即可。使用上述方法可以避免使用Fortran77进行umat的编写。
建模生成job的inp文件之后,可以使用命令行运行该job:
abaqus job=xxx user=xxx.obj
5 单元测试
对一个实体单元进行单轴拉伸测试:

5.1 各向同性
设置材料参数为:
E=210000, v=0.3
Abaqus计算的结果为:

利用子程序计算的结果为:

5.2 横观各向同性
设置材料参数为:

Abaqus计算的结果为:

子程序计算的结果为:

5.3 正交各向异性
设置材料参数为:

Abaqus计算的结果为:

子程序计算的结果为:

免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删
武汉格发信息技术有限公司,格发许可优化管理系统可以帮你评估贵公司软件许可的真实需求,再低成本合规性管理软件许可,帮助贵司提高软件投资回报率,为软件采购、使用提供科学决策依据。支持的软件有: CAD,CAE,PDM,PLM,Catia,Ugnx, AutoCAD, Pro/E, Solidworks ,Hyperworks, Protel,CAXA,OpenWorks LandMark,MATLAB,Enovia,Winchill,TeamCenter,MathCAD,Ansys, Abaqus,ls-dyna, Fluent, MSC,Bentley,License,UG,ug,catia,Dassault Systèmes,AutoDesk,Altair,autocad,PTC,SolidWorks,Ansys,Siemens PLM Software,Paradigm,Mathworks,Borland,AVEVA,ESRI,hP,Solibri,Progman,Leica,Cadence,IBM,SIMULIA,Citrix,Sybase,Schlumberger,MSC Products...