vpsc-vumat子程序的开发与应用文章推荐_Abaqus_UM_理论_VPS

vpsc-vumat子程序的开发与应用文章推荐

原始文献：《A multiscale simulation framework of the accumulative roll bonding process accounting for texture evolution》

DOI：10.1016/j.msea.2015.02.005

在金属材料研究领域，粘塑性自洽模型（Visco-Plastic Self-Consistent, VPSC） 是研究织构演变与力学响应的经典之作。不同于简单的 Taylor 模型（假设所有晶粒应变相同）或 Sachs 模型（假设应力相同），VPSC 将每个晶粒视为嵌入在“等效介质”中的椭球夹杂。它巧妙地平衡了晶粒间的应力与应变分配，既考虑了晶粒形貌的影响，也能精确捕捉由于晶体转动引起的织构演变。

对于从事铝合金、镁合金等具有显著各向异性材料研究的同学来说，VPSC是预测材料在复杂加工路径下表现的有力工具。然而，原生的 VPSC 通常是针对均匀变形设计的，面对实际工程中复杂的几何边界和非均匀变形（如轧制、冲压），它需要一个更强大的载体。

Abaqus 作为有限元分析（FEA）的标杆，擅长处理复杂的边界条件和几何接触。将 VPSC 以 VUMAT（用户材料子程序） 的形式集成进 Abaqus，能实现“1+1 > 2”的效果，例如宏微观耦合：每一个有限元积分点都代表一个多晶集合。有限元计算宏观应变，VPSC 在微观层面计算晶体旋转和硬化，再反馈回宏观应力。非均匀场预测：你不仅能看到工件的整体变形，还能清晰地观察到厚度方向、圆周方向上织构分布的异质性。复杂工艺仿真：只有融入有限元，才能真正模拟非对称轧制等具有复杂应力状态的工艺。

今天推荐的是Prakash 等人在 Materials Science & Engineering A 上发表的经典论文。该论文针对累积叠轧（ARB）中，材料每道次减薄 50%，网格在两三道次后就会严重畸变。此外，层数成倍增加，微观状态如何继承的问题，提出了一种状态变量映射技术。在网格畸变前，通过插值算法将织构（取向）、晶粒形状（变形梯度）等信息转移到新网格。

这保证了材料“记忆”的连续性。同时论文采用了应力驱动的自协调迭代，并引入了两级并行计算（MPI + OpenMP），这在 2026 年依然是非常经典的设计。

作者成功捕捉到了 ARB 厚度方向上的织构梯度（中心 S 组分与表面剪切组分）。

以及形貌的演化特征：

轧制的局部应力状态：

作者的模拟结果表明：ARB 过程中，上一道次的表面（剪切区）在叠轧后进入下一道次的中心，导致织构在厚度方向上不断重新分布和细化。同时“两级并行”比单一并行模式在处理这类复杂多晶模型时具有压倒性的时间优势。

结合作者的理论（尤其是分段线性化和应力驱动的求解思路）我们可以把独立的vpsc子程序编写进abaqus里面，为了避免复杂的雅可比推导，以及适用各种复杂的变形工况，推荐使用abaqus的显式求解器，即vumat程序

以下展示一个使用vpsc-镁合金本构模型，模拟包含1个单元，单元包好100个晶粒在RD方向压缩20%的模拟效果（原始模型参数取自vpsc官方案例,为了减少计算时间使用高应变率进行计算，应力高于vpsc模拟）：

等效塑性应变：

第一个晶粒的累计剪切滑移：

发生孪晶次数;

变形后的形状演化：

来源：我的博士日记

镁合金孪晶建模经典文章推荐

