1，首先利用geom模塊布種子，種子數(shù)為40，模型空間尺寸為10*10*10

2，根據(jù)布種區(qū)域生成voronoi鑲嵌的晶體結(jié)構(gòu)

3，檢查晶體結(jié)構(gòu)是否正確生成

此時文件中會現(xiàn)30.grain.vtr文件，我們可以用paraview可視化多晶板材模型

4，利用damask施加載荷并求解

在此處我們需要兩個個文件，分別是

1. ，tension.load 用來表示載荷的施加

2. ，material.config 用來表示材料的基本屬性（晶體取向，本構(gòu)

模型的選擇）

Tension.load 文件的內(nèi)容為

代碼：fdot 1.0e-3 0 0 0 * 0 0 0 * stress * * * * 0 * * *

0 time 100 incs 200 freq 1

material.config的文件內(nèi)容包含

1，指定分析步的輸出內(nèi)容（包含晶體本身的輸出，和phase的輸出）我們指定的輸出如下圖所示

material.config的文件內(nèi)容包含

2，指定材料本構(gòu)模型，這里彈性模型采用hooke線彈性模型，塑性使用唯象的晶體塑性本構(gòu)模型，晶體結(jié)構(gòu)為Bcc，包含24個滑移系，分別指定馬氏體，鐵素體的對應(yīng)材料參數(shù)

material.config的文件內(nèi)容包含

3，每一個晶粒的取向，采用bunge方式描述，在文件中給出每一個晶粒的歐拉角，以及每一個晶粒種兩相結(jié)構(gòu)所占的比值（這里假設(shè)兩者的占比為0.82，0.18）

完成后保存兩個文件在同一個工作目錄，保證晶粒個數(shù)，以及歐拉角與生成的幾何體的晶體相同

求解進(jìn)行中和求解成功的示意圖如下

5，生成后處理文件夾，指定要輸出的內(nèi)容這里我們指定應(yīng)力代碼：

生成后處理文件夾及成功示意圖

6，對后處理數(shù)據(jù)進(jìn)行細(xì)分，便于輸出完整過程，表征整個演化流程

7，轉(zhuǎn)到后處理文件夾之中

8,生成Vtk文件，便于使用paraview后處理相關(guān)數(shù)據(jù)

9，將vtr格式的文件用paraview打開，如下圖所示，P和F的分布情況

最后,有相關(guān)需求,歡迎通過公眾號聯(lián)系我們.

公眾號:320科技工作室.