-
径向振动压电超声换能器具有辐射面积大、沿半径方向有全指向性、接收灵敏度高、结构简单且性能稳定等优点[1-4], 在超声医学成像[5,6]、超声化学反应[7,8]、超声污水处理[9,10]、水声通信[11,12]等领域被广泛应用. 径向复合压电超声换能器由金属圆环和压电陶瓷圆环复合而成, 金属圆环可以增加压电陶瓷的散热, 提高其功率容量. 为了增大径向复合压电超声换能器输出端的位移幅值, 实现能量集中, 输出端金属圆环多采用变厚度金属圆环. 变厚度金属圆环可置于压电陶瓷环内以实现向内聚能, 也可置于压电陶瓷外从而实现大面积声能辐射. 径向复合压电超声换能器根据金属圆环轴向厚度沿半径变化的函数关系被分成锥型、指数型、幂函数型、阶梯型等[13-16].
对于变厚度环型压电超声换能器的设计和振动分析, 通常使用有限元分析法[17]和等效电路法[18,19]. 等效电路法通过力电类比原理将力学量类比为电学量从而把振动问题简化. 传统理论研究变厚度环型压电超声换能器通常需求解变厚度金属薄圆环径向振动的波动方程. 基于变厚度径向振动的波动方程, 许龙等[20,21]研究了幂函数型环型聚能器的径向振动和阶梯圆环压电超声换能器的径向振动; 刘世清等[22]研究了线性变厚度环型聚能器的径向振动; 王晓宇和林书玉[23]研究了线性变厚度径向变幅器的耦合振动. 但是对于复杂函数变厚度如指数型、悬链线型环型压电超声换能器, 波动方程求解比较困难. 传输矩阵法可以将连续结构分解成一系列微元结构, 一对微元之间通过传输矩阵衔接, 从而将复杂问题简单化. 基于此, 本文提出了用传输矩阵法[24]将任意变厚度金属圆环的径向振动等效近似为
$ N $ 个等厚度圆环径向振动叠加的分析方法, 得到了任意变厚度金属圆环径向振动的共振频率方程和位移放大系数表达式, 以及任意变厚度环型压电超声换能器的共振频率方程. 研究了锥型变厚度、幂函数型变厚度、指数型变厚度、悬链线型变厚度金属圆环的共振频率和位移放大系数与几何尺寸的关系, 以及不同尺寸下锥型变厚度环型压电换能器径向一阶、二阶共振频率和反共振频率. -
如图1所示, 变厚度环型压电超声换能器由变厚度金属内圆环和等厚度压电陶瓷外圆环组成, 变厚度金属圆环的内环半径为
$ {R}_{{\rm{b}}} $ 、外环半径为$ {R}_{{\rm{a}}} $ 、内环半径处轴向厚度为$ {h}_{{\rm{b}}} $ 、外环半径处轴向厚度为$ {h}_{{\rm{a}}} $ , 等厚度压电陶瓷圆环的轴向厚度为$ {h}_{{\rm{a}}} $ 、内环半径为$ {R}_{{\rm{a}}} $ 、外环半径为$ {R}_{{\rm{c}}} $ . 假设变厚度环型压电超声换能器的轴向厚度尺寸远小于径向尺寸, 忽略其纵向振动所带来的耦合效应, 此时环型压电超声换能器做纯径向振动. 环型压电超声换能器的径向振动由金属内圆环径向振动和压电陶瓷外圆环径向振动复合而成, 因此需要分别研究变厚度金属圆环的径向振动和等厚度压电陶瓷圆环的径向振动. -
根据轴向厚度沿半径的变化函数, 变厚度金属薄圆环可以分成锥型、阶梯型、幂函数型等. 为了求解任意变厚度金属薄圆环径向振动的等效电路, 本文采用传输矩阵法, 将变厚度金属薄圆环沿半径方向等分成
$ N $ 个径向厚度为$ \Delta r $ 的薄圆环. 当$ N $ 足够大、薄圆环的径向厚度$ \Delta r $ 的值足够小时, 薄圆环可以看成是轴向厚度不变的等厚度薄圆环. 此时, 变厚度金属薄圆环径向振动的等效电路可以由$ N $ 个等厚度金属薄圆环径向振动的等效电路串联而成.根据文献[25]可知, 等厚度金属薄圆环径向振动的等效电路如图2所示, 图中
图 2 等厚度金属薄圆环径向振动的等效电路图
Figure 2. Equivalent circuit of metal thin circular annular in radial vibration.
$ \begin{split} {Z}_{11}=\;& {\rm{j}}\frac{2{Z}_{2}}{{\rm{\pi }}k{R}_{2}\left[{{\rm{J}}}_{1}\right(k{R}_{1}\left){{\rm{Y}}}_{1}\right(k{R}_{2})-{{\rm{J}}}_{1}(k{R}_{2}\left){{\rm{Y}}}_{1}\right(k{R}_{1})} \\ & \times \left[\frac{{{\rm{J}}}_{1}\left(k{R}_{2}\right){{\rm{Y}}}_{0}\left(k{R}_{1}\right)-{{\rm{J}}}_{0}\left(k{R}_{1}\right){{\rm{Y}}}_{1}\left(k{R}_{2}\right)-{{\rm{J}}}_{1}\left(k{R}_{1}\right){{\rm{Y}}}_{0}\left(k{R}_{1}\right)+{{\rm{J}}}_{0}\left(k{R}_{1}\right){{\rm{Y}}}_{1}\left(k{R}_{1}\right)}{{{\rm{J}}}_{1}\left(k{R}_{1}\right){{\rm{Y}}}_{0}\left(k{R}_{1}\right)-{{\rm{J}}}_{0}\left(k{R}_{1}\right){{\rm{Y}}}_{1}\left(k{R}_{1}\right)}\right] \\ & +{\rm{j}}\frac{2{Z}_{1}(1-\nu )}{{\rm{\pi }}(k{R}_{1}{)}^{2}\left[{{\rm{J}}}_{1}\left(k{R}_{1}\right){{\rm{Y}}}_{0}\left(k{R}_{1}\right)-{{\rm{J}}}_{0}\left(k{R}_{1}\right){{\rm{Y}}}_{1}\left(k{R}_{1}\right)\right]}, \end{split} $ $ \begin{split} {Z}_{12}=\;&{\rm{j}}\frac{2{Z}_{2}}{{\rm{\pi }}k{R}_{2}\left[{{\rm{J}}}_{1}\right(k{R}_{1}\left){{\rm{Y}}}_{1}\right(k{R}_{2})-{{\rm{J}}}_{1}(k{R}_{2}\left){{\rm{Y}}}_{1}\right(k{R}_{1})} \\ &\times \left[\frac{{{\rm{J}}}_{1}\left(k{R}_{1}\right){{\rm{Y}}}_{0}\left(k{R}_{2}\right)-{{\rm{J}}}_{0}\left(k{R}_{2}\right){{\rm{Y}}}_{1}\left(k{R}_{1}\right)-{{\rm{J}}}_{1}\left(k{R}_{2}\right){{\rm{Y}}}_{0}\left(k{R}_{2}\right)+{{\rm{J}}}_{0}\left(k{R}_{2}\right){{\rm{Y}}}_{1}\left(k{R}_{2}\right)}{{{\rm{J}}}_{1}\left(k{R}_{2}\right){{\rm{Y}}}_{0}\left(k{R}_{2}\right)-{{\rm{J}}}_{0}\left(k{R}_{2}\right){{\rm{Y}}}_{1}\left(k{R}_{2}\right)}\right] \\ & -{\rm{j}}\frac{2{Z}_{2}\left(1-\nu \right)}{{\rm{\pi }}(k{R}_{2}{)}^{2}\left[{{\rm{J}}}_{1}\left(k{R}_{2}\right){{\rm{Y}}}_{0}\left(k{R}_{2}\right)-{{\rm{J}}}_{0}\left(k{R}_{2}\right){{\rm{Y}}}_{1}\left(k{R}_{2}\right)\right]}, \end{split} $ $ {Z}_{13}={\rm{j}}\frac{2{Z}_{2}}{{\rm{\pi }}k{R}_{2}\left[{{\rm{J}}}_{1}\right(k{R}_{1}\left){{\rm{Y}}}_{1}\right(k{R}_{2})-{{\rm{J}}}_{1}(k{R}_{2}\left){{\rm{Y}}}_{1}\right(k{R}_{1})}, $ 其中
$ {Z}_{1}=\rho c{S}_{1} $ ,$ {Z}_{2}=\rho c{S}_{2} $ ,$ {S}_{1}=2{\text{π}}{R}_{1}{h}_{1} $ ,$ {S}_{2}= 2{\text{π}}{R}_{2}{h}_{1} $ ,$ {R}_{1} $ ,${R}_{2}{,\; h}_{1}$ 分别是金属薄圆环的外半径、内半径和轴向厚度,$ \rho $ 是金属薄圆环的密度,$\nu$ 是金属薄圆环的泊松比,$ c $ 是金属薄圆环内的声速,$ k $ 是波数,$ {\rm{J}} $ 是贝塞尔函数,$ {\rm{Y}} $ 是诺依曼函数.图2中
$ {F}_{1} $ ,$ {V}_{1} $ ,$ {F}_{2},\; {V}_{2} $ 分别是第一个薄圆环输入端与输出端的径向受力和径向振动速度. 根据等厚度金属薄圆环径向振动的机械等效电路图和基尔霍夫定律可得:$ \left\{\begin{aligned} & {F}_{1}=\left({Z}_{11}+{Z}_{13}\right){V}_{1}-{Z}_{13}{V}_{2}\text{, }\\ & {F}_{2}=-\left({Z}_{12}+{Z}_{13}\right){V}_{2}+{Z}_{13}{V}_{1}. \end{aligned}\right. $ 将输入端的力学参量
$ {F}_{1} $ 和径向振动速度$ {V}_{1} $ 移到等号左边, 输出端的力学参数$ {F}_{2} $ 和径向振动速度$ {V}_{2} $ 移到等号右边, (4)式可以改写为$ \left[\begin{array}{c}{F}_{1}\\ {V}_{1}\end{array}\right]={\boldsymbol{M}}_{1}\left[\begin{array}{c}{F}_{2}\\ {V}_{2}\end{array}\right] , $ 其中
$ {\boldsymbol{M}}_{1}=\left[\begin{array}{cc}1+\dfrac{{Z}_{11}}{{Z}_{13}}& \text{}{Z}_{11}+{Z}_{12}+\dfrac{{Z}_{11}{Z}_{12}}{{Z}_{13}}\\ \dfrac{1}{{Z}_{13}}& 1+\dfrac{{Z}_{12}}{{Z}_{13}}\end{array}\right] . $ 由于变厚度金属薄圆环的径向应力
$ F $ 和径向振动速度$ V $ 连续, 所以第2个薄圆环输入端的径向应力和径向振动速度等于第1个薄圆环输出端的力学参数$ {F}_{2} $ 和径向振动速度$ {V}_{2} $ . 类比第1个薄圆环输入输出端$ F $ ,$ V $ 之间的关系, 第2个薄圆环输入端参数$ {F}_{2} $ 和$ {V}_{2} $ 可以由输出端参数$ {F}_{3} $ 和$ {V}_{3} $ 表示为$ \left[\begin{array}{c}{F}_{2}\\ {V}_{2}\end{array}\right]={\boldsymbol{M}}_{2}\left[\begin{array}{c}{F}_{3}\\ {V}_{3}\end{array}\right]. $ 进而第1个薄圆环输入端的力参数
$ {F}_{1} $ 和振动速度$ {V}_{1} $ 可以由第2个薄圆环输出端参数$ {F}_{3} $ 和$ {V}_{3} $ 表示为$ \left[\begin{array}{c}{F}_{1}\\ {V}_{1}\end{array}\right]={\boldsymbol{M}}_{1}{\boldsymbol{M}}_{2}\left[\begin{array}{c}{F}_{3}\\ {V}_{3}\end{array}\right] . $ 以此类推, 若将变厚度金属薄圆环分成
$ N $ 个等厚度薄圆环串联, 则第1个薄圆环输入端参数$ {F}_{1} $ ,$ {V}_{1} $ 可以由第$ N $ 个金属薄圆环输出端力参数$ {F}_{N+1} $ 和振动速度$ {V}_{N+1} $ 表示为$ \left[\begin{array}{c}{F}_{1}\\ {V}_{1}\end{array}\right]={\boldsymbol{M}}_{1}{\boldsymbol{M}}_{2}\cdot \cdot \cdot {\boldsymbol{M}}_{i}\cdots {\boldsymbol{M}}_{N}\left[\begin{array}{c}{F}_{N+1}\\ {V}_{N+1}\end{array}\right] , $ 其中
$ {\boldsymbol{M}}_{i}=\left[\begin{array}{cc}1+\dfrac{{Z}_{i1}}{{Z}_{i3}}& \text{}{Z}_{i1}+{Z}_{i2}+\dfrac{{Z}_{i1}{Z}_{i2}}{{Z}_{i3}}\\ \dfrac{1}{{Z}_{i3}}& 1+\dfrac{{Z}_{i2}}{{Z}_{i3}}\end{array}\right], $ $ {\boldsymbol{M}}_{1}{\boldsymbol{M}}_{2}\cdots {\boldsymbol{M}}_{i}\cdots {\boldsymbol{M}}_{N}=\left[\begin{array}{cc}{M}_{11}& {M}_{12}\\ {M}_{21}& {M}_{22}\end{array}\right] . $ 由于
$ {F}_{1}={F}_{{\rm{a}}} $ ,$ {V}_{1}={V}_{{\rm{a}}} $ ,$ {F}_{N+1}={F}_{{\rm{b}}} $ ,$ {V}_{N+1}={V}_{{\rm{b}}} $ , 则(9)式可以改写为$ \left[\begin{array}{c}{F}_{{\rm{a}}}\\ {V}_{{\rm{a}}}\end{array}\right]=\left[\begin{array}{cc}{M}_{11}& {M}_{12}\\ {M}_{21}& {M}_{22}\end{array}\right]\left[\begin{array}{c}{F}_{{\rm{b}}}\\ {V}_{{\rm{b}}}\end{array}\right]. $ 将变厚度金属薄圆环内外半径处的力学参量
$ {F}_{{\rm{a}}} $ ,$ {F}_{{\rm{b}}} $ 移至等号左边, 径向振动速度$ {V}_{{\rm{a}}} $ ,$ {V}_{{\rm{b}}} $ 移动到等号右边, (12)式可以改写为$ \left\{\begin{aligned} & {F}_{{\rm{a}}}={Z}_{n11}{V}_{a}+{Z}_{n12}{V}_{{\rm{b}}},\\ &{F}_{{\rm{b}}}={Z}_{n21}{V}_{a}+{Z}_{n22}{V}_{{\rm{b}}},\end{aligned}\right. $ 其中
${Z}_{n11}=\dfrac{{M}_{11}}{{M}_{21}}$ ,${Z}_{n12}={M}_{12}-\dfrac{{M}_{11}{M}_{22}}{{M}_{21}}$ ,${Z}_{n21}= {1}/{{M}_{21}}$ ,${Z}_{n22}=- {{M}_{22}}/{{M}_{21}}$ .根据电力类比原理以及非互易的二端口等效网络知识, 由(13)式可得到变厚度金属薄圆环径向振动的等效电路图如图3所示.
图 3 变厚度金属薄圆环径向振动的等效电路图
Figure 3. Equivalent circuit diagram of metal thin annular with variable thickness in radial vibration.
当变厚度金属薄圆环的输出端自由, 即
$ {F}_{{\rm{b}}}=0 $ 时, 变厚度金属薄圆环径向振动的共振频率方程为$ {Z}_{{\rm{i}}{\rm{n}}}=\frac{{Z}_{n11}{Z}_{n22}-{Z}_{n12}{Z}_{n21}}{{Z}_{n22}}=0. $ 此时, 变厚度金属薄圆环的位移放大系数为
$ {M}^{*}=\left|\frac{1}{{M}_{22}}\right| . $ -
根据任意变厚度金属薄圆环径向振动的等效电路图, 得到其共振频率方程和位移放大系数表达式. 选取N=100, 使用数值计算软件MATLAB根据(14)式和(15)式分别计算①锥型金属圆环、②幂函数型金属圆环、③指数型金属圆环和④悬链线型金属圆环径向振动的共振频率和位移放大系数. 锥型金属圆环的厚度随半径以线性函数规律变化, 幂函数型金属圆环的厚度随半径以幂函数
$ h\left(r\right)=m/{r}^{n} $ 规律变化, 指数型金属圆环的厚度随半径以指数函数$ h\left(r\right)=a{{\rm{e}}}^{-\beta r} $ 规律变化, 悬链线型金属圆环的厚度随半径以悬链线函数$ h\left(r\right)= {h}_{{\rm{a}}}{\rm{c}}{\rm{o}}{\rm{s}}{\rm{h}}\left[\gamma \right({R}_{{\rm{a}}}-r\left)\right] $ 规律变化. 金属圆环的尺寸参数是: 内环半径$ {R}_{{\rm{b}}}=25{\rm{ }}{\rm{m}}{\rm{m}} $ , 外环半径$ {R}_{{\rm{a}}}=50{\rm{ }}{\rm{m}}{\rm{m}} $ , 外环半径处厚度$ {h}_{{\rm{a}}}=10{\rm{ }}{\rm{m}}{\rm{m}} $ , 内环半径处厚度$ {h}_{{\rm{b}}} $ 的取值范围为$ 1—10{\rm{ }}{\rm{m}}{\rm{m}} $ . 选取金属圆环的材料为铝, 其材料参数为: 密度$ \rho =2700{\rm{ }}{\rm{k}}{\rm{g}}/{{\rm{m}}}^{3} $ , 杨氏模量$ E=70~{\rm{G}}{\rm{P}}{\rm{a}} $ , 泊松比$ \upsilon =0.33 $ . 为了与理论计算结果进行对比, 使用有限元软件COMSOL6.0对上述材料尺寸的金属圆环进行仿真模拟. 4种变厚度金属圆环的一阶、二阶径向共振频率和位移放大系数的理论解和数值解如表1和表2所示. 表1和表2中$ f{\rm{和}}{f}^{*} $ 分别表示变厚度金属圆环径向振动共振频率的理论解和数值解,$ {M}^{*} $ ,$ {M}^{**} $ 分别表示变厚度金属圆环位移放大系数的理论解和数值解, 相对误差计算方法为:${\varDelta }_{f}= \left|(f-{f}^{*})/{f}^{*}\right|$ ,${\varDelta }_{{M}^{*}}= \left|({M}^{*}-{M}^{**})/{M}^{**}\right|$ .$ {h}_{{\rm{b}}}/{\rm{m}}{\rm{m}} $ $ {h}_{{\rm{a}}}/{\rm{m}}{\rm{m}} $ $ {f}_{{\rm{r}}1}/{\rm{H}}{\rm{z}} $ ${ {f}^{ {\rm{*} } }_{ {\rm{r} }1} }/{\rm{H} }{\rm{z} }$ ${\varDelta }_{ {f}_{ {\rm{r} }1} }/{\rm{\%} }$ $ {f}_{{\rm{r}}2}/{\rm{H}}{\rm{z}} $ ${ {f}^{*}_{r2} }/{\rm{H} }{\rm{z} }$ ${\varDelta }_{ {f}_{ {\rm{r} }2} }/{\rm{\%} }$ 锥型 6 10 21894 21892 0.01 112380 110790 1.44 幂函数型 5 10 21679 21679 0 113330 111720 1.44 指数型 4 10 21401 21399 0.01 112770 111480 1.16 悬链线型 3 10 20968 20957 0.05 109700 108890 0.74 表 1 变厚度金属圆环径向一阶、二阶共振频率
Table 1. Radial first and second order resonance frequencies of metal rings with variable thickness.
$ {h}_{{\rm{b}}}/{\rm{m}}{\rm{m}} $ $ {h}_{{\rm{a}}}/{\rm{m}}{\rm{m}} $ $ {M}_{{\rm{r}}1}^{{\rm{*}}} $ $ {M}_{{\rm{r}}1}^{{\rm{*}}{\rm{*}}} $ ${\varDelta }_{ {M}_{ {\rm{r} }1}^{*} }$/% $ {M}_{{\rm{r}}2}^{{\rm{*}}} $ $ {M}_{{\rm{r}}2}^{{\rm{*}}{\rm{*}}} $ ${\varDelta }_{ {M}_{ {\rm{r} }2}^{*} }/{\rm{\%} }$ 锥型 6 10 1.1985 1.1987 0.02 1.8304 1.8160 0.79 幂函数型 5 10 1.1997 1.2000 0.02 1.9876 1.9670 1.04 指数型 4 10 1.1992 1.1995 0.02 2.2501 2.2135 1.16 悬链线型 3 10 1.1959 1.1956 0.02 2.8118 2.7359 2.77 表 2 变厚度金属圆环径向一阶、二阶共振位移放大系数
Table 2. Radial first and second order resonance displacement amplification coefficients of metal rings with variable thickness.
4种变厚度金属圆环(①锥型金属圆环、②幂函数型金属圆环、③指数型金属圆环、④悬链线型金属圆环)径向一阶、二阶共振频率和位移放大系数随内外环厚度比
$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的变化曲线如图4和图5所示, 其中-T为理论结果, -FEM为有限元仿真结果.图 4 4种变厚度金属圆环一阶、二阶径向共振频率与
$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的关系 (a)一阶径向共振; (b)二阶径向共振Figure 4. The relationship between the first and second order radial resonance frequencies of four kinds of variable thickness metal rings and thickness ratio
$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ : (a) First-order radial resonance; (b) second-order radial resonance.图 5 4种变厚度环型聚能器一阶和二阶径向共振位移放大系数随
$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的变化曲线 (a)一阶径向共振; (b)二阶径向共振Figure 5. The relationship between the first and second order radial resonance displacement amplification coefficient and thickness ratio
$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ of four kinds of variable thickness metal rings: (a) First-order radial resonance; (b) second-order radial resonance.由图4可知, 4种变厚度金属圆环的一阶共振频率随
$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的增大而减小, 当内环半径处的厚度$ {h}_{{\rm{b}}} $ 减小时, 金属圆环等效半径增大, 共振频率减小. 二阶共振频率随着$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的增大而增大, 径向二阶共振频率受金属圆环等效质量影响较大, 当内环半径处厚度$ {h}_{{\rm{b}}} $ 减小时, 金属圆环等效质量减小, 共振频率向高频偏移. 另外, 4种变厚度金属圆环一阶共振频率的数值解几乎在理论解随$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的变化曲线上, 二阶共振频率的数值解随$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的变化趋势与理论解一致, 且相对误差在2%以内, 这表明传输矩阵法可以用来研究任意变厚度金属圆环的径向振动. 由于径向二阶的声波波长比一阶小, 耦合振动效应较为明显, 使得径向二阶共振频率的相对误差大于径向一阶.由图5可知, 4种变厚度金属圆环(①锥型金属圆环、②幂函数型金属圆环、③指数型金属圆环和④悬链线型金属圆环)一阶共振、二阶共振的位移放大系数随
$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的增大而增大, 这表明内环半径处的厚度$ {h}_{{\rm{b}}} $ 越小, 环型聚能器向内聚能的效果越好. 整体来看, 一阶共振的位移放大系数从大到小依次是锥型、幂函数型、指数型、悬链线型, 二阶共振的位移放大系数从大到小依次是悬链线型、指数型、幂函数型、锥型. 4种聚能器的一阶共振位移放大系数较小, 在1.2左右, 二阶共振位移放大系数较大, 在2—4范围. 这表明位移放大系数不止与厚度随半径的函数变化关系有关, 还与工作模态相关. 另外, 4种聚能器径向一阶、二阶共振位移放大系数的数值解基本与理论解随$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的变化曲线重合, 验证了理论解的正确性. -
本文将变厚度金属圆环的径向振动等效近似为N个等厚金属圆环径向振动的叠加, N的取值大小对计算精度影响较大. 选取锥型变厚度金属圆环的材料为铝, 尺寸参数是: 内环半径
$ {R}_{{\rm{b}}}=25{\rm{ }}{\rm{m}}{\rm{m}} $ , 外环半径$ {R}_{{\rm{a}}}=50{\rm{ }}{\rm{m}}{\rm{m}} $ , 外环半径处厚度$ {h}_{{\rm{a}}}= 10\;{\rm{mm}} $ , 内环半径处厚度$ {h}_{{\rm{b}}}=1~{\rm{m}}{\rm{m}} $ . 锥型变厚度金属圆环的一阶、二阶径向共振频率随N的变化趋势如图6所示, 由图6可知, 当N的取值大于50时, 锥型金属圆环径向一阶、二阶共振频率趋于稳定, 即等厚度金属圆环的径向厚度小于$ 0.5\;{\rm{m}}{\rm{m}} $ 时, 共振频率的计算精度满足要求. -
根据文献[26]可知, 在厚度方向极化的压电陶瓷薄圆环径向振动的等效电路图如图7所示. 图7中
$ {F}_{{\rm{a}}}, {V}_{{\rm{a}}} $ 和$ {F}_{{\rm{c}}}, {V}_{{\rm{c}}} $ 分别代表压电陶瓷内环和外环的径向作用力和径向振动速度,$ {n}_{1} = {\text{π}}{k}_{{\rm{p}}}{R}_{{\rm{c}}}/2 $ ,${n}_{2} = {\text{π}}{k}_{{\rm{p}}}{R}_{{\rm{a}}}/2, {k}_{{\rm{p}}} = {\omega }_{{\rm{p}}}/{c}_{{\rm{p}}}, {\omega }_{{\rm{p}}} = 2{\text{π}}f , {c}_{{\rm{p}}} = [{S}_{11}^{{\rm{E}}}{\rho}_{{\rm{p}}} (1 - {\nu }_{{\rm{p}}}^{2})]^{-1/2}$ ,${\nu }_{{\rm{p}}} = -{S}_{12}^{{\rm{E}}}/{S}_{11}^{{\rm{E}}}$ ,${N}_{31}={{\text{π}}}^{2}{k}_{{\rm{p}}}{R}_{{\rm{a}}}{R}_{c}{d}_{31}/ ({s}_{11}^{{\rm{E}}} + {s}_{12}^{{\rm{E}}})$ ,${C}_{0{\rm{r}}} = ( {{\varepsilon }_{33}^{{\rm{T}}}{S}_{{\rm{p}}}}/{{h}_{{\rm{a}}}}) \{1 - 2{d}_{31}^{2}/[{\varepsilon }_{33}^{{\rm{T}}}({s}_{11}^{{\rm{E}}} + {s}_{12}^{{\rm{E}}})]\}$ ,${S}_{{\rm{p}}} = {\text{π}}({R}_{{\rm{c}}}^{2} - {R}_{{\rm{a}}}^{2}),$ $ {N}_{31} $ 和$ {C}_{0{\rm{r}}} $ 分别代表压电陶瓷圆环径向振动的机电耦合系数和钳定电容,$ {d}_{31} $ 是压电应变常数,$ {s}_{11}^{{\rm{E}}} $ 和$ {s}_{12}^{{\rm{E}}} $ 是恒定电场下的弹性柔顺系数,$ {\varepsilon }_{33}^{{\rm{T}}} $ 是恒定应力下的介电常数. 图7中机械阻抗$ {Z}_{{\rm{P}}1} $ ,$ {Z}_{{\rm{P}}2} $ ,$ {Z}_{{\rm{P}}3} $ 的具体表达式如下:图 7 压电陶瓷圆环径向振动的等效电路图
Figure 7. Equivalent circuit diagram of piezoelectric ceramic annular in radial vibration.
$ \begin{split} {Z}_{{\rm{P}}1}=& \frac{{{\rm{\pi }}}^{2}{\left({k}_{{\rm{p}}}{R}_{{\rm{c}}}\right)}^{2}{Z}_{{\rm{r}}{\rm{a}}}}{4{\rm{j}}}\left[\frac{{{\rm{Y}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{c}}}\right){{\rm{J}}}_{0}\left({k}_{{\rm{p}}}{R}_{{\rm{a}}}\right)-{{\rm{J}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{c}}}\right){{\rm{Y}}}_{0}\left({k}_{{\rm{p}}}{R}_{{\rm{a}}}\right)}{{{\rm{J}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{c}}}\right){{\rm{Y}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{a}}}\right)-{{\rm{J}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{a}}}\right){{\rm{Y}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{c}}}\right)}+\frac{1-{\nu }_{{\rm{p}}}}{{k}_{{\rm{p}}}{R}_{{\rm{a}}}}\right] \\ & -{\rm{j}}\frac{{Z}_{{\rm{r}}{\rm{a}}}}{2}\frac{{\rm{\pi }}{k}_{{\rm{p}}}{R}_{{\rm{c}}}}{{{\rm{J}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{c}}}\right){{\rm{Y}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{a}}}\right)-{{\rm{J}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{a}}}\right){{\rm{Y}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{c}}}\right)}\text{, } \end{split} $ $ \begin{split} {Z}_{{\rm{P}}2}=& \frac{{{\rm{\pi }}}^{2}{\left({k}_{{\rm{p}}}{R}_{{\rm{a}}}\right)}^{2}{Z}_{{\rm{r}}{\rm{c}}}}{4{\rm{j}}}\left[\frac{{{\rm{Y}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{a}}}\right){{\rm{J}}}_{0}\left({k}_{{\rm{p}}}{R}_{{\rm{c}}}\right)-{\rm{J}}\left({k}_{{\rm{p}}}{R}_{{\rm{a}}}\right){{\rm{Y}}}_{0}\left({k}_{{\rm{p}}}{R}_{{\rm{c}}}\right)}{{{\rm{J}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{c}}}\right){{\rm{Y}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{a}}}\right)-{{\rm{J}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{a}}}\right){{\rm{Y}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{c}}}\right)}+\frac{1-{\nu }_{{\rm{p}}}}{{k}_{{\rm{p}}}{R}_{{\rm{c}}}}\right] \\ & -{\rm{j}}\frac{{Z}_{{\rm{r}}{\rm{c}}}}{2}\frac{{\rm{\pi }}{k}_{{\rm{p}}}{R}_{{\rm{a}}}}{{{\rm{J}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{c}}}\right){{\rm{Y}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{a}}}\right)-{{\rm{J}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{a}}}\right){{\rm{Y}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{c}}}\right)}, \end{split} $ $ {Z}_{{\rm{P}}3}={\rm{j}}\frac{{Z}_{{\rm{r}}{\rm{a}}}}{2}\frac{{\rm{\pi }}{k}_{{\rm{p}}}{R}_{{\rm{c}}}}{{{\rm{J}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{c}}}\right){{\rm{Y}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{a}}}\right)-{{\rm{J}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{a}}}\right){{\rm{Y}}}_{1}\left({k}_{{\rm{p}}}{R}_{{\rm{c}}}\right)}\text{, } $ 其中,
$ {Z}_{{\rm{r}}{\rm{a}}}={\rho }_{{\rm{p}}}{c}_{{\rm{p}}}{S}_{{\rm{r}}{\rm{a}}} $ ,$ {Z}_{{\rm{r}}{\rm{c}}}={\rho }_{{\rm{p}}}{c}_{{\rm{p}}}{S}_{{\rm{r}}{\rm{c}}} $ ,$ {S}_{{\rm{r}}{\rm{a}}}=2{\text{π}}{R}_{{\rm{a}}}{h}_{{\rm{a}}} $ ,$ {S}_{{\rm{r}}{\rm{c}}}=2{\text{π}}{R}_{{\rm{c}}}{h}_{{\rm{a}}} $ ,$ {h}_{{\rm{a}}} $ 是压电陶瓷圆环的厚度,$ {k}_{{\rm{p}}} $ 和$ {c}_{{\rm{p}}} $ 分别是压电陶瓷圆环的波数和径向声速.变厚度金属薄圆环与压电陶瓷薄圆环在半径
$ {R}_{{\rm{a}}} $ 处满足径向应力连续和振动速度连续, 因此, 变厚度金属薄圆环径向振动与压电陶瓷薄圆环径向振动的等效电路串联, 可得到环型压电超声换能器径向振动的等效电路图如图8所示.图 8 环型压电超声换能器径向振动的等效电路图
Figure 8. Equivalent circuit diagram of circular piezoelectric ultrasonic transducer in radial vibration.
当环型压电超声换能器内外表面所受径向作用力为0时, 即图8中
$ {F}_{{\rm{b}}}={F}_{{\rm{c}}}=0 $ 时, 环型压电超声换能器电端输入阻抗表达式为$ {Z}_{{\rm{e}}}=\frac{1/\left({\rm{j}}\omega {C}_{0{\rm{r}}}\right)\cdot {Z}_{{\rm{m}}}/{N}_{31}^{2}}{1/\left({\rm{j}}\omega {C}_{0{\rm{r}}}\right)+{Z}_{{\rm{m}}}/{N}_{31}^{2}}=\frac{{Z}_{{\rm{m}}}}{{\rm{j}}\omega {C}_{0{\rm{r}}}{Z}_{{\rm{m}}}+{N}_{31}^{2}}, $ 其中
$ {Z}_{{\rm{m}}}=\frac{{(n}_{1}^{2}{Z}_{{\rm{i}}{\rm{n}}}+{Z}_{{\rm{p}}1}){Z}_{{\rm{p}}2}}{{n}_{1}^{2}{Z}_{{\rm{i}}{\rm{n}}}+{Z}_{{\rm{p}}1}+{Z}_{{\rm{p}}2}}+{Z}_{{\rm{p}}3}. $ 因此环型压电超声换能器径向振动共振频率方程为
$ {Z}_{{\rm{m}}}=0, $ 径向振动反共振频率方程为
$ {Z}_{{\rm{m}}}=\infty . $ -
根据(21)式和(22)式可以求得任意变厚度环型压电超声换能器径向振动的共振频率和反共振频率, 本文以锥型变厚度环型压电超声换能器为例, 锥型压电超声换能器由锥型变厚度金属薄圆环和等厚度压电陶瓷圆环组成. 变厚度金属薄圆环的尺寸参数是: 内环半径
$ {R}_{{\rm{b}}}=20{\rm{ }}{\rm{m}}{\rm{m}} $ , 外环半径${R}_{{\rm{a}}}= 40{\rm{ }}{\rm{m}}{\rm{m}}$ , 外环半径处厚度$ {h}_{{\rm{a}}}=10{\rm{ }}{\rm{m}}{\rm{m}} $ , 内环半径处厚度$ {h}_{{\rm{b}}} $ 的取值范围为$ 1—10{\rm{ }}{\rm{m}}{\rm{m}} $ . 压电陶瓷圆环的尺寸参数是: 内环半径$ {R}_{{\rm{a}}}=40{\rm{ }}{\rm{m}}{\rm{m}} $ , 外环半径$ {R}_{{\rm{c}}}= 45{\rm{ }}{\rm{m}}{\rm{m}} $ , 压电陶瓷圆环的厚度$ {h}_{{\rm{a}}}=10{\rm{ }}{\rm{m}}{\rm{m}} $ . 选取金属圆环的材料为铝, 选取压电陶瓷圆环的材料为PZT-4, 其材料参数是:$ {\rho }_{{\rm{p}}}=7500{\rm{ }}{\rm{k}}{\rm{g}}/{{\rm{m}}}^{3} $ ,${S}_{11}^{{\rm{E}}}=12.3\times {10}^{-12}{{\rm{m}}}^{2}/{\rm{N}}$ ,${S}_{12}^{{\rm{E}}} = -4.05 \times {10}^{-12}{{\rm{m}}}^{2}/{\rm{N}}$ ,${d}_{31}=-1.23\times {10}^{-10}{\rm{C}}/{\rm{N}}$ ,$ {\varepsilon }_{33}^{{\rm{T}}}=11.4946\times {10}^{-9} $ . 为了与理论计算结果进行对比, 使用有限元软件COMSOL6.0对上述材料尺寸的环型压电超声换能器进行仿真模拟, 模拟所得一阶共振、二阶共振的振型图如图9所示. 部分尺寸锥型变厚度环型压电超声换能器一阶、二阶共振频率和反共振频率的理论解和数值解如表3和表4所示.图 9 锥型变厚度环型压电超声换能器一阶、二阶径向共振振型图 (a)一阶径向共振; (b)二阶径向共振
Figure 9. First-order and second-order radial resonance mode shapes of conical variable thickness annular piezoelectric ultrasonic transducer: (a) First-order radial resonance; (b) second-order radial resonance.
$ {h}_{{\rm{b}}}/{\rm{m}}{\rm{m}} $ $ {h}_{{\rm{a}}}/{\rm{m}}{\rm{m}} $ $ {f}_{{\rm{r}}1}/{\rm{H}}{\rm{z}} $ ${ {f}^{*}_{ {\rm{r} }1} }/{\rm{H} }{\rm{z} }$ ${\varDelta }_{ {f}_{ {\rm{r} }1} }/$% $ {f}_{{\rm{a}}1}/{\rm{H}}{\rm{z}} $ ${ {f}^{*}_{ {\rm{a} }1} }/{\rm{H} }{\rm{z} }$ ${\varDelta }_{ {f}_{ {\rm{a} }1} }/$% $ {K}_{{\rm{e}}{\rm{f}}{\rm{f}}1} $ $ {K}_{{\rm{e}}{\rm{f}}{\rm{f}}1}^{{\rm{*}}} $ 9 10 22002 21989 0.06 22276 22269 0.03 0.156 0.158 6 10 21059 21042 0.08 21355 21348 0.03 0.166 0.169 3 10 19886 19856 0.15 20209 20192 0.08 0.178 0.182 表 3 锥型变厚度环型压电换能器一阶共振频率和反共振频率
Table 3. The first-order resonant frequency and anti-resonant frequency of conical variable thickness annular piezoelectric transducer.
$ {h}_{{\rm{b}}}/{\rm{m}}{\rm{m}} $ $ {h}_{{\rm{a}}}/{\rm{m}}{\rm{m}} $ $ {f}_{{\rm{r}}2}/{\rm{H}}{\rm{z}} $ ${ {f}^{*}_{ {\rm{r} }2} }/{\rm{H} }{\rm{z} }$ ${\varDelta }_{ {f}_{ {\rm{r} }2} }/$% $ {f}_{{\rm{a}}2}/{\rm{H}}{\rm{z}} $ ${ {f}^{*}_{ {\rm{a} }2} }/{\rm{H} }{\rm{z} }$ ${\varDelta }_{ {f}_{ {\rm{a} }2} }/$% $ {K}_{{\rm{e}}{\rm{f}}{\rm{f}}2} $ $ {K}_{{\rm{e}}{\rm{f}}{\rm{f}}2}^{{\rm{*}}} $ 9 10 95014 93401 1.73 96069 94814 1.32 0.148 0.172 6 10 98432 96383 2.13 99471 97922 1.58 0.144 0.177 3 10 105005 102020 2.93 106094 103830 2.18 0.143 0.186 表 4 锥型变厚度环型压电换能器二阶共振频率和反共振频率
Table 4. The second-order resonant frequency and anti-resonant frequency of conical variable thickness annular piezoelectric transducer.
表3和表4中
$ f $ 和$ {f}^{*} $ 分别表示锥型变厚度环型压电超换能器径向振动共振频率的理论解和数值解,$ {K}_{{\rm{e}}{\rm{f}}{\rm{f}}} $ 和$ {K}_{{\rm{e}}{\rm{f}}{\rm{f}}}^{{\rm{*}}} $ 分别为其有效机电耦合系数的理论解和数值解, 相对误差$ \varDelta =\left|(f-{f}^{*})/{f}^{*}\right| $ . 由表3和表4可知, 径向一阶振动共振频率的理论解和数值解吻合非常好, 相对误差不超过0.2%, 径向二阶振动共振频率的理论解和数值解吻合较好, 相对误差不超过3%. 这验证了基于传输矩阵法的任意变厚度环型压电超声换能器径向振动共振频率理论解的正确性. 由于径向二阶的声波波长比一阶小, 耦合振动效应较为明显, 二阶共振频率、反共振频率的相对误差大于径向一阶.锥型变厚度环型压电超声换能器径向一阶、二阶共振频率和反共振频率随内外环厚度比
$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的变化曲线如图10所示. 由图10可知, 径向一阶共振频率、反共振频率随着随$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的增大而减小, 径向二阶共振频率、反共振频率随着随$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的增大而增大. 随着$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的增大, 一阶、二阶的共振频率和反共振频率的差值增大, 有效机电耦合系数增大.图 10 锥型变厚度环型压电超声换能器一阶、二阶径向共振频率和反共振频率与
$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的关系 (a)一阶共振和反共振; (b)二阶共振和反共振Figure 10. The relationship between the first and second order radial resonance frequency and the anti-resonance frequency and the thickness ratio
$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ of a conical variable thickness annular piezoelectric ultrasonic transducer: (a) The first-order radial resonance and anti-resonance; (b) the second-order radial resonance and anti-resonance. -
本文基于传输矩阵法将变厚度金属圆环的径向振动转化为
$ N $ 个等厚度金属圆环径向振动的叠加, 得到了任意变厚度金属薄圆环径向振动的等效电路图、共振频率方程和位移放大系数表达式, 分析了锥型、幂函数型、指数型、悬链线型金属圆环的位移放大系数随内外环厚度比$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的变化关系. 研究发现, 整体上一阶共振频率随$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的增大而减小, 二阶共振频率随着$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的增大而增大. 一阶、二阶的位移放大系数均随$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的增大而增大, 且径向一阶共振位移放大系数从大到小依次是锥型、幂函数型、指数型和悬链线型, 径向二阶共振位移放大系数从大到小依次是悬链线型、指数型、幂函数型、锥型. 在此基础上, 推导了由任意变厚度金属圆环和压电陶瓷圆环复合而成的压电超声换能器径向振动的等效电路和共振频率方程, 计算了不同尺寸的环型压电换能器的共振频率和反共振频率, 为了验证理论计算结果的正确性, 使用有限元软件进行仿真模拟, 模拟所得数值解与理论解符合较好, 最大相对误差不超过3%.
-
变厚度环型径向振动压电超声换能器可以实现阻抗变换、能量集中, 具有辐射面积大、全指向性等优点, 在功率超声、水声等领域被广泛应用. 由于求解复杂变厚度金属圆环径向振动的波动方程比较困难, 本文使用传输矩阵法将变厚度金属圆环的径向振动转化为
$ N $ 个等厚度金属圆环径向振动的叠加, 得到了任意变厚度金属薄圆环径向振动的等效电路图、共振频率方程和位移放大系数表达式, 分析了锥型、幂函数型、指数型、悬链线型金属圆环的位移放大系数与几何尺寸的关系. 在此基础上, 推导了由任意变厚度金属圆环和等厚度压电圆环复合而成的压电超声换能器径向振动的等效电路和共振频率方程. 为了验证理论结果的正确性, 使用有限元软件进行仿真, 所得一阶、二阶的共振频率和位移放大系数的数值解与理论解符合较好. 本研究给出了任意变厚度金属圆环径向振动的普适解, 为设计和优化径向压电超声换能器提供了理论指导.The variable thickness annular radial piezoelectric ultrasonic transducer can realize impedance transformation and energy concentration, has the advantages of large radiation area and full directivity, and is widely used in power ultrasound, underwater acoustic and other fields. Because solving complex variable thickness metal ring radial vibration wave equation is more difficult, in this paper, the radial vibration of metal rings with variable thickness is transformed into the superposition of the radial vibrations of N metal rings with equal thickness by using the transfer matrix method. The equivalent circuit diagram, the resonance frequency equation and the expression of the displacement amplification coefficient of the radial vibration of the metal thin ring with arbitrary thickness are obtained. The relationship between the displacement amplification coefficient and the geometric size of the cone, power function, exponential and catenary metal rings is analyzed. On this basis, the equivalent circuit and resonance frequency equation of radial vibration of piezoelectric ultrasonic transducer which is composed of a metal ring with variable thickness and a piezoelectric ring with equal thickness are derived. In order to verify the correctness of the theoretical results, the finite element software is used in simulation, and the numerical solutions of the first and second order resonance frequency and displacement amplification coefficients are in good agreement with the theoretical solutions. In this paper, the universal solution of radial vibration of metal ring with arbitrary variable thickness is given, which provides theoretical guidance for designing and optimizing the radial piezoelectric ultrasonic transducers.-
Keywords:
- annular piezoelectric ultrasonic transducer /
- radial vibration /
- transfer matrix method /
- equivalent circuit
[1] 孙淑珍, 李俊宝 2019 声学学报 44 743
Google Scholar
Sun S Z, Li J B 2019 Acta Acust. 44 743
Google Scholar
[2] 刘世清, 杨先莉, 张志良, 陈赵江 2013 声学学报 38 188
Google Scholar
Liu S Q, Yang X L, Zhang Z L, Chen Z J 2013 Acta Acust. 38 188
Google Scholar
[3] Jin O K, Jung G L 2007 J. Sound Vib. 300 241
Google Scholar
[4] Wang H M, Luo D S 2016 Appl. Math. Model. 40 2549
Google Scholar
[5] Hyun K L, Dabin J 2020 Optik (Stuttg) 222 165480
Google Scholar
[6] Luka P, Duje M, Fran M, Marko S, Marko B, Sven L 2022 IEEE Int. Ultrason. Symp. 124 106737
Google Scholar
[7] Kai X Z, Huang S M, Wu L, Ran T, Peng Y J, Mao Z M, Chen F, Li G R 2019 J. Mater. Sci. Technol. 35 2107
Google Scholar
[8] Feng J J, Liu L, Ma X Hong, Yuan J H, Wang A J 2017 Int. J. Hydrogen Energy 42 2071
Google Scholar
[9] Mohammad J B, Ammad M, Thomas L, Jochen B, Konrad K 2022 Bioresour Technol. 348 126785
Google Scholar
[10] Liu H B, Wang X K, Qin S, Lai W J, Yang X, Xu S Y 2021 Sci. Total Environ. 789 147862
Google Scholar
[11] Li S, Wang Y H, Wu S S, Niu W D, Yang S Q 2022 Appl. Math. Model. 109 455
Google Scholar
[12] Liu Z Y, Miao K, Tan Z M 2022 Appl. Acoust. 187 108497
Google Scholar
[13] 王晓宇, 林书玉 2021 声学学报 46 271
Google Scholar
Wang X Y, Lin S Y 2021 Acta Acust. 46 271
Google Scholar
[14] 刘垚, 李斌 2017 电加工与模具 03 61
Google Scholar
Liu Y, Li B 2017 Electromach. Mould. 03 61
Google Scholar
[15] 李井, 丁艳红, 梁欣 2016 机械设计与制造 7 197
Google Scholar
Li J, Ding Y H, Liang X 2016 Mach. Design Manufact. 7 197
Google Scholar
[16] 巩建辉, 王晨丰, 吴承启 2022 机械工程师 1 151
Google Scholar
Gong J H, Wang C F, Wu C Q 2022 Mech. Engineer. 1 151
Google Scholar
[17] Wang S, Shan J J, Lin S Y 2022 IEEE Int. Ultrason. Symp. 120 106640
Google Scholar
[18] 陈诚, 林书玉 2021 物理学报 70 017701
Google Scholar
Chen C, Lin S Y 2021 Acta Phys. Sin. 70 017701
Google Scholar
[19] Hu L Q, Wang S, Lin S Y 2022 Chin. Phys. B 31 508
Google Scholar
[20] 许龙, 李伟东 2019 声学学报 44 826
Google Scholar
Xu L, Li W D 2019 Acta Acust. 44 826
Google Scholar
[21] 许龙, 范秀梅 2021 应用声学 40 878
Google Scholar
Xu L, Fan X M 2021 J. Appl. Acoust. 40 878
Google Scholar
[22] 刘世清, 林书玉, 郭建中 2006 压电与声光 03 347
Google Scholar
Liu S Q, Lin S Y, Guo J Z 2006 Piezoelect. Acoustoop. 03 347
Google Scholar
[23] 王晓宇, 林书玉 2020 陕西师范大学学报(自然科学版) 48 107
Google Scholar
Wang X Y, Lin S Y 2020 J. Shaanxi Normal Univ. (Nat. Sci. Ed.) 48 107
Google Scholar
[24] 邓婷婷, 傅波 2017 机械工程师 1 36
Google Scholar
Deng T T, Fu B 2017 Mech. Engineer. 1 36
Google Scholar
[25] Lin S Y, Xu L 2012 IEEE Int. Ultrason. Symp. 52 103
Google Scholar
[26] Lin S Y 2007 Sens. Actuators A Phys. 134 505
Google Scholar
-
图 4 4种变厚度金属圆环一阶、二阶径向共振频率与
$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的关系 (a)一阶径向共振; (b)二阶径向共振Fig. 4. The relationship between the first and second order radial resonance frequencies of four kinds of variable thickness metal rings and thickness ratio
$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ : (a) First-order radial resonance; (b) second-order radial resonance.图 5 4种变厚度环型聚能器一阶和二阶径向共振位移放大系数随
$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的变化曲线 (a)一阶径向共振; (b)二阶径向共振Fig. 5. The relationship between the first and second order radial resonance displacement amplification coefficient and thickness ratio
$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ of four kinds of variable thickness metal rings: (a) First-order radial resonance; (b) second-order radial resonance.图 10 锥型变厚度环型压电超声换能器一阶、二阶径向共振频率和反共振频率与
$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ 的关系 (a)一阶共振和反共振; (b)二阶共振和反共振Fig. 10. The relationship between the first and second order radial resonance frequency and the anti-resonance frequency and the thickness ratio
$ {h}_{{\rm{a}}}/{h}_{{\rm{b}}} $ of a conical variable thickness annular piezoelectric ultrasonic transducer: (a) The first-order radial resonance and anti-resonance; (b) the second-order radial resonance and anti-resonance.表 1 变厚度金属圆环径向一阶、二阶共振频率
Table 1. Radial first and second order resonance frequencies of metal rings with variable thickness.
$ {h}_{{\rm{b}}}/{\rm{m}}{\rm{m}} $ $ {h}_{{\rm{a}}}/{\rm{m}}{\rm{m}} $ $ {f}_{{\rm{r}}1}/{\rm{H}}{\rm{z}} $ ${ {f}^{ {\rm{*} } }_{ {\rm{r} }1} }/{\rm{H} }{\rm{z} }$ ${\varDelta }_{ {f}_{ {\rm{r} }1} }/{\rm{\%} }$ $ {f}_{{\rm{r}}2}/{\rm{H}}{\rm{z}} $ ${ {f}^{*}_{r2} }/{\rm{H} }{\rm{z} }$ ${\varDelta }_{ {f}_{ {\rm{r} }2} }/{\rm{\%} }$ 锥型 6 10 21894 21892 0.01 112380 110790 1.44 幂函数型 5 10 21679 21679 0 113330 111720 1.44 指数型 4 10 21401 21399 0.01 112770 111480 1.16 悬链线型 3 10 20968 20957 0.05 109700 108890 0.74 表 2 变厚度金属圆环径向一阶、二阶共振位移放大系数
Table 2. Radial first and second order resonance displacement amplification coefficients of metal rings with variable thickness.
$ {h}_{{\rm{b}}}/{\rm{m}}{\rm{m}} $ $ {h}_{{\rm{a}}}/{\rm{m}}{\rm{m}} $ $ {M}_{{\rm{r}}1}^{{\rm{*}}} $ $ {M}_{{\rm{r}}1}^{{\rm{*}}{\rm{*}}} $ ${\varDelta }_{ {M}_{ {\rm{r} }1}^{*} }$/% $ {M}_{{\rm{r}}2}^{{\rm{*}}} $ $ {M}_{{\rm{r}}2}^{{\rm{*}}{\rm{*}}} $ ${\varDelta }_{ {M}_{ {\rm{r} }2}^{*} }/{\rm{\%} }$ 锥型 6 10 1.1985 1.1987 0.02 1.8304 1.8160 0.79 幂函数型 5 10 1.1997 1.2000 0.02 1.9876 1.9670 1.04 指数型 4 10 1.1992 1.1995 0.02 2.2501 2.2135 1.16 悬链线型 3 10 1.1959 1.1956 0.02 2.8118 2.7359 2.77 表 3 锥型变厚度环型压电换能器一阶共振频率和反共振频率
Table 3. The first-order resonant frequency and anti-resonant frequency of conical variable thickness annular piezoelectric transducer.
$ {h}_{{\rm{b}}}/{\rm{m}}{\rm{m}} $ $ {h}_{{\rm{a}}}/{\rm{m}}{\rm{m}} $ $ {f}_{{\rm{r}}1}/{\rm{H}}{\rm{z}} $ ${ {f}^{*}_{ {\rm{r} }1} }/{\rm{H} }{\rm{z} }$ ${\varDelta }_{ {f}_{ {\rm{r} }1} }/$% $ {f}_{{\rm{a}}1}/{\rm{H}}{\rm{z}} $ ${ {f}^{*}_{ {\rm{a} }1} }/{\rm{H} }{\rm{z} }$ ${\varDelta }_{ {f}_{ {\rm{a} }1} }/$% $ {K}_{{\rm{e}}{\rm{f}}{\rm{f}}1} $ $ {K}_{{\rm{e}}{\rm{f}}{\rm{f}}1}^{{\rm{*}}} $ 9 10 22002 21989 0.06 22276 22269 0.03 0.156 0.158 6 10 21059 21042 0.08 21355 21348 0.03 0.166 0.169 3 10 19886 19856 0.15 20209 20192 0.08 0.178 0.182 表 4 锥型变厚度环型压电换能器二阶共振频率和反共振频率
Table 4. The second-order resonant frequency and anti-resonant frequency of conical variable thickness annular piezoelectric transducer.
$ {h}_{{\rm{b}}}/{\rm{m}}{\rm{m}} $ $ {h}_{{\rm{a}}}/{\rm{m}}{\rm{m}} $ $ {f}_{{\rm{r}}2}/{\rm{H}}{\rm{z}} $ ${ {f}^{*}_{ {\rm{r} }2} }/{\rm{H} }{\rm{z} }$ ${\varDelta }_{ {f}_{ {\rm{r} }2} }/$% $ {f}_{{\rm{a}}2}/{\rm{H}}{\rm{z}} $ ${ {f}^{*}_{ {\rm{a} }2} }/{\rm{H} }{\rm{z} }$ ${\varDelta }_{ {f}_{ {\rm{a} }2} }/$% $ {K}_{{\rm{e}}{\rm{f}}{\rm{f}}2} $ $ {K}_{{\rm{e}}{\rm{f}}{\rm{f}}2}^{{\rm{*}}} $ 9 10 95014 93401 1.73 96069 94814 1.32 0.148 0.172 6 10 98432 96383 2.13 99471 97922 1.58 0.144 0.177 3 10 105005 102020 2.93 106094 103830 2.18 0.143 0.186 -
[1] 孙淑珍, 李俊宝 2019 声学学报 44 743
Google Scholar
Sun S Z, Li J B 2019 Acta Acust. 44 743
Google Scholar
[2] 刘世清, 杨先莉, 张志良, 陈赵江 2013 声学学报 38 188
Google Scholar
Liu S Q, Yang X L, Zhang Z L, Chen Z J 2013 Acta Acust. 38 188
Google Scholar
[3] Jin O K, Jung G L 2007 J. Sound Vib. 300 241
Google Scholar
[4] Wang H M, Luo D S 2016 Appl. Math. Model. 40 2549
Google Scholar
[5] Hyun K L, Dabin J 2020 Optik (Stuttg) 222 165480
Google Scholar
[6] Luka P, Duje M, Fran M, Marko S, Marko B, Sven L 2022 IEEE Int. Ultrason. Symp. 124 106737
Google Scholar
[7] Kai X Z, Huang S M, Wu L, Ran T, Peng Y J, Mao Z M, Chen F, Li G R 2019 J. Mater. Sci. Technol. 35 2107
Google Scholar
[8] Feng J J, Liu L, Ma X Hong, Yuan J H, Wang A J 2017 Int. J. Hydrogen Energy 42 2071
Google Scholar
[9] Mohammad J B, Ammad M, Thomas L, Jochen B, Konrad K 2022 Bioresour Technol. 348 126785
Google Scholar
[10] Liu H B, Wang X K, Qin S, Lai W J, Yang X, Xu S Y 2021 Sci. Total Environ. 789 147862
Google Scholar
[11] Li S, Wang Y H, Wu S S, Niu W D, Yang S Q 2022 Appl. Math. Model. 109 455
Google Scholar
[12] Liu Z Y, Miao K, Tan Z M 2022 Appl. Acoust. 187 108497
Google Scholar
[13] 王晓宇, 林书玉 2021 声学学报 46 271
Google Scholar
Wang X Y, Lin S Y 2021 Acta Acust. 46 271
Google Scholar
[14] 刘垚, 李斌 2017 电加工与模具 03 61
Google Scholar
Liu Y, Li B 2017 Electromach. Mould. 03 61
Google Scholar
[15] 李井, 丁艳红, 梁欣 2016 机械设计与制造 7 197
Google Scholar
Li J, Ding Y H, Liang X 2016 Mach. Design Manufact. 7 197
Google Scholar
[16] 巩建辉, 王晨丰, 吴承启 2022 机械工程师 1 151
Google Scholar
Gong J H, Wang C F, Wu C Q 2022 Mech. Engineer. 1 151
Google Scholar
[17] Wang S, Shan J J, Lin S Y 2022 IEEE Int. Ultrason. Symp. 120 106640
Google Scholar
[18] 陈诚, 林书玉 2021 物理学报 70 017701
Google Scholar
Chen C, Lin S Y 2021 Acta Phys. Sin. 70 017701
Google Scholar
[19] Hu L Q, Wang S, Lin S Y 2022 Chin. Phys. B 31 508
Google Scholar
[20] 许龙, 李伟东 2019 声学学报 44 826
Google Scholar
Xu L, Li W D 2019 Acta Acust. 44 826
Google Scholar
[21] 许龙, 范秀梅 2021 应用声学 40 878
Google Scholar
Xu L, Fan X M 2021 J. Appl. Acoust. 40 878
Google Scholar
[22] 刘世清, 林书玉, 郭建中 2006 压电与声光 03 347
Google Scholar
Liu S Q, Lin S Y, Guo J Z 2006 Piezoelect. Acoustoop. 03 347
Google Scholar
[23] 王晓宇, 林书玉 2020 陕西师范大学学报(自然科学版) 48 107
Google Scholar
Wang X Y, Lin S Y 2020 J. Shaanxi Normal Univ. (Nat. Sci. Ed.) 48 107
Google Scholar
[24] 邓婷婷, 傅波 2017 机械工程师 1 36
Google Scholar
Deng T T, Fu B 2017 Mech. Engineer. 1 36
Google Scholar
[25] Lin S Y, Xu L 2012 IEEE Int. Ultrason. Symp. 52 103
Google Scholar
[26] Lin S Y 2007 Sens. Actuators A Phys. 134 505
Google Scholar
计量
- 文章访问数: 374
- PDF下载量: 22
- 被引次数: 0