基于声辐射模态的水下航行器传感器优化配置研究.pdf

第45卷第16 期2023年8 月舰船科学技术SHIP SCIENCEANDTECHNOLOGYVol.45,No.16Aug.,2023基于声辐射模态的水下航行器传感器优化配置研究汤智胤，邹海龙1，喻柳丁1，李想，沈惠杰1，马锐锋2(1.海军工程大学动力工程学院，湖北武汉430 0 33；2.海军工程大学军用电气科学与技术研究所，湖北武汉430 0 7 2)摘要：本文研究当水下航行器外壳体振动传感器总数确定之后，如何利用这有限个数的传感器得到更多航行器壳体振动信息，从而更准确地评估航行器当前声隐身状态的问题。基于结合声辐射模态理论，给出传感器优化配置方法和步骤，并利用类似水下航行器的椭球壳体结构作为研究对象进行仿真，对该方法进行验证。研究结果表明，针对类椭球结构的水下航行器，低频段壳体传感器应布置在两端及中间位置。关键词：声辐射模态；传感器优化配置；水下航行器中图分类号：TP391文章编号：16 7 2-7 6 49(2 0 2 3)16-0 10 8-0 4Research on underwater vehicle optimal sensor configuration based onTANG Zhi-yin,ZOU Hai-long,YU Liu-ding,LI Xiang,SHEN Hui-jie,MA Rui-feng?(1.College of Power Engineering,Naval University of Engineering,Wuhan 430033,China;2.Institute of Military Electrical Science and Technology,Naval University of Engineering,Wuhan 430033,China)Abstract:The problem of how to arrange the sensors on the shell of the underwater vehicle aimed to assess the vehicleacoustic stealth situation better was researched.With the use of acoustic radiation mode theory,the method and steps of op-timal sensor placement has been put forward.The simulation analysis of a quasi-elliptic shell structure similar to the com-mon underwater vehicle has been given to validate the method in addition.The results showed that aimed at low frequencyband the sensors should be located on the both ends and intermediate position of the structure.Key words:acoustic radiation modes;optimized sensors configuration;underwater vehicle0引言海战场对抗环境下水下航行器的隐身性，特别是声隐身性是其最基本和最重要的技术性能-3。评估水下航行器的隐身安全性，一种常用的方法是依靠航行器壳体表面安装的振、声传感器测量得到。然而在航行器的振、声传感器数量确定之后，如何能利用有限的传感器得到更多信息，从而更精确地评估水下航行器的隐身安全性呢？这就需要进行水下航行器的传感器优化配置研究。对于航行器壳体传感器，通常最主要且数量最多的传感器是壳体加速度传感器4。对于大量的加速度传感器，如何能够合理地将其布置在航行器壳体上，收稿日期：2 0 2 2-0 7-2 2基金项目：湖北省自然科学基金青年基金资助项目(2 0 2 1CFB266)作者简介：汤智胤（19 8 1-），男，博士，副教授，研究方向为船舶辅助机械。文献标识码：Aacoustic radiation modesdoi:10.3404/j.issn.1672-7649.2023.16.022最大程度地采集航行器的振动信息，从而更加精确地进行声隐身安全性评估，需要认真探索。本文利用声辐射模态理论，对航行器壳体加速度传感器优化配置进行研究。1传感器优化配置方法及步骤近年来，国内外学者提出了声辐射模态理论5-8 。声辐射模态及其辐射效率类似于振动问题中的主振形和固有频率。声辐射模态是矢量空间中一组相互正交的基，每组基代表一种可能的辐射形式。每一声辐射模态对应一个独立的辐射效率。声辐射模态的形状由辐射体的几何形状和振动频率决定9 。在中、低频时，声辐射模态的辐射效率随着阶数的增加而快速降第45卷低，所以仅利用前几阶声辐射模态就能近似计算出总声功率。每一阶声辐射模态的声功率是相互独立的。利用声辐射模态理论研究声辐射问题的优势在于消除了结构模态中复杂的耦合项，使研究更简单10-16-18 。若仅针对水下航行器辐射噪声预报问题进行传感器优化配置，其目标是寻找对于辐射声功率最大的前若干阶声辐射模态进行监测，从而进行传感器位置的优化布置。然而对于水下航行器声隐身状态评估而言，并不是辐射噪声能量越大的声辐射模态，对水下航行器声隐身状态贡献也越大。因为水下航行器声隐身状态还取决于对抗方声呐的探测性能和海洋声传播损失12。所以，对于面向水下航行器声隐身状态评估的传感器优化配置，除了考虑声辐射模态的辐射效率外，还需考虑传播损失以及对抗方声呐在不同频率上的不同特性。基于水下航行器声隐身状态评估的艇载传感器优化配置可按照以下步骤进行：1）确认水下航行器频谱中对声隐身状态影响最大的频率成分进行航行器水下航行实验时，测量其辐射噪声，并针对航行器声隐身状态中最敏感的频率成分，与声呐测量和水下声传播相结合，锁定最敏感频率成分。令航行器水下辐射噪声的线谱成分分别为fi=1,2,，n，其对应频率的声强度分别为sl(fi)i=1,2,n；设航行器遭遇最常见的对抗方声呐在这些频率成分上的探测能力分别为ds(fi)i=1,2,n；航行器通常执行任务的海区在这些频率成分上的传播损失分别为tl(fi),i=1,2,.,n。令函数s(f)ds(f)=,.,n.st(fi)=tl(fi)可表示航行器声隐身状态的安全程度，则频率fst则为st(fi)i=1,2,.,n取最大值所对应的频率，即fst=argmax(st(fi),st(f2),.,st(fn)。2）将该频率按照声辐射模态进行分解，得到对声隐身状态贡献最大的前若干阶声辐射模态在进行传感器优化配置之前，可先将传感器均匀地布置在壳体上，作为传感器初步配置方案。这样将航行器振动表面划分成若干面积相等的单元，在某频率成分各单元上的法向速度构成的向量为vER8X1，对应各单元的表面声压为p，则p=Zv其中，Z为每个单元的振动速度与表面声压之间的声阻抗矩阵。而辐射声功率为：汤智胤，等：基于声辐射模态的水下航行器传感器优化配置研究2式中：S为每个单元的面积；Re()表示取实部；上标H表示复共轭转置；R=SR e(Z)为实正定矩阵8 。对R进行特征值分解为R=QQT的形式。因为R为实正定矩阵，所以QH=QT，所以声功率表示为：(5)i=1i=1式中，Q为声辐射模态矩阵，Q=q1,q2,qcl。其中qi，i=1,2,c 为一种可能的相互正交的速度分布，其代表了一种声辐射模态，故各阶声辐射模态下的声功率相互独立。为对角阵，其中第i个对角元素为入i，均为正实数,且i2.c；y=y 1,2,y e T=Q T v,lyi伪第阶声辐射模态的伴随系数。评估共研究前c阶声辐射模态。研究发现，在中低频由前一、二阶辐射模态产生的声功率约占总声功率90%以上，而且在中、低频时声辐射模态形状随频率的变化较小7-8 。由于既考虑了声辐射功率，又考虑了海洋传播损失和对抗方声呐的探测能力，所以对于频率fst上的辐射声贡献越大的声辐射模态，则对于航行器声隐身状态的贡献也越大，需投入尽量多的传感器进行监测；对于S贡献较小的声辐射模态，则对声隐身状态贡献较小，可以忽略。3）结合传感器数量限制，针对前若干阶声辐射模态型进行传感器布置，最有效地监测这些声辐射模态(1)当算出前若干阶对航行器声隐身状态贡献较大的声辐射模态后，针对这些声辐射模态进行传感器优化布置。传感器优化配置时，只需将传感器按照其数量要求，来监测前若干个对航行器声隐身状态贡献较大(2)声辐射模态。对于这种模态已知情况下的传感器优化配置问题，国内外已经研究得较为成熟。其中常用的方法包括有效独立法、Guyan模型缩减法等。有效独立法利用Fisher信息阵使关心的模态向量尽可能线性无关，从而采集到最大的模态反应信息13；Guyan模型缩减法通过刚度或质量子矩阵构成的转换矩阵，把那些对模态反应起主要作用的自由度保留下来作为测点的位置14；Thomas等15 认为选择的测点位置应使所有模态(3)向量的内积具有较小的余弦值。这些优化方法都能起到使测点远离各振型节点的作用，从而利用较少的传感器监测关心的模态。109.1W=ASRe(vHHp)=ASRe(HZv)=vHRv。22(4)C:110利用声辐射模态理论进行面向航行器声隐身状态评估的传感器优化配置，最主要的优点在于消除了结构模态中复杂的耦合项，且在中、低频时的辐射效率随声辐射模态阶数的增加而迅速降低，少数的前几阶声辐射模态对航行器的声隐身状态有很大的贡献16 。而本方法的难点是需事先计算出航行器的声阻抗矩阵，这是实船应用中最主要的问题。2仿真验证取类似水下航行器的椭球壳体结构作为研究对象，椭球模型长轴为5m，两短轴均为1m；结构材料密度p=7850kg/m，结构材料的杨氏模量E=2.110，损耗因子为0.3，外场介质为水，声速co=1500m/s。假设结构的激励同航行器类似存在若干种固定的激励模式，即存在有限个工况。并假定在5种常用工况下，对声隐身状态影响最大的5个频率的辐射噪声分别为：10 Hz（k l=0.42，k 为波数；1为结构的最大尺寸）、8 0 Hz（k l=3.35）、10 0 H z（k l=4.19）、125Hz（k l=5.2 4）、16 0 H z（k l=6.7 0）。利用声辐射模态理论，计算得到这5个频率前4阶声辐射模态形式，如图1 图5所示。这5个频率前4阶声辐射模态对应的特征值如表1所示。可以看出，在低频（kl2元，k为波数；1为结构最大尺寸）时，各阶声辐射模态对应的特征值入；（i=1,2,3,4）随模态阶数的增加而迅速降低，即各阶声辐射模态对声隐身状态的贡献程度随模态阶数的增大迅速减小。特别地，前2 阶模态所对应的特征值要比其他高阶模态的特征值大得多，即对声隐身状态的贡献也大得多。而且，从图1图5可观察到，5个频率的前2阶声辐射模态的形式基本相同，即中间和两端部较强，且为轴向对称。所以，可根据壳体振动加速度传感器数量限制，着重监测前2 阶声辐射模态，即两端和中间部分。当然，在传感器数量较大情况下，除了(a)第1阶声辐射模态(c)第3阶声辐射模态图110 Hz时结构的前4阶声辐射模态Fig.1 The first 4 acoustic radiation modes of the structure in 10 Hz舰船科学技术(a)第1阶声辐射模态(c)第3阶声辐射模态图2 8 0 Hz时结构的前4阶声辐射模态Fig.2 The first 4 acoustic radiation modes of the structure in 80 Hz(a)第1阶声辐射模态(c)第3阶声辐射模态图310 0 Hz时结构的前4阶声辐射模态Fig.3 The first 4 acoustic radiation modes of thestructure in 100 Hz(a)第1阶声辐射模态(b)第2 阶声辐射模态(c)第3阶声辐射模态(d)第4阶声辐射模态图412 5Hz时结构的前4阶声辐射模态Fig.4 The first 4 acoustic radiation modes of thestructurein125Hz前2 阶声辐射模态外，还可由低到高顺序针对其他声(b)第2 阶声辐射模态辐射模态进行传感器优化布置，以达到更好地进行声隐身状态评估的目的。3结语(d)第4阶声辐射模态本文针对在水下航行器壳体振动传感器数量确定的情况下，如何更好地对传感器的位置进行优化配置，以提高水下航行器声隐身状态评估准确度问题展第45卷(b)第2 阶声辐射模态(d)第4阶声辐射模态(b)第2 阶声辐射模态(d)第4阶声辐射模态第45卷(a)第1阶声辐射模态(c)第3阶声辐射模态图516 0 Hz时结构的前4阶声辐射模态Fig.5 The first 4 acoustic radiation modes of thestructure in 160 Hz表1不同频率各阶声辐射模态对应的特征值Tab.1 Eigenvalues corresponding to each acoustic radiation modesin differentfrequency第1阶模态对第2 阶模态对第3阶模态对第4阶模态对计算频率/Hz应特征值入应特征值入应特征值入3应特征值入4102.50580114.490100154 040125200300160259180开研究，提出了基于声辐射模态理论的传感器优化配置方法，给出了初步的优化配置方案，并通过算例分析证明了方法的实用性。由于实际水下航行器的数据难以获取，本文算例以类似水下航行器的椭球壳体结构为研究对象，说明方法的可行性。虽然实际水下航行器结构比算例中所研究的规则椭球体结构要复杂得多，然而在低频时，可类似一个椭球壳体。本方法用于实际航行器传感器优化配置的一个主要难点在于计算航行器的声阻抗矩阵。在实际应用中，可利用有限元边界元方法进行实船建模计算，算出声阻抗矩阵后再利用本方法进行传感器优化布置，以取得更准确的航行器声隐身状态评估结果。参考文献：1 FENG X,YAN L,CHUN W.Review on vibration isolationtechnologyJ.Journal of Physics:Conference Series,2021,1820:15.【2 苏强，王桂波，朱鹏飞，等.国外潜艇声隐身前沿发展综述1.舰船科学技术,2 0 14,36(1):1-9.SU Qiang,WANG Guibo,ZHU Pengfei,et al.A review of thefrontier development of submarine acoustic stealth abroadJ.汤智胤，等：基于声辐射模态的水下航行器传感器优化配置研究(b)第2 阶声辐射模态动与冲击,2 0 0 8,2 7(7):113-117.TANG Zhiyin,JIN Guangwen,HE Lin.Research on rapidevaluation of submarine acoustic stealth statusJ.Vibration(d)第4阶声辐射模态and Shock,2008,27(7):113-117.5 CURNFARE K A.The radiation modes of baffled finiteplate9.Journal Acoustical Society of America,1995,98(3):1570-1580.6 GUO Liang,ZHU Hai-chao,MAO Rong-fu,et al.Sound fieldseparation technique based on acoustic radiation modesJ.Journal of Ship Mechanics,2017,21(6):779-790.【7】苏俊博，朱海潮，毛荣富，等.基于声辐射模态的声场重建中的测点优化方法J.振动与冲击，2 0 17,36(3):145-150.SU Junbo,ZHU Haichao,MAO Rongfu,et al.Optimization13.74.828187713355 5721701993.5873891014342087 224111Ship Science and Technology,2014,36(1):1-9.3何琳.潜艇隐身技术进展.舰船科学技术,2 0 0 6,2 8(S2):9-17.HE Lin.Progress in submarine stealth technologyJ.Ship Sci-ence and Technology,2006,28(S2):9-17.【4汤智胤，金广文，何琳.潜艇声隐身状态快速评估研究.振2.7method of measurement points in sound field reconstruction7130based on acoustic radiation modeJ.Vibration and Shock,16.9962017,36(3):145150.388848聂永发，朱海潮.利用源强密度声辐射模态重建声场.物87154理学报,2 0 14(10):2 56-2 6 7.9聂永发，朱海潮.复杂结构声辐射模态的计算.应用声学，2014,33(6):534540.10鱼海涛，王英民，王奇.利用声辐射模态重构任意目标的散射声场.应用声学,2 0 17,36(3):2 6 4-2 7 5.11 吴吴锦武，原海朋，毛崎波.层合板结构声辐射模态传感器分析与测量J.振动与冲击，2 0 16,35(13):2 6-30.12】汤智胤,姜荣俊，何琳.潜艇声隐身态势评估方法研究1.武汉理工大学学报(交通科学与工程版),2 0 0 7,31(1):17-2 0.13 KAMMER D C.Sensor placements for on-orbit modal identifi-cation of large space structuresJJ.Journal of Guidance,Con-trol,and Dynamics,1991,14(2):252-259.14 GUYAN R I.Reduction of stffness and mass matrics.AIAA Journal,1995,3(2):380.15 THOMAS G C,CLARK R D.A modal test design strategy formodal correlationCJ/In:Proc 13th International Modal Analy-sis Conference.New York,1995:927-933.16】毛崎波，姜哲.通过声辐射模态研究结构噪声的有源控制.江苏理工大学学报(自然科学版),2 0 0 0(4):1-6.17】石炜.矩形薄板的振动与声辐射研究及其控制D.成都：西南交通大学,2 0 10.18】毛崎波，姜哲.通过声辐射模态研究结构声辐射的有源控制.声学学报,2 0 0 1(3):2 7 7-2 8 1.MAO Qibo,JIANG Zhe.Research on active control of struc-tural acoustic radiation through acoustic radiation modalJ.Journal of Acoustics,2001(3):277-281.