基于凌志matlab腳本實現(xiàn)拉伸壓縮，平面應(yīng)變壓縮過程中的織構(gòu)演化

具體操作過程如下：

一，在for文件中加入輸出內(nèi)容如下:

代碼：

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

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

位置如圖：

選擇特定增量步格式化輸出結(jié)果：

      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%），平面應(yīng)變壓縮（軋制變形）（75%），初始取向隨機，取向分布如圖所示：

模型包含512個單元，每個單元代表一個晶粒，并且擁有不同的取向，在abaqus">abaqus運算結(jié)束后，得到對應(yīng)的temp.dat文件。內(nèi)容如下圖所示：

得到數(shù)據(jù)后，打開matlab程序，依次運行matlab腳本，將數(shù)據(jù)導入mtex進行處理，運行順序為data5→ebsd6→Mtex7，注意文件路徑的修改，其中data5負責信息的讀入，ebsd6負責將信息轉(zhuǎn)化為mtex支持的格式。mtex7負責極圖的繪制。

得到的變形后的取向如下圖所示，與大多數(shù)文獻報道結(jié)果一致，重現(xiàn)了典型的織構(gòu)特征分布。從而驗證了腳本的正確性。