首页/文章/ 详情

基于C++相场的枝晶凝固织构仿真

10分钟前浏览0

关键词:相场模拟;枝晶生长;织构预测;CUDA加速;增材制造

一、文章简要介绍

增材制造的凝固速度横跨10的负4次方到10米每秒,晶粒在极端热条件下互相竞争,最终留下的晶体取向决定了构件的各向异性性能,但织构怎么选出来一直缺定量工具。 Northeastern 大学 Karma 团队在 Nature Communications 上给出了一套组合方案:先在激光重熔 AlSi 薄膜的 DTEM 实验里量化织构,再用 C/C++ 写的 CUDA 相场程序在 GPU 上复现枝晶生长竞争,反推出固液界面自由能各向异性的真实参数(DOI:10.1038/s41467-025-66655-2)。对做凝固与织构仿真的工程师来说,这套从实验到相场的闭环很有参考价值。

图1 DTEM拍摄的熔池动态、凝固速度曲线与熔池边界IPF图。

图2 熔池边界晶粒竞争IPF时间序列,不利取向晶粒逐步被淘汰。

二、仿真步骤

步骤一:实验数据获取。用动态透射电镜实时拍摄 Al1Si、Al3Si、Al6Si 薄膜的熔池凝固过程,实测凝固速度 0.5~2 m/s,落在增材制造特征区间;再用 ASTAR 取向成像配合 MATLAB 2024a 的 MTEX 工具箱完成晶粒取向标定,熔池边界按椭圆界面逐帧追踪。

步骤二:织构统计。对每个晶粒建立局部坐标系,把三次对称的三个特征方向投影到凝固方向上求和对齐量,用一千万个随机取向晶粒建立统计分布,再以 Z 评分衡量实测织构偏离随机的程度,显著性阈值取1.96,对应 p 值0.05。

图3 不同硅含量与凝固阶段的实验织构离散极图。

图4 SAST统计结果,随机分布与Z评分显著性判据。

步骤三:定量相场模拟。采用 Karma 定量相场模型描述枝晶生长,固液界面自由能各向异性用两参数形式,参数由分子动力学原子模拟给出;演化方程用 C/C++ 实现并跑在 CUDA 框架上,V100/H200 GPU 单精度计算,典型网格规模 5120×5120。

图5 相场模拟织构与实验Z评分、晶粒数的对比。

三、关键结果与发现

单参数各向异性只偏好立方〈100〉生长方向,模拟出的织构与实验对不上;换成两参数形式后,AlSi 薄膜里观察到的〈110〉主导织构被成功复现,说明即便在每秒数米的快速凝固下,枝晶生长方向仍由界面自由能各向异性主导。晶粒竞争过程也被定量还原:初始29个晶粒随固相分数推进被淘汰到13个、再剩5个,zs[100] 评分越过1.96阈值并逼近完美对齐值,与实验逐段吻合(DOI:10.1038/s41467-025-66655-2)。

图6 单参数各向异性下的模拟织构,与实验不符。

图7 两参数各向异性模拟成功复现实验织构走向。

来源:320科技工作室
增材MATLAB分子动力学
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-11-14
最近编辑:10分钟前
320科技工作室
硕士 | 结构工程师 微❤️ CAE320
获赞 243粉丝 449文章 532课程 0
点赞
收藏
作者推荐

基于Moldflow的薄壁壳件注塑工艺参数多目标优化仿真

关键词:Moldflow;薄壁壳件;注塑工艺;翘曲变形;多目标优化一、文章简要介绍汽车内饰里大量薄壁壳结构,注塑时体收缩和翘曲互相牵制,调参靠试错既慢又不稳。这篇江苏师范大学团队2026年发表在Scientific Reports的工作,用Moldflow仿真样本做数据底座,建立RIME-RF-MOGWO多目标优化框架,把两个质量指标同时压了下去(DOI:10.1038/s41598-026-65446-z)。带弯边、加强筋和开口特征的壳件,壁厚0.9到1.8毫米。二、Moldflow仿真建模零件在UG中建模,壁厚0.9到1.8毫米,最小壁厚与最大投影面积平方根之比只有0.00589,远低于0.045的薄壁判据。模型导入Autodesk Moldflow Insight 2023,双域三角网格离散,过渡区局部加密,最终31302个单元,匹配率94.2%。和86228单元的细网格相比,收缩率差异0.00%、翘曲差异0.04%,精度够、算得快。双域三角网格31302个单元,冷却水路环绕型腔布置。材料选PC/ABS合金,浇口、流道和冷却水路固定为模具配置,只把七个工艺参数当变量:注射压力、注射时间、保压压力、保压时间、熔体温度、模具温度、冷却时间。Moldflow按充填、保压、冷却、翘曲标准序列求解,冷却用瞬态传热,收缩按PVT行为评估,翘曲由残余应力导出。翘曲集中在边角区域,中间大片区域变形较小。全因子7因子5水平要78125次仿真,不现实;改用正交试验加补充采样,共完成95次Moldflow仿真。95个样本平均体收缩率15.45%、平均翘曲0.6000毫米,没有一个组合能同时满足收缩12%到13%、翘曲小于0.3毫米——必须做多目标优化。收缩分布同样不均,边角与浇口远端偏高。三、关键结果与发现随机森林置换重要性加ANOVA筛选显示,注射时间、熔体温度、保压时间贡献最大,分别占44.21%、23.42%、21.50%;注射压力和模具温度贡献为负,固定在80兆帕和70摄氏度,维度降到五个。SMOTE补齐稀疏区样本,RIME优化RF超参数,两个代理模型测试集R平方分别0.93880和0.89401。注射时间权重最高44.21%,压力和模温贡献为负。RIME优化超参数,目标值50代内逐步收敛。MOGWO在五维空间搜出帕累托前沿,代表解回到Moldflow复算:初始候选在边界外推区偏差大,把复算结果当反馈样本迭代修正后,预测与复算明显吻合。最终折中解Run 107:体收缩率12.51%、翘曲0.2962毫米,比原始样本平均值分别降19.02%和50.63%,双双达标。种群逐代逼近低收缩低翘曲区域,菱形为最终前沿。迭代后RF预测与Moldflow复算趋势吻合,Run 107为最终解。来源:320科技工作室

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