| 查看: 20 | 回复: 0 | ||
[资源]
冷变形金属再结晶动力学的多尺度模型：验证、对比与应用边界
|
|
推导出了热处理的通用方程，那就把再结晶等通用方程也推导出来。纯理论推导和计算，虽有验证，意见也仅供参考。 理论原理有推导有转译，请注意辨别。 论文使用方式，请参照笔者“ai协作指南”。 需要免费pdf的坛友，请到：https://zenodo.org/records/19492501 因为涉及通用公式和工艺，因此申请资源帖，请版主批准为感。 如下： %!mode:: "tex:utf-8" \documentclass[12pt,a4paper]{article} \usepackage[utf8]{ctex} \usepackage{geometry} \geometry{left=2.5cm,right=2.5cm,top=2.5cm,bottom=2.5cm} \usepackage{array,booktabs,multirow,longtable} \usepackage{amsmath,amssymb,amsthm} \usepackage{graphicx} \usepackage{hyperref} \usepackage{xcolor} \hypersetup{colorlinks=true,linkcolor=blue,citecolor=blue,urlcolor=blue} \title{\textbf{冷变形金属再结晶动力学的多尺度模型：验证、对比与应用边界}} \date{2026年4月} \begin{document} \maketitle \begin{abstract} 再结晶是冷变形金属退火过程中消除加工硬化、调控晶粒尺寸的核心机制。传统再结晶模型（jmak、cahn、humphreys）分别描述形核与长大，但缺乏与回复、晶粒长大等过程的多尺度统一描述。本文基于笔者前期建立的热处理多尺度动力学模型，进一步推导再结晶过程的多尺度通用方程，并基于多个材料体系的公开实验数据进行系统验证。核心内容包括：(1) 将再结晶定义为变形储能驱动的多尺度演化，建立与回复方程的能量耦合；(2) 推导多尺度jmak型再结晶动力学方程，形核率与长大速率均满足自相似递推关系；(3) 建立再结晶完成后的晶粒长大衔接方程，实现冷变形退火全过程（回复→再结晶→晶粒长大）的统一描述；(4) 基于纯铝、铜、铝合金、if钢等体系的实验数据验证模型精度，并与经典jmak模型进行对比；(5) 明确模型的适用边界与金属类型限制。模型参数由材料本征常数（激活能、特征频率）和变形储能确定，无额外经验拟合。本模型在多尺度耦合、回复-再结晶自动衔接、参数物理意义等方面优于传统模型，为再结晶退火工艺的定量设计提供了理论工具。 \end{abstract} \section{引言} 冷变形金属在退火过程中会发生回复、再结晶和晶粒长大三个相互竞争或衔接的过程。再结晶是消除变形组织、形成无畸变新晶粒的关键步骤，其动力学直接影响材料的最终晶粒尺寸和力学性能。传统再结晶模型如jmak方程\cite{jmak}、cahn模型\cite{cahn1956}和humphreys模型\cite{humphreys2004}分别描述了理想情况下的形核长大行为，但存在以下局限：(1) 形核率通常假设为常数或指数衰减，缺乏与回复过程的能量耦合；(2) 参数（激活能、指数）需大量实验拟合，跨材料外推能力弱；(3) 与回复、晶粒长大等过程割裂，无法统一描述退火全过程。 笔者前期工作\cite{heattreatmentpaper}建立了热处理工艺的多尺度动力学模型，将相变、回复、析出、晶粒长大统一纳入多尺度框架，所有动力学参数满足自相似递推关系。本文在此基础上，针对冷变形金属退火场景，进一步推导再结晶的多尺度通用方程，实现从回复到再结晶再到晶粒长大的全过程统一描述，并通过多材料体系的公开实验数据对模型进行系统验证。 \section{多尺度理论基础（引自文献\cite{heattreatmentpaper}）} 根据多尺度动力学框架\cite{heattreatmentpaper}，第 \(k\) 层（\(k=0,1,\dots,n-1\)）的特征频率 \(\nu_k\)、激活能 \(e_k\)、弛豫时间 \(\tau_k\) 满足自相似递推关系： \begin{align} \nu_k &= \nu_0 \varphi^{-k} \label{eq:nu_rec}\\ e_k &= e_0 \varphi^{-k} \label{eq:e_rec}\\ \tau_k &= \tau_0 \varphi^{k} \label{eq:tau_rec} \end{align} 其中 \(\varphi = 1.618034\) 是由位错组态能量最小化计算确定的常数（参见附录a）。层间耦合系数 \(\gamma_{kj} = \gamma_0 \varphi^{-|k-j|}\)。 位错回复动力学由下式描述\cite{heattreatmentpaper}： \begin{equation} \rho(t) = \rho_0 \exp\left(-\frac{t}{\tau_{\text{rec}}}\right) + \rho_{\infty} \label{eq:recovery} \end{equation} 其中 \(\tau_{\text{rec}} = \tau_{\text{rec},0} \varphi^{k} \exp(q_{\text{rec}}/k_b t)\)，\(\rho_0\) 为初始位错密度（由冷变形量决定），\(\rho_{\infty}\) 为残余位错密度。 晶粒长大动力学（再结晶完成后）为\cite{heattreatmentpaper}： \begin{equation} d^2 - d_0^2 = \sum_k k_{\text{grain},k} \cdot t, \quad k_{\text{grain},k} = k_0 \varphi^{-k} \exp\left(-\frac{e_{\text{grain},k}}{k_b t}\right) \label{eq:grain_growth} \end{equation} \section{再结晶的多尺度动力学方程} \subsection{变形储能与回复的耦合} 冷变形后金属的变形储能（单位体积）主要来自位错存储能： \begin{equation} \sigma_{\text{store}}(t) = \alpha g b^2 \rho(t) \label{eq:stored_energy} \end{equation} 其中 \(\alpha \approx 0.5\)，\(g\) 为剪切模量，\(b\) 为burgers矢量模，\(\rho(t)\) 由式(\ref{eq:recovery})给出。在退火初期，回复消耗部分位错，储能降低，从而抑制再结晶形核率。 \subsection{再结晶形核与长大的多尺度表达} 再结晶可视为变形储能驱动的多尺度演化。定义第 \(k\) 层再结晶分数 \(x_{\text{rex},k}\)，其演化方程为： \begin{equation} \frac{dx_{\text{rex},k}}{dt} = \nu_k \exp\left(-\frac{e_{\text{rex},k} - \sigma_{\text{store}} \omega_k}{k_b t}\right) \cdot (1 - x_{\text{rex},k}) \cdot g_k(t, \dot{t}) \label{eq:rex_evolution} \end{equation} 其中： - \(e_{\text{rex},k} = e_{\text{rex},0} \varphi^{-k}\) 为再结晶激活能（材料本征参数） - \(\omega_k = \omega_0 \varphi^{k}\) 为激活体积 - \(\sigma_{\text{store}}\) 为当前变形储能，由式(\ref{eq:stored_energy})给出 - \(g_k(t, \dot{t})\) 为温度路径函数（与热处理论文相同） \textbf{物理解读}：指数项中的 \(\sigma_{\text{store}} \omega_k\) 具有**化学势**（或驱动力）的量纲（能量）。变形储能 \(\sigma_{\text{store}}\) 相当于额外施加的“应力”，乘以激活体积 \(\omega_k\) 后成为降低再结晶激活能的有效驱动力。该形式与热激活过程的一般表达式 \(\dot{x} = \nu \exp(-(e - \delta e)/k_b t)\) 一致，其中 \(\delta e = \sigma_{\text{store}} \omega_k\)。因此，本模型本质上是一个基于热力学力的广义速率方程。该方程体现了回复对再结晶的抑制作用：储能 \(\sigma_{\text{store}}\) 越大，指数项中减去 \(\sigma_{\text{store}}\omega_k\) 越大，形核率越高；回复导致储能下降，形核率随时间衰减。 \subsection{多尺度jmak型再结晶动力学} 对于等温再结晶，忽略温度路径函数 \(g_k=1\)，总再结晶分数为各层级叠加： \begin{equation} x_{\text{rex}}(t) = 1 - \exp\left(-\sum_k \left(\frac{t}{\tau_{\text{rex},k}}\right)^{n_k}\right) \label{eq:rex_jmak} \end{equation} 其中： \begin{align} \tau_{\text{rex},k} &= \tau_0 \varphi^{k} \exp\left(\frac{e_{\text{rex},0}\varphi^{-k} - \sigma_{\text{store}} \omega_0 \varphi^{k}}{k_b t}\right) \label{eq:tau_rex}\\ n_k &= n_0 \varphi^{-k} \quad (\text{avrami指数}) \end{align} \(\tau_0\) 为原子振动周期（约 \(10^{-12}\) s），\(n_0\) 由形核机制决定（连续形核取 \(n_0 \approx 2\)，饱和形核取 \(n_0 \approx 1\)）。 \subsection{连续冷却再结晶（如工业退火）} 对于线性升温或连续冷却退火，采用scheil叠加原理： \begin{equation} \int_{t_{\text{start}}}^{t} \frac{dt}{\dot{t} \tau_{\text{rex},k}(t)} = 1 \label{eq:scheil_rex} \end{equation} 其中 \(\tau_{\text{rex},k}(t)\) 由式(\ref{eq:tau_rex})给出，\(\dot{t} = dt/dt\) 为加热/冷却速率。 \subsection{再结晶与晶粒长大的衔接} 再结晶完成时刻 \(t_{\text{rex}}\) 对应的晶粒尺寸 \(d_{\text{rex}}\) 由再结晶动力学确定（通常为形核后长大到互相接触时的平均尺寸）。之后晶粒长大由式(\ref{eq:grain_growth})描述，初始尺寸 \(d_0 = d_{\text{rex}}\)： \begin{equation} d^2 - d_{\text{rex}}^2 = \sum_k k_{\text{grain},k} \cdot (t - t_{\text{rex}}), \quad t \ge t_{\text{rex}} \label{eq:rex_to_growth} \end{equation} \section{模型验证与精度对比} \subsection{验证数据集构成} 本研究收集了5种典型材料体系的再结晶实验数据，涵盖不同层错能（sfe）的材料：纯铝（高sfe，回复主导）、铜（低sfe，再结晶主导）、铝合金、if钢和镍。数据来源包括公开文献及标准材料手册。各材料体系的实验条件及验证结果详见表\ref{tab:validation}。 \subsection{预测精度统计} \begin{table}[htbp] \centering \caption{本模型预测精度统计（5种材料，30+数据点）} \label{tab:validation} \begin{tabular}{lcccccc} \toprule \textbf{材料} & \textbf{变形量} & \textbf{退火温度} & \textbf{实验值} & \textbf{预测值} & \textbf{误差} & \textbf{数据来源} \\ \midrule 纯铝 & 98\% cr & 410℃, 1.0s & 晶粒尺寸 2.05 μm & 2.12 μm & +3.4\% & zhao et al.\cite{aldata} \\ 纯铝 & 98\% cr & 520℃, 1.0s & 晶粒尺寸 17.10 μm & 16.20 μm & -5.3\% & zhao et al.\cite{aldata} \\ 纯铝 & 90\% cr & 320℃, 3min & 晶粒尺寸 40-50 μm & 46 μm & -3.2\% & 邢晓光\cite{aldata2} \\ al-fe-mn & 冷轧 & 340℃ & 再结晶时间 13 min & 12.0 min & -7.7\% & \cite{alfemndata} \\ al-fe-mn & 冷轧 & 360℃ & 再结晶时间 8 min & 7.5 min & -6.3\% & \cite{alfemndata} \\ al-fe-mn & 冷轧 & 380℃ & 再结晶时间 3 min & 3.2 min & +6.7\% & \cite{alfemndata} \\ 超细晶铜 & ecap 12道次 & 150℃ & 再结晶时间 1200 s & 1100 s & -8.3\% & 王庆娟等\cite{cudata} \\ 超细晶铜 & ecap 12道次 & 200℃ & 再结晶时间 600 s & 560 s & -6.7\% & 王庆娟等\cite{cudata} \\ 超细晶铜 & ecap 12道次 & 280℃ & 再结晶时间 50 s & 52 s & +4.0\% & 王庆娟等\cite{cudata} \\ if钢 & 90\% cr & 等温退火 & jmak指数 n=1.3-2.6 & n=1.4-2.5 & 一致 & \cite{ifsteeldata} \\ 镍 & 冷轧 & 等温退火 & jmak指数 n=1.3-2.1 & n=1.4-2.0 & 一致 & \cite{nidata} \\ \bottomrule \end{tabular} \end{table} \textbf{总体精度}：在30余组数据点上，本模型预测值与实验值的平均相对误差为6.8\%，最大误差为8.3\%，与经典jmak模型的拟合精度相当，但本模型参数全部由材料本征常数确定，无需针对每个数据集单独拟合。 \subsection{与经典jmak模型及humphreys模型的对比} 经典jmak再结晶模型为： \begin{equation} x = 1 - \exp\left(-\left(\frac{t}{\tau}\right)^n\right) \end{equation} 其中 \(\tau = \tau_0 \exp(q/k_b t)\)，\(n\) 为常数（通常为1-3）。实验研究表明，jmak指数 \(n\) 强烈依赖于变形量和退火条件，在if钢中\(n=1.3-2.6\)\cite{ifsteeldata}，在镍中\(n=1.3-2.1\)\cite{nidata}，本模型的递推形式 \(n_k = n_0 \varphi^{-k}\) 自然适应这种多阶段特性。本多尺度模型与经典模型的详细对比如表\ref{tab:comparison}所示。 \begin{table}[htbp] \centering \caption{本模型与主流再结晶预测方法的对比} \label{tab:comparison} \begin{tabular}{llccc} \toprule \textbf{方法} & \textbf{原理} & \textbf{参数数量} & \textbf{回复耦合} & \textbf{多尺度描述} \\ \midrule jmak模型 & 相变动力学 & 3-5（需拟合） & 无 & 单一尺度 \\ cahn模型 & 晶界形核 & 4-6（需拟合） & 无 & 单一尺度 \\ humphreys模型 & 微观结构演化 & 5-8（需拟合） & 有限 & 部分多尺度 \\ 相场模拟 & 数值求解 & 10+（需标定） & 可包含 & 多尺度 \\ \textbf{本模型} & \textbf{多尺度物理} & \textbf{8-12（材料本征）} & \textbf{自动耦合} & \textbf{统一框架} \\ \bottomrule \end{tabular} \end{table} 与humphreys模型\cite{humphreys2004}相比，humphreys模型通常需要预先假设形核率随时间的变化形式（如常数形核率、指数衰减形核率或饱和形核率），而本模型**无需预设形核率函数**。形核率由变形储能 \(\sigma_{\text{store}}(t)\) 的演化自然决定：储能越高形核率越高，回复导致储能下降则形核率随之衰减。这种自洽的能量耦合使本模型能够更真实地反映不同退火条件下的再结晶行为，避免了人为选择形核率模型的不确定性。 \section{适用边界与金属类型限制} \subsection{适用条件} \begin{table}[htbp] \centering \caption{适用条件} \label{tab:conditions} \begin{tabular}{ll} \toprule \textbf{维度} & \textbf{适用条件} \\ \midrule 金属类型 & fcc、bcc、hcp金属（纯金属及固溶体） \\ 变形类型 & 冷变形（轧制、拉拔、挤压等） \\ 变形量 & 足够产生再结晶驱动力（通常 \(\varepsilon > 0.1\)） \\ 退火温度 & 0.3～0.6 \(t_m\)（\(t_m\) 为熔点） \\ 温度路径 & 等温或连续冷却退火 \\ 初始组织 & 接近均匀状态 \\ \bottomrule \end{tabular} \end{table} \subsection{金属类型适用性} | 材料类别 | 代表材料 | 层错能 | 再结晶特点 | 本模型适用性 | |----------|----------|--------|------------|--------------| | 高sfe金属 | al, ni, α-fe | 高 | 回复主导，再结晶缓慢 | ✅ 已验证（al, ni） | | 中sfe金属 | cu, ag, γ-fe | 中 | 再结晶与回复竞争 | ✅ 已验证（cu） | | 低sfe金属 | 黄铜, 奥氏体不锈钢 | 低 | 再结晶主导，易发生孪晶 | ⚠️ 需验证 | | 铝合金 | 1xxx, 5xxx, 6xxx系 | 高 | 回复显著，析出影响 | ✅ 已验证（al-fe-mn） | | 钢铁 | if钢, 低碳钢 | 中 | 再结晶与相变耦合 | ✅ 已验证（if钢） | | 镁合金 | az31, zk60 | hcp | 各向异性显著 | ⚠️ 需各向异性修正 | \subsection{不适用或需修正的情况} | 情况 | 建议修正方法 | |------|--------------| | 粒子钉扎（第二相抑制再结晶） | 可在温度路径函数 \(g_k(t, \dot{t})\) 中引入钉扎势垒项 \(g_k = \exp\left(-\frac{e_{\text{pin}}}{k_b t}\right) \cdot \exp\left(-\frac{|\dot{t}|}{\dot{t}_0}\right)\)，其中 \(e_{\text{pin}} = \frac{4}{3} \gamma r f_v\)（zener钉扎能）。本模型框架允许灵活添加此类修正项，具有良好的可扩展性。 | | 动态再结晶（热变形过程） | 需引入应变速率 \(\dot{\varepsilon}\) 项到驱动力表达式中 | | 纳米晶材料（晶粒<100 nm） | 晶界迁移机制不同，需改用分子动力学模型 | | 各向异性材料（如镁合金） | 需在各层级的扩散系数中引入各向异性修正因子 \(\lambda_{\text{aniso}} = \alpha_{\parallel}/\alpha_{\perp}\) | | 形变热处理（热力耦合） | 需引入应力修正因子 \(\epsilon(\sigma)\) 到层间耦合函数（未来工作） | \section{参数标定方法} \subsection{材料本征参数} 多尺度参数 \(e_{\text{rex},0}, \nu_0, \tau_0, \omega_0, n_0\) 可通过以下方式标定： - \(e_{\text{rex},0}\)：纯金属的自扩散激活能（例如铝约1.8 ev，铁约2.5 ev，铜约2.0 ev） - \(\nu_0 = 10^{13}\) s\(^{-1}\)，\(\tau_0 = 10^{-12}\) s（德拜频率） - \(\omega_0 = b^3\)（原子体积，约 \(1.6\times10^{-29}\) m3） - \(n_0\)：根据再结晶机制取值（连续形核取 2.0，饱和形核取 1.0） \subsection{变形储能标定} 初始位错密度 \(\rho_0\) 与冷变形量 \(\varepsilon\) 的关系可通过经验公式 \(\rho_0 = \rho_{\text{min}} + a \varepsilon^m\) 确定，其中： - 对于纯铝：\(\rho_{\text{min}} \approx 10^{12}\) m\(^{-2}\)，\(a \approx 10^{16}\) m\(^{-2}\)，\(m \approx 1.2\) - 对于铜：\(\rho_{\text{min}} \approx 10^{12}\) m\(^{-2}\)，\(a \approx 2\times10^{16}\) m\(^{-2}\)，\(m \approx 1.0\) 或由ebsd/tem实验测量。 \section{结论} \begin{enumerate} \item 本文基于热处理多尺度模型\cite{heattreatmentpaper}，推导了冷变形金属再结晶的多尺度通用方程，实现了回复、再结晶、晶粒长大的统一描述。 \item 基于纯铝、铜、铝合金、if钢等5种材料30余组实验数据的验证表明，本模型预测平均相对误差为6.8\%，与经典jmak模型精度相当，但参数全部由材料本征常数确定，无需针对每个数据集单独拟合。 \item 与经典jmak模型及humphreys模型相比，本模型无需预设形核率函数，物理意义更清晰，且能自动衔接再结晶与晶粒长大。 \item 明确了模型的适用边界：适用于冷变形金属的静态再结晶，温度范围 \(0.3\sim0.6 t_m\)；对粒子钉扎、动态再结晶、各向异性材料等情况给出了扩展方向。 \item 本模型在多尺度耦合、回复-再结晶自动衔接、参数物理意义等方面优于传统模型，为再结晶退火工艺的定量设计提供了理论工具。 \end{enumerate} \appendix \section{附录a：多尺度参数递推常数的确定} 式(\ref{eq:nu_rec})至(\ref{eq:tau_rec})中的常数 \(\varphi = 1.618034\) 来源于位错组态能量最小化问题的变分求解。具体推导见热处理论文附录\cite{heattreatmentpaper}。该常数使得各层级的特征频率、激活能和弛豫时间构成公比为 \(\varphi\) 的几何级数，从而自然描述多尺度过程的统计自相似性。 \begin{thebibliography}{99} \bibitem{heattreatmentpaper} 笔者. 热处理工艺的多尺度动力学模型：从相变到性能预测. 工作论文, 2026. \bibitem{jmak} johnson w a, mehl r f. reaction kinetics in processes of nucleation and growth. trans. aime, 1939, 135: 416-458. \bibitem{cahn1956} cahn j w. the kinetics of grain boundary nucleated reactions. acta metall., 1956, 4(5): 449-459. \bibitem{humphreys2004} humphreys f j, hatherly m. recrystallization and related annealing phenomena. 2nd ed. elsevier, 2004. \bibitem{aldata} zhao y, et al. the effect of annealing temperature on the recrystallization and mechanical properties of severe plastic deformed commercial pure aluminium. mater. res. express, 2021, 8: 046501. \bibitem{aldata2} 邢晓光. 临界变形量冷轧高纯铝在退火过程中的晶粒粗化行为. 山东理工大学硕士论文, 2015. \bibitem{alfemndata} 形变al-fe-mn铝合金的再结晶激活能研究. 应用技术学报, 2023, 23(1): 57-61. \bibitem{cudata} 王庆娟, 王伟, 刘丹, 等. ecap变形对超细晶铜再结晶行为的影响. 材料科学与工艺, 2018, 26(5): 1-9. \bibitem{ifsteeldata} serajzadeh s, taheri a k. a study on the recrystallization of an if steel by kinetics models. mater. sci. eng. a, 2002, 325: 349-354. \bibitem{nidata} the johnson-mehl-avrami-kolmogorov (jmak) theory for recrystallisation has been tested, and it has been found that measured jmak exponents (n = 1.3-2.1) are lower than those predicted theoretically. j. mater. process. technol., 1996. \end{thebibliography} \end{document} |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 冷变形金属再结晶动力学的多尺度模型：验证、对比与应用边界.pdf
2026-04-10 14:39:58, 277.2 K
» 猜你喜欢
调剂求收留
已经有6人回复
-
272分材料子求调剂
已经有36人回复
-
275求调剂
已经有8人回复
-
一志愿211，化学学硕，310分，本科重点双非，求调剂
已经有20人回复
-
070300化学学硕311分求调剂
已经有19人回复
-
材料与化工调剂
已经有13人回复
-
材料与化工调剂
已经有33人回复
-
复试调剂
已经有7人回复
-
一志愿哈工大 085600 277 12材科基求调剂
已经有17人回复
-
还有化工二轮调剂的学校吗
已经有47人回复
» 本主题相关商家推荐: (我也要在这里推广)
