1. 正弦摆动焊接熔池仿真到底在算什么
如果你刚接触 Fluent 焊接仿真，大概率会先做一个「直线移动高斯热源」的案例：热源沿 X 方向匀速走，温度场跟着跑，看起来挺像那么回事。但真实焊接工艺里，焊枪往往不是一条直线走到底的，摆动焊、扫描焊、锯齿焊都很常见。摆动的作用是让热影响区更宽、熔池形态更可控、气孔和裂纹倾向更低。问题在于，一旦热源轨迹变成正弦摆动，Fluent 自带的那些固定热源边界条件就不够用了，必须靠 UDF 在每个时间步里实时更新热源中心位置。
这篇内容面向的就是这个场景：用 Fluent + UDF 跑通一个正弦摆动焊接熔池模型，材料用 AISI 1045 中碳钢，工件是三维钢板，热源是移动高斯面热流，运动方式是沿焊接方向前进的同时横向做正弦摆动。关注的结果包括温度场、液相分数、熔池形态、熔池深度变化和流动趋势。热源轨迹的核心公式就两行：
x_center(t) = v * t y_center(t) = A * sin(2πft)其中 v 是焊接速度，A 是摆动幅度，f 是摆动频率。别看公式简单，真正落到 UDF 里要处理的东西不少：面循环、线程指针、时间步获取、高斯热流分布、单位一致性。我试过把这几块拆开单独验证，比一上来就全塞进去要稳得多。
适合谁看：正在做 Fluent 移动热源仿真、焊接熔池 CFD、高斯热源 UDF、摆动焊接仿真、熔池温度场后处理的同学。如果你已经能跑通直线热源，想进一步做摆动轨迹，这篇可以直接跟着操作。如果你连 UDF 编译都没做过，建议先把 Fluent UDF 编译环境跑通再回来。
另外说明一下，仿真过程中会用到一些 AI 辅助工具来生成 UDF 骨架、检查边界条件设置、整理后处理脚本。这些工具如果各自用一套 Key，管理起来很乱，后面会讲怎么用 TaoToken 统一管理调用通道，让仿真辅助工具的配置集中在一处。
2. 前置准备：UDF 编译环境与 TaoToken 统一 Key 通道
2.1 Fluent UDF 编译环境检查
在写 UDF 之前，先确认你的 Fluent 能编译 UDF。Windows 下需要 Visual Studio 的 C++ 编译器，Linux 下需要 GCC。检查方法很简单，在 Fluent 里打开 User-Defined > Functions > Compiled，点 Build，如果弹出一堆编译错误，说明环境没配好。
一个常见的坑是 Fluent 版本和编译器版本不匹配。比如 Fluent 2022 R1 配 VS2019 一般没问题，配 VS2022 有时会报cl.exe not found。解决办法是在 Fluent 启动前设置好环境变量，或者在 Fluent 的 UDF 编译设置里手动指定编译器路径。
# Windows 下检查编译器是否可用 where cl # Linux 下检查 which gcc如果这两条命令找不到可执行文件，先去装编译器，别急着写 UDF。
2.2 为什么仿真辅助工具需要统一 Key 通道
做焊接熔池仿真的时候，我通常会同时用几个 AI 辅助工具：一个用来生成 UDF 骨架，一个用来检查边界条件设置是否合理，一个用来整理后处理脚本。如果每个工具都单独配 Key，时间一长自己都记不清哪个 Key 对应哪个工具，额度用完了也不知道。
TaoToken 的做法是提供一个统一的 API 通道，你只需要在 TaoToken 控制台创建一个 Key，然后各个工具都指向同一个 API 地址就行。这样额度、调用记录、模型切换都在一个地方管理。对于仿真这种需要反复调试的场景，统一通道能省不少事。
具体操作：先到 TaoToken 控制台创建一个 API Key，然后在各个 AI 辅助工具的配置里把 base_url 指向https://taotoken.net/api，把 api_key 填成你创建的那个。模型名称按工具要求填，比如claude-sonnet-4-20250514或者gpt-4o之类的，具体看 TaoToken 文档里支持的模型列表。
如果你只是想让 AI 帮你解释一段 UDF 代码或者检查热源公式，用模型对话就行。如果你要长期做仿真辅助，比如批量生成不同参数的 UDF、自动整理后处理数据，那用 Coding Plan 会更划算，因为调用频率高，按量计费容易超预算。
2.3 仿真文件目录结构建议
在开始写 UDF 之前，建议先把目录结构理清楚，不然后面文件一多容易乱：
sinusoidal_welding/ ├── case/ │ ├── mesh.msh │ └── welding.cas ├── udf/ │ ├── heat_source.c │ └── compile.log ├── post/ │ ├── temperature_contour.py │ └── melt_pool_depth.py └── results/ ├── mp4/ └── curves/UDF 源文件单独放一个目录，编译日志也留着，后面排查编译错误的时候有用。后处理脚本单独放，不要和 case 文件混在一起。
3. 可复制配置：正弦摆动热源 UDF 与动网格参数
3.1 UDF 源码骨架
下面这个 UDF 骨架可以直接复制到heat_source.c里，然后根据你的模型尺寸改参数。核心逻辑是在每个时间步计算热源中心位置，然后对面上的每个网格单元计算高斯热流。
#include "udf.h" /* 焊接工艺参数 */ #define WELD_SPEED 5.0e-3 /* 焊接速度 m/s */ #define OSC_AMPLITUDE 3.0e-3 /* 摆动幅度 m */ #define OSC_FREQ 2.0 /* 摆动频率 Hz */ #define BEAM_RADIUS 2.0e-3 /* 热源有效半径 m */ #define Q_MAX 1.5e8 /* 最大热流密度 W/m^2 */ /* 热源中心位置计算 */ static void get_heat_source_center(real time, real *x_c, real *y_c) { *x_c = WELD_SPEED * time; *y_c = OSC_AMPLITUDE * sin(2.0 * M_PI * OSC_FREQ * time); } /* 高斯热流分布 */ static real gaussian_heat_flux(real r) { return Q_MAX * exp(-3.0 * r * r / (BEAM_RADIUS * BEAM_RADIUS)); } /* 面热源 UDF：在每个时间步更新热流 */ DEFINE_PROFILE(sinusoidal_heat_flux, thread, position) { face_t f; real time = CURRENT_TIME; real x_c, y_c; get_heat_source_center(time, &x_c, &y_c); begin_f_loop(f, thread) { real x[ND_ND]; real r, flux; F_CENTROID(x, f, thread); /* 计算面中心到热源中心的距离 */ r = sqrt((x[0] - x_c) * (x[0] - x_c) + (x[1] - y_c) * (x[1] - y_c)); flux = gaussian_heat_flux(r); F_PROFILE(f, thread, position) = flux; } end_f_loop(f, thread) }这段代码的关键点有三个。第一，CURRENT_TIME获取当前时间步的物理时间，Fluent 在瞬态求解时会自动更新。第二，F_CENTROID获取面网格中心坐标，用来算距离。第三，F_PROFILE把热流值赋给对应的面。注意单位，Fluent 默认是国际单位制，所以速度用 m/s，热流密度用 W/m²。
3.2 动网格参数配置
正弦摆动热源本身不需要动网格，因为热源是作用在固定面上的，只是热流分布随时间变化。但如果你后面要加熔池表面变形或者自由表面，就需要动网格了。这里先讲热源部分的设置。
在 Fluent 里，把热源作用的面单独命名，比如weld_top_surface，然后在 Boundary Conditions 里把该面的 Heat Flux 设为udf sinusoidal_heat_flux。瞬态求解设置里，时间步长建议先设小一点，比如 1e-4 s，总时间根据焊接长度和速度算。
# Fluent TUI 命令示例 define/boundary-conditions/heat-flux weld_top_surface # 选择 udf sinusoidal_heat_flux时间步长和网格尺寸要匹配。如果网格尺寸是 0.5 mm，焊接速度是 5 mm/s，那热源走过一个网格需要 0.1 s。时间步长如果设 1e-4 s，一个网格要 1000 步，计算量很大。可以先从 1e-3 s 开始试，看温度场是否稳定，再决定要不要加密。
3.3 材料参数与相变设置
AISI 1045 中碳钢的关键参数：密度 7850 kg/m³，比热 460 J/(kg·K)，导热系数 50 W/(m·K)，固相线温度 1700 K，液相线温度 1750 K。这些参数在 Fluent 的 Materials 里设置。
液相分数用 Fluent 自带的 Solidification/Melting 模型，设置固相线和液相线温度，Fluent 会自动计算液相分数。如果你要观察熔池形态，液相分数云图比温度云图更直观，因为液相分数直接对应熔池区域。
# 材料参数设置示例 密度: 7850 kg/m3 比热: 460 J/(kg-K) 导热系数: 50 W/(m-K) 固相线温度: 1700 K 液相线温度: 1750 K 熔化潜热: 270000 J/kg4. 验证请求与成功结果：温度场与熔池形态检查
4.1 热源轨迹验证
在跑完整模型之前，先做一个简单的验证：把 UDF 编译好，加载到 Fluent 里，然后只跑几个时间步，看热源中心位置是否按正弦轨迹移动。方法是在 UDF 里加一行打印语句，把每个时间步的 x_c 和 y_c 输出到控制台。
Message0("time=%f, x_c=%f, y_c=%f\n", time, x_c, y_c);如果输出显示 x_c 随时间线性增加，y_c 在正负幅度之间正弦变化，说明轨迹逻辑正确。这一步花不了几分钟，但能避免后面跑了几小时才发现热源没动。
4.2 温度场云图检查
跑完几个时间步后，在 Fluent 后处理里看温度云图。正常情况下，温度高值区应该跟随热源移动，并且因为横向摆动，热影响区不是一条直线，而是呈现周期性摆动轨迹。如果温度场一直固定在某个位置不动，说明 UDF 没生效，检查边界条件里是否正确绑定了 UDF。
一个常见的现象是温度场出现了但数值不对，比如最高温度只有几百度，远低于钢材熔点。这通常是热流密度单位或者数值量级搞错了。检查 Q_MAX 是否用了 W/m²，如果模型尺寸是 mm 级，热源半径也要换算成 m。
4.3 熔池液相分数与深度提取
液相分数云图能直接反映熔池区域。在 Fluent 里创建等值面，液相分数大于 0.5 的区域就是熔池。然后沿焊接方向取中间截面，观察熔池深度随时间的变化。
# 后处理提取熔池深度 1. 创建等值面 liquid_fraction = 0.5 2. 在中间截面 y=0 上提取该等值面的最大深度 3. 导出随时间变化的数据如果熔池深度曲线出现剧烈震荡，通常是时间步长太大或者网格太粗。先把时间步长减半，如果还是震荡，再加密熔池附近网格。
4.4 温度云图 MP4 导出
Fluent 后处理里可以导出瞬态动画。设置好温度云图的显示范围，然后在 Solution Animation 里创建动画序列，导出格式选 MP4。注意帧率不要设太高，否则文件很大，一般 10-15 fps 就够了。
# 动画导出设置 Solution Animation > Create Display Type: Contour Export Format: MP4 Frame Rate: 125. 本篇常见错排查
5.1 UDF 编译报错M_PI未定义
M_PI在标准 C 里不是默认定义的，有些编译器需要手动定义。解决办法是在 UDF 开头加一行：
#ifndef M_PI #define M_PI 3.14159265358979323846 #endif5.2 热源不移动或轨迹不对
先检查CURRENT_TIME是否返回了正确的时间值。如果 Fluent 是稳态求解，CURRENT_TIME始终为 0，热源就不会动。确保求解器设置为瞬态。另外检查DEFINE_PROFILE是否绑定到了正确的面，有时候面选错了，热源作用在了一个不动的面上。
5.3 温度场不收敛或发散
焊接熔池仿真发散的原因很多，最常见的是时间步长太大、网格太粗、热流密度太高。排查顺序：先把热流密度减半，看是否收敛；如果收敛了，说明热流密度量级有问题。如果还是发散，把时间步长减半。如果仍然发散，检查网格质量，尤其是热源附近的网格。
5.4 液相分数不更新
如果液相分数一直是 0 或者一直是 1，检查 Solidification/Melting 模型是否开启，固相线和液相线温度是否设置正确。另外，如果材料参数里没有设置熔化潜热，液相分数计算会不准确。
5.5 后处理动画导出失败
MP4 导出失败通常是编码器问题。Fluent 自带的编码器有时候和系统不兼容，可以先把动画导出为图片序列，然后用其他工具合成 MP4。或者检查导出路径是否有中文或空格，Fluent 对路径比较敏感。
5.6 TaoToken API 调用返回 401
如果 AI 辅助工具调用 TaoToken API 返回 401，先检查 API Key 是否正确复制，有没有多余空格。然后检查 base_url 是否填成了https://taotoken.net/api，注意不要加多余的路径。如果还是 401，到 TaoToken 控制台确认 Key 是否被禁用或者额度是否用完。
6. 仿真辅助工具的调用配置与后续扩展
6.1 用 TaoToken 统一管理仿真辅助工具
做焊接熔池仿真的时候，AI 辅助工具主要用在几个地方：生成 UDF 骨架、检查边界条件设置、整理后处理脚本、解释报错信息。这些工具如果各自用一套 Key，管理起来很麻烦。TaoToken 的统一通道可以让这些工具都指向同一个 API 地址，你只需要在 TaoToken 控制台维护一个 Key。
具体配置：在工具的配置文件里设置base_url = https://taotoken.net/api，api_key = 你的 TaoToken Key。模型名称按工具要求填，比如claude-sonnet-4-20250514。如果你要长期做仿真辅助，建议用 Coding Plan，因为调用频率高，按量计费容易超预算。
6.2 后续扩展方向
当前模型属于工程验证阶段，跑通了热源轨迹、温度场和液相分数。后续可以继续完善：加密熔池附近局部网格、引入温度相关材料参数、对比不同摆动幅度和频率下的熔池形态、加入 Marangoni 驱动流动、尝试 VOF 自由表面模型、输出更完整的熔宽熔深曲线。
尤其是直线焊接和正弦摆动焊接的对比，值得认真做。只有把两种热源放在一起对比，才能清楚看出摆动焊接对熔池宽度、熔深和温度分布的影响。
6.3 分阶段验证思路
焊接熔池仿真最怕一上来就把模型做得太复杂。如果第一步就加入 VOF、自由表面、表面张力、Marangoni 效应、相变、流动、复杂材料参数，很容易不收敛，最后连问题出在哪里都很难判断。
建议分三个阶段：阶段 A 只做热源轨迹 + 温度场 + 液相分数验证，确认模型能跑通、热源逻辑正确、结果趋势可信。阶段 B 加入熔池内部流动、浮力、Marangoni 效应。阶段 C 再加入保护气体区域、VOF 自由表面、表面张力随温度变化、熔池表面变形。每个阶段稳定后再进入下一阶段，效率会高很多。
如果你在配置 UDF 或者排查收敛问题时需要 AI 辅助，可以到 TaoToken 的模型对话里直接问，把报错信息和 UDF 代码贴进去，让它帮你定位问题。如果你要长期做仿真辅助，比如批量生成不同参数的 UDF、自动整理后处理数据，用 Coding Plan 会更合适。API Key 在 TaoToken 控制台的 API Keys 页面创建，接入文档在 doc 页面可以查到详细的 base_url 和模型列表。
