基于声黑洞-Moonie换能器结构的微型化超声手术刀

,  ,  ,  , 

上一篇

下一篇

陈诚, 陈慧琴, 徐春龙, 郭建中, 林书玉. 基于声黑洞-Moonie换能器结构的微型化超声手术刀[J]. 物理学报, 2026, 75(8): 081003. doi: 10.7498/aps.75.20251458
引用本文: 陈诚, 陈慧琴, 徐春龙, 郭建中, 林书玉. 基于声黑洞-Moonie换能器结构的微型化超声手术刀[J]. 物理学报, 2026, 75(8): 081003. doi: 10.7498/aps.75.20251458
Cheng CHEN, Huiqin CHEN, Chunlong XU, Jianzhong GUO, Shuyu LIN. A miniaturized ultrasonic scalpel based on acoustic black hole and Moonie transducer structures[J]. Acta Physica Sinica, 2026, 75(8): 081003. doi: 10.7498/aps.75.20251458
Citation: Cheng CHEN, Huiqin CHEN, Chunlong XU, Jianzhong GUO, Shuyu LIN. A miniaturized ultrasonic scalpel based on acoustic black hole and Moonie transducer structures[J]. Acta Physica Sinica, 2026, 75(8): 081003. doi: 10.7498/aps.75.20251458

基于声黑洞-Moonie换能器结构的微型化超声手术刀

    通讯作者: E-mail: guojz@snnu.edu.cn.;  E-mail: sylin@snnu.edu.cn.

