直角三通加热结构的多物理场耦合数值分析

栾振辉 陈涛

摘 要：以盘式干燥机加热盘为研究对象，为了研究其内部流场、温度场及热应力分布情况，将加热盘简化为三通管结构，提出了流-固-热-力多耦合分析方法，建立有限元模型，联合有限元软件ABAQUS求解4种时间下（0s，30s，60s，90s）加热盘的流场、温度场以及热应力分布。结果表明模型内通入高温烟气，加热盘入口处温度最高，随着时间变化30s之后流体模型内高温烟气分布均匀，温度均匀;固体模型温度场随着时间变化，温度升高，热应力逐渐增大，90s时热应力最大，发生热应力集中现象。

关键词：加热盘;流固耦合;数值模拟;热应力

中图分类号：U173.2 文献标志码：A

文章编号：1672-1098（2021）02-0009-07

收稿日期：2020-11-10

作者简介：栾振辉（1961-），男，安徽凤台人，教授，博士，研究方向：流体传动与控制。

Numerical Analysis of Multiphysics Coupling of Right-angle Three-way Heating Structure

LUAN Zhenhui， CHEN Tao

（School of Mechanical Engineering， Anhui University of Science and Technology， Huainan Anhui232001， China）

Abstract：Taking the heating plate of the disc dryer as the research object， in order to study its internal flow field， temperature field， and thermal stress distribution， the heating plate is simplified to a three-way tube structure， and a fluid-solid-thermal-mechanical multi-coupling analysis method is proposed. The finite element model was established， and the finite element software ABAQUS was combined to solve the flow field， temperature and thermal stress distribution of the heating plate in 4 periods of time （0s， 30s， 60s， 90s）. The results show the followings. The temperature at the inlet of the heating plate is the highest when the high-temperature flue gas introduced into the model， and the high-temperature flue gas is well-distributed and the temperature is uniform in the fluid model after 30s The temperature field of the solid model changes with the time. The stress gradually increases with the increasing temperature， and the thermal stress is the largest after 90s with the thermal stress concentration.

Key words：heating plate;fluid-structure coupling;numerical simulation;thermal stress

盤式干燥机主要由多层固定空心加热圆形载料盘及转耙搅拌装置组成，是一种立式连续的以热传导为主的干燥设备[1]，如图1所示，其广泛应用于干燥、加热、脱水等领域[2-3]。加热盘是盘式干燥机的核心部件之一，其结构为一空心圆盘，分为大小两种，如图2所示（以小加热盘为例）。加热盘工作时受到高温烟气直接作用，导致其局部温度过高且产生较大热应力，极易发生热疲劳破坏。工程实践中，经常出现加热盘过早损坏的情况。为了提高加热盘的使用寿命，有必要对加热盘的受热情况进行分析。

多物理场耦合分析是一种综合考虑多种情况影响下的研究方法。文献[4]建立流固共轭传热模型，对钠泵机械密封摩擦副进行流固热耦合分析，对钠泵机械密封结构设计和改进具有指导意义;文献[5]建立热力耦合模型，仿真分析了制动器的制动过程，联立Manson-Coffin方程预测制动鼓的疲劳寿命;文献[6]建立热-流-固耦合模型，分析了不同工况下对单金属密封端性能的影响;文献[7]建立热力耦合模型，通过预测模型耦合求解，验证了预测模型结果符合实验结果，对工程实践有指导作用;文献[8]建立热力耦合模型，分析了不同工况下镐形截齿的温度分布，为进一步研究截齿应力提供了理论依据;文献[9]建立热力耦合模型，分析了应力分布、热变形和温度场情况，对高强度柴油机活塞结构优化具有指导作用;文献[10]建立流固耦合传热模型，分析了柴油机气缸盖的温度场分布，并通过实验验证了模型的正确性;文献[11]耦合多物理场建立流固耦合传热模型， 并用实验验证了可行性。 国外学者对流固耦合、 热力耦合研究较多， 文献[12]提出了三维弱耦合流体热机械系统的形状和拓扑优化框架; 文献[13]提出了将热边界条件纳入CPFEM模型的方法;文献[14]建立用有限体积法离散质量、动量和能量守恒的宏观控制方程， 模拟具有弯曲边界的流体-固体共轭传热问题; 文献[15]提出采用弱可压缩光滑粒子流体力学（WCSPH）方法求解流体-固体共轭传热（CHT）内流动问题; 文献[16]基于热耦合思想， 采用多目标优化方法对结构进行了优化。

基于国内外学者对流固耦合、热力耦合分析研究较多，而流固耦合研究变形固体在流场作用下的各种方式以及固体变形对流场的影响，热力耦合研究应力场与温度场之间相互作用的过程，通常对这两种方法单独进行研究，而综合研究较少。物理现象通常都不是单独存在的，在工程实践中，流场、温度场及应力场之间是相互作用、相互影响的，需要进行综合考虑。为此，提出一种多物理场耦合分析方法，综合考虑多种情况，相比较单个物理场分析或两个物理场分析更具有实用性。

