《ABAQUS_Fortran二次开发》由会员分享，可在线阅读，更多相关《ABAQUS_Fortran二次开发（71页珍藏版）》请在金锄头文库上搜索。
1、目 录摘 要IABSTRACTII1.绪论11.1.课题的研究背景11.2.本文的研究内容和方法22.基于ABAQUS软件的二次开发32.1.ABAQUS介绍32.2.ABAQUS各模块简介32.3.ABAQUS的二次开发平台52.4.ABAQUS的二次开发语言63.用户材料子程序UMAT83.1.UMAT开发环境设置83.2.UMAT注意事项93.3.UMAT接口的原理103.4.UMAT的使用方法124.材料非线性问题144.1.材料的弹塑性本构关系144.2.非线性有限元算法理论174.3.增量理论常刚度法公式推导204.4.增量理论切线刚度法公式推导215.UMAT程序设计和编码255
2、.1.本构关系描述255.2.常刚度法程序设计275.3.常刚度法程序编码295.4.切线刚度法程序设计325.5.切线刚度法程序编码355.6.程序的调试386.程序验证406.1.问题描述406.2.本构关系416.3.ABAQUS自带材料模型计算416.4.常刚度法的UMAT验证436.5.切线刚度法的UMAT验证456.6.两种算法的比较分析477.结论与展望517.1.结论517.2.展望51致 谢53参考文献54附1：ABAQUS自带弹塑性材料验证的INP文件55附2：用于算法验证的INP文件61摘 要 ABAQUS软件功能强大，特别是能够模拟复杂的非线性问题，它包括了多种材料本构
3、关系及失效准则模型，并具有良好的开放性，提供了若干个用户子程序接口，允许用户以代码的形式来扩展主程序的功能。 本文主要研究了ABAQUS用户子程序UMAT的开发方法，采用FORTRAN语言编制了各向同性硬化材料模型的接口程序，研究该类材料的弹塑性本构关系极其实现方法。 本文紧紧围绕UMAT的二次开发技术，首先对其接口原理做了详细介绍，然后针对非线性有限元增量理论中的常刚度法和切线刚度法的算法理论做了深入的剖析，推导出了常刚度法和切线刚度法的算法理论的具体表达式，然后分别编制了两种算法的UMAT程序，最后建立了一个具体的验算模型，通过与ABAQUS自带弹塑性本构关系的计算结果相比较，验证两者的正
4、确性。 本文还对常刚度法和切线刚度法得算法效率做了对比，得出了在非线性程度较高时切线刚度法效率高于常刚度法的结论。关键字： ABAQUS、UMAT、有限元、材料非线性、FORTRAN、切线刚度ABSTRACT ABAQUS software powerful, especially to simulate complex non-linear problem, which includes a wide range of material constitutive model and failure criteria, and has a good open, providing a numb
5、er of user subroutine interface that allows users to code form to expand the functions of the main program. This paper studies the user subroutine UMAT of ABAQUS development methods, the use of FORTRAN language isotropic hardening material model of the interface program, studied the effects of such
6、material is extremely elastic-plastic constitutive relation method. This article UMAT tightly around the secondary development of technology, the first principle of its interface detail, and then for the theory of nonlinear finite element incremental stiffness of the regular tangent stiffness method
7、 and the theory of algorithms to do an in-depth analysis of deduced a regular tangent stiffness and rigidity of the law of the specific expression of algorithm theory, and then the preparation of the two algorithms, respectively, of the UMAT program, and finally the establishment of a specific model
8、 checking, bringing with ABAQUS elasto-plastic constitutive relation of the calculated results compared to verify the correctness of the two. This article also often stiffness and tangent stiffness method was to do a comparison of algorithm efficiency is obtained when a higher degree in the non-line
9、ar tangent stiffness method more efficient than the conclusions of law often stiffness. KEY WORDS：ABAQUS、UMAT、Finite element、Material nonlinearity、FORTRAN、Tangent stiffness1. 绪论1.1. 课题的研究背景 有限单元法基本思想的提出，可以追溯到克劳夫(R.W.Clough)在1943年的工作1，他第一次尝试应用定义在三角形区域上的分片连续函数和最小位能原理相结合，来求解St. Venant扭转问题。1960年克劳夫进一步处理
10、了平面弹性问题，并第一次提出了“有限单元法”的名称，使人们开始认识了有限单元法的功效。 四十多年来，随着电子计算机的广泛应用和发展，有限单元法的理论和应用都得到迅速的，持续不断的发展，其应用己由弹性力学平面问题扩展到空间问题、板壳问题，由静力学问题扩展到稳定问题、动力问题和波动问题。分析的对象从弹性材料扩展到塑性、粘弹性、粘塑性和复合材料等，从固体力学扩展到流体力学、传热学等连续介质力学领域。在工程分析中的作用已从分析和校核扩展到优化设计并和计算机辅助设计。 利用有限元软件解决工程和科学问题，是有限元理论应用于工程设计和科学研究实践的主要形式。由于工程设计的巨大市场需要，有限元软件的发展是很迅
11、速的，目前常用的大型有限元软件常见的有Sap2000，ADINA，MSC/NASTRAN，MSC Marc，ANSYS，ABAQUS等，这些软件的共同特点是具有丰富的单元库和求解器，强大而可靠的分析功能，人们利用这些软件解决了很多工程建设和工业产品设计中遇到的问题，取得了巨大的经济技术效益。 由于工程问题的千差万别，不同的用户有不同的专业背景和发展方向，通用软件不免在具体的专业方面有所欠缺，针对这些不足，大部分的通用软件都提供了二次开发功能，以帮助用户减少重复性的编程工作、提高开发起点、缩短研发周期、降低开发成本，并能简化后期维护工作，给用户带来很多方便。基于通用软件平台进行开发，是目前研究的
12、一个重要发展方向。 ABAQUS也提供了若干用户子程序（User Subroutines）接口，它是一个功能非常强大且适用的分析工具，与命令行的程序格式相比，用户子程序的限制少得多，从而使用更加灵活方便。针对ABAQUS所提供的本构关系模型种类有限，无法满足工程应用需要的问题，用户子程序中的用户材料子程序（User-defined MateriaMechanical Behavior，简称UMAT）接口可以帮助用户定义自己的材料本构模型和算法，这是ABAQUS的独到之处。由于其操作方便，能被灵活地应用于各个领域中，尤其受到用户的青睐。1.2. 本文的研究内容和方法 ABAQUS中用户材料子程序
13、UMAT的开发主要解决两方面的问题：本构模型的建立和积分算法的选择。 本文主要研究非线性材料的UMAT实现方法，并重点研究其迭代算法部分，目前，用户材料子程序UMAT的迭代算法主要是常刚度法，常刚度法的优点在于算法原理较简单，程序编写较方便，缺点是当遇到复杂非线性材料时，其迭代次数较多，收敛速度也较慢，在这个情况下，本文采取的是一种迭代次数较少且收敛速度较快的切线刚度法，具体就是采用FORTRAN语言编制了基于Von-Mises模型的接口程序,并采用切线刚度算法，通过与ABAQUS自带本构关系计算的结果相比较，验证其正确性。 本文的研究工作紧紧围绕UMAT的二次开发技术，首先根据有限元方法推导
14、材料非线性问题算法的公式，然后参考UMAT接口规范设计程序的算法流程，继而编写出该程序，最后建立一个具体的本构和具体的模型做测试，验证程序的正确性，在这一过程中，调试是一个非常重要的过程，占用了大量的时间，在调试程序时采用了将中间变量输出到文本的方式，这样能明确跟进迭代过程，发现算法或程序的缺陷。 本文采用的本构关系是经过归纳和抽象的，也就是说本文的程序并不仅仅是只针对某个具体模型和问题，而是针对所有符合抽象出的各向同性硬化材料，这样做的好处是能保证程序的通用性和复用性，避免以后的重复劳动，当然，这也是符合ABAQUS软件设计UMAT接口的宗旨的。 2. 基于ABAQUS软件的二次开发2.1.
15、 ABAQUS介绍 ABAQUS是一套功能强大的基于有限元法的工程模拟软件2，其解决问题的 范围从相对简单的线性分析到最富有挑战性的非线性模拟问题。ABAQUS具备十分丰富的、可模拟任意实际形状的单元库。并与之对应拥有各种类型的材料模型库，可以模拟大多数典型工程材料的性能，其中包括金属、橡胶、高分子材料、复合材料、钢筋混凝土、可压缩弹性的泡沫材料以及岩石和土这样的地质材料。 作为通用的模拟分析工具，ABAQUS 不仅能解决结构分析中的问题，还能模拟和研究各种领域中的问题，如热传导、质量扩散、电子元器件的热控制（热一电耦合分析）、声学分析、土壤力学分析（渗流应力耦合分析）和压电介质力学分析。 ABAQUS为用户提供了广泛的功能，且使用起来又十分简明。最复杂的问题也可以很容易地建立模型3。例如复杂的多部件问题可以通过对每个部件定义材料模型和几何形状，然后再把它们组装起来而构成。在大部分模拟分析问题中，甚至在高度非线性问题中，用户也只需要提
