Defining a new state variable

[mjm1985Manu]

Hi DAMASK team,

I am using DAMASK v3.0.1 and would like to define a state variable, D, as a damage indicator for an uncoupled ductile damage criterion. My intention is to delete or remove a voxel when D = 1 (with D = 0 representing an undamaged state and D = 1 representing full damage). For values of D < 1, no softening or stiffness degradation is needed, since this is an uncoupled model.

Could you please let me know whether I need to modify the current source codes (e.g., phase_damage.f90, homogenization_damage.f90, phase_damage_isobrittle.f90) or if it would be possible to define this criterion elsewhere, such as in phase.f90?

Thank you in advance for your guidance.