文章名称：《A constitutive model for hcp materials deforming by slip and twinning: application to magnesium alloy AZ31B》DOI：10.1016/S0749-6419(03)00039-1在镁合金晶体塑性建模的发展脉络里，2003 年 Staroselsky 这篇文章并不是最新的，也不是形式最复杂的。但如果真正从“建立一个能跑、能解释现象、还能继续扩展的模型”这个角度去看，它仍然是一篇非常值得反复阅读的经典工作。这篇文章研究的对象，是 AZ31B 镁合金在室温条件下的塑性变形行为。作者关注的问题非常明确：为什么这种 HCP 结构材料在不同加载方向下，会表现出强烈的不对称性、明显的织构演化，以及非常突出的孪晶效应？换句话说，这篇文章不是简单去拟合一条应力—应变曲线，而是试图回答：镁合金在室温下究竟是靠哪些机制在变形，这些机制又如何共同决定宏观响应。这篇工作的建模核心思想：第一，它非常强调孪晶不是附属机制，而是主导机制之一。在 HCP 镁合金里，单靠常规滑移并不能解释很多室温下的实验现象，尤其是压缩拉伸不对称和织构快速变化。Staroselsky 这篇文章抓住了这一点，把孪晶明确地放到了与滑移同等重要的位置。第二，它采用了一种非常有启发性的处理方式：把孪晶作为一种具有方向性的 pseudo-slip 来处理。这样做的好处是，一方面保留了孪晶“有极性、不可正反完全对称”的物理特征，另一方面又能把它自然嵌入有限变形晶体塑性框架中。这个思想到今天看仍然非常高明，因为它在“物理真实性”和“程序可实现性”之间找到了很好的平衡。第三，这篇文章并没有急着把硬化写得非常复杂。恰恰相反，作者采用了一个相对简洁的 non-hardening 假设，认为材料宏观上表现出来的许多“硬化感”，很大程度上其实来自取向变化和孪晶导致的晶格重排，而不只是每个滑移系的阻力不断增加。这个判断非常重要，因为它提醒我们：在 HCP 镁合金里，单纯盯着“硬化参数”往往是不够的，织构和重取向本身就是塑性响应的重要组成部分。第四，作者还加入了一个各向同性 accommodation 项，用来描述晶界附近非晶体学协调变形的作用。这个处理非常值得重视。很多时候我们做晶体塑性，只把目光放在晶内滑移和孪晶上，却容易忽略多晶材料中晶粒之间并不是天然完全协调的。Staroselsky 这篇文章清楚地认识到：如果不考虑这部分效应，数值计算中的应力水平会偏高，甚至难以合理匹配实验。因此，这个附加项虽然形式上不复杂，但在建模思想上非常成熟。主要考虑的滑移和孪晶如下：拉伸变形的实验于模拟结果对比：压缩变形的模拟和实验结果对比：从结果上看，这篇文章得到的结论也非常有代表性。作者指出，AZ31B 镁合金室温变形的主要载荷承担机制，不是单一滑移，而是 basal、pris matic、pyramidal <a> 滑移与 {10-12}<10-11> 孪晶共同作用的结果。同时，他们认为材料在变形过程中表现出的显著各向异性和应变硬化特征，很大一部分来源于变形诱导的织构演化，而孪晶又是其中最关键的驱动力之一。这个结论今天看依然不过时，因为它直接抓住了镁合金室温塑性的本质。推荐文章的最主要原因是：2003 的这个文章的价值，不在于它把所有机制都做全了，而在于它先把最重要的几个问题讲清楚了——谁在变形、谁在重取向、谁在影响应力水平。对于刚开始做 HCP 晶体塑性的人来说，这种建模路径非常值得学习。孪晶不能只作为“后处理现象”看待，而应该进入本构主框架。如果模型里没有真正把孪晶放进去，那么很多镁合金室温响应就只能停留在表面拟合，难以解释深层原因。一个好模型不仅要能解释，还要能扩展。 Staroselsky 这篇文章的框架虽然简洁，但非常适合作为后续更复杂模型的基础。无论是进一步加入显式硬化、温度效应，还是发展到更复杂的 twin-detwin、parent-child 思路，它都可以作为一个非常扎实的起点。后续很多孪晶模型基于此进行二次开发，因此实现该文章的数值模型对于孪晶的研究非常有帮助：使用文章的公式，讲整体算法集成到abaqus的vumat子程序相对容易，因为不需要推导一致性雅可比。但是率无关模型通常数值稳定性较差。使用该本构模型模拟效果如下：初始RVE模型：沿着X方式施加拉伸变形，变形结束后应力分布：变形后的孪晶体积分数：阈值孪晶体积分数（文章中推荐0.3-1.0的随机值）：变形结束是否发生孪生变形：来源：我的博士日记