hypela2 - subroutine isotropic hardening based on perturbation

傲慢与偏见

本子程序基于MARC的子程序Hypela2编写的。对于不同的材料本构，其应力和切线刚度矩阵都不同，因此需要修改本程序的代码。为便于后续程序的编写，对程序作了模块化设计。
程序包含：
1. 应力迭代子程序1（基于应力矢量）
2. 应力迭代子程序2（基于应力张量）
3. tangent modulus算法1（解析方法），更正了以前的一个typo error。
4. Tangent moudlus算法2（摄动理论，数值解法）
5. 应力应变摄动子程序
6. 状态变量结果输出子程序plotv.
7. 矢量张量转换子程序。
本程序特色在于，切线刚度矩阵与材料本构无关，因此理论上可以应用于任意材料本构。对于复杂材料本构，切线刚度矩阵难以得到解析解的情况下，基于摄动理论的数值解法是很有价值的。

peiss8888

谢谢分享！