气动充放气系统二维非定常流场数值模拟

气动充放气系统二维非定常流场数值模拟朱冬;杨庆俊;包钢【摘要】以典型气动系统为研究对象，考虑系统内部气体、元件和环境空气之间的耦合传热，将系统内部流场、元件和外部空气作为整体进行研究，建立气动充放气系统的二维模型，采用湍流κ-ε两方程模型，利用有限体积法(FVM)和动网格技术对充放气过程进行计算，得到系统内气体状态参数。

在充放气过程中，气体速度最大值出现在放气管出口处，可达到超声速范围；气缸中温度梯度主要存在于内壁面处。

在系统元件中，进气节流口壁面温度的上升幅度最大，放气节流口壁面及紧邻的气管壁面温度降幅最大。

通过试验验证本方法的可行性，为非定常流场的计算和分析提供参考。

%Considering the coupling heat transfer of pneumatic components, air inside the system and ambient air, the flow field inside system, components and air outside system of typical pneumatic system treated as research object were investigated as an entirety. And the two-dimension model of pneumatic charging and discharging system was putted forward. The unsteady flow field was calculated by finite volume method (FVM) and dynamic mesh with the turbulence k-ε model. The results show that the air's maximal velocity appears at the outlet of the discharging pipe, which is faster than speed of sound in the charging and discharging process, and the temperature gradient inside cylinder exists nearby the cylinder wall. Among the components, the temperature of charging speed controller rises most, and the temperature of discharging speed controller and the neighboring pipe reduces most. The calculation results are consistent with those of experiment, which validates theproposed method is feasible. It illustrates a new way for calculation and analysis of unsteady flow field.【期刊名称】《中南大学学报（自然科学版）》【年(卷),期】2011(042)006【总页数】6页(P1663-1668)【关键词】气动系统;充放气;非定常流;有限体积法【作者】朱冬;杨庆俊;包钢【作者单位】哈尔滨工业大学气动技术中心,黑龙江哈尔滨,150001;哈尔滨工业大学气动技术中心,黑龙江哈尔滨,150001;哈尔滨工业大学气动技术中心,黑龙江哈尔滨,150001【正文语种】中文【中图分类】TH138随着气动系统在现代工业中越来越广泛的应用，人们对气动系统运行过程中压缩空气的状态变化越来越关注[1-3]。

系统内气体状态包括层流和湍流等多种形式，且气体与系统元件、元件与外部空气之间存在着非稳态传热[4]。

同时，许多气动元件的结构复杂，其内部流场的测量变得困难有时甚至不可能实现[5]，因此，对气动流场进行数值计算就显得尤为必要。

许多学者对气动系统流场进行了研究，如：李军等[6]采用有限差分法对气动充放气系统进行了一维流场计算，对气缸或容腔中的气体采用集中参数法，将节流口进行等熵简化处理；李玉军等[7]将压模型和压力回复模型引入到一维流场中用来处理节流口边界条件，计算得到的流场不能反映气管、气缸等元件的径向参数分布以及节流口对管路流场的节流特性，难以满足实际需要。

杨丽红[8]对固定容积容器的放气过程进行了二维仿真研究，计算时考虑容器壁和内部气体的传热，得到了放气时放气口的流场分布。

本文作者建立了气动充放气系统的二维模型，考虑系统内气体、元件和外部气体之间的热传递，采用湍流k-ε两方程模型，运用有限体积法(FVM)和动网格技术[9-11]对充放气过程进行了计算。

将计算结果与试验结果和绝热模型计算结果进行对比，表明本方法与试验结果具有较好的一致性，证明了本方法的正确性。

这一方法可以用于计算出气缸内部和节流口处的流场，为研究节流口节流特性及系统流场参数分布提供详细的数据和参考。

1 气动充放气系统模型的建立图1所示为典型气动回路。

气动回路包括气源、电磁阀、管路、节流阀、气缸和负载等。

气缸在换向阀的控制下进行充放气。

通过进气回路向气缸进气腔充气，通过放气回路从气缸放气腔向系统外放气。

活塞未运动时，属于固定容积的充放气；活塞运动时，进气腔和排气腔的容积随时间不断发生变化，此时属于变容积充放气。

充气时，进气腔气体受到压缩，温度升高，进气腔侧元件包括气缸壁、节流阀和气管等和内部气体进行热交换而温度升高；元件外部空气温度也由于元件温度的升高而升高。

放气过程中，放气腔气体由于膨胀而温度降低，放气腔侧元件和内部气体进行热交换而温度降低；元件外部空气温度也由于元件温度的降低而降低。

图1 典型气动回路Fig.1 Typical pneumatic circuit1.1 数学模型的建立气体在充放气系统中的流动属于可压缩湍流。

