CN112579967B - 一种应用抛物方程模型求解水下声场简正波模态的方法 - Google Patents
一种应用抛物方程模型求解水下声场简正波模态的方法 Download PDFInfo
- Publication number
- CN112579967B CN112579967B CN202011398003.7A CN202011398003A CN112579967B CN 112579967 B CN112579967 B CN 112579967B CN 202011398003 A CN202011398003 A CN 202011398003A CN 112579967 B CN112579967 B CN 112579967B
- Authority
- CN
- China
- Prior art keywords
- sound field
- eigenvalue
- mode
- equation model
- parabolic equation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 239000011159 matrix material Substances 0.000 claims abstract description 37
- 230000005540 biological transmission Effects 0.000 claims abstract description 12
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 7
- 238000004134 energy conservation Methods 0.000 claims abstract description 3
- 238000004364 calculation method Methods 0.000 claims description 15
- 238000012546 transfer Methods 0.000 claims description 9
- 239000013598 vector Substances 0.000 claims description 5
- 230000002441 reversible effect Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Operations Research (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种应用抛物方程模型求解水下声场简正波模态的方法,步骤如下:将能量守恒水声抛物方程模型中的多步矩阵运算简化为单步传递矩阵运算形式,并对传递矩阵进行特征值分解,根据传递矩阵的准对称性,获得一组相互正交的声场传播模态;通过与简正波理论声场表达式进行对比,得到基于抛物方程模型的简正波本征值和本征模态的近似表达式,实现声场简正波模态快速求解。本发明的有益效果为:本发明提出的简正波本征值和本征向量计算的新方法计算结果与简正波模型得出的结果吻合良好,同时新方法可以根据网格划分,灵活计算简正波的模态。
Description
技术领域
本发明涉及水声传播建模与海洋声场分析领域,主要是一种应用抛物方程模型求解水下声场简正波模态的方法。
背景技术
水声传播模型作为水声科学研究的基础组成部分,是分析水声传播规律及物理机理的重要工具,为水声设备提供基本的理论依据和技术参考,对于改善声呐的工作性能具有非常重大的现实意义。水声传播理论研究的核心问题是将实际声波传播的物理问题抽象成数学问题,即求解满足一定边界条件的波动方程问题。根据波动方程的求解方法和近似条件的不同,多种水声传播理论被陆续提出,极大地推动了声场计算和分析手段的进步。现有的水下声场传播模型主要包括简正波理论、抛物方程方法、射线模型、快速场(波数积分)方法、有限差分和有限元方法等;同时,还衍生出各种理论方法的混合模型,例如射线-简正波、简正波- 抛物方程、射线-简正波-抛物方程等声传播模型。
在处理与距离有关声传播问题时,基于简正波理论发展起来的耦合简正波模型虽然可以有效计算,但声场计算效率大打折扣。相比于其它大多数声传播理论,抛物方程模型具有快速、灵活的优点,可有效的解决与距离有关环境下的声传播问题,但其作为一种数值计算方法,难以有效分析声传播问题的内在物理机理和声场干涉结构,阻碍了该方法的进一步发展。
发明内容
本发明的目的在于克服现有技术存在的不足,而提供一种应用抛物方程模型求解水下声场简正波模态的方法。
本发明的目的是通过如下技术方案来完成的。一种应用抛物方程模型求解水下声场简正波模态的方法,步骤如下:将能量守恒水声抛物方程模型中的多步矩阵运算简化为单步传递矩阵运算形式,并对传递矩阵进行特征值分解,根据传递矩阵的准对称性,获得一组相互正交的声场传播模态;通过与简正波理论声场表达式进行对比,得到基于抛物方程模型的简正波本征值和本征模态的近似表达式,实现声场简正波模态快速求解。
本发明的有益效果为:本发明提出的简正波本征值和本征向量计算的新方法计算结果与简正波模型得出的结果吻合良好,同时新方法可以根据网格划分,灵活计算简正波的模态。
附图说明
图1为抛物方程模型计算简正波本征值和本征向量实现流程图。
图2为海洋波导环境参数分布图。
图3是声场模态本征函数结果比较示意图。
具体实施方式
下面将结合附图对本发明做详细的介绍:
本发明公开了一种应用抛物方程模型求解水下声场简正波模态的方法,步骤如下:首先,根据水声抛物方程模型递推计算方程,将多步矩阵计算简化为单步传递矩阵计算,然后对计算得到的传递矩阵进行特征值分解。根据矩阵的对称性,将抛物方程模型声场计算表达式改写成由特征值和特征向量表示的形式。最后,通过与简正波理论声场模态表达式对比,获得简正波本征值和本征向量的表达式,实现简正波模态的求解。
图1为本专利的实现过程图,采用传统抛物方程模型,计算得到传递矩阵,然后进行特征值分解,将抛物方程模型声场计算表达式改写成由特征值和特征向量表示的形式。最后,通过与简正波理论声场模态表达式对比,计算得到简正波本征值和本征函数。
图2为与距离无关海洋(水平)波导示意图,图中标出了海洋地声参数及声源信息。
图3是分别采用抛物方程模型和简正波理论对水平波导声传播模态本征函数求解,结果表明,采用抛物方程模型计算得到的模态本征函数与简正波模型计算结果吻合很好。
表1是分别采用抛物方程模型和简正波理论对水平波导声传播模态本征值求解。采用抛物方程模型计算得到的模态本征值的实部和虚部与简正波模型计算结果都高度一致,证明新方法的正确性。
本发明的具体实施方式如下:
(1)在二维柱坐标系下,水下声场满足的单向抛物方程模型声场计算公式为:
对上式中的算子采用Galerkin方法进行离散,可以获得如下的多步矩阵方程的形式:
RlUl=SlUl-1(l=1,2,...,L),
(2)抛物方程模型计算声场的过程是由上一个水平距离处的不同深度声场值向量向下一个水平距离处的不同深度声场值向量递推求解的过程,从最直观的角度看,可以写成如下的形式:
其中,矩阵ΥN×N为传递矩阵。
(3)在计算Rl和Sl的过程中发现,两个矩阵的各个元素总体上关于主对角线对称,即矩阵为准对称阵,满足:
传递矩阵Υ也为准对称阵,即:
传递矩阵为可逆矩阵,因此可以进行特征值分解,即:
Υ=ΨΓΨ-1,
Ψ=[φ1,φ2,φ3,...,φN],
Γ=diag(λ1,λ2,λ3,...,λN),
其中,列向量φn(n=1,2,...,N)为传递矩阵Υ的特征值λn(n=1,2,...,N)对应的特征向量。特征值与特征向量满足如下的关系:
Υφn=λnφn.
由于传递矩阵为准对称矩阵,而对称矩阵的特征向量之间相互正交,因此可以得出,特征向量组{φn},(n=1,2,...,N)为一组完备的准正交基,即特征向量之间内积满足:
传递矩阵对声场的作用可以分解成传递矩阵对声场各正交分量独立作用的叠加.
(4)对于与距离无关声传播问题,每个水平步长递推求解过程中,传递矩阵保持不变,因此声场的计算公式可以表示成:
U(mΔr)=Υm-1U(Δr).
采用特征值和特征向量可以表示成:
抛物方程模型的模态可以视为简正波模态的均匀离散采样形式,但由于抛物方程模型可以有效计算海底的声场,因此两者之间存在一定的差别。变量U(mΔr)经过反推可得水下声场声压表达式可以写成:
(5)通过与简正波理论声场模态表达式对比,可以得到简正波本征值和本征模态的表达式:
表1声场模态本征值计算结果比较
模态阶数 | 简正波理论 | 抛物方程方法 |
1 | 0.8376169-0.0000003i | 0.8376173-0.0000003i |
2 | 0.8371932-0.0000010i | 0.8371948-0.0000010i |
3 | 0.8364865-0.0000023i | 0.8364900-0.0000023i |
4 | 0.8354959-0.0000040i | 0.8355020-0.0000041i |
5 | 0.8342201-0.0000063i | 0.8342292-0.0000063i |
6 | 0.8326576-0.0000089i | 0.8326700-0.0000090i |
7 | 0.8308066-0.0000120i | 0.8308221-0.0000122i |
8 | 0.8286647-0.0000155i | 0.8286829-0.0000157i |
9 | 0.8262295-0.0000193i | 0.8262495-0.0000198i |
可以理解的是,对本领域技术人员来说,对本发明的技术方案及发明构思加以等同替换或改变都应属于本发明所附的权利要求的保护范围。
Claims (1)
1.一种应用抛物方程模型求解水下声场简正波模态的方法,其特征在于:步骤如下:将能量守恒水声抛物方程模型中的多步矩阵运算简化为单步传递矩阵运算形式,并对传递矩阵进行特征值分解,根据传递矩阵的准对称性,获得一组相互正交的声场传播模态;通过与简正波理论声场表达式进行对比,得到基于抛物方程模型的简正波本征值和本征模态的近似表达式,实现声场简正波模态快速求解;
具体步骤如下:
(1)在海洋环境参数水平方位变换缓慢的情况下,基于高阶Padé近似的抛物方程模型的水下声场计算公式获得如下的多步矩阵方程的形式:
RlUl=SlUl-1(l=1,2,...,L),
(2)从数学的角度看,等维列向量之间的变换通过一次等维可逆方阵实现,该方式视为向量变换的最简形式:
(3)传递矩阵为可逆矩阵,进行特征值分解,即:
Υ=ΨΓΨ-1,
Ψ=[φ1,φ2,φ3,...,φN],
Γ=diag(λ1,λ2,λ3,...,λN),
(4)抛物方程模型的模态视为简正波模态的均匀离散采样形式,U(mΔr)经过反推得水下声场声压表达式写成:
(5)通过与简正波理论声场模态表达式对比,得到简正波本征值和本征模态的表达式:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011398003.7A CN112579967B (zh) | 2020-12-04 | 2020-12-04 | 一种应用抛物方程模型求解水下声场简正波模态的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011398003.7A CN112579967B (zh) | 2020-12-04 | 2020-12-04 | 一种应用抛物方程模型求解水下声场简正波模态的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112579967A CN112579967A (zh) | 2021-03-30 |
CN112579967B true CN112579967B (zh) | 2024-01-05 |
Family
ID=75127497
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011398003.7A Active CN112579967B (zh) | 2020-12-04 | 2020-12-04 | 一种应用抛物方程模型求解水下声场简正波模态的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112579967B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117725762B (zh) * | 2024-02-06 | 2024-05-24 | 西北工业大学 | 一种基于时域谱展开的水下瞬态声场仿真方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107016159A (zh) * | 2017-02-28 | 2017-08-04 | 浙江海洋大学 | 本征值确定方法及装置 |
CN107885934A (zh) * | 2017-11-07 | 2018-04-06 | 哈尔滨工程大学 | 基于耦合fem‑pe的海洋信道下弹性结构声辐射预报方法 |
CN110135052A (zh) * | 2019-05-12 | 2019-08-16 | 哈尔滨工程大学 | 浅海信道下弹性结构辐射声场的计算方法 |
-
2020
- 2020-12-04 CN CN202011398003.7A patent/CN112579967B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107016159A (zh) * | 2017-02-28 | 2017-08-04 | 浙江海洋大学 | 本征值确定方法及装置 |
CN107885934A (zh) * | 2017-11-07 | 2018-04-06 | 哈尔滨工程大学 | 基于耦合fem‑pe的海洋信道下弹性结构声辐射预报方法 |
CN110135052A (zh) * | 2019-05-12 | 2019-08-16 | 哈尔滨工程大学 | 浅海信道下弹性结构辐射声场的计算方法 |
Non-Patent Citations (2)
Title |
---|
一种水平变化可穿透波导中声传播问题的耦合简正波方法;杨春梅;骆文于;张仁和;秦继兴;;物理学报(第09期);全文 * |
基于WKBZ理论的耦合简正波-抛物方程理论;彭朝晖, 李风华;中国科学A辑(第02期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112579967A (zh) | 2021-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107885934B (zh) | 基于耦合fem-pe的海洋信道下弹性结构声辐射预报方法 | |
CN103941220B (zh) | 一种基于稀疏重构的网格外目标波达方向估计方法 | |
CN107100611B (zh) | 一种地质勘探钻进过程双层智能钻速建模方法 | |
CN112285788B (zh) | 一种基于电磁波动方程的cpml吸收边界条件加载方法 | |
CN101604019B (zh) | 一种海洋环境与声场不确实性的表征和传递的快速计算方法 | |
Benedetti et al. | A fast dual boundary element method for 3D anisotropic crack problems | |
CN107844627A (zh) | 一种仅输出时变结构模态参数贝叶斯估计方法 | |
KR20170056687A (ko) | 시퀀스 재귀 필터링 3차원 변분(3d-var) 기반의 실측 해양 환경 데이터 동화방법 | |
CN113050075B (zh) | 一种基于扩散映射的水下声源匹配场定位方法 | |
CN112579967B (zh) | 一种应用抛物方程模型求解水下声场简正波模态的方法 | |
CN109933877B (zh) | 代数多重网格三维变分数据同化方法 | |
CN116304647A (zh) | 一种基于瞬态响应的浮式结构频响函数识别方法 | |
CN112014791A (zh) | 存在阵列误差的阵列pca-bp算法的近场源定位方法 | |
CN115906303A (zh) | 一种基于机器学习的平面微波滤波器设计方法及装置 | |
CN117148460A (zh) | 一种接地长线源瞬变电磁法的高效反演方法 | |
CN117473734A (zh) | 一种基于模态展开的封闭空间瞬态声场仿真方法 | |
CN109164488B (zh) | 一种梯形网格有限差分地震波场模拟方法 | |
CN114662040A (zh) | 一种周期起伏海底声场简正波模态计算方法 | |
CN113935166A (zh) | 一种浅海信道下三维结构声场预报方法、系统及存储介质 | |
Shang et al. | Stochastic dynamic analysis of acoustic-structural coupled systems under non-stationary random excitations | |
CN115859729A (zh) | 基于深度学习的水下井口系统疲劳损伤评估方法及装置 | |
CN115688307A (zh) | 考虑土体强度不确定性的隧道掌子面失稳风险评估方法 | |
CN110738003B (zh) | 一种重型拖拉机pto壳体时变可靠性分析方法 | |
CN111428342B (zh) | 一种基于频域谱分解的随机动载荷识别方法 | |
CN111259578B (zh) | 基于半解析法和边界元法预测直筒笼式阀门振动噪声方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |