基于Fluent薄钢板对焊过程数值模拟研究.docVIP

基于Fluent薄钢板对焊过程数值模拟研究 　　摘要：金属平板的焊接过程是一个非常复杂的非稳态传热过程。对该物理过程进行简化的基础上，建立了数学模型。利用Fluent软件进行了数字模拟，数据结果显示了金属平板内温度场的分布，弥补了实验手段的不足。 　　关键词：钢板焊接；温度场；Fluent软件；数值模拟 　　中图分类号：TG174.44文献标识码：A文章编号：1006-8937（2011）16-0091-02 　　 　　焊接过程的热传导是一个非常复杂的非稳态过程。傅里叶定律与能量守恒定律导出的导热微分方程和定解条件构成了焊接导热问题的完整的数学描写。但由于实际问题各种各样，边界条件十分复杂，利用解析的方法很难求解这类微分方程。随着计算机内存和并行技术的发展，通过计算机获得导热问题的数值解的方法迅速发展。对物理问题进行数值求解的基本思路是把原来在时间，空间坐标系中连续的物理量的场，用有限个离散点上的值的集合来代替，通过求解一定方法建立起来的关于这些值代数方程，来获得离散点数上被求物理量的值。计算流体力学（Computational fluid Dynamics，CFD）就是其中一种有效的研究传热问题的数值模拟方法。 　　1Fluent软件介绍 　　Fluent 软件是目前市场上最流行的CFD 软件。利用Fluent软件可以计算平板焊接过程的温度场和揭示传热规律。只要涉及流体，热传递及化学反应问题等工程问题，都可以用Fluent软件求解。Fluent 程序包括以下几个部分：Gambit用于建立几何结构和网格的生成；Fluent用于传热流动等模拟计算的求解器；prePDF用于模拟PDF 燃烧过程；TGrid用于从现有的边界网格生成体网格。 　　Fluent 求解的思路和步骤是首先根据实际情况抽象和建立物理模型,然后选择求解模型,设置边界条件和初始化条件,设置迭代和控制参数,进行计算,经后处理和分析得出结论,或者根据实验结果重新调整参数进行再计算直至达到可接受的结果。 　　2两块薄钢板焊接过程热传导的数学模型 　　为对两块平板的对接焊过程（图1）进行数值计算，对其物理过程做以下简化处理： 　　钢板中的温度仅是x及时间τ的函数。 　　焊枪的热源作用在钢板上时钢板吸收的热流密度q（x）=qmexp（-3x2/re2）。rm为电弧有效加热半径；qm为最大热流密度。 　　钢板上下表面的散热可用q=h（t-tf）计算，侧壁绝热。 　　焊接热源为非熔化极热源，并假设由焊丝填充的焊道部分在焊接前就已存在了。由于此部分仅占整个焊件极小部分，故对整个导热过程影响不大。 　　钢板的物性为常数，周围环境温度不变，钢板的对流换热系数不变。 　　因此可以把图1所示的问题简化成图2所示的一维非稳态导热问题。 　　3FLUENT 的分析及数值模拟结果 　　①几何模型及网格划分。文章对长度H为12 cm,厚度δ为1 cm的两块相同的钢板对接焊过程进行模拟。由对称可知只需对一块钢板进行模拟。利用Fluent 自带的前处理软件Gambit来进行建模和网格划分。为了更好分辨热源附近的温度场采用了非均分网格,钢板的计算网格如图3 所示。 　　②计算参数的设置。将在Gambit中画好的钢板网格导入Fluent软件中，选择二维单精度压力基、隐式非稳态求解器，激活Energy Equation 选项。输入材料的物性参数为：密度7 800 kg/m3 ，导热系数为41.9 m2/s,比热容为670 J/（kg?K）,对流换热系数为12.6 W/(m2?K).电弧有效加热半径为0.71 cm,最大热流密度为1.68×108 W/m2。由于内热源是温度和距离的函数，故使用用户自定义函数（User-Defined Function）在文本编辑器中用C语言编辑内热源，热流密度作用时间为1s。钢板四周边界条为绝热。设置残差为1×10-6,钢板和环境初始温度为298 K。时间步长的选择不能太大,也不宜太小，文章时间步长为0.01 s。 　　③模拟结果及分析。图4为5 s时钢板内的温度分布图。从图中可以看出在高温热源的作用下，开始时时钢板的电弧加热斑点处温度可高达2 700 K左右,随着离电弧加热斑点距离的增加，温度急剧下降，在电弧有效加热半径之外的钢板的温度基本不变。由于向长度方向导热，并在对流换热作用下，5 s时电弧加热斑点处钢板温度已降到1 350 K左右，钢板升温的区域。 　　4结语 　　钢板的焊接过程是一个非常复杂的非稳态过程，必须抓住主要影响因素，忽略次要因素，对薄钢板焊接的物理过程进行了简化，建立了一维非稳态数学模型，利用Fluent软件进行了数值模型，得到的结果与文献的结果基本相吻合，因此模拟结果是正确的，利用Fluent软件分析钢板