对于Damask3.0版本感兴趣的同学可以扫描下面图片加入知识星球交流学习。

Damask2.0.3版本是最后一个和Abaqus有接口的版本，在Damask2.0.3的官网中已经进行了声明，目前最新的版本已经不支持和Abaqus联合使用。但是Damask和abaqus联合使用仍然是学习晶体塑性有限元方法不错的工具。我曾使用过Abaqus联合damask平台，这里我对使用过程中个人的一点经验进行简单的介绍，希望各位读者在研究过程中少一些技术上的障碍。

这么做的好处是：微结构演化/晶体塑性细节交给 DAMASK，复杂几何与边界条件交给成熟的 FEM。老用户要注意：官方早在 v2 时代就发出弃用 Abaqus 支持的提醒——2.0.3 是最后一个支持 Abaqus 的版本，之后不再维护。

前面介绍了damask2.03版本和abaqus配合使用，以及damask3.0版本和marc配合使用的案例介绍，当前推文的主要介绍如何把damask子程序正确移植到windows平台下并被abaqus正确调用，把 DAMASK 跑在 Windows 上这件事，看起来像是在“折腾编译环境”，但它真正的价值并不是让程序“能跑”，而是让晶体塑性与位错密度这类微观机制模型真正进入一个更高效、更可复现、更贴近工程实际的工作流

四、培训内容（1），abaqus调用damask实现FCC织构演化模拟------以多晶铁铝为例子（2），abaqus调用damask实现BCC织构演化模拟------以多晶铁素体为例子（3），abaqus调用damask实现HCP织构演化模拟------以多晶镁为例子（4）FCC，BCC，HCP多晶局部应力应变场模拟，状态变量说明，初始取向赋予，后处理取向提取

四、培训内容（1），abaqus调用damask实现FCC织构演化模拟------以多晶铁铝为例子（2），abaqus调用damask实现BCC织构演化模拟------以多晶铁素体为例子（3），abaqus调用damask实现HCP织构演化模拟------以多晶镁为例子（4）FCC，BCC，HCP多晶局部应力应变场模拟，状态变量说明，初始取向赋予

检查晶体结构是否正确生成 代码：geom_check 30grains.geom 利用damask施加载荷并求解 4，代码：DAMASK_spectral --geometry 30grains.geom --load tension.load 在此处我们需要两个个文件（分别设置用于求解），分别是 （1） tension.load------- 用来表示载荷的施加 （2

1，首先利用geom布种子，种子数为40，模型空间尺寸为64*64*1（表示二维平面问题）代码： 2，根据布种区域生成voronoi镶嵌的晶体结构 代码： 3 检查结构是否正确生成 4，利用damask施加载荷并求解 代码：fdot 1.0e-3 0 0 0 * 00 0 *stress * * ** 0

，分别为：不用系统对应的硬化表示为：作者考虑了3组滑移+1组孪晶共24个系统，这里说明一点的是，在孪生区域对应的滑移包含的个数为6*18，即每一组滑移都会因为特定的孪晶方向而旋转。

1，首先利用geom模块生成包含30个晶粒，模型空间尺寸为32*32*2代码：2，根据布种区域生成voronoi镶嵌的晶体结构代码：3，检查晶体结构是否正确生成 代码：此时文件中会现20.grain.vtr文件，我们可以用paraview可视化多晶板材模型4，利用damask施加载荷并求解代码：在此处我们需要两个个文件

且DAMASK同时支持两种方法，可以通过使用类似的模型进行对比，研究两类方法计算结果的差异。这里以DREAM 3D　为建模软件，建立两类模型，模拟DP双相钢在周期性边界体条件下整体的应力应变响应以及织构演化情况。

1，首先利用geom模块布种子，种子数为40，模型空间尺寸为10*10*102，根据布种区域生成voronoi镶嵌的晶体结构3，检查晶体结构是否正确生成此时文件中会现30.grain.vtr文件，我们可以用paraview可视化多晶板材模型4，利用damask施加载荷并求解在此处我们需要两个个文件

相关做法完美的集中到damask3.0版本里面，然而需要指出的是：DAMASK/谱方法更偏向规则网格与RVE范式，而工程里经常需要：任意几何与复杂边界（非周期、接触、局部细化等），以及不同工艺路径（多道次、换向、局部约束），Abaqus CPFEM（UMAT/VUMAT）在这些方面更“通用”，所以把“remesh + 状态变量映射”做成一套工作流，就能把大变形晶体塑性更稳地推进到更高压缩/更大应变阶段

相比于6.5X版本的Dream 3D，Dream 3D 7.0版本已从头开始完全重写，代码库的重写产生了一个新的底层库，新的用户界面：新版本的重大更新是引入了可视化界面（效果类似于Paraview）而且可以对数值进行转化和相关后处理，非常适合Abaqus或damask的后处理：支持导入的数据结构包括：支持多种颜色渲染，渲染效果如图：同时新版本特别适合分析

Abaqus调用damask实现轧制变形中FCC，BCC织构演化分析 案例实操一 1，使用abaqus建立20*20*20（mm）的立方块 2，对立方块进行单元划分共包含1000个单元 3，假设每个单元代表一个单独的晶粒，通过脚本随机赋予每个单元材料属性 4，施加对应的边界提交（60%的下压量） 5，提交与后处理材料数据 包含1000个晶粒的有限元模型

*degree,cs,ss);plotPDF(ori,Miller({1,0,0},{1,1,0},{1,1,1},cs),'all') %绘制极图plotIPDF(ori,[vector3d.X,vector3d.Y,vector3d.Z],'all')%绘制反极图下图是EVOCD平台拉伸教程的模拟结果，与本教程的结果一致。

还有如果有同学需要damask和neper、ABAQUS等，请参考我的另一篇帖子——基于linux虚拟机的晶体塑性有限元计算平台（已安装damask和neper等）》http://www.jishulink.com/content/post/1287159【VMware】如何打开已有虚拟机：https://jingyan.baidu.com/article/6181c3e0c44ebd152ef15332

如考虑孪生和滑移两者的贡献，需修改Huang基于滑移原理的子程序，或者换DAMASK平台进行模拟， 链接里是DAMASK平台模拟镁合金的变形。图5.3 DAMASK平台模拟镁合金变形的参数设置下面，将采用Huang的基于滑移原理晶体塑性子程序，开展HCP晶格材料以滑移为主的塑性变形模拟。