CN114330060A - 一种周期排列任意数量大矩形开口声传递损失的计算方法 - Google Patents
一种周期排列任意数量大矩形开口声传递损失的计算方法 Download PDFInfo
- Publication number
- CN114330060A CN114330060A CN202111616789.XA CN202111616789A CN114330060A CN 114330060 A CN114330060 A CN 114330060A CN 202111616789 A CN202111616789 A CN 202111616789A CN 114330060 A CN114330060 A CN 114330060A
- Authority
- CN
- China
- Prior art keywords
- rectangular opening
- formula
- wave
- sound
- rectangular
- 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.)
- Granted
Links
- 230000005540 biological transmission Effects 0.000 title claims abstract description 86
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000004364 calculation method Methods 0.000 claims abstract description 37
- 238000013016 damping Methods 0.000 claims abstract description 31
- 230000005855 radiation Effects 0.000 claims abstract description 31
- 230000000694 effects Effects 0.000 claims abstract description 30
- 230000008878 coupling Effects 0.000 claims abstract description 19
- 238000010168 coupling process Methods 0.000 claims abstract description 19
- 238000005859 coupling reaction Methods 0.000 claims abstract description 19
- 239000002245 particle Substances 0.000 claims abstract description 19
- 230000014509 gene expression Effects 0.000 claims abstract description 6
- 239000011148 porous material Substances 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 7
- 230000001902 propagating effect Effects 0.000 claims description 7
- 239000011541 reaction mixture Substances 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 230000000149 penetrating effect Effects 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims 1
- 238000012546 transfer Methods 0.000 description 9
- 238000011160 research Methods 0.000 description 5
- 230000000737 periodic effect Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 3
- 230000001808 coupling effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- 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
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种周期排列任意数量大矩形开口声传递损失的计算方法，首先定义多矩形开口的坐标系，然后建立多矩形开口耦合的互辐射阻抗表达式，接着建立多矩形开口两侧截面处的力平衡式，然后计算考虑粘滞阻尼效应的开口内特征阻抗和特征波数表达式，进而建立开口两端的声传递矩阵，基于力平衡式和声传递矩阵，建立各开口声波出射口处质点振速与入射声压关系式，最后求得耦合声传递损失。本发明可以实现考虑开口间互耦合和粘滞阻尼效应的任意数量大矩形开口传递损失的准确计算。
Description
技术领域
本发明涉及物理专业中噪声类领域，尤其涉及一种能够获得同时考虑粘滞阻尼效应和孔间耦合效应的多矩形开口声传递率及声传递损失的计算方法。
背景技术
周期排列大尺寸矩形开口或矩形狭缝结构在日常生活中随处可见的，如百叶窗的叶片形成的周期排列的矩形狭缝结构等。目前已有研究团队针对单一的矩形开口声传递特性提供了快速的简便计算方法，但是对于百叶窗形成的周期性排列的大尺寸狭长矩形开口而言，还需考虑热粘滞损失以及开口间耦合辐射带来的影响。相似的研究在周期排列的矩形开口中已有进行，如航空发动机声衬、墙体材料护面层甚至手机扬声器的出声口等。这些矩形开口往往单个尺寸不大并且距离较近，多数结构中孔的大小一致。对于周期排列等尺寸小矩形开口，如微穿孔板，国内外学者及研究团队如马大猷院士，进行过深入且大量的研究。这些研究往往将单个矩形开口等效成声阻抗，然后采用等效电路法计算其传递阻抗，进而求得其吸声系数。但是，对于矩形开口阵列而言，其耦合辐射阻抗计算和声传递特性解析框架仍然较为复杂且难以求解。
发明内容
本发明为克服上述现有技术的不足之处，提出一种周期排列任意数量大矩形开口声传递损失的计算方法，具体是一种可以获得考虑开口间互耦合和粘滞阻尼效应的任意数量多矩形开口传递损失的方法，以期能够灵活获取等尺寸、任意数量、周期排列且考虑粘滞阻尼效应的多矩形开口声传递率及声传递损失计算方法。
本发明为解决技术问题采用如下技术方案：
本发明获得一种考虑开口间互耦合和粘滞阻尼效应的任意数量多矩形开口传递损失计算方法，具体为一种周期排列任意数量大矩形开口声传递损失的计算方法，对于贯穿壁面的多个矩形孔，入射侧矩形开口宽度为2ai1、2ai2、......、2ain，长度为2bi1、2bi2、......、2bin，面积为Si1、Si2、......、Sin，出射侧矩形开口宽度为2ao1、2ao2、......、2aon，长度为2bo1、2bo2、......、2bon，面积为So1、So2、......、Son，任意两个矩形开口的中心距为r12、r13、......、rmn，所述计算方法按如下步骤进行：
步骤a、定义坐标系
以其中第h个矩形开口在入射端横截面的几何中心为坐标原点，以垂直于矩形开口声波入射侧截面并朝向矩形开口声波透射侧的方向为z轴正方向，以平行于矩形开口声波入射侧截面的长度方向一侧为y轴正方向，以平行于矩形开口声波入射侧截面的宽度方向一侧为x轴正方向，建立观察坐标系进行研究。
首先假定平面波斜入射到矩形开口，入射声压pi以θi，
的角度入射到第h个矩形开口，
为声波入射到第h个矩形开口过程中产生的反射声压，而
为在此开口入射端产生的辐射声压，
和
分别为矩形开口入射侧内部声压及其质点振速。
为此开口透射端的辐射声压，
和
为此开口内部透射端口声压及其质点振速。
步骤b、计算任意两个矩形开口的互辐射阻抗
按式(1)计算获得任意两个矩形开口的互辐射阻抗
式(1)中j为虚数单位，ρ0＝1.2kg.m-3kg·m-3为空气密度，c0＝340m·s-1为空气中的声速，
为第i个矩形开口的特征波数，r为第i个矩形开口与第h个矩形开口的中心距，根据模态正交性，只考虑m＝p，n＝q的情况，
表示第i个矩形开口与第h个矩形开口在z＝0处的互辐射阻抗，当i＝h时，
为第h个矩形开口在声波入射端的自辐射阻抗；
步骤c、计算矩形开口声波入射侧截面处的力
第h个矩形开口入射端辐射声压的导致的力
可表达为：
式(2)中，
为第i个矩形开口入射端产生的辐射声压，S0h为入射侧第h个矩形开口的横截面积，
为第h个矩形开口内第(m，n)阶高阶波所对应的截面模态振型；
按式(3)计算获得第h个矩形开口声波入射侧截面处的力平衡式如下：
式(3)中，
为入射声压，
为声波入射到第h个矩形开口过程中产生的反射声压，
为第h个矩形开口入射侧内部声压；
其中，第h个矩形开口的内部声压
及内部质点振速
的计算公式为：
式(4)和(5)中，
和
为第(m，n)阶高阶波的模态声压和质点振速，
为第h个矩形开口在z轴方向上传播的第(m，n)阶高阶波的波数，
和
分别表示孔隙内传播介质中的特征阻抗和特征波数；当需考虑粘滞阻尼效应时，应代入考虑粘滞阻尼效应的Zc和kc，即：
式(5.1)中，特征阻抗Zc和特征波数kc的表达式为：
式(6)和(7)中，j为虚数单位，S0为入射侧矩形开口的横截面积，
Y′＝jωC(ω)S0 (9)
式(8)和(9)中，j为虚数单位，S0为入射侧矩形开口的横截面积，ω为角频率，有：
式(10)和(11)中，γ＝1.4为介质的比热比，ρ0＝1.2kg.m-3kg·m-3为空气密度，c0＝340m·s-1为空气中的声速，μ＝1.86×10-5Pa·s为剪切粘度，Cp＝1010J·kg-1·K-1为恒温比热，λ＝0.026W·m-1·K-1为导热系数，此外，对于横截面为矩形的孔隙或者开口，横截面形状影响系数Fr(η)的计算公式为：
在式(12)中，j为虚数单位，ω为角频率；在实际计算时，求和上限s和t可通过收敛性分析获得，且有：
然后把
分别代入式(12)，再把计算结果分别代入式(10)、式(11)，得到对应的复密度ρ(ω)和复压缩系数函数C(ω)；
步骤d、计算矩形开口声波入射侧截面和出射侧截面的声传递矩阵
按式(3)计算获得第h个矩形开口第(m，n)阶高阶波的模态声压
和质点振速
在声波入射端截面(z＝0处)和透射端截面(z＝l处)的声传递矩阵为：
式中Ah、Bh、Ch、Dh分别为第h个矩形开口声波入射端截面和透射端截面的声传递矩阵的四个元素；Ah、Bh、Ch、Dh中包含矩形开口内声传播介质的特征阻抗Zf和特征波数kf，当需考虑粘滞阻尼效应时，应代入考虑粘滞阻尼效应的Zc和kc；
为声波入射侧截面处第h个矩形开口内侧的声压，
为声波出射侧截面处第h个矩形开口内侧的声压，
为声波出射侧截面第h个矩形开口内的质点振速，
为声波出射侧截面第h个矩形开口内的质点振速；
步骤e、计算矩形开口内部透射端口声压
计算获得第h个矩形开口声波出射侧截面处的力平衡式，
式(16)中Slh为透射端第h个矩形开口的横截面积，
为第i个矩形开口与第h个矩形开口在z＝l处(即声波透射端)的互辐射阻抗，当i＝h，
为第h个矩形开口在声波透射端的自辐射阻抗，
为在z＝l处，即声波透射端第h个矩形开口内的质点振速；
因为互模态辐射阻抗可以忽略不计，式(9)可简化为：
由上式可得
式(18)中
步骤f、计算矩形开口声传递率
按式(20)计算获得第n个矩形开口的入射声功率，声源的声功率为：
联立式(4)、(5)、(20)，即可求得透射声功率：
计算获得全体以
的角度入射到H个矩形开口的声波耦合声传递率，
为获得式(22)中
的值，需计算多个
通过联立式(2)、(3)、(4)、(15)并将对应的矩形开口的截面尺寸代入，多次重复上述计算即可求得全部的
最终可以按式(23)计算获得考虑粘滞阻尼效应的多矩形开口耦合声传递损失，
TL＝-10lg(τ) (23)。
进一步的，步骤c中第h个矩形开口的内部质点振速
在不考虑热粘性损失时，那么，Zf＝ρ0c0，kf＝ω/c0，即
同时，步骤d的式(15)中，Ah、Bh、Ch、Dh中包含矩形开口内声传播介质的特征阻抗Zf和特征波数kf，在不考虑热粘性损失的情形下，Zf＝ρ0c0，kf＝ω/c0。
更进一步的，步骤b的式(1)中，r的值为远远大于第i个矩形开口与第h个矩形开口宽度之和的2倍的任意值。
与已有技术相比，本发明有益效果体现在：
1、本申请所述方法可实现考虑孔间耦合和粘滞阻尼效应的多矩形开口声传递损失的准确计算：考虑孔间耦合体现在在式(16)中，当考虑孔间耦合时代入了式(1)所得的互辐射阻抗
考虑粘滞阻尼效应体现在式(5)、式(15)中，当考虑粘滞阻尼效应时孔隙内传播介质中的特征阻抗和特征波数应代入考虑粘滞阻尼效应的Zc和kc；
2、本申请所述方法方法可实现任意数量多矩形开口的耦合声传递损失计算：任意数量多矩形开口体现在本方法适用的矩形开口数量H的取值范围为2至无限大，即式(2)、式(18)和式(22)中的H可以代入2至无限大中的任意整数，即可得到任意数量多矩形开口的耦合声传递损失计算结果；
3、本申请所述方法方法可实现任意间距排列多矩形开口耦合声传递损失计算：任意间距排列多矩形开口体现在式(1)中的第i个矩形开口与第h个矩形开口的中心距r，r的值可以为远远大于第h个矩形开口与第i个矩形开口宽度之和的2倍的任意值；
4、本申请所述方法可实现等尺寸或不等尺寸多矩形开口耦合传递损失计算，等尺寸或不等尺寸的多矩形开口是指第h个矩形开口的入射侧截面长度2aih入射侧截面宽度2bih、出射侧长度2aoh以及出射侧截面尺寸长度2boh和第i个矩形开口相对应的截面尺寸可以相等或不相等，即这两个矩形开口的入射侧截面面积Sin、出射侧截面面积Son可以相等或不相等，在计算中只需将对应的矩形开口的截面尺寸代入式(1)、式(2)、式(3)、式(14)、式(16)、式(17)、式(18)、式(19)、式(20)、式(21)、式(22)。
5、本申请所述方法中的大矩形开口是指在计算该矩形开口声传递特性时须将开口内部高阶波声能传递纳入考虑范围内的矩形开口，体现于将第h个矩形开口内第(m，n)阶高阶波声能所对应的截面模态振型
代入式(2)、式(3)、式(4)、式(5)、式(16)、式(17)、式(19)以及在式(4)、式(5)中代入第(m，n)阶高阶波的模态声压
第(m，n)阶高阶波的质点振速
第h个矩形开口在z轴方向上传播的第(m，n)阶高阶波的波数
故本专利既可用于只考虑平面波(只取m＝n＝0)的矩形开口，也可用于考虑高阶波(m、n为0到无限大的任意正整数)的大矩形开口。
综上，本申请可以获得同时考虑开口间互耦合和粘滞阻尼效应的任意数量多矩形开口传递损失的方法，并能够灵活获取等尺寸、任意数量、周期排列且考虑粘滞阻尼效应的多矩形开口声传递率及声传递损失计算方法。
附图说明
图1为本发明所述任意数量多矩形开口示意图；
图2为本发明实施例中，考虑粘滞阻尼效应和孔间耦合效应的周期排列的十五个等尺寸矩形开口传递损失的声传递损失计算验证结果图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例1
参见图1，本申请公开的一种周期排列任意数量大矩形开口声传递损失的计算方法中，考虑开口间互耦合和粘滞阻尼效应的任意数量多矩形开口传递损失计算方法，对于贯穿壁面的多个矩形孔，处在壁面一侧的是矩形开口声波入射侧，处在壁面另一侧的是矩形开口声波出射侧，入射侧矩形开口宽度为2ai1、2ai2、......、2ain，长度为2bi1、2bi2、......、2bin，面积为Si1、Si2、......、Sin，出射侧矩形开口宽度为2ao1、2ao2、......、2aon，长度为2bo1、2bo2、......、2bion，面积为So1、So2、......、Son，任意两个矩形开口的中心距为r12、r13、......、rmn，其特征是所述计算方法按如下步骤进行：
步骤a、定义坐标系
以其中第h个矩形开口在入射端横截面的几何中心为坐标原点，以垂直于矩形开口声波入射侧截面并朝向矩形开口声波透射侧的方向为z轴正方向，以平行于矩形开口声波入射侧截面的长度方向一侧为y轴正方向，以平行于矩形开口声波入射侧截面的宽度方向一侧为x轴正方向，建立观察坐标系进行研究。
首先假定平面波斜入射到矩形开口，入射声压pi以θi，
的角度入射到第h个矩形开口，
为声波入射到第h个矩形开口过程中产生的反射声压，而
为在此开口入射端产生的辐射声压，
和
分别为矩形开口入射侧内部声压及其质点振速。
为此开口透射端的辐射声压，
和
为此开口内部透射端口声压及其质点振速。
步骤b、计算任意两个矩形开口的互辐射阻抗
按式(1)计算获得任意两个矩形开口的互辐射阻抗
式(1)中j为虚数单位，ρ0＝1.2kg.m-3kg·m-3为空气密度，c0＝340m·s-1为空气中的声速，
为第i个矩形开口的特征波数，r为第i个矩形开口与第h个矩形开口的中心距，步骤b的式(1)中，r的值为远远大于第i个矩形开口与第h个矩形开口宽度之和的2倍的任意值。根据模态正交性，只考虑m＝p，n＝q的情况，
表示第i个矩形开口与第h个矩形开口在z＝0处的互辐射阻抗，当i＝h时，
为第h个矩形开口在声波入射端的自辐射阻抗；
步骤c、计算矩形开口声波入射侧截面处的力
第h个矩形开口入射端辐射声压的导致的力
可表达为：
式(2)中，
为第i个矩形开口入射端产生的辐射声压，S0h为入射侧第h个矩形开口的横截面积，
为第h个矩形开口内第(m，n)阶高阶波所对应的截面模态振型；
按式(3)计算获得第h个矩形开口声波入射侧截面处的力平衡式如下：
式(3)中，
为入射声压，
为声波入射到第h个矩形开口过程中产生的反射声压，
为第h个矩形开口入射侧内部声压；
其中，第h个矩形开口的内部声压
及内部质点振速
的计算公式为：
式(4)和(5)中，
和
为第(m，n)阶高阶波的模态声压和质点振速，
为第h个矩形开口在z轴方向上传播的第(m，n)阶高阶波的波数，
和
分别表示孔隙内传播介质中的特征阻抗和特征波数；当需考虑粘滞阻尼效应时，应代入考虑粘滞阻尼效应的Zc和kc，即：
式(5.1)中，特征阻抗Zc和特征波数kc的表达式为：
式(6)和(7)中，j为虚数单位，S0为入射侧矩形开口的横截面积，
Y′＝jωC(ω)S0 (9)
式(8)和(9)中，j为虚数单位，S0为入射侧矩形开口的横截面积，ω为角频率，有：
式(10)和(11)中，γ＝1.4为介质的比热比，ρ0＝1.2kg.m-3kg·m-3为空气密度，c0＝340m·s-1为空气中的声速，μ＝1.86×10-5Pa·s为剪切粘度，Cp＝1010J·kg-1·K-1为恒温比热，λ＝0.026W·m-1·K-1为导热系数，此外，对于横截面为矩形的孔隙或者开口，横截面形状影响系数Fr(η)的计算公式为：
在式(12)中，j为虚数单位，ω为角频率；在实际计算时，求和上限s和t可通过收敛性分析获得，且有：
然后把
分别代入式(12)，再把计算结果分别代入式(10)、式(11)，得到对应的复密度ρ(ω)和复压缩系数函数C(ω)；
步骤d、计算矩形开口声波入射侧截面和出射侧截面的声传递矩阵
按式(3)计算获得第h个矩形开口第(m，n)阶高阶波的模态声压
和质点振速
在声波入射端截面(z＝0处)和透射端截面(z＝l处)的声传递矩阵为：
式中Ah、Bh、Ch、Dh分别为第h个矩形开口声波入射端截面和透射端截面的声传递矩阵的四个元素；Ah、Bh、Ch、Dh中包含矩形开口内声传播介质的特征阻抗Zf和特征波数kf，当需考虑粘滞阻尼效应时，应代入考虑粘滞阻尼效应的Zc和kc；
为声波入射侧截面处第h个矩形开口内侧的声压，
为声波出射侧截面处第h个矩形开口内侧的声压，
为声波出射侧截面第h个矩形开口内的质点振速，
为声波出射侧截面第h个矩形开口内的质点振速；
步骤e、计算矩形开口内部透射端口声压
计算获得第h个矩形开口声波出射侧截面处的力平衡式，
式(16)中Slh为透射端第h个矩形开口的横截面积，
为第i个矩形开口与第h个矩形开口在z＝l处(即声波透射端)的互辐射阻抗，当i＝h，
为第h个矩形开口在声波透射端的自辐射阻抗，
为在z＝l处，即声波透射端第h个矩形开口内的质点振速；
因为互模态辐射阻抗可以忽略不计，式(9)可简化为：
由上式可得
式(18)中
步骤f、计算矩形开口声传递率
按式(20)计算获得第n个矩形开口的入射声功率，声源的声功率为：
联立式(4)、(5)、(20)，即可求得透射声功率：
计算获得全体以
的角度入射到H个矩形开口的声波耦合声传递率，
为获得式(22)中
的值，需计算多个
通过联立式(2)、(3)、(4)、(15)并将对应的矩形开口的截面尺寸代入，多次重复上述计算即可求得全部的
最终可以按式(23)计算获得考虑粘滞阻尼效应的多矩形开口耦合声传递损失，
TL＝-10lg(τ) (23)。
实施例2
同实施例1不同之处在于，步骤c中第h个矩形开口的内部质点振速
在不考虑热粘性损失时，那么，Zf＝ρ0c0，kf＝ω/c0，即
式中
和
为第(m，n)阶高阶波的模态声压和质点振速，
为第h个矩形开口在z轴方向上传播的第(m，n)阶高阶波的波数，
和
分别表示孔隙内传播介质中的特征阻抗和特征波数，ρ0＝1.2kg.m-3为空气密度，ω为角频率。同时，步骤d中式(15)，Ah、Bh、Ch、Dh中包含矩形开口内声传播介质的特征阻抗ZF和特征波数kf，在计算中也要代入ZF＝ρ0c0，kF＝ω/c0。
方法的检验：
为了验证本实施例中考虑开口间互耦合和粘滞阻尼效应的任意数量多矩形开口传递损失计算方法，对壁面上宽度2a为2mm、长度2b为50mm、深度l为5mm、中心距r为3mm的周期排列的十五个等尺寸矩形开口传递损失的声传递损失计算验证结果与通过声学有限元法仿真得到的结果进行对比。
图2所示，本发明方法计算的结果与通过声学有限元法仿真得到的结果几乎相等，图2中实线表示通过声学有限元法仿真得到的结果，虚线表示本发明所述声传递损失计算公式的计算结果，从图2中可以看出，两条声传递损失曲线的谷完全重合，仅在部分频率范围内，声传递损失的值有偏差，但偏差仍在可以忽略的范围内。
本施例表明，本发明方法能够很好地计算出考虑开口间互耦合和粘滞阻尼效应的任意数量多矩形开口传递损失。
尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种周期排列任意数量大矩形开口声传递损失的计算方法，对于贯穿壁面的多个大矩形开口，处在壁面一侧的是矩形开口声波入射侧，处在壁面另一侧的是矩形开口声波出射侧，入射侧矩形开口宽度为2ai1、2ai2、......、2ain，长度为2bi1、2bi2、......、2bin，面积为Si1、Si2、......、Sin，出射侧矩形开口宽度为2ao1、2ao2、......、2aon，长度为2bo1、2bo2、......、2bon，面积为So1、So2、......、Son，任意两个矩形开口的中心距为r12、r13、......、rmn，其特征在于，所述计算方法按如下步骤进行：
步骤a、定义坐标系
以其中第h个矩形开口在入射端横截面的几何中心为坐标原点，以垂直于矩形开口声波入射侧截面并朝向矩形开口声波透射侧的方向为z轴正方向，以平行于矩形开口声波入射侧截面的长度方向一侧为y轴正方向，以平行于矩形开口声波入射侧截面的宽度方向一侧为x轴正方向，建立观察坐标系；
首先假定平面波斜入射到矩形开口，入射声压pi以θi，
的角度入射到第h个矩形开口，
为声波入射到第h个矩形开口过程中产生的反射声压，而
为在此开口入射端产生的辐射声压，
和
分别为矩形开口入射侧内部声压及其质点振速，
为此开口透射端的辐射声压，
和
为此开口内部透射端口声压及其质点振速；
步骤b、计算任意两个矩形开口的互辐射阻抗
按式(1)计算获得任意两个矩形开口的互辐射阻抗
式(1)中j为虚数单位，ρ0＝1.2kg.m-3kg·m-3为空气密度，c0＝340m·s-1为空
气中的声速，
为第i个矩形开口的特征波数，r为第i个矩形开口与第h个矩形开口的中心距；根据模态正交性，只考虑m＝p，n＝q的情况，
表示第i个矩形开口与第h个矩形开口在z＝0处的互辐射阻抗，当i＝h时，
为第h个矩形开口在声波入射端的自辐射阻抗；
步骤c、计算矩形开口声波入射侧截面处的力
第h个矩形开口入射端辐射声压导致的力
可表达为：
式(2)中，H为矩形开口数量，
为第i个矩形开口入射端产生的辐射声压，S0h为入射侧第h个矩形开口的横截面积，
为第h个矩形开口内第(m，n)阶高阶波所对应的截面模态振型；
按式(3)计算获得第h个矩形开口声波入射侧截面处的力平衡式如下：
式(3)中，
为入射声压，
为声波入射到第h个矩形开口过程中产生的反射声压，
为第h个矩形开口入射侧内部声压；
其中，第h个矩形开口的内部声压
及内部质点振速
的计算公式为：
式(4)和(5)中，
和
为第(m，n)阶高阶波的模态声压和质点振速，
为第h个矩形开口在z轴方向上传播的第(m，n)阶高阶波的波数，
和
分别表示孔隙内传播介质中的特征阻抗和特征波数；当需考虑粘滞阻尼效应时，应代入考虑粘滞阻尼效应的Zc和kc，即：
式(5.1)中，特征阻抗Zc和特征波数kc的表达式为：
式(6)和(7)中，j为虚数单位，S0为入射侧矩形开口的横截面积，
Y′＝jωC(ω)S0 (9)
式(8)和(9)中，j为虚数单位，S0为入射侧矩形开口的横截面积，ω为角频率，
式(10)和(11)中，γ＝1.4为介质的比热比，ρ0＝1.2kg.m-3kg·m-3为空气密度，c0＝340m·s-1为空气中的声速，μ＝1.86×10-5Pa·s为剪切粘度，Cp＝1010J·kg-1·K-1为恒温比热，λ＝0.026W·m-1·K-1为导热系数，此外，对于横截面为矩形的孔隙或者开口，横截面形状影响系数Fr(η)的计算公式为：
在式(12)中，j为虚数单位，ω为角频率；在实际计算时，求和上限s和t可通过收敛性分析获得，且有：
然后把
分别代入式(12)，再把计算结果分别代入式(10)、式(11)，得到对应的复密度ρ(ω)和复压缩系数函数C(ω)；
步骤d、计算矩形开口声波入射侧截面和出射侧截面的声传递矩阵
按式(3)计算获得第h个矩形开口第(m，n)阶高阶波的模态声压
和质点振速
在声波入射端截面和透射端截面的声传递矩阵为：
式中Ah、Bh、Ch、Dh分别为第h个矩形开口声波入射端截面和透射端截面的声传递矩阵的四个元素；Ah、Bh、Ch、Dh中包含矩形开口内声传播介质的特征阻抗Zf和特征波数kf，当需考虑粘滞阻尼效应时，应代入考虑粘滞阻尼效应的Zc和kc；
为声波入射侧截面处第h个矩形开口内侧的声压，
为声波出射侧截面处第h个矩形开口内侧的声压，
为声波出射侧截面第h个矩形开口内的质点振速，
为声波出射侧截面第h个矩形开口内的质点振速；
步骤e、计算矩形开口内部透射端口声压
计算获得第h个矩形开口声波出射侧截面处的力平衡式，
式(16)中Slh为透射端第h个矩形开口的横截面积，
为第i个矩形开口与第h个矩形开口在z＝l处的互辐射阻抗，当i＝h，
为第h个矩形开口在声波透射端的自辐射阻抗，
为在z＝l处，即声波透射端第h个矩形开口内的质点振速；
因为互模态辐射阻抗可以忽略不计，式(9)可简化为：
由上式可得
式(18)中
步骤f、计算矩形开口声传递率
按式(20)计算获得第n个矩形开口的入射声功率，声源的声功率为：
联立式(4)、(5)、(20)，即可求得透射声功率：
计算获得全体以θi，
的角度入射到H个矩形开口的声波耦合声传递率，
为获得式(22)中
的值，需计算多个
通过联立式(2)、(3)、(4)、(15)并将对应的矩形开口的截面尺寸代入，多次重复上述计算即可求得全部的
最终可以按式(23)计算获得考虑粘滞阻尼效应的多矩形开口耦合声传递损失，即：
TL＝-10lg(τ) (23)。
2.根据权利要求1所述的一种周期排列任意数量大矩形开口声传递损失的计算方法，其特征在于，在不考虑热粘性损失时，特征阻抗Zf＝ρ0c0，特征波数kf＝ω/c0，步骤c中第h个矩形开口的内部质点振速
的计算公式为：
式中，
和
为第(m，n)阶高阶波的模态声压和质点振速，
为第h个矩形开口在z轴方向上传播的第(m，n)阶高阶波的波数，
和
分别表示孔隙内传播介质中的特征阻抗和特征波数，ρ0＝1.2kg.m-3为空气密度，ω为角频率；同时，步骤d的式(15)中，Ah、Bh、Ch、Dh中包含矩形开口内声传播介质的特征阻抗Zf和特征波数kf，在不考虑热粘性损失的情形下，在计算中也要代入Zf＝ρ0c0，kf＝ω/c0。
3.根据权利要求1所述的一种周期排列任意数量大矩形开口声传递损失的计算方法，其特征在于，所述矩形开口数量H的取值范围为2至无限大，即式(2)、式(18)和式(22)中的H可以代入2至无限大中的任意整数，进而得到任意数量大矩形开口的耦合声传递损失计算结果。
4.根据权利要求1所述的一种周期排列任意数量大矩形开口声传递损失的计算方法，其特征在于，步骤b的式(1)中，r的值为大于第i个矩形开口与第h个矩形开口宽度之和的2倍的任意值，进而实现任意间距排列多矩形开口耦合声传递损失计算。
5.根据权利要求1所述的一种周期排列任意数量大矩形开口声传递损失的计算方法，其特征在于，第h个矩形开口的入射侧截面长度aih、入射侧截面宽度2bih、出射侧长度2aoh以及出射侧截面尺寸长度2boh均和第i个矩形开口相对应的截面尺寸可以相等或不相等，在将对应的矩形开口的截面尺寸代入式(1)、式(2)、式(3)、式(14)、式(16)、式(17)、式(18)、式(19)、式(20)、式(21)、式(22)中即可实现等尺寸或不等尺寸多矩形开口耦合传递损失计算。
6.根据权利要求1所述的一种周期排列任意数量大矩形开口声传递损失的计算方法，将第h个矩形开口内第(m，n)阶高阶波声能所对应的截面模态振型
代入式(2)、式(3)、式(4)、式(5)、式(16)、式(17)、式(19)中，并在式(4)、式(5)中代入第(m，n)阶高阶波的模态声压
第(m，n)阶高阶波的质点振速
第h个矩形开口在z轴方向上传播的第(m，n)阶高阶波的波数
可得大矩形开口的耦合声传递损失计算结果，所述大矩形开口是指在计算该矩形开口声传递特性时须将开口内部高阶波声能传递纳入考虑范围内的矩形开口。
7.根据权利要求6所述的一种周期排列任意数量大矩形开口声传递损失的计算方法，当m＝n＝0时，可得平面波的矩形开口的声传递损失计算结果。
8.根据权利要求6所述的一种周期排列任意数量大矩形开口声传递损失的计算方法，当m、n为0到无限大的任意正整数时，可得高阶波的大矩形开口的声传递损失计算结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111616789.XA CN114330060B (zh) | 2021-12-27 | 2021-12-27 | 一种周期排列任意数量大矩形开口声传递损失的计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111616789.XA CN114330060B (zh) | 2021-12-27 | 2021-12-27 | 一种周期排列任意数量大矩形开口声传递损失的计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114330060A true CN114330060A (zh) | 2022-04-12 |
| CN114330060B CN114330060B (zh) | 2024-02-13 |
Family
ID=81014288
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111616789.XA Active CN114330060B (zh) | 2021-12-27 | 2021-12-27 | 一种周期排列任意数量大矩形开口声传递损失的计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114330060B (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004125662A (ja) * | 2002-10-03 | 2004-04-22 | Ono Sokki Co Ltd | 透過損失測定方法、透過損失測定装置、および透過損失測定プログラム |
| CN104198584A (zh) * | 2014-09-18 | 2014-12-10 | 合肥工业大学 | 一种获得圆形开口声传递率及声传递损失的方法 |
| CN104198032A (zh) * | 2014-08-14 | 2014-12-10 | 合肥工业大学 | 一种矩形开口声传递率及声传递损失计算方法 |
-
2021
- 2021-12-27 CN CN202111616789.XA patent/CN114330060B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004125662A (ja) * | 2002-10-03 | 2004-04-22 | Ono Sokki Co Ltd | 透過損失測定方法、透過損失測定装置、および透過損失測定プログラム |
| CN104198032A (zh) * | 2014-08-14 | 2014-12-10 | 合肥工业大学 | 一种矩形开口声传递率及声传递损失计算方法 |
| CN104198584A (zh) * | 2014-09-18 | 2014-12-10 | 合肥工业大学 | 一种获得圆形开口声传递率及声传递损失的方法 |
Non-Patent Citations (1)
| Title |
|---|
| 杨亮;季振林;: "阻性管道消声性能预测的边界元数值配点混合方法", 振动工程学报, no. 03, 15 June 2016 (2016-06-15) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114330060B (zh) | 2024-02-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bucciarelli et al. | A multilayer microperforated panel prototype for broadband sound absorption at low frequencies | |
| Yu et al. | Hybrid silencers with micro-perforated panels and internal partitions | |
| Aurégan et al. | Low frequency sound attenuation in a flow duct using a thin slow sound material | |
| CN112071296B (zh) | 隔音装置 | |
| CN108255781B (zh) | 一种高超声速目标表面动态温度建模方法 | |
| WO2018153318A1 (zh) | 一种降噪装置和机柜 | |
| Yu et al. | Modeling vibroacoustic systems involving cascade open cavities and micro-perforated panels | |
| Morrow | Point-to-point correlation of sound pressures in reverberation chambers | |
| CN114330060A (zh) | 一种周期排列任意数量大矩形开口声传递损失的计算方法 | |
| Selamet et al. | Acoustic attenuation performance of circular expansion chambers with single-inlet and double-outlet | |
| Gao et al. | Broadband ventilated sound insulation in a highly sparse acoustic meta-insulator array | |
| Lee et al. | On the acoustic absorption of multi-layer absorbers with different inner structures | |
| Sakamoto et al. | Theoretical analysis of the sound-absorption coefficient and transmission loss for longitudinal clearances among the close-packed cylinders (Three kinds of estimation and experiment) | |
| Jiang et al. | BEM modeling of mufflers with diesel particulate filters and catalytic converters | |
| Farooqui et al. | Validation of low frequency noise attenuation using locally resonant patches | |
| CN110807288A (zh) | 可定制宽频的高效通风吸声器有限元模拟及演示验证方法 | |
| Sakamoto et al. | Theoretical analysis and measurement of sound transmission loss in louver elements with a sound attenuating function using a Helmholtz resonator array | |
| Sakamoto et al. | Study of arrayed side branch silencer for ventilation door (Theoretical analysis and measurement of transmission loss) | |
| CN108417127B (zh) | 一种声学超颖表面实验装置 | |
| Rämmal et al. | Sound reflection at an open end of a circular duct exhausting hot gas | |
| Zhou et al. | Acoustic properties of multiple cavity resonance liner for absorbing higher-order duct modes | |
| Goto et al. | Sound attenuation of periodic micro-perforated chamber mufflers using the spectral transfer matrix method | |
| Yu et al. | Sound attenuation using duct silencers with micro-perforated panel absorbers | |
| Simon | Long elastic open neck acoustic resonator in flow | |
| Zhao et al. | Experimental study of broadband noise attenuation using variable locally reacting impedance |
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 |