A miniaturized ultrasonic scalpel based on acoustic black hole and Moonie transducer structures

    Corresponding authors: E-mail: guojz@snnu.edu.cn.;  E-mail: sylin@snnu.edu.cn.
  • 摘要: 本文提出了一种基于声黑洞(acoustic black hole, ABH)与Moonie换能器结构的微型化超声手术刀, 旨在利用ABH结构的能量局域聚集效应以及Moonie换能器的位移放大特性有效解决超声手术刀在微型化过程中振动性能衰减问题. 该手术刀由声黑洞轮廓金属端帽、工作刀头与纵向极化压电圆环组成, 并实现了多重振动模式转换. 利用等效电路法建立了ABH-Moonie型手术刀多模态耦合振动的解析理论模型, 实现了快速频率预测. 采用有限元法分析了ABH-Moonie型手术刀的耦合振动特性, 并探究了声黑洞参数和负载对振动性能的影响. 相较于普通Moonie结构手术刀(m < 2), ABH-Moonie型手术刀最大振幅输出显著提升. 在负载条件下, 手术刀仍能保持良好的振动性能, 包括较小的频率漂移以及振动减幅. 最后加工了手术刀样机并进行实验测试, 实测手术刀的设计模态能够被有效激发, 频率误差仅为1.5%, 充分验证了ABH-Moonie复合增效结构的可行性以及解析理论模型的准确性. 本研究可以为微型化高性能压电器件的设计优化提供全新视角.
  • 加载中
  • 图 1  ABH-Moonie型手术刀的结构示意图 (a)组成部件; (b)声黑洞轮廓设计; (c)几何尺寸

    Figure 1.  Structure diagram of the ABH-Moonie type scalpel: (a) Components; (b) ABH-profile design; (c) geometric dimensions.

    图 2  ABH-Moonie型手术刀耦合振动的机电等效电路

    Figure 2.  Electromechanical equivalent circuit of the ABH-Moonie type scalpel in coupled vibration.

    图 3  电导纳-频率响应曲线

    Figure 3.  Electrical admittance-frequency response curve.

    图 4  手术刀的耦合振动模态 (a) 3D剖面图; (b) 2D剖面图

    Figure 4.  Coupled vibration mode of the ABH-Moonie type scalpel: (a) 3D cross-sectional view; (b) 2D cross-sectional view.

    图 5  共振/反共振频率、有效机电耦合系数随声黑洞参数m变化曲线

    Figure 5.  Curves of resonance/antiresonance frequencies and effective electromechanical coupling coefficients varying with the ABH parameter m.

    图 6  归一化最大振动位移随声黑洞参数m变化曲线

    Figure 6.  Curve of normalized maximum vibration displacement varying with the ABH parameter m.

    图 7  ABH-Moonie型手术刀切割仿体组织的有限元建模

    Figure 7.  Finite element model of the ABH-Moonie type scalpel for cutting phantom tissue.

    图 8  仿体负载对电导纳频响特性的影响曲线

    Figure 8.  Effect curves of phantom load on the frequency response of electrical admittance.

    图 9  仿体负载对刀头振动位移的影响曲线

    Figure 9.  Effect curves of phantom load on the vibration displacement of the blade.

    图 10  ABH-Moonie型手术刀样机

    Figure 10.  ABH-Moonie type scalpel prototype.

    图 11  ABH-Moonie型手术刀的阻抗分析实验

    Figure 11.  Impedance analysis experiment of the ABH-Moonie type scalpel.

    图 12  测量得到的ABH-Moonie型手术刀的电阻抗频响及振动相位曲线

    Figure 12.  Measured electrical impedance frequency response and vibration phase curves of the ABH-Moonie type scalpel.

    图 13  ABH-Moonie型手术刀的模态测量实验

    Figure 13.  Mode measurement experiment of the ABH-Moonie type scalpel.

    图 14  测量得到的ABH-Moonie型手术刀的耦合振动模态

    Figure 14.  Measured coupled vibration mode of the ABH-Moonie type scalpel.

    表 1  ABH-Moonie型手术刀的材料参数

    Table 1.  Material parameters of the ABH-Moonie type scalpel.

    参数数值参数数值
    $ \rho /(\text{kg}{\cdot} {\text{m}}^{-}{}^{3}) $2700$ c_{13}^{\text{E}}/({\text{m}}^{2}{\cdot} {\text{N}}^{\text{-1}}) $7.43×1010
    $ E\text{/}(\text{N}{\cdot} {\text{m}}^{-}{}^{2}) $7×1010$ c_{33}^{\text{E}}/({\text{m}}^{2}{\cdot} {\text{N}}^{\text{-1}}) $11.5×1010
    $ \sigma $0.33$ {e}_{31}\text{/}(\text{N}{\cdot} \text{m}{\cdot} {\text{V}}^{\text{-1}}) $-5.2
    $ {\rho }_{0}/(\text{kg}{\cdot} {\text{m}}^{-}{}^{3}) $7500$ {e}_{33}\text{/}(\text{N}{\cdot} \text{m}{\cdot} {\text{V}}^{\text{-1}}) $15.1
    $ c_{11}^{\text{E}}/({\text{m}}^{2}{\cdot} {\text{N}}^{\text{-1}}) $13.9×1010$ \varepsilon _{33}^{\text{S}}\text{/}(\text{C}{\cdot} {\text{m}}^{-1}) $5.69×10–9
    $ c_{12}^{\text{E}}/({\text{m}}^{2}{\cdot} {\text{N}}^{\text{-1}}) $7.78×1010
    下载: 导出CSV
  • [1] Sun D, Zhou Z Y, Liu Y H, Shen W Z 1997 IEEE Trans. Biomed. Eng. 44 462 doi: 10.1109/10.581934
    [2] Schafer M E 2023 IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 70 147 doi: 10.1109/TUFFC.2022.3225519
    [3] O'Daly B J, Morris E, Gavin G P, O'Byrne J M, McGuinness G B 2008 J. Mater. Process. Technol. 200 38 doi: 10.1016/j.jmatprotec.2007.11.041
    [4] Wang S, Shan J J, Lin S Y 2024 Chin. Phys. B. 33 104302 doi: 10.1088/1674-1056/ad6a0c
    [5] Chen Q, Yao M, Li L, Liu X N, Liu Y Q, Wang F F, Tang Y X, Zhang Q Z, Lin D, Luo H S, Zhao X Y 2025 Sens. Actuator Phys. 382 116149 doi: 10.1016/j.sna.2024.116149
    [6] Li X, Stritch T, Manley K, Lucas M 2021 IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 68 2543 doi: 10.1109/TUFFC.2021.3065207
    [7] 胡负稷, 魏翔宇, 王海波, 周佳骏, 冯耿超 2023 应用声学 5 1053 Hu F J, Wei X Y, Wang H B, Zhou J J, Feng G C 2023 Deng J, Gao N S, Chen X 2023 J. Appl. Acoust. 5 1053
    [8] Hofmann F, Haeberlin A, Brot S D, Stahel A, Keppner H, Burger J 2023 Ultrasonics. 130 106927 doi: 10.1016/j.ultras.2023.106927
    [9] Kuang Y, Sadiq M, Cochran S, Cochran S, Huang Z H 2015 IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 62 1957 doi: 10.1109/TUFFC.2015.007300
    [10] Bejarano F, Feeney A, Wallace R, Simpson H, Lucas M 2016 Ultrasonics. 72 24 doi: 10.1016/j.ultras.2016.07.004
    [11] Dogan A, Uchino K, Newnham R E 1997 IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 44 597 doi: 10.1109/58.658312
    [12] Narayanan, M, Schwartz, R W 2010 J Electroceram. 24 205 doi: 10.1007/s10832-008-9559-2
    [13] Ochoa P, Pons J L, Villegas M, Fernandez J F 2006 Sens. Actuator Phys. 132 63 doi: 10.1016/j.sna.2006.05.031
    [14] Gao N S, Guo X Y, Deng J, Cheng B Z, Hou H 2021 Appl. Acoust. 173 107694 doi: 10.1016/j.apacoust.2020.107694
    [15] Lee J Y, Jeon W 2017 J. Acoust. Soc. Am. 141 1437 doi: 10.1121/1.4976687
    [16] Mi Y Z, Zhai W, Cheng L, Xi C Y, Yu X 2021 Appl. Phys. Lett. 118 114101 doi: 10.1063/5.0042514
    [17] 朱浩 李俊宝 葛晓辉 李德鹏 2025 物理学报 74 194301 doi: 10.7498/aps.74.20250741 Zhu H, Li J B, Ge X H, Li D P 2025 Acta Phys. Sin. 74 194301 doi: 10.7498/aps.74.20250741
    [18] Zhang L F, Tang X, Qin Z Y, Chu F L 2022 Appl. Phys. Lett. 121 013902 doi: 10.1063/5.0089382
    [19] Liu P Z, Tian Z H, Yang K C, Naquin T D, Hao N J, Huang H Y, Chen J Y, Ma Q X, Bachman H, Zhang P, Xu X H, Hu J H, Huang T J 2022 Sci. Adv. 8 eabm2592 doi: 10.1126/sciadv.abm2592
    [20] Liu P Z, Huang H Y, Wang X, Tang Q, Qi X M, Su S F, Xiang Z H, Hu J H 2023 Mech. Syst. Signal Process. 188 109991 doi: 10.1016/j.ymssp.2022.109991
    [21] 王怡, 陈诚, 林书玉 2025 物理学报 74 044303 doi: 10.7498/aps.74.20241326 Wang Y, Chen C, Lin S Y 2025 Acta Phys. Sin. 74 044303 doi: 10.7498/aps.74.20241326
    [22] 高南沙, 张智成, 王谦, 郭鑫羽, 陈克安, 侯宏 2022 科学通报 67 1203 doi: 10.1360/TB-2021-0439 Gao N S, Zhang Z C, Wang Q, Guo X Y, Chen K A, Hou H 2022 Chin. Sci. Bull. 67 1203 doi: 10.1360/TB-2021-0439
    [23] Pelat A, Gautier F, Conlon S C, Semperlotti F 2020 J. Sound. Vib. 476 115316 doi: 10.1016/j.jsv.2020.115316
    [24] Zhang S, Zhou J W, Ma Q X, Huang H Y, Li J W, Yang Y M, Hu J H 2026 Int. J. Mech. Sci. 308 110972
    [25] Jiang H, Ren W B, Ye W X, Lin S Y 2026 Appl. Acoust. 242 111100 doi: 10.1016/j.apacoust.2025.111100
    [26] Chen C, Tang Y F, Ren W B, Wang Y, Guo J Z, Lin S Y 2024 Ultrasonics 143 107417 doi: 10.1016/j.ultras.2024.107417
    [27] 陈诚, 刘洋, 郭建中, 林书玉 2024 中国科学: 物理学 力学 天文学 54 124311 doi: 10.1360/SSPMA-2024-0159 Chen C, Liu Y, Guo J Z, Lin S Y 2024 Sin, -Phys, Mech, Astron. 54 124311 doi: 10.1360/SSPMA-2024-0159
    [28] Chen C, Liu Y, Wang C H, Guo J Z, Lin S Y 2024 Ultrason. Sonochem. 111 107106 doi: 10.1016/j.ultsonch.2024.107106
    [29] Wang Y, Chen C, Ren W B, Chen H Q, Lin S Y 2025 Sens. Actuator A Phys. 395 117042 doi: 10.1016/j.sna.2025.117042
    [30] Chen H Q, Wang Y, Jiang H, Chen C, Lin S Y 2026 Sens. Actuator A Phys. 400 117500 doi: 10.1016/j.sna.2026.117500
    [31] Zhang X L, Liang B 2018 Appl. Acoust. 129 284 doi: 10.1016/j.apacoust.2017.08.009
    [32] 张小丽, 林书玉, 付志强, 王勇 2013 物理学报 62 034301 doi: 10.7498/aps.62.034301 Zhang X L, Lin S Y, Fu Z Q, Wang Y 2013 Acta Phys. Sin. 62 034301 doi: 10.7498/aps.62.034301
    [33] Feng F L, Shen J Z, Deng J J 2006 Ultrasonics 44 e723 doi: 10.1016/j.ultras.2006.05.083
    [34] Chen C, Guo J Z, Lin S Y 2024 Ultrasonics 138 107248 doi: 10.1016/j.ultras.2024.107248
  • 加载中
