matlab刚性球散射声场编程,基于小球反射法的超声脉冲波声场特性研究

超声换能器在超声检测、超声医疗和超声制造等系统中有着广泛的应用。由于各种超声技术都是利用换能器产生的振动能量与物质相互作用,从而产生各种物理、化学和生物等效应[1],所以各系统超声作用的效果和声场的性质密切相关,因此掌握换能器声场分布的信息对衡量换能器的性能至关重要。目前,国内外测量换能器声场的方法主要有辐射力法、水听器法、光纤检测法、光学检测法[2]和小球反射法[3]。虽然相对于前4种方法基于小球反射法的测量系统测量精度较低,但其结构简单,成本低廉,方便使用,且能正确反映出声场声压分布趋势,因此能满足一般工业应用的需要。本文推导了小球散射超声换能器脉冲波声场的指向性函数计算公式,并由此建立了小球反射的超声换能器脉冲波声场声压分布理论模型,计算了由不同尺寸的小球反射的声场声压分布。并通过实验验证了理论模型的正确性。1小球散射声场指向性分析计算1.1小球散射单频连续波声场的指向性分析计算超声直入射式换能器的圆形晶片可看作理想的活塞声源[4-7],在其辐射声场的远场区域内,换能器辐射的超声波可近似看作平面波[8]。如图1所示,位于换能器A的远场区域中的半径a的刚性小球B散射由换能器激发的平面波,引入球坐标系r,取小球的球心为球坐标的原点,X轴为极轴,整个问题关于X轴对称,与无关。在球坐标系下,平面入射波声压为pI(x,t)=p0ei(t-kx)=p0eite-ikrcos(1)式中I为入射波;p0为换能器表面的声压;为角频率;t为时间;r为离小球圆心的距离;k为波数;x为换能器表面到小球圆心的距离。图1小球散射示意图因小球表面的振动分布与无关,所以小球散射的声场是轴对称的,散射波是各阶轴对称的球面波的叠加,可简化为pS(r,,t)=l=0BlPl(cos)hl(2)(kr)eit(2)式中Bl为常数,由边界条件决定;Pl(cos)为勒让德函数;h(l2)(kr)为第二类l阶球汉克尔函数;S为散射波。为了使入射平面波与小球散射波在形式上可比较,将入射平面波用球函数叠加,式(1)可表示为pI(r,,t)=p0eitl=0(-i)l(2l+1).Pl(cos)jl(kr)(3)式中jl(kr)为l阶球贝塞尔函数。在刚性小球的球面上(r=a),满足垂直振速为零的边界条件,于是i0r(pI+pS)|r=a=0(4)将式(2)、(3)代入式(4)中可得Bl=-(-j)l(2l+1)dd(ka)[jl(ka)]dd(ka)[h2l(ka)]p0(5)根据指向性函数的定义,由式(2)、(3)和(5)可得,散射声场的指向性函数[1]D()=pp0S=k1al=0ejl+21(-j)l.(2l+1)dd[[hjl2l((kkaa))]]//dd((kkaa)).Pl(cos)(6)在实际计算中,第90阶的球贝塞尔函数的导数已很接近于零,所以在MATLAB中编写程序对指向性函数进行计算时,只要计算到第100阶就足够了,在ka=1、5、10和50时分别计算小球散射单频稳态声场的指向性函数,得到最大值归一化后的指向性分布如图2所示,当声频率较低、小球较小时,小球面向声源一侧的散射较均匀。当声频率较高、小球较大时,开始出现花瓣状散射,且频率越高花瓣越多,同时,小球背向声源一侧有较强的散射波。图2小球散射单频稳态声场指向性1.2小球散射脉冲声场的指向性分析计算在实际超声检测中广泛使用的是脉冲波,所以测量换能器脉冲声场的特性才有实际的意义。换能器脉冲声场特性和单频稳态声场特性不同,相应地小球对脉冲声场的散射声场特性与对单频稳态声场的散射声场特性也有差别。根据声场线性叠加原理,利用傅