I want to implement a constitutive model for ice in compression which is developed using damage mechanics. Damage mechanics in the model utilizes Schapery's approach, which in turn derived from non-linear viscoelastic theory.I also want to implement an explicit solver which can capture high-speed dynamic events accurately.Could you help me understand what it will take to implement this material model, what SfePy can already do, and what I will need to contribute to  for making this project work?

 

 This project is being carried out to know how easy is to use SfePy for the research that is going on in our university and to contribute for open source softwares like SfePy. The results will also be validated by comparing with the experimental.

 Please refer the attached paper for further details about the material model.