搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于传输矩阵法的任意变厚度环型压电超声换能器

董宜雷 陈诚 林书玉

引用本文:
Citation:

基于传输矩阵法的任意变厚度环型压电超声换能器

董宜雷, 陈诚, 林书玉

Arbitrary variable thickness annular piezoelectric ultrasonic transducer based on transfer matrix method

Dong Yi-Lei, Chen Cheng, Lin Shu-Yu
PDF
HTML
导出引用
  • 变厚度环型径向振动压电超声换能器可以实现阻抗变换、能量集中, 具有辐射面积大、全指向性等优点, 在功率超声、水声等领域被广泛应用. 由于求解复杂变厚度金属圆环径向振动的波动方程比较困难, 本文使用传输矩阵法将变厚度金属圆环的径向振动转化为$ 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.
      通信作者: 林书玉, sylin@snnu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11674206, 11874253, 12174240)资助的课题.
      Corresponding author: Lin Shu-Yu, sylin@snnu.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 11674206, 11874253, 12174240).
    [1]

    孙淑珍, 李俊宝 2019 声学学报 44 743Google Scholar

    Sun S Z, Li J B 2019 Acta Acust. 44 743Google Scholar

    [2]

    刘世清, 杨先莉, 张志良, 陈赵江 2013 声学学报 38 188Google Scholar

    Liu S Q, Yang X L, Zhang Z L, Chen Z J 2013 Acta Acust. 38 188Google Scholar

    [3]

    Jin O K, Jung G L 2007 J. Sound Vib. 300 241Google Scholar

    [4]

    Wang H M, Luo D S 2016 Appl. Math. Model. 40 2549Google Scholar

    [5]

    Hyun K L, Dabin J 2020 Optik (Stuttg) 222 165480Google Scholar

    [6]

    Luka P, Duje M, Fran M, Marko S, Marko B, Sven L 2022 IEEE Int. Ultrason. Symp. 124 106737Google 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 2107Google Scholar

    [8]

    Feng J J, Liu L, Ma X Hong, Yuan J H, Wang A J 2017 Int. J. Hydrogen Energy 42 2071Google Scholar

    [9]

    Mohammad J B, Ammad M, Thomas L, Jochen B, Konrad K 2022 Bioresour Technol. 348 126785Google Scholar

    [10]

    Liu H B, Wang X K, Qin S, Lai W J, Yang X, Xu S Y 2021 Sci. Total Environ. 789 147862Google Scholar

    [11]

    Li S, Wang Y H, Wu S S, Niu W D, Yang S Q 2022 Appl. Math. Model. 109 455Google Scholar

    [12]

    Liu Z Y, Miao K, Tan Z M 2022 Appl. Acoust. 187 108497Google Scholar

    [13]

    王晓宇, 林书玉 2021 声学学报 46 271Google Scholar

    Wang X Y, Lin S Y 2021 Acta Acust. 46 271Google Scholar

    [14]

    刘垚, 李斌 2017 电加工与模具 03 61Google Scholar

    Liu Y, Li B 2017 Electromach. Mould. 03 61Google Scholar

    [15]

    李井, 丁艳红, 梁欣 2016 机械设计与制造 7 197Google Scholar

    Li J, Ding Y H, Liang X 2016 Mach. Design Manufact. 7 197Google Scholar

    [16]

    巩建辉, 王晨丰, 吴承启 2022 机械工程师 1 151Google Scholar

    Gong J H, Wang C F, Wu C Q 2022 Mech. Engineer. 1 151Google Scholar

    [17]

    Wang S, Shan J J, Lin S Y 2022 IEEE Int. Ultrason. Symp. 120 106640Google Scholar

    [18]

    陈诚, 林书玉 2021 物理学报 70 017701Google Scholar

    Chen C, Lin S Y 2021 Acta Phys. Sin. 70 017701Google Scholar

    [19]

    Hu L Q, Wang S, Lin S Y 2022 Chin. Phys. B 31 508Google Scholar

    [20]

    许龙, 李伟东 2019 声学学报 44 826Google Scholar

    Xu L, Li W D 2019 Acta Acust. 44 826Google Scholar

    [21]

    许龙, 范秀梅 2021 应用声学 40 878Google Scholar

    Xu L, Fan X M 2021 J. Appl. Acoust. 40 878Google Scholar

    [22]

    刘世清, 林书玉, 郭建中 2006 压电与声光 03 347Google Scholar

    Liu S Q, Lin S Y, Guo J Z 2006 Piezoelect. Acoustoop. 03 347Google Scholar

    [23]

    王晓宇, 林书玉 2020 陕西师范大学学报(自然科学版) 48 107Google Scholar

    Wang X Y, Lin S Y 2020 J. Shaanxi Normal Univ. (Nat. Sci. Ed.) 48 107Google Scholar

    [24]

    邓婷婷, 傅波 2017 机械工程师 1 36Google Scholar

    Deng T T, Fu B 2017 Mech. Engineer. 1 36Google Scholar

    [25]

    Lin S Y, Xu L 2012 IEEE Int. Ultrason. Symp. 52 103Google Scholar

    [26]

    Lin S Y 2007 Sens. Actuators A Phys. 134 505Google Scholar

  • 图 1  变厚度环型压电超声换能器的纵向截面示意图

    Fig. 1.  Schematic diagram of longitudinal section of variable thickness annular piezoelectric ultrasonic transducer.

    图 2  等厚度金属薄圆环径向振动的等效电路图

    Fig. 2.  Equivalent circuit of metal thin circular annular in radial vibration.

    图 3  变厚度金属薄圆环径向振动的等效电路图

    Fig. 3.  Equivalent circuit diagram of metal thin annular with variable thickness in radial vibration.

    图 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.

    图 6  锥型金属圆环的一阶、二阶径向共振频率与N的关系 (a)一阶径向共振; (b)二阶径向共振

    Fig. 6.  The relationship between the first and second order radial resonant frequencies of conical metal rings and N: (a) First-order radial resonance; (b) second-order radial resonance.

    图 7  压电陶瓷圆环径向振动的等效电路图

    Fig. 7.  Equivalent circuit diagram of piezoelectric ceramic annular in radial vibration.

    图 8  环型压电超声换能器径向振动的等效电路图

    Fig. 8.  Equivalent circuit diagram of circular piezoelectric ultrasonic transducer in radial vibration.

    图 9  锥型变厚度环型压电超声换能器一阶、二阶径向共振振型图 (a)一阶径向共振; (b)二阶径向共振

    Fig. 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.

    图 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{\%} }$
    锥型61021894218920.011123801107901.44
    幂函数型510216792167901133301117201.44
    指数型41021401213990.011127701114801.16
    悬链线型31020968209570.051097001088900.74
    下载: 导出CSV

    表 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{\%} }$
    锥型6101.19851.19870.021.83041.81600.79
    幂函数型5101.19971.20000.021.98761.96701.04
    指数型4101.19921.19950.022.25012.21351.16
    悬链线型3101.19591.19560.022.81182.73592.77
    下载: 导出CSV

    表 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{*}}} $
    91022002219890.0622276222690.030.1560.158
    61021059210420.0821355213480.030.1660.169
    31019886198560.1520209201920.080.1780.182
    下载: 导出CSV

    表 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{*}}} $
    91095014934011.7396069948141.320.1480.172
    61098432963832.1399471979221.580.1440.177
    3101050051020202.931060941038302.180.1430.186
    下载: 导出CSV
  • [1]

    孙淑珍, 李俊宝 2019 声学学报 44 743Google Scholar

    Sun S Z, Li J B 2019 Acta Acust. 44 743Google Scholar

    [2]

    刘世清, 杨先莉, 张志良, 陈赵江 2013 声学学报 38 188Google Scholar

    Liu S Q, Yang X L, Zhang Z L, Chen Z J 2013 Acta Acust. 38 188Google Scholar

    [3]

    Jin O K, Jung G L 2007 J. Sound Vib. 300 241Google Scholar

    [4]

    Wang H M, Luo D S 2016 Appl. Math. Model. 40 2549Google Scholar

    [5]

    Hyun K L, Dabin J 2020 Optik (Stuttg) 222 165480Google Scholar

    [6]

    Luka P, Duje M, Fran M, Marko S, Marko B, Sven L 2022 IEEE Int. Ultrason. Symp. 124 106737Google 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 2107Google Scholar

    [8]

    Feng J J, Liu L, Ma X Hong, Yuan J H, Wang A J 2017 Int. J. Hydrogen Energy 42 2071Google Scholar

    [9]

    Mohammad J B, Ammad M, Thomas L, Jochen B, Konrad K 2022 Bioresour Technol. 348 126785Google Scholar

    [10]

    Liu H B, Wang X K, Qin S, Lai W J, Yang X, Xu S Y 2021 Sci. Total Environ. 789 147862Google Scholar

    [11]

    Li S, Wang Y H, Wu S S, Niu W D, Yang S Q 2022 Appl. Math. Model. 109 455Google Scholar

    [12]

    Liu Z Y, Miao K, Tan Z M 2022 Appl. Acoust. 187 108497Google Scholar

    [13]

    王晓宇, 林书玉 2021 声学学报 46 271Google Scholar

    Wang X Y, Lin S Y 2021 Acta Acust. 46 271Google Scholar

    [14]

    刘垚, 李斌 2017 电加工与模具 03 61Google Scholar

    Liu Y, Li B 2017 Electromach. Mould. 03 61Google Scholar

    [15]

    李井, 丁艳红, 梁欣 2016 机械设计与制造 7 197Google Scholar

    Li J, Ding Y H, Liang X 2016 Mach. Design Manufact. 7 197Google Scholar

    [16]

    巩建辉, 王晨丰, 吴承启 2022 机械工程师 1 151Google Scholar

    Gong J H, Wang C F, Wu C Q 2022 Mech. Engineer. 1 151Google Scholar

    [17]

    Wang S, Shan J J, Lin S Y 2022 IEEE Int. Ultrason. Symp. 120 106640Google Scholar

    [18]

    陈诚, 林书玉 2021 物理学报 70 017701Google Scholar

    Chen C, Lin S Y 2021 Acta Phys. Sin. 70 017701Google Scholar

    [19]

    Hu L Q, Wang S, Lin S Y 2022 Chin. Phys. B 31 508Google Scholar

    [20]

    许龙, 李伟东 2019 声学学报 44 826Google Scholar

    Xu L, Li W D 2019 Acta Acust. 44 826Google Scholar

    [21]

    许龙, 范秀梅 2021 应用声学 40 878Google Scholar

    Xu L, Fan X M 2021 J. Appl. Acoust. 40 878Google Scholar

    [22]

    刘世清, 林书玉, 郭建中 2006 压电与声光 03 347Google Scholar

    Liu S Q, Lin S Y, Guo J Z 2006 Piezoelect. Acoustoop. 03 347Google Scholar

    [23]

    王晓宇, 林书玉 2020 陕西师范大学学报(自然科学版) 48 107Google Scholar

    Wang X Y, Lin S Y 2020 J. Shaanxi Normal Univ. (Nat. Sci. Ed.) 48 107Google Scholar

    [24]

    邓婷婷, 傅波 2017 机械工程师 1 36Google Scholar

    Deng T T, Fu B 2017 Mech. Engineer. 1 36Google Scholar

    [25]

    Lin S Y, Xu L 2012 IEEE Int. Ultrason. Symp. 52 103Google Scholar

    [26]

    Lin S Y 2007 Sens. Actuators A Phys. 134 505Google Scholar

  • [1] 张陶然, 莫润阳, 胡静, 陈时, 王成会, 郭建中. 黏弹介质包裹的液体腔中气泡的动力学分析. 物理学报, 2021, 70(12): 124301. doi: 10.7498/aps.70.20201876
    [2] 陈诚, 林书玉. 基于2-2型压电复合材料的新型宽频带径向振动超声换能器. 物理学报, 2021, 70(1): 017701. doi: 10.7498/aps.70.20201352
    [3] 李宇涵, 邓联文, 罗衡, 贺龙辉, 贺君, 徐运超, 黄生祥. 双层螺旋环超表面复合吸波体等效电路模型及微波损耗机制. 物理学报, 2019, 68(9): 095201. doi: 10.7498/aps.68.20181960
    [4] 易红霞, 肖刘, 苏小保. 传输矩阵法在行波管内部反射引起的增益波动计算中的应用. 物理学报, 2016, 65(12): 128401. doi: 10.7498/aps.65.128401
    [5] 陈姝媛, 阮存军, 王勇. 带状注速调管多间隙扩展互作用输出腔等效电路的研究. 物理学报, 2014, 63(2): 028402. doi: 10.7498/aps.63.028402
    [6] 吴超, 吕绪良, 曾朝阳, 贾其. 基于阻抗模拟的等效电磁参数研究. 物理学报, 2013, 62(5): 054101. doi: 10.7498/aps.62.054101
    [7] 王秀芝, 高劲松, 徐念喜. 利用等效电路模型快速分析加载集总元件的微型化频率选择表面. 物理学报, 2013, 62(20): 207301. doi: 10.7498/aps.62.207301
    [8] 尹彬, 柏云龙, 齐艳辉, 冯素春, 简水生. 拉锥型啁啾光纤光栅滤波器的研究. 物理学报, 2013, 62(21): 214213. doi: 10.7498/aps.62.214213
    [9] 胡丰伟, 包伯成, 武花干, 王春丽. 荷控忆阻器等效电路分析模型及其电路特性研究. 物理学报, 2013, 62(21): 218401. doi: 10.7498/aps.62.218401
    [10] 张小丽, 林书玉, 付志强, 王勇. 弯曲振动薄圆盘的共振频率和等效电路参数研究. 物理学报, 2013, 62(3): 034301. doi: 10.7498/aps.62.034301
    [11] 白春江, 李建清, 胡玉禄, 杨中海, 李斌. 利用等效电路模型计算耦合腔行波管注-波互作用. 物理学报, 2012, 61(17): 178401. doi: 10.7498/aps.61.178401
    [12] 包伯成, 胡文, 许建平, 刘中, 邹凌. 忆阻混沌电路的分析与实现. 物理学报, 2011, 60(12): 120502. doi: 10.7498/aps.60.120502
    [13] 罗小光, 何济洲. 摇摆型棘齿热电子隧穿制冷. 物理学报, 2011, 60(9): 090506. doi: 10.7498/aps.60.090506
    [14] 於陆勒, 盛政明, 张 杰. 均匀等离子体光栅的色散特性研究. 物理学报, 2008, 57(10): 6457-6464. doi: 10.7498/aps.57.6457
    [15] 童 凯, 崔卫卫, 汪梅婷, 李志全. 一维缺陷光子晶体温度的测量. 物理学报, 2008, 57(2): 762-766. doi: 10.7498/aps.57.762
    [16] 田 赫, 掌蕴东, 王 号, 邱 巍, 王 楠, 袁 萍. 微环耦合谐振光波导中的色散控制模型与数值仿真. 物理学报, 2008, 57(10): 6400-6403. doi: 10.7498/aps.57.6400
    [17] 田 赫, 掌蕴东, 王 号, 邱 巍, 王 楠, 袁 萍. 光脉冲在微环耦合谐振光波导中传输线性特性的数值仿真. 物理学报, 2008, 57(11): 7012-7016. doi: 10.7498/aps.57.7012
    [18] 郝保良, 刘濮鲲, 唐昌建. 二维非正交坐标斜方格金属光子带隙结构. 物理学报, 2006, 55(4): 1862-1867. doi: 10.7498/aps.55.1862
    [19] 李洪奇. 介观压电石英晶体等效电路的量子化. 物理学报, 2005, 54(3): 1361-1365. doi: 10.7498/aps.54.1361
    [20] 王均宏. 脉冲电压电流沿偶极天线传播过程的等效电路法分析. 物理学报, 2000, 49(9): 1696-1701. doi: 10.7498/aps.49.1696
