CN113051663B - 一种管路系统振动噪声评估装置及评估方法 - Google Patents
一种管路系统振动噪声评估装置及评估方法 Download PDFInfo
- Publication number
- CN113051663B CN113051663B CN202110276704.1A CN202110276704A CN113051663B CN 113051663 B CN113051663 B CN 113051663B CN 202110276704 A CN202110276704 A CN 202110276704A CN 113051663 B CN113051663 B CN 113051663B
- Authority
- CN
- China
- Prior art keywords
- vibration
- noise
- pipeline system
- substructure
- pump
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/10—Noise analysis or noise optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开一种管路系统振动噪声评估装置及评估方法,所述系评估方法包括:S1、构建管路系统模型,包括振动噪声源子模型、振动噪声响应子模型;S2、基于噪声源子模型,测量泵、阀运行时的振动、噪声数据;S3、获取管路系统中各子结构的阻抗特性;S4、基于阻抗特性,分别求解各子结构的振动、噪声响应;S5、基于泵、阀的振动、噪声数据,以及各子结构的振动、噪声响应,获取泵、阀到各待评估子结构的振动、噪声传递函数;S6、改变工况,测得泵、阀的振动、噪声数据,通过振动、噪声传递函数计算待评估子结构的振动加速度和声压。本发明能够对不同工况下管路系统中各子结构振动、噪声进行快速准确评估,实施过程简单。
Description
技术领域
本发明涉及船舶减振降噪技术领域,特别是涉及一种管路系统振动噪声评估装置及评估方法。
背景技术
随着舰船技术的迅猛发展,对舰船的隐蔽性要求越来越高。通海管路系统作为舰船的重要组成部分,在其运行过程中,产生的水噪声随着管路向外辐射,极大降低了舰船的隐蔽性。因此,对舰船管路系统的振动噪声水平进行评估,为管路系统的低噪声设计提供技术支撑显得尤为重要。
发明内容
本发明的目的是提供一种管路系统振动噪声评估装置及评估方法,以解决现有技术中存在的技术问题,能够对不同工况下管路系统中各子结构振动、噪声进行快速准确评估,实施过程简单。
为实现上述目的,本发明提供了如下方案:本发明提供一种管路系统振动噪声评估装置,包括模型构建模块、噪声源获取模块、阻抗获取模块、振动噪声响应获取模块、振动噪声评估模块;
所述模型构建模块用于构建管路系统模型;所述管路系统模型包括振动噪声源子模型、振动噪声响应子模型;
所述噪声源获取模块基于所述振动噪声源子模型获取泵、阀的振动、噪声数据;
所述阻抗获取模块基于所述振动噪声源子模型、振动噪声响应子模型,分别获取管路系统中各子结构的阻抗特性;
所述振动噪声响应获取模块根据所述阻抗获取模块获取的管路系统中各子结构的阻抗特性,分别求解管路系统中各子结构的振动、噪声响应;
所述振动噪声评估模块根据所述噪声源获取模块所获取的泵、阀的振动、噪声数据、所述振动噪声响应获取模块所获取的管路系统中各子结构的振动、噪声响应,获取泵、阀到各待评估子结构的振动、噪声传递函数,基于不同工况下的泵、阀的振动、噪声数据,以及所述振动、噪声传递函数,完成管路系统中各子结构振动和噪声的评估。
优选地,所述振动噪声源子模型基于管路系统中的泵、阀结构进行构建;所述振动噪声响应子模型基于管路系统中的其他待评估子结构进行构建。
优选地,所述振动噪声响应获取模块采用阻抗-导纳法求解管路系统中各子结构的振动、噪声响应。
本发明还提供一种管路系统振动噪声评估方法,包括如下步骤:
S1、基于管路系统中的各子结构,构建管路系统模型;所述管路系统模型包括振动噪声源子模型、振动噪声响应子模型;
S2、基于噪声源子模型,测量泵、阀运行时的振动、噪声数据;
S3、获取管路系统中各子结构的阻抗特性;
S4、基于管路系统中各子结构的阻抗特性,分别求解管路系统中各子结构的振动、噪声响应;
S5、基于泵、阀的振动、噪声数据,以及管路系统中各子结构的振动、噪声响应,获取泵、阀到各待评估子结构的振动、噪声传递函数;
S6、改变工况,重复步骤S2测得泵、阀的振动、噪声数据,通过步骤S5中得到的振动、噪声传递函数计算待评估子结构的振动加速度和声压,完成管路系统中各子结构的振动、噪声评估。
优选地,所述步骤S1中,所述振动噪声源子模型基于管路系统中的泵、阀结构进行构建;所述振动噪声响应子模型基于管路系统中的其他待评估子结构进行构建。
优选地,所述步骤S2中,所述噪声数据的测量采用两声源法、两负载法中的一种。
优选地,所述步骤S3中,分别采用实验测试方法、理论计算方法、有限元分析方法获取管路系统中各子结构的阻抗特性。
优选地,所述步骤S4中,采用阻抗-导纳法分别求解管路系统中各子结构的振动、噪声响应。
本发明公开了以下技术效果:
(1)本发明将泵、阀作为作为振动、噪声源,通过阻抗/声阻抗描述振动、噪声源与管路系统中其他子结构的振动、噪声传递特性,实现了管路系统的简化,且通过振动、噪声源与管路系统中其他子结构的振动、噪声传递特性,能够通过测试不同工况下泵、阀的振动、噪声数据,实现不同工况下管路系统中各子结构振动、噪声的快速准确评估,无需求解复杂的方程组,实施过程简单;
(2)本发明在获取管路系统各子结构的阻抗特性时,综合采用了实验测量、理论计算和有限元分析,相比与只采用实验测量的方法,本发明实施简单,同时保证了评估结果的准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明管路系统振动噪声评估装置结构示意图;
图2为本发明管路系统振动噪声评估方法流程图;
图3为本发明实施例中两声源法测量噪声数据的有源双端口理论模型结构示意图;
图4为本发明实施例中通过两声源法测泵的噪声数据的结构示意图;
图5为本发明实施例中管路系统振动与声预报的阻抗-导纳综合法示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
参照图1所示,本实施例提供一种管路系统振动噪声评估装置,包括:模型构建模块、噪声源获取模块、阻抗获取模块、振动噪声响应获取模块、振动噪声评估模块;
所述模型构建模块用于构建管路系统模型;所述管路系统模型包括振动噪声源子模型、振动噪声响应子模型;所述振动噪声源子模型基于管路系统中的泵、阀结构进行构建;所述振动噪声响应子模型基于管路系统中的其他待评估子结构进行构建,例如,管段、挠性软管、消声器、支脚;
所述噪声源获取模块基于所述振动噪声源子模型获取泵、阀的振动、噪声数据;
所述阻抗获取模块基于所述振动噪声源子模型、振动噪声响应子模型,分别获取管路系统中各子结构的阻抗特性;
所述振动噪声响应获取模块根据所述阻抗获取模块获取的管路系统中各子结构的阻抗特性,采用阻抗-导纳法分别求解管路系统中各子结构的振动、噪声响应;
所述振动噪声评估模块根据所述噪声源获取模块所获取的泵、阀的振动、噪声数据、所述振动噪声响应获取模块所获取的管路系统中各子结构的振动、噪声响应,获取泵、阀到各待评估子结构的振动、噪声传递函数,基于不同工况下的泵、阀的振动、噪声数据,以及所述振动、噪声传递函数,完成管路系统中各子结构振动和噪声的评估。
参照图2所示,本实施例提供一种管路系统振动噪声评估方法,具体包括如下步骤:
S1、基于管路系统中的各子结构,构建管路系统模型;
所述管路系统模型包括振动噪声源子模型、振动噪声响应子模型;所述振动噪声源子模型基于管路系统中的泵、阀结构进行构建;所述振动噪声响应子模型基于管路系统中的其他待评估子结构进行构建,例如,管段、挠性软管、消声器、支脚;通过将泵、阀作为振动源和噪声源,用阻抗/声阻抗描述管路系统中其他子结构的振动、噪声传递特性,以达到简化管路系统、获取理论模型的目的。
S2、基于噪声源子模型,测量泵、阀运行时的振动、噪声数据;
本实施例中,振动数据的测量参考ISO9611;噪声数据的测量采用两声源法、两负载法中的一种;
两声源法测量噪声数据的具体方法包括:
在声学元件尺寸远小于声波波长的情况下,将噪声源看作为位于泵的进出口两端口之间的一个点,以进、出端口声压P和体积速度Q为状态变量,建立有源双端口理论模型,如图3所示;
图3中,Pi和Qi为声源入口端声压和体积振速,Po和Qo为声源出口端声压和体积振速,Psi和Pso为声源在入口和出口的声压,用阻抗矩阵表示为:
本实施例中,通过两声源法测泵的噪声数据的方法具体包括:
首先,在待测泵的两端均布置两个水听器和一个外加声源,如图4所示;
最后,根据测得的8个声压信号,通过阻抗矩阵计算泵的被动源参数,即声阻抗;
其中,n=1,2,c表示介质声速,ρ表示介质密度,f表示频率,j表示虚数单位,l1、l2、l3和l4分别表示四个水听器距泵的入口/出口的距离,Si、So表示进出口横截面积。由此可以计算出泵的被动源特性参数进一步关闭外加声源,开启待测泵,测得4个水听器声压信号并结合计算得到的泵的声阻抗矩阵,可计算得到泵的主动源特性参数Psi和Pso:
至此,泵的主动源特性参数和被动源特性参数均通过测量计算获得。
S3、获取管路系统中各子结构的阻抗特性;
本实施例中,对于泵、阀结构,采用实验测试方法获取阻抗特性;对于直管结构,采用理论计算方法获取阻抗特性;对于弯管结构,采用三维有限元软件计算得到阻抗特性。采用理论计算方法获取直管结构的阻抗特性如下式所示:
S4、基于管路系统中各子结构的阻抗特性,采用阻抗-导纳法分别求解管路系统中各子结构的振动、噪声响应。
如图5所示,管路系统包括10个节点,分别为1~10号节点。
首先仅考虑管路中的振动传递,将管路系统划分为一系列的管路单元;其中与1、6、9~10号节点连接的弯管、直管和挠性接管都是用阻抗矩阵来表示其振动传递特性,在这些节点可以直接通过阻抗矩阵的相互叠加获得总阻抗矩阵。这一过程同有限元法中的总刚度矩阵的组装过程类似。由于水泵的传递特性是用导纳矩阵表示,与其相连的单元一般是用阻抗矩阵表示,水泵所在的2~5号节点不能通过阻抗矩阵的相互叠加来获得总阻抗矩阵。在2~5号节点,水泵与管道、隔振器之间的相互作用力是平衡的:
其中:Z22、Z21、Z55、Z56分别是2号节点、5号节点处,管路单元的输入与传递阻抗矩阵;v1、v2、v5、v6分别为1、2、5、6号节点的速度。
其中,z3,z4分别为3号节点、4号节点(即水泵机脚隔振器)的输入阻抗,v3、v4分别为3、4号节点的速度。
根据上式得到:
其中,I为6阶单位矩阵,上标L、R分别表示节点两端单元;vsp(p=2,3,4,5)表示节点p的自由速度;ypq(p,q=2,3,4,5)表示节点p与节点q之间导纳矩阵的元素;z66、z65分别表示6号节点处管路单元的输入与传递阻抗矩阵。矩阵中略去的部分如果遇到泵、阀单元对应节点,重复应用这种方法进行阻抗、导纳矩阵的组装,否则可以直接应用阻抗矩阵相加进行组装。这种方法将泵的导纳矩阵与相邻单元的阻抗矩阵融合起来,形成了以各节点速度为未知量,以泵的自由速度为激励源的线性代数方程组。求解该方程组就可以求得各个节点的振动响应。
从声传播的角度,由于泵、阀及所有的管路单元的声传递特性都用声阻抗矩阵表示,管路系统的声阻抗矩阵可以很方便地由各节点的声压连续条件得到,如下式所示:
以水泵的噪声源特性为基础得到管路声传播计算式:
式中,Ps2、Ps5分别为节点2、节点5处的内部噪声源;apq(p,q=1,2,3,4,5,6)表示节点p与节点q间阻抗矩阵元素;vf1、vf2、vf5、vf6分别为节点1、节点2、节点5、节点6处的流体扰动速度;L、R分别表示节点两端单元。
上面对结构振动与声分别应用了阻抗—导纳综合法。事实上二者在传播过程中是互相耦合的,上述推导忽略了与泵连接的弯管和挠性接管的振动与声的耦合阻抗。考虑耦合阻抗,弯管和挠性接管节点的力、声压及速度有如下关系:
其中:Fi、Fj分别节点i、节点j的力,Pi、Pj分别为节点i、节点j的声压;cij、dij分别是结构与声的耦合阻抗;aij是声阻抗。水泵的两个法兰所在的2、5号节点振动与声压分别应力平衡和声压连续,得到考虑耦合作用的振动与声传播计算表达式:
同理,求解方程组就可以求得各个节点的噪声响应。
S5、基于泵、阀的振动、噪声数据,以及管路系统中各子结构的振动、噪声响应,获取泵、阀到各待评估子结构的振动、噪声传递函数H振动、H噪声;其中,
式中,H表示传递函数,a表示振动加速度;p表示声压。
S6、改变工况,重复步骤S2测得泵、阀的振动、噪声数据,即振动加速度a′泵/阀和声压p′泵/阀,通过步骤S5中得到的振动、噪声传递函数H振动、H噪声,计算待评估子结构的振动加速度a′评估点和声压p′评估点,完成管路系统中各子结构的振动、噪声评估:
a′评估点=H振动*a′泵/阀
p′评估点=H噪声*p′泵/阀。
以上所述的实施例仅是对本发明的优选方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。
Claims (8)
1.一种管路系统振动噪声评估装置,其特征在于,包括模型构建模块、噪声源获取模块、阻抗获取模块、振动噪声响应获取模块、振动噪声评估模块;
所述模型构建模块用于构建管路系统模型;所述管路系统模型包括振动噪声源子模型、振动噪声响应子模型;
所述噪声源获取模块基于所述振动噪声源子模型获取泵、阀的振动、噪声数据;
所述阻抗获取模块基于所述振动噪声源子模型、振动噪声响应子模型,分别获取管路系统中各子结构的阻抗特性;
所述振动噪声响应获取模块根据所述阻抗获取模块获取的管路系统中各子结构的阻抗特性,分别求解管路系统中各子结构的振动、噪声响应;
所述振动噪声评估模块根据所述噪声源获取模块所获取的泵、阀的振动、噪声数据、所述振动噪声响应获取模块所获取的管路系统中各子结构的振动、噪声响应,获取泵、阀到各待评估子结构的振动、噪声传递函数,基于不同工况下的泵、阀的振动、噪声数据,以及所述振动、噪声传递函数,完成管路系统中各子结构振动和噪声的评估。
2.根据权利要求1所述的管路系统振动噪声评估装置,其特征在于,所述振动噪声源子模型基于管路系统中的泵、阀结构进行构建;所述振动噪声响应子模型基于管路系统中的其他待评估子结构进行构建。
3.根据权利要求1所述的管路系统振动噪声评估装置,其特征在于,所述振动噪声响应获取模块采用阻抗-导纳法求解管路系统中各子结构的振动、噪声响应。
4.一种采用如权利要求1-3任一项所述的管路系统振动噪声评估装置的评估方法,其特征在于,包括如下步骤:
S1、基于管路系统中的各子结构,构建管路系统模型;所述管路系统模型包括振动噪声源子模型、振动噪声响应子模型;
S2、基于噪声源子模型,测量泵、阀运行时的振动、噪声数据;
S3、获取管路系统中各子结构的阻抗特性;
S4、基于管路系统中各子结构的阻抗特性,分别求解管路系统中各子结构的振动、噪声响应;
S5、基于泵、阀的振动、噪声数据,以及管路系统中各子结构的振动、噪声响应,获取泵、阀到各待评估子结构的振动、噪声传递函数;
S6、改变工况,重复步骤S2测得泵、阀的振动、噪声数据,通过步骤S5中得到的振动、噪声传递函数计算待评估子结构的振动加速度和声压,完成管路系统中各子结构的振动、噪声评估。
5.根据权利要求4所述的管路系统振动噪声评估方法,其特征在于,所述步骤S1中,所述振动噪声源子模型基于管路系统中的泵、阀结构进行构建;所述振动噪声响应子模型基于管路系统中的其他待评估子结构进行构建。
6.根据权利要求4所述的管路系统振动噪声评估方法,其特征在于,所述步骤S2中,所述噪声数据的测量采用两声源法、两负载法中的一种。
7.根据权利要求4所述的管路系统振动噪声评估方法,其特征在于,所述步骤S3中,分别采用实验测试方法、理论计算方法、有限元分析方法获取管路系统中各子结构的阻抗特性。
8.根据权利要求4所述的管路系统振动噪声评估方法,其特征在于,所述步骤S4中,采用阻抗-导纳法分别求解管路系统中各子结构的振动、噪声响应。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110276704.1A CN113051663B (zh) | 2021-03-15 | 2021-03-15 | 一种管路系统振动噪声评估装置及评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110276704.1A CN113051663B (zh) | 2021-03-15 | 2021-03-15 | 一种管路系统振动噪声评估装置及评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113051663A CN113051663A (zh) | 2021-06-29 |
CN113051663B true CN113051663B (zh) | 2021-11-23 |
Family
ID=76512212
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110276704.1A Active CN113051663B (zh) | 2021-03-15 | 2021-03-15 | 一种管路系统振动噪声评估装置及评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113051663B (zh) |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6697741B2 (en) * | 2002-01-31 | 2004-02-24 | Honda Giken Kogyo Kabushiki Kaisha | Method and system for evaluating and monitoring hydraulic pump noise levels |
US7570989B2 (en) * | 2004-11-22 | 2009-08-04 | Cardiodynamics International Corporation | Method and apparatus for signal assessment including event rejection |
CN103616227B (zh) * | 2013-12-09 | 2015-12-02 | 中国船舶重工集团公司第七〇五研究所 | 管路消声器降噪效果评估装置及评定方法 |
CN110282074B (zh) * | 2019-07-15 | 2020-08-04 | 中国船舶科学研究中心(中国船舶重工集团公司第七0二研究所) | 船舶冷却管路系统变工况低噪声配置调控方法 |
-
2021
- 2021-03-15 CN CN202110276704.1A patent/CN113051663B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113051663A (zh) | 2021-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chandler-Wilde et al. | Efficient calculation of the Green function for acoustic propagation above a homogeneous impedance plane | |
To et al. | A transient testing technique for the determination of matrix parameters of acoustic systems, I: Theory and principles | |
Singh et al. | Mathematical modeling of multicylinder compressor discharge system interactions | |
Fang et al. | Acoustic attenuation analysis of silencers with multi-chamber by using coupling method based on subdomain division technique | |
Trabelsi et al. | Passive and active acoustic properties of a diaphragm at low Mach number: Experimental procedure and numerical simulation | |
CN113051663B (zh) | 一种管路系统振动噪声评估装置及评估方法 | |
Mimani et al. | Acoustic end-correction in a flow-reversal end chamber muffler: A semi-analytical approach | |
Kulkarni et al. | Effect of extended inlet and outlet lengths on transmission loss of double expansion chamber reactive muffler | |
Mimani et al. | Transverse plane-wave analysis of short elliptical end-chamber and expansion-chamber mufflers | |
Du et al. | Acoustic performance of a water muffler | |
Shah et al. | Experimental study on the acoustic properties of perforates under flow using three-port technique | |
US6752018B2 (en) | Method and apparatus for characterizing an acoustic impedance | |
Schulz et al. | Impedance eduction of acoustic liners based on four different levels of physical modeling | |
Okasha et al. | Analysis of pipeline networks using two-ports | |
Chen et al. | A three-dimensional straightforward method for liner impedance eduction in uniform grazing flow | |
Byrne | Calculating the acoustical properties of plane uniform constructions | |
Hornikx et al. | The extended Fourier pseudospectral time-domain (PSTD) method for fluid media with discontinuous properties | |
CN110566471B (zh) | 一种基于功能参数的便捷式泵声学性能获取方法 | |
Bugaru et al. | Transfer matrix method for a single-chamber Mufflers | |
Larsson | Fluid power attenuators analysis, measurements and performance optimization | |
Ingle et al. | Experimental investigation of a double expansion chamber reactive muffler for stationary diesel engine | |
Kuribayashi et al. | Experimental Investigation on the Helmholtz Hydraulic Silencer With Flat Cylindrical Vessel Configuration | |
Chapnik et al. | Noise reduction using finite-length flexible segments | |
Tuozzo et al. | OpenPulse-Open source code for numerical modelling of low-frequency acoustically induced vibration in gas pipeline systems | |
CN111337093B (zh) | 一种超声波流量测量方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |