AbaqusUmat子程序_非线性弹性

成竹

各位大神好，正在学习编程中。我咨询的问题如下：
在Abaqus的子程序Umat中编写非线性弹性的程序，用的是Fortran语言。我现在有应力应变的方程，但是求解一致切线刚度矩阵需要求导，一致切线刚度矩阵是6*6的矩阵，目前我不会求导（因为维数多，若是一维问题就没啥问题），也不会在Fortran中实现这个程序的编写。麻烦大神指导或者推荐参考资料。盼回复。谢谢。