一种针对浅海分层弹性海底的浅海声场计算方法技术

本发明专利技术公开一种针对浅海分层弹性海底的浅海声场计算方法，属于水声工程技术领域，计算方法基于快速场（FFP）声场计算方法和点声源作用下模拟运算了多层弹性海底模型建立了针对浅海分层弹性海底的浅海声场计算方法，本发明专利技术所构建的多层弹性海底模型在弹性介质半空间层和均匀立体介质层间还具有弹性介质层，弥补了现有技术中基于半无限海底模型的缺陷，利用该方法反演得到的声传播损失与实际测量值具有更高的吻合度；结合快速场计算浅海声场，计算方法快速准确，既包含了离散模式部分，又包含了连续部分。

A method of calculating shallow sea sound field for shallow sea layered elastic seabed

The invention discloses a method for calculation of shallow water in shallow layered elastic seabed, which belongs to the technical field of underwater acoustic engineering, calculation method based on fast field (FFP) field calculation method and effect of point source under the simulation calculation method of multilayer elastic seabed model established in shallow water in shallow layered elastic seabed, multilayer elastic model of the submarine the invention of the flexible dielectric layer in elastic half space medium layer and a uniform three-dimensional dielectric layer, defects in semi infinite seabed model based on existing technology, sound propagation and actual measurement by using the inversion method to the value of having a higher degree of agreement; combined with the rapid shallow water field calculation, calculation the method is rapid and accurate, which includes not only the discrete model, and contains a continuous part.

