在用黃永剛子程序huang.for做晶體塑性有限元分析時，提交計算后log文件顯示奇異矩陣，模型的位移邊界條件用一個子模型驅(qū)動，相同的模型不使用umat可以計算，只要加上umat就會有奇異矩陣。