图( 14) 表( 1)
计量
  • 文章访问数:  319
  • HTML全文浏览数:  319
  • PDF下载数:  4
  • 施引文献:  0
出版历程
  • 收稿日期:  2025-10-24
  • 刊出日期:  2026-04-20

基于声黑洞-Moonie换能器结构的微型化超声手术刀

    通讯作者: E-mail: guojz@snnu.edu.cn.; 
    通讯作者: E-mail: sylin@snnu.edu.cn.
  • 1. 长安大学理学院, 西安 710064
  • 2. 陕西师范大学, 陕西省超声重点实验室, 西安 710119

摘要: 本文提出了一种基于声黑洞(acoustic black hole, ABH)与Moonie换能器结构的微型化超声手术刀, 旨在利用ABH结构的能量局域聚集效应以及Moonie换能器的位移放大特性有效解决超声手术刀在微型化过程中振动性能衰减问题. 该手术刀由声黑洞轮廓金属端帽、工作刀头与纵向极化压电圆环组成, 并实现了多重振动模式转换. 利用等效电路法建立了ABH-Moonie型手术刀多模态耦合振动的解析理论模型, 实现了快速频率预测. 采用有限元法分析了ABH-Moonie型手术刀的耦合振动特性, 并探究了声黑洞参数和负载对振动性能的影响. 相较于普通Moonie结构手术刀(m < 2), ABH-Moonie型手术刀最大振幅输出显著提升. 在负载条件下, 手术刀仍能保持良好的振动性能, 包括较小的频率漂移以及振动减幅. 最后加工了手术刀样机并进行实验测试, 实测手术刀的设计模态能够被有效激发, 频率误差仅为1.5%, 充分验证了ABH-Moonie复合增效结构的可行性以及解析理论模型的准确性. 本研究可以为微型化高性能压电器件的设计优化提供全新视角.

