尝试使用单个元素使用 Huang 的 UMAT 进行 运行 模拟时出错
Error when trying to run simulation with Huang's UMAT using a single element
我正在尝试使用Huang的UMAT (http://www.columbia.edu/~jk2079/Kysar_Research_Laboratory/Single_Crystal_UMAT.html) in a cube element in Abaqus that I built following this tutorial: https://www.youtube.com/watch?v=9GQkNKI09ow.
每次我尝试 运行 模拟时,Abaqus/Standard 分析都会因错误而退出。消息文件不存在,在监视器的警告选项卡中,我可以看到以下内容:“部分控件选项上的参数 hourglass = enhanced 与这些元素相关: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
非常感谢您的帮助。
我解决了这个问题,我想分享我是如何解决的。沙漏只是一个警告,您可以使用 C3D20R 元素解决。错误是由于脚本中的 STOP 一词引起的。通过将其替换为 CALL XIT,UMAT 可以工作(https://abaqus-docs.mit.edu/2017/English/SIMACAESUBRefMap/simasub-c-xit.htm). This applies to Abaqus2019 linked with Fortran according to the steps in https://www.researchgate.net/publication/349991987_Linking_ABAQUS_20192020_and_Intel_oneAPI_Base_Toolkit_FORTRAN_Compiler。
我正在尝试使用Huang的UMAT (http://www.columbia.edu/~jk2079/Kysar_Research_Laboratory/Single_Crystal_UMAT.html) in a cube element in Abaqus that I built following this tutorial: https://www.youtube.com/watch?v=9GQkNKI09ow.
每次我尝试 运行 模拟时,Abaqus/Standard 分析都会因错误而退出。消息文件不存在,在监视器的警告选项卡中,我可以看到以下内容:“部分控件选项上的参数 hourglass = enhanced 与这些元素相关: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
非常感谢您的帮助。
我解决了这个问题,我想分享我是如何解决的。沙漏只是一个警告,您可以使用 C3D20R 元素解决。错误是由于脚本中的 STOP 一词引起的。通过将其替换为 CALL XIT,UMAT 可以工作(https://abaqus-docs.mit.edu/2017/English/SIMACAESUBRefMap/simasub-c-xit.htm). This applies to Abaqus2019 linked with Fortran according to the steps in https://www.researchgate.net/publication/349991987_Linking_ABAQUS_20192020_and_Intel_oneAPI_Base_Toolkit_FORTRAN_Compiler。