【技术实现步骤摘要】
一种针对浅海分层弹性海底的浅海声场计算方法
本专利技术属于水声工程领域，尤其是涉及一种针对浅海分层弹性海底的浅海声场计算方法。
技术介绍
在水声领域的研究中，声场建模和海底声参数测量一直都是研究的热点。随着水声工程技术的不断发展，在国内外水声工作者的共同努力下，针对水平分层模型下的各种浅海声场计算方法，如：简正波理论、射线理论、抛物方程、快速场理论等日趋完善。在声场建模的基础上，采用声学反演方法可以快速、有效的测定环境参数，与传统测量方法相比节省了大量的人力、物力，因而也越来越受到水声研究者的关注。在已有研究中，对具有弹性海底的浅海波导中的声场计算研究多是基于半无限海底模型展讨论。考虑在实际海洋环境中，海底往往是有一层沉积层和半无限弹性体构成的两层介质，在部分地区甚至需考虑为多层介质，为满足浅海声场实际建模需要，需建立针对浅海分层弹性海底的浅海声场计算方法。
技术实现思路
本专利技术的目的在于提供一种针对浅海分层弹性海底的浅海声场计算方法，该方法弥补了现有技术中基于半无限海底模型的缺陷。本专利技术的另一目的在于提供载有一种搭载有针对浅海分层弹性海底的浅海声场计算方法的设备，装载该移动硬盘的主机散热、除尘效果好，异物无法进入机箱内。本专利技术为实现上述目的所采取的技术方案为：一种针对浅海分层弹性海底的浅海声场计算方法，计算方法基于快速场(FFP)声场计算方法和点声源作用下模拟运算了多层弹性海底模型建立了针对浅海分层弹性海底的浅海声场计算方法。多层弹性海底模型为：点声源位于厚度为h1，密度及声速分别为ρ1、c1的均匀流体介质层中，z＝0为海面，海面向下z为正值，声源深度为z0，第二层为各向同性的弹性介质，第三层为各向同性的弹性介质半空间，其密度、压缩波声速和剪切波声速分别为ρ2、cp2、cs2和ρ3、cp3、cs3，流体层中变量用位移势函数φ1表示，第二层弹性海底中势函数有标量势函数φ2和矢量势函数第三层中势函数有标量势函数φ3和矢量势函数声源为单频谐和点声源，则各势函数满足的波动方程分别是：其中代表控制方程；k1＝ω/c1代表水中波数；k2＝ω/cp2代表第二层弹性体中压缩波波数；χ2＝ω/cs2代表第二层弹性体中剪切波波数；k3＝ω/cp3代表半无限空间弹性体中压缩波波数；χ3＝ω/cs3代表半无限空间弹性体中剪切波波数；当矢量势函数中只有一个分量ψ2θ和ψ3θ不为零时，此时ψ2θ、ψ3θ满足的波动方程为：将上述各势函数表示成Fourier-Bessel积分形式：得到关于Z1(z,ξ)、Z2(z,ξ)、Z3(z,ξ)和f2θ(z,ξ)、f3θ(z,ξ)的常微分方程。根据海面的压力释放条件和弹性海底无穷远辐射条件选择，Z1(z,ξ)、Z2(z,ξ)、Z3(z,ξ)和f2θ(z,ξ)、f3θ(z,ξ)分别如下：其中在柱坐标系下，r,z方向的位移分别用q,w表示，于是根据点源条件、z＝h1处的流体/弹性边界条件(法向位移连续、法向应力连续和切向应力为零)和z＝h2处的弹性/弹性边界条件(法向位移、法向应力、切向位移和切向应力均连续)可建立9个方程，整理为矩阵形式可表示为：其中：a11＝sinβ1z0；a12＝-sinβ1z0；a13＝-cosβ1z0；a14＝a15＝a16＝a17＝a18＝a19＝0；a21＝cosβ1z0；a22＝-cosβ1z0；a23＝sinβ1z0；a24＝a25＝a26＝a27＝a28＝a29＝0；a31＝0；a32＝β1cosβ1h1；a33＝-β1sinβ1h1；a38＝a39＝0；a41＝0；a42＝sinβ1h1；a43＝cosβ1h1；a48＝a49＝0；a58＝a59＝0；a61＝a62＝a63＝0；a71＝a72＝a73＝0，a81＝a82＝a83＝0；a91＝a92＝a93＝0；利用矩阵运算可以求解得到A、B、C、D1、D2、E1、E2、F、G各参数在仿真条件下的具体数值，根据求解得到的各层势函数系数，基于快速场理论可将其离散化为：即可以求解得到多层弹性海底模型中各层势函数数值解，通过势函数与各声场物理量间的计算关系，可求解得到各层声场，上述技术方法既包含了离散模式，又包含了连续模式。该模型为三层模型，设置科学，合理，对多层弹性海底模型的计算具有很好的指导意义，所构建的多层弹性海底模型在弹性介质半空间层和均匀立体介质层间还具有弹性介质层，弥补了现有技术中基于半无限海底模型的缺陷，利用该方法反演得到的声传播损失与实际测量值具有更高的吻合度；结合快速场计算浅海声场，计算方法快速准确，既包含了离散模式部分，又包含了连续部分。搭载有针对浅海分层弹性海底的浅海声场计算方法的设备，包括移动硬盘，移动硬盘设在主机箱体上，主机箱体内部设有循环散热管，循环散热管内填充有散热剂，且在循环散热管端部连接有电磁阀，电磁阀与控制器连接，主机箱体侧面不规则布有三个散热孔，散热孔内置散热风扇，散热孔顶面连接有防尘罩，循环散热管的端部可与外接水箱连接，并且利用散热剂来对主机内部进行散热，为更好的扩大主机箱体内部的散热效果，保护主机内部元件，特别是载有针对浅海分层弹性海底的浅海声场计算方法的移动硬盘。为优化上述技术方案，采取的措施还包括：散热风扇内倾斜环绕设置扇叶，相邻扇叶的边缘中部通过弯形支撑板连接，各扇叶连接有弧形扰流片，扰流片的弧面与扇叶旋转方向相反，扰流片高度等于扇叶两端面的高度差距离，扰流片上开设有半月形扰流孔，扰流孔与扰流片竖直面的倾斜角为15°-35°，设计的散热风扇将散热冷气吹到主机箱各表面对其散热，还可出去部件表面的灰尘杂质保证主机运行的稳定性、安全性，散热风扇的扇叶之间通过弯形支撑板的连接有效提升扇叶整体稳定效果，并防止扇叶在运行过程中有坚硬物体如铅笔、螺丝等物体进入机箱内损坏机箱内部部件或扇叶影响散热风扇的散热效果，通过在扇叶表面设置扰流板和带有倾斜角的半月形扰流孔可扰乱扇热风扇旋转带动的气流，使气流更容易到达主机箱体各部位，对各部件散热，且散热均匀，并且在扰乱散热风扇的同时扰流板自身也产生了与散热风扇相反的气流，该气流在扰乱散热风扇旋转带出的气流的同时可利用自身产生的气流对与扰流板位置对应的部件散热。与现有技术相比，本专利技术的优点在于：本专利技术公开了一种针对浅海分层弹性海底的浅海声场计算方法，所构建的多层弹性海底模型在弹性介质半空间层和均匀立体介质层间还具有弹性介质层，弥补了现有技术中基于半无限海底模型的缺陷，利用该方法反演得到的声传播损失与实际测量值具有更高的吻合度；结合快速场计算浅海声场，计算方法快速准确，既包含了离散模式部分，又包含了连续部分。附图说明图1为具有两层弹性海底的浅海声场环境模型示意图；图2为实施例1中表1仿真条件下流体层中声压传播损失曲线仿真结果；图3为实施例1中表2仿真条件下流体层中声压传播损失曲线仿真结果；图4为一种针对浅海分层弹性海底的浅海声场计算设备的结构示意图；图5为散热风扇的结构示意图；图6为扰流叶结构示意图。附图标记说明：1主机箱体；2散热风扇；201扇叶；202扰流叶；203支撑板；204扰流孔；3防尘罩；4散热孔；5循环散热管；6电磁阀；7控制器；8移动硬盘。具体实施方式以下结合实施例和附图对本专利技术作进一步详细描述：实施例1：一种针对浅海分层弹性本文档来自技高网...

