[论文解读] Computational modeling of microstructure

Mitchell Luskin|arXiv (Cornell University)|May 1, 2003

Shape Memory Alloy Transformations参考文献 37被引用 11

一句话总结

本文提出了一套计算框架，用于通过具有细尺度振荡的能量最小化变形的数值分析，模拟马氏体和铁磁性晶体中的微结构。该框架引入了一种成核算法，通过在每个时间步以概率方式引入奥氏体或马氏体相，实现准静态相变模拟，并通过稳定量（如Young测度和非线性积分）的误差估计进行了验证。

ABSTRACT

Many materials such as martensitic or ferromagnetic crystals are observed to be in metastable states exhibiting a fine-scale, structured spatial oscillation called microstructure; and hysteresis is observed as the temperature, boundary forces, or external magnetic field changes. We have developed a numerical analysis of microstructure and used this theory to construct numerical methods that have been used to compute approximations to the deformation of crystals with microstructure.

研究动机与目标

• 开发一种对马氏体和铁磁性晶体中微结构的数值分析方法，以捕捉细尺度振荡变形。
• 构建在网格细化下变形梯度点态不稳定的情况下仍保持稳定的数值方法，用于计算微结构。
• 通过在温度或载荷变化期间实现亚稳相的成核，模拟相变的准静态演化。
• 确保在网格细化下，体积分数（Young测度）和非线性积分等稳定量的收敛性。
• 通过在时间离散框架中引入温度依赖的成核概率，模拟相变中的滞后行为。

提出的方法

• 将能量泛函表述为框架无关弹性能量密度 $φ(F,\theta)$ 的积分，在低温下于马氏体变体 $\mathcal{U}$ 上进行最小化。
• 对 $\phi(F,\theta)$ 施加二次增长和大变形增长条件，以确保能量有限且变形具有一致连续性。
• 在理想化分析中，将允许的变形空间 $\mathcal{A}$ 定义为具有狄利克雷边界条件的形式，并忽略载荷和界面能。
• 使用有限元空间 $\mathcal{P}_1(\tau) \times \mathcal{P}_0(\tau)$ 在三角剖分区域上逼近变形和相指标 $b$。
• 采用时间离散格式，每个时间步的变形通过从成核初始迭代出发的Polak-Ribière共轭梯度法计算。
• 引入一种基于温度依赖概率函数 $P(\theta)$ 的随机成核策略，根据局部状态和温度在单元上初始化奥氏体或马氏体相。

实验结果

研究问题

• RQ1 当变形梯度在网格细化下不稳定时，如何构建微结构的稳定数值逼近？
• RQ2 何种数值框架能够实现对马氏体晶体在不同温度和载荷下相变滞后行为的精确计算？
• RQ3 如何在时间演化过程中避免陷入局部能量极小值，从而模拟亚稳态？
• RQ4 何种条件可确保在微结构数值逼近中，Young测度和非线性积分等稳定量的收敛性？
• RQ5 在准静态相变过程中，如何以物理解释的方式建模新相的成核？

主要发现

• 即使在点态梯度于网格细化下不稳定的情况下，仍为稳定量（如非线性积分 $\int_\Omega f(x,\nabla y(x))\,dx$ 和局部体积分数，即Young测度）建立了误差估计。
• 该数值方法通过能量最小化变形的有限元逼近成功捕捉了马氏体变体间具有振荡梯度的微结构形成。
• 成核算法通过在每个时间步基于温度和局部状态以概率方式引入奥氏体或马氏体相，实现了相变的模拟。
• 该方法再现了马氏体相变中的滞后行为，包括加热过程中高温相的成核和冷却过程中低温相的成核。
• 通过在“帐篷”构型下对CuAlNi薄膜的模拟验证了该方法，其中通过基于梯度的准则更新相指标 $b^{[0]}(x_K,t_\ell)$ 实现成核。
• 该算法可针对升温与降温过程采用不同的概率函数 $P(\theta)$，并可在特定温度区间内禁止逆相变。