风扇气动噪声分析

仿真有限元分析优化设计代做加微信15527309811咨询报价，仅需100元起步

 

 

仿真有限元分析优化设计CAE代做咨询电话02759102656，

加微信15527309811咨询报价，仅需100元起步

 

风扇气动噪声分析服务利用仿真工具进行风扇/风机流动特性和噪声特性计算。风扇可以用于发动机的冷却等很多场景，合理的风扇设计将极大地提高风扇的效率，但由于管道风扇内部流动非常复杂，通过理论计算对其流动进行定性分析十分困难，风洞试验虽然可以得到其流动参数和噪声特性，但也无法对流场内部的流动细节进行描述。因而进行仿真计算是一种有效的设计计算方法。

 

管道风扇气动噪声分析：

1.1介绍

风扇可以用于发动机的冷却等很多场景，合理的风扇设计将极大地提高风扇的效率，但由于管道风扇内部流动非常复杂，通过理论计算对其流动进行定性分析十分困难，风洞试验虽然可以得到其流动参数和噪声特性，但也无法对流场内部的流动细节进行描述。

演示如何利用Fluent进行风扇流动特性和噪声特性计算。

1.2几何建模和流场计算域建立

风扇外径为384mm，轮毂直径为140mm，轮毂比为0.365，8扇叶均匀分布，外流场建模充分考虑到进气试验标准，入口区长度至少为入口处管道直径的六倍；而出口区的长度则应保证至少为出口位置管道直径的十倍；至于旋转流体区，是指包含了风扇本体以及周围流场的圆柱体区域，应当保证其尺寸尽量靠近风扇叶片的直径，最终风扇模型和外流场模型分别如下图所示。

 

 

 

 

 

 

 

1.3模型网格的划分

网格生成作为仿真计算中的关键环节，其结果直接控制了后续计算过程的效率与精度。为了保证划分结果的质量，应选择合适的网格尺寸，防止太疏或太密的网格产生，在流量梯度较大的流动区域内，应当尽量提高网格质量（高细密度，较小的歪斜度）；至于梯度小的区域可以在保证精度的基础上适当较少网格数目。

本案例旋转流体区由于包含了风扇本体且流动情况最为复杂，为了保证足够的计算精度，该区域网格尺寸最小。管道区网格尺寸较旋转区略大，最终划分结果如下图。

 

 

 

 

1.4边界条件设定与旋转模型选取

完成网格生成后需进行边界条件的设置。在流动的计算过程需要设定的边界条件包括：

（1） 流动入口条件：根据吸气试验的要求将流动入口设置为压力边界条件，其中入口处压值定义为大气压力，且气体沿轴线方向流动；

（2） 流动出口条件：根据吸气试验的要求将流动入口设置为压力边界条件，出口压力值定义为 0，即出口处没有外界的作用；

（3） 壁面边界条件：主要为通流区的管壁表面。

对风扇旋转运动的仿真则是通过 MRF 模型来实现的。Fluent 中常用的多运动坐标系模型包括： SMM（滑移网格模型）,MPM(混合面模型)以及 MRF（多重参考系模型）。考虑到风扇中气体运动属于定常流动，所以选择计算量相对较少的 MRF 基准。

作为旋转机械仿真中最常使用的模型，MRF 模型计算思路在于：将算法区间分成数个运动相互独立的子区间，先在各子部分间对流场方程进行求解，通过各部分间的交界面完成流场信息的传递。

作为 CFD 模型中唯一运动的旋转流体域，将其边界条件设定为 Fluid（流动域）。在 Fluent 软件中选择 MRF，并且将风机转速定义为坐标系的转速。

 

 

 

 

1.5计算方程选择与仿真参数设置

对于风扇内部的稳态流动，采用定常计算模型进行仿真，且计算过程中不考虑重力的影响。利用 SIMPLE 方程完成速度与压力的解耦，将湍流模型定义为 RNG k-epsilon;的双方程模型。

在 Fluent 软件对参数进行设定时，根据实际工况将流体材料定义为空气且认为风扇内部流体不可压缩；由于流动过程中没有热能的交换所以不对能量守恒方程进行求解，只考虑流体连续方程以及动量方程。对于控制方程中的湍动能耗散项以及动量项等使用二阶迎风的离散格式，在迭代过程中使用欠松驰因子以加速收敛。

1.6风扇流场计算结果分析

用Fluent软件对转速为2000rpm的风扇进行计算，得到包括速度矢量图、压力云图结果如下所示。

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2风扇气动噪声分析

 

2.1噪声分析步骤

