基于lingzhi-matlab脚本与huang-umat实现变形过程中的织构演化预测

基于凌志matlab脚本实现拉伸压缩，平面应变压缩过程中的织构演化

具体操作过程如下：

一，在for文件中加入输出内容如下:

代码：

  FILE1='C:\Users\Desktop\VPSC\cp-huang-al\texture2\temp.dat'

      OPEN(91,FILE=FILE1,STATUS='UNKNOWN',ACCESS='APPEND')

位置如图：

选择特定增量步格式化输出结果：

      if(kinc==34) then

      write (91,100)  kinc, noel, npt, coords(1),coords(2),coords(3),

     2            (statev(IDNOR),IDNOR=3*NSLPTL+1,3*NSLPTL+3),

     3            (statev(IDDIR),IDDIR=6*NSLPTL+1,6*NSLPTL+3)

 100   FORMAT (1X,I3,I5,I2,F10.4,F10.4,F10.4,F30.22,F30.22,F30.22,

     1                      F30.22,F30.22,F30.22) 

      endif

      close(91)

分别模拟简单拉伸（75%），简单压缩（75%），平面应变压缩（轧制变形）（75%），初始取向随机，取向分布如图所示：

模型包含512个单元，每个单元代表一个晶粒，并且拥有不同的取向，在abaqus">abaqus运算结束后，得到对应的temp.dat文件。内容如下图所示：

得到数据后，打开matlab程序，依次运行matlab脚本，将数据导入mtex进行处理，运行顺序为data5→ebsd6→Mtex7，注意文件路径的修改，其中data5负责信息的读入，ebsd6负责将信息转化为mtex支持的格式。mtex7负责极图的绘制。

得到的变形后的取向如下图所示，与大多数文献报道结果一致，重现了典型的织构特征分布。从而验证了脚本的正确性。