Your objective is to implement a plasticty return algorithm for the case of a Drucker-Prager model with isotropic linear hardening plasticity. You can use the provided Mathematica code as a template to modify or otherwise as a resource for creating your own implementation.
Mathematica Template Code (.nb)
Mathematica Template Code (.pdf)
Use resonable material properties and produce convincing plots that assist in the verification of your model.
Submit a working version of your code and sufficient verification evidence to Canvas by Friday, December 4, 11:00 AM.