离心泵流动诱导噪声的数值预测.docxVIP

PAGE 1 PAGE 1 离心泵流动诱导噪声的数值预测 上述方法均只考虑了由水动力非定常性引起的噪声，并没有很好地考虑由于结构机械振动引起的辐射噪声。叶轮及蜗壳内部产生的噪声往外辐射时，必定要引起蜗壳结构的振动，因此结构与声场相互之间的耦合作用不能忽视。本文通过流体、固体结构及声场的耦合方法来对一离心泵的流淌诱导噪声进行数值求解，期望能为水泵流淌诱导噪声辐射特性的研究供应设计参考。  1离心泵设计参数及计算模型  本文的研究对象是一个中等比转速离心泵，其比转速n=137.6，如图1所示。水泵的主要设计尺寸及性能参数在表1中列出。 图1离心泵计算模型表1离心泵主要设计尺寸及性能参数 2流场及声场计算方法  2.1流场掌握方程  在本文中由于泵内流体为水，并且忽视其温度的变化，因此可以将流体看成不可压缩、绝热的流体来进行计算，能量方程在计算中可以不考虑。连续性方程和动量方程可表示为：  其中，ρ为流体密度，u和p分别表示为速度和压力，t为时间，μ和S分别为动力粘度和源项。  上述方法均只考虑了由水动力非定常性引起的噪声，并没有很好地考虑由于结构机械振动引起的辐射噪声。叶轮及蜗壳内部产生的噪声往外辐射时，必定要引起蜗壳结构的振动，因此结构与声场相互之间的耦合作用不能忽视。本文通过流体、固体结构及声场的耦合方法来对一离心泵的流淌诱导噪声进行数值求解，期望能为水泵流淌诱导噪声辐射特性的研究供应设计参考。  1离心泵设计参数及计算模型  本文的研究对象是一个中等比转速离心泵，其比转速n=137.6，如图1所示。水泵的主要设计尺寸及性能参数在表1中列出。 图1离心泵计算模型表1离心泵主要设计尺寸及性能参数 2流场及声场计算方法  2.1流场掌握方程  在本文中由于泵内流体为水，并且忽视其温度的变化，因此可以将流体看成不可压缩、绝热的流体来进行计算，能量方程在计算中可以不考虑。连续性方程和动量方程可表示为：  其中，ρ为流体密度，u和p分别表示为速度和压力，t为时间，μ和S分别为动力粘度和源项。  式中，K为叶轮旋转一周所需计算的步数，在本文中取为360，即每个时间步下叶轮转过1°，N为叶轮转速,依据表1中的水泵设计参数可以计算得△t=5.747×10-5s。当计算中各物理量最大残差小于10-4时，认为在各虚拟时间步内计算收敛。实际上，由于人耳对频率敏感的范围集中在50～8KHz之间，而依据后续声学计算中的频域响应范围最大值fmax的计算公式：  由之前计算出来的时间步长可计算到的最大响应频率为8700Hz，因此该时间步长既可以满意声学计算频域的要求又不会大幅度增加计算周期。经数值试验发觉，当叶轮转过7～10圈之后，各监测点上的数值变化已呈较为明显的周期性分布，此时将蜗壳内壁面上的压力脉动在下一阶段的非定常计算时将其输出，作为后续声场计算的载荷激励源。  2.3声场计算方法 图3声振耦合计算流程 鉴于声学边界元方法可以很好地解释复杂空间内流淌噪声的辐射、反射、透射以及声场与结构振动的耦合作用，本文的声场分析采用LMSVirtual.lab中的声学边界元耦合模块来完成。图3所示为声振耦合计算流程，首先通过非定常流场计算输出五个完整旋转周期的蜗壳壁面声压数据作为声场计算声源，之后通过对结构进行模态分析获得壳体结构上的模态响应，再在LMSVirtual.lab声学求解软件中将上述二者的计算结果结合声学边界元模型进行耦合求解，获得壳体表面的振速，最终将壳体表面振速作为边界条件计算得到外场辐射声压及指向性等信息。  本文通过CFD结合声振耦合求解方法针对一离心泵所产生的流淌诱导噪声在泵壳外部的辐射状况进行了数值预估。蜗壳内壁面上的非定常压力脉动监测值显示出了蜗舌四周区域的脉动幅值较大，而在远离蜗舌位置处的脉动相对变化较小，在叶片扫掠频率及其谐频下的压力脉动要明显高于其它频率。  2.4壳体结构的模态响应分析  模态是结构在特定固有频率下的振动特性，它可以预估结构在此频段内受到各种振源作用下的实际振动响应。本文采用壳单元并基于有限元方法来进行模态分析，图4为壳体简化模型，其中壳体壁厚为12mm，结构材料为钢，约束位置在叶轮轮毂下侧，如图中圆圈细线所示的位置，设该接触面在各方向上的位移为零。 图4泵壳结构及约束位置 通过对结构进行模态计算可获得各阶固有频率下的振动特性。由于之前已有很多试验已经证明，在离心泵中主要的压力脉动集中在基频以及前4～5阶谐频上，因此相应频率范围内的结构振动特性应予以重视，表2所示为壳体前6阶固有频率值。 表2壳体结构