用户子程序UMAT设置

ABAQUS用户子程序UMAT的设置在property模块进行。CreatmaterialTGeneralTUserMaterial，输入MechanicalConstants,即为子程序中需要输入的参数，如下图所示。Module:三PrcpertyU&erC'efintdFieldUsedJ^utpirtVariables11MaterialLjstrDenselyDeprvsrRe-guhrizotiNarnecMaterial-1Drsnipiion：Corartants：Mod令EdrtMaterialNamt:Miartenal-1De-Ecnpton:Module:三PrcpertyU&erC'efintdFieldUsedJ^utpirtVariables11MaterialLjstrDenselyDeprvsrRe-guhrizotiNarnecMaterial-1Drsnipiion：Corartants：Mod令EdrtMaterialNamt:Miartenal-1De-Ecnpton:Material口皀haviorsUserMeteriageneralechanicalIhermalEU&erMaterialUsermete^ial0p«:MechsmcalM*dhani»lC&nsUintG112输入于程序参数二IUseunsymmetrkmate-~ialstiffnetsDanaCreatmaterialTGeneralTDepvar，在这里面设置状态变量的个数，与子程序中变量STATEV(NSTATEV)有关，如下图所示。MaterialISehavion?UcAifMdRatdlIDepvar|[lEgularizj'JDn山刖Ma^ialUserDefinedField晒t0：UwxOutputV^ridblacplkit3加ILMaterialISehavion?UcAifMdRatdlIDepvar|[lEgularizj'JDn山刖Ma^ialUserDefinedField晒t0：UwxOutputV^ridblacplkit3加IL豆腐日由1hl^diankalIhermalflecfirkaModule::Pr口parlyNumberofEDlutinn-deptstattvariablMiVariablenuribertontrallingelementdeledon[Abaqus/E^plkitonly)iMsSeirialBehaviors。嘶「 输入状态吻量个数0?而用户子程序的格式可以通过搜索ABAQUSDocumentation，在该文档中可以找到子程序UMAT的书写方法以及相关变量的含义。SUBROUTINEUMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT,DRPLDE,DRPLDT,STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,CELENT,DFGRDO,DFGRD1,NOEL,NPT,LAYER,KSPT,JSTEP,KINC)CINCLUDE'ABA_PARAM.INC'CHARACTER*80CMNAMEDIMENSIONSTRESS(NTENS),STATEV(NSTATV),DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),PR0PS(NPR0PS),C00RDS(3),DR0T(3,3),DFGRD0(3,3),DFGRDl(3,3),JSTEP(4)usercodingtodefineDDSDDE,STRESS,STATEV,SSE,SPD,SCDand,ifnecessary,RPL,DDSDDT,DRPLDE,DRPLDT,PNEWDTRETURNEND必须在子程序UMAT中提供材料本构的雅可比（Jacobian）矩阵（即变量DDSDDE），即应力增量对应变增量的变化率。而状态变量STATEV在不需要的时候可以不必定义。SSE,SPD,SCD是一些能量计算相关的变量，对仿真结果没有影响，在不输出能量时可以不必定义。推荐文章《ABAQUS子程序UMAT里弹塑性本构的实现》,里面有关于子程序UMAT的基本原理以及实例展示。同时在ABAQUSDocumentation里《AbaqusUserSubroutinesReferenceGuide》中UMAT部分也有UMAT子程序中变量含义、使用条件等的详细讲解，且其中给出了一个线粘弹性材料的UMAT子程序，但相关内容均为英文，理解有一定困难。建议两部分内容结合起来理解。在提交计算时，需要选择所编辑的子程序，如下图所示：Name:tension2_2Model:tension2Analysisproduct:Abaqus/StandardDescription:SubmissionGeneralMemoryParallelizationPrecisiorPreprocessorPrintoutZUPrintanechooftheinputdataPrint