删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

薄膜与Helmholtz腔耦合结构低频带隙

本站小编 Free考研考试/2021-12-29

摘要:设计了一种含薄膜壁的Helmholtz型声子晶体, 该结构利用了空气和薄膜的耦合振动, 一方面将刚性壁转变为柔性壁, 降低了一阶振动时的等效刚度, 使第一带隙下限分别低于同参数下的普通Helmholtz型声子晶体和薄膜, 另一方面基于局域共振原理, 由于薄膜的出现和腔口空气通道长度的增加, 使得结构在低频范围内存在多个振动模态, 从而将原有一个带隙扩展为多个带隙. 将该结构带隙上下限分别等效为环形系统和串联系统, 用传递矩阵法和有限元法两种方法计算了其低频带隙范围, 两种方法结果吻合良好. 通过调整参数对带隙调控规律进行了进一步分析, 结果显示, 在低频范围内, 既可以通过改变与腔口空气通道或薄膜相关的参数, 在保证其中某些带隙变化不大的情况下, 单独调整其他带隙; 也可以通过调整内外腔体积, 对所有带隙进行调控.
关键词: Helmholtz共振腔/
低频带隙/
噪声控制/
声子晶体

English Abstract


--> --> -->
自声子晶体提出以来[1], 其一直是声学领域国内外****研究重点之一[27], 另一方面, Helmholtz共振腔这种声学器件在声学超材料中有着广泛的应用[812]. 近年来, 利用Helmholtz共振原理构建的声子晶体及超材料已逐渐由简单叠加结构[13,14]逐渐发展为长开口[15]、多腔[16]、多开口[17,18]、嵌套[19]、多层复合[20]等结构.
由于空气密度的限制, 仅利用空气的振动往往无法进一步提升声子晶体的性能, 与此同时, 薄膜作为一种轻质材料, 其在低频方面也具有较好的隔声性能[2123]. 能否利用两种结构的耦合, 构建出低频隔声性能更好的声子晶体, 就成为值得探讨的问题. 实际上, 近年来已有关于含薄膜的Helmholtz腔仿真研究[24,25], 腔体与薄膜耦合[2628]以及薄膜与穿孔板-腔耦合[29]等结构的研究出现. 但薄膜与Helmholtz腔耦合结构的理论计算与声子晶体带隙研究仍较少.
本文在之前Helmholtz腔与固/固型声子晶体的耦合研究[30]基础上, 设计了一种含薄膜壁的Helmholtz型声子晶体, 对其带隙机理进行了详细分析, 用传递矩阵法(transfer matrix method, TMM)和有限单元法(finite element method, FEM)计算了其低频带隙上下限. 该结构第一带隙下限分别低于同参数下的普通Helmholtz型声子晶体和薄膜, 且质量小于同尺寸传统Helmholtz型声子晶体, 进一步提高了Helmholtz腔在小尺寸、轻结构下控制大波长的能力, 提高了其在工程上的应用价值.
带薄膜壁的Helmholtz结构横截面如图1所示, 其晶格常数为a, 腔体框架边长为l, 腔壁厚度为b, 在腔右侧由悬臂梁形成“W”型开口, 其空气通道总长度为l1 = n × (lb) + b, 宽度为s, 其中n为悬臂梁的个数. 将腔体左侧壁更换为厚度为br的硅橡胶薄膜, 并在其上黏附有长hs、厚bs的铝制质量块, 薄膜受到y方向张力T的作用. 因框架材料一般为金属, 其声阻抗一般在空气的105倍以上, 薄膜材料的103倍以上, 对带隙的影响较小, 故将其设定为固定约束状态.
图 1 带薄膜壁的Helmholtz结构横截面
Figure1. Cross section of Helmholtz resonator structure with a membrane wall.

a = 53 mm, l = 50 mm, b = 1 mm, n = 2 (l1 = 99 mm), s = 1 mm, br = 1 mm, hs = 5 mm, bs = 1 mm, T = 1 × 106 N/m2, 先将其按照第一布里渊区进行扫描, 再将其沿纵向对3个元胞结构进行串联, 分别计算得出其在1700 Hz以下的结构能带图和隔声量曲线如图2所示. 从图2可以看出, 其在1700 Hz以下存在3个完全带隙(灰色区域), 分别为88.40—119.06 Hz, 302.09—533.03 Hz和772.31—891.44 Hz (各带隙起止点已在图中标出), 与此同时出现了多个平直带; 其对应隔声量曲线分别在各带隙下限处出现了40 dB以上的隔声峰. 若将薄膜也设定为固定约束状态, 则该结构变为普通二维Helmholtz结构, 用同样的方法计算得出的结构能带图和隔声量曲线如图3所示, 其在1700 Hz以下范围只存在1个完全带隙(116.60—318.34 Hz), 最大隔声峰为36 dB. 同时, 通过FEM计算得出其在相同条件下的薄膜基频为240.59 Hz.
图 2 带薄膜壁的Helmholtz结构 (a) 带隙图; (b) 隔声曲线
Figure2. Band diagram (a) and transmission spectrum (b) of the Helmholtz resonator structure with a membrane wall.

