CN109783969A - 一种应用传递矩阵计算不同厚度均匀材料声学特性的方法 - Google Patents
一种应用传递矩阵计算不同厚度均匀材料声学特性的方法 Download PDFInfo
- Publication number
- CN109783969A CN109783969A CN201910075965.XA CN201910075965A CN109783969A CN 109783969 A CN109783969 A CN 109783969A CN 201910075965 A CN201910075965 A CN 201910075965A CN 109783969 A CN109783969 A CN 109783969A
- Authority
- CN
- China
- Prior art keywords
- acoustic
- thickness
- matrix
- transmission matrix
- homogeneous material
- 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.)
- Pending
Links
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明涉及一种应用传递矩阵由已知厚度材料推导该材料在不同厚度情况下声学特性的方法,所述方法包括:通过实验测得均匀材料特定厚度下的声传递矩阵,进而计算该材料不同厚度的声传递矩阵,最后计算预测该材料不同厚度下的吸声系数和隔声量。材料的厚度对其声学特性有重要影响,同样的材料,厚度不同也会表现出不同的声学特性。本发明提出的方法可以从实测得到的材料声传递矩阵出发,计算该材料不同厚度下的吸声系数和隔声量,特别适合计算预测泡沫海绵、毡类和纤维类等局部作用材料的声学特性。本发明能根据噪声治理的实际需要,实现材料声学特性的优化设计,避免材料选择和制备的盲目,为噪声污染治理提供有力的技术支撑。
Description
技术领域
本发明涉及不同厚度材料声学特性预测的技术,具体地说涉及一种应用传递矩阵由已知厚度材料推导该材料在不同厚度情况下声学特性的方法。
背景技术
声学材料的吸声和隔声量等声学特性在噪声治理工程中有着重要应用,准确掌握材料的声学特性是对其进行合理有效应用的关键,这就需要一种能够对材料声学特性进行准确测量和预测的方法。
在传统技术中,通常采用混响室法和管测法测量材料声学特性。混响室法中声波的入射方式与实际工程较为接近,但测量结果的复现性差,材料的面积会对测量结果有影响,并且测量代价大,不易操作;管测法可以测量小样本材料的吸声系数和隔声量,测量方便,但只能测量现有样本的声学特性,如果材料的厚度发生改变,其声学特性也会发生改变,需要制作相应厚度的测试样本进行测量。声学材料的厚度对其声学特性有重要影响,材料在不同厚度情况下会表现出不同的声学特性。如果完全通过实测方法获得满足所需声学效果的材料,需要为测量制作不同厚度的材料测试样本,效率极低。
为了高效搞清材料在不同厚度情况下的声学特性,研究不同厚度材料声学特性的计算预测方法很有必要。如果掌握不同厚度材料声学特性的计算预测方法,就能根据噪声治理的实际需要,实现材料声学特性的优化设计,避免材料选择和制备的盲目,为噪声污染治理工作提供有力的技术支撑,减少企业资源消耗,提高企业效益。
发明内容
为了解决背景技术存在的不同厚度均匀材料下声学特性的测量计算问题,本发明提供了一种应用传递矩阵计算不同厚度均匀材料声学特性的方法,能够对不同厚度均匀材料的声学特性进行预测计算。
本发明解决其技术问题所采用的技术方案是:一种应用传递矩阵由已知厚度均匀材料推导该均匀材料在不同厚度情况下声学特性的方法,其特征在于包括以下步骤:
a,获取均匀材料的参考厚度。
其中,均匀材料是指满足均匀性这一特性的材料。材料的均匀性是指在组成材料的空间内各个不同的点处的材料性质是相同的,如密度等。具体地,获取均匀材料的参考厚度,并将该参考厚度作为参考样本,以保证后续计算的准确性。当然,也可以将参考厚度理解为已知厚度。
b,获取所述均匀材料的参考厚度对应的参考声传递矩阵。
其中,声传递矩阵是通过矩阵来描述声音在传递过程中输入与输出之间的关系。均匀材料的参考厚度与声传递矩阵之间存在对应关系,该对应关系可以是一一对应关系,也可以是多对一或者一对多的对应关系。
具体地,在获得均匀材料的参考厚度之后,获取该均匀材料的参考厚度对应的参考声传递矩阵。参考声传递矩阵也可以作为参考样本。
c,获取所述均匀材料的目标厚度。
其中,目标厚度是指所需厚度。
d,根据所述参考厚度、参考声传递矩阵、目标厚度以及下述公式(1),可通过矩阵的幂运算获得所述均匀材料在目标厚度下的目标声传递矩阵,目标厚度与参考厚度的比值为矩阵的幂,
其中,为参考声传递矩阵,d为均匀材料的目标厚度,d0为均匀材料的参考厚度,[Td]为目标声传递矩阵。
具体地,在获得参考厚度、参考声传递矩阵和目标厚度之后,将其代入上述公式(1)中,便可得到均匀材料的目标厚度对应的目标声传递矩阵,即[Td]。
e,根据所述目标声传递矩阵,得到所述均匀材料的目标厚度所对应的声学特性,所述声学特性包括吸声系数和隔声量。
具体地,在获得该均匀材料的目标声传递矩阵[Td]的基础上,可以依据相应的声学特性公式来获得该均匀材料在实际声学边界条件下的吸声系数和隔声量。
其中参考声传递矩阵的获取方式,包括以下步骤:
(a)以均匀材料两个表面的声压P和声速V为声传递矩阵端口参量,设一个表面的位置坐标为x=0,另一表面位置坐标为x=d,d为材料厚度,两表面声压和声速之间可建立如下关系:
通过两种或两种以上不同的端口状态参量可以求出声传递矩阵中的四个参数T11、T12、T21、T22,这四个参数反映特定材料本身固有的声学特性,与材料的声学边界条件无关;
(b)以材料两个表面的前后四个声波幅值A、B、C、D为声传递矩阵端口参量,其中A、B为声源管中入射波和反射波的声压幅值,C、D为接收管中入射波和反射波的声压幅值。四个声波幅值之间可建立如下关系:
通过阻抗管的末端设置两种不同的边界负载条件,可以求出声传递矩阵中的四个参数S11、S12、S21、S22;
通过[S]矩阵和[T]矩阵之间的转换关系,将[S]矩阵转换为[T]矩阵。
本发明的有益效果是:可以通过测量特定厚度均匀材料的声传递矩阵,计算该均匀材料在不同厚度情况下的声传递矩阵,再由声传递矩阵计算均匀材料的声学特性,只需测量材料一种厚度的样本,就能通过计算得到该材料各种厚度下的声学特性,可以根据噪声治理的实际需求指导材料声学特性的设计,使材料声学特性设计成为可能。
附图说明
图1为一个实施例提供的测量材料传递矩阵参数的原理图;
图1中:1、2、3及4为传感器放置位置;5为被测材料;A、B为1、2位置处的声压幅值;C、D为3、4位置处的声压幅值;
图2一个实施例提供的测试系统连接示意图;
图3一个实施例提供的计算原理流程图;
图4一个实施例提供的预测与实测50mm泡沫海绵的吸声系数进行对比;
图4中:预测结果为在刚性壁条件下,通过测量厚度25mm泡沫海绵的声传递矩阵,应用本发明的方法计算预测厚度50mm泡沫海绵的吸声系数;
图5一个实施例提供的预测结果与实测50mm泡沫海绵的隔声量进行对比;
图5中:预测结果为在刚性壁条件下,通过测量厚度25mm泡沫海绵的声传递矩阵,应用本发明的方法计算预测厚度50mm泡沫海绵的隔声量。
具体实施方式
下面参照附图和具体实施例对本发明做进一步详细说明。本发明可以用许多不同的形式来实施,并不局限于此处描述的实施例,提供这些实施例是为了详尽地且完全的公开本发明,并且向所属技术领域的技术人员充分传达本发明的范围。
本申请一实施例提供了一种应用传递矩阵计算不同厚度均匀材料声学特性的方法,该方法包括以下步骤:
a,获取均匀材料的参考厚度。
其中,均匀材料是指满足均匀性这一特性的材料。材料的均匀性是指在组成材料的空间内各个不同的点处的材料性质是相同的,如密度等。具体地,获取均匀材料的参考厚度,并将该参考厚度作为参考样本,以保证后续计算的准确性。当然,也可以将参考厚度理解为已知厚度。
b,获取所述均匀材料的参考厚度对应的参考声传递矩阵。
其中,声传递矩阵是通过矩阵来描述声音在传递过程中输入与输出之间的关系。均匀材料的参考厚度与声传递矩阵之间存在对应关系,该对应关系可以是一一对应关系,也可以是多对一或者一对多的对应关系。
具体地,在获得均匀材料的参考厚度之后,获取该均匀材料的参考厚度对应的参考声传递矩阵。参考声传递矩阵也可以作为参考样本。
c,获取所述均匀材料的目标厚度。
其中,目标厚度是指所需厚度。
d,根据所述参考厚度、参考声传递矩阵、目标厚度以及下述公式(1),可通过矩阵的幂运算获得所述均匀材料在目标厚度下的目标声传递矩阵,目标厚度与参考厚度的比值为矩阵的幂,
其中,为参考声传递矩阵,d为均匀材料的目标厚度,d0为均匀材料的参考厚度,[Td]为目标声传递矩阵。
具体地,在获得参考厚度、参考声传递矩阵和目标厚度之后,将其代入上述公式(1)中,便可得到均匀材料的目标厚度对应的目标声传递矩阵,即[Td]。
e,根据所述目标声传递矩阵,得到所述均匀材料的目标厚度所对应的声学特性,所述声学特性包括吸声系数和隔声量。
具体地,在获得该均匀材料的目标声传递矩阵[Td]的基础上,可以依据相应的声学特性公式来获得该均匀材料在实际声学边界条件下的吸声系数和隔声量。
其中参考声传递矩阵的获取方式,包括以下步骤:
1、使用阻抗管进行测量,将待测均匀材料样本加工为符合阻抗管截面尺寸的形状,按照图1中的方式将待测样本置于上、下游中间测量管段。
2、按照图2所示的测量系统在上、下游处各安置2个,总共4个幅值和相位匹配的传感器,测量过程中改变阻抗管管道末端的声阻抗,测量两种末端状态下的上、下游管段声压。实测获取特定厚度均匀材料声传递矩阵,测量有以下两种方法,实际测量时选择其一即可。
(a)通过传感器测量得到上、下游的声压,获得4个测点之间的传递函数。以均匀材料两个表面的声压P和声速V为声传递矩阵端口参量,设一个表面的位置坐标为x=0,另一表面位置坐标为x=d,d为材料厚度,两表面声压和声速之间可建立如下关系:
因为均匀材料声学特性在测量过程中保持不变,均匀材料任一表面端口参量声压P和声速V的改变都会导致另一表面端口参量的改变,通过两种或两种以上不同的端口状态参量可以求出声传递矩阵中的四个参数T11、T12、T21、T22,这四个参数反映特定材料本身固有的声学特性,与材料的声学边界条件无关;
(b)以材料两个表面的前后四个声波幅值A、B、C、D为声传递矩阵端口参量,其中A、B为声源管中入射波和反射波的声压幅值,C、D为接收管中入射波和反射波的声压幅值。四个声波幅值之间可建立如下关系:
通过阻抗管的末端设置两种不同的边界负载条件,可以求出声传递矩阵中的四个参数S11、S12、S21、S22;
3、通过[S]矩阵和[T]矩阵之间的转换关系,将[S]矩阵转换为[T]矩阵。
通过具体数值来对本申请的发明构思做出进一步的说明。具体而言,本实施例按照图3所示计算原理流程,通过对25mm厚的泡沫海绵进行吸声系数与隔声量的测量与计算,从而预测50mm厚的泡沫海绵的声学特性。
在本实施例中,提供了一种声学特性计算装置,该装置使用阻抗管进行上、下游通道组装,传声器1与2、3与4的位置间距皆为0.1m,2与3的位置间距为0.35m。
进一步地,基于上述装置,应用声传递矩阵计算不同厚度均匀材料声学特性的方法,包括以下步骤:
步骤1:使用阻抗管测量,将待测样本加工为符合阻抗管尺寸的形状,置于上、下游中间管段。
步骤2:在上、下游处各放置两个幅值和相位匹配过的传感器,测量材料上、下游处的声压,改变阻抗管管道末端的声阻抗,测量材料在两种末端状态下的上、下游声压。
步骤3:由1、2位置及3、4位置测量的声压结果,分别得到系数A、B、C、及D的表达式如下:
其中是波数,c0是空气中声速。
均匀材料两个表面的前后四个声波幅值A、B、C、D之间可建立如下关系:
因为测量过程中均匀材料声学特性保持不变,在阻抗管末端设置两种不同的声学边界,求出声传递矩阵中的四个参数S11、S12、S21、S22。
由系数A到D可获得材料两表面的声压与声速信息:
声音穿过被测均匀材料时,由均匀材料本身的性质,前后表面上的声压和声速存在下述关系:
其中d为均匀材料本身厚度,通过两种末端状态下的测量结果,可解得反映均匀材料本身固有特性的参数T11、T12、T21、T22,即获得均匀材料的传递矩阵。
步骤4:建立[T]矩阵和[S]矩阵的对应关系:
步骤5:通过测量得到厚度为d0均匀材料的传递矩阵设厚度为d相同均匀材料的声传递矩阵为[Td],则有:
步骤6:已知厚度为d均匀材料的传递矩阵为[T],从[T]出发,可计算出该厚度均匀材料的吸声系数、隔声量等声学特性参数。
不同厚度泡沫海绵前后表面上的声压、声速存在如下关系:
均匀材料与刚性壁接触的表面处介质质点速度为0,所以:
P0=T11Px,
V0=T21Px。
均匀材料前表面的声阻抗率求解公式:
通过以上各式可得均匀材料的吸声系数:
α=1-|rp|2
从图4预测与实测50mm泡沫海绵的吸声系数对比中可以看出预测与实测曲线有较好的吻合度。
声压透射系数:
隔声量计算公式:
从图5预测与实测50mm泡沫海绵的隔声系数对比中可以看出预测与实测曲线有较好的吻合度。测量材料一种厚度的样本,通过计算得到该材料各种厚度下的声学特性具有较好的准确性。
Claims (2)
1.一种应用传递矩阵计算不同厚度均匀材料声学特性的方法,其特征在于,包括以下步骤:
a,获取均匀材料的参考厚度;
b,获取所述参考厚度的均匀材料对应的参考声传递矩阵;
c,获取所述均匀材料的目标厚度;
d,根据所述参考厚度、参考声传递矩阵、目标厚度以及下述公式(1),可通过矩阵的幂运算获得所述均匀材料在目标厚度下的目标声传递矩阵,目标厚度与参考厚度的比值为矩阵的幂:
其中,为参考声传递矩阵,d为均匀材料的目标厚度,d0为均匀材料的参考厚度,[Td]为目标声传递矩阵;
e,根据所述目标声传递矩阵,得到所述均匀材料在实际声学边界条件下所对应的声学特性,所述声学特性包括吸声系数和隔声量。
2.根据权利要求1所述的一种应用传递矩阵计算不同厚度均匀材料声学特性的方法,其特征在于,所述参考声传递矩阵的获取方式,包括以下步骤:
(a)以均匀材料两个表面的声压P和声速V为声传递矩阵端口参量,设一个表面的位置坐标为x=0,另一表面位置坐标为x=d,d为材料厚度,两表面声压和声速之间可建立如下关系:
通过两种或两种以上不同的端口状态参量可以求出声传递矩阵中的四个参数T11、T12、T21、T22,这四个参数反映特定材料本身固有的声学特性,与材料的声学边界条件无关;
(b)以材料两个表面的前后四个声波幅值A、B、C、D为声传递矩阵端口参量,其中A、B为声源管中入射波和反射波的声压幅值,C、D为接收管中入射波和反射波的声压幅值。四个声波幅值之间可建立如下关系:
通过阻抗管的末端设置两种不同的边界负载条件,可以求出声传递矩阵中的四个参数S11、S12、S21、S22;
通过[S]矩阵和[T]矩阵之间的转换关系,将[S]矩阵转换为[T]矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910075965.XA CN109783969A (zh) | 2019-01-26 | 2019-01-26 | 一种应用传递矩阵计算不同厚度均匀材料声学特性的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910075965.XA CN109783969A (zh) | 2019-01-26 | 2019-01-26 | 一种应用传递矩阵计算不同厚度均匀材料声学特性的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109783969A true CN109783969A (zh) | 2019-05-21 |
Family
ID=66501404
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910075965.XA Pending CN109783969A (zh) | 2019-01-26 | 2019-01-26 | 一种应用传递矩阵计算不同厚度均匀材料声学特性的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109783969A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111007152A (zh) * | 2019-12-24 | 2020-04-14 | 无锡吉兴汽车声学部件科技有限公司 | 一种声学性能综合评估方法及系统 |
CN114371224A (zh) * | 2021-12-03 | 2022-04-19 | 华南理工大学 | 一种aba结构随机入射吸声系数的测试方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1742048A1 (en) * | 2005-06-08 | 2007-01-10 | Rieter Technologies AG | Acoustic device for determining parameters of a porous material |
CN103091401A (zh) * | 2013-01-15 | 2013-05-08 | 西北工业大学 | 一种测量声学材料特征阻抗和复波数的装置和测量方法 |
CN103115966A (zh) * | 2013-01-15 | 2013-05-22 | 西北工业大学 | 一种采用脉冲法测量声学材料隔声量的装置和测量方法 |
CN104198032A (zh) * | 2014-08-14 | 2014-12-10 | 合肥工业大学 | 一种矩形开口声传递率及声传递损失计算方法 |
CN104535647A (zh) * | 2014-11-26 | 2015-04-22 | 华晨汽车集团控股有限公司 | 一种多层材料吸隔声性能预测装置及方法 |
-
2019
- 2019-01-26 CN CN201910075965.XA patent/CN109783969A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1742048A1 (en) * | 2005-06-08 | 2007-01-10 | Rieter Technologies AG | Acoustic device for determining parameters of a porous material |
CN103091401A (zh) * | 2013-01-15 | 2013-05-08 | 西北工业大学 | 一种测量声学材料特征阻抗和复波数的装置和测量方法 |
CN103115966A (zh) * | 2013-01-15 | 2013-05-22 | 西北工业大学 | 一种采用脉冲法测量声学材料隔声量的装置和测量方法 |
CN104198032A (zh) * | 2014-08-14 | 2014-12-10 | 合肥工业大学 | 一种矩形开口声传递率及声传递损失计算方法 |
CN104535647A (zh) * | 2014-11-26 | 2015-04-22 | 华晨汽车集团控股有限公司 | 一种多层材料吸隔声性能预测装置及方法 |
Non-Patent Citations (2)
Title |
---|
冯涛 等: ""用传递矩阵法计算多层泡沫材料的法向吸声系数"", 《机械设计与制造》 * |
杨梦露 等: ""空气背衬条件下泡沫材料法向吸声系数混合计算方法研究"", 《噪声与振动控制》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111007152A (zh) * | 2019-12-24 | 2020-04-14 | 无锡吉兴汽车声学部件科技有限公司 | 一种声学性能综合评估方法及系统 |
CN114371224A (zh) * | 2021-12-03 | 2022-04-19 | 华南理工大学 | 一种aba结构随机入射吸声系数的测试方法 |
CN114371224B (zh) * | 2021-12-03 | 2023-08-22 | 华南理工大学 | 一种aba结构随机入射吸声系数的测试方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2020101234A4 (en) | An Optimization Method for Excitation Parameters of Ultrasonic Infrared Thermography Crack Nondestructive Testing | |
Howard et al. | Detectability of corrosion damage with circumferential guided waves in reflection and transmission | |
CN105651215B (zh) | 一种超声声速未知条件下的涂层厚度测量方法 | |
Zhang et al. | Efficient immersion imaging of components with nonplanar surfaces | |
Yuan et al. | Phased array guided wave propagation in curved plates | |
CN106770664A (zh) | 一种基于全聚焦成像算法改进边缘缺陷检测的方法 | |
CN106680376A (zh) | 一种基于三维声强阵列的隔声测量系统与测量方法 | |
Jones et al. | Effects of flow profile on educed acoustic liner impedance | |
CN106767580B (zh) | 一种确定复合材料层压结构中缺陷铺层深度的超声方法 | |
Benstock et al. | The influence of surface roughness on ultrasonic thickness measurements | |
CN109696480B (zh) | 一种基于改进时间反转算法的玻璃纤维复合材料声发射源定位成像方法 | |
CN104297346A (zh) | 超声平板导波的金属板材无损检测系统及其检测方法 | |
JP2009063372A (ja) | 空中超音波探傷装置及びその方法 | |
CN109783969A (zh) | 一种应用传递矩阵计算不同厚度均匀材料声学特性的方法 | |
CN105353043A (zh) | 基于abaqus的金属薄板微裂纹时间反转定位方法 | |
Mu et al. | Defect circumferential sizing by using long range ultrasonic guided wave focusing techniques in pipe | |
Nocke et al. | Impedance deduction from broad-band, point-source measurements at grazing incidence | |
CN110208383A (zh) | 一种基于反转路径差信号的板结构兰姆波无参考成像方法 | |
CN112432998B (zh) | 一种带有声腔结构的橡胶板粘接缺陷超声波无损检测方法 | |
CN206563731U (zh) | 一种基于三维声强阵列的隔声测量系统 | |
Bezdek et al. | A coupled finite-element, boundary-integral method for simulating ultrasonic flowmeters | |
CN105548356A (zh) | 基于t模态导波检测带环焊缝小弯曲半径弯管缺陷的方法 | |
CN102867099A (zh) | 一种矩形超声换能器声场仿真方法 | |
CN104793056B (zh) | 一种测量水声压电换能器在非消声水池中平均辐射阻抗的方法 | |
CN113219053B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190521 |