• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 15
    • 下载费用:30 金币  

    重庆时时彩万能后二: 弹性叶片螺旋桨流噪声预测方法.pdf

    关 键 词:
    弹性 叶片 螺旋桨 噪声 预测 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201611060235.5

    申请日:

    2016.11.23

    公开号:

    CN106777542A

    公开日:

    2017.05.31

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G06F 17/50申请日:20161123|||公开
    IPC分类号: G06F17/50 主分类号: G06F17/50
    申请人: 大连理工大学
    发明人: 邹丽; 鲁建华; 孙铁志; 宗智; 姜胜超; 于游
    地址: 116024 辽宁省大连市高新园区凌工路2号
    优先权:
    专利代理机构: 大连东方专利代理有限责任公司 21212 代理人: 阎昱辰;李洪福
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201611060235.5

    授权公告号:

    |||

    法律状态公告日:

    2017.06.23|||2017.05.31

    法律状态类型:

    实质审查的生效|||公开

    摘要

    本发明公开了一种弹性叶片螺旋桨流噪声预测方法,包括如下步骤:—建立待分析的螺旋桨和流体域的全尺寸模型,该模型中包括随螺旋将转动的旋转域和不发生旋转的静止域;网格化所述的全尺寸模型;设置噪声的采样点;为网格化的全尺寸模型设置流场和声场的求解参数;—为所述的全尺寸模型的更新动态网格;设置固体部分螺旋桨的材料属性。网格划分和耦合面参数;设置流固耦合计算过程参数和流固耦合控制方程;其中:为应力;n为法向向量;d为位移;下标f和s分别代表流体与固体;进行流固耦合计算,直至算法收敛;获取所述监测点处流场数据,根据FW??H方程计算监测点处的声压级;依据快速傅里叶变换FFT功能得到声场的功率谱密度曲线。

    权利要求书

    1.一种弹性叶片螺旋桨流噪声预测方法,其特征在于包括如下步骤:
    —建立待分析的螺旋桨和流体域的全尺寸模型,该模型中包括随螺旋将转动的旋转域
    和不发生旋转的静止域;网格化所述的全尺寸模型;设置噪声的采样点;为网格化的全尺寸
    模型设置流场和声场的求解参数;
    —为所述的全尺寸模型的更新动态网格;设置固体部分螺旋桨的材料属性。网格划分
    和耦合面参数;设置流固耦合计算过程参数和流固耦合控制方程;
    <mfenced open = "" close = "}"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&tau;</mi> <mi>f</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>n</mi> <mi>f</mi> </msub> <mo>=</mo> <msub> <mi>&tau;</mi> <mi>s</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>n</mi> <mi>s</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>d</mi> <mi>f</mi> </msub> <mo>=</mo> <msub> <mi>d</mi> <mi>s</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced>
    其中:τ为应力;n为法向向量;d为位移;下标f和s分别代表流体与固体;
    —进行流固耦合计算,直至算法收敛;获取所述监测点处流场数据,根据FW-H方程计算
    监测点处的声压级;依据快速傅里叶变换FFT功能得到声场的功率谱密度曲线,对所得结果
    进行声学分析,完成弹性叶片螺旋桨流噪声的预测。
    2.根据权利要求1所述的弹性叶片螺旋桨流噪声预测方法,其特征还在于在计算弹性
    模型步骤前还具有验证建立的全尺寸模型正确性的步骤:
    —计算假定为刚性的螺旋桨声场和流场的数据;
    —通过对比得到的刚性螺旋桨与试验的声场和流场数值,验证所述全尺寸模型的正确
    性;
    —若满足正确性要求,则继续计算弹性桨流场和声场数据。
    3.根据权利要求1所述的弹性叶片螺旋桨流噪声预测方法,其特征还在于根据螺旋桨
    投影原理及坐标转换公式完成螺旋桨叶切面二维平面坐标至三维空间坐标的转换,进而完
    成螺旋桨叶切面轮廓曲线的绘制;
    <mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>X</mi> </mtd> </mtr> <mtr> <mtd> <mi>Y</mi> </mtd> </mtr> <mtr> <mtd> <mi>Z</mi> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>Z</mi> <mn>1</mn> </msub> <mi>sin</mi> <mi>&phi;</mi> <mo>+</mo> <msub> <mi>X</mi> <mn>1</mn> </msub> <mi>cos</mi> <mi>&phi;</mi> <mo>+</mo> <mi>L</mi> <mi> </mi> <mi>sin</mi> <mi>&phi;</mi> <mo>-</mo> <mi>r</mi> <mi> </mi> <mi>tan</mi> <mi>&theta;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>r</mi> <mi> </mi> <mi>cos</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>Z</mi> <mn>1</mn> </msub> <mi>cos</mi> <mi>&phi;</mi> <mo>-</mo> <msub> <mi>X</mi> <mn>1</mn> </msub> <mi>sin</mi> <mi>&phi;</mi> <mo>+</mo> <mi>L</mi> <mi> </mi> <mi>cos</mi> <mi>&phi;</mi> </mrow> <mi>r</mi> </mfrac> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>r</mi> <mi> </mi> <mi>sin</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>Z</mi> <mn>1</mn> </msub> <mi>cos</mi> <mi>&phi;</mi> <mo>-</mo> <msub> <mi>X</mi> <mn>1</mn> </msub> <mi>sin</mi> <mi>&phi;</mi> <mo>+</mo> <mi>L</mi> <mi> </mi> <mi>cos</mi> <mi>&phi;</mi> </mrow> <mi>r</mi> </mfrac> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
    其中:φ为螺距角;θ为纵倾角;L为最大厚度线和参考线之间距离;r为叶切面半径;X1、
    Y1、Z1为局部坐标系下的坐标值;X,Y,Z为全局坐标系下坐标值。
    4.根据权利要求1所述的弹性叶片螺旋桨流噪声预测方法,其特征还在于所述的流固
    耦合的控制方程为满足连续性方程和纳维-斯托克斯方程N-S方程,方程函数形式如下:
    <mrow> <mfrac> <mrow> <mo>&part;</mo> <mi>&rho;</mi> </mrow> <mrow> <mo>&part;</mo> <mi>t</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <mo>&part;</mo> <mrow> <mo>(</mo> <msub> <mi>&rho;u</mi> <mi>x</mi> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <mo>&part;</mo> <mrow> <mo>(</mo> <msub> <mi>&rho;u</mi> <mi>y</mi> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <mo>&part;</mo> <mrow> <mo>(</mo> <msub> <mi>&rho;u</mi> <mi>z</mi> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&part;</mo> <mi>z</mi> </mrow> </mfrac> <mo>=</mo> <mn>0</mn> </mrow>
    <mfenced open = "" close = "}"> <mtable> <mtr> <mtd> <mrow> <mi>&rho;</mi> <mfrac> <mrow> <msub> <mi>du</mi> <mi>x</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mo>&part;</mo> <mi>p</mi> </mrow> <mrow> <mo>&part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>+</mo> <mi>&rho;</mi> <mi>X</mi> <mo>+</mo> <msub> <mi>&mu;&Delta;u</mi> <mi>x</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&rho;</mi> <mfrac> <mrow> <msub> <mi>du</mi> <mi>y</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mo>&part;</mo> <mi>p</mi> </mrow> <mrow> <mo>&part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>+</mo> <mi>&rho;</mi> <mi>Y</mi> <mo>+</mo> <msub> <mi>&mu;&Delta;u</mi> <mi>y</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&rho;</mi> <mfrac> <mrow> <msub> <mi>du</mi> <mi>z</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mo>&part;</mo> <mi>p</mi> </mrow> <mrow> <mo>&part;</mo> <mi>z</mi> </mrow> </mfrac> <mo>+</mo> <mi>&rho;</mi> <mi>Z</mi> <mo>+</mo> <msub> <mi>&mu;&Delta;u</mi> <mi>z</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced>
    其中:ρ为流体密度;t为时间;ux、uy、uz分别为速度在空间直角坐标系x、y、z轴的分量;p
    为压力;X、Y、Z分别为外力在x、y、z方向上的分量;μ是流体动力粘性系数;Δ是拉普拉斯算
    子。
    5.根据权利要求4所述的弹性叶片螺旋桨流噪声预测方法,其特征还在于流固耦合计
    算过程中湍流场计算采用RNGk-ε湍流模型:
    <mfenced open = "" close = "}"> <mtable> <mtr> <mtd> <mrow> <mi>&rho;</mi> <mfrac> <mrow> <mi>D</mi> <mi>k</mi> </mrow> <mrow> <mi>D</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>y</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <msub> <mi>&alpha;</mi> <mi>k</mi> </msub> <msub> <mi>&mu;</mi> <mrow> <mi>e</mi> <mi>f</mi> <mi>f</mi> </mrow> </msub> <mfrac> <mrow> <mo>&part;</mo> <mi>k</mi> </mrow> <mrow> <mo>&part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>G</mi> <mi>k</mi> </msub> <mo>+</mo> <msub> <mi>G</mi> <mi>b</mi> </msub> <mo>-</mo> <mi>&rho;</mi> <mi>&epsiv;</mi> <mo>-</mo> <msub> <mi>Y</mi> <mi>M</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&rho;</mi> <mfrac> <mrow> <mi>D</mi> <mi>&epsiv;</mi> </mrow> <mrow> <mi>D</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>y</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <msub> <mi>&alpha;</mi> <mi>&epsiv;</mi> </msub> <msub> <mi>&mu;</mi> <mrow> <mi>e</mi> <mi>f</mi> <mi>f</mi> </mrow> </msub> <mfrac> <mrow> <mo>&part;</mo> <mi>&epsiv;</mi> </mrow> <mrow> <mo>&part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>C</mi> <mrow> <mn>1</mn> <mi>&epsiv;</mi> </mrow> </msub> <mfrac> <mi>&epsiv;</mi> <mi>k</mi> </mfrac> <mrow> <mo>(</mo> <msub> <mi>G</mi> <mi>k</mi> </msub> <mo>+</mo> <msub> <mi>C</mi> <mrow> <mn>3</mn> <mi>&epsiv;</mi> </mrow> </msub> <msub> <mi>G</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>C</mi> <mrow> <mn>2</mn> <mi>&epsiv;</mi> </mrow> </msub> <mi>&rho;</mi> <mfrac> <msup> <mi>&epsiv;</mi> <mn>2</mn> </msup> <mi>k</mi> </mfrac> <mo>-</mo> <mi>R</mi> </mrow> </mtd> </mtr> </mtable> </mfenced>
    其中:ρ为流体密度;为随体导数;k为湍流动能;ε为湍流耗散率;μ为流体动力粘性系
    数;αk、αε为湍流普朗特数的倒数;μeff和R为修正参数;Gk和Gb分别为层流速度梯度和浮力引
    起的湍流动能;YM为可压缩流体湍流扩张贡献量;C1ε、C2ε、C3ε为经验常数。
    6.根据权利要求1所述的弹性叶片螺旋桨流噪声预测方法,其特征还在于:所述的求解
    参数至少包括流场部分湍流模型选择、流体密度、旋转轴、转速以及边界条件;
    声学??椴糠种辽侔ㄔ冻∶芏?即流体水的密度、水中声速以及水中的参考声压;
    求解部分至少包括设置求解算法、收敛残差以及时间步长。
    7.根据权利要求1所述的弹性叶片螺旋桨流噪声预测方法,其特征还在于所述的声学
    分析至少包括预报螺旋桨流噪声的声压级大??;
    通过计算得到流噪声声压级大小随螺旋桨进速系数J的变化规律,进而得到螺旋桨流
    噪声;
    <mrow> <mi>J</mi> <mo>=</mo> <mfrac> <mi>v</mi> <mrow> <mi>n</mi> <mi>D</mi> </mrow> </mfrac> </mrow>
    其中v为螺旋桨速度,计算中为进口水的流速,n为转速,D为螺旋桨直径。

    说明书

    弹性叶片螺旋桨流噪声预测方法

    技术领域

    本发明涉及一种弹性叶片螺旋桨流噪声的预测方法,设计专利分类号G06计算;推
    算;计数G06F电数字数据处理G06F17/00特别适用于特定功能的数字计算设备或数据处理
    设备或数据处理方法G06F17/50计算机辅助设计。

    背景技术

    螺旋桨噪声是潜艇噪声的主要部分,按生成机理可分为桨叶振动噪声、空泡噪声、
    流噪声,其中桨叶振动噪声和空泡噪声是螺旋桨噪声的主要部分,而且随着螺旋桨优化设
    计理论的逐步成熟,桨叶的振动噪声和空泡噪声都能够通过优化螺旋桨设计得到一定程度
    的解决。螺旋桨流噪声是指湍流与螺旋桨相互作用而产生压力脉动及叶梢处涡脱落等而产
    生的噪声,其产生的原因主要是水流与桨的相互作用。对于潜艇来说,其低速航行时,螺旋
    桨流噪声所占比重虽然很小,但它却对声呐信噪比产生很大的影响。螺旋桨流噪声会增加
    声呐的背景噪声,影响声呐对目标声信号的接收,从而影响了潜艇的作战能力。潜艇在高速
    航行时,噪声主要为流噪声,从而降低了潜艇的隐蔽性,增加了暴露的风险。因此,通过对螺
    旋桨流噪声的研究来优化螺旋桨设计,减小螺旋桨噪声,对于提升潜艇的性能有重要意义。

    由于螺旋桨流噪声试验水听器安装困难、试验中存在声反射、透射等问题,因此国
    内外对螺旋桨流噪声的研究主要以数值模拟为主,如Seol H将势流理论与声类比方法结
    合,通过数值方法预报螺旋桨的空泡和非空泡噪声;谢建波等人在Lighthill声类比理论和
    FW-H方程的基础上,应用带平均流效果的格林函数对螺旋桨辐射线谱噪声进行了理论和数
    值分析;龚京风等人运用大涡模拟和K-FWH方程相结合的方法,对螺旋桨水动力噪声进行了
    预报,证明了大涡模拟结合K-FWH模型进行螺旋桨噪声预报的可行性。

    但是在过去的研究中,学者们主要将螺旋桨视为刚性体或者只进行流体-结构的
    单向耦合,即先计算流场,然后将压力结果传递到结构中,并不考虑流体-结构之间的双向
    耦合作用。随着复合材料的发展及人们对柔性表面降噪研究的深入,通过改变螺旋桨叶片
    的弹性来减小螺旋桨流噪声成为科研人员的新课题。采用弹性叶片的螺旋桨必然不能再视
    为刚性体,这就要求在数值模拟的过程中考虑流体-结构之间的双向耦合作用,耦合后的流
    场计算完毕之后进行流噪声的分析。

    发明内容

    本发明针对以上问题的提出,而研制的一种弹性叶片螺旋桨流噪声预测方法,包
    括如下步骤:

    —建立待分析的螺旋桨和流体域的全尺寸模型,该模型中包括随螺旋将转动的旋
    转域和不发生旋转的静止域;网格化所述的全尺寸模型;设置噪声的采样点;为网格化的全
    尺寸模型设置流场和声场的求解参数;

    —为所述的全尺寸模型的更新动态网格;设置固体部分螺旋桨的材料属性。网格
    划分和耦合面参数;设置流固耦合计算过程参数和流固耦合控制方程;


    其中:τ为应力;n为法向向量;d为位移;下标f和s分别代表流体与固体;

    —进行流固耦合计算,直至算法收敛;获取所述监测点处流场数据,根据FW-H方程
    计算监测点处的声压级;依据快速傅里叶变换FFT功能得到声场的功率谱密度曲线,对所得
    结果进行声学分析,完成弹性叶片螺旋桨流噪声的预测。

    作为优选的实施方式,计算弹性模型步骤前还具有验证建立的全尺寸模型正确性
    的步骤:

    —计算假定为刚性的螺旋桨声场和流场的数据;

    —通过对比得到的刚性螺旋桨与试验的声场和流场数值,验证所述全尺寸模型的
    正确性;

    —若满足正确性要求,则继续计算弹性桨流场和声场数据。

    作为优选的实施方式,所述的螺旋桨投影原理及坐标转换公式完成螺旋桨叶切面
    二维平面坐标至三维空间坐标的转换,进而完成螺旋桨叶切面轮廓曲线的绘制;


    其中:φ为螺距角;θ为纵倾角;L为最大厚度线和参考线之间距离;r为叶切面半
    径;X1、Y1、Z1为局部坐标系下的坐标值;X,Y,Z为全局坐标系下坐标值。

    作为优选的实施方式,所述的流固耦合的控制方程为满足连续性方程和纳维-斯
    托克斯方程N-S方程,方程函数形式如下:



    其中:ρ为流体密度;t为时间;ux、uy、uz分别为速度在空间直角坐标系x、y、z轴的分
    量;p为压力;X、Y、Z分别为外力在x、y、z方向上的分量;μ是流体动力粘性系数;Δ是拉普拉
    斯算子。

    更进一步的,流固耦合计算过程中湍流场计算采用RNGk-ε湍流模型:


    其中:ρ为流体密度;为随体导数;k为湍流动能;ε为湍流耗散率;μ为流体动力粘
    性系数;αk、αε为湍流普朗特数的倒数;μeff和R为修正参数;Gk和Gb分别为层流速度梯度和浮
    力引起的湍流动能;YM为可压缩流体湍流扩张贡献量;C1ε、C2ε、C3ε为经验常数。

    作为优选的实施方式,所述的求解参数至少包括流场部分湍流模型选择、流体密
    度、旋转轴、转速以及边界条件;

    声学??椴糠种辽侔ㄔ冻∶芏?即流体水的密度、水中声速以及水中的参考声
    压;

    求解部分至少包括设置求解算法、收敛残差以及时间步长。

    作为优选的实施方式,所述的声学分析至少包括预报螺旋桨流噪声的声压级大
    ??;

    通过计算得到流噪声声压级大小随螺旋桨进速系数J的变化规律,进而得到螺旋
    桨流噪声;


    其中v为螺旋桨速度,计算中为进口水的流速,n为转速,D为螺旋桨直径。

    附图说明

    为了更清楚的说明本发明的实施例或现有技术的技术方案,下面将对实施例或现
    有技术描述中所需要使用的附图做一简单地介绍,显而易见地,下面描述中的附图仅仅是
    本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可
    以根据这些附图获得其他的附图。

    图1为本发明主要流程图

    图2为螺旋桨投影原理示意图

    图3为螺旋桨三维模型示意图

    图4为旋转域及静止域示意图

    图5为所计算流体域尺寸示意图

    图6为流体域与结构域网格划分示意图

    图7为螺旋桨敞水性能曲线计算值与试验值对比图

    图8为螺旋桨流噪声随进速系数变化曲线图

    具体实施方式

    为使本发明的实施例的目的、技术方案和优点更加清楚,下面结合本发明实施例
    中的附图,对本发明实施例中的技术方案进行清楚完整的描述:

    如图1-6所示:

    步骤一几何建模:本步骤中采用的软件为Solidworks,根据图2所示的螺旋桨投影
    原理示意图及坐标转换公式:


    其中:φ为螺距角;θ为纵倾角;L为最大厚度线和参考线之间距离;r为叶切面半
    径;X1、Y1、Z1为局部坐标系下的坐标值;X,Y,Z为全局坐标系下坐标值。

    可以将DTMB 4119螺旋桨叶切面二维平面坐标转换成三维空间坐标,从而在
    Solidworks中得到螺旋桨叶切面轮廓曲线,

    然后以各叶切面导边和随边线作为放样曲线进行曲面放样。再将所形成的曲面进
    行修补,得到封闭的曲面,随后将封闭的曲面转换为实体得到一个叶片模型。接下来建立桨
    毂模型、对叶片进行阵列操作,完成螺旋桨的建模,如图3所示。

    对螺旋桨进行数值模拟属于典型的旋转机械问题,因此对流体域的建模采用旋转
    机械的一般建模方法,即建立静止域和旋转域,

    其中旋转域是在模拟中随着螺旋桨一起转动的一部分,静止域不发生旋转。流体
    域的建模主要是在SolidWorks中分别作出与静止域和旋转域直径相同的圆柱体。然后在静
    止域圆柱上进行布尔运算切除掉旋转域圆柱,形成最终的静止域、在旋转域圆柱上进行布
    尔运算切除掉螺旋桨,形成最终的旋转域,如图4所示。将旋转域和静止域进行装配操作,得
    到最终的计算流体域,尺寸如图5所示。将所建立的螺旋桨和流体域模型以.stp文件导出,
    以供后续使用。

    在得到模型之后,还具有验证正确性的步骤,在本实施例中,计算所得螺旋桨敞水
    性能曲线(如图7,包括推力系数、转矩系数、效率)与试验值进行对比,实施例中计算结果相
    对误差小于5%。一般来说,数值模拟误差达5%以内已经较精确,所以可以判定模型是正确
    的。

    步骤二流体域网格划分:本步骤运用ANSYS Workbench中meshing??槎粤魈逵蚪?br />行网格划分?;质逼蒙栉狢FD分析;

    考虑到流固耦合分析在流体计算时需要进行动网格的设置,同时本身模型的结构
    形式较复杂,为了使所划分的网格更好表达螺旋桨形状,采用非结构化四面体网格对流体
    域进行划分,螺旋桨处的网格适当加密,尤其指在如叶梢部分流态变化比较大的地方网格
    要比其他地方密,以保持计算精度。

    流体域网格划分如图6所示。

    步骤三设置CFD计算软件:此步骤中CFD计算满足质量守恒与动量守恒定律,即满
    足连续性方程和N-S方程,因此基本控制方程为:



    其中:ρ为流体密度;t为时间;ux、uy、uz分别为速度在空间直角坐标系x、y、z轴的分
    量;p为压力;X、Y、Z分别为外力在x、y、z方向上的分量;μ是流体动力粘性系数;Δ是拉普拉
    斯算子。

    首先选择瞬态计算,然后在ANSYS Fluent软件的Model选项Viscous中选择RNGk-ε
    湍流模型,其余保持默认。

    RNGk-ε模型在有漩涡的湍流场计算中有更高的精度,其方程为:


    其中:ρ为流体密度;为随体导数;k为湍流动能;ε为湍流耗散率;μ为流体动力粘
    性系数;αk、αε为湍流普朗特数的倒数;μeff和R为修正参数;Gk和Gb分别为层流速度梯度和浮
    力引起的湍流动能;YM为可压缩流体湍流扩张贡献量;C1ε、C2ε、C3ε为经验常数。

    Cell Zone Condition中将旋转域静止域设为流体,旋转域勾选MRF(即多重参考
    系)选项,根据所研究工况对旋转轴、转速等参数进行设置。

    边界条件通过Boundary Condition选项进行设置本发明采用速度进口、压力出
    口,进口速度根据实际工况取值,若实际工况为非均匀流,则需要采用UDF进行设置,出口压
    力设为0,进口、出口处湍流参数根据一般经验选择湍流强度和湍流粘度比,均取2%。旋转
    域与静止域相连的面设为interface,静止域的圆柱面设为静止无滑移壁面,桨叶和桨毂设
    置为无滑移壁面。

    对声场进行处理需要在Model选项中打开Fluent的声学???,即FW-H模型,并进行
    设置。

    声类比理论是Lighthill在研究湍流声激发中将N-S方程变形成非齐次波动方程,
    将流体力学和声学联系在一起而形成的,经过后来的学者改进形成了FW-H方程,其基本形
    式如下:


    其中:为监测点位置;t为时间;pT为厚度噪声;pL为载荷噪声;pQ为四极子噪声。

    流体密度按Fluent数据库中液体水的密度998.2kg/m3设置、水中声速设为1500m/
    s、自由流速度即为进口速度、水中的参考声压设为1×10-6Pa、声源设为桨叶。设置检测点坐
    标,作为后面流噪声采样点,以便预测此点的螺旋桨流噪声。

    本发明的实例中为了与其他文献进行对比是在距桨轴中心轴向距离分别为0、
    0.5R、0.7R,径向距离分别为0.3R、0.7R、R、1.2R、2R处共设置15个采样点来监测噪声。

    以上设置完毕根据实际需要设置Fluent算法、收敛残差、时间步长等参数,进行初
    始的刚性桨计算,得到流场与声场数据与试验结果对比,有利于选择合适的算法及验证模
    型的精确性以便后续计算。本发明实例计算中得到的螺旋桨敞水性能曲线与试验值对比如
    图6所示。

    步骤四动网格设置:在弹性桨计算时,桨叶处于变形状态,会导致流场域的形状发
    生变化,因此在流固耦合第一次计算完毕后需要对流体域的网格进行更新。ANSYS Fluent
    提供了网格弹簧光顺、网格重构和动态层网格三种网格更新方式,其中动态层网格应用于
    六面体网格,而模拟中流体域均为四面体网格,因此只需要运用网格弹簧光顺和网格重构
    两种方式进行网格更新。首先在动网格选项中激活网格弹簧光顺和网格重构功能,弹簧常
    数设为0.6、边界松弛因子设为0.6,其他参数先采用默认值,根据计算的结果进行逐步调
    整。参数设置完毕之后选择动网格区域,计算中主要考虑弹性叶片螺旋桨,所以将桨叶及周
    围旋转域设为动网格区域,其中桨叶设为System Coupling,旋转域设置为Deforming。

    步骤五结构域设置:首先在Transient Structural??橹懈菟扑懵菪暗牟?br />料设置材料属性。然后打开meshing,偏好设为mechanical分析,选择自动划分得到满足计
    算要求的网格如图6所示。接着设置桨叶为流固耦合面、结构的约束以及时间步长。时间步
    长设置为与流体部分计算相同,流固耦合面与Fluent中动网格区域相对应设为桨叶,设置
    桨毂为固支端。

    结构部分只考虑弹性变形,因此控制方程为:

    σ=E·ε

    其中:σ为应力;E为弹性模量;ε为应变。

    步骤六流固耦合设置:System Coupling??榭梢允迪至魈逵虢峁辜扑愕氖萁?br />换,需在此??橹薪旭詈仙柚?,主要是对时间步长、数据交换面及求解顺序进行设置。时
    间步长设置要满足使流体、结构、耦合时间步长完全对应保证耦合计算的准确性。将流体和
    结构求解器中的桨叶设为数据交换面,在流固耦合的过程中进行数据传递。计算中的变形
    过程为螺旋桨在流体中运动时受力导致桨叶变形,随后桨叶的变形又导致流体域的形状发
    生变化,因此求解顺序应该是先流体后固体,即Fluent设为1,Transient Structural设为
    2。设置完毕后,点击Update开始双向流固耦合计算直至收敛。

    一般情况下流固耦合的收敛残差都是动态平衡的,收敛条件可以参考残差但并不
    完全依赖于残差。流固耦合计算由于进行网格更新,所以收敛残差是波动的。在一个时间步
    长内,收敛残差由大变小,下一个时间步长内残差由大变小,残差最大值小于前一个时间
    步,因此随着时间的增长,如果一直是这个规律,收敛残差的最大值会逐渐变小,最终达到
    动态平衡。

    本发明判断达到收敛条件要保证两点:1.残差达到动态平衡;2.在每个时间步内
    的最后一次迭代,计算界面上默认显示的物理量数值基本一致,说明继续算下去的结果是
    基本不变的。

    流固耦合部分应满足流体与结构相交处各变量相等或守恒,本发明的计算不考虑
    温度,因此流固耦合控制方程为:


    其中:τ为应力;n为法向向量;d为位移;下标f和s分别代表流体与固体。

    步骤七声场分析:再次回到ANSYS Fluent???,对声学部分进行处理,读取监测点
    处流场数据,根据FW-H方程计算监测点处的声压级,同时可以根据ANSYS Fluent内置的快
    速傅里叶变换(FFT)功能得到声场的功率谱密度曲线,对所得结果进行声学分析,完成对弹
    性叶片螺旋桨流噪声的预测。

    本发明实例中计算了不同工况下螺旋桨流噪声声压级大小,并找到了流噪声声压
    级大小随螺旋桨进速系数J:


    其中v为螺旋桨速度,计算中为进口水的流速,n为转速,D为螺旋桨直径)的变化规
    律:一般的,随进速系数不断增大,螺旋桨流噪声不断减小。

    以上所述,仅为本发明较佳的具体实施方式,但本发明的?;し段Р⒉痪窒抻诖?,
    任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其
    发明构思加以等同替换或改变,都应涵盖在本发明的?;し段е?。

    关于本文
    本文标题:弹性叶片螺旋桨流噪声预测方法.pdf
    链接地址://www.4mum.com.cn/p-6020970.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 3d组选包胆 白小姐网站六肖中特 广西快3免费计划软件下载 lbkeno规律 幸运28稳赚技巧 北京pk10分析软件免费 吉林快三追号计划软件手机版 飞艇计划稳 时时走势图软件 双色球开奖各省数据 小胆神3b独胆 竞彩自由过关投注 北京pk10规律大小 全网最早36码无错特围网址 时时彩怎么玩稳赚 极速赛车开奖结果