1,首先利用geom模块生成包含30个晶粒,模型空间尺寸为32*32*2
代码:
2,根据布种区域生成voronoi镶嵌的晶体结构
代码:
3,检查晶体结构是否正确生成 代码:
此时文件中会现20.grain.vtr文件,我们可以用
paraview可视化多晶板材模型
4,利用damask施加载荷并求解
代码:
在此处我们需要两个个文件,分别是
,tension.load 用来表示载荷的施加
,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科技工作室