CN111651925B - 空调管路低频噪声声源识别方法及系统 - Google Patents
空调管路低频噪声声源识别方法及系统 Download PDFInfo
- Publication number
- CN111651925B CN111651925B CN202010504995.0A CN202010504995A CN111651925B CN 111651925 B CN111651925 B CN 111651925B CN 202010504995 A CN202010504995 A CN 202010504995A CN 111651925 B CN111651925 B CN 111651925B
- Authority
- CN
- China
- Prior art keywords
- compressor
- pipeline
- frequency
- sound source
- low
- 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
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000008878 coupling Effects 0.000 claims abstract description 23
- 238000010168 coupling process Methods 0.000 claims abstract description 23
- 238000005859 coupling reaction Methods 0.000 claims abstract description 23
- 239000007787 solid Substances 0.000 claims abstract description 23
- 239000003507 refrigerant Substances 0.000 claims description 12
- 238000004378 air conditioning Methods 0.000 claims description 10
- 230000001133 acceleration Effects 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000006835 compression Effects 0.000 claims description 4
- 238000007906 compression Methods 0.000 claims description 4
- 238000013507 mapping Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 5
- 230000005284 excitation Effects 0.000 description 4
- 238000001816 cooling Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- 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
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L17/00—Speaker identification or verification techniques
- G10L17/26—Recognition of special voice characteristics, e.g. for use in lie detectors; Recognition of animal voices
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Acoustics & Sound (AREA)
- Health & Medical Sciences (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Remote Sensing (AREA)
- Multimedia (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明涉及空调领域,具体涉及一种空调管路低频噪声声源识别方法。能够解决单一结构场仿真无法准确识别管路低频噪声声源位置的问题,同时极大地提高了低频噪声声源识别的准确度。技术方案概括为,计算压缩机旋转惯性力矩,对压缩旋转惯性力矩进行变换,获得各个频率点下旋转惯性力矩的幅值,生成压缩机管路系统的流固耦合有限元模型,将各个频率点下旋转惯性力矩的幅值加载到压缩机管路系统的流固耦合有限元模型上,计算得到压缩机管路系统中各段管路的频率响应曲线,对比各段管路的频率响应曲线的幅值大小来确定低频噪声声源位置。
Description
技术领域
本发明涉及空调领域,具体涉及一种空调管路低频噪声声源识别方法及系统。
背景技术
低频噪音由于其声波波长长,具有很强穿透性,因而在空调行业内,空调器的管路低频噪音问题是行业内的难点。空调器的管路低频噪音大多数情况下是由管路振动引起,因此,目前解决管路低频噪音的手段都是通过降低管路振动。比如,在管路上添加减震块或者腻子,导致管路成本增加。
目前,针对解决低频管路噪音问题,普遍的做法都是以单一结构场的方式,通过压缩机管路系统的固有频率和振型以及管路的响应去进行预判,而这种预判在很大程度上存在误差。比如,空调系统在某一频率点下运行,制冷模式无低频噪音问题,而制热模式存在噪音问题。若仅用单一结构场的仿真方式,无法准确识别和预判产生低频噪声声源位置,也就无法有目的性地去优化管路,进而也就无法实现缩短设计周期的目的。
现有技术如CN110986264A公开的一种空调低频共振噪声识别控制方法及空调,在空调管路和壳体上布置振动传感器以及在外机布置噪声传感器进行数据采集,通过给出噪声和振动的阈值来识别及判断低频共振噪声。此种方式依赖于振动传感器以及噪声传感器采集数据的准确度,又因为低频噪声共振,极易出现误差。
发明内容
本发明的目的是提供一种空调管路低频噪声声源识别方法及系统,能够解决单一结构场仿真无法准确识别管路低频噪声声源位置的问题,同时极大地提高了低频噪声声源识别的准确度。
本发明采取如下技术方案实现上述目的,空调管路低频噪声声源识别方法,包括以下步骤:
步骤(1)、计算压缩机旋转惯性力矩;
步骤(2)、对压缩机旋转惯性力矩进行FFT(Fast Fourier Transformation)变换,获得各个频率点下压缩机旋转惯性力矩的幅值;
步骤(3)、生成压缩机管路系统的流固耦合有限元模型;
步骤(4)、将各个频率点下旋转惯性力矩的幅值加载到压缩机管路系统的流固耦合有限元模型上得到各个管路的频率响应曲线;
步骤(5)、根据各个管路的频率响应曲线来识别低频噪声声源位置。
进一步的是,在步骤(1)中,所述计算压缩机旋转惯性力矩包括:压缩机旋转惯性力矩为时域数据,设置压缩机旋转惯性力矩为M,则M=Md-Mg,其中,Md为驱动力矩,Mg为阻力矩。
进一步的是,所述压缩机转轴压缩制冷剂的有效功率Pd=G(h2-h1),其中G为流量质量,h2为制冷剂在排气压力和温度确定条件下的焓值,h1为制冷剂在吸气压力和温度确定条件下的焓值。
进一步的是,在步骤(2)中,所述获得各个频率点下压缩机旋转惯性力矩的幅值包括:对压缩机旋转惯性力矩进行FFT变换,获得各个频率点下的频谱成分,频谱成分中包括各个频率点下的压缩机旋转惯性力矩的幅值信息。
进一步的是,在步骤(3)中,所述确定压缩机管路系统的流固耦合有限元模型包括:先进行压缩机管路内流场分析,以质量流量进口、压力出口为边界条件,将计算得到的流场结果映射到压缩机管路系统的结构场有限元模型中,进而得到压缩机管路系统的流固耦合有限元模型。
进一步的是,所述频率响应曲线为加速度频率响应曲线或速度频率响应曲线。
进一步的是,在步骤(5)中,所述低频噪声声源位置识别包括:通过对比各段管路的频率响应曲线的幅值大小来确定低频噪声声源位置,该段管路频率响应曲线的幅值越大,表示该段管路为低频噪声声源的概率越高。
空调管路低频噪声声源识别系统,包括:
力矩计算模块,用于计算压缩机旋转惯性力矩;
FFT变换模块,用于对压缩机旋转惯性力矩进行FFT变换,获得各个频率点下压缩机旋转惯性力矩的幅值;
流固耦合模块,用于生成压缩机管路系统的流固耦合有限元模型;
载荷加载模块,用于将各个频率点下旋转惯性力矩的幅值加载到压缩机管路系统的流固耦合有限元模型上得到各个管路的频率响应曲线;
声源识别模块,用于根据各个管路的频率响应曲线来识别低频噪声声源位置。
本发明采用压缩机管路系统的流固耦合有限元模型取代压缩机管路系统的单一场结构模型,解决了单一结构场仿真无法准确识别管路低频噪声声源位置的问题;通过对压缩旋转惯性力矩进行FFT变换,获得各个频率点下旋转惯性力矩的幅值,然后将各个频率点下旋转惯性力矩的幅值加载到压缩机管路系统的流固耦合有限元模型上,计算得到压缩机管路系统中各段管路的频率响应曲线,对比各段管路的频率响应曲线的幅值大小来确定低频噪声声源位置,该段管路频率响应曲线的幅值越大,表示该段管路为低频噪声声源的概率越高;本发明不依赖于振动传感器以及噪声传感器等数据采集设备,通过对各段管路频率响应曲线进行对比分析得出低频噪声声源位置,极大地提高了低频噪声声源识别的准确度。
附图说明
图1是本发明空调管路低频噪声声源识别方法的方法流程图。
图2是本发明一种实施例压缩机管路系统中各管段加速度曲线数据对比图。
图3是本发明一种实施例压缩机管路系统中吸气管噪声优化前后的数据对比图。
具体实施方式
本发明空调管路低频噪声声源识别方法,其方法流程图如图1,包括以下步骤:
步骤101:计算压缩机旋转惯性力矩;
步骤102:对压缩机旋转惯性力矩进行FFT变换,获得各个频率点下压缩机旋转惯性力矩的幅值;
步骤103:生成压缩机管路系统的流固耦合有限元模型;
步骤104:将各个频率点下旋转惯性力矩的幅值加载到压缩机管路系统的流固耦合有限元模型上得到各个管路的频率响应曲线;
步骤105:根据各个管路的频率响应曲线来识别低频噪声声源位置。
步骤101中,计算压缩机旋转惯性力矩的具体实施步骤包括:压缩机旋转惯性力矩为时域数据,设置压缩机旋转惯性力矩为M,则M=Md-Mg,其中,Md为驱动力矩,Mg为阻力矩。
其中驱动力矩其中Pd为压缩机转轴压缩制冷剂的有效功率,f为压缩机的运行频率,压缩机转轴压缩制冷剂的有效功率Pd=G(h2-h1),G为流量质量,h2为制冷剂在排气压力和温度确定条件下的焓值,h1为制冷剂在吸气压力和温度确定条件下的焓值;阻力矩其中R为压缩机偏心转子半径,ε为偏心转子偏心率,h为压缩机气缸高度,θ为压缩机转子转角,pθ为压缩机压缩腔内压力,ps为压缩机吸气腔内压力,也等于压缩机吸气口压力。
步骤102中,获得各个频率点下压缩机旋转惯性力矩的幅值的具体实施手段包括:对压缩机旋转惯性力矩进行FFT变换,获得各个频率点下的频谱成分,频谱成分中包括各个频率点下的压缩机旋转惯性力矩的幅值信息。
具体地,当压缩机运行频率为fHz,经过FFT变换,可得到fHz(1倍频)压缩机旋转惯性力矩的幅值,2fHz(2倍频)压缩机旋转惯性力矩的幅值,……nfHz(n倍频)压缩机旋转惯性力矩的幅值。压缩机为变频压缩机,运行频率范围为10~100Hz。
步骤103中,确定压缩机管路系统的流固耦合有限元模型的具体实施步骤包括:先进行压缩机管路内流场分析,以质量流量进口、压力出口为边界条件,将计算得到的流场结果映射到压缩机管路系统的结构场有限元模型中,进而得到压缩机管路系统的流固耦合有限元模型。
步骤104中,将各个频率点下旋转惯性力矩的幅值加载到压缩机管路系统的流固耦合有限元模型上得到各个管路的频率响应曲线。具体地,可以是:将1倍频压缩机旋转惯性力矩的幅值(假设10Hz~100Hz)加载后,得到1倍频载荷激励下的频响曲线;将2倍频压缩机旋转惯性力矩的幅值(20Hz~200Hz)加载,得2倍频载荷激励下的频响曲线;……将n倍频压缩机旋转惯性力矩的幅值(10nHz~100nHz)加载,得到n倍频载荷激励下的频响曲线。
其中,频率响应曲线为加速度频率响应曲线或速度频率响应曲线。
步骤105中,低频噪声声源位置识别的具体实施步骤包括:通过对比各段管路的频率响应曲线的幅值大小来确定低频噪声声源位置,该段管路频率响应曲线的幅值越大,表示该段管路为低频噪声声源的概率越高。具体的,若在压缩机80Hz运行下管路产生160Hz低频噪声,则对比各段管路在2倍频载荷激励条件下的响应幅值,响应幅值越大,对产生低频噪声的贡献量越大,即该管段是引起该低频噪声的主因。
本发明的一种实施例,压缩机在55-65Hz范围内存在55-65Hz基频传递音,将1倍频下压缩机旋转惯性力矩的最大幅值加载到压缩机管路系统的流固耦合有限元模型上得到如图2所示的吸气管段、排气管段、冷凝器连接管段和截止阀连接管段的加速度响应曲线图,吸气管段加速度幅值大,因此判断出吸气管段是引起该基频噪声的主因,通过优化吸气管段,该基频噪声得以消除,优化后的吸气管段加速度频响曲线与原吸气管段加速度频响曲线对比如图3。
综上所述,本发明解决了单一结构场仿真无法准确识别管路低频噪声声源位置的问题,同时极大地提高了低频噪声声源识别的准确度,进一步能够针对性的对管路进行优化。
Claims (9)
1.空调管路低频噪声声源识别方法,其特征在于,包括以下步骤:
步骤(1)、计算压缩机旋转惯性力矩;
步骤(2)、对压缩机旋转惯性力矩进行FFT变换,获得各个频率点下压缩机旋转惯性力矩的幅值;
步骤(3)、生成压缩机管路系统的流固耦合有限元模型;
步骤(4)、将各个频率点下旋转惯性力矩的幅值加载到压缩机管路系统的流固耦合有限元模型上得到各个管路的频率响应曲线;
步骤(5)、根据各个管路的频率响应曲线来识别低频噪声声源位置,具体包括:通过对比各段管路的频率响应曲线的幅值大小来确定低频噪声声源位置,该段管路频率响应曲线的幅值越大,表示该段管路为低频噪声声源的概率越高。
6.根据权利要求1所述的空调管路低频噪声声源识别方法,其特征在于,在步骤(2)中,所述获得各个频率点下压缩机旋转惯性力矩的幅值包括:对压缩机旋转惯性力矩进行FFT变换,获得各个频率点下的频谱成分,频谱成分中包括各个频率点下的压缩机旋转惯性力矩的幅值信息。
7.根据权利要求1所述的空调管路低频噪声声源识别方法,其特征在于,在步骤(3)中,所述生成压缩机管路系统的流固耦合有限元模型包括:先进行压缩机管路内流场分析,以质量流量进口、压力出口为边界条件,将计算得到的流场结果映射到压缩机管路系统的结构场有限元模型中,进而得到压缩机管路系统的流固耦合有限元模型。
8.根据权利要求1所述的空调管路低频噪声声源识别方法,其特征在于,在步骤(4)中,所述频率响应曲线为加速度频率响应曲线或速度频率响应曲线。
9.空调管路低频噪声声源识别系统,其特征在于,包括:
力矩计算模块,用于计算压缩机旋转惯性力矩;
FFT变换模块,用于对压缩机旋转惯性力矩进行FFT变换,获得各个频率点下压缩机旋转惯性力矩的幅值;
流固耦合模块,用于生成压缩机管路系统的流固耦合有限元模型;
载荷加载模块,用于将各个频率点下旋转惯性力矩的幅值加载到压缩机管路系统的流固耦合有限元模型上得到各个管路的频率响应曲线;
声源识别模块,用于根据各个管路的频率响应曲线来识别低频噪声声源位置,具体包括:通过对比各段管路的频率响应曲线的幅值大小来确定低频噪声声源位置,该段管路频率响应曲线的幅值越大,表示该段管路为低频噪声声源的概率越高。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010504995.0A CN111651925B (zh) | 2020-06-05 | 2020-06-05 | 空调管路低频噪声声源识别方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010504995.0A CN111651925B (zh) | 2020-06-05 | 2020-06-05 | 空调管路低频噪声声源识别方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111651925A CN111651925A (zh) | 2020-09-11 |
CN111651925B true CN111651925B (zh) | 2022-09-16 |
Family
ID=72347139
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010504995.0A Active CN111651925B (zh) | 2020-06-05 | 2020-06-05 | 空调管路低频噪声声源识别方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111651925B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113218049B (zh) * | 2021-05-12 | 2022-04-01 | 四川长虹空调有限公司 | 一种变频空调器压缩机与管路快速匹配方法 |
CN113239489B (zh) * | 2021-05-12 | 2022-12-09 | 四川长虹空调有限公司 | 变频单转子压缩机振动的仿真评价方法 |
CN114152331A (zh) * | 2021-12-02 | 2022-03-08 | 华能山东石岛湾核电有限公司 | 一种高温气冷堆一回路压力容器响声监测定位系统及方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110986264A (zh) * | 2019-11-21 | 2020-04-10 | 珠海格力电器股份有限公司 | 一种空调低频共振噪声识别控制方法及空调 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AT504424B1 (de) * | 2008-02-08 | 2009-05-15 | Avl List Gmbh | Verfahren zur bestimmung der anteile einzelner übertragungswege |
CN105279327A (zh) * | 2015-10-16 | 2016-01-27 | 奇瑞汽车股份有限公司 | 一种汽车车内振动和噪声的计算机仿真预测方法 |
CN107132052B (zh) * | 2017-05-05 | 2020-12-11 | 广州汽车集团股份有限公司 | 电动车辆车内振动、噪声预测方法和装置 |
CN108896258B (zh) * | 2018-03-30 | 2020-05-26 | 四川长虹空调有限公司 | 一种变频压缩机振动载荷计算方法 |
CN109871596A (zh) * | 2018-12-27 | 2019-06-11 | 武汉理工大学 | 新能源汽车电驱动桥振动噪声的主要贡献结构预测方法 |
CN110287569A (zh) * | 2019-06-18 | 2019-09-27 | 珠海格力电器股份有限公司 | 空调管路设计的分析方法和装置 |
-
2020
- 2020-06-05 CN CN202010504995.0A patent/CN111651925B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110986264A (zh) * | 2019-11-21 | 2020-04-10 | 珠海格力电器股份有限公司 | 一种空调低频共振噪声识别控制方法及空调 |
Also Published As
Publication number | Publication date |
---|---|
CN111651925A (zh) | 2020-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111651925B (zh) | 空调管路低频噪声声源识别方法及系统 | |
US11378088B2 (en) | Control system for centrifugal compressor | |
JP6436785B2 (ja) | 空気調和装置 | |
CN1223728A (zh) | 利用模糊逻辑的离心式冷却装置的变速控制 | |
CN109408963B (zh) | 空调器不同环境工况下的配管振动响应函数关系拟合方法 | |
Aprea et al. | Determination of the compressor optimal working conditions | |
CN103649546A (zh) | 压缩机间隙控制 | |
CN110541811B (zh) | 一种快速获取空调压缩机低频转矩补偿角度的方法 | |
CN111639455B (zh) | 变频空调压缩机配管振动仿真方法 | |
US9816742B2 (en) | Variable frequency drive apparatuses, systems, and methods and controls for same | |
Jang et al. | Investigation on noise of rotary compressors using fluid-structure interaction | |
CN117432647A (zh) | 一种大功率高压比离心压缩机性能测试系统及其测试方法 | |
CN109162926B (zh) | 具有声音控制系统的压缩机 | |
CN105026855B (zh) | 对叶轮失速引起的控制不稳进行检测和恢复的方法和系统 | |
Xue et al. | Modeling and experimental investigation of a variable speed drive water source heat pump | |
CN109145498B (zh) | 一种压缩机系统的管路气流脉动计算方法 | |
CN113959049A (zh) | 一种高负荷降载的噪音改善控制方法及变频空调系统 | |
CN114440395A (zh) | 空调降噪控制方法、装置、空调及存储介质 | |
Willie | Conjugate heat transfer simulation and cooling optimization of the flow inside a screw compressor | |
CN113779728A (zh) | 一种适用于电动空压机的疲劳寿命加速测试方法 | |
Lee et al. | Development of capacity modulation compressor based on a two-stage rotary compressor–Part ΙΙ: Performance experiments and P–V analysis | |
CN118246293B (zh) | 一种管路生命周期优化设备及方法 | |
KR100315959B1 (ko) | 2단 원심압축기의 성능시험장치 | |
Mao et al. | Low-Frequency Discontinuous Noise Diagnosis and Reduction of Air Conditioner Outdoor Unit | |
CN111102695B (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 |