? ? ? ? 為了將晶體塑性模擬結(jié)果與EBSD實(shí)驗(yàn)得到的反極圖map (IPF?map) 相對(duì)比，在晶體塑性有限元模擬完成后，需要利用每個(gè)高斯點(diǎn)的歐拉角和高斯點(diǎn)坐標(biāo)，畫(huà)出模擬的 IPF map。這就需要利用模擬所得歐拉角求出對(duì)應(yīng)的rgb值，然后用rgb顏色畫(huà)圖。在DAMASK中有相關(guān)的介紹，但是由于本人用abaqus UMAT進(jìn)行晶體塑性模擬，暫時(shí)沒(méi)有想到用abaqus python解決上述問(wèn)題的方法，所以這里先從.odb文件中將每個(gè)高斯點(diǎn)的坐標(biāo)及歐拉角提取出來(lái)，借助mtex生成每個(gè)歐拉角對(duì)應(yīng)的rgb，然后做截面，這樣就可以得到想要的截面上的ipf map了。

? ? ? ? 由于EBSD得到的也是試樣一個(gè)截面處上的取向信息，那么這樣對(duì)比起來(lái)應(yīng)該是比較方便的。同時(shí)由于matlab可以利用方程定義截面，所以甚至可以畫(huà)出一個(gè)曲面截面上的IPF map，這一點(diǎn)還未在其他地方見(jiàn)到過(guò)，是利用matlab進(jìn)行abaqus后處理的一個(gè)小的優(yōu)勢(shì)。

step -1 :? 從abaqus .odb文件中提取高斯點(diǎn)坐標(biāo)和歐拉角到csv文件，一個(gè)簡(jiǎn)單地python腳本

step -2?: 將高斯點(diǎn)坐標(biāo)及歐拉角讀到matlab里面進(jìn)行處理(b站代碼塊不支持matlab，直接貼圖吧不然很亂)

上面只是一個(gè)腳本，沒(méi)有寫(xiě)成函數(shù)，用起來(lái)可能后有一些問(wèn)題。自己做一些記錄也給可能需要的人提供個(gè)參考。