1，首先利用geom布種子，種子數(shù)為40，模型空間尺寸為64*64*1（表示二維平面問題）代碼：

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

代碼：

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

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

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

0time  200  incs400 freq  1，

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

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

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

模型的選擇）

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

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

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

2，指定材料本構(gòu)模型，這里彈性模型采用hooke線彈性模型，塑性使用唯象的晶體塑性本構(gòu)模型，晶體結(jié)構(gòu)為Hcp，不同基面包含不同個數(shù)的滑移系，及對應(yīng)的參數(shù)，要分別給出，同時需要給出孿晶的對應(yīng)參數(shù)，其余參數(shù)含義見官網(wǎng)

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

代碼：

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

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

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

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

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

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

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

最后，工作室damask和abaqus聯(lián)合仿真培訓(xùn)火熱進行中，需要的小伙伴歡迎來聊。

Damask和abaqus晶體塑性聯(lián)合仿真培訓(xùn)通知