-
基于弯曲波在声黑洞(acoustic black hole, ABH)结构中振幅不断增大的特性, 提出了一种新型径向夹心式径-弯复合换能器, 该换能器由径向夹心式圆环换能器与外围的环形ABH结构组成. ABH结构的存在实现了换能器径向振动与弯曲振动之间的转换, 提高了换能器的声辐射性能. 利用几何声学的方法建立了ABH结构弯曲振动的解析模型, 给出了其弯曲振动的本征频率, 并结合有限元方法讨论了换能器机电转换性能随尺寸变化的关系. 通过有限元方法给出了该换能器在空气中的辐射声压场、辐射声强以及辐射指向性, 仿真结果表明, ABH结构的存在能够改善换能器弯曲振动的机电转换性能, 提高换能器的声辐射性能, 使换能器呈现出一定的辐射指向性. 最后通过实验对换能器样机的电阻抗特性以及振动模态进行测量, 实验结果与仿真相符合.A novel radial sandwich radial bending composite transducer is proposed based on the increasing amplitude of bending waves in the structure of an acoustic black hole (ABH). The transducer consists of a radial sandwich circular ring transducer and an outer acoustic black hole structure. The existence of the ABH structure realizes the conversion between radial vibration and bending vibration of the transducer, and improves the acoustic radiation performance of the transducer. Based on the one-dimensional acoustic black hole bending vibration theory, the bending vibration equation of the circular ABH structure is given. An analytical model of the ABH structure bending vibration is established by using the method of geometrical-acoustics, and the eigen frequency of its bending vibration is given. The relationship between the electromechanical conversion performance and the size of the transducer is discussed by using the finite element method, and the results show that when the eigen frequency of the bending vibration of the ABH structure is close to the resonant frequency of the transducer, the electromechanical coupling coefficient of the transducer reaches a maximum value. The radiation sound pressure field, radiation sound intensity, and radiation directionality of the transducer in the air are also analyzed by the finite element method. The results indicate that the presence of the ABH structure can improve the electromechanical conversion performance of the bending vibration of the transducer, enhance the acoustic radiation performance of the transducer, make the transducer exhibit radiation directionality and the lower order of the bending vibration mode of the ABH structure, the stronger radiation directionality of the transducer towards the left side and the right side. Finally, in order to verify the accuracy of the simulation and theory, we made a prototype transducer and carried out the experimental test of transducer electrical impedance and vibration mode. The measured results are in good agreement with the simulation results. The displacement distribution of the radiation surface of the transducer at resonance frequency is measured, which verifies the bending vibration mode of the ABH structure. This type of transducer is expected to be used as acoustic projector in the air or water.
-
Keywords:
- acoustic black hole /
- transducer /
- radial vibration /
- bending vibration
1. 引 言
径向夹心式圆环换能器由金属圆环和用压电材料制成的圆环复合而成. 其性能稳定, 具有声辐射面积大、接收全指向性等优点, 被广泛应用于水声技术、地质勘探、超声化学反应等各种超声技术中[1–3].
声黑洞是一种奇特声学现象, 由Mironov[4]在1988年首先提出, 可简述为弯曲波在楔形结构中传播时, 其波速将随板厚减小而减小, 而振幅逐渐增大. 理想情况下, 当板厚按照二次以上幂指数减小到零时, 即可实现弯曲波的零反射, 因此这种结构被称为声黑洞. 随后Krylov等[5–7]将其应用于一维梁中, 并提出利用这种结构来操纵弯曲波, 实现减振降噪等应用.
当弯曲波在声黑洞结构中传播时, 其振幅随厚度减小而增大的特性使得弯曲波的能量汇聚在结构尖端处并以声能的形式向外辐射出去. 声黑洞结构的这种特性被Remillieux等[8,9]提出可以用来改善换能器与空气之间的阻抗匹配, 提高换能器在空气中的辐射性能.
压电复合圆环型换能器除了径向振动模态还具有弯曲振动模态, 相对于径向振动模态, 弯曲振动模态的振动位移较大, 因此更易于空气等介质相匹配. 然而关于弯曲振动模态的研究比较少, 直到近十年来才有学者做过相关方面的研究. Aronov[10]在2013年用能量法对复合圆环型换能器的弯曲振动进行了研究, 给出了其弯曲振动模态的等效机电电路图以及有效机电耦合系数的表达式, Aronov[11]同时也研究了开槽后压电圆环的弯曲振动模态. 吴德林等[12]在Aronov[10,11]的基础上对压电复合圆环型换能器弯曲振动的四极子模态进行了研究. 相对于径向振动, 复合圆环型换能器弯曲振动模态的机电转换性能较差, 这种缺点在一定程度上限制了其弯曲振动模态的应用.
本文借助振动模态转换的方法来改善圆环型换能器弯曲振动模态机电转换性能. 模态转换是一种借助多个模态进行换能器设计的方法, 能够利用两种不同振动模态的优点从而满足换能器的设计性能要求, 被广泛应用于换能器的设计中[13–15]. 本文的主要设计思路是将一种环形的ABH结构应用于换能器结构设计中, 相比较于一维ABH梁的平直结构, 环形ABH能减小结构尺寸, 同时也可以与夹心式圆环换能器复合在一起, 将换能器的径向振动转换为弯曲振动, 这样既提高了换能器弯曲振动模态的机电转换性能又利用弯曲振动大位移的特性增强了换能器的声辐射能力.
2. ABH结构的弯曲振动分析
图1(a)为所研究的径向夹心式径-弯复合换能器俯视图, 该换能器由一个径向夹心式圆环换能器与外围环组成. 径向夹心式圆环换能器则由两个金属环与其中间的压电环夹心而成, 压电环的极化方向沿径向.
图1(b)为截取外围环1/4部分的示意图, 其由厚度均匀部分与变厚度ABH部分组成, 厚度均匀部分与径向夹心式圆环换能器耦合在一起, 将径向夹心式圆环换能器的径向振动转换为ABH部分的弯曲振动. 外围环的厚度随弧度变化, 厚度表达式可以写为
$$ h\left(\theta \right)=\left\{\begin{aligned} &\frac{{h}_{1}-{h}_{0}}{{\left(m{\mathrm{\pi }}/2\right)}^{p}}{\theta }^{p}+{h}_{0},&& 0\leqslant \theta < m{\mathrm{\pi }}/2,\\ &{h}_{1},&& m{\mathrm{\pi }}/2\leqslant \theta \leqslant {\mathrm{\pi }}/2,\end{aligned}\right. $$ (1) 其中p为环厚度随弧度变化的幂次, $ {h}_{0} $为ABH部分$ \theta =0 $处的最小厚度, 而$ {h}_{1} $则为ABH部分的最大厚度.
图1(c) 为该换能器的振动模态示意图, 出于观察方便, 将图中换能器的振动位移进行放大, 换能器的实际振动位移小于图中示意.
为了更方便地研究该换能器的径-弯复合振动模态, 将该换能器的振动模态分成两部分来研究: ABH部分的弯曲振动模态和径向夹心式圆环换能器的径向振动模态. 由于径向夹心式圆环换能器的径向振动模态已有较多的研究[16–19], 因此我们主要讨论ABH部分的弯曲振动模态.
关于ABH部分的研究, 可以借助一维ABH梁理论, 一维ABH结构中弯曲振动控制方程可以写成如下形式[5]:
$$ {\nabla }^{2}\left(D\left(\theta \right){\nabla }^{2}w\right)-{\omega }^{2}\rho h\left(\theta \right)w=0, $$ (2) 其中$ {\nabla }^{2} $为拉普拉斯算子、w表示挠度、ρ为密度、$ \omega $为角频率、$ \text{D}\left(\theta \right)={Eh}^{3}\left(\theta \right)/[12(1-{\nu }^{2})] $为抗弯刚度, ν和E分别为材料的泊松比和杨氏模量. $ \omega = 2{\mathrm{\pi }}f $ 为角频率.
在二维极坐标下ABH部分的振动控制方程(2)可以写为
$$ \frac{1}{{r\left(\theta \right)}^{2}}\frac{{\partial }^{2}}{\partial {\theta }^{2}}\left(\frac{D\left(\theta \right)}{{r\left(\theta \right)}^{2}}\frac{{\partial }^{2}w}{\partial {\theta }^{2}}\right)-{\omega }^{2}\rho h\left(\theta \right)w=0. $$ (3) 选取ABH部分中面处为参考位置, 则r$ \left(\theta \right) $=$ {r}_{4}+ {h}_{1}-h\left(\theta \right)/2 $. 通过几何声学近似的方法, 可以将该方程的挠度解表示为如下复数形式:
$$ \begin{array}{c}w\left(\theta \right)=A\left(\theta \right){{\mathrm{e}}}^{{\mathrm{i}}\varPhi \left(\theta \right)},\end{array} $$ (4) 其中A($ \theta $)和$ \varPhi \left(\theta \right)={k}_{{\mathrm{p}}}s\left(\theta \right) $ 分别为弯曲波随弧度变化的幅值和累积相位, 其中$ {k}_{{\mathrm{p}}}=\omega /{c}_{{\mathrm{p}}} $为厚度均匀环中的波数, $ {c}_{{\mathrm{p}}}=2{{c}_{{\mathrm{t}}}(1-{c}_{{\mathrm{t}}}^{2}/{c}_{{\mathrm{l}}}^{2})}^{1/2} $ 为弯曲波的相速度, $ {c}_{{\mathrm{l}}} $和$ {c}_{{\mathrm{t}}} $分别为材料中纵波和横波的波速. 将(4)式代入方程(3)之中, 略去$ {k}_{{\mathrm{p}}} $ 的1, 2, 3次项, 可以得到弯曲波的程函方程:
$$ { |\nabla s (\theta)|}^{4}={k}^{4}/{k}_{{\mathrm{p}}}^{4}={R}^{4}{\rho {{\omega }}}^{2}{\varepsilon (\theta )}^{2}h (\theta)/(D(\theta){k}_{{\mathrm{p}}}^{4}), $$ (5) $ \nabla $为梯度算子, $ R = {r}_{4} + {h}_{1} $, $ \varepsilon (\theta) = (R - h(\theta))/R, \; k\left(\theta \right)={R\omega }^{1/2}{\varepsilon\left(\theta \right)}^{1/2}{\left(\rho h\left(\theta \right)/D\left(\theta \right)\right)}^{1/4} $, 方程(5)的一个解可以表示为[7]
$$ s (\theta) = \frac{1}{{k}_{{\mathrm{p}}}}\int {R{{\omega }}}^{1/2}{{\varepsilon (\theta )}^{1/2} (\rho h (\theta)/D (\theta))}^{1/4}{\mathrm{d}}\theta . $$ (6) 因此弯曲波在ABH部分中的累计相位$ \varPhi \left(\theta \right)= {k}_{{\mathrm{p}}}s\left(\theta \right) $可以表示为
$$ \varPhi =\int _{0}^{\theta }k\left(\theta \right){\mathrm{d}}\theta . $$ (7) 弯曲波在ABH部分的累积相位$ \varPhi \left(\theta \right) $可写为由$ \theta =0 $处到$ \theta = $m $ {\mathrm{\pi }}/2 $的积分表达式:
$$ \varPhi ={\int }_{0}^{m{\mathrm{\pi }}/2}k\left(\theta \right){\mathrm{d}}\theta . $$ (8) (8)式中积分的下限为ABH部分的最薄处, 也就是$ \theta =0 $处, 积分上限对应ABH部分的最厚处, 即$ \theta =m{\mathrm{\pi }}/2 $. 当变厚度部分中弯曲波的累计相位为[20]
$$ \begin{array}{c}\varPhi =2n\pi ,\end{array} $$ (9) 此时ABH部分可以看作为对应n倍弯曲波波长的振子, 其中n的取值为$ N/2 $, N取正整数. 对方程(9)求解, 即得到ABH部分n阶弯曲振动模态所对应的本征频率. 需要注意的是得出上面的结果需要满足几何声学近似条件[7]:
$$ \left|\frac{1}{{k}^{2}}\frac{{\mathrm{d}}k}{{\mathrm{d}}\theta }\right|\ll 1. $$ (10) 对于本文所研究的ABH结构, 代入铝的材料参数以及$ k\left(\theta \right) $, (10)式可以化简为
$$ \left|\frac{1}{{k}^{2}}\frac{{\mathrm{d}}k}{{\mathrm{d}}\theta }\right| < 19.73\frac{p\sqrt{{{\Delta }}h}}{({r}_{4}-{h}_{1})\sqrt{{({\theta }_{1})}^{p}\omega }}{\theta }^{\tfrac{p}{2}-1}, $$ (11) 其中$ h(\theta ) = {{\Delta }}h{\theta }^{p}/{\theta }_{1}^{p} + {h}_{0} $, $ {{\Delta }}h = {h}_{1} - {h}_{0} $, $ {\theta }_{1} $=$ m{\mathrm{\pi }}/2 $. 因为$ 0\leqslant \theta \leqslant {\theta }_{1} $; 将(11)式进一步化简得到:
$$ 19.73\frac{p\sqrt{{{\Delta }}h}}{({r}_{4}-{h}_{1})\sqrt{{({\theta }_{1})}^{p}\omega }}{\theta }^{\tfrac{p}{2}-1} \leqslant 19.73\frac{p\sqrt{{{\Delta }}h}}{({r}_{4}-{h}_{1}){\theta }_{1}\sqrt{\omega }}. $$ (12) 由(12)式可知, 对几何声学近似条件影响最大的是弯曲振动的频率, 本文研究的结构尺寸以及频率范围(>18 kHz)满足该近似条件.
3. 换能器的有效机电耦合系数分析
当改变ABH部分的m参数也即变厚度部分的长度时, ABH部分弯曲振动某一n阶本征频率与该换能器的谐振频率将在某一m值处相等, 径向夹心式圆环换能器与外围环的耦合最强点将出现在该m值附近. 此时径向夹心式圆环换能器的径向振动与ABH部分n阶弯曲振动的耦合达到最强, 换能器在该点的辐射性能也最佳.
换能器的机电耦合系数是描述换能器在能量转换过程中能量相互耦合程度的一个物理量, 在这里我们用换能器的有效机电耦合系数来描述ABH部分与夹心式圆环换能器的耦合程度. 当二者之间的耦合达到最强时, 换能器的有效机电耦合系数最大.
换能器的有效机电耦合系数$ {k}_{{\mathrm{e}}{\mathrm{f}}{\mathrm{f}}} $定义为[21] $ {k}_{{\mathrm{e}}{\mathrm{f}}{\mathrm{f}}}^{2}= ({f}_{{\mathrm{n}}}^{2}-{f}_{{\mathrm{m}}}^{2})/{f}_{{\mathrm{n}}}^{2} $ , 其中$ {f}_{{\mathrm{m}}} $为最大导纳频率, 即回路中电流最大时所对应的换能器谐振频率, $ {f}_{{\mathrm{n}}} $为最小导纳频率, 即回路中电流最小时所对应的频率.
本文设计了3种不同尺寸的径向夹心式径-弯复合换能器(换能器的具体尺寸见表1), 并对换能器有效机电耦合系数随m的变化关系进行讨论. 径向夹心式径-弯复合换能器的有效机电耦合系数 $ {k}_{{\mathrm{e}}{\mathrm{f}}{\mathrm{f}}} $ 通过COMSOL Multiphysics有限元仿真给出.
换能器编号 $ {r}_{1}/{\mathrm{m}}{\mathrm{m}} $ $ {r}_{2}/{\mathrm{m}}{\mathrm{m}} $ $ {r}_{3}/{\mathrm{m}}{\mathrm{m}} $ $ {r}_{4}/{\mathrm{m}}{\mathrm{m}} $ p $ {h}_{1}/{\mathrm{m}}{\mathrm{m}} $ $ {h}_{0}/{\mathrm{m}}{\mathrm{m}} $ a 30 35 40 45 3 5 1.0 b 30 35 40 45 2 5 1.5 c 30 35 40 45 2 5 2.0 d 30 35 40 45 — — — 表1给出的是几种不同换能器的尺寸, 其中a, b, c为径向夹心式径-弯复合换能器, d为夹心式圆环换能器.
图2给出了a, b, c三种换能器在不同m下的有效机电耦合系数 $ {k}_{{\mathrm{e}}{\mathrm{f}}{\mathrm{f}}} $, 径向夹心式径-弯复合换能器的谐振频率以及ABH部分所对应的n阶弯曲振动本征频率. 换能器的谐振频率以及有效机电耦合系数$ {k}_{{\mathrm{eff}}} $由有限元仿真给出, ABH部分n阶弯曲振动本征频率由求解方程(9)给出. 从图2可以看到, 随着m的变化换能器的$ {k}_{{\mathrm{e}}{\mathrm{f}}{\mathrm{f}}} $将会出现一个最大值. 换能器a的$ {k}_{{\mathrm{e}}{\mathrm{f}}{\mathrm{f}}} $ 最大值约为0.289, 换能器b的$ {k}_{{\mathrm{e}}{\mathrm{f}}{\mathrm{f}}} $ 最大值约为0.289, 换能器c的$ {k}_{{\mathrm{e}}{\mathrm{f}}{\mathrm{f}}} $ 最大值约为0.291.
Aronov[10]给出了圆环换能器弯曲振动模态的有效机电耦合系数, 其得到的弯曲模态的有效机电耦合系数$ {k}_{{\mathrm{e}}{\mathrm{f}}{\mathrm{f}}} \approx 0.6{k}_{31} $. 本文使用的压电材料为PZT-4, 其$ {k}_{31} $为0.33, 因此其弯曲模态所对应的$ {k}_{{\mathrm{e}}{\mathrm{f}}{\mathrm{f}}} $ $ \approx $ 0.198. 在吴德林等[12]的论文中, 其复合圆环换能器弯曲振动四极子模态的最大 $ {k}_{{\mathrm{e}}{\mathrm{f}}{\mathrm{f}}} $不超过0.6$ {k}_{31} $, 而从图2中可以看到, 我们所设计的3种径向夹心式径-弯复合换能器的最大$ {k}_{{\mathrm{e}}{\mathrm{f}}{\mathrm{f}}} $都大于复合圆环弯曲振动模态的有效机电耦合系数$ 0.6{k}_{31} $. 这充分说明借助径向振动所激发的弯曲振动的机电转换性能优于复合圆环换能器的弯曲振动模态.
设换能器谐振频率与ABH部分n阶本征频率随m的变化曲线的交点为$ {m}_{1} $, 径向夹心式径-弯复合换能器最大机电耦合系数所对应的m为$ {m}_{2} $. 从图2中可以看出, 径向夹心式径-弯复合换能器最大机电耦合系数所对应的$ {m}_{2} $值处于换能器谐振频率与ABH部分弯曲振动谐振频率交点值$ {m}_{1} $附近. 设$ {m}_{1} $与$ {m}_{2} $ 之间的绝对差值$ \varDelta =\left|{m}_{2}-{m}_{1}\right| $, 从图2中可以看出, 换能器a的$ {\varDelta }_{{\mathrm{a}}}=0.004 $、换能器b的$ {\varDelta }_{{\mathrm{b}}}=0.02 $ 、换能器c的$ {\varDelta }_{{\mathrm{c}}}=0.02 $.
这种寻找曲线交点的方法一方面有助于分析出与夹心式圆环换能器径向振动所耦合的ABH部分弯曲振动模态的阶数, 图2分别标注出了ABH部分所对应弯曲振动模态的阶次, 图2(a)中n = 2, 而图2(b), (c)中n = 1.5. 除此之外我们也可以通过计算换能器的谐振频率给出其频率随尺寸参数m变化的曲线, 并求解方程(9)给出ABH部分本征频率随尺寸参数m变化的曲线, 然后根据两条曲线的交点确定耦合最强时m值的大致范围, 从而可以大大减少仿真的计算量. 因为换能器谐振频率及ABH部分n阶弯曲振动本征频率与m的关系接近于线性, 所以可以借助少量的计算进行曲线拟合来找到曲线交点的位置, 并进一步通过有限元法确定耦合最强点所对应的m值.
4. 径向夹心式径-弯复合换能器的声场分析
本节给出了a, b, c三种不同型号径向夹心式径-弯复合换能器其最大导纳频率$ {f}_{{\mathrm{m}}} $所对应的空气中辐射声压分布、声强分布以及指向性图案. 在仿真的过程中采用二维平面建模. 仿真中使用的压电材料为PZT-4, 压电材料的极化方向沿径向, 压电环内外两壁面施加的有效电压差为1 V. 仿真中所用铝材料的杨氏模量和泊松比分别为$ E= 70\;{\mathrm{G}}{\mathrm{P}}{\mathrm{a}},\;\nu =0.33 $, 密度$ \rho =2700\;{\mathrm{k}}{\mathrm{g}}/{{\mathrm{m}}}^{3} $.
4.1 声压分布
图3给出的是几种换能器在空气中的辐射声压分布以及对应的谐振频率$ {f}_{{\mathrm{m}}} $, 其中a, b, c为表1给出的3种径向夹心式径-弯复合换能器, 3种换能器的m尺寸参数为图2中机电耦合系数最大时所对应的值, 图3(d)则为夹心式圆环换能器的辐射声压图.
由上述分析可知, 对ABH部分的尺寸进行合理优化设计, 使得ABH结构与夹心式圆环换能器的耦合达到最强, 可以提高换能器的声辐射能力.
4.2 声强分布
图4中a, b和 c为3种具有不同尺寸径-弯复合换能器的辐射声强分布图, 图4(d)为夹心式圆环换能器的辐射声强分布图, 图中声强的单位为W/m2. 从图4可以看到, 具有ABH结构的径向夹心式径-弯复合换能器其最大辐射声强相较于圆环型换能器有大幅度的提高. 其中以换能器b的辐射能力最强. 换能器a, c的最大辐射声强也有不同程度的提高. 相较于图3, 图4可以反映出声能量的辐射方向, 由于ABH部分处于换能器的左右两侧, 因此3种换能器的声辐射主要集中在换能器的两侧.
4.3 指向性分析
图5给出的是距离换能器中心1 m处4种换能器的指向性图. 图5中a, b和 c为图3中3种不同尺寸径-弯复合换能器的指向性图案, d为夹心式圆环型换能器的指向性图案. 相较于夹心式圆环型换能器的全指向性, 径-弯复合换能器的指向性图案出现较多小旁瓣, 且左右两侧的旁瓣较为尖锐. ABH结构的存在使得换能器向左右两侧的辐射能力得到提升. 相比图5中a型号换能器, b, c型号换能器指向性图案中左右两侧瓣更为尖锐. b, c型号换能器与ABH部分n = 1.5阶弯曲振动相耦合, 而a型号换能器则与n = 2阶弯曲振动相耦合. 由仿真结果和上述分析可知, 当换能器与ABH部分较低阶弯曲振动相耦合时, 其左右两侧的辐射指向性更强.
5. 实验验证
为了验证基于声黑洞理论设计径-弯复合换能器的实际可行性, 制作了如图6所示的换能器样机, 其中金属内外环材料为铝, 压电陶瓷为径向极化的PZT-4圆环, 换能器的尺寸参数见表2. 利用阻抗分析和激光测振实验对换能器样机的电参数特性以及振动模态进行实验分析.
尺寸参数 r1/mm r2/mm r3/mm r4/mm p h1/mm h0/mm m 22 26 30 34 2 4 1 0.65 5.1 换能器阻抗分析实验
图7为有限元仿真所得到的换能器阻抗随频率的变化曲线, 从图中可知换能器的仿真谐振频率频率$ {f}_{{\mathrm{F}}{\mathrm{E}}{\mathrm{M}}} $约为23.682 kHz. 为了验证仿真结果, 采用WK6500B精密阻抗分析仪对换能器的阻抗进行了测量. 从图8(b)可以看出, 阻抗分析仪测得换能器的谐振频率$ {f}_{{\mathrm{E}}{\mathrm{X}}} $约为23.323 kHz, 实验测量结果与仿真的误差$ \varDelta =\left|{f}_{{\mathrm{F}}{\mathrm{E}}{\mathrm{M}}}-{f}_{{\mathrm{E}}{\mathrm{X}}}\right|/{f}_{{\mathrm{E}}{\mathrm{X}}}= $1.54%. 误差的主要来源有: 1) 在理论分析和数值模拟中, 使用了标准的材料参数, 而换能器实物的材料参数与仿真参数有一定差距; 2) 在仿真中忽略了换能器的机械/介电损耗, 但是在实验测量中不能忽略; 3) 由于实物加工精度以及换能器装配的原因, 使得仿真模型与实物有一定的区别. 同时由于换能器样机的谐振频率低于仿真值, 使得换能器另一个弯曲模态(对应图7中第一个谷点)也被激发出, 因此在23 kHz附近出现两个谐振点.
5.2 换能器模态分析实验
图9为利用有限元仿真测得的换能器在谐振频率23.682 kH时的振型. 为进一步验证仿真结果, 使用图10(a)的Polytec PSV-400全场扫描式激光振动测量系统对换能器侧边ABH环部分的振型进行测量并与仿真结果进行对比.
图10(b) 为利用激光测振仪对换能器侧边ABH环位移振幅频响的测量结果, 实验中换能器两端施加的是电压幅值为1 V的扫频电信号. 测量结果中, 换能器位移共振峰所对应的频率为23.305 kHz, 与阻抗分析仪所测的谐振频率值接近. 从图10(b)的振型测量结果可以看出, 换能器侧边ABH环的振型与图9有限元仿真结果的振型相一致, 测量结果中换能器侧边ABH环厚度最小处的位移振幅最大, 这也符合弯曲波在声黑洞结构中振幅随壁厚减小而不断增大的特性. 综上可知, 借助ABH结构将换能器径向振动转换为ABH部分弯曲振动的方法是切实可行的.
6. 结 论
本文利用环形ABH结构设计了一种新型径向夹心式径-弯复合换能器, 实现了换能器径向振动与弯曲振动之间的转换, 并通过实验对换能器的径-弯复合模态进行验证, 利用阻抗分析仪测得的换能器谐振频率与仿真的误差为1.54%, 利用激光测振仪测得的换能器的振型与仿真结果一致. 从理论计算与仿真结果可以得到以下几点结论.
1) 当换能器总体谐振频率与ABH部分某一阶弯曲振动本征频率接近时, 换能器径-弯复合模态的有效机电耦合系数最大.
2) 换能器径-弯复合模态的最大有效机电耦合系数值是大于复合压电圆环型换能器的弯曲振动模态的.
3) 借助ABH部分将夹心式圆环换能器径向振动转换为弯曲振动的径-弯复合模态可以提高换能器的声辐射能力.
4) 不同于径向夹心式圆环换能器的无指向性, 具有ABH结构的径向夹心式径-弯复合换能器表现出一定的指向性, ABH结构提高了换能器向左右两侧的辐射能力, 并且与换能器相耦合的ABH部分的弯曲振动阶数越低, 其左右两侧的辐射指向性越强.
本研究不仅证明在换能器设计中加入声黑洞结构的可行性, 也为径-弯复合换能器的设计提供了参考模型.
[1] 胡银丰 2017 声学与电子工程 39 25
Hu Y F 2017 Acoust Electr Eng. 39 25
[2] 曹雪砷, 王易敏, 王申 2015 测井技术 39 160Google Scholar
Cao X S, Wang Y M, Wang S 2015 Well Logging Technol. 39 160Google Scholar
[3] 刘世清, 麻磊磊 2020 陕西师范大学学报(自然科学版) 48 60
Liu S Q, Ma L L 2020 J. Shaanxi Normal Univ. (Nat. Sci. Ed. ). 48 60
[4] Mironov M A 1988 Sov. Phys. Acoustics. 34 318
[5] Krylov V V 1989 Sov. Phys. Acoustics. 35 176
[6] Krylov V V 1990 Mosc. U. Phys. B+. 45 65
[7] Krylov V V, Tilman F 2004 J. Sound Vib. 274 605Google Scholar
[8] Remillieux M C, Anderson B E, Le Bas P Y, Ulrich T J 2014 Ultrasonics 54 1409Google Scholar
[9] Anderson B E, Remillieux M C, Le Bas P Y, Ulrich T J, Pieczonka L 2015 Ultrasonics 63 141Google Scholar
[10] Aronov B S 2013 J. Acoust. Soc. Am. 134 1021Google Scholar
[11] Aronov B S 2013 J. Acoust. Soc. Am. 133 3875Google Scholar
[12] 吴德林, 戴郁郁, 陈浩, 张海澜, 王秀明 2020 声学学报 45 103Google Scholar
Wu D L, Dai Y Y, Chen H, Zhang H L, Wang X M 2020 Acta Acustica. 45 103Google Scholar
[13] 林书玉, 鲜小军 2014 陕西师范大学学报(自然科学版) 42 31
Lin S Y, Xian X J 2014 J. Shaanxi Normal Univ. (Nat. Sci. Ed. ). 42 31
[14] 许龙, 林书玉 2012 声学学报 37 408Google Scholar
Xu L, Lin S Y 2012 Acta Acustica. 37 408Google Scholar
[15] 李凤鸣, 刘世清, 许龙, 张海岛, 曾小梅, 陈赵江 2023 中国科学: 物理学 力学 天文学 53 198
Li F M, Liu S Q, Xu L, Zhang H D, Zeng X M, Chen Z J 2023 Sci. Sin. Phys. Mech. Astron. 53 198
[16] 董宜雷, 陈诚, 林书玉 2023 物理学报 72 054304Google Scholar
Dong Y L, Chen C, Lin S Y 2023 Acta Phys. Sin. 72 054304Google Scholar
[17] 高炳山, 林书玉 2010 应用声学 29 217Google Scholar
Gao B S, Lin S Y 2010 J. Appl. Acoust. 29 217Google Scholar
[18] 胡静, 林书玉, 王成会, 莫润阳 2015 西北大学学报(自然科学版) 45 709
Hu J, Lin S Y, Wang C H, Mo R Y 2015 J. Northwest Univ. (Nat. Sci. Ed.) 45 709
[19] 王帅军, 林书玉 2011 陕西师范大学学报(自然科学版) 39 23
Wang S, Lin S 2011 J. Shaanxi Normal Univ. (Nat. Sci. Ed.) 39 23
[20] Krylov V V 2020 J. Sound Vib. 468 115100Google Scholar
[21] 林书玉 2004 超声换能器的原理及设计 (北京: 科学出版社) 第20页
Lin S Y 2004 The Theory and Design of Ultrasonic Transducers (Beijing: Science Press) p20
-
图 1 径向夹心式径-弯复合换能器示意图 (a) 换能器总体示意图; (b) 外围环1/4部分的示意图; (c) 换能器振动模态示意图
Fig. 1. The schematic diagram of a composite radial bending transducer: (a) The overall schematic diagram of the transducer; (b) a schematic diagram of the 1/4 part of the outer ring; (c) a schematic diagram of the vibration mode of the transducer.
图 2 不同m下换能器的有效机电耦合系数、谐振频率以及ABH部分的弯曲振动本征频率 (a) a型号换能器; (b) b型号换能器; (c) c型号换能器
Fig. 2. Effective electro-mechanical coupling coefficient, resonant frequency of transducer and eigen frequency of bending vibration of ABH part with different m: (a) Model a transducer; (b) model b transducer; (c) model c transducer.
图 3 四种不同换能器在空气中的辐射声压分布图 (a) a型号换能器, fm = 18.766 kHz; (b) b型号换能器, fm = 18.829 kHz; (c) c型号换能器, fm = 18.841 kHz; (d) d型号换能器, fm = 19.01 kHz
Fig. 3. The radiation sound pressure distribution of four different transducers in the air: (a) Model a transducer, fm = 18.766 kHz; (b) model b transducer, fm = 18.829 kHz; (c) model c transducer, fm = 18.841 kHz; (d) model d transducer, fm = 19.01 kHz.
图 4 四种不同换能器在空气中的辐射声强分布图 (a) a型号换能器, fm = 18.766 kHz; (b) b型号换能器, fm = 18.829 kHz; (c) c型号换能器, fm = 18.841 kHz; (d) d型号换能器, fm = 19.01 kHz
Fig. 4. The radiation sound intensity distribution of four different transducers in the air: (a) Model a transducer, fm = 18.766 kHz; (b) model b transducer, fm = 18.829 kHz; (c) model c transducer, fm = 18.841 kHz; (d) model d transducer, fm = 19.01 kHz.
图 5 四种换能器的声场指向性图 (a) a型号换能器, fm = 18.766 kHz; (b) b型号换能器, fm = 18.829 kHz; (c) c型号换能器, fm = 18.841 kHz; (d) d型号换能器, fm = 19.01 kHz
Fig. 5. Four types of transducer sound field directionality graphics: (a) Model a transducer, fm = 18.766 kHz; (b) model b transducer, fm = 18.829 kHz; (c) model c transducer, fm = 18.841 kHz; (d) model d transducer, fm = 19.01 kHz.
表 1 几种不同换能器的尺寸
Table 1. The sizes of several different transducers.
换能器编号 $ {r}_{1}/{\mathrm{m}}{\mathrm{m}} $ $ {r}_{2}/{\mathrm{m}}{\mathrm{m}} $ $ {r}_{3}/{\mathrm{m}}{\mathrm{m}} $ $ {r}_{4}/{\mathrm{m}}{\mathrm{m}} $ p $ {h}_{1}/{\mathrm{m}}{\mathrm{m}} $ $ {h}_{0}/{\mathrm{m}}{\mathrm{m}} $ a 30 35 40 45 3 5 1.0 b 30 35 40 45 2 5 1.5 c 30 35 40 45 2 5 2.0 d 30 35 40 45 — — — 表 2 换能器尺寸参数
Table 2. Size parameters of the transducer.
尺寸参数 r1/mm r2/mm r3/mm r4/mm p h1/mm h0/mm m 22 26 30 34 2 4 1 0.65 -
[1] 胡银丰 2017 声学与电子工程 39 25
Hu Y F 2017 Acoust Electr Eng. 39 25
[2] 曹雪砷, 王易敏, 王申 2015 测井技术 39 160Google Scholar
Cao X S, Wang Y M, Wang S 2015 Well Logging Technol. 39 160Google Scholar
[3] 刘世清, 麻磊磊 2020 陕西师范大学学报(自然科学版) 48 60
Liu S Q, Ma L L 2020 J. Shaanxi Normal Univ. (Nat. Sci. Ed. ). 48 60
[4] Mironov M A 1988 Sov. Phys. Acoustics. 34 318
[5] Krylov V V 1989 Sov. Phys. Acoustics. 35 176
[6] Krylov V V 1990 Mosc. U. Phys. B+. 45 65
[7] Krylov V V, Tilman F 2004 J. Sound Vib. 274 605Google Scholar
[8] Remillieux M C, Anderson B E, Le Bas P Y, Ulrich T J 2014 Ultrasonics 54 1409Google Scholar
[9] Anderson B E, Remillieux M C, Le Bas P Y, Ulrich T J, Pieczonka L 2015 Ultrasonics 63 141Google Scholar
[10] Aronov B S 2013 J. Acoust. Soc. Am. 134 1021Google Scholar
[11] Aronov B S 2013 J. Acoust. Soc. Am. 133 3875Google Scholar
[12] 吴德林, 戴郁郁, 陈浩, 张海澜, 王秀明 2020 声学学报 45 103Google Scholar
Wu D L, Dai Y Y, Chen H, Zhang H L, Wang X M 2020 Acta Acustica. 45 103Google Scholar
[13] 林书玉, 鲜小军 2014 陕西师范大学学报(自然科学版) 42 31
Lin S Y, Xian X J 2014 J. Shaanxi Normal Univ. (Nat. Sci. Ed. ). 42 31
[14] 许龙, 林书玉 2012 声学学报 37 408Google Scholar
Xu L, Lin S Y 2012 Acta Acustica. 37 408Google Scholar
[15] 李凤鸣, 刘世清, 许龙, 张海岛, 曾小梅, 陈赵江 2023 中国科学: 物理学 力学 天文学 53 198
Li F M, Liu S Q, Xu L, Zhang H D, Zeng X M, Chen Z J 2023 Sci. Sin. Phys. Mech. Astron. 53 198
[16] 董宜雷, 陈诚, 林书玉 2023 物理学报 72 054304Google Scholar
Dong Y L, Chen C, Lin S Y 2023 Acta Phys. Sin. 72 054304Google Scholar
[17] 高炳山, 林书玉 2010 应用声学 29 217Google Scholar
Gao B S, Lin S Y 2010 J. Appl. Acoust. 29 217Google Scholar
[18] 胡静, 林书玉, 王成会, 莫润阳 2015 西北大学学报(自然科学版) 45 709
Hu J, Lin S Y, Wang C H, Mo R Y 2015 J. Northwest Univ. (Nat. Sci. Ed.) 45 709
[19] 王帅军, 林书玉 2011 陕西师范大学学报(自然科学版) 39 23
Wang S, Lin S 2011 J. Shaanxi Normal Univ. (Nat. Sci. Ed.) 39 23
[20] Krylov V V 2020 J. Sound Vib. 468 115100Google Scholar
[21] 林书玉 2004 超声换能器的原理及设计 (北京: 科学出版社) 第20页
Lin S Y 2004 The Theory and Design of Ultrasonic Transducers (Beijing: Science Press) p20
期刊类型引用(1)
1. 陈诚,刘洋,郭建中,林书玉. 弯曲声黑洞结构纵弯耦合振动超声手术刀. 中国科学:物理学 力学 天文学. 2024(12): 97-110 . 百度学术 其他类型引用(0)
计量
- 文章访问数: 2238
- PDF下载量: 60
- 被引次数: 1