二维控制方程如下。

(1) 连续性方程[12]：式中：ρ为空气的密度；x为轴向坐标；r为径向坐标；t为时间；ux为轴向速度；ur为径向速度。

(2) 轴向和径向的动量守恒方程：(3) 能量方程：式中：T为气体热力学温度；λ为导热系数；cp为气体的定压比热容；ST为广义源项。

(4) 状态方程：式中：R为摩尔气体常数。

(5) 活塞运动方程：式中：m为活塞和负载质量；p1为气缸无杆腔绝对压力；p2为气缸有杆腔绝对压力；pa为大气绝对压力；A1为气缸无杆腔侧活塞面积；A2为气缸有杆腔侧活塞面积；f为摩擦力；F1为负载力。

(6) k-ε两方程湍流模型[13]：其中：为湍流脉动动能；ε为脉动动能耗散率；µt为湍流黏度系数。

方程中常数取值为：c1=1.44，c2=1.92，cµ=0.09，σk=1.0，σε=1.3。

1.2 几何模型及网格划分图1所示的气动系统具有复杂的三维结构。

为了减少计算量，本文将气缸简化为圆柱体，将节流阀简化为节流口，气缸的进气口和出气口分别挪移到气缸两底面中心，如图2所示。

模型中的气源区域与进气管相连，为系统提供恒压和恒温的气体。

放气管出口与大气区域相连。

图2 简化的气动回路模型Fig.2 Simple model of pneumatic circuit简化后的气动回路进一步简化为以水平中心线为对称轴的二维模型。

利用 Fluent软件的前处理软件Gambit建立二维模型，如图2所示。

采用四边形网格和三角形网格对模型进行分区域网格化。

气管和节流口等处采用较密的三角形网格，气管壁等处采用四边形网格。

模型网格如图3所示(由于尺寸限制只显示进气节流口处网格)。

图3 模型网格Fig.3 Mesh of model2 初始条件和边界条件的设定初始时模型中气源区、气缸放气腔、放气节流口与放气管压力均为供气压力，其他区域为大气压力。

边界条件包括进气口、放气口和壁面等。

进气口为压力入口边界，压力为供气压力，温度为供气温度；放气口为压力出口边界，压力为大气压，温度为大气环境温度。

给定各壁面温度，包括气管壁面、节流阀壁面和气缸壁面等。

活塞的壁面设为对流传热。

管壁、节流口壁和气缸壁处为耦合壁面，计算内部气体、壁面和外部气体之间的热传递。

近壁面区采用标准壁面函数。

3 计算结果与分析3.1 计算结果与试验对比为了验证本方法的正确性，依照表1中参数建立模型，利用Fluent软件进行计算。

气缸无杆腔为进气腔，有杆腔为放气腔。

用 PISO法对压力场和速度场进行求解。

运动方程通过自定义函数实现，实时读取活塞两侧的压力，根据式(6)计算得到活塞的速度并返回到模型中。

活塞以相应速度运动使得进气腔逐渐增大而排气腔逐渐减小。

利用有限差分法建立相同参数的绝热模型进行计算，同时，按照图1和表1建立气动回路进行充放气试验，测量气缸两腔压力。

进气腔压力如图4所示，放气腔压力如图5所示。

从图4和5可以看出：与绝热模型的计算结果相比，本方法的计算曲线与试验曲线具有更好的一致性，表明了本方法的可行性和正确性。

整个充放气过程约为0.55 s。

初始时活塞静止，充气腔和放气腔分别为定容充气和定容放气。

在活塞伸出过程中，充气腔压力先升高后又降低，放气腔压力单调降低；充放气0.33 s时，活塞杆完全伸出，活塞速度瞬时突变为零，充气腔压力逐渐增大为气源压力；放气腔压力迅速减小最终稳定为大气压。

表1 元件参数Table 1 Parameters of components序号元件或参数型号或数值1 气管 TU0805，管内径5 mm，外径8 mm，进气管长3.2 m，放气管长2 m 2 节流阀 AS3000，进气口、放气口节流面积均为12.3 mm2 3 气缸CDM2B40-300A 4 负载/N 9.1 kg，水平加载5 压力传感器 PSE510-R06，量程0～1 MPa 6 供气压力/MPa 0.6 7 供气温度/K环境温度/K 293 293图4 气缸进气腔压力曲线Fig.4 Pressure curves of charging chamber in cylinder图5 气缸放气腔压力曲线Fig.5 Pressure curves of discharging chamber in cylinder3.2 计算结果分析图 6所示为充放气过程中放气管轴线的压力分布。

放气管长2 m，与放气节流口连接处为起点，与大气区域相连处即放气管出口为末端。

气体经节流口进入放气管，由于射流的作用压力有明显下降。