I've done a model of a simple [0/90/45/-45]_S plate under uniform traction Nx=100.

When I use a Elastic Material with Fail Stress I get the correct failure index calculated by max. stress and Tsai-Wu.

When I use a Elastic Material with Hashin Damage (no Fail Stress values given, but the same values given as Hashin parameters), I get failure indexes about 100 times smaller, and for some layers I get zero. This is obviously wrong.

I must note that I omit to provide the Damage evolution properties, on purppose because the manual says,... quote...

21.3.2. ... If you define a damage initiation model without defining an associated evolution law, the initiation criteria will affect only output. Thus, you can use these criteria to evaluate the propensity of the material to undergo damage without modeling the damage process. end quote.

Also, note that the Fail stresses in compression XC and YC must be entered as negative in Fail Stress, but must be enteres as positive in Hashin Damage.

Does anyone know if the Hashin Failure works in 6.10-2?... or have any idea what I might be doing wrong.

I add part of the .inp file in case you can figure this out:

model here

** Section: Section-1*Shell Section, elset=_PickedSet2, composite, symmetric1.25, 3, Material-Hashin, 0., a1.25, 3, Material-Hashin, 90., b1.25, 3, Material-Hashin, 45., c1.25, 3, Material-Hashin, -45., d*End Instance** 

more cards here

** MATERIALS** note in *shell section that I am not using this material *Material, name=Material-Fail*Elastic, type=ENGINEERING CONSTANTS133860., 7706., 7706., 0.301, 0.301, 0.396, 4306., 4306. 2760.,*Fail Stress1830.,-1096.,   57., -228.,   71.,   -1.,    0.**

** note in *shell section that I am using this material *Material, name=Material-Hashin*Damage Initiation, criterion=HASHIN, alpha=1.1830., 1096.,   57.,  228.,   71., 21.48*Elastic, type=ENGINEERING CONSTANTS133860., 7706., 7706., 0.301, 0.301, 0.396, 4306., 4306. 2760.,**

** BOUNDARY CONDITIONS** ** Name: BC-1 Type: Symmetry/Antisymmetry/Encastre*Boundary_PickedSet5, XSYMM** Name: BC-2 Type: Symmetry/Antisymmetry/Encastre*Boundary_PickedSet6, YSYMM** ----------------------------------------------------------------** ** STEP: Step-1** *Step, name=Step-1*Static1., 1., 1e-05, 1.** ** LOADS** ** Name: Load-1   Type: Shell edge load*Dsload_PickedSurf4, EDNOR, -100.** ** OUTPUT REQUESTS** *Restart, write, frequency=0** ** FIELD OUTPUT: F-Output-1** *Output, field*Node OutputCF, PHILSM, PSILSM, RF, U*Element Output, directions=YES1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16CFAILURE, DAMAGEC, DAMAGEFC, DAMAGEFT, DAMAGEMC, DAMAGEMT, DAMAGESHR, DAMAGET, DMICRT, E, ERPRATIO, HSNFCCRT, HSNFTCRT, HSNMCCRT, HSNMTCRT, JKS, SDEG, SHRRATIO*Element Output, directions=YES17, 18, 19, 20, 21, 22, 23, 24CFAILURE, DAMAGEC, DAMAGEFC, DAMAGEFT, DAMAGEMC, DAMAGEMT, DAMAGESHR, DAMAGET, DMICRT, E, ERPRATIO, HSNFCCRT, HSNFTCRT, HSNMCCRT, HSNMTCRT, JKS, SDEG, SHRRATIO*Contact OutputBDSTAT, CRSTS, CSDMG, CSMAXSCRT, CSMAXUCRT, CSQUADSCRT, CSQUADUCRT, DBS, DBSF, DBT, EFENRRTR, ENRRT, OPENBC** ** HISTORY OUTPUT: H-Output-1** *Output, history, variable=PRESELECT*End Step