CN112951188B - 一种有源式微穿孔板吸声器及提升其低频吸声性能的方法 - Google Patents
一种有源式微穿孔板吸声器及提升其低频吸声性能的方法 Download PDFInfo
- Publication number
- CN112951188B CN112951188B CN202110115539.1A CN202110115539A CN112951188B CN 112951188 B CN112951188 B CN 112951188B CN 202110115539 A CN202110115539 A CN 202110115539A CN 112951188 B CN112951188 B CN 112951188B
- Authority
- CN
- China
- Prior art keywords
- sound
- microperforated panel
- active
- panel
- absorber
- 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 81
- 239000006096 absorbing agent Substances 0.000 title claims abstract description 80
- 238000010521 absorption reaction Methods 0.000 title claims abstract description 70
- 230000014509 gene expression Effects 0.000 claims description 84
- 238000006073 displacement reaction Methods 0.000 claims description 52
- 230000008569 process Effects 0.000 claims description 26
- 230000005284 excitation Effects 0.000 claims description 21
- 239000011159 matrix material Substances 0.000 claims description 15
- 230000004044 response Effects 0.000 claims description 13
- 238000013016 damping Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 11
- 239000000463 material Substances 0.000 claims description 10
- 239000002245 particle Substances 0.000 claims description 9
- 239000011148 porous material Substances 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 229910052782 aluminium Inorganic materials 0.000 claims description 4
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 claims description 4
- 239000013598 vector Substances 0.000 claims description 4
- 230000002068 genetic effect Effects 0.000 claims description 3
- 238000010845 search algorithm Methods 0.000 claims description 3
- 101100184723 Homo sapiens PMPCA gene Proteins 0.000 abstract description 73
- 102100025321 Mitochondrial-processing peptidase subunit alpha Human genes 0.000 abstract description 73
- 230000001276 controlling effect Effects 0.000 abstract description 3
- 230000001105 regulatory effect Effects 0.000 abstract description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 239000011358 absorbing material Substances 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000002657 fibrous material Substances 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 239000012774 insulation material Substances 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/16—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/162—Selection of materials
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/16—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/172—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using resonance effects
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K2210/00—Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
- G10K2210/30—Means
- G10K2210/301—Computational
- G10K2210/3036—Modes, e.g. vibrational or spatial modes
-
- 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
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
Abstract
本发明提供了一种有源式微穿孔板吸声器及提升其低频吸声性能的方法，解决采用“有源”方式提高MPPA低频吸声性能时存在难以获得理想的活塞式控制面源、且一维声场条件也过于简化的弊端，导致难以进行广泛工程应用的问题。本发明将传统微穿孔板吸声器的刚性空腔后壁用点力控制的弹性平板代替形成有源式微穿孔板吸声器，通过点力控制弹性平板的结构振动来调控有源式微穿孔板吸声器腔内的声场分布，进而调整有源式微穿孔板吸声器表面的阻抗特性，使其与空气介质的特征阻抗相匹配，进而达到显著提升微穿孔板吸声器低频吸声性能的目的。
Description
技术领域
本发明属于噪声控制技术领域，具体涉及一种有源式微穿孔板吸声器及提升其低频吸声性能的方法。
背景技术
随着社会发展及生活水平不断提高，噪声污染作为一项重要的环境问题越来越受到人们的重视，噪声控制对人对物均是一项非常重要的工作。传统的噪声控制方式是通过添加吸声与隔声材料，通过吸收并耗散噪声能量而达到控制噪声的目的。微穿孔板是该类材料中的典型代表，微穿孔板可由各种基材(如金属、玻璃、木材等)穿微孔而制成，更加环保、经久耐用，适用于环境苛刻的特殊场合；因此，微穿孔板吸声材料具有广阔的应用前景。
微穿孔板吸声器(Micro-perforated panel absorber,简称MPPA)是微穿孔吸声材料在工程中的典型应用；MPPA由微穿孔板及其后部空腔组合而成。微穿孔板具有较大的声阻及较小的声抗，因而具有良好的吸声性能，由后部空腔产生的共鸣吸声效应，显著增大了微孔中空气介质与板结构的相对速度，使得MPPA在较宽的频带内具有优越的吸声性能。然而，对于特性参数的MPPA，其有效的吸声频带(吸声系数大于0.5的频带)只能限定在特定的范围内。
为了有效拓宽MPPA的吸声频带，提出了许多改进设计措施，例如：将多个微穿孔板平行的插入MPPA的空腔内部形成多层MPPA，或者将具有不同腔深的多个MPPA并行排布形成大面积MPPA阵列，此类方法将不同MPPA的有效吸声频带相互叠加，显著拓宽了结构总的吸声频带；亦或，通过改进MPPA的结构，实现单个MPPA有效吸声频带的拓宽，如用不规则的空腔代替传统矩形空腔来增强MPPA内的声-振耦合效应，来拓宽有效吸声频带；或空腔内填充各向异性的纤维材料、减小微穿孔板的孔隙直径与板厚、将微穿孔板与声学超材料及Helmholtz共鸣器结合等方法均能显著拓宽MPPA的有效吸声频带。
上述方法均能有效拓宽MPPA的中高频吸声性能，但对低频吸声性能的改善非常有限，而这在工程实际中却有很迫切的应用需求。
MPPA的空腔越深，有效吸声频带越趋于低频段，但太深的空腔深度会导致MPPA体积庞大，使其无法工程应用。因而，有学者提出了“无源”和“有源”两类方法对此进行改进设计。其中，“无源”方法主要是在MPPA内部引入一些“共鸣结构”来拓宽低频吸声频带，如Helmholtz共鸣器、分流扬声器及机械阻抗板等，该方法对低频段的拓宽非常有限。“有源”方法则是通过在MPPA内部引入有源控制技术来大范围拓展低频吸声频带，是拓宽MPPA低频吸声频带的有效方法。
然而，现有“有源”方法仅针对小尺寸MPPA、且在一维声场环境下，提出将MPPA的空腔后部用活塞式可变阻抗面来代替，通过控制空腔声场来提升MPPA的低频吸声性能；该方法在设定MPPA空腔内部为平面声场的条件下，提出了“声压释放”与“阻抗匹配”两种策略来提升低频吸声性能。由于实际中难以获得理想的活塞式控制面源，且一维声场条件也过于简化。这些弊端使得该技术目前难以得到广泛的工程应用，严重阻碍了该技术的发展。
发明内容
本发明的目的在于解决采用“有源”方式提高MPPA低频吸声性能时存在难以获得理想的活塞式控制面源、且一维声场条件也过于简化的弊端，导致难以进行广泛工程应用的问题，而提供一种有源式微穿孔板吸声器及提升其低频吸声性能的方法。
为实现上述目的，本发明所提供的技术解决方案是：
一种有源式微穿孔板吸声器，其特殊之处在于，由微穿孔板及其后部空腔组合而成，所述后部空腔的后壁为点力控制的弹性平板。
同时，本发明还提供了一种提升上述有源式微穿孔板吸声器低频吸声性能的方法，其特殊之处在于，包括以下步骤：
1)有源式微穿孔板吸声器空腔内声场的声压表达
根据有源式微穿孔板吸声器的几何模型参数、材料参数和初级激励与次级激励参数，对有源式微穿孔板吸声器空腔内的声场进行建模计算，推导出有源式微穿孔板吸声器空腔内微穿孔板和弹性平板表面的声压分布表达式；
2)有源式微穿孔板吸声器入射侧声场的声压表达
根据有源式微穿孔板吸声器入射侧界面上满足的速度连续边界条件，推导入射侧声场中反射声波的声压表达式，从而获得入射侧声场的总声压表达式；
3)有源式微穿孔板吸声器中微穿孔板和弹性平板振动位移的求解
根据模态叠加原理及微穿孔板与弹性平板振动位移满足的波动方程，推导出微穿孔板与弹性平板位移模态幅值满足的矩阵方程，联立两个矩阵，获得微穿孔板与弹性平板位移模态幅值的表达式；
4)最优控制力的求解及控制后吸声系数的计算
将有源式微穿孔板吸声器的吸声系数作为最优的控制目标，通过优化搜索算法获得使吸声系数最大的最优次级控制力幅值，进而可以获得控制后有源式微穿孔板吸声器系统总的振动响应及控制后的吸声系数。
进一步地，步骤1)中，所述有源式微穿孔板吸声器的几何模型参数包括微穿孔板的长度a、宽度b、厚度hM，弹性平板的厚度h2，空腔深度D，微穿孔板的孔隙直径d以及孔隙率σ；
所述有源式微穿孔板吸声器的材料参数包括微穿孔板的密度ρM、弹性平板的密度ρ2、微穿孔板的杨氏模量EM、弹性平板的杨氏模量E2、微穿孔板的泊松比σM、弹性平板的泊松比σ2、微穿孔板的模态阻尼比ξM,mn、弹性平板的模态阻尼比ξ2,mn、空气的密度ρ0、空气的声速c0；
所述有源式微穿孔板吸声器的初级激励与次级激励参数包括初级激励为垂直入射的平面波时，平面波的幅值为p0；次级控制点力的布置位置为(xs,ys)＝(0.05,0.05)时，控制力的幅值Fs。
进一步地，步骤1)的具体步骤为：
1.1)有源式微穿孔板吸声器空腔内声场响应求解
在初级垂直入射平面波及次级控制力的激励下，有源MPPA空腔内产生声场声压。根据空腔内声速度势Φ(x,y,z,t)满足的波动方程，同时利用模态叠加原理，以及空腔内六个壁面上满足的速度连续的边界条件，可推导获得有源MPPA空腔内声速度势的表达式；
具体是：在初级垂直入射平面波及次级控制力的激励下，有源式微穿孔板吸声器空腔内产生声场声压；根据空腔内声速度势Φ(x,y,z,t)满足的齐次波动方程，空腔内四个刚性侧壁处，存在空气质点振速为0的边界条件，利用模态叠加原理，结合格林第二公式以及模态函数的正交性，推导获得有源式微穿孔板吸声器空腔内声速度势的表达式为：
式中，ω为角频率，U与W为x与y方向上声模态序数的上限值，(u,w)为声模态序数，c0为空气的声速；
其中，系数变量Luw与Nuw取决于空腔内声速度势在微穿孔板与弹性平板两侧满足的如下速度连续的边界条件：
式中，D为空腔深度，v2(x,y,t)为弹性平板的质点振速，v1(x,y,t)为临近微穿孔板表面的质点振速，取微穿孔板振速和微孔中空气介质振速的空间平均值，具体表示为：
式中，vM(x,y)为微穿孔板的振速，P为入射侧的总声压，PD为腔内微穿孔板表面声压，σ为微穿孔板的孔隙率，z0为微穿孔板的孔隙阻抗；
根据边界条件(2)与(3)，结合空腔声模态函数的正交性获得声速度势表达式中的系数变量Luw与Nuw的表达式如下：
式中，声模态函数
系数
将式(5)与式(6)带入式(1)中，获得有源式微穿孔板吸声器空腔内任意点声速度势的表达式；
1.2)有源式微穿孔板吸声器腔内微穿孔板表面的声压表达
根据空腔内声速度势与声压的关系，同时利用声模态函数的正交性，推导获得有源MPPA腔内靠近微穿孔板表面的声压表达式；
具体是：根据空腔内声速度势与声压的关系式获得腔内微穿孔板表面声压PD的表达式为：
式中，
由于式(7)的右端仍然含有自身PD，因此，将式(7)左右两端同乘以声模态函数Ψuw(x,y)，并利用模态函数的正交性，获得如下等式：
将式(8)带入式(7)中，获得腔内微穿孔板表面声压PD的表达式：
式中，
1.3)有源式微穿孔板吸声器腔内弹性平板表面的声压表达
根据空腔内声速度势与声压的关系，同时利用声模态函数的正交性，推导获得有源MPPA腔内靠近弹性平板表面的声压表达式；
具体是：根据空腔内声速度势与声压的关系式并利用式(8)，推导获得有源式微穿孔板吸声器腔内弹性平板表面声压P0的表达式：
式中，
进一步地，步骤2)的具体步骤为：
2.1)有源式微穿孔板吸声器入射侧反射声波的声速度势求解
根据有源MPPA入射侧声场在有源MPPA表面满足的速度连续的边界条件，推导获得反射波声速度势的表达式；
具体是：根据反射声波在入射侧管道内满足的齐次波动方程，结合四边刚性壁、入口处为开口的边界条件，获得反射声波声速度势Φr(x,y,z,ω)的表达式为：
式中，系数变量Ruw取决于入射侧表面满足的如下速度连续的边界条件：
式中，Φi为入射平面波的声速度势，k0＝ω/c0为声波数；
根据式(12)，结合声模态函数的正交性，获得系数变量Ruw的表达式：
将式(13)带入式(11)中，获得反射声波的声速度势的表达式；
2.2)有源式微穿孔板吸声器入射侧反射声波的声压表达
根据入射侧反射声波声速度势与反射声波声压的关系，同时利用声模态函数的正交性，推导获得反射波声压的表达式；
具体是：根据入射侧反射声波声速度势与反射声波声压的关系，获得入射侧反射声波的声压表达式：
式中，
式(14)中入射侧总声压P＝Pi+Pr；
式(14)两边同乘以模态函数Ψuw(x,y)，并利用模态函数的正交性，获得如下的表达式：
将式(15)带入式(14)中，获得反射声波声压的表达式：
式中，
其中，Pi为入射波在有源式微穿孔板吸声器表面处的值，
进一步地，步骤3)的具体步骤为：
3.1)微穿孔板振动响应建模
根据微穿孔板振动位移满足的波动方程，同时利用模态叠加原理及模态函数的正交性，推导获得微穿孔板的位移模态幅值Amn满足的方程式：
式中，系数
ωmn与ξmn为第(m,n)阶模态的共振频率与模态阻尼比，
ρM与hM分别为微穿孔板的密度与厚度，
为微穿孔板的模态振型函数，
其中，系数
(m,n)为振动模态的模态序数；
将式(17)中微穿孔板和弹性平板的位移vM(x,y)与v2(x,y)进行模态展开，经化解可获得微穿孔板与弹性平板位移模态幅值Amn与Bmn满足的方程：
定义如下变量：
则M×N个未知的微穿孔板位移模态幅值与弹性平板的位移模态幅值满足如下矩阵方程：
C1A+D1B＝V (19)
式中，A＝[A11,L,Amn,L,AMN]T，B＝[B11,L,Bmn,L,BMN]T，V＝[V11,L,Vmn,L,VMN]T；
3.2)弹性平板振动响应建模
根据弹性平板振动位移满足的波动方程，同时利用模态叠加原理及模态函数的正交性，推导获得微穿孔板与弹性平板的位移模态幅值Amn与Bmn满足如下的方程：
式中，
ω2,mn与ξ2,mn为弹性平板第(m,n)阶模态的共振频率与模态阻尼比，
ρ2与h2为弹性平板的密度与厚度；
Qs,mn为广义次级模态力，
定义如下的变量：
则M×N个未知的微穿孔板与弹性平板位移模态幅值满足如下矩阵方程：
C2A+D2B＝V2+FsQs (21)
式中，V2＝[V2,11,L,V2,mn,L,V2,MN]T；
3.3)两种板的位移模态幅值求解
将微穿孔板与弹性平板位移模态幅值满足的矩阵方程(20)与(21)联立，获得如下的矩阵方程组：
进而获得微穿孔板与弹性平板的位移模态幅值为：
进一步地，步骤4)的具体步骤为：
4.1)最优控制力的求解
将有源式微穿孔板吸声器的吸声系数作为最优的控制目标，表示为：
式中，Πinc与Πabs分别为入射的总声能及有源式微穿孔板吸声器吸收的声能，SMPP为微穿孔板的面积；
采用遗传算法求得使吸声系数α最大时的最优控制力幅值Fs,opt；
4.2)控制后有源式微穿孔板吸声器吸声系数的计算
获得最优控制力幅值Fs,opt后，带入式(23)便可获得最优控制下的位移模态幅值矢量A与B；然后根据吸声系数的计算公式(24)，获得控制后有源MPPA的吸声系数。
进一步地，所述微穿孔板和弹性平板的材质均为铝。
本发明的优点是：
1.本发明将传统微穿孔板吸声器的刚性空腔后壁用点力控制的弹性平板代替形成有源式微穿孔板吸声器，通过点力控制弹性平板的结构振动来调控有源式微穿孔板吸声器腔内的声场分布，进而调整有源式微穿孔板吸声器表面的阻抗特性，使其与空气介质的特征阻抗相匹配，进而达到显著提升微穿孔板吸声器低频吸声性能的目的。
2.本发明方法将有源控制技术引入传统的微穿孔板吸声器中形成有源式微穿孔板吸声器(也称“有源MPPA”)，不仅可以在更宽的低频范围内提升MPPA的吸声性能，同时该结构用点力控制的弹性平板作为次级控制源、且无腔内一维平面声场的限制条件，使得该方法构建的“有源MPPA”更易于工程实现。
附图说明
图1为有源式微穿孔板吸声器模型示意图；
图2为有源式微穿孔板吸声器微穿孔板尺寸为0.2m×0.2m时控制前后吸声系数对比图；
图3为有源式微穿孔板吸声器微穿孔板尺寸为0.2m×0.3m时控制前后吸声系数对比图。
具体实施方式
以下结合附图和具体实施例对本发明的内容作进一步的详细描述：
一种基于结构振动声辐射控制构建的提升有源式微穿孔板吸声器低频吸声性能的方法，其中有源式微穿孔板吸声器(以下统称有源MPPA)由微穿孔板及其后部空腔组合而成，后部空腔的后壁为点力控制的弹性平板，且微穿孔板和弹性平板的材质均为铝，该方法具体包括以下步骤：
1)有源MPPA空腔内声场的声压表达
根据有源MPPA的几何模型参数、材料参数和初级激励与次级激励参数，对有源MPPA空腔内的声场进行建模计算，推导出有源MPPA空腔内微穿孔板和弹性平板表面的声压分布表达式；
上述几何模型参数包括微穿孔板的长度a、宽度b、厚度hM，弹性平板的厚度h2，空腔深度D，微穿孔板的孔隙直径d以及孔隙率σ；
上述材料参数包括微穿孔板的密度ρM、弹性平板的密度ρ2、微穿孔板的杨氏模量EM、弹性平板的杨氏模量E2、微穿孔板的泊松比σM、弹性平板的泊松比σ2、微穿孔板的模态阻尼比ξM,mn、弹性平板的模态阻尼比ξ2,mn、空气的密度ρ0、空气的声速c0；
上述初级激励与次级激励参数包括初级激励为垂直入射的平面波时，平面波的幅值为p0；次级控制点力的布置位置为(xs,ys)＝(0.05,0.05)时，控制力的幅值Fs；
该步骤的过程是：首先，建立有源MPPA空腔中声速度势满足的波动方程，根据空腔各表面处满足的速度连续边界条件推导出空腔内声速度势的表达式；其次，根据空腔内声速度势与声压的关系式，同时利用声模态函数的正交性，推导出空腔内微穿孔板表面及弹性平板表面的声压表达式。具体过程如下：
1.1)有源MPPA空腔内声场响应求解
如图1所示，在垂直入射平面波激励下，部分声能量被有源MPPA吸收进入空腔内，同时，腔内声场也受到点力控制的弹性平板辐射声压的调控作用。声能量的吸收量取决于有源MPPA表面的阻抗与空气特征阻抗的匹配程度。空腔内的声速度势Φ(x,y,z,t)满足齐次波动方程，在腔内四个刚性侧壁处，存在空气质点振速为0的边界条件。基于此，根据模态叠加原理，同时结合格林第二公式以及模态函数的正交性可得腔内声速度势的表达式为：
式中，ω为角频率，U与W为x与y方向上声模态序数的上限值，(u,w)为声模态序数，c0为空气的声速；
其中，未知的系数变量Luw与Nuw取决于空腔内声速度势在微穿孔板与弹性平板两侧满足的如下速度连续的边界条件：
式中，D为空腔深度，v2(x,y,t)为弹性平板的质点振速，v1(x,y,t)为临近微穿孔板表面的质点振速，可认为是微穿孔板振速和微孔中空气介质振速的空间平均值，具体表示为：
式中，vM(x,y)为微穿孔板的振速，P为入射侧的总声压，PD为腔内微穿孔板表面声压，σ为微穿孔板的孔隙率，z0为微穿孔板的孔隙阻抗；
根据边界条件(2)与(3)，结合空腔声模态函数的正交性获得声速度势表达式中的系数变量Luw与Nuw的表达式如下：
式中，声模态函数
系数
将式(5)与式(6)带入式(1)中，获得有源MPPA空腔内任意点声速度势的表达式；
1.2)有源MPPA腔内微穿孔板表面的声压表达
根据空腔内声速度势与声压的关系式获得腔内微穿孔板表面(z＝-D)声压PD的表达式为：
式中，
由于式(7)的右端仍然含有自身PD，因此，将式(7)左右两端同乘以声模态函数Ψuw(x,y)，并利用模态函数的正交性，获得如下等式：
将式(8)带入式(7)中，获得腔内微穿孔板表面声压PD的表达式：
式中，
1.3)有源MPPA腔内靠近弹性平板表面的声压表达
根据空腔内声速度势与声压的关系式并利用式(8)，推导获得有源MPPA腔内弹性平板表面(z＝0)声压P0的表达式：
式中，
2)有源MPPA入射侧声场的声压表达
根据有源MPPA入射侧界面上满足的速度连续边界条件，推导入射侧声场中反射声波的声压表达式，从而获得入射侧声场的总声压表达式；
2.1)有源MPPA入射侧反射声波的声速度势求解
根据反射声波在入射侧管道内满足的齐次波动方程，结合四边刚性壁、入口处为开口的边界条件，获得反射声波声速度势Φr(x,y,z,ω)的表达式为：
式中，未知系数变量Ruw取决于入射侧表面满足的如下速度连续的边界条件：
式中，Φi为入射平面波的声速度势，k0＝ω/c0为声波数；
根据式(12)，结合声模态函数的正交性，获得系数变量Ruw的表达式：
将式(13)带入式(11)中，获得反射声波的声速度势的表达式；
2.2)有源MPPA入射侧反射声波的声压表达
根据入射侧反射声波声速度势与反射声波声压的关系，获得入射侧(z＝-D)反射声波的声压表达式：
式中，
式(14)中入射侧总声压P＝Pi+Pr；
式(14)两边同乘以模态函数Ψuw(x,y)，并利用模态函数的正交性，获得如下的表达式：
将式(15)带入式(14)中，获得反射声波声压的表达式：
式中，
其中，Pi为入射波在有源MPPA表面处的值，
3)有源MPPA中微穿孔板和弹性平板振动位移的求解
根据模态叠加原理及微穿孔板与弹性平板振动位移满足的波动方程，推导出微穿孔板与弹性平板位移模态幅值满足的矩阵方程，联立两个矩阵，获得微穿孔板与弹性平板位移模态幅值的表达式；
3.1)微穿孔板振动响应建模
由附图1可知，微穿孔板受到入射侧的总声压P与空腔内微穿孔板表面声压PD的作用。根据微穿孔板振动位移满足的波动方程，同时利用模态叠加原理及模态函数的正交性，推导获得微穿孔板的位移模态幅值Amn满足的方程式：
式中，系数
ωmn与ξmn为第(m,n)阶模态的共振频率与模态阻尼比，
ρM与hM分别为微穿孔板的密度与厚度，
为微穿孔板的模态振型函数，
其中，系数
(m,n)为振动模态的模态序数；
将式(17)中微穿孔板和弹性平板的位移vM(x,y)与v2(x,y)进行模态展开，经化解可获得微穿孔板与弹性平板位移模态幅值Amn与Bmn满足的方程：
定义如下变量：
则M×N个未知的微穿孔板位移模态幅值与弹性平板的位移模态幅值满足如下矩阵方程：
C1A+D1B＝V (19)
式中，A＝[A11,L,Amn,L,AMN]T，B＝[B11,L,Bmn,L,BMN]T，V＝[V11,L,Vmn,L,VMN]T；
3.2)弹性平板振动响应建模
由附图1可知，弹性平板受到空腔内声压P0及次级控制力fs＝Fsδ(x-xs,y-ys)的作用。同理，根据弹性平板振动位移满足的波动方程，结合模态叠加原理及模态函数的正交性，获得微穿孔板与弹性平板的位移模态幅值Amn与Bmn满足如下的方程：
式中，
ω2,mn与ξ2,mn为弹性平板第(m,n)阶模态的共振频率与模态阻尼比，
ρ2与h2为弹性平板的密度与厚度；
Qs,mn为广义次级模态力，
如果进一步定义如下的变量：
则M×N个未知的微穿孔板与弹性平板位移模态幅值满足如下矩阵方程：
C2A+D2B＝V2+FsQs (21)
式中，V2＝[V2,11,L,V2,mn,L,V2,MN]T；
3.3)两种板的位移模态幅值求解
将微穿孔板与弹性平板位移模态幅值满足的矩阵方程(20)与(21)联立，获得如下的矩阵方程组：
进而获得微穿孔板与弹性平板的位移模态幅值为：
4)最优控制力的求解及控制后吸声系数的计算
将有源MPPA的吸声系数作为最优的控制目标，通过优化搜索算法获得使吸声系数最大的最优次级控制力幅值，进而可以获得控制后有源MPPA系统总的振动响应及控制后的吸声系数。
4.1)最优控制力的求解
对于有源控制系统，理论上最优的控制目标函数为有源MPPA的吸声系数，其可表示为：
式中，Πinc与Πabs分别为入射的总声能及有源MPPA吸收的声能，SMPP为微穿孔板的面积；由式(23)求得微穿孔板与弹性平板的位移模态幅值矢量A与B的表达式后，入射侧的总声压P(x,y)与表面振速v1(x,y)均可求得。由于A与B均为次级控制力幅值Fs的函数，因此，吸声系数α也为Fs的函数。针对上述较为复杂的隐式函数，可采用优化搜索算法如遗传算法求得使吸声系数α最大时的最优控制力幅值Fs,opt，就可求解出控制后有源MPPA的吸声系数。工程实现时，可构建恰当的误差传感策略获取与吸声系数非常相关的误差信号，通过自适应控制算法(Fx-LMS)进行实时跟踪控制。
4.2)控制后有源MPPA吸声系数的计算
获得最优控制力幅值Fs,opt后，带入式(23)便可获得最优控制下的位移模态幅值矢量A与B，即整个系统的振动响应即可获得；然后根据吸声系数的计算公式(24)，获得控制后有源MPPA的吸声系数。
为了验证该方法，本发明还通过以下仿真实验进一步说明：
一、确定有源MPPA的参数
如附图1所示，本实施例中微穿孔板与弹性平板的材质均为铝材，次级单点力源fs作用于弹性平板(0.05m，0.05m)位置处。有源MPPA的几何参数、材料属性参数及初级激励参数如表1～表3所示。仿真实验中微穿孔板和弹性平板的尺寸分别取①和②两种情形。
表1 有源MPPA的几何参数
表2 有源MPPA的材料属性参数
表3 初级激励参数
二、计算有源MPPA控制前后的吸声系数
按照以上参数计算有源MPPA控制前后的吸声系数。
三、仿真结果及分析
微穿孔板和弹性平板的尺寸分别为①和②两种情形时，计算获得有源MPPA控制前后的吸声系数，如附图2与3所示。从计算结果可以看出，控制后，MPPA在很宽的低频段内吸声系数几乎达到最优值1，充分说明该有源MPPA在低频段吸声性能的优越性。
综上，该发明不仅可以大幅提升传统MPPA低频段宽带的吸声性能，同时系统易于实现、且并没有MPPA空腔内要存在一维声场的限制条件，更加适合工程应用。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明公开的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本发明的保护范围之内。
Claims (7)
1.一种提升有源式微穿孔板吸声器低频吸声性能的方法，其特征在于：
所述有源式微穿孔板吸声器由微穿孔板及其后部空腔组合而成，所述后部空腔的后壁为点力控制的弹性平板；
具体提升方法包括以下步骤：
1)有源式微穿孔板吸声器空腔内声场的声压表达
根据有源式微穿孔板吸声器的几何模型参数、材料参数和初级激励与次级激励参数，对有源式微穿孔板吸声器空腔内的声场进行建模计算，推导出有源式微穿孔板吸声器空腔内微穿孔板和弹性平板表面的声压分布表达式；
2)有源式微穿孔板吸声器入射侧声场的声压表达
根据有源式微穿孔板吸声器入射侧界面上满足的速度连续边界条件，推导入射侧声场中反射声波的声压表达式，从而获得入射侧声场的总声压表达式；
3)有源式微穿孔板吸声器中微穿孔板和弹性平板振动位移的求解
根据模态叠加原理及微穿孔板与弹性平板振动位移满足的波动方程，推导出微穿孔板与弹性平板位移模态幅值满足的矩阵方程，联立两个矩阵，获得微穿孔板与弹性平板位移模态幅值的表达式；
4)最优控制力的求解及控制后吸声系数的计算
将有源式微穿孔板吸声器的吸声系数作为最优的控制目标，通过优化搜索算法获得使吸声系数最大的最优次级控制力幅值，进而可以获得控制后有源式微穿孔板吸声器系统总的振动响应及控制后的吸声系数。
2.根据权利要求1所述方法，其特征在于：
步骤1)中，所述有源式微穿孔板吸声器的几何模型参数包括微穿孔板的长度a、宽度b、厚度hM，弹性平板的厚度h2，空腔深度D，微穿孔板的孔隙直径d以及孔隙率σ；
所述有源式微穿孔板吸声器的材料参数包括微穿孔板的密度ρM、弹性平板的密度ρ2、微穿孔板的杨氏模量EM、弹性平板的杨氏模量E2、微穿孔板的泊松比σM、弹性平板的泊松比σ2、微穿孔板的模态阻尼比ξM,mn、弹性平板的模态阻尼比ξ2,mn、空气的密度ρ0、空气的声速c0；
所述有源式微穿孔板吸声器的初级激励与次级激励参数包括初级激励为垂直入射的平面波时，平面波的幅值为p0；次级控制点力的布置位置为(xs,ys)＝(0.05,0.05)时，控制力的幅值Fs。
3.根据权利要求2所述方法，其特征在于，步骤1)的具体步骤为：
1.1)有源式微穿孔板吸声器空腔内声场响应求解
在初级垂直入射平面波及次级控制力的激励下，有源式微穿孔板吸声器空腔内产生声场声压；根据空腔内声速度势Φ(x,y,z,t)满足的齐次波动方程，空腔内四个刚性侧壁处，存在空气质点振速为0的边界条件，利用模态叠加原理，结合格林第二公式以及模态函数的正交性，推导获得有源式微穿孔板吸声器空腔内声速度势的表达式为：
式中，ω为角频率，U与W为x与y方向上声模态序数的上限值，(u,w)为声模态序数，c0为空气的声速；
其中，系数变量Luw与Nuw取决于空腔内声速度势在微穿孔板与弹性平板两侧满足的如下速度连续的边界条件：
式中，D为空腔深度，v2(x,y,t)为弹性平板的质点振速，v1(x,y,t)为临近微穿孔板表面的质点振速，取微穿孔板振速和微孔中空气介质振速的空间平均值，具体表示为：
式中，vM(x,y)为微穿孔板的振速，P为入射侧的总声压，PD为腔内微穿孔板表面声压，σ为微穿孔板的孔隙率，z0为微穿孔板的孔隙阻抗；
根据边界条件(2)与(3)，结合空腔声模态函数的正交性获得声速度势表达式中的系数变量Luw与Nuw的表达式如下：
式中，声模态函数
系数
将式(5)与式(6)带入式(1)中，获得有源式微穿孔板吸声器空腔内任意点声速度势的表达式；
1.2)有源式微穿孔板吸声器腔内微穿孔板表面的声压表达
根据空腔内声速度势与声压的关系式获得腔内微穿孔板表面声压PD的表达式为：
式中，
由于式(7)的右端仍然含有自身PD，因此，将式(7)左右两端同乘以声模态函数Ψuw(x,y)，并利用模态函数的正交性，获得如下等式：
将式(8)带入式(7)中，获得腔内微穿孔板表面声压PD的表达式：
式中，
1.3)有源式微穿孔板吸声器腔内弹性平板表面的声压表达
根据空腔内声速度势与声压的关系式并利用式(8)，推导获得有源式微穿孔板吸声器腔内弹性平板表面声压P0的表达式：
式中，
4.根据权利要求3所述方法，其特征在于，步骤2)的具体步骤为：
2.1)有源式微穿孔板吸声器入射侧反射声波的声速度势求解
根据反射声波在入射侧管道内满足的齐次波动方程，结合四边刚性壁、入口处为开口的边界条件，获得反射声波声速度势Φr(x,y,z,ω)的表达式为：
式中，系数变量Ruw取决于入射侧表面满足的如下速度连续的边界条件：
式中，Φi为入射平面波的声速度势，k0＝ω/c0为声波数；
根据式(12)，结合声模态函数的正交性，获得系数变量Ruw的表达式：
将式(13)带入式(11)中，获得反射声波的声速度势的表达式；
2.2)有源式微穿孔板吸声器入射侧反射声波的声压表达
根据入射侧反射声波声速度势与反射声波声压的关系，获得入射侧反射声波的声压表达式：
式中，
式(14)中入射侧总声压P＝Pi+Pr；
式(14)两边同乘以模态函数Ψuw(x,y)，并利用模态函数的正交性，获得如下的表达式：
将式(15)带入式(14)中，获得反射声波声压的表达式：
式中，
其中，Pi为入射波在有源式微穿孔板吸声器表面处的值，
5.根据权利要求4所述方法，其特征在于，步骤3)的具体步骤为：
3.1)微穿孔板振动响应建模
根据微穿孔板振动位移满足的波动方程，同时利用模态叠加原理及模态函数的正交性，推导获得微穿孔板的位移模态幅值Amn满足的方程式：
式中，系数
ωmn与ξmn为第(m,n)阶模态的共振频率与模态阻尼比，
ρM与hM分别为微穿孔板的密度与厚度，
为微穿孔板的模态振型函数，
其中，系数
(m,n)为振动模态的模态序数；
将式(17)中微穿孔板和弹性平板的位移vM(x,y)与v2(x,y)进行模态展开，经化解可获得微穿孔板与弹性平板位移模态幅值Amn与Bmn满足的方程：
定义如下变量：
则M×N个未知的微穿孔板位移模态幅值与弹性平板的位移模态幅值满足如下矩阵方程：
C1A+D1B＝V (19)
式中，A＝[A11,…,Amn,…,AMN]T，B＝[B11,…,Bmn,…,BMN]T，V＝[V11,…,Vmn,…,VMN]T；
3.2)弹性平板振动响应建模
根据弹性平板振动位移满足的波动方程，同时利用模态叠加原理及模态函数的正交性，推导获得微穿孔板与弹性平板的位移模态幅值Amn与Bmn满足如下的方程：
式中，
ω2,mn与ξ2,mn为弹性平板第(m,n)阶模态的共振频率与模态阻尼比，
ρ2与h2为弹性平板的密度与厚度；
Qs,mn为广义次级模态力，
定义如下的变量：
则M×N个未知的微穿孔板与弹性平板位移模态幅值满足如下矩阵方程：
C2A+D2B＝V2+FsQs (21)
式中，
3.3)两种板的位移模态幅值求解
将微穿孔板与弹性平板位移模态幅值满足的矩阵方程(20)与(21)联立，获得如下的矩阵方程组：
进而获得微穿孔板与弹性平板的位移模态幅值为：
6.根据权利要求5所述方法，其特征在于，步骤4)的具体步骤为：
4.1)最优控制力的求解
将有源式微穿孔板吸声器的吸声系数作为最优的控制目标，表示为：
式中，Πinc与Πabs分别为入射的总声能及有源式微穿孔板吸声器吸收的声能，SMPP为微穿孔板的面积；
采用遗传算法求得使吸声系数α最大时的最优控制力幅值Fs,opt；
4.2)控制后有源式微穿孔板吸声器吸声系数的计算
获得最优控制力幅值Fs,opt后，带入式(23)获得最优控制下的位移模态幅值矢量A与B；然后根据吸声系数的计算公式(24)，获得控制后有源式微穿孔板吸声器的吸声系数。
7.根据权利要求1-6任一所述方法，其特征在于：
所述微穿孔板和弹性平板的材质均为铝。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110115539.1A CN112951188B (zh) | 2021-01-28 | 2021-01-28 | 一种有源式微穿孔板吸声器及提升其低频吸声性能的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110115539.1A CN112951188B (zh) | 2021-01-28 | 2021-01-28 | 一种有源式微穿孔板吸声器及提升其低频吸声性能的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112951188A CN112951188A (zh) | 2021-06-11 |
| CN112951188B true CN112951188B (zh) | 2024-05-03 |
Family
ID=76239396
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110115539.1A Active CN112951188B (zh) | 2021-01-28 | 2021-01-28 | 一种有源式微穿孔板吸声器及提升其低频吸声性能的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112951188B (zh) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1048248A (zh) * | 1990-04-21 | 1991-01-02 | 郑成仁 | 振动吸收式消音器 |
| DE19717266C1 (de) * | 1997-03-15 | 1998-04-09 | Fraunhofer Ges Forschung | Einrichtung für Schallabsorber mit mikroperforierten Bauteilen |
| JP2007011034A (ja) * | 2005-06-30 | 2007-01-18 | Kajima Corp | 微細穿孔板利用の吸音構造及び吸音材 |
| DE102009005577A1 (de) * | 2009-01-21 | 2010-07-29 | Thomas Friedrich | Thermoaktive Fertigteilplatte mit integrierter Funktionsplatte |
| CN202268160U (zh) * | 2011-10-12 | 2012-06-06 | 中国科学院合肥物质科学研究院 | 一种自适应微穿孔板吸声器 |
| CN102543061A (zh) * | 2012-01-20 | 2012-07-04 | 江苏大学 | 薄膜机械阻抗与微穿孔板声阻抗结合的宽频吸声结构 |
| CN102646414A (zh) * | 2012-05-14 | 2012-08-22 | 南京大学 | 基于微穿孔和腔内共振系统的组合吸声结构 |
| CN103559877A (zh) * | 2013-07-17 | 2014-02-05 | 南京大学 | 一种基于分流扬声器和微穿孔板的复合吸声结构 |
| CN104715749A (zh) * | 2015-03-17 | 2015-06-17 | 中国科学院合肥物质科学研究院 | 基于自适应微穿孔板吸声器的声阻抗调节装置及调节方法 |
| CN105427853A (zh) * | 2015-10-30 | 2016-03-23 | 东南大学 | 宽频带微穿孔板吸声体及其性能预测方法、结构设计方法 |
| CN105609264A (zh) * | 2016-01-19 | 2016-05-25 | 江西理工大学 | 一种采用微穿孔板的非晶合金油浸式变压器 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2256722A1 (de) * | 2009-05-29 | 2010-12-01 | Akusik & Innovation GmbH | Schalldämpfendes bzw. -absorbierendes Material |
-
2021
- 2021-01-28 CN CN202110115539.1A patent/CN112951188B/zh active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1048248A (zh) * | 1990-04-21 | 1991-01-02 | 郑成仁 | 振动吸收式消音器 |
| DE19717266C1 (de) * | 1997-03-15 | 1998-04-09 | Fraunhofer Ges Forschung | Einrichtung für Schallabsorber mit mikroperforierten Bauteilen |
| JP2007011034A (ja) * | 2005-06-30 | 2007-01-18 | Kajima Corp | 微細穿孔板利用の吸音構造及び吸音材 |
| DE102009005577A1 (de) * | 2009-01-21 | 2010-07-29 | Thomas Friedrich | Thermoaktive Fertigteilplatte mit integrierter Funktionsplatte |
| CN202268160U (zh) * | 2011-10-12 | 2012-06-06 | 中国科学院合肥物质科学研究院 | 一种自适应微穿孔板吸声器 |
| CN102543061A (zh) * | 2012-01-20 | 2012-07-04 | 江苏大学 | 薄膜机械阻抗与微穿孔板声阻抗结合的宽频吸声结构 |
| CN102646414A (zh) * | 2012-05-14 | 2012-08-22 | 南京大学 | 基于微穿孔和腔内共振系统的组合吸声结构 |
| CN103559877A (zh) * | 2013-07-17 | 2014-02-05 | 南京大学 | 一种基于分流扬声器和微穿孔板的复合吸声结构 |
| CN104715749A (zh) * | 2015-03-17 | 2015-06-17 | 中国科学院合肥物质科学研究院 | 基于自适应微穿孔板吸声器的声阻抗调节装置及调节方法 |
| CN105427853A (zh) * | 2015-10-30 | 2016-03-23 | 东南大学 | 宽频带微穿孔板吸声体及其性能预测方法、结构设计方法 |
| CN105609264A (zh) * | 2016-01-19 | 2016-05-25 | 江西理工大学 | 一种采用微穿孔板的非晶合金油浸式变压器 |
Non-Patent Citations (1)
| Title |
|---|
| 有源消声的理论进展与发展趋向;陈克安;《噪声与振动控制》;19920831(第第4期期);第28-33页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112951188A (zh) | 2021-06-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gao et al. | Design, fabrication and sound absorption test of composite porous metamaterial with embedding I-plates into porous polyurethane sponge | |
| JP6970880B2 (ja) | 音響メタマテリアル騒音制御法およびダクトシステムにおける装置 | |
| US20240046908A1 (en) | Thin-layer low-frequency underwater sound insulation metamaterial | |
| CN112687252B (zh) | 调控夹层板结构低频隔声性能的装置及其参数确定方法 | |
| Sakagami et al. | Absorption characteristics of a space absorber using a microperforated panel and a permeable membrane | |
| CN112951188B (zh) | 一种有源式微穿孔板吸声器及提升其低频吸声性能的方法 | |
| CN116011122A (zh) | 一种计算周期异形管道超材料声传输的计算方法 | |
| Que et al. | Tunable low frequency band gaps and sound transmission loss of a lever-type metamaterial plate | |
| CN110135062B (zh) | 圆形多孔热防护材料的优化设计方法 | |
| Huang | A theory of reactive control of low-frequency duct noise | |
| Liu et al. | Analysis of an ultra-low frequency and ultra-broadband phononic crystals silencer with small size | |
| Sugahara | A 3D-printed sound-absorbing material based on multiple resonator-like unit cells for low and middle frequencies | |
| CN116343736A (zh) | 一种提升大尺寸微穿孔板吸声器斜入射低频吸声性能的有源控制方法 | |
| CN112528399B (zh) | 一种基于能量集中汇聚的水下平台壁板振动宽频控制方法 | |
| CN210597696U (zh) | 一种吸声复合结构单元及具有其的吸声体阵列 | |
| Goto et al. | Sound attenuation of periodic micro-perforated chamber mufflers using the spectral transfer matrix method | |
| CN116825067A (zh) | 一种实现水下低频宽带全吸收的薄式声学超表面 | |
| CN114242027B (zh) | 一种基于多孔材料的复合吸声超构表面 | |
| CN217108003U (zh) | 应用于水下航行器低频减振降噪的局域共振型平板结构 | |
| Zhang | Vibration and acoustic radiation control of a panel with piezoelectric oscillators. | |
| Bobrovnitskiĭ | Theory of the new high-efficiency absorbing and nonscattering coatings | |
| Aygun | Vibroacoustic analysis of composite thin fiberglass plate | |
| Liu et al. | A Hybrid Numerical Analytical Method for Calculating Hydrodynamic Self-Noise of Ribbed Acoustic Window | |
| CN116453494A (zh) | 一种实现水下低频宽带吸收的封闭式耐压水声吸收器 | |
| Wang et al. | Research on Ultrabroadband Acoustic Absorbers Based on Slow-wave Metamaterials |
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 |