使用单个元素运行黄的UMAT仿真时出错

Question

我试图在Abaqus中的多维数据集元素中使用黄的UMAT (http://www.columbia.edu/~jk2079/Kysar_Research_Laboratory/Single_Crystal_UMAT.html)，我在下面的教程中构建了这个元素：https://www.youtube.com/watch?v=9GQkNKI09ow. 

每次我试图运行模拟时，Abaqus/Standard都会出现一个错误。消息文件不存在，在监视器的警告选项卡中，我可以读到以下内容：“参数沙漏=增强的区段控件选项与这些元素相关: c3d8r、cax4r、cgax4r、cpeg4r、cpe4r、cps4r、m3d4r、s4r、sc8r及其混合、热和压力元素(如果适用)。它也适用于所有类型的修改的三角形和四面体元素。如果该特性仅应用于这些元素类型，则可以忽略此警告”。

下面的信息也出现在我第一次尝试运行模拟时："1元素被定义为零小时玻璃刚度。您可以使用*沙漏刚度或更改元素类型。元素集ErrElemZeroHourGlassStiffness已经识别了这些元素“。

这个问题似乎与沙漏参数有关，但我还没有找到解决办法。有谁能帮我了解一下是什么问题吗？我的目标是运行这个模拟。

*Heading
** Job name: Job-1 Model name: Job-1
** Generated by: Abaqus/CAE 2019
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=CUBE-1
*Node
      1,           1.,           1.,           1.
      2,           1.,           0.,           1.
      3,           1.,           1.,           0.
      4,           1.,           0.,           0.
      5,           0.,           1.,           1.
      6,           0.,           0.,           1.
      7,           0.,           1.,           0.
      8,           0.,           0.,           0.
*Element, type=C3D8R
1, 5, 6, 8, 7, 1, 2, 4, 3
*Nset, nset=SET-1, generate
 1,  8,  1
*Elset, elset=SET-1
 1,
** Section: Section-1-SET-1
*Solid Section, elset=SET-1, controls=EC-1, material=MATERIAL-1
,
*End Part
**  
**
** ASSEMBLY
**
*Assembly, name=Assembly
**  
*Instance, name=CUBE-1, part=CUBE-1
*End Instance
**  
*Nset, nset=SET-1, instance=CUBE-1, generate
 5,  8,  1
*Elset, elset=SET-1, instance=CUBE-1
 1,
*Nset, nset=SET-2, instance=CUBE-1
 3, 4, 7, 8
*Elset, elset=SET-2, instance=CUBE-1
 1,
*Nset, nset=SET-3, instance=CUBE-1, generate
 2,  8,  2
*Elset, elset=SET-3, instance=CUBE-1
 1,
*Nset, nset=SET-4, instance=CUBE-1, generate
 1,  4,  1
*Elset, elset=SET-4, instance=CUBE-1
 1,
*End Assembly
** 
** ELEMENT CONTROLS
** 
*Section Controls, name=EC-1, hourglass=ENHANCED
1., 1., 1.
** 
** MATERIALS
** 
** COMMENTS FROM *DEPVAR
** =====================
** 
**   ITRATN , ITRMAX , GAMERR ,
**     --   ,   --   ,   --   ,
**  ITRATN:  parameter determining whether iteration method is used to
**    solve increments of stresses and state variables in terms of
**    strain increments
** 
**    ITRATN=0. --- no iteration
**    otherwise --- iteration
** 
**  ITRMAX:  maximum number of iterations
** 
**  GAMERR:  absolute error of shear strains in slip systems
** 
** 
*Material, name=MATERIAL-1
*Depvar
    125,
*User Material, constants=160, unsymm
168400.,121400., 75400.,     0.,     0.,     0.,     0.,     0.
     0.,     0.,     0.,     0.,     0.,     0.,     0.,     0.
     0.,     0.,     0.,     0.,     0.,     0.,     0.,     0.
     1.,     0.,     0.,     0.,     0.,     0.,     0.,     0.
     1.,     1.,     1.,     1.,     1.,     0.,     0.,     0.
     0.,     0.,     0.,     0.,     0.,     0.,     0.,     0.
     0.,     0.,     0.,     0.,     0.,     0.,     0.,     0.
    -1.,     0.,     1.,     0.,     0.,     1.,     0.,     0.
     0.,     1.,     0.,     0.,     1.,     0.,     0.,     0.
    10.,  0.001,     0.,     0.,     0.,     0.,     0.,     0.
     0.,     0.,     0.,     0.,     0.,     0.,     0.,     0.
     0.,     0.,     0.,     0.,     0.,     0.,     0.,     0.
  541.5,  109.5,   60.8,     0.,     0.,     0.,     0.,     0.
     1.,     1.,     0.,     0.,     0.,     0.,     0.,     0.
     0.,     0.,     0.,     0.,     0.,     0.,     0.,     0.
     0.,     0.,     0.,     0.,     0.,     0.,     0.,     0.
     0.,     0.,     0.,     0.,     0.,     0.,     0.,     0.
     0.,     0.,     0.,     0.,     0.,     0.,     0.,     0.
    0.5,     1.,     0.,     0.,     0.,     0.,     0.,     0.
     1.,    10.,  1e-05,     0.,     0.,     0.,     0.,     0.
** ----------------------------------------------------------------
** 
** STEP: Step-1
** 
*Step, name=Step-1, nlgeom=YES
*Static
0.5, 1., 1e-05, 1.
** 
** BOUNDARY CONDITIONS
** 
** Name: Disp-BC-1 Type: Symmetry/Antisymmetry/Encastre
*Boundary
SET-1, XSYMM
** Name: Disp-BC-2 Type: Symmetry/Antisymmetry/Encastre
*Boundary
SET-2, ZSYMM
** Name: Disp-BC-3 Type: Symmetry/Antisymmetry/Encastre
*Boundary
SET-3, YSYMM
** Name: Disp-BC-4 Type: Displacement/Rotation
*Boundary
SET-4, 1, 1, 0.001
** 
** OUTPUT REQUESTS
** 
*Restart, write, frequency=0
** 
** FIELD OUTPUT: F-Output-1
** 
*Output, field, variable=PRESELECT
** 
** HISTORY OUTPUT: H-Output-1
** 
*Output, history, variable=PRESELECT
*End Step

非常感谢你的帮助。

Answer 1

我解决了这个问题，我想和大家分享一下我是如何做到的。沙漏只是一个警告，您可以使用C3D20R元素解决这个问题。此错误是由于脚本中的停止字所致。通过将其替换为调用XIT，UMAT工作(https://abaqus-docs.mit.edu/2017/English/SIMACAESUBRefMap/simasub-c-xit.htm)。这适用于根据Abaqus2019中的步骤与Fortran链接的https://www.researchgate.net/publication/349991987_Linking_ABAQUS_20192020_and_Intel_oneAPI_Base_Toolkit_FORTRAN_Compiler。