DAMASK实例1–二维多晶铝

1.创建多晶模型

输入：

seeds_fromRandom -N 20 --grid 64 64 1 >20grains.seeds

上述命令随机生成20个seeds，64*64划分

geom_fromVoronoiTessellation --grid 64 64 1 20grains.seeds

上述命令通过20个种子生成几何模型20grains.geom

geom_check 20grains.geom

上述命令生成可视的vtr文件，通过Paraview，上方选项如下图所示：

则可以看到生成的晶体模型：

2.创建加载文件及材料属性文件

创建加载文件tension.load,如下图所示：

在（1,1,1)方向施加力，总时间200，增量步400，每1步均输出结果保存。

在材料属性文件material.config中，定义20个晶粒，只有一个相，且输出我们需要的变量：f:应力，p:应变，如下图所示：

Al的材料参数如下所示：（并未展示完全，详情见资料链接）

3.计算求解及后处理

接下来输入下述命令即可进行计算求解：

DAMASK_spectral --geometry 20grains.geom --load tension.load

对得到的数据进行后处理，输出应力应变，输入如下命令：

postResults 20grains_tension.spectralOut --cr f,p

可以看到多了一个PostProc文件夹，在里面存在一个txt格式的文件，存放在每个晶粒九个方向的应力应变如下图所示：

如果想输出更多所需结果，例如晶粒取向，织构等，则可输入下述命令：

postResults 20grains_tension.spectralOut --cr f,p,orientation,texture,grainrotation --separation x,y,z --increments --range 400 400 1 --split

得到20grains_tension_inc400.txt文件，如下图所示:

对计算结果进行可视化后处理，输入下述命令：

vtk_rectilinearGrid 20grains_tension_inc400.txt

可以得到vtr文件，将计算输出结果添加到该文件，输入：

vtk_addRectilinearGridData --data f,p,orientation,texture,grainrotation --vtk 20grains_tension_inc400_pos\(cell\).vtr 20grains_tension_inc400.txt

可以得到最终可视化结果的应力（f）应变（p），如下图所示（图中左上角标黄）：

相关资料链接：百度云 提取码：94fo

Fischerstardust于HIT(SZ)