因此，针对加热盘的传热过程及热应力，本文拟用流-固-热-力多耦合分析方法，探究其内部流场、温度场、热应力分布情况。鉴于加热盘传热过程比较复杂，本研究拟将加热盘简化为三通结构进行研究，通过建模及仿真分析研究其内部情况，最后通过实物模型进行实验室验证。

1 理论模型

1.1 流体域方程

加热盘内高温烟气流动遵循质量守恒定律、动量守恒定律、能量守恒定律[17]。

ρft+xjρfvi=0（1）

ρfvit+xiρfvivj=-pxi+xiμ·vixi（2）

ρEt+SymbolQC@·vρE+p=SymbolQC@·keffSymbolQC@T-∑jhjJj+τeff·v+Sh（3）

式中：v为流体速度;ρf为流体密度;下标i，j=1，2，3，表示3个方向的分量，t为时间;p为流体压力;μ为动力粘度;ρ为流体的密度，E为流体微团的总能，τeff为有效应力，hj为组分j的焓，keff为有效导热系数，Jj为组分j的扩散通量，Sh为包括化学反应热及其他体积热的源项。

1.2 固体域方程

Msd··+Csd·+Ksd+τs=0（4）

式中：Ms为固体单元质量矩阵，d为固体单元的位移矢量，Cs为固体单元阻尼矩阵，Ks为固体单元刚度矩阵，τs为固体单元受到的应力。

固体区域增加了由温差引起的热变形为

fT=αT·SymbolQC@T（5）

式中：αT为与温度相关的热膨胀系数。

1.3 流-固耦合界面方程

流体和固体区域采用直接流-固耦合方法，实现固体传热和流体传热的耦合，得到所需要流体流场以及流体和固体温度场，满足能量连续性条件温度和热流密度相等[18]。控制方程为

Tf=Ts（6）

qf=-λfTfn=-λsTsn=qs（7）

式中：Tf为流体温度，λf为流体导热系数，Ts为固体温度，λs为固体导热系数，qf为流体侧热流密度，qs为固体侧热流密度，n为为流-固交界面法向量。

2 有限元建模

为简化计算，将图2所示的加热盘的进口端理想化为一个直角三通模型，一个入口和两个出口的截面尺寸均为40mm×40mm，入口长度为200mm，两个出口相距400mm，壁厚为6mm。高温烟气由入口端流入，两个出口端流出。通过FSI（Fluid-Structure-Interaction）耦合ABAQUS软件的传热模块（Heat transfer）和流体模块（CFD），模拟高温烟气在流动过程中与三通模型之间的传热过程。

首先，采用DC3D8单元对直角三通实体模型划分网格，经试算确定网格尺寸为10mm，为有效表征三通模型壁厚方向的温度梯度，壁厚划分三层单元，如图3所示。直角三通模型的内表面定义为流固耦合界面，外表面施加表面膜约束实现对流换热，空气的自然对流换热系数为10W/（m2·K），环境温度为20℃。直角三通实体模型的材料为不锈鋼，相关热力学参数如表1所示。

其次，基于直角三通实体模型，建立三维流场模型如图3所示。采用FC3DB单元对其进行离散，单元网格尺寸与直角三通实体模型相一致，即为10mm。流场模型的入口流速为3.5m/s，入口温度为800℃，出口表压为0MPa。流场模型的外表面（除去一个入口和两个出口）定义为流固耦合界面。流场模型通过这个流固耦合界面实现烟气与直角三通模型内表面之间的传热，烟气、模型的相关力学参数如表1所示。

最后，将直角三通实体模型计算的瞬态温度场导入ABAQUS/Standard模块，采用顺序热应力耦合分析技术计算直角三通实体模型响应的热应力。在该静力学计算模块，直角三通实体模型由C3D8R单元进行离散，网格尺寸与传热模型一样。固定直角三通实体模型的底面，以限制其刚体位移，详细的计算流程如图4所示。

3 结果分析

根据高温烟气的进口流速， 计算其雷诺数为1 201.5，故可判断为层流状态。图5给出了直角三通内部烟气流动90s时刻的速度场分布，显然可见烟气的最大流速位于进口直边段的中心位置，两端出口直边段的流速十分均匀。图6展示了烟气在流动0s、30s、60s和90s时刻的温度分布。从图6中可以看到，在0s时刻，高温烟气开始进口端向内部流动，30s后高温烟气的温度场基本保持不变，这是因为烟气从进口端流到出口端只要0.12s（t=（0.22+0.2）/3.5），30s后烟气自身的温度分布基本保持稳定了。

图7所示为高温烟气与直角三通接触面沿路径的温度分布情况。由图7可见，高温烟气外表面（流固耦合界面）上的温度整体上沿其流动方向逐渐减小（即进口端温度最高，出口端温度最低），但在三通内部烟气流速变向的位置温度产生的波动。随着时间的增加，烟气同一位置上的温度逐渐升高，说明在流固传热过程中，三通实体模型的温度逐渐升高，如图8所示。在三通模型内部的流固耦合界面上选取三个位置，分别记为A、B和C，其中A点位于进口端，B点位于拐角处，C点位于出口端。A、B和C三个位置处的温度随时间的变化如图9所示，显然可见，三个位置处的温度都随时间逐渐增大，但A点温度升高的速率最大，B点次之，C点最小，三点位置处的温度变化与图8相一致。

