晶体塑形计算蠕变变形

umat编写，最好是基于huang的umat本构。我想用基于单晶的蠕变通过改变不同单元的晶体取向获得多晶的蠕变本构。