I want to model progressive impact damage in laminated composite, I see in papers, they usually use VUMAT subroutine, but there are other subroutines to model damage. Is there anyone who knows this problem? And is there any other way to model progressive damage?