在 Fluent 中对于风机噪声的仿真是分为两个部分先后完成的：

（1） 首先使用大涡模拟模型（LES）对风扇流场中的瞬态控制方程求解获得流场的动态稳定值，通过计算结果得到风扇的噪声源（即风扇叶片上的动态载荷）；

（2） 接下来则是通过求解 FW-H 模型的方法对风机载荷进行分析并得到噪声值。

2.2瞬态流场仿真边界条件设定

声场仿真过程中由于其 CFD 模型与流场极为相似因此不再另行建立模型，而是对原有流场模型的边界条件进行修改。由于噪声特性的仿真属于非定常计算，虽然同样将旋转流体域设为唯一的运动区域，但是改用滑移网格模型对风扇的动叶片与静止区域进行耦合以保证瞬态计算的精确度。

在控制方程的离散过程中使用PISO 算法代替原来的 SIMPLE 方法，相比较而言 PISO 算法在原有“预测-修正”方法的基础上添加了一个再修正过程，对原有计算结果进行了二次改进，有效的提高了计算精度与方程的收敛速度。至于迭代过程中参数的设置，将时间步长设定为0.0001s，而计算的截止频率取 6000Hz，在每个时间步长内计算 40 次，迭代次数为 1000。

2.3气动噪声边界条件设定与后处理

将2.2节中的计算结果与 FW-H 方程相结合在叶片表面使用二重积分获得随需要的声压值信号，在进行噪声参数设定时，以风扇本体为噪声源，而监测点则按照 GB/T2888-2008《风机和罗茨鼓风机噪声测量方法》中的规定设置，取风机前 1m 处噪声结果作为分析。

最终可以得到1m处噪声值计算结果为78dB(A)，完成声场计算过程后得到的数据是时域信号，还需使用Fluent 软件后处理功能中的快速傅立叶变换模块（FFT）完成时频转换获得声压级频谱图，如下所示。

 

 

 

 

 

 

 

风机气动噪声分析：

1模型建立

首先对模型进行几何处理，模型只包含两个流体域，一个为风机内部的旋转域，一个为风机外部的静止域，入口定义为inlet，出口定义为outlet，机器人定义为jiqiren，其余壁面定义为wall，如下图。

 

 

 

 

以meshing模式打开fluent对流场网格进行网格划分，应用Watertight Geometry流程划分网格，首先将第1节中建立的几何导入进来，之后将interface接触的交界面定义为200mm，接下来进行面网格划分，设置面网格最小尺寸为30mm，最大尺寸为1000mm，网格增长率为1.1，为了对曲面进行精确捕捉，设置Curvature和Proximity方法对曲率和微小结构进行识别，将inlet设置为速度入口边界条件，out设置为压力出口边界条件，其余边界均暂时设置为wall边界，添加边界层，第一层高度为0.1mm，共3层，增长率为1.1，最后对体网格进行划分，网格选择为四面体网格，最大网格尺寸为1000mm，最终网格划分结果如下图。

 

 

 

 

 

 

 

2仿真与后处理

入口速度分别设置为3m/s和12m/s，风机转速为24rpm。

仿真设置时，首先选择压力基稳态求解器，湍流模型选择为SST k-w模型，入口速度确定为3m/s和12m/s，外场壁面设置为sym边界，设置迭代次数为2000开始计算，计算收敛后，转变为瞬态，sst k-w模型，时间步长0.0001s，共计算50步进行瞬态计算，直到结果收敛即可。

3m/s结果分析

速度云图

 

 

 

 

涡云图

 

 

 

 

湍流云图

 

 

 

 

12m/s结果分析

速度云图

 

 

 

 

涡云图

 

 

 

 

湍流云图

 

 

 

 

3气动噪声模拟

瞬态计算完成后激活FWH气动噪声模型，噪声源设置为风机，共设置三个接受位置，坐标分别为（-1，0）、（0.5，0）和（1，0），共计算250步，计算过程的残差图如下。

 

 

 

 

4气动噪声结果分析

3m/s噪声结果分析

监测点1（-1，0）声压图：27.43db

 

 

 

 

监测点2（0.5，0）声压图：29.53db

 

 

 

 

监测点3（1，0）声压图：26.88db

 

 

 

 

12m/s噪声结果分析

监测点1声压图：67.42db

 

 

 

 

监测点2声压图：77.90db

 

 

 

 

监测点3声压图：66.81db

 

 

 

 

结论：

可见，风机声压值都不大，随着来流风速的增加，监测点1（-1，0）、监测点2（0.5，0）和监测点3的声压值均大幅增加。