计量
  • 文章访问数:  375
  • PDF下载量:  22
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-11-03
  • 修回日期:  2022-11-26
  • 上网日期:  2022-12-17
  • 刊出日期:  2023-03-05

基于传输矩阵法的任意变厚度环型压电超声换能器

  • 陕西师范大学物理学与信息技术学院, 陕西省超声重点实验室, 西安 710119
  • 通信作者: 林书玉, sylin@snnu.edu.cn
    基金项目: 国家自然科学基金(批准号: 11674206, 11874253, 12174240)资助的课题.

摘要: 变厚度环型径向振动压电超声换能器可以实现阻抗变换、能量集中, 具有辐射面积大、全指向性等优点, 在功率超声、水声等领域被广泛应用. 由于求解复杂变厚度金属圆环径向振动的波动方程比较困难, 本文使用传输矩阵法将变厚度金属圆环的径向振动转化为$ N $个等厚度金属圆环径向振动的叠加, 得到了任意变厚度金属薄圆环径向振动的等效电路图、共振频率方程和位移放大系数表达式, 分析了锥型、幂函数型、指数型、悬链线型金属圆环的位移放大系数与几何尺寸的关系. 在此基础上, 推导了由任意变厚度金属圆环和等厚度压电圆环复合而成的压电超声换能器径向振动的等效电路和共振频率方程. 为了验证理论结果的正确性, 使用有限元软件进行仿真, 所得一阶、二阶的共振频率和位移放大系数的数值解与理论解符合较好. 本研究给出了任意变厚度金属圆环径向振动的普适解, 为设计和优化径向压电超声换能器提供了理论指导.

