低率相关性的晶体塑性模型，不同迭代方案计算时间的差异性

常规的唯象晶体塑性模型的流动方程通常使用幂律形式：

其中m为率相关系数，对于较小的m值，如≤0.01，整体的响应结果被认为接近率无关响应，然而该参数显著影响积分效率，对于不同的迭代方案，其对稳定性的影响也不仅相同，这里尝试进行简单对比，对比指标和总计算增量步数和计算时间（所有程序均使用单核计算）：

所有迭代方案使用相同的硬化模型和相同材料参数，并对包含200个晶粒的多晶模型进行20%的拉伸变形模拟。如下图所示：

（1）对于以弹性变形梯度和塑性速度梯度为迭代变量的写法：

计算时间：

增量步数：

（2）以PK2应力和滑移系统当前强度为迭代变量的双重迭代全隐式迭代方案：

计算时间：

增量步数：

（3）以PK2应力和滑移系统当前强度为迭代变量的双重迭代半隐式迭代方案：

计算时间：

增量步数：

（4）以PK2应力和滑移系统当前强度为迭代变量的单次迭代求解方程组全隐式迭代方案：

计算时间：

增量步数：

（5）以滑移系剪切应变为迭代变量的迭代方案：

无法完成模拟的收敛！！！

增量步数：

（6）以柯西应力为迭代变量的迭代方案：

计算时间：

增量步数：

（7）以偏应力为迭代变量的迭代方案：

计算时间：

增量步数：

模拟得到的效果图：

来源：我的博士日记

晶体塑性每日文章推荐（十八）

文章doi：10.1016/j.ijplas.2009.10.004上一篇推文介绍了基于L2范数最小化计算GND的推文，基于类似的思想可以实现率无关数值模型的构建，作者的创作思想就是利用塑性变形过程中最大能量耗散原理，将单晶屈服问题视作约束优化问题，其中约束就是每一个滑移系统的屈服函数，并将传统率无关晶体塑性模型中的数值奇异问题，通过创建约束条件的组合进行优化分析得到率无关的晶体塑性数值解。该率无关的本构计算方法计算效率相对较高，且数值稳定性很好，与以往研究和实验结果具有良好的一致性。此外也有很多率相关模型使用奇异值分解进行数值求解。作者的思想弹塑性问题通常被定义为约束优化问题，旨在寻找给定应变增量的最佳应力张量和内部变量。在这样的问题中，目标函数是基于最大耗散原理定义的，约束是屈服函数。示意图和公式为：以塑性变形率方程为切入点：λ为一致性参数在单晶中，整体塑性变形是多个滑移系统上滑移的结果。在晶体塑性问题中，变形由多个屈服面定义，这些屈服面不平滑相交，屈服函数的数量取决于晶体中滑移系的数量。假设施密德定律对于单晶的塑性变形是有效的，那么对于任何滑移系统，屈服函数可以定义为引入约束：在计算该方程过程中经常出现雅可比的非正定性问题（可以采用奇异值分解）同时确定滑移系统开动也存在计算量过大的问题作者这里引入组合约束思想，则原始求解模型转化为：在此基础上将多个约束转化为等效的一个约束： h（x）取值为：移植到所研究的晶体塑性模型中等效的约束为：将等效的约束作为单晶的屈服面：根据流动法则，计算得到塑性变形率其中剪切应变率作者所提出的模型对于不同ρ时的屈服面形状如下并且根据以往研究金属的层错能对屈服面有影响，因此将ρ与层错能进行关联,其中Γ是材料的 SFE，G是剪切模量，b是伯格斯矢量的大小。作者的模型在vumat子程序中进行实现，为了与以往文献对比，作者的硬化模型使用方程为：模拟了fcc-al以及fcc-coppor，材料参数为（ρ=80）模拟得到的效果可以看到作者提供的数值模型具有良好的精度，同时大大减少了计算时间，并克服了由于滑移的多重性而导致的奇异性问题。来源：我的博士日记