点声源入射下水下弹性球壳声散射分析.pdf
《点声源入射下水下弹性球壳声散射分析.pdf》由会员分享,可在线阅读,更多相关《点声源入射下水下弹性球壳声散射分析.pdf(10页珍藏版)》请在咨信网上搜索。
1、点声源入射下水下弹性球壳声散射分析郑金焱,陈美霞,董文凯(华中科技大学船舶与海洋工程学院,湖北武汉,430074)摘要:为了研究水下目标的声散射问题,在经典 Rayleigh 简正级数解析解的基础上,结合格林函数在球坐标系中的表达式,计算了点声源球面波与平面波的折算因子;基于平面波入射下弹性球壳的散射声场,完成了点声源入射下的弹性球壳散射声场求解,并进行了有限元验证。以此为基础,计算并分析了典型球状目标的目标强度;从目标强度的角度出发,对点声源与球心的距离进行量化,提出球面波可近似为平面波处理的条件。研究结果表明:当点声源在球体表面附近时,散射体的目标强度曲线峰值相对于平面波入射下会发生偏移;
2、此外在无因次频率 ka50 时,球面波可近似为平面波入射,点声源和球心的临界距离与球壳半径、球壳厚度有关。研究结果可为水下目标的探测、识别及水声测量实验提供理论参考。关键词:点声源;弹性球壳;声散射特性;目标强度中图分类号:TJ630.1;U666.7文献标识码:A文章编号:2096-3920(2023)05-0725-10DOI:10.11993/j.issn.2096-3920.2022-0053AcousticScatteringAnalysisofUnderwaterElasticSphericalShellunderPointSoundSourceIncidenceZHENG Jin
3、yan,CHEN Meixia,DONG Wenkai(SchoolofNavalArchitectureandOceanEngineering,HuazhongUniversityofScienceandTechnology,Wuhan430074,China)Abstract:Inordertostudytheacousticscatteringofunderwatertargets,theconversionfactorofplanewaveandsphericalwaveofpointsoundsourcewascalculatedonthebasisoftheanalyticalso
4、lutionoftheclassicalRayleighnormalseries,combinedwiththeexpressionofGreensfunctioninthesphericalcoordinatesystem.Basedonthescatteredsoundfieldoftheelasticsphericalshellundertheincidenceofplanewave,thescatteredsoundfieldoftheelasticsphericalshellundertheincidenceofpointsoundsourcewassolved,andfinitee
5、lementverificationwascarriedout.Onthisbasis,thetargetstrengthoftypicalsphericaltargetswerecalculatedandanalyzed.Fromtheperspectiveoftargetstrength,thedistancefromthecenterofthesphereandthepointsoundsourcewasquantified,andtheconditionthatthesphericalwavewasapproximatelyequivalenttotheplanewavewasprop
6、osed.Theresultsshowthatthepeakvalueofthetargetstrengthcurveofthespherewillshiftrelativetotheincidenceofplanewavewhenthepointsoundsourceisnearthesphere.Inaddition,thesphericalwavecanbeapproximatedtotheplanewaveincidencewiththedimensionlessfrequencyka50,andthedistancebetweenthepointsoundsourceandthesp
7、hericalcenterisrelatedtotheradiusandthicknessofthesphericalshell.Theresultscanprovideatheoreticalreferenceforthedetection,identification,andunderwateracousticmeasurementexperimentsofunderwatertargets.Keywords:pointsoundsource;elasticsphericalshell;acousticscatteringcharacteristic;targetstrength收稿日期:
8、2022-08-26;修回日期:2022-10-28.基金项目:国家自然科学基金项目资助(52071152).作者简介:郑金焱(2000-),男,在读硕士,主要研究方向为振动与噪声控制.第31卷第5期水下无人系统学报Vol.31 No.52023年10月JOURNALOFUNMANNEDUNDERSEASYSTEMSOct.2023引用格式 郑金焱,陈美霞,董文凯.点声源入射下水下弹性球壳声散射分析 J.水下无人系统学报,2023,31(5):725-734.水下无人系统学报sxwrxtxb.xml-7250引言水下球形目标的声散射特性一直是水声技术的研究热点:降低结构的散射声场能有效提高潜艇
9、的隐身性能;研究目标的回波强度,可以提升主动声呐的远程探测距离。因此,研究目标的声散射特性具有重要的理论价值和广泛的工程应用背景。国内外围绕该问题在声散射理论、声散射计算方法以及实验模型等方面进行了广泛而细致的研究。声波在水下传播时遇到障碍物时会产生散射。Strutt 等1提出了当球体尺寸远小于声波波长时散射声场的近似解。Morse2推导了任意频率的平面波入射下刚性球的散射理论解。Li3和刘国利4等研究了平面波斜入射情况下无限长弹性圆柱壳体的共振散射问题。对于散射问题,积分方程法也可得到严格解,用分离变量法求解 Helmholtz 方程可以计算已知表面振速物体的辐射声场,并得到精确解5。Faw
10、cett6给出不同流体填充下的弹性壳体散射声场解。朱凤芹7采用分离变量法推导出平面波入射下内部填充流体的弹性壳体散射声场表达式,利用数值方法分析了球体、球壳和柱壳的声散射频率特性并进行实验验证。近年来,随着声散射研究理论和计算机技术的快速发展和不断深入,庄悦等8基于 Boit 理论,建立了流体多孔球的散射模型,并采用数值的方法进行计算,简单讨论了不同沙子孔隙率和半径下沙粒的背向散射情况。Godin 等9利用无穷级数展开的形式,研究了低频球面声波对球形障碍物的散射问题。梁国龙等10从工程实际出发,用数值法分析了 4 种情况下弹性球壳散射声场对矢量水听器测向影响的规律。李运思等11运用边界元法开展
11、了刚性椭球颗粒的声散射数值研究,并采取一种声散射函数近似的计算方法研究了液体颗粒的散射。徐慧等12应用声学仿真软件,对球壳阵列下水下目标的隔声性能进行了相关研究。在以往的研究中,散射体一般都为球形,Juan等13将研究方向进行拓展,用解析方法给出了长椭球和扁椭球在平面波作用下的精确解,对散射系数使用算法进行简化求取,并进行了数值讨论。Li 等14基于声散射理论,对具有多层介质覆盖的弹性球壳声散射进行研究计算,发现在低频时隐身效果显著提高。张健等15利用无穷级数中函数的近似形式,推导出了水下低频球面声波入射下阻抗球的散射声压渐进解的公式。何嘉华等16借助镜像源方法以及等效源原理,研究了阻抗半空间
12、边界下二维圆柱的声散射问题。周彦玲等17基于经典 Rayleigh 简正级数解的方法和有限元数值方法对水下球体目标进行计算,研究了球体散射的相位特性。吴杰等18针对敷设空腔覆盖层的弹性球壳进行仿真研究,得出敷设空腔覆盖层能够有效降低球体的目标强度,发现在一定范围内增大空腔覆盖层的物理参数可以有效增强结构的隐身性能。王明升等19通过声散射理论,采用分波序列的方法分析了液体球及弹性球在 Bessel 波束作用下声辐射力的变化规律,并给出不同介质组合及壳层厚度对粒子的声辐射力的影响。臧雨宸等20通过理论推导给出柱面波入射下多层球的轴向声辐射力表达式,并进行了相应的数值计算。前人所做的工作证明了研究规
13、则目标的声散射问题具有工程意义,并且平面波入射下各种规则目标的声散射研究已经相当充分。然而在实际水声环境中,声源的情况复杂(如各种不规则形状的活塞声源、阵列声源),将其简化处理为平面波并不准确。点声源辐射的是球面波,在近场与平面波的特性有较大的区别,进而导致目标的散射声场特性差异较大。求取点声源入射下的散射声场,对复杂声源入射下的声散射研究工作具有借鉴意义。文中首先从理论上推导出平面波与点声源球面波的关系,然后基于平面波入射下弹性球壳的散射声场,运用叠加法推导出在球面波激发下球体散射声压的精确解。同时,利用有限元软件验证理论结果的准确性。然后计算了球面波入射下刚性球、弹性球和弹性球壳的目标强度
14、。以目标强度为出发点,对点声源与球心的距离进行量化,提出球面波近似为平面波入射的条件,并进行了相关讨论,最后以球体和球壳为例对球面波入射时的散射特性进行了分析。文中研究扩充了复杂声源入射下规则目标声散射特性的内容。2023年10月水下无人系统学报第31卷726JournalofUnmannedUnderseaSystemssxwrxtxb.xml-1理论推导1.1折算因子(r,)图 1 所示为球坐标系,其中坐标原点为球心 O,点声源距离球心矢径为 r1,声压测量点距球心矢径为 r,点声源位于 z 轴上,平面波从 z 轴方向入射。平面波 P1xyz(r,)Orr1图1球面入射几何模型Fig.1G
15、eometrymodelofsphericalincidence截取球体的 xOz 剖面进行分析,由三维无限大空间的 Green 函数在球坐标系中的展开式,利用加法定理可以将点声源发散的球面波展开i=eik|rr1|rr1|=ikn=0(2n+1)Pncos()jn(kr)h(1)n(kr1),0rr1jn(kr1)h(1)n(kr),0r1r(1)式中:Pn为 n 阶 Legnedre 函数;jn为 n 阶球 Bessel函数;hn(1)为第 1 类球 Hankel 函数;k 为流体自由波数。为简化分析,取单位点声源距离球心 1m 处声压为 1Pa,并略去时间因子 exp(it)。由于所研究
16、问题为轴对称,所以式(1)中的球谐函数与方位角 无关。取 rr1的入射波势函数i=ikn=0(2n+1)Pn(cos)jn(kr)h(1)n(kr1)(2)Pn(cos)=(1)nPn(cos)由于,所以入射波势函数化简为i=i2n+1kn=0(2n+1)jn(kr)h(1)n(kr1)Pn(cos)=n=0injn(kr)Pn(cos)(3)in=i2n+1k(2n+1)h(1)n(kr1)式中,。将入射平面波按球面波分解,势函数表达为i=eikrcos=n=0in(2n+1)jn(kr)Pn(cos)=n=0injn(kr)Pn(cos)(4)in=in(2n+1)其中:。将两者入射波势函
17、数逐项相除,就可将球面波和平面波联系起来,得到折算因子Qn=inin=in+1kh(1)n(kr1),n=0,1,2(5)1.2平面波入射下弹性球壳声散射对于如图 2 所示的弹性球壳,整个空间被分为 3 个区域:区域 1 是外部流体空间,其密度、声速分别为 1、c1;区域 2 是弹性球壳,其外半径为 a,内半径为 b,球壳密度为 0,纵波波速为 CL,剪切波速为 CT;区域 3 是壳内流体介质,其密度和声速分别为 2、c2。入射平面波 P1Ps1c12c2ab区域 1区域 2区域 3(r,)zxO图2弹性球壳及入射平面波Fig.2Elasticsphericalshellandincident
18、planewave考虑一单位振幅的平面波沿 z 轴入射,将声波按球面波分解可表示为P1=eik1rcos=n=0in(2n+1)jn(k1r)Pn(cos)(6)式中,k1=/c1为壳体外部区域波数。同理将散射声场 Ps和内部声场 P2写成球面波展开的形式Ps=n=0in(2n+1)h(1)n(k1r)Pn(cos)Bn(7)P2=n=0in(2n+1)jn(k2r)Pn(cos)Gn(8)式中:k2=/c2为壳体内部区域波数;Bn和 Gn为待定系数。2023年10月郑金焱,等:点声源入射下水下弹性球壳声散射分析第5期水下无人系统学报sxwrxtxb.xml-727在弹性介质中引入标量势和向量
19、势,考虑到球体的轴对称性,标量势和向量势 可表示为=n=0Cnjn(kLr)+Dnnn(kLr)Pn(cos)(9)=n=0Enjn(kTr)+Fnnn(kTr)Pn(cos)(10)式中:Cn、Dn、En和Fn为待定系数;nn为球Neumann函数。根据球坐标中应力与应变的关系,应用壳体表面上应力与振速连续的边界条件,可以列出确定上述各待定系数的矩阵方程DX=A(11)式中:D 为 66 阶矩阵;X 为待求系数矩阵。D 和A 中相应元素在文献7 给出,根据Cramer法则,Bn的解为Bn=fndn(12)式中:fn是用 A 代替 D 中的第 1 列所得到的行列式;dn为方程的系数行列式。若球
20、壳内部为真空,则去掉行列式 fn、dn中第 5 行第 6 列即可。若球壳内半径为 0,弹性球壳退化为弹性球,此时弹性介质中的势函数去掉球 Neumann 函数,矩阵 D 以及矩阵 A 变为 33 阶的行列式,相应的系数 Bn也发生变化,矩阵元素在这里不列出。1.3球面波入射下弹性球壳声散射(r,)取球坐标,坐标原点在球心 O,球壳相关参数与平面波入射时一致,点声源距离球心距离为 r1,同样截取模型的 xOz 剖面,建立如图 3 所示的数学模型。Psab(r1,)xzO1c12c2(r,)图3弹性球壳及入射球面波Fig.3Elasticsphericalshellandincidentspher
21、icalwave球面波的入射声压取 r46,n=50 时和 n=65 时目标强度曲线不一致,说明 n=50 并不满足收敛性要求,而n=65 和 68 时曲线保持不变,说明目标强度曲线已经收敛,即 n=65 即可满足要求。为保证计算结果的准确性以及计算速度的快速性,后文计算结果涉及 ka50 的计算时,均取 n=65。2.2仿真验证在 COMSOLMultiphysics 软件中,建立弹性球壳模型如图 5 所示。球壳水域PML 层图5弹性球壳几何建模Fig.5Geometricmodelingofelasticsphericalshell采用二维轴对称方法建立该模型,基本参数如下:1)单位点声源
22、激励,入射频率为 10kHz;2)点声源与球心距离为 4m,球壳材料为铝,纵波波速CL=6040ms1,横波波速 CT=3040ms1,密度0=2700kgm3;3)球壳外部流体介质为水,声速c=1500ms1,密度=1000kgm3,为了模拟无反射边界条件,在水域最外层设置完美匹配层(per-fectmatchedlayer,PML)。验证方案:保证球壳外半径 a=0.5m 不变,改变内半径 b 和球壳内部填充情况,计算距离球心10m 的散射声压。为了方便曲线更加直观显示,对散射声压取级,即 Lp=20lg(Ps/P0),P0=106Pa,并与前述的理论解进行对比,如图 6 所示。图 6 中
23、多组方案的理论计算结果与有限元(finiteelementmethod,FEM)结果对比显示,两者几乎完全吻合。由此证明,该理论方法用于计算水01020304050302520151050TS/dBkan=50n=65n=68图4不同截断数下 ka50 时距离球心 1000m 处目标强度Fig.4Target strength at a distance of 1000 m fromspherical center with ka50 under different trun-cationnumbers2023年10月郑金焱,等:点声源入射下水下弹性球壳声散射分析第5期水下无人系统学报sxwr
24、xtxb.xml-729下弹性球壳的散射特性是可行且有效的。3数值计算与讨论3.1目标强度分析3.1.1各类球形目标的频响曲线在工程中,弹性球壳通常内部不充水,鉴于此,文中考虑的弹性球壳内部均为真空。在前述定义下,先计算了相同尺寸的刚性球、铝球和铝球壳的目标强度频响曲线。球体参数为外半径 0.5m,壳厚 0.02m,点声源距离球心 4m,外部流体介质为水,如图 7 所示。05 00010 00015 00020 00025 0005040302010010TS/dBf/Hz刚性球弹性球弹性球壳图7各类球形目标的目标强度Fig.7Targetstrengthofvarioussphericalt
25、argets对比刚性球和弹性球的目标强度曲线,可以发现,在低频时,实心弹性球的散射主要贡献为刚性散射,弹性散射的贡献相对较少,随着频率的增加,弹性散射的贡献便体现出来,而且随着 ka 的不断增大,弹性散射的贡献更突出;对比刚性球和弹性球壳的频响曲线,可以看到,弹性球壳的目标强度曲线更加曲折,随着 ka 的增加,共振峰值越来越密集。而弹性球与弹性球壳对比显示,由于球壳刚度更低,球壳的弹性振动更容易被激发,在低频时就已经被激发引起共振峰。3.1.2弹性球目标强度由于弹性球壳散射特性复杂,这里首先对弹性球进行分析,讨论点声源与球心距离改变时的目标强度曲线与平面波激励下的关系。图 8 为 a=0.5m
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 声源 入射 下水 弹性 球壳声 散射 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。