图 3 普通Helmholtz结构的(a)带隙图和(b)隔声曲线
Figure3. Band diagram (a) and transmission spectrum (b) of the ordinary Helmholtz resonator structure.

通过以上分析可以发现, 将Helmholtz型声子晶体的一个刚性壁换为带分布质量的张紧膜后, 其低频隔声性能得到了提升. 具体表现为: 第一带隙下限得到进一步降低, 且同时低于同条件下的普通Helmholtz结构和薄膜结构; 出现了新的隔声峰, 且高度高于原有结构; 虽然第一带隙的宽度减小, 但在低频范围内出现了新的带隙, 使得总带隙宽度得到提升.
为研究薄膜与Helmholtz腔的耦合作用, 取该结构前3个带隙起止点处声压场和薄膜振型进行分析, 如图4所示, 其中纯灰色部分为薄膜振型(位移经放大), 右侧图例为声压场数值, 单位为Pa.
图 4 (a) 模态A (88.40 Hz)、(b) 模态B (119.06 Hz)、(c) 模态C (302.09 Hz)、(d) 模态D (533.03 Hz)、(e) 模态E (772.31 Hz)、(f) 模态F (891.44 Hz) 的薄膜振型和声场压力图
Figure4. Vibration mode of the membrane and sound pressure distribution diagrams of point A (88.40 Hz) (a), B (119.06 Hz) (b), C (302.09 Hz) (c), D (533.03 Hz) (d), E (772.31 Hz) (e), and F (891.44 Hz) (f).

图4可以看出, 在模态A, C, E处, 结构声压场变化规律完全相同, 均为内腔声压最大, 并通过腔口空气通道逐渐过渡到外腔. 外腔左右两部分声压呈反对称分布, 其中薄膜侧为正, 腔口侧为负, 且这种差异随着带隙阶数的增大而增强, 但外腔声压和均为零. 此时声波被完全局域在内腔中, 振动与外腔无关, 与其对应于带隙下限相匹配. 而在模态B, D, F处, 结构声压场分布与前述相反, 内腔声压最小, 且为负值, 通过腔口空气通道过渡至外腔, 外腔声压最大. 此时振动与内腔外腔都有关, 声波可以在腔外传播, 对应于带隙上限.
由于膜的振动是各阶主振型叠加的结果, 通过振型图仅能推断某阶主振型占主要地位, 在后续分析中, 将占主要地位的某阶主振型称为其某阶振动. 从振型图可以看出, 随着频率的升高, 薄膜振动逐渐由低阶转向高阶, 但在带隙上下限处均没有发现反对称振型(该种振动模态下薄膜上下位移呈反对称分布, 平均位移为零)的参与.
对于出现多个平直带的原因, 与之前研究得出的结论相同[30], 是由薄膜的反对称振型造成的, 这里不再进行讨论.
另外, 对于模态A, 可以看出膜与腔口空气做同向振动, 这样实际上减小了内腔空气弹簧刚度, 导致第一带隙下限下降; 与此类似, 模态B中内腔空气弹簧刚度增大, 外腔减小, 但由于两者体积变化比例不同, 其总体刚度是减小的, 导致其第一带隙上限也会下降.
经过以上分析可看出, 对于该结构在1700 Hz以下产生的多个带隙, 其不同带隙上限或下限处声压场分布规律均是相同的, 只是薄膜振动模态不同, 但各带隙上限和下限的声压场分布规律不同.
在此对上下限分别构建等效系统, 如图5所示, 其中X1表示薄膜平均位移; X2表示腔口通道内空气质心位移; N1N2分别为内腔、外腔对薄膜的总压力, 采用TMM与连续体振动相结合的方法进行计算.
图 5 (a) 带隙下限系统示意图; (b) 带隙上限系统示意图
Figure5. (a) System corresponding to starting frequency of band gaps; (b) system corresponding to cut-off frequency of band gaps.