English Abstract

    • 超声手术刀作为现代外科手术中的重要工具, 以其精确切割、同步止血和对周围组织损伤小等优点得到了广泛应用[1–5]. 随着微创手术、智能手术机器人和精准医疗的快速发展, 对超声手术刀的微型化需求日益迫切. 然而, 传统Langevin型超声手术刀在微型化过程中面临一个根本性矛盾: 尺寸减小通常导致振幅下降和声能辐射效率降低, 进而影响手术切割效果和止血能力[6]. 因此, 探索微型化超声手术刀的性能强化方案是目前该领域的重要研究课题之一.

      为克服传统Langevin型超声手术刀的设计局限, 研究者们相继提出了多种微型化超声手术刀设计方案. 例如, Hu等[7]基于刻槽设计提出了一种纵弯耦合振动小型化超声手术刀, 通过耦合弯振模态带来横向振幅增益, 使得手术刀的振动能力得到了显著提升; Hofmann等[8]提出一种平面钛基微型化超声手术刀, 该手术刀采用d31振动模式在10.2 W功率驱动下使切割幅度达到25.3 μm; Kuang等[9]基于PMN-PT压电单晶贴片设计提出了一种高性能超声微刀, 并对其大功率驱动下的振动性能与非线性特征进行了系统研究; Lucas等[10]基于Cymbal弯张换能器结构提出了一种小型化超声骨刀, 该骨刀仅凭单个压电陶瓷片即可实现与传统超声骨刀相当的超声振幅. 这些研究虽技术路线各异, 但均聚焦于结构优化以突破性能瓶颈, 其中基于弯张型换能器的设计理念尤为关键. 而作为弯张型换能器的另一典型代表, Moonie换能器由宾夕法尼亚大学材料研究室于20世纪80年代开发, 其独特的设计结构能将压电陶瓷的径向振动转换为轴向弯曲振动, 从而输出更大的位移振幅. 凭借这一优势, Moonie换能器已在能量收集、压电驱动等多个领域获得广泛应用[11–13]. 因此, 基于Moonie换能器结构对微型化超声手术刀进行性能优化将是一条有效的技术路径.

      另一方面, 声黑洞结构作为一种新兴功能结构在多个领域展现出广阔前景, 包括减振降噪[14–16]、信号放大[17]、能量收集[18]以及粒子操控等[19–21]. 该结构的核心机理在于利用特定的厚度分布设计(通常为幂函数形式), 实现对弹性体内弯曲波的调控, 从而在特定区域形成能量聚集带, 实现对声波的有效捕获与聚焦[22,23]. 值得注意的是, 声黑洞因其独特的薄壁结构, 近年来已迅速扩展到微型机电系统设计领域. 例如, Zhang等[24]利用径向声黑洞结构激发声流进行高效冷却, 在低功耗下获得高能量利用系数, 为微型器件的散热问题提供了新思路. Jiang等[25]通过设计双声黑洞轮廓综合提升了微型换能器的空气辐射声压与电声转换效率. 这些不同领域的进展共同表明, 声黑洞结构在应对微型器件能量密度不足、转换效率低等共性挑战方面具有独特优势, 也为本研究提供了跨领域依据与原理支撑.

      在前期工作中, 我们率先将ABH结构引入功率超声振动系统设计中, 开展了一系列开创性研究, 包括: 1)提出了多款集成ABH结构的新型超声手术刀, 显著增强了刀头的振幅和侧向声能辐射[26,27]; 2)开发了新型浸入式声黑洞化学反应器, 有效挖掘了探头的空化潜力, 提升了声处理范围与效率[28]; 3)提出了二维声黑洞结构超声镊子辐射器, 在增强辐射的同时实现了对微粒的多功能精准操控[29,30]. 这些研究充分验证了ABH结构在功率超声器件中强化振动、促进空化与提升声操控能力的综合性能优势, 也为构建高能量密度的振动系统提供了全新思路.

      在上述研究背景下, 本文的核心工作是继续将声黑洞结构的能量局域聚焦效应与Moonie换能器的位移放大特性相结合, 在有限空间内强化局部振动能量密度, 从而有效解决超声手术刀在微型化过程中振动性能衰减问题. 主要研究内容包括: 首先, 建立了ABH-Moonie型手术刀的多模耦合振动机电等效电路模型; 其次, 采用有限元法分析了手术刀的频率响应与耦合振动特性, 并探究了声黑洞参数变化和负载施加对手术刀性能的影响; 最后, 加工了手术刀样机, 并开展了阻抗分析和振动测量实验.

    • 本文提出的ABH-Moonie型手术刀几何结构如图1(a), (b)所示, ABH-Moonie型手术刀由纵向极化压电陶瓷圆环、金属圆环、ABH圆板和工作刀头组成, ABH圆板的厚度h随着半径r在$ 0 \lt r \lt {r}_{\text{ABH}} $处按照幂函数规律变化, 且m ≥ 2, h(r)的表达式为

      在图1(c)中, r0, R1分别是压电陶瓷的内半径和工作刀头的半径, r1是ABH圆板、金属圆环和压电陶瓷的外半径, rABH, h0和h1是二维声黑洞轮廓的半径、初始厚度和最大厚度, h2和h3分别是金属圆环和压电陶瓷的厚度, L1是工作刀头的长度.

      ABH-Moonie型手术刀的基本设计原理如下: 纵向极化的压电陶瓷圆环通过逆压电效应和泊松效应产生径向-纵向耦合振动, 激励金属圆环产生纵向振动, 进而经由环形接触面驱动具有二维ABH结构的圆板产生弯曲振动. 在声黑洞效应作用下, 圆板中的弯曲振动振幅自边缘向中心逐渐增强, 并在中心区域达到最大值. 理论上, 当h0趋近于零时, 弯曲波幅值可趋于无穷大. 此外, 由于结构刚度的降低, ABH圆板的弯曲共振频率相应下降, 有利于在较低频超声范围内实现结构复合, 并提升振动系统的能量转换效率. 工作刀头与ABH圆板中心相连接, 从而将该处的弯曲振动转换为刀头的纵向振动. 工作刀头的设计也可进行灵活调整, 为了满足超声牙科、骨科的特殊形状刀头要求, 图1(a)中的工作刀头可进行套扣处理, 从而连接其他异形刀头.

    • 由于ABH圆板的半径远大于厚度, 在自由边界条件下, 弯曲振动ABH圆板的等效集中参数可由等厚度薄圆盘弯曲振动的集中参数近似得出, 在圆盘的等效柔顺常数Cn和等效质量Mn的表达式中, 用ABH圆板的等效厚度heq代替圆盘的厚度h, 即可近似求出ABH圆板的等效质量Mn-ABH和等效柔顺常数Cn-ABH.

      基于薄板经典理论, 在自由边界条件下, 等厚度圆形薄板弯曲振动的横向位移ζ(r, t)可以根据线弹性理论写成[31,32]:

      其中, j是虚数单位, A和B是由边界条件决定的两个常数, J0(kr)和I0(kr)分别是零阶第一类贝塞尔函数和零阶第一类修正贝塞尔函数. 在与圆板相关的参数中, $k^4 = \rho h\omega^2/D$, $D = Eh^3/[12(1 - \sigma^2)] $, k为波数, ρ为密度, h为厚度, ω为角频率, E为杨氏模量, σ为泊松比.

      圆板弯曲振动的弯矩和剪切力可以表示为

      考虑自由边界条件时, (Mr)a = 0, (Qr)a = 0, 其中a是圆板的半径:

      第n阶弯曲振动的位移分布和振动速度的本征函数可表示为

      考虑圆盘微小的体积元dv = 2πrhρdr, 整个圆板的动能Ekn和势能Epn可表示为

      圆板中心处的动能$ E_{\text{k}n}^{\prime} $和势能$ E_{\text{p}n}^{\prime} $为

      选取中心处为动能和势能的参考点, 则有:

      因此, 圆板弯曲振动的等效质量Mn和等效弹性系数Kn可以表达成:

      其中m1是弯曲振动圆盘的质量. 代入边界条件(5)式和(6)式后, (17)式和(18)式可以写成:

      在此基础上即可推导ABH圆板的等效质量Mn-ABH和等效柔顺系数Cn-ABH. ABH圆板的厚度表达式如下:

      体积和等效厚度为

      将(24)式得出的heq代替h代入(19)式, (20)式和(21)式就能求出ABH圆板弯曲振动的等效集中参数Cn-ABH, Mn-ABH, 进而求得圆板弯曲振动的等效机械阻抗可以表达成:

      最后可得圆板弯曲振动频率方程为

      为了验证上述推导的准确性, 利用有限元法对其频率计算结果进行了初步验证. 将ABH圆板的几何尺寸取为: r1 = 15 mm, rABH = 12 mm, h0 = 1 mm, h1 = 3 mm, m = 2.5, 材料选择为铝, 密度为2700 kg/m3, 杨氏模量为7×1010 N/m2, 泊松比为0.33. 根据(26)式计算得到的ABH圆板设计模态的弯曲共振频率为62345 Hz, 有限元分析所得对应模态的共振频率为63130 Hz, 二者相对误差仅为1.3%, 证明了等效厚度近似方法的可行性. 需要说明的是, 该方法仅适用于薄板模型, 且要求厚度变化相对平缓.

    • 基于传输线理论、力和振速的连续性条件, 建立了ABH-Moonie型手术刀耦合振动的机电等效电路模型, 如图2所示. 红色部分代表纵向极化压电陶瓷圆环纵-径耦合振动的多端机电等效电路, 其中Zz1和Zz2为纵向振动机械阻抗, Zr1, Zr2和Zr3为径向振动机械阻抗, V和I为施加的激励电压与电流, C0为等效电容, nz和nr为等效机械变压器, Zc0为纵-径向弹性耦合阻抗[33]. 灰色部分是金属圆环纵向振动的T型网络, 其中Z1和Z2为等效机械阻抗. 蓝色部分代表ABH圆板弯曲振动的等效集中参数电路, Cn-ABH为等效柔顺常数, Mn-ABH为等效质量. 黄色部分代表工作刀头纵向振动的T型网络, 其中Z3和Z4为等效机械阻抗. 图2中各个元件具体表达式如下[34]:

      其中, kPr = ω/cPr为压电陶瓷径向波数, $k_{\rm Pz} = \omega/c_{\rm Pz} $为压电陶瓷纵向波数, $c_{\rm Pr} = (c_{11}^{\text{E}}/\rho_0)^{1/2} $为压电陶瓷径向声速, $c_{\rm Pr} =( c_{33}^{\text{E}}/\rho_0)^{1/2} $为压电陶瓷纵向声速, ρ0为压电陶瓷密度. 进一步化简电路可得

      利用网孔电流法, 可得出该电路的基尔霍夫电压定律方程如下:

      根据(39)式—(42)式, 求解可得到v1, v2, v3和v4的表达式. 最后可以得到ABH-Moonie型手术刀的等效输入电导纳为

      基于(43)式, 可得ABH-Moonie型手术刀的共振频率方程为

      求解上式即可得到ABH-Moonie型手术刀的共振频率.

    • 基于商业软件COMSOL Multiphysics 5.4的固体力学和静电模块, 利用有限元法对ABH-Moonie型手术刀的动力学特性进行了数值建模与探究. ABH-Moonie型手术刀几何参数R1, r0, rABH, r1, h0, h1, h2, h3, L1分别为1, 3, 12, 15, 1, 3, 2, 3, 5 mm. 其余材料参数如表1所示. 建立了二维轴对称仿真模型以计算手术刀电导纳频率响应特性以及耦合振动模态. 为了确保计算精度, 网格划分选择了极细化. 激励电压为1 V, 并且忽略负载和损耗的影响. 另外, 利用MATLAB R2018b数值计算软件对等效电路法得到的频率方程(45)进行求解. 等效电路法与有限元法的频率计算范围均为10—100 kHz, 计算步长为1 Hz.

      图3显示了通过等效电路法和有限元法获得的ABH-Moonie型手术刀的电导纳频率响应曲线. 对于目标设计模态, ABH-Moonie型手术刀的理论计算共振频率为66645 Hz, 略高于有限元计算的64357 Hz, 误差仅为3.5%. 误差主要源于理论分析模型中的边界条件近似连续, 并且忽略了其他非设计振动模式的耦合效应. 相比有限元法, 解析理论模型的建立不仅实现了快速频率预估, 更可以用来辅助辨别耦合振动的主导因素, 从而更准确把握设计原理和潜在物理机制.

      接着, 对ABH-Moonie型手术刀在64357 Hz共振频率下的耦合振动模态进行分析, 如图4所示. 可以看出, 底部压电陶瓷的径向振动较为明显. 受泊松效应影响, 压电陶瓷同时产生纵向振动, 进而激发金属圆环作纵向振动. 金属圆环进一步激励ABH圆板产生弯曲振动, 该圆板呈现出两个振动节圆, 且弯曲振幅在中心位置达到最大值. 最终, ABH圆板驱动位于中心的工作刀头沿轴向振动, 完成了径向-纵向-弯曲-纵向的多重振动模式转换. 此外, 得益于ABH结构与Moonie换能器复合结构的协同增效, 振动能量被聚焦在轴心位置, 从而使得刀头整体位移显著增大.

    • 本节继续探究声黑洞参数对ABH-Moonie型手术刀振动性能的影响. 作为声黑洞结构的核心参数, 截面厚度变化指数m是调控振幅大小的关键[22]. 因此, 选取声黑洞参数m作为自变量, 分别探究其与共振/反共振频率(fr, fa)、有效机电耦合系数keff以及刀头最大输出振幅的变化关系. 其中有效机电耦合系数被定义为

      由图5可知, 随着声黑洞参数m的增大, ABH-Moonie型手术刀的共振/反共振频率、有效机电耦合系数keff均呈下降的趋势, 其中在m < 2.5时, 有效机电耦合系数keff的下降趋势较快, 在m > 2.5时, 有效机电耦合系数keff的变化趋势趋于平缓, 幅值保持在0.354附近. 接着, 提取了刀头尖端最大振动位移, 并分析了振幅大小与声黑洞参数m的变化关系. 由图6可见, 当m ≥ 2时, ABH-Moonie型手术刀的最大输出振幅明显高于非声黑洞设计的普通Moonie结构手术刀. 以m分别为1.5和2.5的两种结构为例进行比较, 结果显示, 相较于普通Moonie结构(m = 1.5), ABH-Moonie型手术刀的最大输出振幅提高了约45%. 综上分析可得, 声黑洞设计虽然会小幅降低机电耦合系数, 但其带来的振幅显著提升对微型化超声手术刀的性能补偿尤为关键. 另外, 共振频率的降低亦有利于器件小型化或匹配特定的生物力学频率.

      此外, 随着参数m的增大, ABH-Moonie型手术刀刀头的最大振幅呈现先迅速上升后逐渐下降的趋势. 由此可见, 在声黑洞设计范围(m ≥ 2)内, 最大振幅输出与参数m并非简单的单调递增关系, 其主要原因在于截断厚度(即初始厚度h0)的影响. 在实际加工过程中, 手术刀的截断厚度不可能为零, 而非零截断厚度会导致弯曲波在边界处发生反射, 反射波与入射波相互干涉, 从而影响振幅大小. 因此, 在给定截断厚度h0后, 选取适配的参数m将是实现手术刀更大振幅输出的关键.

    • 为了模拟ABH-Moonie型手术刀的实际工况, 进一步研究了仿体组织负载对手术刀振动特性的影响. 建立了有限元模型以分析手术刀作用于仿体组织时的振动特性, 如图7所示. 刀头与仿体组织模型的轴向接触长度为3 mm. 所选取的仿体组织模型为线弹性材料, 忽略了黏弹性等非线性因素的影响. 仿体组织的具体材料与尺寸参数如下: 密度为1000 kg/m3, 泊松比为0.4955, 等效弹性模量为25 kPa, 长度为20 mm, 宽度为3 mm, 高度为5 mm.

      图8显示了仿体组织负载对ABH-Moonie型手术刀电导纳频率响应特性的影响. 当手术刀作用于仿体组织时, 其共振频率从无负载时的64357 Hz下降至64355 Hz, 频率漂移仅为2 Hz, 这在一定程度上将降低动态跟踪频率和电学匹配的难度.

      图9显示了仿体组织负载对手术刀振动位移幅值的影响. 可以看出, 负载的存在会降低手术刀的最大输出幅度, 最大减幅仅为7.4%, 满足实际工程应用需求. 总体来说, 对于仿体组织负载的影响, ABH-Moonie型手术刀仍能保持良好的振动性能, 具体表现为较小的频率漂移与振幅减幅.

    • 图10所示为加工ABH-Moonie型手术刀样机, 并在小信号激励下对其电学与振动特性进行实验测量. 手术刀的工作刀头、ABH圆板、金属圆环为一体化加工, 材料为铝, 压电陶瓷材料选用了PZT-4. 另外, 采用环氧树脂导电胶对金属和压电陶瓷进行粘接. 在粘接过程中通过对厚度和固化压力的控制, 确保样机的结构强度和功率测试的稳定性. 手术刀的径向最大尺寸仅为30 mm, 轴向最大尺寸仅为13 mm.

    • 利用WK6500 B精密阻抗分析仪测量了ABH-Moonie型手术刀样机的电学特性, 如图11所示. 频率扫描范围为60—75 kHz, 激励电压为1 V, 测量点数选择为400. 图12显示了ABH-Moonie型手术刀样机的电阻抗频率响应和相位曲线, 测得其共振频率为65639 Hz, 与理论设计结果的误差小于1.5%, 与仿真设计结果的误差小于2%, 验证了解析理论与仿真模型的准确性. 误差的主要来源如下: 1)理论建模中存在近似处理, 包括ABH圆板的等效高度近似以及各部件的边界条件连续近似; 2)实际使用的材料参数与标准材料参数存在差异; 3)实验测量中不可忽略机械损耗、介电损耗以及负载效应; 4)手术刀的加工采用胶粘工艺, 胶体间的微小气泡难于完全去除导致频率误差出现.

    • 使用全场激光扫描振动测量系统(Polytec PSV-400)对ABH-Moonie型手术刀样机的振动模态进行了测量, 如图13所示, 频率扫描范围为60—75 kHz, 激励电压为1 V. 图14显示了测量得到的ABH-Moonie型手术刀的耦合振动模态, 对应的共振频率分别为65672 Hz. 二维ABH结构设计使得微型手术刀的弯曲振幅在短半径圆板上累积强放大, 能量密度在中心位置达到最大, 验证了复合增效结构的实际可行性.

    • 本文通过融合ABH结构的能量局域聚集效应以及Moonie换能器的位移放大机制, 提出了一种新型多模态耦合振动微型化超声手术刀, 以满足现代外科手术对器械微型化与高性能的迫切需求. 研究采用解析建模、数值仿真与实验测试相结合的方法, 获得如下主要结论.

      1)建立的机电等效电路模型可准确预测手术刀频率特性, 为器件的性能优化与参数设计提供了理论依据.

      2)手术刀实现了径向-纵向-弯曲-纵向的多重振动模式转换, 仿真与实测数据吻合良好. ABH-Moonie复合结构设计使得手术刀的振动位移显著提高.

      3)声黑洞结构参数对手术刀的机电转换效率及最大输出振幅具有显著影响, 实际应用中需结合具体场景需求予以合理选择.

      4)在负载条件下, ABH-Moonie型手术刀仍能保持良好的振动性能, 包括较小的频率漂移以及振动减幅.

      5)ABH-Moonie换能器复合增效结构在强化微型化压电器件性能方面具有巨大潜力, 未来可应用于微创手术器械、微纳操控、精密驱动等多个领域.

    参考文献 (34)

目录

/

返回文章
返回