1，首先利用geom布种子，种子数为40，模型空间尺寸为64*64*1（表示二维平面问题）代码：

2，根据布种区域生成voronoi镶嵌的晶体结构

代码：

3 检查结构是否正确生成

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

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

0time  200  incs400 freq  1，

在此处我们需要两个个文件，分别是

            （1） ，tension.load 用来表示载荷的施加

             （2） ， material.config------- 用来表示材料的基本属性（晶体取向，本构

模型的选择）

material.config的文件内容包含

1，指定分析步的输出内容（包含晶体本身的输出，和phase的输出）我们指定的输出如下图所示

material.config的文件内容包含

2，指定材料本构模型，这里弹性模型采用hooke线弹性模型，塑性使用唯象的晶体塑性本构模型，晶体结构为Hcp，不同基面包含不同个数的滑移系，及对应的参数，要分别给出，同时需要给出孪晶的对应参数，其余参数含义见官网

完成后保存两个文件在同一个工作目录，保证晶粒个数，以及欧拉角与生成的几何体的晶体相同

代码：

求解进行中和求解成功的示意图如下

5，生成后处理文件夹，指定要输出的内容这里我们指定应力代码：

生成后处理文件夹及成功示意图

6 对后处理数据进行细分，便于输出完整过程，表征整个演化流程

7 转到后处理文件夹之中

8,生成Vtk文件，便于使用paraview后处理相关数据

9，将vtr格式的文件用paraview打开，如下图所示，f和p的分布的分布情况