将腔口通道内空气视为均质弹性杆, 其传递矩阵[31]
${{ U}_1} = \left[ {\begin{array}{*{20}{c}} {\cos \dfrac{\omega }{c}{l_1}}&{\dfrac{{ - \sin \dfrac{\omega }{c}{l_1}}}{{\omega c{\rho _{\rm{a}}}s}}} \\ {\omega c{\rho _{\rm{a}}}s\sin \dfrac{\omega }{c}{l_1}}&{\cos \dfrac{\omega }{c}{l_1}} \end{array}} \right],$
其中m1 = al1, ρa为空气密度, c为在空气中的声速, ω为角频率. 设内外腔内空气在振动过程中压强均匀, 则内腔空气传递矩阵为
${{ U}_2} = \left[ {\begin{array}{*{20}{c}} {\dfrac{{{l_{\rm{r}}}}}{s}}&{ - \dfrac{{{V_2}}}{{{\rho _{\rm{a}}}{c^2}s{l_{\rm{r}}}}}} \\ 0&{\dfrac{s}{{{l_{\rm{r}}}}}} \end{array}} \right],$
外腔空气的传递矩阵为
${{ U}_4} = \left[ {\begin{array}{*{20}{c}} {\dfrac{s}{{{l_{\rm{r}}}}}}&{ - \dfrac{{{V_4}}}{{{\rho _{\rm{a}}}{c^2}s{l_{\rm{r}}}}}} \\ 0&{\dfrac{{{l_{\rm{r}}}}}{s}} \end{array}} \right],$
其中V2, V4分别表示内腔和外腔体积; l r = l – 2b, 为薄膜的长度.
对于薄膜的纵向振动, 采用Rayleigh-Ritz法[32]求解, 同时考虑张力和弹性模量的影响, 其强迫振动方程为
${E_{\rm{r}}}{J_{\rm{r}}}\frac{{{\partial ^4}y}}{{\partial {x^4}}} + T\frac{{{\partial ^2}y}}{{\partial {x^2}}} + {\rho _{\rm{r}}}{b_{\rm{r}}}\frac{{{\partial ^2}y}}{{\partial {t^2}}} = p(x,t),$
其中Erρr分别表示硅橡胶弹性模量和密度; Jr为截面对中性轴的惯性矩; p为薄膜收到的外力.
取基础函数为φi = 1 – cos(2πnx/l r), 这种取法计算简便, 但舍弃了反共振振型, 故不能计算出平直带的振动频率. 考虑薄膜上质量块的分布作用, 此时其等效刚度矩阵和等效质量矩阵中各元素为
$\begin{split} {K_{ij}} \; &= {E_{\rm{r}}}{J_{\rm{r}}}\int_0^{{l_{\rm{r}}}} {{{\varphi }''_i}} {{\varphi}''_j}{\rm{d}}x + T{b_{\rm{r}}}\int_0^{{l_r}} {{{\varphi}'_i}{{\varphi}'_j}{\rm{d}}x} \\& = \left\{ {\begin{aligned}& {\frac{{8{E_{\rm{r}}}{J_{\rm{r}}}{i^4}{{\text{π}}^4}}}{{{l_{\rm{r}}}^3}} + \frac{{2T{b_{\rm{r}}}{i^2}{{\text{π }}^2}}}{{{l_{\rm{r}}}}},}\quad\quad\quad{i = j,} \\ & {0,}\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad {i \ne j,} \end{aligned}} \right. \end{split}$
${M_{ij}} = {\rho _{\rm{r}}}{b_{\rm{r}}}\int_0^{{l_{\rm{r}}}} {{\varphi _i}{\varphi _j}{\rm{d}}x} + {\rho _{\rm{s}}}{b_{\rm{s}}}\int_{\frac{{{l_{\rm{r}}} - {h_{\rm{s}}}}}{2}}^{\frac{{{l_{\rm{r}}} + {h_{\rm{s}}}}}{2}} {{\varphi _i}{\varphi _j}{\rm{d}}x},$
由此可解出特征值矩阵Λ和其对应的特征向量矩阵A.
振动时, 薄膜在空气的作用下, 相当于受到周期性均布激振力的作用, 设均布力为p sinωt, 则正则广义力为
${Q_i}\left( t \right) = p\sin \omega t\int_0^{{l_{\rm{r}}}} {{\varphi _i}{\rm{d}}x} = {l_{\rm{r}}}p\sin \omega t.$
此时方程(4)可写为
$\ddot \xi + {{\varLambda}} \xi = {{{A}}^{\rm{T}}}{{Q}}\left( t \right),$
其解为
${\xi _i} = \frac{1}{{{\varLambda _{ii}} - {\omega ^2}}}{B_i}p\sin \omega t,$
其中BiATQ(t)的第i个元素除以psinωt后的结果. 令η = Aξ, 则薄膜在空气作用下的稳态响应为
$y(x,t) = p\sin \omega t\sum\limits_{i = 1}^n {{\varphi _i}} {\eta _i},$
则其在振动过程中对空气造成的最大体积改变量为
$\begin{split} pV\left( \omega \right) \; &= p\int_0^{{l_{\rm{r}}}} {\sum\limits_{i = 1}^n {{\varphi _i}} {\eta _i}{\rm{d}}x} \\ &= p\sum\limits_{i = 1}^n {\left[ {\frac{1}{{{\varOmega _i}^2 - {\omega ^2}}}\sum\limits_{j = 1}^n {{{\left( {{\alpha _{ij}}{l_{\rm{r}}}} \right)}^2}} } \right]},\end{split}$
其中αijA中各元素, Ωi为矩阵Λ中对角线元素.
通过各传递矩阵及(11)式, 可分别对带隙上下限对应的系统进行求解.
对带隙下限, 设传递顺序为内腔-腔口, 则有
${{{U}}_{{\rm{down}}}} = {{ U}_1}{{ U}_2},$
$\left[ {\begin{array}{*{20}{c}} {{X_2}} \\ 0 \end{array}} \right] = {{{U}}_{{\rm{down}}}}\left[ {\begin{array}{*{20}{c}} {{X_1}} \\ {{N_1}} \end{array}} \right].$
并且, 根据薄膜平均位移在TMM和Rayleigh-Ritz法下计算结果应相同, 可得
${X_1}{l_{\rm{r}}} = \frac{{{N_1}}}{{{l_{\rm{r}}}}}V\left( \omega \right).$
联立(13)和(14)式可得
$\begin{split} & \left| {\begin{array}{*{20}{c}} {U_{{\rm{down}}}^{{\rm{1,1}}}}&{U_{{\rm{down}}}^{1,2}}&{ - 1} \\ {U_{{\rm{down}}}^{2,1}}&{U_{{\rm{down}}}^{{\rm{2,2}}}}&0 \\ { - {l_{\rm{r}}}}&{ - V\left( \omega \right)/{l_{\rm{r}}}}&0 \end{array}} \right|= \\ & \frac{{{V_2}\omega \sin \dfrac{{{l_1}\omega }}{c} + V(\omega ){c^2}{\rho _{\rm{a}}}\omega \sin \dfrac{{{l_1}\omega }}{c} - cs\cos \dfrac{{{l_1}\omega }}{c}}}{c}\\ & = 0,\end{split} $
其中$ U_{\rm down}^{i,j}$表示Udown中各元素.
同样, 对带隙上限, 设传递顺序为内腔-腔口-外腔, 则有
${{{U}}_{{\rm{up}}}} = {{ U}_4}{{ U}_1}{{ U}_2},$
$\left[ {\begin{array}{*{20}{c}} {{X_1}} \\ {{N_2}} \end{array}} \right] = {{{U}}_{{\rm{up}}}}\left[ {\begin{array}{*{20}{c}} {{X_1}} \\ {{N_1}} \end{array}} \right],$
${X_1}{l_{\rm{r}}} = \frac{{{N_2} - {N_1}}}{{{l_{\rm{r}}}}}V\left( \omega \right).$
联立(17)和(18)式可得
$\begin{split} & \left| {\begin{array}{*{20}{c}} {U_{{\rm{up}}}^{{\rm{1,1}}} - 1}&{U_{{\rm{up}}}^{{\rm{1,2}}}}&0 \\ {U_{{\rm{up}}}^{{\rm{2,1}}}}&{U_{{\rm{up}}}^{{\rm{2,2}}}}&{ - 1} \\ { - {l_{\rm{r}}}}&{ - V\left( \omega \right)/{l_{\rm{r}}}}&{V\left( \omega \right)/{l_{\rm{r}}}} \end{array}} \right|=0, \\ & {V_2}{V_4}{\omega ^2} + [\left( {{V_2} + {V_4}} \right)V(\omega ){c^2}{\rho _{\rm{a}}}{\omega ^2} - {c^2}{s^2}]\sin \frac{{{l_1}\omega }}{c} \\ & - \left[ {\left( {{V_2} + {V_4}} \right)cs\omega + 2V(\omega ){c^3}{\rho _{\rm{a}}}s\omega } \right]\cos \frac{{{l_1}\omega }}{c} \\ & + V(\omega ){c^3}{\rho _{\rm{a}}}s\omega \left( {{{\sin }^2}\frac{{{l_1}\omega }}{c} + {{\cos }^2}\frac{{{l_1}\omega }}{c}} \right) \\ &+ V(\omega ){c^3}{\rho _{\rm{a}}}s\omega = 0,\\[-17pt]\end{split} $
其中$ U_{\rm up}^{i,j}$ 表示Uup中各元素.
根据(15)与(19)式可分别计算出带隙下限与带隙上限的频率. 从这两式可以看出, 在低频范围内, 式子左端取得零值的主要因素有式中的V(ω)和其他含ω项, 其中V(ω)与薄膜的振动模态有关, 其他含ω项均来自于弹性杆传递矩阵. 这说明随着频率的增大, 每当薄膜或腔口空气的振动模态发生改变时, 都将出现一个新的带隙, 亦即产生了一种新的局域共振模态. 另外, 由于两者的耦合性及内外腔空气的作用, 带隙的上下限将不会出现在原固有频率处, 而是发生一定的偏移.
为验证理论计算方法的适用性, 并进一步研究带隙形成规律, 用FEM和TMM计算其带隙上下限频率随参数改变的变化情况, 并取其前三个带隙进行分析. 在此, 除研究的参数外, 其余参数与上文中相同. 从(15)和(19)式可以看出, 影响该结构带隙的主要参数有V2, V4, l1, s以及与薄膜相关的函数V(ω), 在此仅有针对性的选取个别参数进行分析, 其中误差项是以FEM所得结果作为真实值计算得出.
从上文理论计算可发现, 薄膜受到张力T、重物质量及分布、薄膜长度l r等因素会并且只会影响薄膜相关的函数V(ω), 在此首先选取了重物长度l s作为变量进行分析, 结果如表1所列. 同时, 作为比较, 采用FEM计算了同条件下附加金属片薄膜的纵向振动固有频率(不含反共振频率), 结果如表2所列. 从表1表2可以看出, 增大l s, 薄膜一阶固有频率下降, 二阶固有频率增大, 而结构第二、第三带隙变化趋势与之完全相同, 变化幅度也很接近, 而结构第一带隙向低频方向移动, 但变化幅度较小. 该现象说明此结构在1700 Hz以下新出现的第二、三带隙分别是由于薄膜出现了前两阶振动模态引起的. 而第一带隙由于仍然对应于腔口空气的振动, 通过增大l s的方式增加等效质量是一种间接的调控方式, 对于该带隙的优化效果并不理想.
l s/10–3 m468101214
1阶固有频率252.4237.9226.8218.2211.6206.6
2阶固有频率751.3782.8794.0814.8843.4879.1


表2薄膜附加金属片长度l s对薄膜固有频率的影响
Table2.Effect of the parameter l s on natural frequency of membrane.

ls/10–3 m第一带隙下限第一带隙上限第二带隙下限第二带隙上限第三带隙下限第三带隙上限
FEM
TMM
误差/%FEM
TMM
误差/%FEM
TMM
误差/%FEM
TMM
误差/%FEM
TMM
误差/%FEM
TMM
误差/%
489.23.7121.44.4314.22.8557.70.6790.72.2900.31.6
92.5126.7323.0561.2808.3914.7
688.93.6120.14.2297.71.8526.8–0.8789.10.4905.20.8
92.2125.1303.1522.8792.4912.2
888.73.5118.93.9285.10.9508.3–2.1795.6–1.1916.1–0.5
91.9123.5287.6497.5786.8912.0
1088.63.4117.93.5275.6–0.1497.9–3.5815.9–3.7932.5–2.3
91.6122.0275.4480.5786.0910.8
1288.53.1117.23.0268.4–1.0492.6–4.8843.8–6.9953.9–5.2
91.3120.7265.6469.1785.3904.4
1488.52.8116.72.3262.9–2.0490.5–5.9878.5–11966.3–8.2
91.0119.4257.6461.5779.9887.3


表1薄膜附加金属片长度l s对低频带隙的影响
Table1.Effect of the parameter l s on low-frequency band gaps.

另外, 当l s较小时, 两种计算结果接近, 但当l s > 10 × 10–3 m后, 误差开始显著增大, 这是由于在用Rayleigh-Ritz法对薄膜进行处理时, 仅通过(6)式对分布质量进行了处理, 而忽略了附加金属片对薄膜等效刚度的影响. 随着l s增大, 这种影响逐渐增大, 导致了误差不断增大.
表3显示的是薄膜张力对带隙的影响, 可以看出, 随着薄膜张力的增大, 其各带隙上下限均有增大的趋势, 但第二、三带隙的增长幅度大于第一带隙, 这与上文中所得出第二、三带隙对应于薄膜的振动模态产生和改变相一致. 对于第一带隙, 由于其对应的是腔口空气的振动模态, 可在分析时忽略薄膜质量的影响, 此时随着张力增大, 结构趋向于刚性壁, 其带隙上下限逐渐与无薄膜结构接近. 当张力增大到108 N/m2后, 其第一带隙上下限已基本与无薄膜结构一致, 且1700 Hz以下已无其他完整带隙.
T/106 N·m–1第一带隙下限第一带隙上限第二带隙下限第二带隙上限第三带隙下限第三带隙上限
FEM
TMM
误差/%FEM
TMM
误差/%FEM
TMM
误差/%FEM
TMM
误差/%FEM
TMM
误差/%FEM
TMM
误差/%
0.574.33.489.32.8259.20.7440.3–1.2574.10.1770.10.5
76.891.8261.0434.9574.4774.3
1.596.33.4143.14.6345.53.2589.60.2952.12.21035.11.7
99.6149.7356.7591.0973.01053.1
2.5103.32.9174.64.6415.04.2648.00.71217.42.81274.12.3
106.4182.7432.6652.41251.01303.2
3.5106.82.7196.44.4474.94.8691.41.11434.13.01478.62.5
109.6205.0497.6698.81477.61516.1
4.5108.82.5212.54.1528.25.1729.61.41621.63.21642.53.3
111.5221.3555.2740.11673.41696.5
10113.02.1257.83.0757.35.9907.92.71645.95.81741.83.1
115.3265.6801.7932.81740.61796.2
100116.21.8311.01.21654.15.21737.33.12270.89.02375.36.1
118.3314.91740.41791.62475.72520.3


表3薄膜张力T对低频带隙的影响
Table3.Effect of the parameter T on low-frequency band gaps.

腔口空气通道长度l1对低频带隙的影响如表4所列, 可以看出, 随着l1的增大, 第一带隙上下限均向低频方向移动, 而第二带隙下限变化不大, 这与上文提出的对应关系相符合.
l1/mm第一带隙下限第一带隙上限第二带隙下限第二带隙上限第三带隙下限第三带隙上限
FEM
TMM
误差/%FEM
TMM
误差/%FEM
TMM
误差/%FEM
TMM
误差/%FEM
TMM
误差/%FEM
TMM
误差/%
9988.44.4119.15.7302.13.4533.01.3772.33.4891.42.4
92.3125.9312.5540.0798.6912.9
14874.53.4101.55.1301.23.5513.31.1772.63.4873.72.4
77.1106.7311.7519.1798.9894.4
19766.02.989.94.7301.93.5500.01.1772.43.3836.01.9
67.994.1312.5505.3798.2851.8
24660.12.681.44.5303.23.5488.51.0697.61.2734.10.9
61.785.1313.7493.3705.9740.6
29555.82.575.04.3304.73.4475.70.9587.80.7637.50.7
57.278.2315.2479.9591.9642.3
34452.52.469.84.2306.33.4458.40.7507.50.6558.93.1
53.772.7316.6461.5510.8576.4


表4腔口空气通道长度l1对低频带隙的影响
Table4.Effect of the parameter l1 on low-frequency band gaps.

但第二带隙上限也不断下降, 特别是当l1大于246 mm后, 第三带隙上下限急剧下降. 从带隙图分析发现, 腔口空气二阶振动对应的带隙(第四带隙)随着l1的增大不断向低频方向移动, 压缩了第三带隙及第二带隙上限. 直至l1 = 295 mm后, 腔口空气二阶振动对应的带隙下降到薄膜二阶振动对应带隙以下, 成为第三带隙, 如图6所示, 腔口空气表现为中间压强最大, 两端最小. 该现象说明随着l1的增大, 腔口空气在1700 Hz以下范围内的振动模态增多, 固有频率下降. 实际上, 当l1 = 344 mm时, 该结构在1700 Hz以下已有6个带隙, 其分别对应于腔口空气一阶振动、薄膜一阶振动、腔口空气二阶振动、薄膜二阶振动、腔口空气三阶振动和腔口空气四阶振动.
图 6 l1 = 295 mm时(a) 第三带隙下限和(b) 第三带隙上限的声场压力图
Figure6. Sound pressure distribution diagrams at starting frequency (a) and cutoff frequency (b) of the 3th band gap when l1 = 295 mm.

由(15)式可以看出, 外腔体积V4不影响带隙下限, 而由(19)式可以看出, 对带隙上限, 内腔体积V2与外腔体积V4的作用完全相同, 故只对内腔体积V2进行分析, 如表5所列. 从表5可以看出, 随着V2的增大, 各带隙上下限均向低频方向移动, 这是由于腔体积增加会减小其等效刚度. 另外, 在此减小V2的方式是在内腔中增加刚性填充物, 这种方法会使得内腔形状不规则, 腔内声压不均匀, 导致误差上升. 但即便刚性填充物占内腔比例达到66% (此时内腔体积为7.07 × 10–4 m3), 最大误差仍较小, 说明本文采用的理论计算方法也适用于其他较为复杂结构.
V2/10–4 m3第一带隙下限第一带隙上限第二带隙下限第二带隙上限第三带隙下限第三带隙上限
FEM
TMM
误差/%FEM
TMM
误差/%FEM
TMM
误差/%FEM
TMM
误差/%FEM
TMM
误差/%FEM
TMM
误差/%
7.07108.56.3120.55.8406.91.0566.21.0809.91.9953.11.2
115.3127.6411.0571.8825.3964.2
10.57103.34.8121.14.8361.41.6558.3–0.1796.51.9925.41.4
108.3126.8367.2557.8811.4938.7
14.0798.04.2121.04.4335.01.9550.0–0.1790.91.8913.11.4
102.2126.4341.3549.4804.9925.8
17.5793.33.9120.94.3317.72.1544.4–0.1788.31.6906.51.3
96.9126.1324.4543.9801.1918.1
19.3291.13.8120.84.3311.12.2542.3–0.1787.41.6904.11.2
94.5126.0317.9541.8799.7915.3


表5内腔体积V2对低频带隙的影响
Table5.Effect of the parameter V2 on low-frequency band gaps.

从整体上看, 该种带薄膜壁的Helmholtz结构可变参数很多, 且各参数对不同带隙的影响程度不尽相同. 因此, 在低频范围内, 既可以通过改变与腔口空气通道或薄膜相关的参数, 在保证其中某些带隙变化不大的情况下, 单独调整其他带隙; 也可以通过调整内外腔体积, 对所有带隙进行调控.
本文设计了一种含薄膜壁的Helmholtz型声子晶体, 建立了系统等效模型, 通过研究发现: 1)该结构第一带隙下限分别低于同参数下的普通Helmholtz型声子晶体和薄膜, 且在1700 Hz以下范围内将原有一个带隙扩展为多个带隙; 2)利用TMM与Rayleigh-Ritz法相结合的方式可以较为精确地计算出其带隙上下限; 3)该结构的各带隙分别对应于腔口空气和薄膜的各阶振型, 因此可以通过改变与腔口空气通道或薄膜相关的参数, 在保证其中某些带隙变化不大的情况下, 单独调整其他带隙; 与此同时, 内外腔体积对所有带隙均有影响, 故可通过调整其体积对所有带隙进行调控. 这些结论对构建Helmholtz腔与薄膜耦合声学超材料具有指导意义, 有利于推动低频隔声技术的发展.
相关话题/结构 计算 质量 系统 材料

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 薄膜底面Helmholtz腔声学超材料的隔声性能
    摘要:针对低频噪声的隔离问题,设计了一种薄膜底面Helmholtz腔声学超材料,该超材料由薄膜底面Helmholtz腔附加质量单元构成.使用有限元法,计算了超材料在20—1200Hz频段内的传输损失曲线与各阶共振频率,并在实验中验证了数值计算的真实性.研究结果表明,超材料在20—1200Hz频段内隔 ...
    本站小编 Free考研考试 2021-12-29
  • 塑料类高分子聚合物材料水中目标声学参数反演
    摘要:塑料类高分子聚合物材料作为3D打印领域的一类主要材料,在水下声学模型和结构成型中的应用越来越广泛,这类新材料的声学参数直接影响着3D打印的水下声学模型和结构的声学性能.基于水中目标声散射的Rayleigh简正级数解及塑料类高分子聚合物材料实心球中亚音速Rayleigh波低频共振机理分析,获取了 ...
    本站小编 Free考研考试 2021-12-29
  • 利用海洋环境噪声空间特性估计浅海海底分层结构及地声参数
    摘要:海洋环境噪声场中包含了海洋中的诸多信息,海底地声参数是影响海洋环境噪声场空间分布的主要因素之一.对于不同的海底分层结构,海底反射损失会根据沉积层厚度和各层声速呈现出不同的临界角和干涉条纹结构.本文利用Harrison能流理论,从理想反射系数出发,分别考虑了声速、密度、衰减系数、沉积层厚度等几种 ...
    本站小编 Free考研考试 2021-12-29
  • 全固态磁制冷系统物理模型的研究进展
    摘要:磁制冷是一种节能环保的制冷技术,具有广阔的应用前景.目前,基于主动磁回热循环的磁制冷系统被广泛研究并诞生了多个原型制冷机.然而,这些系统主要采用流体换热,导致系统存在工作频率低、回热损失大、子部件设计复杂等问题,使得制冷机成本升高和效率降低.针对上述问题和难点,引入固态传热增强机制和全固态磁制 ...
    本站小编 Free考研考试 2021-12-29
  • Li离子电池负极材料石墨炔在B, N掺杂调控下的储Li性能优化
    摘要:一种理想的Li离子电池负极材料需要具有较高的储Li容量和较低的体积膨胀比.本文应用密度泛函理论研究了二维多孔石墨炔在B,N原子掺杂调控后作为Li离子电池负极材料时的储Li性能.计算结果表明,B在石墨炔结构中的掺杂可以增强Li与石墨炔之间的吸附作用,储Li容量可以增加到2061.62mAh/g, ...
    本站小编 Free考研考试 2021-12-29
  • V-4Cr-4Ti/Ti复合材料界面的辐照损伤特性研究
    摘要:采用透射电子显微镜和纳米压痕技术研究了高能Fe10+离子注入后V-4Cr-4Ti/Ti复合材料界面辐照硬化、辐照缺陷以及微观结构的变化.根据样品成分和微观结构的特征,将该复合材料分成V-4Cr-4Ti基体、界面I、界面II和Ti基体四个区域.纳米压痕结果表明,辐照后样品均发生了辐照硬化,但界面 ...
    本站小编 Free考研考试 2021-12-29
  • 缺陷离子调控对BiFeO<sub>3</sub>-BaTiO<sub>3</sub>基钙钛矿材料的铁电光伏特性影响
    摘要:BiFeO3-BaTiO3铁电材料具有优异的压电和铁电性能,近年来受到广泛的关注.该材料既保持了BiFeO3体系高的自发极化强度Ps的优点,也克服了BiFeO3体系难以合成纯钙钛矿相等缺点,被认为是非常有前景的铁电、压电以及光伏材料.本文采取传统固相法制备了Bi(Fe0.96Mg0.02–xT ...
    本站小编 Free考研考试 2021-12-29
  • 结构光照明技术在二维激光诱导荧光成像去杂散光中的应用
    摘要:背景杂散光信号的干扰制约了激光片光成像技术的发展,本文将结构光照明技术应用到激光片光成像测量中来消除杂散光的干扰.介绍了基于结构光照明技术的工作原理和实验测量系统,基于Matlab软件理论分析了相位移动结构光照明技术具有完全消除杂散光的作用,并针对稳定罗丹明B溶液进行了二维激光诱导荧光成像实验 ...
    本站小编 Free考研考试 2021-12-29
  • 大质量转动沃尔夫-拉叶星的形成及内部核合成研究
    摘要:沃尔夫-拉叶(Wolf-Rayetstars,WR)星是一类非常特殊的恒星,具有强烈的星风损失,造成氢包层丢失,仅具有裸露的氦核.WR星被认为是Ib/Ic型超新星的前身星,研究WR星的形成及内部核合成具有重要意义.根据转动恒星的角动量转移和元素扩散方程,研究了影响WR星结构与演化的各种物理因素 ...
    本站小编 Free考研考试 2021-12-29
  • 一种基于压缩感知和多维混沌系统的多过程图像加密方案
    摘要:随着计算机科学的快速发展,信息的存储和传播常常在各类计算机硬件以及多种网络之间进行,传统的信息加密方案已逐渐不再适用.因此,基于计算机的信息加密算法近年来逐步成为研究热点.通过结合小波包变换、压缩感知、混沌系统等理论,一种基于压缩感知和多维混沌系统的多过程图像加密方案被提出.该加密方案实现了针 ...
    本站小编 Free考研考试 2021-12-29