-
塑料类高分子材料甲基丙烯酸甲酯-亚克力(PMMA)圆柱中亚音速Rayleigh波低频隧穿共振可引起反向散射增强, 在低频标准散射体设计等领域具有重要应用价值. 提出一种微弱形变的规则波纹表面结构, 可实现水中PMMA圆柱反向散射低频共振频率的无源调控. 利用微扰法推导了水中微弱形变规则波纹圆柱反向散射低频共振频率偏移的近似解, 讨论了波纹微扰系数、周期对规则波纹圆柱共振频率偏移的影响规律. 基于Rayleigh波相位匹配方法分析了低频共振频率偏移的机理. 研究表明: 微弱形变规则波纹圆柱中亚音速Rayleigh波沿微弱形变波纹表面传播, 与光滑圆柱体相比, 传播路径的改变引起Rayleigh波传播相位变化, 导致了Rayleigh波低频共振频率发生偏移. 最后开展了微弱形变规则波纹圆柱体声散射特性水池实验, 获取了其反向散射共振频率, 明显观察到了规则波纹圆柱共振频率偏移现象, 与理论预报结果吻合较好.
-
关键词:
- 规则波纹 /
- 圆柱 /
- 亚音速Rayleigh波 /
- 共振散射
Backscattering enhancement associated with subsonic Rayleigh wave resonance of a polymethlmethacrylate (PMMA) cylinder is observed at low frequencies in water, which suggests that the PMMA cylinders may have essential applications in the low-frequency standard scatterer design. A slightly deformed surface structure with regular corrugation is presented to manipulate the low-frequency backscattering resonance frequency of PMMA cylinder passively. Using the perturbation method, the approximate resonance frequency shift is derived for an infinite slightly deformed cylinder with regular corrugations. Based on the phase matching of Rayleigh waves, the mechanism of low-frequency resonance frequency shift is revealed. Studies show that compared with a bare cylinder, a small boundary deformation can cause the propagation path of Rayleigh waves to change, namely, the Rayleigh waves propagate along the slightly deformed surface with regular corrugations. The modified propagation path can cause the propagation phase to change, which brings about the low-frequency backscattering resonance frequency shift of a PMMA cylinder. Furthermore, how the resonance frequency shifts with the deformation coefficient and period of the corrugation is discussed in detail. The backscattering resonance frequency of the regular corrugated cylinder shifts to low frequency with the increase of ξ under the condition of the deformation coefficient ξ > 0, but shifts to high frequency with the decrease of ξ at ξ < 0, and the resonance frequency shift increases with the increase of absolute value of deformation coefficient. When corrugation period m < 10, the phase variation with m is too small, so the corrugation period has little effect on the resonance frequency shift. Finally, acoustic scattering experiment of the regular corrugated cylinder is conducted in the tank. The resonance frequency shift is obviously observed in the experiment, which is in good agreement with the theoretical prediction. Hence, the characteristics of backscattering enhancement associated with subsonic Rayleigh wave and the resonance frequency shifts make the PMMA deformed cylinder have potential applications such as in standard scatter design and identification using “AcoustiCode”.-
Keywords:
- regular corrugation /
- cylinder /
- subsonic Rayleigh wave /
- resonance scattering
1. 引 言
在光学领域, Ge等[1,2]利用微扰理论研究了边界微弱形变光学微腔的共振特性, 研究表明微腔边界形变导致不同共振模态之间耦合, 使得其共振频率发生偏移, 并在很大程度改变腔内射线动力学特性和输出方向, 且光波的输出方向对边界变形极其灵敏, 在光束控制和高分辨率检测方面具有重要应用前景. 在声学领域, Fawcett[3]用数值方法计算了外径不变内径微弱形变水中球壳声散射频域特征, 结果显示边界的微弱形变可以使得由于球壳中传播的a0波与刚性背景耦合形成的中频吻合共振[4]得到抑制, 但未对抑制中频吻合共振的机理进行研究; 提出了一种时域Kirchhoff/衍射混合算法[5], 此算法只适用于微弱形变刚性目标高频声散射的计算. Antonio等[6]提出一种数值曲面积分方法, 可以用来计算不规则表面或随机粗糙表面几何目标的声散射问题, 但此方法只适用于刚性目标. 本文基于微扰法分析了微弱形变弹性圆柱体低频声散射特征, 提出了一种利用微弱形变规则波纹表面结构调控塑料类高分子材料甲基丙烯酸甲酯-亚克力(PMMA)圆柱反向散射低频共振频率的方法.
塑料类高分子聚合物材料PMMA是3D打印领域的一类主要材料, 其密度与水接近, 剪切波速度小于水中声速[7], 近年来在声学领域逐渐受到关注. Hefner和Marston[8]通过理论和实验研究了水中PMMA球体低频反向散射增强现象, 研究表明由于亚音速Rayleigh波共振使得PMMA球体目标低频反向散射大幅度增强, 并指出这类PMMA目标在水中低频标准散射体设计等方面具有重要应用价值. PMMA目标低频共振频率的人工调控, 是设计满足不同频率需求的声散射标准体, 是水中被动声学标记体的基础[9,10]. 本文提出了一种应用于水中PMMA圆柱体低频共振频率调控的规则波纹表面, 第2节基于微扰法推导了无限长微弱形变规则波纹圆柱体声散射共振频率偏移的近似解, 分析了其低频共振频率偏移的机理, 即规则波纹表面导致圆柱中亚音速Rayleigh波传播路径改变, 从而引起传播相位变化, 使得低频共振频率发生偏移. 第3节讨论波纹微扰系数、周期等波纹参数对水中PMMA规则波纹圆柱低频共振散射频率偏移的影响规律. 第4节进行了光滑PMMA圆柱体和规则波纹PMMA圆柱体(相同波纹周期、不同波纹微扰系数)声散射特性水池试验, 获取了其反向散射共振特性, 观察到共振频率偏移特征, 同时验证了第2节规则波纹圆柱共振频率偏移近似公式的有效性和正确性.
2. 无限长微弱变形表面圆柱体声散射
2.1 微扰理论
无限长微弱形变规则波纹表面圆柱体示意如图1所示.
若圆柱表面存在微小变形[3], 即
R(a,φ)=a+εf(φ), (1) 其中, a为未变形圆柱半径, ε为波纹微扰系数,
f(φ)=sin2(mφ),m=0,1,2,3,⋯ , m是[0,π] 的波纹周期. 则波纹圆柱圆周长度S为S(ε,m)=∫2π0[a+εsin2(mφ)]2+[2mεsin(mφ)cos(mφ)]2dφ. (2) 由(2)式可以看出, 微弱形变规则波纹圆柱圆周长度与未变形圆柱体半径a、波纹微扰系数ε和波纹周期m有关.
假设
|εf(φ)|≪a , 入射平面波沿着x轴正向垂直入射到无限长波纹圆柱, 设散射声场具有如下形式[11,12]:pscat(r,φ)=∞∑n=0εnin[bnH(1)n(kr)cos(nφ)+ε∑p≠nαbpH(1)p(kr)cos(pφ)], (3) 这里忽略时间因子e–iωt, 其中k = ω/c是入射平面波波数, ω为角频率, c为水中声速, r为接收点到目标中心的距离,
φ 为入射平面波与x轴正向的夹角,εn={1,n=0,2,else, bn是散射系数,H(1)n(kr) 为第一类汉克尔函数. 则微弱形变规则波纹圆柱外部总声场为pt(r,φ)=∞∑n=0εnin{[Jn(kr)+bnH(1)n(kr)]cos(nφ)+ε∑p≠nαbpH(1)p(kr)cos(pφ)}. (4) 微弱形变规则波纹圆柱内部声场:
Φ(r,φ)=∞∑n=0εnin[cnJn(kdr)cos(nφ)+ε∑p≠nαcpJp(kdr)cos(pφ)],Ψ(r,φ)=∞∑n=0εnin[enJn(ksr)sin(nφ)+ε∑p≠nαepJp(ksr)sin(pφ)], (5) 其中, Jn是n阶Bessel函数, αbp, αcp和αep是待定的耦合系数, kd和ks分别为纵波波数和横波波数, p为耦合阶次.
根据边界条件:
τrr(a+εf(φ),φ)=−pt(a+εf(φ),φ), ur(2)(a+εf(φ),φ)=ur(1)(a+εf(φ),φ), τrφ(a+εf(φ),φ)=0, 其中,
τrr 是结构中的法向应力,u(1)r 和u(2)r 分别为结构和水中的法向位移,τrφ 是结构中的切向应力. 基于微扰理论将(6)式按照ε的一阶泰勒级数展开, 则有τrr(a,φ)+pt(a,φ)=−εf(φ)[∂τrr∂r+∂pt∂r](a,φ), u(2)r(a,φ)−u(1)r(a,φ)=−εf(φ)[∂u(2)r∂r−∂u(1)r∂r](a,φ), τrφ(a,φ)=−εf(φ)∂τrφ∂r(a,φ), 由(7)式, 得
[Dd11Dd12Dd13000Dd21Dd22Dd230000Dd32Dd33000Dd41Dd42Dd43Dd44Dd45Dd46Dd51Dd52Dd53Dd54Dd55Dd560Dd62Dd630Dd65Dd66]×(bncnenαbpαcpαep)=(A1∗A2∗0A∗4A∗50), (8) 根据克莱姆法则得到:
bn=BdnDdn,αbp=AbpnDdn, (9) 其中
Bdn ,Abpn 和Ddn 为6 × 6的行列式, 行列式具体元素见附录A.无限长规则波纹圆柱远场反向散射形态函数[13]表示为
f(x)=2√iπx∞∑n=0εn[Bdn(x)Ddn(x)(−1)n+ε∑p≠nAbpn(x)Ddn(x)(−1)p]. (10) 令
x=kR 为无限长微弱形变规则波纹圆柱特征方程Ddn(x)=0 的根. 当ε=0 时, (10)式为无限长光滑圆柱远场反向散射形态函数, 此时通过文献[14]的方法求解特征方程得到无限长光滑圆柱特征值的实部, 记为x0 . 根据共振散射理论[15], 无限长圆柱特征值的实部对应共振峰的频率, 则无限长圆柱共振频率f0=x0c2πa , 对应的相速度[16]:cph=2πf0a/n,n=0,1,2,3,⋅⋅⋅. (11) 当
ε≠0 时, 无限长微弱形变规则波纹圆柱共振频率记为f. 将x在x0 附近按照ε的一阶展开为x=x0+εx1. (12) 将(6)式在
x0 附近按照ε的一阶展开, 具体推导过程见附录A, 经过一系列数学化简得x1=−x0Ann, (13) 其中
Ann=εnπ∫π0f(φ)acos(nφ)cos(nφ)dφ. 则微弱形变规则波纹圆柱特征值偏移量为
Δx=x−x0=−εx0Ann. (14) 对应的微弱形变规则波纹圆柱共振频率偏移量可表示为
Δf=f−f0=−εf0Ann. (15) 无限长微弱形变规则波纹圆柱共振频率近似为
f=f0−εf0Ann. (16) 定义微弱形变规则波纹圆柱归一化波纹微扰系数ξ = ε/a, 利用(10)式计算半径a = 0.05 m材料为PMMA的无限长光滑圆柱和波纹参数m = 6, ξ = –3%的无限长微弱形变规则波纹圆柱反向散射形态函数幅频特性, 所用材料参数见表1. 计算频率范围为1—18 kHz, 频率间隔为5 Hz, 所对应无因次频率ka的范围为0.20—4.19, ka较小, 文中所研究频段为低频.
表 1 计算所用材料参数Table 1. Material parameters used in the calculations.材料 密度/(kg·m–3) 纵波波速/(m·s–1) 剪切波波速/(m·s–1) PMMA 1190 2690 1340 水 1000 1500 — 图2(a)对比了无限长光滑圆柱(黑色实线)和微弱形变规则波纹圆柱(红色点划线)反向散射形态函数幅频特性, 研究发现边界微弱变形使得PMMA圆柱共振频率发生了偏移. 图2(b)是图2(a)蓝色虚线框中共振峰的局部放大图, 无限长光滑圆柱对应的共振峰频率f0 = 6745 Hz, 微弱形变规则波纹圆柱的共振峰频率f1 = 6845 Hz与利用(16)式计算的微弱形规则波纹圆柱共振频率f2 = 6846.2 Hz基本一致. 因此, 可以用(16)式近似估计无限长微弱形变规则波纹圆柱的共振频率.
2.2 微弱形变圆柱低频共振机理
根据(11)式计算得到无限长PMMA圆柱Rayleigh波相速度频散曲线, 如图3所示. 从图3可以看出, Rayleigh波相速度(黑色实线)随着频率的增大逐渐趋于常数
cR≈1060 m/s, 小于水中声速, 为亚音速波. 文献[8,17]利用射线近似理论详细描述了PMMA球体中亚音速Rayleigh波对低频反向散射增强的贡献及传播路径. 对于平面波垂直入射无限长PMMA圆柱, Rayleigh波沿圆柱表面周向传播, 如图4(a)所示. 亚音速Rayleigh波在临界角θl=π/2 的B1点被激发, 在B1, B2点(箭头处)耦合沿圆柱表面周向传播, D1, D2点处耦合从D1点出射进入散射声场. 在焦散线上, Rayleigh波的投影切向速度等于水中的声速. 声场从渐消到传播的转变发生在半径b=ac/cR 附近. 在a<r<b 区域为渐消波区域, Rayleigh波通过耦合隧穿此渐消波区域再辐射形成反向散射声场. 因为入射平面波对称激发, 所以环绕波成对出现, 它们的绕行方向相反, 从D1到B1方向传播的波这里没有显示.当圆柱表面具有微弱形变时, Rayleigh波传播过程与光滑圆柱相同, 但传播路径发生改变, 如图4(b)所示. Rayleigh波沿规则波纹表面传播, 由微弱形变表面所确定的焦散半径
bd=[a+εsin2(mφ)]c/cR 发生变化. Rayleigh波传播路径的改变引起传播相位变化, 由环绕波分析理论[18]可知, Rayleigh波沿规则波纹圆柱表面周向传播的圆周数v=0,1,⋯ , 在接收点形成程差为S的等间隔波列, 这些波列的相位差为Δφ(f,ξ,m)=ωS(ξ,m)cR=2πfS(ξ,m)cR. (17) 若这些波列的相位差是
2π 的M倍(M为整数), 同相叠加形成共振, 由此得到共振条件Δφ=2Mπ,M=1,2,3,⋯, (18) cR=S(ξ,m)Mf. (19) 当
ε=0 时,cR=2πaf0/M , 无限长PMMA圆柱中Rayleigh波频散曲线与曲线2πaf0/M (图3中蓝色点划线)的交点所对应的横坐标为f′0 , 无限长PMMA圆柱的共振频率为f0 , 见表2. 当ε≠0 时, 利用(16)式求得的规则波纹圆柱共振频率为fξ , 无限长PMMA圆柱中Rayleigh波频散曲线与曲线Sf/M (图3中红色虚线)的交点所对应的横坐标为f′ξ , 这里m = 6,ξ=5% . 定义由共振条件确定的共振频率相对误差为δf=|f−f′|/f×100% , 由图5可知对于光滑圆柱和规则波纹圆柱, 由共振条件确定的共振频率相对误差均在0.3%以下. 因此, 利用Rayleigh波共振条件可以近似估计无限长PMMA圆柱反向散射低频共振频率. 同时, 这也证明了规则波纹PMMA圆柱中Rayleigh波速度不变, 但Rayleigh波传播路径发生改变引起传播相位的变化, 使得PMMA圆柱低频共振频率发生偏移, 从而实现了对PMMA圆柱反向散射低频共振频率的调控.表 2 无限长PMMA圆柱和规则波纹圆柱共振频率Table 2. Resonance frequency of PMMA bare and corrugated cylinderl f0/Hz f′0/Hz fξ /Hz f′ξ/Hz 1 6419 6420 6259 6240 2 9953 9955 9704 9700 3 13383 13385 13048 13050 4 16790 16790 16370 16375 3. 波纹参数对低频共振频率的影响
以上研究表明, 边界微弱形变可以使得PMMA圆柱反向散射低频共振频率发生偏移. 由(16)式很直观地看出无限长规则波纹圆柱反向散射低频共振频率与波纹微扰系数、波纹周期有关, 基于微扰法分析边界变形结构参数如波纹微扰系数、波纹周期对无限长PMMA圆柱反向散射低频共振频率的影响. 图6(a)为规则波纹圆柱反向散射波纹微扰系数-频率谱, 颜色表示形态函数幅值, 横轴是频率, 单位为kHz, 纵轴是归一化波纹微扰系数ξ范围–10%—10%, 间隔1%. 图中黑色虚线为根据(16)式获取的各个共振峰频率随波纹微扰系数变化曲线, 当归一化波纹微扰系数ξ > 0时, 波纹在
[a,a+ε] 区域内起伏, 规则波纹圆柱反向散射共振频率随着ξ的增大向低频偏移, 且频率偏移量随着ξ的增大而增大; 当归一化波纹微扰系数ξ < 0时, 波纹在[a+ε,a] 区域内起伏, 规则波纹圆柱反向散射共振频率随着ξ的减小向高频偏移, 且频率偏移量随着ξ的减小而增大. 因此可以通过改变波纹微扰系数对水中PMMA圆柱反向散射低频共振频率进行无源调控. 图6(b)为规则波纹圆柱反向散射频率-波纹周期谱, 颜色表示形态函数幅值, 横轴是频率, 单位kHz, 纵轴是波纹周期m范围2—6, 这里ξ = –3%. 图中黑色虚线为根据(16)式获取的共振峰频率随波纹周期变化曲线, 可见当波纹周期比较小时, 波纹周期对无限长规则波纹圆柱反向散射共振频率偏移影响较小. 不同波纹周期所对应相位变化量为ΔΦm=Δφ(f0,ξ,m)−Δφ(f0,ξ,1) , 图7给出共振频率f0 = 6419 Hz, ξ = –3%, m取值0—20所对应的相位变化量ΔΦm , 可见当m < 10时, m变化引起的相位变化较小, 不足以对共振频率偏移产生重要作用.4. 实 验
4.1 实验布置
针对PMMA规则波纹圆柱开展声散射特性水池实验, 此次实验是在中国科学院声学研究所北海站消声水池完成. 3个实验模型(一个光滑圆柱, 两个规则波纹圆柱, 波纹周期m = 6, 归一化波纹微扰系数分别为ξ = –3%和ξ = –10%)均为亚克力材料, 密度ρ = 1150 kg/m3, 半径a = 0.05 m, 长度L = 0.5 m. 为了方便吊放, 模型两端各打4个孔, 实验模型垂直置于水池中绕中轴线转动. 发射换能器采用刚性连接吊放位置距离目标3.65 m, 水听器型号为B&K8103, 吊放位置距离目标2.11 m. 发射换能器、水听器和目标三者处于同一深度, 距离水面3.87 m, 实验模型及具体布放如图8和图9所示. 为排除市电干扰, 测量过程中所有设备全程处于UPS直流供电状态.
4.2 实验数据分析
发射信号为线性调频信号, 频率为5—15 kHz, 脉宽0.8 ms, 发射周期为500 ms. 将水听器接收到的时域信号进行宽带信号处理[17], 截取直达波(即发射阵直接到达水听器的信号)时域信号
pi(t) 和回波时域信号ps(t) 进行傅里叶变换得到其频谱Pi(f) 和Ps(f) , 再将回波信号频谱除以直达波信号频谱并对距离进行修正获得3个实验模型反向散射形态函数幅频特性曲线, 如图10所示. 实验中在所研究频段明显观察到圆柱反向散射形态函数3个共振峰(黑色线), 微弱形变规则波纹圆柱反向散射共振峰明显的向高频偏移, 且ξ = –10%的变形圆柱(蓝色虚线)共振峰频率偏移量大于ξ = –3%的变形圆柱(红色点划线)共振峰偏移量, 与第3节理论分析结果一致.从图10中3个实验模型的反向散射形态函数中提取共振峰频率, 光滑圆柱共振峰频率记为f0, 规则波纹圆柱共振峰频率表示为
f′ξ , 利用(16)式获取的共振峰频率为fξ , 详细数值见表3. 这里定义实验获取的微弱形变规则波纹圆柱共振峰频率相对误差为表 3 实验模型共振峰频率Table 3. Resonance frequencies of experimental objects.l f0/Hz fξ1/Hz f′ξ1 /Hz fξ2/Hz f′ξ2/Hz 1 6649 6748.7 6852 6981.4 6966 2 8643 8772.6 9098 9075.1 9294 3 12410 12596.0 12996 13031.0 13192 δξ=|f′ξl−fξl|fξl×100%. (20) 为了更为直观地观察各个共振频率, 以图11(a), (b)柱状图形式呈现, 分别对应归一化波纹微扰系数ξ = –3%和ξ = –10%规则波纹的共振频率. 实验获取微弱形变规则波纹圆柱共振峰频率与近似(16)式所估算的结果吻合较好, 相对误差均在5%以内(见图11(c)). 因此, 可以用(16)式快速估计微弱形变规则波纹圆柱的共振频率, 从而实现对PMMA圆柱反向散射低频共振频率的调控.
图12为实验获取的微弱形变规则波纹圆柱的频率-角度谱, 横轴表示入射角度, 纵轴表示频率, 单位为kHz, 颜色表示形态函数幅值, 黑色虚线表示共振峰频率. 从图12可以看出: 微弱形变规则波纹圆柱低频共振频率在各个角度都发生了偏移, 且在各个方位角偏移量基本相同. 主要的原因是归一化微扰系数ξ = –3%和ξ = –10%所对应的微扰系数分别为ε = 0.0015 m和ε = 0.005 m, 在所研究频段内Rayleigh波波长λ > 0.07 m (所对应频率为f = 15 kHz), 此时Rayleigh波波长远大于波纹微扰系数, 对散射起主要作用的仍然是低频Rayleigh波共振, 波纹多镜面反射的影响尚未清晰显示出来[19], 因此微弱形变规则波纹PMMA圆柱低频共振频率偏移量与入射方向无关.
5. 结 论
本文提出了一种微弱形变的规则波纹表面结构, 可实现PMMA圆柱体反向散射低频共振频率的无源调控. 利用微扰法推导了无限长微弱形变规则波纹圆柱共振频率的近似解, 研究发现由于边界微弱形变Rayleigh波传播路径发生改变, 沿着规则波纹表面传播, 传播路径的改变引起传播相位的变化, 从而导致共振峰频率发生偏移. 讨论了波纹微扰系数、周期对规则波纹圆柱共振频率偏移的影响规律: 当归一化波纹微扰系数ξ > 0时, 规则波纹圆柱反向散射共振频率随着ξ的增大向低频偏移, 当ξ < 0时, 规则波纹圆柱反向散射共振频率随着ξ的减小而向高频偏移, 共振峰频率偏移量随着波纹微扰系数绝对值的增大而增大; 当波纹周期m < 10时, 波纹周期对共振频率偏移影响较小. 最后开展了规则波纹圆柱体声散射特性水池实验, 获取了不同波纹微扰系数规则波纹圆柱体的反向散射共振频率, 明显观察到了规则波纹圆柱共振峰频率偏移现象, 与理论预报结果吻合较好. PMMA圆柱低频亚音速Rayleigh波隧穿共振引起反向散射幅度增强和规则波纹结构引起低频共振频率偏移的特征, 使得规则波纹PMMA圆柱在标准散射体设计、利用“声条形码”进行身份识别方面具有潜在的应用价值.
此外, 文中基于微扰法推导的用于计算无限长规则波纹圆柱声散射的解, 只适用于微弱形变, 即波纹微扰系数和波纹周期较小的情况, 此时忽略了由于法向变化引起的微扰, 只考虑了函数值引起的微扰, 是一种近似解. 若要对本问题或相关问题进行深入探讨, 获取更为精确的解, 可参考文献[20]中规则形腔中简正模式的微扰近似方法.
附录A
对于正横方向入射无限长规则波纹圆柱体, 行列式
Ddn ,Bdn 和Abpn 分别为Ddn=|Dd11Dd12Dd13000Dd21Dd22Dd230000Dd32Dd33000Dd41Dd42Dd43Dd44Dd45Dd46Dd51Dd52Dd53Dd54Dd55Dd560Dd62Dd630Dd65Dd66|, Bdn=|A∗1Dd12Dd13000A∗2Dd22Dd230000Dd32Dd33000A∗4Dd42Dd43Dd44Dd45Dd46A∗5Dd52Dd53Dd54Dd55Dd560Dd62Dd630Dd65Dd66|, Abpn=|Dd11Dd12Dd13A∗100Dd21Dd22Dd23A∗2000Dd32Dd33000Dd41Dd42Dd43A∗4Dd45Dd46Dd51Dd52Dd53A∗5Dd55Dd560Dd62Dd630Dd65Dd66|, A∗1=−[A1+ε(A′1−2A1)Ann],A∗2=−[A2+ε(A′2−2A2)Ann],A∗4=−(A′1−2A1)Anp,A∗5=−(A′2−2A2)Anp, Dd11=d11+ε(d′11−2d11)Ann,Dd12=d12+ε(d′12−2d12)Ann,Dd13=d13+ε(d′13−2d13)Ann,Dd21=d21+ε(d′21−d21)Ann,Dd22=d22+ε(d′22−d22)Ann,Dd23=d23+ε(d′23−d23)Ann,Dd32=d32+ε(d′32−2d32)A1nn,Dd33=d33+ε(d′33−2d33)A1nn,Dd41=(d′11−2d11)Anp,Dd42=(d′12−2d12)Anp,Dd43=(d′13−2d13)Anp,Dd44=d11p+ε(d′11p−2d11p)App,Dd45=d12p+ε(d′12p−2d12p)App,Dd46=d13p+ε(d′13p−2d13p)App,Dd51=(d′21−d21)Anp,Dd52=(d′22−d22)Anp,Dd53=(d′23−d23)Anp,Dd54=d21p+ε(d′21p−d21p)App,Dd55=d22p+ε(d′22p−d22p)App,Dd55=d22p+ε(d′22p−d22p)App,Dd62=(d′32−2d32)A1np,Dd63=(d′33−2d33)A1np,Dd65=d32p+ε(d′32p−2d32p)A1pp,Dd66=d33p+ε(d′33p−2d33p)A1pp, A1=ρ1ρ21ρ1ω2k2sa2Jn(ka), A2=−kaρ1ω2J′n(ka),A′1=ρ1ρ21ρ1ω2[2k2sa2Jn(ka)+kak2sa2J′n(ka)],A′2=−kaρ1ω2J′n(ka)−k2a2ρω2J′′n(ka), d11=ρ1ρ21ρ1ω2k2sa2H(1)n(ka),d12=(2n2−k2sa2)Jn(kda)−2kdaJ′n(kda),d13=2n[ksaJ′n(ksa)−Jn(ksa)],d′11=ρ1ρ21ρ1ω2[2k2sa2H(1)n(ka)+kak2sa2H(1)n′(ka)],d′12=−2k2sa2Jn(kda)+kda(2n2−k2sa2−2)J′n(kda)−2k2da2J″n(kda),d′13=2nk2sa2J′′n(ksa), d21=−kaρ1ω2H(1)n′(ka),d22=kdaJ′n(kda), d23=nJn(ksa),d′21=−kaρ1ω2H(1)′n(ka)−k2a2ρ1ω2H(1)″n(ka),d′22=kdaJ′n(kda)+k2da2J″n(kda), d′23=nksaJ′n(ksa), d32=2n[Jn(kda)−kdaJ′n(kda)],d33=2ksaJ′n(ksa)+(k2sa2−2n2)Jn(ksa),d′32=−2nk2da2J″n(kda),d′33=2k2sa2J″n(ksa)+2k2sa2Jn(ksa)+ksa[(ksa)2−2n2+2]J′n(ksa), A1nn=εnπ∫π0f(φ)asin(nφ)sin(nφ)dφ,A1pp=εpπ∫π0f(φ)asin(pφ)sin(pφ)dφ,App=εpπ∫π0f(φ)acos(pφ)cos(pφ)dφ, 其中,
d11p ,d12p ,d13p ,d21p ,d22p ,d23p ,d32p ,d33p 和d′11p ,d′12p ,d′13p ,d′21p ,d′22p ,d′23p ,d′32p ,d′33p 分别等于d11 ,d12 ,d13 ,d21 ,d22 ,d23 ,d32 ,d33 和d′11 ,d′12 ,d′13 ,d′21 ,d′22 ,d′23 ,d′32 ,d′33 中将阶次n替换成p的值.J′n 和J″n 分别代表贝塞尔函数的一阶和二阶导数.将
kdR=c1cdx ,ksR=c1csx 代入(6a)式, 化简得τrrn(x)+ptn(x)=μk2({ρ1ρ21ρ1ω2c21c2sH(1)n(x)Bn[(2n2x2−c21c2s)Jn(c1cdx)−2c1cd1xJ′n(c1cdx)]Cn+2n[c1cs1xJ′n(c1csx)−1x2Jn(c1csx)]En+ρ1ρ21ρ1ω2c21c2sJn(x)}cos(nφ)+ε∑p≠n{ρ1ρ21ρ1ω2c21c2sH(1)p(x)αbp[(2p2x2−c21c2s)Jp(c1cdx)−2c1cd1xJ′p(c1cdx)]αcp+2p[c1cs1xJ′p(c1csx)−1x2Jp(c1csx)]αep}cos(pφ)), 将(A1)式中x在x0处按照ε的一阶展开得
τrrn(x0)+ptn(x0)=μk2[εx1{ρ1ρ21ρ1ω2c21c2sH(1)′n(x0)Bn+[c1cd(2n2x02−c21c2s+2x02)J′n(c1cdx0)−4n2x30Jn(c1cdx0)−2c12cd21x0J″n(c1cdx0)]Cn+2n[c12cs21x0J″n(c1csx0)−c1cs2x02J′n(c1csx0)+2x03Jn(c1csx0)]En+ρ1ρ21ρ1ω2c21c2sJ′n(x0)}cos(nφ)+ε(∑p≠n{ρ1ρ21ρ1ω2c21c2sH(1)p(x0)αbp[(2p2x20−c21c2s)Jp(c1cdx0)−2c1cd1x0J′p(c1cdx0)]αcp+2p[c1cs1x0J′p(c1csx0)−1x20Jp(c1csx0)]αep})cos(pφ)], ∂τrrn∂r+∂ptn∂r=μk3[(ρ1ρ21ρ1ω2c21c2sH(1)′n(x)Bn+[c1cd(2n21x2−c21c2s+21x2)J′n(c1cdx)−2c21c2d1xJ″n(c1cdx)−4n21x3Jn(c1cdx)]Cn+{2nc21c2s1xJ″n(c1csx)−4n[c1cs1x2J′n(c1csx)−1x3Jn(c1csx)]}En+ρ1ρ21ρ1ω2c21c2sJ′n(x))cos(nφ)+ε∑p≠n(ρ1ρ21ρ1ω2c21c2sH(1)′p(x)αbp+{c1cd(2p21x2−c21c2s+21x2)J′p(c1cdx)−2c21c2d1xJ″p(c1cdx)−4p21x3Jp(c1cdx)}αcp+{2pc21c2s1xJ″p(c1csx)−4p[c1cs1x2J′p(c1csx)−1x3Jp(c1csx)]}αep)cos(pφ)], τrrn(x0)+ptn(x0)=−εf(φ)[∂τrrn∂r(x0)+∂ptn∂r(x0)]. [1] Ge L, Song Q H, Redding B, Eberspacher A, Wiersig J, Cao H 2013 Phys. Rev. A. 88 043801
Google Scholar
[2] Ge L, Song Q H, Redding B, Cao H 2013 H Phys. Rev. A. 87 023833.1
Google Scholar
[3] Fawcett J A 2016 IEEE J. Ocean. Eng. 41 682
Google Scholar
[4] Zhang L G, Sun N H, Marston P L 1992 J. Acoust. Soc. Am. 91 1862
Google Scholar
[5] Fawcett J A 2001 J. Acoust. Soc. Am. 109 1312
Google Scholar
[6] Guel-Tapia J A, Villa-Villa F, Mendoza-Suarez A, Perez-Aguilar H 2016 Archives of Acoustics 41 461
Google Scholar
[7] Hartmann B, Jarzynski B 1974 J. Acoust. Soc. Am. 56 1469
Google Scholar
[8] Hefner B T, Marston P L 2000 J. Acoust. Soc. Am. 107 1930
Google Scholar
[9] Satish A, Trivett D, Sabra K G 2020 J. Acoust. Soc. Am. 147 EL517
Google Scholar
[10] Srivastava P, Nichols B, Sabra K G 2017 J. Acoust. Soc. Am. 142 EL573
Google Scholar
[11] Doolittle R D, Überall H, Uginčius P 1968 J. Acoust. Soc. Am. 43 1
Google Scholar
[12] Dubertrand R, Bogomolny E, Djellali N, Lebental M, Schmit C 2008 Phys. Rev. A. 77 013804.1
Google Scholar
[13] Donald B, Gaunaurd G C 1983 J. Acoust. Soc. Am. 73 1448
Google Scholar
[14] 龙云亮, 文希理, 谢处方 1994 数值计算与计算机应用 2 88
Google Scholar
You Y L, Wen X L, Xie C F 1994 Journal on Numerical Methods and Computer Applications 2 88
Google Scholar
[15] 汤渭霖 1995 声学学报 06 456
Google Scholar
Tang W L 1995 ACTA. ACUSTICA. 06 456
Google Scholar
[16] Marston P L, Sun N H 1995 J. Acoust. Soc. Am. 97 777
Google Scholar
[17] 周彦玲, 范军, 王斌 2019 物理学报 68 214301
Google Scholar
Zhou Y L, Fan J, Wang B 2019 Acta Phys. Sin. 68 214301
Google Scholar
[18] 汤渭霖, 范军, 马忠诚 著 2018 水中目标声散射 (北京: 科学出版社) 第99—104页
Tang W L, Fan J, Ma Z C 2018 Acoustic Scattering of Underwater Targets (Beijing: Science Press) pp99–104 (in Chinese)
[19] 彭茜蕤, 周彦玲, 范军 2018 声学技术 37 528
Peng X R, Zhou Y L, Fan J 2018 Technical Acoustics 37 528
[20] 程建春 著 2019 声学原理 (北京: 科学出版社) 第615—619页
Cheng J C 2019 Acoustical Principle (Beijing: Science Press) pp615–619 (in Chinese)
-
表 1 计算所用材料参数
Table 1. Material parameters used in the calculations.
材料 密度/(kg·m–3) 纵波波速/(m·s–1) 剪切波波速/(m·s–1) PMMA 1190 2690 1340 水 1000 1500 — 表 2 无限长PMMA圆柱和规则波纹圆柱共振频率
Table 2. Resonance frequency of PMMA bare and corrugated cylinder
l f0/Hz f′0/Hz fξ /Hz f′ξ/Hz 1 6419 6420 6259 6240 2 9953 9955 9704 9700 3 13383 13385 13048 13050 4 16790 16790 16370 16375 表 3 实验模型共振峰频率
Table 3. Resonance frequencies of experimental objects.
l f0/Hz fξ1/Hz f′ξ1 /Hz fξ2/Hz f′ξ2/Hz 1 6649 6748.7 6852 6981.4 6966 2 8643 8772.6 9098 9075.1 9294 3 12410 12596.0 12996 13031.0 13192 -
[1] Ge L, Song Q H, Redding B, Eberspacher A, Wiersig J, Cao H 2013 Phys. Rev. A. 88 043801
Google Scholar
[2] Ge L, Song Q H, Redding B, Cao H 2013 H Phys. Rev. A. 87 023833.1
Google Scholar
[3] Fawcett J A 2016 IEEE J. Ocean. Eng. 41 682
Google Scholar
[4] Zhang L G, Sun N H, Marston P L 1992 J. Acoust. Soc. Am. 91 1862
Google Scholar
[5] Fawcett J A 2001 J. Acoust. Soc. Am. 109 1312
Google Scholar
[6] Guel-Tapia J A, Villa-Villa F, Mendoza-Suarez A, Perez-Aguilar H 2016 Archives of Acoustics 41 461
Google Scholar
[7] Hartmann B, Jarzynski B 1974 J. Acoust. Soc. Am. 56 1469
Google Scholar
[8] Hefner B T, Marston P L 2000 J. Acoust. Soc. Am. 107 1930
Google Scholar
[9] Satish A, Trivett D, Sabra K G 2020 J. Acoust. Soc. Am. 147 EL517
Google Scholar
[10] Srivastava P, Nichols B, Sabra K G 2017 J. Acoust. Soc. Am. 142 EL573
Google Scholar
[11] Doolittle R D, Überall H, Uginčius P 1968 J. Acoust. Soc. Am. 43 1
Google Scholar
[12] Dubertrand R, Bogomolny E, Djellali N, Lebental M, Schmit C 2008 Phys. Rev. A. 77 013804.1
Google Scholar
[13] Donald B, Gaunaurd G C 1983 J. Acoust. Soc. Am. 73 1448
Google Scholar
[14] 龙云亮, 文希理, 谢处方 1994 数值计算与计算机应用 2 88
Google Scholar
You Y L, Wen X L, Xie C F 1994 Journal on Numerical Methods and Computer Applications 2 88
Google Scholar
[15] 汤渭霖 1995 声学学报 06 456
Google Scholar
Tang W L 1995 ACTA. ACUSTICA. 06 456
Google Scholar
[16] Marston P L, Sun N H 1995 J. Acoust. Soc. Am. 97 777
Google Scholar
[17] 周彦玲, 范军, 王斌 2019 物理学报 68 214301
Google Scholar
Zhou Y L, Fan J, Wang B 2019 Acta Phys. Sin. 68 214301
Google Scholar
[18] 汤渭霖, 范军, 马忠诚 著 2018 水中目标声散射 (北京: 科学出版社) 第99—104页
Tang W L, Fan J, Ma Z C 2018 Acoustic Scattering of Underwater Targets (Beijing: Science Press) pp99–104 (in Chinese)
[19] 彭茜蕤, 周彦玲, 范军 2018 声学技术 37 528
Peng X R, Zhou Y L, Fan J 2018 Technical Acoustics 37 528
[20] 程建春 著 2019 声学原理 (北京: 科学出版社) 第615—619页
Cheng J C 2019 Acoustical Principle (Beijing: Science Press) pp615–619 (in Chinese)
计量
- 文章访问数: 5531
- PDF下载量: 61