针对这一问题，中科院声学所超声学实验室博士生臧雨宸与其导师林伟军研究员等人利用部分波级数展开法和镜像原理，结合球函数的加法公式，成功推导了阻抗边界附近粘弹性球壳在Bessel涡旋声场中的声辐射力矩表达式，并详细分析了球壳与边界间距离、边界声压反射系数、球壳的相对厚度与声束半锥角等因素对声辐射力矩的影响。
相关成果2021年6月在线发表于国际声学期刊 Journal of Sound and Vibration 。
图1 Bessel涡旋声束入射到阻抗边界附近的粘弹性球壳示意图（图/中科院声学所）
计算结果表明，随着边界声压反射系数的增大，声辐射力矩的峰值也随之增大，但峰值位置不变。事实上，这些峰值对应着球壳本身不同阶的本征模式。此外，边界的引入使在球壳与边界之间产生了准驻波场，声辐射力矩随球壳与边界距离的变化呈现周期变化的规律，其变化的周期与声波波长有关。
图2 一阶Bessel涡旋声束对不同声压反射系数下粘弹性PE球壳的声辐射力矩函数（图/中科院声学所）
球壳相对厚度的改变会影响球壳本身的本征模式，进而影响最终的声辐射力矩大小。结果显示，随着球壳厚度不断减小(内外径之比不断增大)，声辐射力矩的峰值会显著增大，且峰值的位置发生移动。但当球壳非常薄(内外径之比为0.99)时，声辐射力的峰值又会随之减小。因此，在实际的声操控中，需要设计合适的球壳的厚度从而获取需要的声辐射力矩。此外，声束的半锥角范围也需要合理选择，当半锥角为零时，声束将退化为平面波，从而无法产生声辐射力矩。
图3 一阶Bessel涡旋声束对不同厚度的粘弹性PE球壳的声辐射力矩函数（图/中科院声学所）
本研究得到国家自然科学基金(11604361和11904384)和中科院青年创新促进会(2019024)的资助。
关键词：
声辐射力矩；Bessel涡旋声束；粘弹性球壳；阻抗边界；声操控
参考文献：
ZANG Yuchen, LIN Weijun, ZHENG Yinfei, SU Chang, WU Pengfei, CHANG Qin, Acoustic radiation torque of a Bessel vortex wave on a viscoelastic spherical shell nearby an impedance boundary. Journal of Sound and Vibration. 2021. DOI: 10.1016/j.jsv.2021.116261.
论文链接：
https://doi.org/10.1016/j.jsv.2021.116261
