目前在研究循环内聚力模型的UMAT开发，但是在开发过程中遇到几个问题（算法用matlab尝试可行）
  第一个疑问：作为循环内聚力模型，一般把第一个循环的加载阶段视为单调载荷，包含了弹性阶段的应力线性变化和软化阶段的应力衰减，相对于循环载荷条件来说，各个积分点的应力更新公式不同，所以需要把单调和循环两个加载条件区分。在编写UMAT的过程中，打算加入一个与单元应变相关的状态变量作为单调载荷与循环载荷的判据，但是如何在逻辑上实现（对于同一个材料点，每一次调用时可能会将上述状态变量重复赋值，不能很清晰的判断加载状态）？

  第二个疑问，对于三维内聚力单元，只能输出S33,S13,S23三个方向的应力。在UMAT中是否只用更新上述三个方向的应力，且这三个应力是否用stress(1),stress(2),stress(3)就可以对应上。

  第三个疑问，与第二个问题类似，损伤演化方程中的应力应变为有效值，对于有效位移 dlta=sqrt(strant(1)**2+strant(2)**2+strant(3)**2)，三个应变值是否应该对应法向应变和对应的两个剪应变？ 主程序传递的stran(2) stran(3) 就是对应的两个剪应变吗？

  希望有过类似经验的老师或同行看到了能给我指点迷津，不胜感激！

目前在研究循环内聚力模型的UMAT开发，但是在开发过程中遇到几个问题（算法用matlab尝试可行）
  第一个疑问：作为循环内聚力模型，一般把第一个循环的加载阶段视为单调载荷，包含了弹性阶段的应力线性变化和软化阶段的应力衰减，相对于循环载荷条件来说，各个积分点的应力更新公式不同，所以需要把单调和循环两个加载条件区分。在编写UMAT的过程中，打算加入一个与单元应变相关的状态变量作为单调载荷与循环载荷的判据，但是如何在逻辑上实现（对于同一个材料点，每一次调用时可能会将上述状态变量重复赋值，不能很清晰的判断加载状态）？

  第二个疑问，对于三维内聚力单元，只能输出S33,S13,S23三个方向的应力。在UMAT中是否只用更新上述三个方向的应力，且这三个应力是否用stress(1),stress(2),stress(3)就可以对应上。

  第三个疑问，与第二个问题类似，损伤演化方程中的应力应变为有效值，对于有效位移 dlta=sqrt(strant(1)**2+strant(2)**2+strant(3)**2)，三个应变值是否应该对应法向应变和对应的两个剪应变？ 主程序传递的stran(2) stran(3) 就是对应的两个剪应变吗？

  希望有过类似经验的老师或同行看到了能给我指点迷津，不胜感激！