In the case of simulation of an iterative process (Let's say we have 30 iterations) in which elements' material properties need to be changed after each iteration based on the results (eg. resulting strain field) of the previous iteration, which method is preferred and why?

Thanks in advance.

PS1: An example of such a process would be a simulation of the 30 days of bone healing by 30 iterations.
PS2: Answers can also be added to the RG: https://www.researchgate.net/post/Abaqus_scripting...