English Abstract

    • 径向振动压电超声换能器具有辐射面积大、沿半径方向有全指向性、接收灵敏度高、结构简单且性能稳定等优点[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}}} $. 假设变厚度环型压电超声换能器的轴向厚度尺寸远小于径向尺寸, 忽略其纵向振动所带来的耦合效应, 此时环型压电超声换能器做纯径向振动. 环型压电超声换能器的径向振动由金属内圆环径向振动和压电陶瓷外圆环径向振动复合而成, 因此需要分别研究变厚度金属圆环的径向振动和等厚度压电陶瓷圆环的径向振动.

      图  1  变厚度环型压电超声换能器的纵向截面示意图

      Figure 1.  Schematic diagram of longitudinal section of variable thickness annular piezoelectric ultrasonic transducer.

    • 根据轴向厚度沿半径的变化函数, 变厚度金属薄圆环可以分成锥型、阶梯型、幂函数型等. 为了求解任意变厚度金属薄圆环径向振动的等效电路, 本文采用传输矩阵法, 将变厚度金属薄圆环沿半径方向等分成$ 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{\%} }$
      锥型61021894218920.011123801107901.44
      幂函数型510216792167901133301117201.44
      指数型41021401213990.011127701114801.16
      悬链线型31020968209570.051097001088900.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{\%} }$
      锥型6101.19851.19870.021.83041.81600.79
      幂函数型5101.19971.20000.021.98761.96701.04
      指数型4101.19921.19950.022.25012.21351.16
      悬链线型3101.19591.19560.022.81182.73592.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}} $时, 共振频率的计算精度满足要求.

      图  6  锥型金属圆环的一阶、二阶径向共振频率与N的关系 (a)一阶径向共振; (b)二阶径向共振

      Figure 6.  The relationship between the first and second order radial resonant frequencies of conical metal rings and N: (a) First-order radial resonance; (b) second-order radial resonance.

    • 根据文献[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{*}}} $
      91022002219890.0622276222690.030.1560.158
      61021059210420.0821355213480.030.1660.169
      31019886198560.1520209201920.080.1780.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{*}}} $
      91095014934011.7396069948141.320.1480.172
      61098432963832.1399471979221.580.1440.177
      3101050051020202.931060941038302.180.1430.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%.

参考文献 (26)

目录

    /

    返回文章
    返回