局部材料坐标下的应变/应力

我们正在进行简单的双层复合材料的模拟，它经历了广泛的屈曲和后屈曲不稳定。模拟是用ABAQUS显式完成的，以捕获高度非线性的几何形状。一切都运行得很好，但是与我们的分析理论相比，我们需要在每个元素上使用Green应变。实体单元的ABAQUS只给出NE或LE投影到全局坐标系上。当然，主应变是可用的，但是我们感兴趣的是局部旋转模式，特别是应变映射到局部共旋转参考坐标系。ABAQUS显然可以访问元素的每个积分点的变形梯度张量F(非对称的3x3拉伸比矩阵)，我如何将参考坐标系的单位向量(x=[1 0 0]， y=[0 1 0]， z=[0 0 1])映射到变形状态下的新单位向量x'， y'， z'。X ' = X。y' = y。每个元素的Fz' = z ?我读过这是可能的VUMAT子程序，但这些似乎只用于特殊的本构关系;系统是线性弹性的。 Is there a way without evoking the subroutines that the above mapping can be done? Thank you for any help you can offer.