将图8中的温度场导入ABAQUS静力学模块中，进行顺序热力耦合计算，得到三通实体模型的热应力分布如图10所示。显然可见，与图8相一致，随着烟气流动时间的增加，直角三通实体模型响应的热应力逐渐增大。由于直角三通的底面固定，限制了温度变化引起的变形，因此直角三通底部的热应力较大。在直角三通底部建立一条路径，沿该路径分布的热应力如图11所示。由图11可见，随着高温烟气流动时间的增加，路径上相同位置响应的热应力逐渐增大，与图10相一致;对于相同的高温烟气流动时间，沿图11（a）分布的热应力呈现近似对称分布，进口端和出口端的热应力较小，拐角处的热应力最大，这与实际加热盘隔板处易产生裂纹的现象相吻合，其主要可归结两方面原因：（1）进口端和出口端有一面是自由的，对温度变化引起的热变形的约束较弱;（2）由图7（b）可见拐角处的温度产生波动，产生较大热变形，进而诱导较大的热应力。

4 结论

（1）高温烟气外表面上的温度整体上沿其流动方向逐渐减小，但在三通内部，随着时间的增加，烟气同一位置上的温度逐渐升高。

（2）随着高温烟气流动时间的增加，三通实体模型响应的热应力逐渐增大。

上述结论与实际加热盘过早失效的现象相一致，这为下一步加热盘的结构优化奠定了基础。

参考文献：

[1] 刘广文. 盘式干燥机[J].农药， 2008（5）：388-389.

[2] 贾海华. 盘式干燥机原理与脱水动力学的研究[D].北京：煤炭科学研究总院，2003.

[3] 刘小龙，蒋家羚，刘宝庆. 盘式连续干燥器的传质传热模型[J].化工学报，2004（7）：1 047-1 052.

[4] 谷继品，陈侃，刘伟，等.一回路主循环钠泵机械密封热流耦合分析[J].液压气动与密封，2020，40（7）：36-42.

[5] 于振华，计三有，何奎.单向自增力式制动器热力耦合分析及疲劳寿命预测[J].机械传动，2020，44（10）：141-147.

[6] 马艺，倪洋，孟祥铠，等.牙轮钻头用单金属密封热-流-固耦合模型与性能分析[J].中国机械工程，2020，31（19）：2 295-2 303.

[7] 余建航，屈治国，李澎，等.逆流湿式冷却塔宏观预测耦合模型及其热力性能研究[J].热能动力工程，2020，35（6）：183-190.

[8] 庞国辉，王义亮.镐形截齿截岩过程的温度场研究[J].煤矿机械，2020，41（6）：29-32.

[9] 蔡翌辉，刘畅，郭豪杰，等.高强度柴油机活塞温度与应力场有限元分析及其结构改进[J].内燃机工程，2020，41（1）：1-11.

[10] 魏丹，庞华庭，段敏，等.柴油机缸盖流固耦合传热分析与实验验证[J].机械设计与制造，2019（12）：67-70.

[11] 高宏亮，李建统，尚芳，等.小型风冷汽油机流固耦合传热模拟计算研究[J].小型内燃机与车辆技术，2020，49（1）：54-62.

[12] F FEPPON，G ALLAIRE，C DAPOGNY，et al.Topology optimization of thermal fluid-structure systems using body-fitted meshes and parallel computing[J].Journal of Computational Physics，2020，417：156-161.

[13] DANIELS Connolly，CHRISTOPHER P KOHAR，WAQAS MUHAMMAD，et al.A coupled thermomechanical crystal plasticity model applied to Quenched and Partitioned steel[J].International Journal of Plasticity，2020， 133：274-279.

[14] YANG L M，SHU C，CHEN Z，et al.Three-dimensional lattice Boltzmann flux solver for simulation of fluid-solid conjugate heat transfer problems with curved boundary[J].Physical review E， 2020， 101（5-1）：053309.

[15] K C Ng，Y L Ng，T W H SHEU，et al.Fluid-Solid Conjugate Heat Transfer Modelling using Weakly Compressible Smoothed Particle Hydrodynamics[J].International Journal of Mechanical Sciences， 2018：362-365.

[16] ROOZBEH GHANI， DAVOOD IRANSHAHI. Conceptual comparison of four configurations in the thermal coupling of ammonia synthesis and 2-butanol dehydrogenation[J].Applied Thermal Engineering， 2019， 154：238-250.

[17] 陈汇龙，李同，赵斌娟，等.动压型机械密封内流场及性能的流固热耦合[J].排灌机械工程学报，2017，35（6）：502-507.

[18] 杨岑玉，胡晓，李雅文，等.流固耦合固体相变蓄热建模與分析方法[J].储能科学与技术，2019，8（3）：538-543.

（责任编辑：李 丽）