1，首先利用geom模塊生成包含30個晶粒，模型空間尺寸為32*32*2

代碼：

2，根據布種區域生成voronoi鑲嵌的晶體結構

代碼：

3，檢查晶體結構是否正確生成 代碼：

此時文件中會現20.grain.vtr文件，我們可以用

paraview可視化多晶板材模型

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

代碼：

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

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

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

模型的選擇）

Tension.load 文件的內容為

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

0 time 200 incs 400 freq 1，

material.config的文件內容包含

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

material.config的文件內容包含

2，指定材料本構模型，這里彈性模型采用hooke線彈性模型，塑性使用唯象的晶體塑性本構模型，晶體結構為Fcc，包含12個滑移系，其余參數含義見官網

material.config的文件內容包含

3，每一個晶粒的取向，采用bunge方式描述，在文件中給出每一個晶粒的歐拉角

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

代碼：

求解進行中和求解成功的示意圖如下

5，生成后處理文件夾，指定要輸出的內容這里我們指定應力代碼：

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

6，對后處理數據進行細分，便于輸出完整過程，表征整個演化流程

7，轉到后處理文件夾之中

8,生成Vtk文件，便于使用paraview后處理相關數據

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

最后,有相關需求歡迎通過公眾號聯系我們

公眾號:320科技工作室