太阳集团通风降温系统
电 话：0579-81328720>
传 真：0579-81328720
联系人
售前咨询:13388660553
技术指导:18858318765
售后服务:15068216608
地址：上海 金华 嘉兴 襄阳
厂房降温风机轴流通风机叶片模态仿真及其对气动噪声的影响后向式
摘要：利用有限元模型对叶轮模态进行了计算,判断分析了各阶模态振型对气动噪声的影响程度 , 求解中利用了ANSYS 的模态循环对称功能 , 同时分析了旋转软化、应力强化对叶轮真实运转状况下模态频率的影响。
Modal Simulation of Axial Fan Blade and The Effect of That on Aerodynamic Noise
Abstract:Calculation is carried out on impeller modal using the finite element model.Ｔhe effect of modal vibration in various stage on aerodynamic noise is judged and analyzed ,the effect of rotating softness and stress intensification on the modal frequency under the impeller in real operation condition is also analyzed at the same time.
一、引言
轴流通风机当其叶片较薄以及过度前掠,重心偏离叶根截面中心时,较高转速造成的离心力和不稳定进气流造成的叶片升力的变化,很容易激发叶片振动。同时由于流固耦合,还可能造成叶片的驰振 ,使叶片提前疲劳损坏,降低风机效率,并产生较大的气动噪声。
在叶轮设计时有必要对其振动模态进行计算,但叶片叶身曲面复杂,用经典理论无法求解 ,因此必须借用有限元模型来计算。ANSYS 是当今比较有名的有限元分析软件之一,具有多种物理场的求解功能 , 可以很方便地进行模态分析;大型 CAD 工程软件 UniGraphics 具有丰富的曲面造型功能 , 非常适合于叶轮等具有复杂曲面实体的造型,建好的实体模型导入ANSYS 即可进行模态分析。
二、叶轮 CAD 模型建立和接口导入
1. 叶轮基本参数
轴流通风机为整体注塑 ABS 塑料叶轮, 叶片数为4,叶片较宽,叶片呈前掠状。工作转速为 860 r/min,轮毂直径为0.147m,叶轮外径为0.42m 。
2. 几何模型建立
通过三坐标测量仪测量得到叶片表面型值点 , 将点阵连接成曲面 , 并利用软件 U G 的曲面剪裁和缝合功能 , 将叶片的曲面连接起来。一旦所有曲面被缝合就自动生成以各曲面为边界的实体。
叶轮为循环对称结构 , 为加快有限元分析过程 , 利用 ANSYS 的循环对称分析功能 , 对一个90 °基本扇区进行求解。建模时使全局坐标系的Z 轴与叶轮旋转轴线对应 , 建立完整叶轮模型,然后用过轮毂轴线两个相互夹角为 90 °的两个平面切出 1/4 的叶轮模型 (图 1) 。
3. 导入几何模型
能够将UG模型导入ANSYS 的方法有 3 种,其中基于直接的模型数据交换的两种是:一是通过标准的数据接口将 CAD 模型数据转入分析工程;另外是通过ANSYS为UG提供的专用接口直接读入UG的prt 文件;第三种借助UG的GFEMFEA。
这里采取第二种方法 , 在功能菜单中点击File → Import → U G, 再选取零件文件即可。
三、预处理和求解
1. 输入材料物理参数
输入ABS材料的物理性能参数:密度为1.2×10-6 g/mm3 ,弹性模量为2.3MPa,泊松比为0.38。
2. 选择单元类型
叶轮表面为变厚度复杂曲面 , 选用 10 节点的四面体单元 solid92 , 该单元采用二次位移模式 , 非常适合对形状不规则的实体划分有限元模型。
为了对基本扇区的两个间隔相对 90 °的轮毂的剖面划分网格 , 还选择了一种二维单元:MESH200 单元 , 并设定单元形状参数为“ triangle with 6 nodes ” (MESH200 单元是专门用来划分网格 , 提供网格占位功能 , 不参加单元运算 ) 。
3. 划分网格
先用 MESH200 三角形平面单元划分轮毂上的两个剖面的一个面上的网格 , 然后通过MSHCOPY 命令将该面上的网格拷贝到另外一个剖面上 ( 完成后的网格单元如图 2) 。 对整个模型用 solid92 单元分网格(完成后的网格单元如图 3) 。
4. 边界条件加载
叶轮通过轮毂的轴心线装配在电机轴线上 , 叶轮除转动外 , 其它运动都被限制。所以使轮毂圆柱装配面的有限元节点 X , Z 方向的自由度得到约束 ( 在全局柱面坐标系中 ) 。
有限单元的节点都有一个坐标系与载荷方向对应 , 在通常情况下节点坐标系与全局迪卡尔坐标系对应 , 必须首先用 NROTAT 命令转换节点坐标到柱面坐标系中 , 然后再在节点上加载位移约束。
5. 循环对称处理
循环对称模态求解是 ANSYS 对循环对称结构提供的一种特殊简化模态求解方法 , 在求解前有一些特殊的预处理。
首先 , 需要选择叶轮上下两个剖面上的节点并建立两个组集 , 取名为“Low”和“High”。其次运行 CYCGEN 的宏在基本扇区上建立第二个扇区 , 模态分析就是通过这两个扇区完成的 , 如果不带参数运行这个命令 , 它将内部耦合和约束方程也拷贝到第二个扇区上 ; 如果运行 CYCGEN ‘ , LOAD ’的命令 , 则会把负载也拷贝到第二个扇区上。这里运行 CYCGEN,‘LOAD’。
6. 求解
选用 Block Lanczos 为求解方法 , 设置求解频率范围为 20Hz 到 200Hz 。对于求解循环对称模态 ,ANSYS 也提供了专用的求解宏指令(不可直接用 solve 命令),该指令格式为：
CYCSOL , NDMIN , NDMAX , NSECTOR , LOW
各参数含义如下 :
NDMIN、NDMAX: 计算的上下节径范围, NDMIN 最小为0 ,NDMAX 对偶数最大可取 n/2 , 对奇数最大可取 ( n - 1) / 2 。
NSECTOR : 循环对称的扇区数,这里为4。LOW: 较低角度剖面上节点构成的组集名称。
该命令对应菜单路径 :Main Menu > Solution > Modal Cyclic Sym。
输入CYCSOL ,0 ,2 ,4 ,LOW进行求解。
四、仿真结果与试验结果的对照
试验是在叶轮安装状况下通过锤击法进行的 ,支架刚度很大,忽略支架的影响,认为结构在 20～200Hz 范围内的振动模态频率由叶轮决定。试验得到的叶轮模态频率值为 58.17Hz,83. 38Hz ,88. 69Hz ,154. 8Hz ;仿真得到的模态频率值约为 62Hz、80Hz、88Hz 和152. 2Hz。试验检测得出的频率与仿真结果对应关系较好,因此相互得到了验证。由于振型试验比较复杂 ,所以没有进一步作振型试验,后面将利用仿真的数据来观察振型。
五、考虑预应力和旋转软化
在真实状况下叶轮是运动的 ,由于离心力和气动载荷的影响 ,叶轮产生拉伸变形,模态有可能与静止状况有很大不同 ,所以必须予以考虑。
影响旋转件频率变化的一种原因是由于离心力对叶片运动产生的预应力的影响 ,造成了叶轮刚度的增大 ,使运行状况下模态频率升高。
另一种原因 :旋转软化,旋转软化使模态频率降低。其原理可以用一个简单的弹簧 - 质量旋转工程说明 (如图4) ,弹簧垂直于旋转轴,当弹簧刚度很高而旋转加速度很小时 ,认为弹簧变形很小。
忽略弹簧变形对质量块向心加速度的影响 ,建立如下平衡方程:
kx = Mωs２ｒ （１）
式中 k ———弹簧刚度
x ———离开平衡位置的距离
ωs ———旋转角速度
r ———质点自由位置相对于转轴的半径
但是如果弹簧刚度不够 ,同时旋转速度又很大 ,由于离心力的影响使弹簧产生较大位移,而该位移同时又使质点离心运动的半径加大,这时的平衡方程写为：
kx=Mωs２(r + x) ( 2 )
如果仍然用 (1) 式的形式表示的话,其平衡方程可以写为：
(k-Mωs２)x=Mωs２r
施加如图 4 所表示载荷时,其振动方程可写为：
Mx-(k-Mωs２)x=f(t)
因此刚度由 k 变为(k-Mωs２),即相当于旋转软化作用 ,旋转速度越高,旋转物体密度越大 ,这种软化作用也就越明显。 应力刚化使模态频率升高 ,旋转软化使模态频率偏低 ,通常应力刚化的作用偏大,所以同时考虑两种因素影响 ,使运转状况下模态频率比静止状况下模态频率偏高。
为了获得真实状况与静止状况下模态的差别 ,又进行了一次模态有限元分析,步骤是在三、四步骤后给叶轮施加一个转动角速度,打开预应力开关 ,选择分析类型为静应力分析,并进行一次静应力分析。然后选择分析类型为模态分析 ,并保证预应力开关为打开状态,同时打开旋转软化选项 ,下面同三、五以后的步骤。
计算结果各振型对应模态频率变化不到1Hz ,因此该叶轮可以采用静止状况下的模态来代替运动状态下的模态 ,从考虑问题的周全性上考虑 ,对应力刚化、旋转软化验证是必要的。
六、振型和对气动噪声影响的分析
为了观察振型 ,用命令Expand 并输入参数4 扩展成整个叶轮以观察振型(菜单路径:Main Menu > General Postprocessing > Expand sector) 。
(1) 一阶振动频率为61.5Hz ,叶片表现为沿径向扭摆 ,变形最大位置在叶根处,整个叶轮振型表现为 1、3叶片摆动时,2、4不动,1、3叶片反向扭摆 (如图5a) 。
(2)二阶振动频率为62Hz,叶片表现为沿径向线扭摆,变形最大位置在叶根处,整个叶轮振型表现为1 、3叶片向前扭摆,2、4叶片向后扭摆(如图5b) 。
(3) 三阶振动频率为62.5Hz ,叶片表现为沿径向线扭摆,变形最大位置在叶根处,整个叶轮振型表现为 4 个叶片以相同形式同向扭摆(如图5c) 。
(4) 四阶振动频率为80.3Hz ,叶片仍然表现为沿径向线扭摆 ,但叶片型面上有弯曲现象出现 ,弯曲最大现象出现在较大叶片半径处,整个叶轮振型表现为1、3叶片向前扭摆,2、4叶片向后扭摆(如图5d)。
(5) 五阶振动频率为80.5Hz ,叶片表现为沿径向线扭摆 ,叶片型面上有弯曲现象出现,弯曲最大现象出现在较大叶片半径处 ,整个叶轮振型表现为1、3叶片反向扭摆,2、4不动(如图5e) 。
(6) 六阶振动频率为87.6Hz ,叶片表现为沿径向线扭摆 ,叶片型面上有弯曲现象出现,弯曲最大现象出现在较大叶片半径处 ,整个叶轮振型表现为 4 个叶片以相同形式同向扭摆(如图5f) 。
(7) 七阶振动频率为152.2Hz,1、3叶片的运动形式以叶片型面上的弯曲波为主 ,弯曲波有两条节线 ,且弯曲波最严重发生在叶片外周处,1、3叶片运动形式相差180°相位;2、4叶片基本不动 ,但在前掠的叶尖处有少量翘曲;对面的叶片变化相位差 180°(如图5g) 。
(8) 八阶振动频率为152.6Hz ,叶片表现4个叶片型面上都出现弯曲波 ,且整个叶轮相对的两个叶片振动情况相同 ,而相邻的叶片振动情况相差 180°相位(如图5h) 。
可以看出叶轮模态振动形式主要是由 4 个叶片周向不同振动组合形式 ,造成这种原因主要是前掠叶片刚度远小于轮毂的刚度 ,即叶片“软” ,而轮毂“硬”;叶片低频振型主要以叶片整体扭摆为主 ,而高频主要以叶片弯曲波为主。从对噪声的影响来看认为前六阶振型影响较大 ,因为由于它的摆对流场有较大影响,造成了叶片进气攻角的变化 ,从而形成叶片表面升力的波动 ,最严重的情况会产生驰振,产生很大气动噪声和效率较大的降低。
七、结论
通过有限元对叶轮模态进行了模态分析,考虑了旋转软化、应力强化对叶轮真实运行状况下模态频率的影响,发现与静止状态差别不大 ,同时分析结果与试验吻合较好。通过振型分析 ,认为较低频率的前六阶振动对气动噪声影响较大 ,为设计低噪声风机提供了借鉴。
下一步将计算叶轮流场 ,得到叶片流场受力并对叶轮作谐波分析 ,通过气动声学公式预测噪声的大小 ,以期望获得叶轮振动对气动噪声的具体量值。
摘要：为了研究短叶片对后向式离心通风机性能的影响，对某风机厂的JMF№4.3A型离心风机加短叶片的叶轮模型进行了数值研究。数值方法采用RNGk-ε湍流模型和SIMPLEC算法。主要研究了长短叶片数、短叶片长度、短叶片周向位置以及短叶片的安装角度对后向式离心通风机性能的影响。为方便比较，计算结果做成了无因次曲线的形式。结果表明：对后向式离心通风机而言，短叶片长度影响风机的做功能力和效率；短叶片的周向位置则对压力和效率均有较大影响,特别是对效率的影响最大；而短叶片的安装角度则对风机的性能影响较小。
关键词：数值模拟；离心通风机；叶片；性能
0引言
国内外都有采用加短叶片的离心通风机，特别是目前对高效率的追求，使得它的应用更为广泛。但国内主要是应用长短叶片离心风机，而对其进行工程而深入的研究还较少，尤其对于采用长短叶片的后向式离心通风机研究更为少见。采用长短叶片的目的是为了减小叶片进口处的堵塞；减小叶道扩张角以减少分离损失；减小滑移系数而提高做功能力。长短叶片叶轮中加短叶片的形式多种多样，但概括起来可以归结为四个方面的问题：1）长短叶片数对风机性能的影响；2）短叶片的长度对风机性能的影响；3）短叶片的周向位置对风机性能的影响；4）短叶片的安装角度对风机性能的影响[1]。本文采用数值模拟手段研究这四个方面的问题。
本文采用目前较为常用的Fluent商用软件包对长短叶片后向离心通风机短叶片的设计规律进行研究，以供工程设计参考。
Fluent软件包主要包括：前处理器Gambit、求解器Fluent、后处理器Fluent或Tecplot。Gambit具有前处理建模和划分网格的功能。Fluent求解是Fluent的核心部分，其数值方法是基于有限体积法[2]的。
1模型的验证
模型采用江阴市精亚风机有限公司的JMF№4.3A型后向离心风机，其叶片长度为159mm，空气动力略图见图1。
图1JMF№4.3A型离心风机空气动力略图
计算模型用三维制图软件建模，然后用Gambit的分块网格功能划分网格，单元总数约为80多万。Fluent使用分离求解器求解全三维的N-S方程；采用RNGk-ε湍流模型、标准壁面函数和SIMPLEC算法进行模拟计算。可动区域则使用MRF（MovingReferenceFrame）模型，转速为2900r/min。叶轮为相对于可动区域静止的可动壁面，其它壁面为无滑移边界。边界条件为速度进口，压力出口。计算收敛且稳定后提取计算结果[3]。
为方便比较，把提取的全压和效率用三阶多项式拟合，做成了无因次曲线的形式。图2中的横坐标为流量系数Φ，主纵坐标为风机内效率η，次纵坐标为全压系数Ψ。其表达式分别如下：
式中：zqv为风机流量，m3/s；D2为叶轮外径，ｍ；u2为叶轮外缘圆周速度，ｍ/s；ptF为风机全压，Pa；ρ为气体密度，kg/m3；Pin为风机内功率，kW。
为使结论具有普遍意义，定义短叶片相对长度ε=短叶片长度/原叶片长度。本文所研究的风机原型是短叶片相对长度ε＝3/4，长短叶片数各为6的叶轮模型，此模型已通过江苏省组织的科技鉴定与新产品新技术鉴定，并投入批量生产。
图2是计算结果与测试值的比较。由图2可以看出，计算得到的全压系数和效率曲线与测试值有所偏差（全压系数偏差约为2%，效率偏差约为3%），但总体趋势较为吻合。这说明本文所采用的计算模型和边界条件用于JMF№4.3A型后向离心风机的模拟计算是可行的。出现这种偏差的原因可能是：由于各种条件的限制，模拟计算所使用的风机叶轮内流动损失模型不能完全符合实际的流动过程，而只是与实际流动过程近似。
图2与原型计算结果的比较
2模拟计算研究
本文通过四个方面分别进行研究。FLUENT计算模型与验证实验所采用的模型一致，但短叶片结构参数有所不同。
2.1叶片数目的影响
为研究不同叶片数目对风机性能的影响，本文采用了长短叶片数各为5、6和8的计算模型。其中叶片数为5时只做了ε=3/4的叶轮模型的模拟计算，而6叶片和8叶片则分别做了ε=3/4和ε=2/3的两种模型的模拟计算[4]。其计算结果见图3和图4。
图3ε=3/4模型的风机特性曲线 图4ε=2/3模型的风机特性曲线
由图3可以看出，叶片数为5时，其压力和效率明显降低。这是因为叶片数过少时，叶片的做功能力明显降低，流道过宽而导致分离损失加大，从而使得效率也大幅降低。叶片数为8时，其压力比6叶片数高很多，但其效率却较低。为此而研究了叶片相对长度ε=2/3，叶片数为6和8的计算模型，计算结果见图4。
图4表明，8叶片数的压力比6叶片数时高，而效率同样也是降低的。这是由于叶片数的增加，也增大了叶片的做功能力，但做功能力的增加，被其所带来的流动损失和摩擦损失所抵消，因而其效率反而有所降低。
2.2短叶片长度的影响
参考文献[1]的工程经验，一号轮为带有不同长度的短叶片，长短叶片数为6或8，周向均布，见图5。计算结果特性曲线见图6和图7。
图66叶片数不同长度短叶片时风机特性曲线图78叶片数不同长度短叶片时风机特性曲线
图6和图7都表明，随着短叶片长度的增加，叶轮的全压系数会增大，即做功能力增强，但增加幅度有所减小。图6中，随着短叶片长度的减小，最高效率点向小流量方向偏移，而总体效率只是在相对长度ε=2/3时仍然较高。而图7中，效率会随着短叶片长度的减小而降低，最高效率点也相对一致。
因此，随着短叶片长度的增加，压力会升高，即叶轮做功能力增大，但效率的变化趋势与叶片数有一定关系。鉴于上面的计算结果分析，此后以长短叶片数各为6的叶轮模型为研究对象。
2.3短叶片周向位置的影响
二号轮为短叶片周向位置不同的叶轮，见图8（P表示压力面，S表示吸力面），图中只给出了3/10τ、4/10τ和1/2τ（τ为栅距）的情况，其余类推。图9和图10分别给出了计算全压和效率特性曲线。由图可以看出，当短叶片偏离中心面时（1/2τ面），其压力和效率都降低，只是效率降低的幅度较大。即使是对称于中心面的位置，靠近吸力面处的压力和效率都比靠近压力面的稍高。
图8二号轮 图9短叶片不同周向位置的风机特性曲线1
图10短叶片不同周向位置的风机特性曲线2 图11三号轮
2.4短叶片安装角度的影响
三号轮为短叶片安装角不同的叶轮，短叶片距离压力面τ/2，其相对长度为1/2，相对安装角分别为-8°、0°和8°，见图11。图12给出了三号轮的全压系数和效率曲线。
由图可以看出，改变短叶片的相对安装角，对风机压力影响较大，无论是相对安装角为正还是为负，压力都低于相对安装角为0°时的压力，特别是相对安装角为正时，压力降低的幅度较大；相对安装角对效率的影响较小，但改变了最高效率点，使其向小流量方向移动。
图12改变短叶片相对安装角时风机特性曲线 图13方案比较选择
2.5方案比较与选择
根据前面四个方面的研究，为提取最佳的方案，对各最优方案进行了比较，其特性曲线见图13。进行比较的叶轮模型分别为叶片数为8，短叶片相对长度为3/4的叶轮，以及叶片数为6，短叶片相对长度分别为3/4和2/3的叶轮，各叶片周向均布。
由图可以看出，当叶片数为8，短叶片相对长度为3/4时，其全压系数最大，即这种叶轮的做功能力最强，但效率较低。而当叶片数为6，短叶片相对长度为2/3时，其效率最高，且最高效率点有向小流量方向移动的趋势，压力比短叶相对长度为3/4的叶轮稍低。基于工程实际需要，可在这三种叶轮模型中进行选取。
3 结论
1）后向离心通风机加入短叶后，叶片总数与原设计的最佳叶片数相同或稍多为宜，不可减少。
2）加入的短叶片长度越长，其全压越高，但升高的幅度减少，效率变化则无明显趋势。
3）在一定叶片数下，短叶片长度的减小，会使最高效率点向小流量方向移动。
4）当短叶片偏离中心面时，其压力和效率都会降低，只是效率降低的幅度较大，即使是对称于中心面的位置，靠近吸力面处的压力和效率都比靠近压力面的稍高。
5）改变短叶片的相对安装角，对压力影响较大，特别是相对安装角为正时，压力降低的幅度很大；相对安装角对效率的影响较小，但改变了最高效率点，使其向小流量方向移动。
参 考 文 献
[1]桂利川，常鸿
太阳集团是水帘生产厂家|环保空调生产厂家|屋顶风机厂家|，太阳集团承接规划：猪场降温|车间降温|厂房降温|猪场通风|车间通风|厂房通风|屋顶排风机|屋顶排热|厂房通风降温|车间通风降温|通风换气排热降温工程|屋顶风机安装|负压风机安装|水帘安装|环保空调安装|通风设备安装|通风降温设备|通风系统安装案例|通风降温系统|屋顶通风机|屋顶排风系统
相关的主题文章:
