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

基于2-2型压电复合材料的新型宽频带径向振动超声换能器

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

摘要:本文提出了一种基于2-2型压电复合材料的新型宽频带径向振动超声换能器, 它主要由内金属圆环和外压电陶瓷复合材料圆环组成. 首先利用Newnham串并联理论和均匀场理论推导了2-2型压电复合材料的等效参数; 其次利用解析法得到了金属圆环和径向极化压电复合陶瓷圆环径向振动的机电等效电路; 最后得到了换能器的六端机电等效电路, 从而得到了换能器的频率方程. 接着分析了换能器共振频率和反共振频率以及有效机电耦合系数与几何尺寸、两相体积占比的关系, 采用仿真软件对新型换能器的径向振动进行了数值模拟. 结果表明, 利用解析法得到的共振频率和反共振频率与数值模拟结果吻合较好. 此外, 对换能器在水下的辐射声场进行了仿真研究, 结果表明新型复合材料径向换能器相比传统纯陶瓷径向换能器, 发射电压响应幅值更大, 工作带宽提高接近一倍, 声匹配更佳.
关键词: 复合材料等效参数/
径向振动/
带宽/
声匹配

English Abstract


--> --> -->
径向振动换能器具有辐射面积大、辐射效率高、径向辐射均匀, 作用范围广等优点, 因此被广泛应用于水声、超声降解以及声化学等超声液体处理技术领域[1]. 径向换能器种类主要有以下几种: 圆管型压电陶瓷换能器[2]、复合型圆管压电换能器[3]和径向复合压电超声换能器等[4,5]. 另一方面, 2-2型压电复合材料是现在研究最具体、发展前景最好的压电复合材料之一. 这种新型材料相较于传统纯压电陶瓷, 具有低阻抗、低机械品质因数、频带宽等优点[6,7].
本文提出了一种基于2-2型压电复合材料的新型径向振动超声换能器, 如图1所示, 它是由内金属圆环和外压电复合陶瓷圆环组成. 首先推导了2-2型压电复合材料的等效参数, 并利用解析法推导了金属圆环和径向极化压电陶瓷复合材料圆环径向振动的机电等效电路和换能器的频率方程. 接着采用仿真软件对新型换能器的径向振动性能进行数值仿真模拟, 并对结果进行分析. 最后研究了新型径向换能器较传统径向纯陶瓷换能器在工作带宽和声匹配上的性能提升.
图 1 新型径向振动复合材料圆环换能器
Figure1. New radial vibration composite ring transducer.

图2所示, 取z方向为极化方向, 压电相为横观各向同性体, 聚合物相为各向同性体, 聚合物相无压电效应.
图 2 2-2型压电复合材料
Figure2. 2-2 piezoelectric composite.

压电相本构关系为[8,9]
$\begin{split}\left[ {\begin{array}{*{20}{c}} {T_1^{\rm{c}}} \\ {T_2^{\rm{c}}} \\ {T_3^{\rm{c}}} \\ {T_4^{\rm{c}}} \\ {T_5^{\rm{c}}} \\ {T_6^{\rm{c}}} \end{array}} \right] =\;& \left[ {\begin{array}{*{20}{c}} {C_{11}^{\rm{c}}}&{C_{12}^{\rm{c}}}&{C_{13}^{\rm{c}}}&0&0&0 \\ {C_{12}^{\rm{c}}}&{C_{11}^{\rm{c}}}&{C_{13}^{\rm{c}}}&0&0&0 \\ {C_{13}^{\rm{c}}}&{C_{13}^{\rm{c}}}&{C_{33}^{\rm{c}}}&0&0&0 \\ 0&0&0&{C_{44}^{\rm{c}}}&0&0 \\ 0&0&0&0&{C_{44}^{\rm{c}}}&0 \\ 0&0&0&0&0&{C_{66}^{\rm{c}}} \end{array}} \right] \\&\times \!\left[\! {\begin{array}{*{20}{c}} {S_1^{\rm{c}}} \\ {S_2^{\rm{c}}} \\ {S_3^{\rm{c}}} \\ {S_4^{\rm{c}}} \\ {S_5^{\rm{c}}} \\ {S_6^{\rm{c}}} \end{array}}\! \right]\! -\! \left[\! {\begin{array}{*{20}{c}} 0&0&{{e_{31}}} \\ 0&0&{{e_{31}}} \\ 0&0&{{e_{33}}} \\ 0&{{e_{15}}}&0 \\ {{e_{15}}}&0&0 \\ 0&0&0 \end{array}}\! \right]\! \times\! \left[ \!{\begin{array}{*{20}{c}} {E_1^{\rm{c}}} \\ {E_2^{\rm{c}}} \\ {E_3^{\rm{c}}} \end{array}}\! \right], \end{split}$
$\begin{split}\left[ {\begin{array}{*{20}{c}} {D_1^{\rm{c}}} \\ {D_2^{\rm{c}}} \\ {D_3^{\rm{c}}} \end{array}} \right] =\;& \left[ {\begin{array}{*{20}{c}} 0&0&0&0&{{e_{15}}}&0 \\ 0&0&0&{{e_{15}}}&0&0 \\ {{e_{31}}}&{{e_{31}}}&{{e_{33}}}&0&0&0 \end{array}} \right]\\&\times \!\left[ \!{\begin{array}{*{20}{c}} {S_1^{\rm{c}}} \\ {S_2^{\rm{c}}} \\ {S_3^{\rm{c}}} \\ {S_4^{\rm{c}}} \\ {S_5^{\rm{c}}} \\ {S_6^{\rm{c}}} \end{array}}\! \right] \!+\! \left[ \!{\begin{array}{*{20}{c}} {\varepsilon _{11}^{\rm{c}}}&0&0 \\ 0&{\varepsilon _{11}^{\rm{c}}}&0 \\ 0&0&{\varepsilon _{33}^{\rm{c}}} \end{array}} \!\right] \!\times\! \left[ \!{\begin{array}{*{20}{c}} {E_1^{\rm{c}}} \\ {E_2^{\rm{c}}} \\ {E_3^{\rm{c}}} \end{array}} \!\right], \end{split}$
其中$T_i^{\rm{c}}\;(i \!=\! 1, 2, 3, \cdots 6)$为应力, $S_i^{\rm{c}}\;(i = 1, 2, 3, \cdots 6)$为应变, $E_{\rm{i}}^{\rm{c}}\;(i \!=\! 1, 2, 3)$为电场强度分量, $D_i^{\rm{c}}(i \!=\! 1, 2, 3)$为电位移分量, $C_{ij}^{\rm{c}}(i, j \!=\! 1, 2, 3, \cdots 6)$为常电场刚度系数, $e_{{\rm{ij}}}^{\rm{c}}(i, j = 1, 2, 3, \cdots 6)$为常应变压电应力系数, $\varepsilon _{ij}^{\rm{c}}(i, j = 1, 2, 3)$为常应变介电系数, 上标c代表压电相, 上标p代表聚合物相, 上标cp代表复合材料.
聚合物相本构关系为
$\left[\! {\begin{array}{*{20}{c}} {T_1^{\rm{p}}} \\ {T_2^{\rm{p}}} \\ {T_3^{\rm{p}}} \\ {T_4^{\rm{p}}} \\ {T_5^{\rm{p}}} \\ {T_6^{\rm{p}}} \end{array}}\! \right] \!\!=\!\!\left[\!\! {\begin{array}{*{20}{c}} {C_{11}^{\rm{p}}}\!&\!{C_{12}^{\rm{p}}}\!&\!{C_{{\rm{12}}}^{\rm{p}}}\!&\!0\!&\!0\!&\!0\! \\ {C_{12}^{\rm{p}}}\!&\!{C_{11}^{\rm{p}}}\!&\!{C_{{\rm{12}}}^{\rm{p}}}\!&\!0\!&\!0\!&\!0\! \\ {C_{{\rm{12}}}^{\rm{p}}}\!&\!{C_{{\rm{12}}}^{\rm{p}}}\!&\!{C_{{\rm{11}}}^{\rm{p}}}\!&\!0\!&\!0\!&\!0\! \\ \!0\!&\!0\!&\!0\!&\!{C_{{\rm{66}}}^{\rm{p}}}\!&\!0\!&\!0\! \\ \!0\!&\!0\!&\!0\!&\!0\!&\!{C_{{\rm{66}}}^{\rm{p}}}\!&\!0\! \\ \!0\!&\!0\!&\!0\!&\!0\!&\!0\!&\!{C_{{\rm{66}}}^{\rm{p}}} \!\end{array}}\!\! \right] \!\times\! \left[\! {\begin{array}{*{20}{c}} {S_1^{\rm{p}}} \\ {S_2^{\rm{p}}} \\ {S_3^{\rm{p}}} \\ {S_4^{\rm{p}}} \\ {S_5^{\rm{p}}} \\ {S_6^{\rm{p}}} \end{array}} \!\right], $
$\left[ {\begin{array}{*{20}{c}} {D_1^{\rm{p}}} \\ {D_2^{\rm{p}}} \\ {D_3^{\rm{p}}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {\varepsilon _{11}^{\rm{p}}}&0&0 \\ 0&{\varepsilon _{11}^{\rm{p}}}&0 \\ 0&0&{\varepsilon _{11}^{\rm{p}}} \end{array}} \right] \times \left[ {\begin{array}{*{20}{c}} {E_1^{\rm{p}}} \\ {E_2^{\rm{p}}} \\ {E_3^{\rm{p}}} \end{array}} \right].$
该复合材料主向(z方向)拉伸时, 不会引起剪切变形的耦合效应, 故可将${T_4} = {T_5} = {T_6} = 0$, ${E_1} = {E_2} = 0$处理, 根据Newnham串并联理论[10]可以写出各个方向应力和应变关系式, 其中${v_{\rm c}}$代表压电相体积分数, ${v_{\rm p}}$代表聚合物相体积分数, ${v_{\rm c}} + {v_{\rm p}} = 1$.
x方向, 压电相和聚合物相串联, 则各相应力相等, 应变为各相应变之和, 即:
$T_1^{\rm{cp}} = T_1^{\rm{c}} = T_1^{\rm{p}},\;S_1^{\rm{cp}} = {v_{\rm{c}}}S_1^{\rm{c}} + {v_{\rm{p}}}S_1^{\rm{p}}.$
y方向, 压电相和聚合物相并联, 则各相应变相等, 应力为各相应力之和, 即:
$T_{\rm{2}}^{\rm{cp}} = {v_{\rm{c}}}T_{\rm{2}}^{\rm{c}} + {v_p}T_{\rm{2}}^{\rm{p}},\;S_2^{\rm{cp}} = S_2^{\rm{c}} = S_2^{\rm{p}}.$
z方向, 压电相和聚合物相并联, 则各相应变相等, 应力为各相应力之和, 电场强度分量相等, 电位移分量为各相电位移分量之和即:
$\begin{split}&T_{\rm{3}}^{\rm{cp}} = {v_{\rm{c}}}T_{\rm{3}}^{\rm{c}} + {v_{\rm{p}}}T_{\rm{3}}^{\rm{p}},\;S_3^{\rm{cp}} = S_3^{\rm{c}} = S_3^{\rm{p}},\;\\&E_3^{\rm{cp}} = E_3^{\rm{c}} = E_3^{\rm{p}},\;D_3^{\rm{cp}} = {v_{\rm{c}}}D_3^{\rm{c}} + {v_p}D_3^{\rm{p}}.\end{split}$
取两相中相等的场量作为独立变量建立本构关系[11], 在压电相中, 取$T_1^{\rm{c}}$, $S_2^{\rm{c}}$, $S_3^{\rm{c}}$, $E_3^{\rm{c}}$为独立变量, 本构方程可表示为
$\left[ {\begin{array}{*{20}{c}} {S_1^{\rm{c}}} \\ {T_2^{\rm{c}}} \\ {T_3^{\rm{c}}} \\ {D_3^{\rm{c}}} \end{array}} \right] = {A^{\rm{c}}} \times \left[ {\begin{array}{*{20}{c}} {T_1^{\rm{c}}} \\ {S_2^{\rm{c}}} \\ {S_3^{\rm{c}}} \\ {E_3^{\rm{c}}} \end{array}} \right], $
${A^{\rm{c}}} = \left| {\begin{array}{*{20}{c}} {\dfrac{1}{{C_{11}^{\rm{c}}}}}&{ - \dfrac{{C_{12}^{\rm{c}}}}{{C_{11}^{\rm{c}}}}}&{ - \dfrac{{C_{13}^{\rm{c}}}}{{C_{11}^{\rm{c}}}}}&{\dfrac{{e_{31}^{\rm{c}}}}{{C_{11}^{\rm{c}}}}} \\ {\dfrac{{C_{12}^{\rm{c}}}}{{C_{11}^{\rm{c}}}}}&{C_{11}^{\rm{c}} - \dfrac{{C{{_{12}^{\rm{c}}}^2}}}{{C_{11}^{\rm{c}}}}}&{C_{13}^{\rm{c}} - \dfrac{{C_{12}^{\rm{c}}C_{13}^{\rm{c}}}}{{C_{11}^{\rm{c}}}}}&{\dfrac{{C_{12}^{\rm{c}}}}{{C_{11}^{\rm{c}}}}e_{31}^{\rm{c}} - e_{31}^{\rm{c}}} \\ {\dfrac{{C_{13}^{\rm{c}}}}{{C_{11}^{\rm{c}}}}}&{C_{13}^{\rm{c}} - \dfrac{{C_{12}^{\rm{c}}C_{13}^{\rm{c}}}}{{C_{11}^{\rm{c}}}}}&{C_{33}^{\rm{c}} - \dfrac{{C{{_{13}^{\rm{c}}}^2}}}{{C_{11}^{\rm{c}}}}}&{\dfrac{{C_{13}^{\rm{c}}}}{{C_{11}^{\rm{c}}}}e_{31}^{\rm{c}} - e_{33}^{\rm{c}}} \\ {\dfrac{{e_{31}^{\rm{c}}}}{{C_{11}^{\rm{c}}}}}&{e_{31}^{\rm{c}} - \dfrac{{C_{12}^{\rm{c}}}}{{C_{11}^{\rm{c}}}}e_{31}^{\rm{c}}}&{e_{33}^{\rm{c}} - \dfrac{{C_{13}^{\rm{c}}}}{{C_{11}^{\rm{c}}}}e_{31}^{\rm{c}}}&{\dfrac{{e{{_{31}^{\rm{c}}}^2}}}{{C_{11}^{\rm{c}}}} + \varepsilon _{33}^{\rm{c}}} \end{array}} \right| .$
在聚合物相中, 取$T_1^{\rm{p}}$, $S_2^{\rm{p}}$, $S_3^{\rm{p}}$, $E_3^{\rm{p}}$为独立变量, 本构方程可表示为
$\left[ {\begin{array}{*{20}{c}} {S_1^{\rm{p}}} \\ {T_2^{\rm{p}}} \\ {T_3^{\rm{p}}} \\ {D_3^{\rm{p}}} \end{array}} \right] = {A^{\rm{p}}} \times \left[ {\begin{array}{*{20}{c}} {T_1^{\rm{p}}} \\ {S_2^{\rm{p}}} \\ {S_3^{\rm{p}}} \\ {E_3^{\rm{p}}} \end{array}} \right], $
${A^{\rm{p}}} = \left| {\begin{array}{*{20}{c}} {\dfrac{1}{{C_{11}^{\rm{p}}}}}&{ - \dfrac{{C_{12}^{\rm{p}}}}{{C_{11}^{\rm{p}}}}}&{ - \dfrac{{C_{12}^{\rm{p}}}}{{C_{11}^{\rm{p}}}}}&0 \\ {\dfrac{{C_{12}^{\rm{p}}}}{{C_{11}^{\rm{p}}}}}&{C_{11}^{\rm{p}} - \dfrac{{C{{_{12}^{\rm{p}}}^2}}}{{C_{11}^{\rm{p}}}}}&{C_{12}^{\rm{p}} - \dfrac{{C{{_{12}^{\rm{p}}}^2}}}{{C_{11}^{\rm{p}}}}}&0 \\ {\dfrac{{C_{12}^{\rm{p}}}}{{C_{11}^{\rm{p}}}}}&{C_{12}^{\rm{p}} - \dfrac{{C{{_{12}^{\rm{p}}}^2}}}{{C_{11}^{\rm{p}}}}}&{C_{11}^{\rm{p}} - \dfrac{{C{{_{12}^{\rm{p}}}^2}}}{{C_{11}^{\rm{p}}}}}&0 \\ 0&0&0&{\varepsilon _{11}^{\rm{p}}} \end{array}} \right| .$
在复合材料中, 取$T_1^{\rm{cp}}$, $S_2^{\rm{cp}}$, $S_3^{\rm{cp}}$, $E_3^{\rm{cp}}$为独立变量, 本构方程可表示为
$\begin{split}&\left[ {\begin{array}{*{20}{c}} {S_1^{\rm{cp}}} \\ {T_2^{\rm{cp}}} \\ {T_3^{\rm{cp}}} \\ {D_3^{\rm{cp}}}\end{array}} \right] = {v_{\rm c}}\left[ {\begin{array}{*{20}{c}} {S_1^{\rm{c}}} \\ {T_2^{\rm{c}}} \\ {T_3^{\rm{c}}} \\ {D_3^{\rm{c}}}\end{array}} \right] + {v_{\rm p}}\left[ {\begin{array}{*{20}{c}} {S_1^{\rm{p}}} \\ {T_2^{\rm{p}}} \\ {T_3^{\rm{p}}} \\ {D_3^{\rm{p}}}\end{array}} \right] = {v_{\rm c}}{A^{\rm{c}}}\left[ {\begin{array}{*{20}{c}} {T_1^{\rm{c}}} \\ {S_2^{\rm{c}}} \\ {S_3^{\rm{c}}} \\ {E_3^{\rm{c}}}\end{array}} \right] \\&+ {v_{\rm p}}{A^{\rm{p}}} \times \left[ {\begin{array}{*{20}{c}} {T_1^{\rm{p}}} \\ {S_2^{\rm{p}}} \\ {S_3^{\rm{p}}} \\ {E_3^{\rm{p}}}\end{array}} \right] = ({v_{\rm c}}{A^{\rm{c}}} + {v_{\rm p}}{A^{\rm{p}}}) \times \left[ {\begin{array}{*{20}{c}} {T_1^{\rm{cp}}} \\ {S_2^{\rm{cp}}} \\ {S_3^{\rm{cp}}} \\ {E_3^{\rm{cp}}}\end{array}} \right], \end{split}$
压电复合材料本构方程可以表示为
$\left[ {\begin{array}{*{20}{c}} {T_1^{\rm{cp}}} \\ {T_2^{\rm{cp}}} \\ {T_3^{\rm{cp}}} \\ {D_3^{\rm{cp}}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {C_{11}^{\rm{cp}}}&{C_{12}^{\rm{cp}}}&{C_{13}^{\rm{cp}}}&{ - e_{31}^{\rm{cp}}} \\ {C_{12}^{\rm{cp}}}&{C_{22}^{\rm{cp}}}&{C_{23}^{\rm{cp}}}&{ - e_{32}^{\rm{cp}}} \\ {C_{13}^{\rm{cp}}}&{C_{23}^{\rm{cp}}}&{C_{33}^{\rm{cp}}}&{ - e_{33}^{\rm{cp}}} \\ {e_{31}^{\rm{cp}}}&{e_{32}^{\rm{cp}}}&{e_{33}^{\rm{cp}}}&{\varepsilon _{33}^{\rm{cp}}} \end{array}} \right] \times \left[ {\begin{array}{*{20}{c}} {S_1^{\rm{cp}}} \\ {S_2^{\rm{cp}}} \\ {S_3^{\rm{cp}}} \\ {E_3^{\rm{cp}}} \end{array}} \right], $
$\begin{split}&C_{11}^{\rm{cp}} = \dfrac{{C_{11}^{\rm{c}}C_{11}^{\rm{p}}}}{{{v_{\rm{c}}}C_{11}^{\rm{p}} + {v_{\rm{p}}}C_{11}^{\rm{c}}}},\;C_{12}^{\rm{cp}} = \dfrac{{{v_{\rm{c}}}C_{12}^{\rm{c}}C_{11}^{\rm{p}} + {v_{\rm{p}}}C_{12}^{\rm{p}}C_{11}^{\rm{c}}}}{{{v_{\rm{c}}}C_{11}^{\rm{p}} + {v_{\rm{p}}}C_{11}^{\rm{c}}}},\;\\&C_{13}^{\rm{cp}} = \dfrac{{{v_{\rm{c}}}C_{13}^{\rm{c}}C_{11}^{\rm{p}} + {v_{\rm{p}}}C_{12}^{\rm{p}}C_{11}^{\rm{c}}}}{{{v_{\rm{c}}}C_{11}^{\rm{p}} + {v_{\rm{p}}}C_{11}^{\rm{c}}}},\\ &C_{{\rm{22}}}^{\rm{cp}} = {v_{\rm{c}}}\left(C_{11}^{\rm{c}} - \frac{{C{{_{12}^{\rm{c}}}^2}}}{{C_{11}^{\rm{c}}}}\right) + {v_{\rm{p}}}\left(C_{11}^{\rm{p}} - \frac{{C{{_{12}^{\rm{p}}}^2}}}{{C_{11}^{\rm{p}}}}\right)\\&\;\;\qquad+ \frac{{{{({v_{\rm{c}}}C_{12}^{\rm{c}}C_{11}^{\rm{p}} + {v_{\rm{p}}}C_{12}^{\rm{p}}C_{11}^{\rm{c}})}^2}}}{{C_{11}^{\rm{c}}C_{11}^{\rm{p}}({v_{\rm{c}}}C_{11}^{\rm{p}} + {v_{\rm{p}}}C_{11}^{\rm{c}})}},\\ &C_{{\rm{23}}}^{\rm{cp}} = {v_{\rm{c}}}\left(C_{{\rm{13}}}^{\rm{c}} - \frac{{C_{12}^{\rm{c}}C_{1{\rm{3}}}^{\rm{c}}}}{{C_{11}^{\rm{c}}}}\right) + {v_{\rm{p}}}\left(C_{{\rm{12}}}^{\rm{p}} - \frac{{C{{_{12}^{\rm{p}}}^2}}}{{C_{11}^{\rm{p}}}}\right) \\&\qquad\;+ \!\frac{{({v_{\rm{c}}}C_{{\rm{13}}}^{\rm{c}}C_{11}^{\rm{p}} \!+\! {v_{\rm{p}}}C_{12}^{\rm{p}}C_{11}^{\rm{c}})\!({v_{\rm{c}}}C_{{\rm{12}}}^{\rm{c}}C_{11}^{\rm{p}} \!+\! {v_{\rm{p}}}C_{12}^{\rm{p}}C_{11}^{\rm{c}})}}{{C_{11}^{\rm{c}}C_{11}^{\rm{p}}({v_{\rm{c}}}C_{11}^{\rm{p}} + {v_{\rm{p}}}C_{11}^{\rm{c}})}},\\ &C_{{\rm{33}}}^{\rm{cp}} = {v_{\rm{c}}}\left(C_{{\rm{33}}}^{\rm{c}} - \frac{{C{{_{{\rm{13}}}^{\rm{c}}}^2}}}{{C_{11}^{\rm{c}}}}\right) + {v_{\rm{p}}}\left(C_{11}^{\rm{p}} - \frac{{C{{_{12}^{\rm{p}}}^2}}}{{C_{11}^{\rm{p}}}}\right)\\&\qquad\;\;+ \frac{{{{({v_{\rm{c}}}C_{{\rm{13}}}^{\rm{c}}C_{11}^{\rm{p}} + {v_{\rm{p}}}C_{12}^{\rm{p}}C_{11}^{\rm{c}})}^2}}}{{C_{11}^{\rm{c}}C_{11}^{\rm{p}}({v_{\rm{c}}}C_{11}^{\rm{p}} + {v_{\rm{p}}}C_{11}^{\rm{c}})}},\\ &e_{31}^{\rm{cp}} = \frac{{{v_{\rm{c}}}C_{11}^{\rm{p}}e_{31}^{\rm{c}}}}{{{v_{\rm{c}}}C_{11}^{\rm{p}} + {v_{\rm{p}}}C_{11}^{\rm{c}}}},\;e_{{\rm{32}}}^{\rm{cp}} = {v_{\rm{c}}}\left(e_{{\rm{31}}}^{\rm{c}} - \frac{{C_{{\rm{12}}}^{\rm{c}}e_{{\rm{31}}}^{\rm{c}}}}{{C_{11}^{\rm{c}}}}\right) \\&\qquad\;\;+ \frac{{{v_{\rm{c}}}e_{{\rm{31}}}^{\rm{c}}({v_{\rm{c}}}C_{{\rm{12}}}^{\rm{c}}C_{11}^{\rm{p}} + {v_{\rm{p}}}C_{12}^{\rm{p}}C_{11}^{\rm{c}})}}{{C_{11}^{\rm{c}}({v_{\rm{c}}}C_{11}^{\rm{p}} + {v_{\rm{p}}}C_{11}^{\rm{c}})}},\\ &e_{{\rm{33}}}^{\rm{cp}} =\! {v_{\rm{c}}}\!\left(\!e_{{\rm{33}}}^{\rm{c}} \!-\! \frac{{C_{{\rm{13}}}^{\rm{c}}e_{{\rm{31}}}^{\rm{c}}}}{{C_{11}^{\rm{c}}}}\!\right) \!+\! \frac{{{v_{\rm{c}}}e_{{\rm{31}}}^{\rm{c}}({v_{\rm{c}}}C_{{\rm{13}}}^{\rm{c}}C_{11}^{\rm{p}} \!+\! {v_{\rm{p}}}C_{12}^{\rm{p}}C_{11}^{\rm{c}})}}{{C_{11}^{\rm{c}}({v_{\rm{c}}}C_{11}^{\rm{p}} + {v_{\rm{p}}}C_{11}^{\rm{c}})}},\\ &\varepsilon _{{\rm{33}}}^{\rm{cp}} = {v_{\rm{c}}}\left(\!\varepsilon _{{\rm{33}}}^{\rm{c}}\! +\! \frac{{e{{_{{\rm{31}}}^{\rm{c}}}^{\rm{2}}}}}{{C_{11}^{\rm{c}}}}\!\right) \!+\! {v_{\rm{p}}}\varepsilon _{{\rm{11}}}^{\rm{p}} \!-\! \frac{{v_{\rm{c}}^2C_{11}^{\rm{p}}e{{_{{\rm{31}}}^{\rm{c}}}^2}}}{{C_{11}^{\rm{c}}({v_{\rm{c}}}C_{11}^{\rm{p}} + {v_{\rm{p}}}C_{11}^{\rm{c}})}} . \end{split}$

2
3.1.2-2型压电陶瓷复合材料圆环径向振动的等效电路
-->图3是一个径向极化的压电陶瓷复合材料圆环, 其中b, ch分别是圆环内外半径和厚度. 压电本构方程可以表示为
图 3 压电陶瓷复合材料圆环
Figure3. Piezoelectric ceramic composite ring.

${T_\theta } = C_{11}^{\rm{cp}}{S_\theta } + C_{12}^{\rm{cp}}{S_z} + C_{13}^{\rm{cp}}{S_r} - e_{31}^{\rm{cp}}{E_r}, $
${T_Z} = C_{12}^{\rm{pc}}{S_\theta } + C_{22}^{\rm{pc}}{S_z} + C_{23}^{\rm{pc}}{S_r} - e_{32}^{\rm{pc}}{E_r}, $
${T_r} = C_{13}^{\rm{cp}}{S_\theta } + C_{23}^{\rm{cp}}{S_z} + C_{33}^{\rm{cp}}{S_r} - e_{33}^{\rm{cp}}{E_r}, $
${D_r} = e_{31}^{\rm{cp}}{S_\theta } + e_{32}^{\rm{cp}}{S_z} + e_{33}^{\rm{cp}}{S_r} + \varepsilon _{33}^{\rm{cp}}{E_r}, $
其中${T_\theta }$, ${T_z}$, ${T_r}$分别表示切向应力, 纵向应力和径向应力; ${S_\theta }$, ${S_z}$, ${S_r}$分别表示切向应变, 纵向应变和径向应变; ${E_r} = {E_r}(r, t)$是径向激励电场; ${{{D}}_r} = {{{D}}_r}(r, t)$是电位移矢量; $C_{ij}^{\rm{cp}}$, $e_{ij}^{\rm{cp}}$, $\varepsilon _{ij}^{\rm{cp}}$分别代表压电陶瓷复合圆环的等效弹性常数、压电常数以及介电常数. 由于换能器厚度远小于径向尺寸, 因此换能器的振动可以看成平面应力问题[12], 即${T_z} = 0$. 将(13)式变形为
${S_Z} = \frac{{e_{32}^{\rm{cp}}{E_r} - C_{12}^{\rm{cp}}{S_\theta } - C_{23}^{\rm{cp}}{S_r}}}{{C_{22}^{\rm{pc}}}}.$
将(16)式分别代入(12)式、(14)式和(15)式得[13]
$\begin{split}{T_\theta } =\;& \left[ {C_{11}^{\rm{cp}} - \frac{{{{(C_{12}^{\rm{cp}})}^2}}}{{C_{22}^{\rm{cp}}}}} \right]{S_\theta } + \left( {C_{13}^{\rm{cp}} - \frac{{C_{12}^{\rm{cp}}C_{23}^{\rm{cp}}}}{{C_{22}^{\rm{cp}}}}} \right){S_r} \\&- \left( {e_{31}^{\rm{cp}} - \frac{{C_{12}^{\rm{cp}}e_{32}^{\rm{cp}}}}{{C_{22}^{\rm{cp}}}}} \right){E_r},\\[-15pt] \end{split}$
$\begin{split}{T_r} =\;& \left[ {C_{13}^{\rm{cp}} - \frac{{C_{12}^{\rm{cp}}C_{23}^{\rm{cp}}}}{{C_{22}^{\rm{cp}}}}} \right]{S_\theta } + \left( {C_{33}^{\rm{cp}} - \frac{{{{(C_{23}^{\rm{cp}})}^2}}}{{C_{22}^{\rm{cp}}}}} \right){S_r} \\&- \left( {e_{33}^{\rm{cp}} - \frac{{C_{23}^{\rm{cp}}e_{32}^{\rm{cp}}}}{{C_{22}^{\rm{cp}}}}} \right){E_r},\\[-15pt] \end{split}$
$\begin{split}{D_r} =\;& \left[ {e_{31}^{\rm{cp}} - \frac{{e_{32}^{\rm{cp}}C_{12}^{\rm{cp}}}}{{C_{22}^{\rm{cp}}}}} \right]{S_\theta } + \left( {e_{33}^{\rm{cp}} - \frac{{e_{32}^{\rm{cp}}C_{23}^{\rm{cp}}}}{{C_{22}^{\rm{cp}}}}} \right){S_r}\\&+ \left( {\varepsilon _{33}^{\rm{cp}} - \frac{{e_{32}^{\rm{cp}}e_{32}^{\rm{cp}}}}{{C_{22}^{\rm{cp}}}}} \right){E_r}.\\[-15pt]\end{split}$
压电陶瓷复合材料圆环径向振动的波动方程及应变和位移关系可以表示为
$\rho \frac{{{\partial ^{\rm{2}}}{\xi _r}}}{{\partial {t^{\rm{2}}}}} = \frac{{\partial {T_r}}}{{\partial r}} + \frac{{{T_r} - {T_\theta }}}{r}, $
$ {S}_{r}=\frac{\partial {\xi }_{r}}{\partial r},\;{S_\theta } = \frac{{{\xi _r}}}{r}.$
电场强度幅值可以表示为
${E_{r0}} = B{}_3{S_\theta } + {B_5}{S_r} + \frac{{{L_3}}}{r}\frac{{C_{22}^{{\rm{cp}}}}}{{\varepsilon _{33}^{{\rm{cp}}}C_{22}^{{\rm{cp}}} + {{(e_{32}^{{\rm{cp}}})}^2}}}.$
式中, ${L_3}$是一个常数, 可由电边界条件决定,
${B_3} = \frac{{e_{32}^{{\rm{cp}}}C_{12}^{{\rm{cp}}} - e_{31}^{{\rm{cp}}}C_{22}^{{\rm{cp}}}}}{{\varepsilon _{33}^{{\rm{cp}}}C_{22}^{{\rm{cp}}} + {{(e_{32}^{{\rm{cp}}})}^2}}},\;{B_5} = \frac{{e_{32}^{{\rm{cp}}}C_{23}^{{\rm{cp}}} - e_{33}^{{\rm{cp}}}C_{22}^{{\rm{cp}}}}}{{\varepsilon _{33}^{{\rm{cp}}}C_{22}^{{\rm{cp}}} + {{(e_{32}^{{\rm{cp}}})}^2}}}.$
将(22)式代入(17)式和(18)式得
${T_\theta } = {B_1}{S_\theta } + {B_2}{S_r} + {B_3}\frac{{{L_3}}}{r}, $
${T_r} = {B_2}{S_\theta } + {B_4}{S_r} + {B_5}\frac{{{L_3}}}{r}.$
式中,
$\begin{split}{B_1} =\;& \left( {C_{11}^{\rm{cp}} - \frac{{C_{12}^{\rm{cp}}C_{12}^{\rm{cp}}}}{{C_{22}^{\rm{cp}}}}} \right) - \left( {e_{31}^{\rm{cp}} - \frac{{C_{12}^{\rm{cp}}e_{32}^{\rm{cp}}}}{{C_{22}^{\rm{cp}}}}} \right)\\&\times\frac{{e_{32}^{\rm{cp}}C_{12}^{\rm{cp}} - e_{31}^{\rm{cp}}C_{22}^{\rm{cp}}}}{{\varepsilon _{33}^{\rm{cp}}C_{22}^{\rm{cp}} + {{(e_{32}^{\rm{cp}})}^2}}}, \end{split}$
$\begin{split}{B_2} =\;& \left( {C_{13}^{\rm{cp}} - \frac{{C_{12}^{\rm{cp}}C_{23}^{\rm{cp}}}}{{C_{22}^{\rm{cp}}}}} \right) - \left( {e_{31}^{\rm{cp}} - \frac{{C_{12}^{\rm{cp}}e_{32}^{\rm{cp}}}}{{C_{22}^{\rm{cp}}}}} \right)\\&\times\frac{{e_{32}^{\rm{cp}}C_{23}^{\rm{cp}} - e_{33}^{\rm{cp}}C_{22}^{\rm{cp}}}}{{\varepsilon _{33}^{\rm{cp}}C_{22}^{\rm{cp}} + {{(e_{32}^{\rm{cp}})}^2}}}, \end{split}$
$\begin{split}{B_4} =\;& \left( {C_{33}^{\rm{cp}} - \frac{{C_{23}^{\rm{cp}}C_{23}^{\rm{cp}}}}{{C_{22}^{\rm{cp}}}}} \right) - \left( {e_{33}^{\rm{cp}} - \frac{{C_{23}^{\rm{cp}}e_{32}^{\rm{cp}}}}{{C_{22}^{\rm{cp}}}}} \right)\\&\times\frac{{e_{32}^{\rm{cp}}C_{23}^{\rm{cp}} - e_{33}^{\rm{cp}}C_{22}^{\rm{cp}}}}{{\varepsilon _{33}^{\rm{cp}}C_{22}^{\rm{cp}} + {{(e_{32}^{\rm{cp}})}^2}}} .\end{split}$
将(21)式、(23)式和(24)式代入(20)式并分离时间变量得
$\frac{{{{\rm{d}}^2}{\xi _{r0}}}}{{{\rm{d}}{r^2}}} + \frac{1}{r}\frac{{{\rm{d}}{\xi _{r0}}}}{{{\rm{d}}r}} + ({k^2} - \frac{{{v^2}}}{{{r^2}}}){\xi _{r0}} - {B'_3}\frac{{{L_3}}}{{{r^2}}} = 0, $
式中,
$\begin{split}&{v^2} = {B_1}/{B_4},\;{B'_3} = {B_1}/{B_3},\;k = \omega /{V_3},\;\\&V_3^2 = {B_4}/\rho,\;\rho = {v_{\rm{c}}}{\rho _{\rm{c}}} + {v_{\rm{p}}}{\rho _{\rm{p}}},\end{split}$
$\omega $是角频率, $\rho $, ${\rho _c}$${\rho _p}$分别是压电复合材料等效密度, 压电相密度和聚合物相密度, $k$${V_3}$分别是波数和径向极化压电陶瓷复合材料圆环中径向振动的传播速度. 方程(25)是广义贝塞尔方程, 其解为
${\xi _{r0}} = {A_1}{J_v}(kr) + {A_2}{Y_v}(kr) + {B'_3}{L_3}{s_{ - 1,v}}(kr), $
式中, ${J_v}(kr)$, ${Y_v}(kr)$$v$阶的第一类和第二类贝塞尔函数, ${s_{ - 1, v}}(kr)$ 是第一类隆梅尔函数, ${A_1}$, ${A_2}$是两个常数. 根据(26)式可以得出径向振动压电陶瓷复合材料圆环的径向振动速度为
${v_{r0}} = {\rm{j}}\omega [{A_1}{J_v}(kr) + {A_2}{Y_v}(kr) + {B'_3}{L_3}{s_{ - 1,v}}(kr)].$
根据速度边界条件: $v\left| {_{r = b}} \right. = {v_b}$, $v\left| {_{r = c}} \right. = - {v_c}$可得
${A_1} = - \frac{{{\rm{j}}{Y_v}(kc)}}{{\omega {\tau _1}}}{v_b} - \frac{{{\rm{j}}{Y_v}(kb)}}{{\omega {\tau _1}}}{v_c} + \frac{{{\tau _2}}}{{{\tau _1}}}{L_3}, $
${A_2} = \frac{{{\rm{j}}{Y_v}(kc)}}{{\omega {\tau _1}}}{v_b} + \frac{{{\rm{j}}{Y_v}(kb)}}{{\omega {\tau _1}}}{v_c} + \frac{{{\tau _3}}}{{{\tau _1}}}{L_3}.$
式中, ${\tau _1}$, ${\tau _2}$${\tau _3}$是三个引入的常数, 其具体表达式为
${\tau _1} = {J_v}(kb){Y_v}(kc) - {J_v}(kb){Y_v}(kc), $
${\tau _2} = {B'_3}[{Y_v}(kb){s_{ - 1,v}}(kc) - {Y_v}(kc){s_{ - 1,v}}(kb)], $
${\tau _3} = {B'_3}[{J_v}(kc){s_{ - 1,v}}(kb) - {J_v}(kb){s_{ - 1,v}}(kc)] .$
压电复合陶瓷圆环内外表面间的电压可表示为
${V_{r0}} = - \int_b^c {{E_{r0}}} {\rm{d}}r = - \frac{{{\rm{j}}kb{X_b}}}{\omega }{v_b} - \frac{{{\rm j}kc{X_c}}}{\omega }{v_c} - \frac{{{L_3}P}}{{{\tau _1}}}.$
式中,
${X_b} = {B_3}\left\{ (v - 2){s_{ - 2,v - 1}}(kb) - \frac{{[{J_{v - 1}}(kb){Y_v}(kc) - {J_v}(kc){Y_{v - 1}}(kb)]{s_{ - 1,v}}(kb)}}{{{\tau _1}}} - \frac{{2{s_{ - 1,v}}(kc)}}{{\pi kb{\tau _1}}}\right\} + \frac{{{B_5}}}{{kb}},$
${X_c} = {B_3}\left\{ (v - 2){s_{ - 2,v - 1}}(kc) + \frac{{[{J_{v - 1}}(kc){Y_v}(kb) - {J_v}(kb){Y_{v - 1}}(kc)]{s_{ - 1,v}}(kc)}}{{{\tau _1}}} + \frac{{2{s_{ - 1,v}}(kb)}}{{\pi kc{\tau _1}}}\right\} + \frac{{{B_5}}}{{kc}},$
$\begin{split}P=\;&{\tau }_{1}\left\{\frac{{B}_{3}{B}^{\prime }_{3}({H}_{c}-{H}_{b})}{\rm{2}{v}^{2}(4-{v}^{2})}+{B}_{5}{B}^{\prime }_{3}{}_{\rm{1}}[{s}_{-1,v}(kc)-{s}_{-1,v}(kb)]+\left(\frac{{C}_{22}^{\rm{cp}}}{{\epsilon }_{33}^{\rm{cp}}{C}_{22}^{\rm{cp}}+{e}_{32}^{\rm{cp}}{e}_{32}^{\rm{cp}}}-\frac{{B}_{3}{B}^{\prime }_{3}}{{v}^{2}}\right)(\mathrm{ln}c-\mathrm{ln}b)\right\}\\&+{\tau }_{2}\{{B}_{3}[(v-2)({J}_{1c}-{J}_{1b})+{J}_{2b}-{J}_{2c}]+{B}_{5}[{J}_{v}(kc)-{J}_{v}(kb)]\\&+{\tau }_{3}\{{B}_{3}[(v-2)({Y}_{1c}-{Y}_{1b})+{Y}_{2b}-{Y}_{2c}]+{B}_{5}[{Y}_{v}(kc)-{Y}_{v}(kb)], \end{split}$
${J_{1r}} = kr{J_v}(kr){s_{ - 2,v - 1}}(kr),$
${J_{2r}} = kr{J_{v - 1}}(kr){s_{ - 1,v}}(kr), $
${Y_{1r}} = kr{Y_v}(kr){s_{ - 2,v - 1}}(kr), $
${Y_{2r}} = kr{Y_{v - 1}}(kr){s_{ - 1,v}}(kr), $
${H_r} = {(kr)^2}_2{F_3}\!\left(\! {[1,1],\left[2,2 + \frac{1}{2}v,2 - \frac{1}{2}v\right],} \right.\left. { \!-\! \frac{1}{4}{{(kr)}^2}} \!\right) .$
根据压电复合陶瓷径向力边界条件: ${F_b} \!= \!- {T_b}{S_b}$, ${F_c} = - {T_c}{S_c}$, 其中${S_b} = 2\pi bh$, ${S_c} = 2\pi ch$, 化简可得
$ \begin{split}{F}^{\prime }_{b}=\;&c{X}_{c}{F}_{b}=\bigg\{-\frac{{n}^{2}}{{\rm{j}}\omega {C}_{0}}\\&-{Z}_{b}\frac{{\rm{j}}{c}^{2}{X}_{c}^{2}[{J}_{v}(kc){Y}_{v-1}(kb)-{J}_{v-1}(kb){Y}_{v}(kc)]}{{\tau }_{1}}\\ &+\frac{4{\rm{j}}\pi h{c}^{2}{X}_{c}^{2}({B}_{2}-{B}_{4}v)}{\omega }\bigg\}{{v}^{\prime }}_{b}\\&+\left(-\frac{{n}^{2}}{{\rm{j}}\omega {C}_{0}}-{Z}_{b}\frac{2{\rm{j}}b{X}_{c}^{}{M}_{b}^{}}{\pi k{\tau }_{1}}\right){{v}^{\prime }}_{c}+n{V}_{r0},\\[-13pt] \end{split} $
$ \begin{split}{F}^{\prime}_{c}=\;&b{X}_{b}{F}_{c}=\bigg\{-\frac{{n}^{2}}{{\rm{j}}\omega {C}_{0}}\\&-{Z}_{c}\frac{{\rm{j}}{b}^{2}{X}_{b}^{2}[{J}_{v}(kb){Y}_{v-1}(kc)-{J}_{v-1}(kc){Y}_{v}(kb)]}{{\tau }_{1}}\\ &+\frac{4{\rm{j}}\pi h{b}^{2}{X}_{b}^{2}({B}_{4}v-{B}_{2})}{\omega }\bigg\}{v}^{\prime }_{c}\\&+\left(-\frac{{n}^{2}}{{\rm{j}}\omega {C}_{0}}-{Z}_{c}\frac{2{\rm{j}}c{X}_{c}^{}{M}_{b}^{}}{\pi k{\tau }_{1}}\right){v}^{\prime }_{b}+n{V}_{r0}.\\[-13pt] \end{split}$
式中,
$ \begin{split}&{Z_b} = \frac{{k{B_4}{S_b}}}{\omega } = \rho {V_3}{S_b},\;{Z_c} = \frac{{k{B_4}{S_c}}}{\omega } = \rho {V_3}{S_c},\;\\&{n_1} = c{X_c},\;{n_2} = b{X_b}, \end{split}$
$ \begin{split}&{v'_b} = \frac{{{v_b}}}{{{n_2}}},\;{v'_c} = \frac{{{v_c}}}{{{n_1}}},\;n = 2\pi h{\tau _1}kbc{X_b}{X_c}/P,\;\\&{C_0} = 2\pi h{\tau _1}/P, \end{split}$
$n$${C_0}$分别是机电转换系数和钳定电容. 流过压电陶瓷复合材料圆环的电流可以由下式得出:
$ \begin{split}{I_r} =\;& - \frac{{{\rm{d}}Q}}{{{\rm{d}}t}} = - {\rm{j}}\omega \int_{ - h/2}^{h/2} {\int_0^{2\pi } {{L_3}{\rm{d}}\theta dz} } = - {\rm{j}}\omega 2\pi h{L_3} \\=\;& - n({v'_b} + {v'_c}) + {\rm{j}}\omega {C_0}{V_{r0}}.\\[-10pt] \end{split}$
基于以上推导, 可以得到径向极化压电陶瓷复合材料圆环径向振动的六端机电等效电路如图4所示.
图 4 径向极化压电陶瓷复合圆环径向振动的机电等效电路
Figure4. Electromechanical equivalent circuit of radial vibration of a piezoelectric ceramic composite ring with radial polarization.

图4中, 各个部分阻抗的表达式如下:
$\begin{split} {Z}_{\rm{1}}=\;&-{Z}_{b}\frac{{\rm{j}}{c}^{2}{X}_{c}^{2}[{J}_{v}(kc){Y}_{v-1}(kb)-{J}_{v-1}(kb){Y}_{v}(kc)]}{{\tau }_{1}}\\&+\frac{{2}{\rm{j}}\pi h{c}^{2}{X}_{c}^{2}({B}_{2}-{B}_{4}v)}{\omega }+{Z}_{b}\frac{2{\rm{j}}c{X}_{c}^{}{X}_{b}^{}}{\pi k{\tau }_{1}},\\[-15pt] \end{split}$
$\begin{split} {Z}_{\rm{2}}=\;&-{Z}_{c}\frac{{\rm{j}}{b}^{2}{X}_{b}^{2}[{J}_{v}(kb){Y}_{v-1}(kc)-{J}_{v-1}(kc){Y}_{v}(kb)]}{{\tau }_{1}}\\&+\frac{{2}{\rm{j}}\pi h{b}^{2}{X}_{b}^{2}({B}_{4}v-{B}_{2})}{\omega }+{Z}_{c}\frac{2{\rm{j}}b{X}_{c}^{}{X}_{b}^{}}{\pi k{\tau }_{1}},\\[-15pt] \end{split}$
${Z_{\rm{3}}} = - {Z_b}\frac{{2{\rm{j}}cX_c^{}X_b^{}}}{{\pi k{\tau _1}}} = - {Z_c}\frac{{2{\rm{j}}bX_c^{}X_b^{}}}{{\pi k{\tau _1}}}.$

2
3.2.金属圆环径向振动的等效电路
-->金属圆环内外半径, 厚度分别为a, b, h, 其径向振动的四端机电等效电路如图5所示.
图 5 径向振动金属圆环的机电等效电路
Figure5. Electromechanical equivalent circuit of radial vibrating metal ring.

图5中, 各个阻抗的表达式如下[13]:
$\begin{split} {Z_{1{\rm{m}}}} =\;& {\rm{j}}\frac{{2{Z_{0a}}}}{{\pi {k_0}a[{J_1}\left( {{k_0}b} \right){Y_1}\left( {{k_0}a} \right) - {J_1}\left( {{k_0}a} \right){Y_1}\left( {{k_0}b} \right)]}} \\ & \times \frac{{{J_1}\left( {{k_0}b} \right){Y_0}\left( {{k_0}a} \right) - {J_0}\left( {{k_0}a} \right){Y_1}\left( {{k_0}b} \right) - {J_1}\left( {{k_0}a} \right){Y_0}\left( {{k_0}a} \right) + {J_0}\left( {{k_0}a} \right){Y_1}\left( {{k_0}a} \right)}}{{{J_1}\left( {{k_0}a} \right){Y_0}\left( {{k_0}a} \right) - {J_0}\left( {{k_0}a} \right){Y_1}\left( {{k_0}a} \right)}} \\ &- {\rm{j}}\frac{{2{Z_{0a}}(1 - {v_0})}}{{\pi {{({k_0}a)}^2}[{J_1}\left( {{k_0}a} \right){Y_0}\left( {{k_0}a} \right) - {J_0}\left( {{k_0}a} \right){Y_1}\left( {{k_0}a} \right)]}}, \end{split} $
$\begin{split} {Z_{2{\rm{m}}}} = \;&{\rm{j}}\frac{{2{Z_{0a}}}}{{\pi {k_0}a[{J_1}\left( {{k_0}b} \right){Y_1}\left( {{k_0}a} \right) - {J_1}\left( {{k_0}a} \right){Y_1}\left( {{k_0}b} \right)]}} \\ &\times \frac{{{J_1}\left( {{k_0}a} \right){Y_0}\left( {{k_0}b} \right) - {J_0}\left( {{k_0}b} \right){Y_1}\left( {{k_0}a} \right) - {J_1}\left( {{k_0}b} \right){Y_0}\left( {{k_0}b} \right) + {J_0}\left( {{k_0}b} \right){Y_1}\left( {{k_0}b} \right)}}{{{J_1}\left( {{k_0}b} \right){Y_0}\left( {{k_0}b} \right) - {J_0}\left( {{k_0}b} \right){Y_1}\left( {{k_0}b} \right)}} \\ &+ {\rm{j}}\frac{{2{Z_{0b}}(1 - {v_0})}}{{\pi {{({k_0}b)}^2}[{J_1}\left( {{k_0}b} \right){Y_0}\left( {{k_0}b} \right) - {J_0}\left( {{k_0}b} \right){Y_1}\left( {{k_0}b} \right)]}}, \end{split} $
$\begin{split} {Z_{{\rm{3m}}}} =\;& \frac{{ - {\rm{2}}{\rm{j}}{Z_{0a}}}}{{{k_0}a\pi \left[ {{J_1}\left( {{k_0}b} \right){Y_1}\left( {{k_0}a} \right) - {J_1}\left( {{k_0}a} \right){Y_1}\left( {{k_0}b} \right)} \right]}} = \frac{{ - {\rm{2}}{\rm{j}}{Z_{0b}}}}{{{k_0}b\pi \left[ {{J_1}\left( {{k_0}a} \right){Y_1}\left( {{k_0}b} \right) - {J_1}\left( {{k_0}b} \right){Y_1}\left( {{k_0}a} \right)} \right]}} \end{split}. $
式中,
$\begin{split}&{S_a} = 2\pi ah,\;{S_b} = 2\pi bh,\;{Z_{0a}} = {\rho _0}{v_r}{S_a},\;{Z_{0b}} = {\rho _0}{v_r}{S_b},\; {k_0} = \omega /{v_r},\;v_r^2 = \frac{{{E_0}}}{{{\rho _0}\left( {1 - v_0^2} \right)}},\\[-15pt]\end{split}$
${\rho _0}$, ${E_0}$${v_0}$分别是金属圆环的密度、杨氏模量和泊松比, ${k_0}$${v_{\rm{r}}}$分别是金属圆环中波数和径向振动的传播速度.
2
3.3.新型径向复合换能器的等效电路及频率方程
-->利用力及振速的连续性条件, 可以得出径向振动压电陶瓷复合材料圆环换能器的机电等效电路, 如图6所示.
图 6 径向振动压电陶瓷复合材料圆环换能器的机电等效电路
Figure6. Electromechanical equivalent circuit of radial vibration piezoelectric ceramic composite toroidal transducer.

图6${Z_{1{\rm{L}}}}$${Z_{2{\rm{L}}}}$分别代表新型换能器内外表面的径向负载阻抗, 当忽略换能器内外负载阻抗, 即${Z_{1{\rm{L}}}}$${Z_{2{\rm{L}}}}$都为0时, 由图6可以得到整个换能器的机械阻抗为[14]
${Z_{\rm{m}}} = {Z_3} + \frac{{{Z_1}({Z_2} + n_2^2{Z_{{\rm{om}}}})}}{{{Z_1} + {Z_2} + n_2^2{Z_{{\rm{om}}}}}}, $
其中压电复合陶瓷圆环机械阻抗为
${Z_{{\rm{om}}}} = {Z_{1{\rm{m}}}} + \frac{{{Z_{2{\rm{m}}}}{Z_{3{\rm{m}}}}}}{{{Z_{2{\rm{m}}}} + {Z_{3{\rm{m}}}}}}, $
换能器的输入电阻抗为
${Z_e} = \frac{V}{I} = \frac{{{n^2} - {\rm j}\omega {C_0}{Z_{\rm{m}}}}}{{{\omega ^2}C_0^2{Z_{\rm{m}}}}}, $
则换能器共振频率方程为
${n^2} - {\rm j} \omega {C_0}{Z_{\rm{m}}} = {\rm{0}}, $
反共振频率方程为
${\omega ^2}C_0^2{Z_{\rm{m}}} = {\rm{0}}.$

2
4.1.换能器振动性能与几何尺寸、两相体积占比的关系
-->通过求解频率方程(43)式和(44)式可以得出换能器的共振及反共振频率, 从而可以对换能器的共振频率与其几何尺寸[15]以及复合材料两相体积占比[16]之间的关系进行进一步的研究. 金属圆环的材料为不锈钢, 其材料参数为[4]
${\rho _0} = $7800 kg/m3, $E = 2.09$ 进行11 N/m2, ${v_0} = 0.28$.
压电陶瓷复合材料圆环的压电相为PZT-4, 其材料参数为
${\rho _c} = 7500$ kg/m3, $C_{11}^{\rm{c}} = $13.9 × 1010 N/m2, $C_{12}^{\rm{c}} = $7.78 × 1010 N/m2, $C_{13}^{\rm{c}} = $7.43 × 10 10 N/m2, $C_{33}^{\rm{c}} = $11.5 × 1010 N/m2, $e_{31}^{\rm{c}} = $–5.2 Nm/V, $e_{33}^{\rm{c}} = $15.1 Nm/V, $\varepsilon _{33}^{\rm{c}} = $5.62 C/m.
聚合物相为环氧树脂, 其材料参数为[17]:
${\rho _0} = $1050 kg/m3, $C_{11}^{\rm{p}} = $8 × 1010 N/m2, $C_{12}^{\rm{c}} = $4.4 × 1010 N/m2, $\varepsilon _{11}^{\rm{c}} = $0.037 × 10–9 C/m.
换能器的几何尺寸分别为:$a \!=\!$6 mm, $c \!=\!$26 mm, $h = $6 mm, $a \leqslant b \leqslant c$.
图7(a)图7(b)是换能器一阶共振频率和反共振频率与其半径比$\tau = \dfrac{{b - a}}{{c - a}}$之间的理论关系. 如图7(a)所示, 同一两相体积占比(${v_{\rm p}}/{v_{\rm c}} = 0.25$)情况下, 当半径比增加时, 新型复合换能器的一阶共振频率和反共振频率呈下降趋势. 用fr, fa分别表示换能器的共振频率和反共振频率, 换能器有效机电耦合系数可以表示为${k_{\rm eff}} = \sqrt {1 - {{\left( {{{{f_r}}}/{{{f_a}}}} \right)}^2}}$. 如图7(b)所示, 同一两相占比(${v_{\rm p}}/{v_{\rm c}} = 0.25$)情况下, 当半径比增加时, 新型复合换能器的有效机电耦合系数呈先增大后减小的趋势, 半径比在0.35左右, 有效机电耦合系数取到最大值. 另外, 如图7(c)图7(d)所示, 同一半径比($\tau = 0.6$)情况下, 当两相占比增加时, 新型复合换能器的一阶共振频率和反共振频率呈上升趋势, 而有效机电耦合系数呈下降趋势. 因此在实际换能器设计中, 应充分考虑尺寸与两相占比对换能器振动性能的影响, 为了得到较高的机电转换效率, 换能器半径比设计应尽量接近于0.35, 聚合物相占比越高虽然会导致机电转换效率降低, 但同时也能带来更好的声匹配能力, 因此在换能器设计中选择较低聚合物相占比即可.
图 7 (a) 换能器的一阶共振频率和反共振频率与几何尺寸之间的关系; (b) 换能器的有效机电耦合系数与几何尺寸之间的关系; (c) 换能器的一阶共振频率和反共振频率与两相占比之间的关系; (d) 换能器的有效机电耦合系数与两相占比之间的关系
Figure7. (a) Relationship between the first-order resonance frequency and anti-resonance frequency of the transducer and the geometrical size; (b) relationship between the effective electromechanical coupling coefficient and geometric dimensions of the transducer; (c) relationship between the first-order resonance frequency and anti-resonance frequency of the transducer and the proportion of the two phases; (d) relationship between the effective electromechanical coupling coefficient of the transducer and the proportion of the two phases.

2
4.2.换能器振动性能的仿真模拟
-->为验证换能器的解析理论, 利用仿真软件(COMSOL Multiphysics 5.4)对换能器的振动模态及共振频率进行了数值模拟. 结果如表1所示, fr, fa分别表示换能器一阶径向振动共振频率、反共振频率的理论计算结果, fr1, fa1分别表示换能器一阶径向振动数值模拟结果, 误差${A_1} = \left| {{f_{\rm{r}}} - {f_{{\rm{r}}1}}} \right|/{f_{{\rm{r}}1}}$, ${A_2} = \left| {{f_{\rm{a}}} - {f_{{\rm{a}}1}}} \right|/{f_{{\rm{a}}1}}$. 而误差主要来源于在解析理论中, 假设换能器的厚度远小于其径向尺寸, 理想情况应该是无限小, 而在数值模拟过程中换能器尺寸则是有限的, 此外在复合材料等效参数推导中并没有考虑材料形状对结构参数的影响, 因此会产生一定误差. 由表1可以看出, 数值模拟得到的换能器在不同半径比或不同两相占比的情况下(空气中), 一阶径向振动的共振频率、反共振频率与根据解析法得出的结果吻合较好, 从而验证了新型换能器理论设计的正确性, 图8表1不同半径比或不同两相占比的换能器共振时的振动模态图像.
a/mmb/mmc/mmvp/vch/mmfr/Hzfa/Hzfr1/Hzfa1/HzA1/%A2/%
616240.429651780.754216.253198.053941.02.663.08
816240.429647428.949485.347552.047741.02.593.65
616240.250650047.052761.252675.052941.04.990.34
816240.250645949.448257.847901.048299.04.070.08


表1新型径向复合材料换能器共振频率的理论及数值模拟结果(一阶径向振动)
Table1.Theoretical and numerical simulation results of the resonance frequency of the new radial composite transducer(first-order radial vibration).

图 8 (a) 换能器共振频率为53198 Hz时的振动模态; (b) 换能器共振频率为47552 Hz时的振动模态; (c) 换能器共振频率为52675 Hz时的振动模态; (d) 换能器共振频率为47901 Hz时的振动模态
Figure8. (a) The vibration mode of the transducer when the resonance frequency is 53198 Hz; (b) the vibration mode of the transducer when the resonance frequency is 47552 Hz; (c) the vibration mode of the transducer when the resonance frequency is 52675 Hz; (d) the vibration mode of the transducer when the resonance frequency is 47901 Hz.

2
4.3.相比传统径向换能器的性能提升
-->为了验证新型换能器相比传统径向换能器的性能提升, 利用仿真软件(COMSOL Multiphysics 5.4)对换能器在水下的辐射声场进行了仿真研究[18], 假设换能器处在无限大水域内工作. 在53到65 kHz范围内, 同一尺寸($a = $6 mm, $b = $16 mm, $c = $24 mm, $h = $6 mm, ${v_{\rm p}}/{v_{\rm c}} = 0.429$)的新型换能器的发射电压响应相比传统径向换能器幅值更大, 其原因在于采用压电复合材料以后, 换能器的声阻抗匹配得到了改善. 另外, 以中心频率对应发射电压响应左右下降3 dB取换能器工作带宽[19,20], 由数值模拟结果可知该尺寸下新型径向换能器在水中振动的中心频率在53082 Hz, 振动模态如图9(a)所示, 传统径向换能器的中心频率在55356 Hz, 振动模态如图9(c)所示, 则由图10可见, 新型换能器工作带宽为6000 Hz左右, 即为图10中两蓝线间隔, 而传统径向换能器工作带宽为3000 Hz左右, 即为图10中两绿线间隔, 两者相差接近一倍, 因此新型换能器相比传统径向换能器性能上有很大的提升, 是一种针对传统换能器性能提升的改进方案. 另外在图10中可以看到新型换能器和传统换能器在68和73 kHz左右出现峰值, 因此对利用仿真软件对两种换能器在对应峰值频率附近寻找特征频率, 其振动模态分别如图9(b)图9(d)所示, 由图可见此时两种换能器的振动模态非常复杂, 存在径向振动和弯曲振动的耦合振动, 且振动位移较大, 因此在水中发射电压响应出现了一个峰值, 但由于振动模态复杂, 并非本文所研究模态.
图 9 (a) 换能器共振频率为53082 Hz时的振动模态; (b) 换能器共振频率为68070 Hz时的振动模态; (c) 换能器共振频率为55356 Hz时的振动模态; (d) 换能器共振频率为73747 Hz时的振动模态
Figure9. (a) The vibration mode of the transducer when the resonance frequency is 53082 Hz; (b) the vibration mode of the transducer when the resonance frequency is 68070 Hz; (c) the vibration mode of the transducer when the resonance frequency is 55356 Hz; (d) the vibration mode of the transducer when the resonance frequency is 73747 Hz.

图 10 新型复合材料径向换能器与传统纯陶瓷径向换能器的发射电压响应曲线
Figure10. The emission voltage response curves of the new radial transducer and the traditional radial transducer.

本文研究了一种新型径向振动压电陶瓷复合材料圆环换能器, 推导了2-2型压电复合材料的等效参数, 得到了径向极化压电陶瓷复合材料圆环换能器的机电等效电路及其频率方程. 利用仿真软件对换能器的径向振动性能进行了数值模拟, 结果表明数值模拟得到的换能器一阶径向振动的共振频率、反共振频率与用解析法得出的结果吻合较好. 主要结论主要有以下几点:
1) 推导了2-2型压电复合材料的等效参数, 并将其应用于径向换能器理论设计中;
2) 同一两相占比(${v_{\rm{p}}}/{v_{\rm{c}}} = 0.25$)情况下, 当半径比增加时, 新型复合换能器的一阶共振频率和反共振频率呈下降趋势. 同一两相占比(${v_{\rm{p}}}/{v_{\rm{c}}} = 0.25$)情况下, 当半径比增加时, 复合换能器的有效机电耦合系数呈先增大后减小的趋势, 半径比在0.35左右, 有效机电耦合系数取到最大值. 另外, 同一半径比($\tau = 0.6$)情况下, 当两相占比增加时, 复合换能器的一阶共振频率和反共振频率呈上升趋势, 而有效机电耦合系数呈下降趋势;
3) 在53到65 kHz范围内, 新型复合材料径向换能器相比传统纯陶瓷径向换能器, 发射电压响应幅值更大, 工作带宽提高接近一倍, 声匹配更佳.
相关话题/机电 金属 工作 电压 软件

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 二维过渡金属硫属化合物相变方法的研究进展
    摘要:近些年来,二维过渡金属硫属化合物因其具有独特的原子结构和能带特征而备受关注.过渡金属硫属化合物不仅存在稳定相,也存在一些亚稳相,这些亚稳相因其独特的电学和光学特性逐渐成为研究的热点.为了获得这些亚稳相,一种能够使过渡金属硫属化合物在两种晶相之间可控相变的方法必不可少.本文首先对过渡金属硫属化合 ...
    本站小编 Free考研考试 2021-12-29
  • 多狄拉克锥的二维CrPSe<sub>3</sub>的半金属铁磁性与电子结构<span style="color:#ff0000;">(已撤稿)&
    摘要:采用第一性原理方法和紧束缚模型研究了二维层状结构过渡金属硫族磷化物CrPSe3的磁学特性和电子结构.二维CrPSe3晶格呈现六角蜂窝状结构,声子谱证明其具备较好的动力学稳定性.二维CrPSe3的电子结构呈现出自旋无能隙磁性半金属特性.蒙特卡罗方法模拟出二维CrPSe3的居里温度为224K,分子 ...
    本站小编 Free考研考试 2021-12-29
  • 金属光栅异常透射增强黑磷烯法拉第旋转的理论研究
    摘要:黑磷是一种具有直接带隙的二维材料,其较宽的带隙填补了石墨烯和二维过渡金属硫化物之间的带隙空白,其特殊的褶皱状晶体结构导致了其独特的面内各向异性,使其具有了独特的力电磁响应特性.本文基于单层黑磷设计了一种金光栅/黑磷/硅的混合等离子体结构的磁光器件.通过金属光栅诱导异常透射显著增强透射率,同时通 ...
    本站小编 Free考研考试 2021-12-29
  • 金属锂在固态电池中的沉积机理、策略及表征
    摘要:全固态金属锂电池的能量密度有望达到现有商业化锂离子电池的2—5倍,且有可能从本质上解决现有液态电解质锂离子电池的安全性问题.要想实现全固态金属锂电池这一颠覆性技术,克服金属锂/固态电解质界面存在的副反应严重、界面接触差、锂枝晶生长等一系列挑战至关重要.本文探讨了金属锂在有机、无机固态电解质中的 ...
    本站小编 Free考研考试 2021-12-29
  • 实用化条件下金属锂负极失效的研究
    摘要:金属锂因其高理论比容量和低电极电势,被视为下一代高比能二次电池理想的负极材料之一.然而,其表面不稳定的固液界面膜及不均匀锂沉积等问题严重限制着其实际应用.目前,金属锂负极的研究大多采用温和的实验条件,这对于理解负极表界面的物理化学性质和揭示锂沉积规律等基础研究具有重要意义.但是,超薄锂负极(< ...
    本站小编 Free考研考试 2021-12-29
  • 全固态金属锂电池负极界面问题及解决策略
    摘要:全固态金属锂电池有望提高当前商用锂离子电池的安全性及能量密度,被广泛认为是下一代电池的重要研发方向.其中的负极-电解质界面与电池性能紧密相连.本文将该界面存在的问题划分为静态及动态两方面,静态问题包括化学不稳定及物理接触差,体现在电池循环前的巨大阻抗,动态问题包括枝晶生长及孔洞形成,体现在电池 ...
    本站小编 Free考研考试 2021-12-29
  • 三重简并拓扑半金属MoP中超快圆偏振光产生和调控光生热电流
    摘要:超快激光可以用来产生和调控拓扑量子材料中的拓扑相变和自旋极化电流,这些光诱导产生的新奇物性源于材料中受到体系对称性保护的线性色散能带结构的简并节点.作为一种同时具有两重和三重简并节点的拓扑半金属,磷化钼(MoP)是一类非常独特的拓扑半金属体系.本文初步探索了三重简并拓扑半金属晶体MoP中产生自 ...
    本站小编 Free考研考试 2021-12-29
  • 难熔金属钒熔化行为的局域原子结构模拟与分析
    摘要:采用经典分子动力学(moleculardynamics,MD)方法,模拟了16000个钒原子在5种不同熔化速率(γ1=1×1011K/s,γ2=1×1012K/s,γ3=1×1013K/s,γ4=1×1014K/s与γ5=1×1015K/s)下原子结构的熔化行为.结果表明:不同熔化速率对难熔金 ...
    本站小编 Free考研考试 2021-12-29
  • 非贵金属表面增强拉曼散射基底的研究进展
    摘要:表面增强拉曼散射(surface-enhancedRamanscattering,SERS)在分析检测领域中具有重要地位,然而随着其不断发展,贵金属SERS基底在实际应用中受到限制.基于C,Ti,Zn,Cu,Mo,W等非贵金属纳米材料的SERS基底相比于贵金属基底具有更优异的经济性、稳定性、选 ...
    本站小编 Free考研考试 2021-12-29
  • 过渡金属(Cr, Mn, Fe, Co)掺杂对TiO<sub>2</sub>磁性影响的第一性原理研究
    摘要:关于过渡金属掺杂TiO2是否会产生室温铁磁性及其磁性的来源存在争议,为了解决此问题,本文采用基于密度泛函理论下的GGA+U方法对体系Ti0.875X0.125O2(X=Cr,Mn,Fe,Co)的磁学性质及光学性质进行了第一性原理的研究.首先计算了铁磁和反铁磁的基态能量,比较后推测出铁磁态为它们 ...
    本站小编 Free考研考试 2021-12-29