GND/Gradient Initialization Silently Fails in Multi-Phase Models When Crystal Elements Do Not Start from Label 1

[Akashgiitk]

Dear Sir, ,
i found this unusual error/bug while using the code.
In multi-phase or two-phase ABAQUS models (e.g., elastic inclusion+ crystal plasticity matrix), the UMAT is only called for crystal phase elements. If ABAQUS assigns those crystal elements labels that do not start from 1 (e.g., elastic elements are 1–99, crystal elements are 100–199), the code produces three compounding errors:

1. numel is inflated to the highest element label (199) instead of the true crystal element count (100).
2. The GND gradient operator initialization condition (sum(ip_init) == numel*numpt) is never satisfied, so GNDs are silently never computed for the entire analysis.
3. All global state variable arrays (statev_gmatinv, statev_tauc, ipcoords, etc.) are over-allocated.

It is advisable to number the crystal elements beginning from 1 for the subset on which the GND computation is intended to be carried out.

Thankyou
Akash