【技术保护点】
一种针对浅海分层弹性海底的浅海声场计算方法，其特征在于：所述计算方法基于快速场(FFP)声场计算方法和点声源作用下模拟运算了多层弹性海底模型建立了针对浅海分层弹性海底的浅海声场计算方法。

【技术特征摘要】
1.一种针对浅海分层弹性海底的浅海声场计算方法，其特征在于：所述计算方法基于快速场(FFP)声场计算方法和点声源作用下模拟运算了多层弹性海底模型建立了针对浅海分层弹性海底的浅海声场计算方法。2.根据权利要求1所述一种针对浅海分层弹性海底的浅海声场计算方法，其特征在于：所述多层弹性海底模型为：点声源位于厚度为h1，密度及声速分别为ρ1、c1的均匀流体介质层中，z＝0为海面，海面向下z为正值，声源深度为z0，第二层为各向同性的弹性介质，第三层为各向同性的弹性介质半空间，其密度、压缩波声速和剪切波声速分别为ρ2、cp2、cs2和ρ3、cp3、cs3，流体层中变量用位移势函数φ1表示，第二层弹性海底中势函数有标量势函数φ2和矢量势函数第三层中势函数有标量势函数φ3和矢量势函数声源为单频谐和点声源，则各势函数满足的波动方程分别是：其中代表控制方程；k1＝ω/c1代表水中波数；k2＝ω/cp2代表第二层弹性体中压缩波波数；χ2＝ω/cs2代表第二层弹性体中剪切波波数；k3＝ω/cp3代表半无限空间弹性体中压缩波波数；χ3＝ω/cs3代表半无限空间弹性体中剪切波波数；当矢量势函数中只有一个分量ψ2θ和ψ3θ不为零时，此时ψ2θ、ψ3θ满足的波动方程为：将上述各势函数表示成Fourier-Bessel积分形式：得到关于Z1(z,ξ)、Z2(z,ξ)、Z3(z,ξ)和f2θ(z,ξ)、f3θ(z,ξ)的常微分方程。3.根据权利要求2所述一种针对浅海分层弹性海底的浅海声场计算方法，其特征在于：根据海面的压力释放条件和弹性海底无穷远辐射条件选择，Z1(z,ξ)、Z2(z,ξ)、Z3(z,ξ)和f2θ(z,ξ)、f3θ(z,ξ)分别如下：其中在柱坐标系下，r,z方向的位移分别用q,w表示，于是4.根据权利要求3所述一种针对浅海分层弹性海底的浅海声场计算方法，其特征在于：根据点源条件、z＝h1处的流体/弹性边界条件(法向位移连续、法向应力连续和切向应力为零)和z＝h2处的弹性/弹性边界条件(法向位移、法向应力、切向位移和切向应力均连续)可建立9个方程，整理为矩阵形式可表示为：其中：a11＝...

【专利技术属性】
技术研发人员：祝捍皓，郑红，林建民，汤云峰，郑广学，
申请(专利权)人：浙江海洋大学，
类型：发明
国别省市：浙江,33