在用黄永刚子程序huang.for做晶体塑性有限元分析时，提交计算后log文件显示奇异矩阵，模型的位移边界条件用一个子模型驱动，相同的模型不使用umat可以计算，只要加上umat就会有奇异矩阵。