Hello!

Please indicate in which aspects of modeling and analyzing the plasticity of HCP crystals I made a mistake. Please specify exactly what the problems may be in the modeling and analysis of crystalloplasticity in HCP. Every time you try to start, the error "The L-2 norm of the residual force overflowed. This may be caused by PRED,ON or birth/dead elements. Please use PRED,OFF and try again." appears. I've already worked hard on this task. I am attaching the APDL code below. The model used is a 0.02 mm cube - single-element and single-grain. The options "Large deformationand "Substeps" are enabled in the analysis settings with a minimum step of 100 and a maximum of 50,000.

Data source: Introduction to Crystal Plasticity Modeling in Ansys Mechanical

https://innovationspace.ansys.com/product/intro-to-crystal-plasticity-modeling-in-ansys-mechanical/?template=rocky

The APDL ANSYS code

/prep7

! Begin writing grain materials

*SET,matid,1

TB,ELAS,matid,,9,

TBDATA,1,70000.0,70000.0,70000.0,50724.637681159424,50724.637681159424,50724.637681159424

TBDATA,7,0.38,0.38,0.38

TB,CTE,matid,

TBDATA,1,1.28e-05,1.28e-05,1.28e-05

TB,PLAS,matid

TB,XTAL,matid,,3,ORIE

TBDATA,1,156.57144,0.0,0.0

TB,XTAL,matid,,1,NSLFAM

TBDATA,1,5

TB,XTAL,matid,,1,FORM

TBDATA,1,1

TB,XTAL,matid,,10,XPARAM

TBDATA,1,2,0,1,30,0,1

TBDATA,7,1,1,1,1

TB,XTAL,matid,1,18,HARD

TBTEMP, 298.0,

TBDATA,1,220.0,220.0,220.0,220.0,220.0,1750.0

TBDATA,7,1750.0,1750.0,1750.0,1750.0,400.0,400.0

TBDATA,13,400.0,400.0,400.0,2,0,1.4

TB,XTAL,matid,1,3,FLHCP

TBTEMP, 298.0,

TBDATA,1,1732000.0,0.005,1

CMSEL,S,GID1

EMODIF,ALL,MAT,matid

ALLSEL,ALL

allsel,all

TREF, 298.0

TUNIF, 298.0

/solu

pred,off

solv

fini