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

显式与隐式有限元方法在晶体塑性分析中的比较与最佳实践


引言

晶体塑性分析作为模拟金属和金属基材料介观尺度下大变形与应变局部化的重要工具，其数值计算涉及高度非线性的积分过程。文章《Comparison of the implicit and explicit finite element methods using crystal plasticity》详细探讨了在晶体塑性模型中，显式方法与隐式方法的性能、适用条件及其在实际应用中的策略选择。该文不仅提供了两种求解器的根本区别概述，还通过实验对比和案例分析，揭示了在不同的工程条件下，哪种方法更优。

显式与隐式方法区别

显式与隐式求解器基于不同的数值积分方法。在隐式方法中，通常采用牛顿迭代或类似的方法求解时间连续方程，要求计算收敛到精确解。这一过程虽可靠准确，但在接触问题等情况下，具有较高的时间步选择限制。相比之下，显式方法基于欧拉向前传播法，每一时间步可以直接计算出结果，因此更适合大规模并行计算和处理复杂接触条件，特别是当计算稳定性和效率成为关键因素时。




计算限制与优化

显式方法计算过程中面临的临界稳定时间问题，基于Le（单元特征长度）与Cd（膨胀波速）来定义时间步大小。增加Le或提高材料密度将导致更长的稳定时间，进而提高计算效率。在准静态分析中，显式求解器的使用可能引入额外的时间缩放，以加快计算过程，关键在于保持动能与总能量的比例不超过5%，从而确保本构响应的准确性和稳定性。

跨程序接口与性能提升

作者采用黄的亚弹性框架，将CUBIT隐式程序接口与ABAQUS/explicit显式程序结合，实现了显式与隐式积分过程的一致性。通过合理设计子程序接口，如通过VUMAT（用于ABAQUS/explicit）提供矢量化接口，来优化数据传递，提高了多核处理器计算效率。此外，基于时间增量分析，准确计算并调整初始时间增量以匹配材料应力响应，确保了响应的准确性与计算效率的平衡。通过质量缩放策略的引入，在保证物理行为准确性的前提下，显着提高了计算效率。

实验案例与案例分析

文章通过一系列实验案例比较了二维与三维变形条件、接触问题以及多晶体分析下的显式与隐式性能。依据具体问题的特点，结果显示了显式方法在处理复杂接触条件、滑动模拟及多核计算环境下的显著优势。同时，通过对等效应力分布、累计剪切应变及滑移系统当前强度等指标的对比分析，客观呈现了不同方法在特定问题场景下的性能差异。

研究结论与建议

显式方法在接触和力学非线性较强烈的场合展现出明显优势，特别是在大尺度变形条件下，其并行计算效率及对接触问题的稳定性较好地体现了其应用价值。然而，对于复杂参量下的计算，显式方法可能需要额外的时间缩放，这在一定程度上限制了其纯粹的“隐式”精度。因此，研究人员和工程师在决定采用显式或隐式求解器时，应综合考虑计算效率、性能需求以及物理问题的具体特性，适时引入优化策略，如质量缩放、合理选择时间步大小等，以实现在特定应用领域内的最优化解决方案。

武汉格发信息技术有限公司，格发许可优化管理系统可以帮你评估贵公司软件许可的真实需求，再低成本合规性管理软件许可,帮助贵司提高软件投资回报率，为软件采购、使用提供科学决策依据。支持的软件有: CAD,CAE,PDM,PLM,Catia,Ugnx, AutoCAD, Pro/E, Solidworks ,Hyperworks, Protel,CAXA,OpenWorks LandMark,MATLAB,Enovia,Winchill,TeamCenter,MathCAD,Ansys, Abaqus,ls-dyna, Fluent, MSC,Bentley,License,UG,ug,catia,Dassault Systèmes,AutoDesk,Altair,autocad,PTC,SolidWorks,Ansys,Siemens PLM Software,Paradigm,Mathworks,Borland,AVEVA,ESRI,hP,Solibri,Progman,Leica,Cadence,IBM,SIMULIA,Citrix,Sybase,Schlumberger,MSC Products...