MATLAB代码实现凝固相场模拟——凝固模型，各向异性枝晶生长纯物质凝固模型，激光增材制造...

vQFQJbUiJ

261人浏览 · 2026-04-04 19:45:00

vQFQJbUiJ · 2026-04-04 19:45:00 发布

MATLAB代码实现凝固相场模拟——凝固模型，各向异性枝晶生长纯物质凝固模型，激光增材制造，选择性激光熔融，SLM，凝固，铸造，焊接等等多种耦合场景提供讲解video~提供讲解video~提供讲解video~ 相场原理推导数学过程的视频讲解和参考文献资料以及模型文件作为在相场模拟领域深耕多年的过来人，跟大家分享些心得体会[拥抱]相场模拟总的来说难度大，系统复杂网络上没有有针对性的，可实践的编程那么作为初学者，步入这个领域就非常困难，有门槛不管怎么看paper，没有上手的编程实践，根本不可能掌握[帐篷] 楼主也是从新手过来[拥抱]，经过多少个无休止的探索，多少个不眠之夜[度假]，总算对这个领域有了一点心得体会所以总结出了这一套编程代码，以及凝固相场的原理推导讲解，从金兹堡-郎道自由能出发，教会大家如何定向，精准的掌握相场理论希望大家都能有收获~

枝晶在金属凝固过程中舞动的模样，总能让我想起显微镜下冰晶生长的优雅轨迹。十年前第一次在实验室看到Al-Cu合金枝晶照片时，我就被这种自然形成的数学之美深深吸引。相场法作为模拟这类现象的利器，其编程实现却让无数初学者在矩阵的迷宫里兜兜转转——今天我们就用MATLAB撕开这层神秘面纱。

先看这段各向异性处理的核心代码：

function anisotropy = calc_anisotropy(theta, epsilon, mode)
    % theta: 界面法向角度
    % epsilon: 各向异性强度
    % mode=1 六重对称 / mode=2 四重对称
    if mode == 1
        delta = epsilon * cos(6*theta); 
    else
        delta = epsilon * cos(4*theta); 
    end
    anisotropy = 1 + delta;
end

参数epsilon控制着枝晶尖端分叉的剧烈程度，就像调节显微镜的焦距——当epsilon=0.04时，你会看到典型的六瓣雪花状生长；调到0.1以上，尖端开始出现肉眼可见的震荡。这种角度依赖的界面能变化，正是枝晶呈现规则对称性的物理本质。

相场方程求解的时空离散化是另一个难点。来看时间步进的核心循环：

for n = 1:max_steps
    % 计算相场梯度
    [phi_x, phi_y] = gradient(phi, dx);
    
    % 各向异性计算
    theta = atan2(phi_y, phi_x);
    A = calc_anisotropy(theta, 0.04, 1); 
    
    % 构建相场方程右端项
    term1 = A.^2 .* (lap_phi - 2*(phi_x.^2 + phi_y.^2).*d2A_dtheta2);
    rhs = mobility * (term1 - df_dphi);
    
    % 显式时间推进
    phi = phi + dt * rhs;
    
    % 实时可视化
    if mod(n,50)==0
        imagesc(phi); axis equal off; drawnow
    end
end

这里采用显式欧拉法虽然简单，但要注意CFL条件——就像煮粥时火候的控制，时间步长dt必须小于dx²/(2*D)，否则数值震荡会让你的枝晶像沸腾的粥一样炸锅。曾有个博士生因为忽略这个条件，模拟出的枝晶长得像外星生物，在组会上被调侃了整整三个月。

当处理激光熔池快速凝固时，需要在模型中耦合温度场：

% 激光热源项
laser_source = P/(2*pi*r^2) * exp(-((X-x0).^2 + (Y-y0).^2)/(2*r^2));

% 热传导方程
kappa = 1e-5; % 热扩散系数
T = T + dt*(kappa*del2(T) + laser_source - latent_heat*dphi_dt);

这个热源项像数码画笔，在金属粉末床上绘制熔池轨迹。参数r控制着光斑尺寸，当设置为50μm时，模拟出的单道熔池形貌与我们的SLM实验观测结果误差小于5%。有趣的是，调整latent_heat参数时，会发现凝固收缩导致的凹陷深度变化——这解释了为什么实际焊接中需要精确控制冷却速率。