CN104778321B - 一种复杂形状声学空间的分析预报方法 - Google Patents
一种复杂形状声学空间的分析预报方法 Download PDFInfo
- Publication number
- CN104778321B CN104778321B CN201510175544.6A CN201510175544A CN104778321B CN 104778321 B CN104778321 B CN 104778321B CN 201510175544 A CN201510175544 A CN 201510175544A CN 104778321 B CN104778321 B CN 104778321B
- Authority
- CN
- China
- Prior art keywords
- acoustic space
- acoustic
- matrix
- coupling
- subspace
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000004458 analytical method Methods 0.000 title claims abstract description 16
- 239000011159 matrix material Substances 0.000 claims abstract description 52
- 230000008878 coupling Effects 0.000 claims abstract description 39
- 238000010168 coupling process Methods 0.000 claims abstract description 39
- 238000005859 coupling reaction Methods 0.000 claims abstract description 39
- 238000010606 normalization Methods 0.000 claims description 12
- 230000010354 integration Effects 0.000 abstract 1
- 238000009795 derivation Methods 0.000 description 4
- 238000009792 diffusion process Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000003542 behavioural effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005381 potential energy Methods 0.000 description 1
Abstract
本发明公开了一种复杂形状声学空间的分析预报方法。将复杂形状声学空间对象划分为N个子空间，N>2，每相邻的两个子空间通过耦合面耦合；选择每一个子空间的声压函数，确定每一个子空间的刚度矩阵和质量矩阵；确定声学空间的耦合能量，得到总体耦合刚度；对所有N个子空间的质量矩阵、刚度矩阵以及总体耦合刚度矩阵分别进行整合排列，形成声学空间的总体质量矩阵和刚度矩阵，得到声学空间的特性方程；根据声学空间的特性方程，得到声学预报信息。本发明具有收敛速度快、适用范围广的优点。
Description
技术领域
本发明属于复杂声学空间的声学分析方法，尤其涉及一种复杂形状声学空间的分析预报方法。
背景技术
封闭空间声场的理论研究回答的是给定了形状且已知了墙壁材料的房间内声波是怎么传播的这样一个物理方面的问题，程建春，田静，创新与和谐:中国声学进展.科学出版社，2008年。对封闭空间声场的研究方法通常分为两类，即几何声学和波动声学(XiangN.,Summers J.E.,Acoustics in coupled rooms:modelling and dataanalysis.International Symposium on Room Acoustics,Seville,2007.)。基于声音沿直线传播这一假设的研究方法统称为几何声学，包括统计声学模型，Summers J.E.,TorresR.,Shimizud Y.,Statistical-acoustics models of energy decay in systems ofcoupled rooms and their relation to geometrical acoustics.Journal of theAcoustical Society of America,2004,116:958–969.，扩散方程模型Jing Y.,Xiang N.,A modified diffusion equation for room-acoustic prediction(L).Journal of theAcoustical Society of America,2007,121:3284–3287；Valeau V.,Picaut J.,HodgsonM.,On the use of a diffusion equation for room-acoustic prediction.Journal ofthe Acoustical Society of America,2006,119:1504–1513等。但是，声音的本质是一种波动现象。这样，将声音看作声线或者是粒子的直线传播的几何声学无法揭示声场的声学特性细节，且这一假设限制了几何声学在高频的发展。于是，一些学者从波动的角度来研究室内声学。莫尔斯在声学中引入简正波理论，求得了室内声场的简正波解Morse P.M.,Vibration and sound.McGraw-Hill,1936。然而，在将波动理论拓展至具有不规则几何形状和不均匀界面声学特性的一般性声学问题时，由于声场几何形状和边界条件无法精确地用数学公式表示，所以只能提出形式解，无法得到精确的解析解：程建春，田静，创新与和谐，中国声学进展.北京：科学出版社，2008年。
发明内容
本发明的目的是提供一种计算速度快、适用范围广的，一种复杂形状声学空间的分析预报方法。
本发明是通过以下技术方案实现的：
一种复杂形状声学空间的分析预报方法，包括以下步骤，
步骤一：将复杂形状声学空间对象划分为N个子空间，N>2，每相邻的两个子空间通过耦合面耦合；
步骤二：选择每一个子空间的声压函数，确定每一个子空间的刚度矩阵和质量矩阵；
步骤三：确定声学空间的耦合能量，得到总体耦合刚度；
步骤四：对所有N个子空间的质量矩阵、刚度矩阵以及总体耦合刚度矩阵分别进行整合排列，形成声学空间的总体质量矩阵和刚度矩阵，得到声学空间的特性方程；
步骤五：根据声学空间的特性方程，得到声学预报信息，声学预报信息包括声学空间的频率和模态。
本发明一种复杂形状声学空间的分析预报方法，还可以包括：
1、确定声学空间的耦合能量为：
其中ρ0为声场介质密度，pl(pm)为的子空间l(m)的声压，n为外法线方向，Slm为子空间l和m间耦合面的面积，ω为系统圆频率，令耦合声能量表达式Wc对声压函数归一化坐标每一个分量Pl的导数等于零，l＝1,2,...,N，得到整个声学空间的总的耦合刚度矩阵Kc。
2、声学空间的特性方程为：
(Kz-ω2Mz)P＝0
其中Kz和Mz分别为整个声学空间的刚度矩阵和质量矩阵，P为全局归一化坐标向量，其可以表示成P＝[P1,P2,...,PN]T。
有益效果：
本发明提出了一种复杂形状声学空间的建模分析方法，用于解决复杂形状声场的声学特性问题。其核心是将复杂形状声学空间分解为子空间声场，并构建子空间的声场耦合公式，依据变分原理求得耦合声场的声学特性参数。本发明没有上述方法中的假设限制，具有更广泛的适用范围。目前国内外还没有基于声场分解的复杂形状声场分析方法的技术出现及相关的文献报道。
本发明突破了局限于简单规则声场的约束，适用于具有复杂形状的声学空间声学问题的预报分析。与有限元方法相比，本发明的复杂声学空间分析预报方法物理意义明晰，且此方法无需任何网格划分并有着收敛快速，计算所需资源少等优点。有限元方法作为一种成熟技术，在结构、声学等多个领域已经得到广泛的应用。由于计算模型的大型化和复杂化，以及计算机计算能力的限制，导致有限元方法存在着计算误差以及计算效率低下等弊端，对于某些声学问题的中高频结果甚至无法求解。通过本发明提出的复杂形状声学空间分析预报方法，既可以解决复杂形状声学空间声学问题，又可以节省大量的计算成本。
附图说明
图1为本发明的流程示意图；
图2是本发明实施方式1的复杂声学空间示意图；
图3为复杂声学空间系统固有频率求解结果表。
具体实施方式
下面将结合附图对本发明做进一步详细说明。
本发明为复杂形状声学空间分析方法，包括以下步骤：
(1)根据复杂声学空间对象的几何特点，将其划分为N个子空间，N>2，对于一般复杂声学空间其预报精度有赖于所划分子空间个数，根据实际情况复杂声学空间划分得越为精细，其预报精度越高，故取N>2(N至少大于2)，相邻的两子空间通过耦合面耦合；
(2)对每一个子空间选择声压函数，依据变分法原理将声压表达式中的归一化向量坐标作为未知量，并形成对应于每一个子空间的刚度矩阵和质量矩阵；
(3)对每相邻的两个子空间之间运用声场耦合公式，将耦合能量转化为总体耦合刚度矩阵；
(4)对子空间的质量矩阵、刚度矩阵以及总体耦合刚度矩阵分别进行排列，形成整个声学空间的总体质量矩阵和刚度矩阵；
(5)得到整个复杂形状声学空间的特性方程，建立基于matlab的求解器即可输出此复杂形状声学空间的频率和模态等声学预报信息。
结合图1，详细介绍本发明的步骤：
1、建立求解域，根据声学空间的几何特点将其划分为多个子空间。对子空间编号，并设相邻的两子空间通过耦合界面耦合；
2、选择子空间声压函数表达式，并对声学子空间进行满足声压函数要求的坐标变换；
3、在变换后的子空间中，依据变分法原理建立子空间拉格朗日泛函；
4、以声压函数表达式中的归一化坐标为未知量，依据能量原理计算得到对应于每一个子空间的刚度矩阵和质量矩阵；
5、在耦合面上，建立声场耦合的能量表达式；
6、依据建立的声场耦合公式，按照能量原理计算得到整个声学空间的总体耦合刚度矩阵；
7、按照全局归一化坐标，对对应于每一个子空间的刚度矩阵和质量矩阵以及对应于整个声学空间的总体耦合刚度矩阵进行排列，形成整个声学空间的总体质量矩阵和刚度矩阵；
8、求解得到的线性方程组，即声学空间的特征方程，得到每一个子空间声压表达式中的归一化坐标；
9、通过计算结果得到其他信息，如固有频率、声压模态，也可进一步进行外部激励作用下的声压响应计算。
结合图2，复杂声学空间划分为N(N＝5)个子声场，分别由腔1，腔2，腔3，腔4和腔5表示。
依据声场特性,假设第s个子空间声压函数在其变换后的局部坐标系内表达式如下：
其中Ps为声压表达式的归一化坐标向量，为声压函数的容许函数向量，可以由正交多项式级数构成。
对每一个子空间，其拉格朗日泛函表达式为Ls＝Us-Ts-Ws，其中Us为子空间s的声场势能，Ts为子空间s的声场动能，Ws为子空间s边界壁面上做的功。在频率为ω时，子空间的声场势能和声场动能表达式如下：
其中ρ0为声场介质密度，c0为声速，Vs为子空间s体积，为梯度算子。
在边界面上，表面声压和空气质点振速在表面法向的分量un有如下关系：
于是在这些边界面上做的总功为：
其中Sl为s子空间的某一边界面面积。如果这些边界为声学刚性边界，那么有从而有Ws＝0。
令拉格朗日泛函对声压函数的归一化坐标取极值，可以得到声学子空间的线性特性方程组，以子空间s为例，其可进一步写为矩阵表达形式：
(Ks-ω2Ms)Ps＝0。 (6)
其中Ks为子空间s的刚度矩阵，Ms为子空间s的质量矩阵。
在耦合界面上，相邻的两个子声场p-和p+通过耦合界面耦合。假定p-为正向，于是耦合界面区域受到的总声压为：
Δp＝p--p+。 (7)
相邻的两声场通过耦合面的能量交换可以写为：
式(8)即为耦合面上的声场耦合能量表达式。由式(8)可以得到在声学空间的所有耦合面上，其总的耦合能量表达式Wc为：
其中pl(pm)为的子空间l(m)的声压，n为外法线方向，Slm为子空间l和m之间耦合面的面积。
令耦合能量表达式对声压函数的归一化坐标取极值，得到总体耦合刚度矩阵Kc。
按照全局归一化坐标，对子空间的刚度矩阵和质量矩阵以及整个声学空间的总体耦合刚度矩阵进行排列，得到总声学空间的特性方程，以总质量矩阵和刚度矩阵形式表示为：
(Kz-ω2Mz)P＝0。 (10)
这里Kz和Mz分别为整个声学空间的刚度矩阵和质量矩阵，其具体形式为：
P为全局归一化坐标向量，其可以表示成P＝[P1,P2,...,PN]T。
求解式(10)即可得到相应声学空间的声学特性预报结果。
根据以上求解步骤，对一个复杂的声学空间进行了验算，如图2所示。这一声学空间由具有倾斜平面和任意曲面的多个子声场构成，因而具有代表性。根据此声学空间的几何形状特点，可以将其划分为5个子空间，相邻子空间之间通过耦合界面联系在一起。图中给出了声场的几何参数，其具体值为：Lx0＝0.5m，Lx1＝Lx2＝Lx4＝Lx5＝1m，Lx3＝2m，Ly0＝1m，Ly1＝Ly5＝Ly6＝1.5m，Ly2＝Ly4＝2m，Ly3＝2.5m，Lz＝2.5m。其中，腔3的曲面函数为f(x,y)＝-(x-3)2/2+2.5。空气密度为ρ0＝1.21kg/m3，声速为c0＝340m/s。取每个子空间声场声压函数三个方向的截断级数分别为7-7-7，9-9-9和11-11-11。
应用本发明预报的系统固有频率结果与Ansys Acoustic30单元仿真结果比较如图3所示。通过本发明在以上例子的实施结果可以看出，本发明可用于耦合声场问题的求解预报。新发明的预报方法求解复杂声学空间声学特性问题的结果与Ansys仿真结果吻合良好，验证了本发明在复杂声学空间声学特性问题预报的正确性。另外，从表中可以看到，本发明预报结具有较快的收敛速度。
Claims (2)
1.一种复杂形状声学空间的分析预报方法，其特征在于：包括以下步骤，
步骤一：将复杂形状声学空间对象划分为N个子空间，N>2，每相邻的两个子空间通过耦合面耦合；
步骤二：选择每一个子空间的声压函数，确定每一个子空间的刚度矩阵和质量矩阵；
步骤三：确定声学空间的耦合能量，得到总体耦合刚度；
步骤四：对所有N个子空间的质量矩阵、刚度矩阵以及总体耦合刚度矩阵分别进行整合排列，形成声学空间的总体质量矩阵和刚度矩阵，得到声学空间的特性方程；
步骤五：根据声学空间的特性方程，得到声学预报信息，声学预报信息包括声学空间的频率和模态；
所述的确定声学空间的耦合能量为：
其中ρ0为声场介质密度，pl为子空间l的声压，pm为子空间m的声压，n为外法线方向，Slm为子空间l和m间耦合面的面积，ω为系统圆频率，令耦合声能量表达式Wc对声压函数归一化坐标每一个分量Pl的导数等于零，l＝1,2,...,N，得到整个声学空间的总的耦合刚度矩阵Kc。
2.根据权利要求1所述的一种复杂形状声学空间的分析预报方法，其特征在于：所述的声学空间的特性方程为：
(Kz-ω2Mz)P＝0
其中Kz和Mz分别为整个声学空间的刚度矩阵和质量矩阵，P为全局归一化坐标向量，其表示成P＝[P1,P2,...,PN]T。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510175544.6A CN104778321B (zh) | 2015-04-14 | 2015-04-14 | 一种复杂形状声学空间的分析预报方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510175544.6A CN104778321B (zh) | 2015-04-14 | 2015-04-14 | 一种复杂形状声学空间的分析预报方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104778321A CN104778321A (zh) | 2015-07-15 |
| CN104778321B true CN104778321B (zh) | 2017-11-17 |
Family
ID=53619782
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510175544.6A Expired - Fee Related CN104778321B (zh) | 2015-04-14 | 2015-04-14 | 一种复杂形状声学空间的分析预报方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104778321B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109883532B (zh) * | 2019-03-14 | 2021-04-30 | 哈尔滨工程大学 | 一种声源识别与声场预报方法 |
| CN113641954B (zh) * | 2021-07-20 | 2022-04-05 | 中国科学院声学研究所 | 一种复杂海洋环境中的三维声场快速预报方法及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1472514A (zh) * | 2003-06-19 | 2004-02-04 | 上海交通大学 | 快速噪声诊断的声场重构方法 |
| JP2004333224A (ja) * | 2003-05-02 | 2004-11-25 | Shimizu Corp | ３次元音響波動解析方法 |
| CN101902680A (zh) * | 2010-07-05 | 2010-12-01 | 南京大学 | 基于有源控制的虚拟耦合空间的设计方法 |
| CN104036140A (zh) * | 2014-06-13 | 2014-09-10 | 中国医学科学院生物医学工程研究所 | 一种用于声学不均匀媒介的磁声耦合成像声压求解方法 |
-
2015
- 2015-04-14 CN CN201510175544.6A patent/CN104778321B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004333224A (ja) * | 2003-05-02 | 2004-11-25 | Shimizu Corp | ３次元音響波動解析方法 |
| CN1472514A (zh) * | 2003-06-19 | 2004-02-04 | 上海交通大学 | 快速噪声诊断的声场重构方法 |
| CN101902680A (zh) * | 2010-07-05 | 2010-12-01 | 南京大学 | 基于有源控制的虚拟耦合空间的设计方法 |
| CN104036140A (zh) * | 2014-06-13 | 2014-09-10 | 中国医学科学院生物医学工程研究所 | 一种用于声学不均匀媒介的磁声耦合成像声压求解方法 |
Non-Patent Citations (2)
| Title |
|---|
| 声场与结构耦合系统的模态分析与灵敏度计算;高剑;《中国优秀博硕士学位论文全文数据库 (硕士) 工程科技Ⅱ辑》;20021215(第2期);11-12 * |
| 轿车车室内噪声的仿真分析;马天飞等;《数字化设计》;20041231;48-51 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104778321A (zh) | 2015-07-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Song | The scaled boundary finite element method in structural dynamics | |
| CN110457790B (zh) | 用于结构变形分析的近场动力学非连续伽辽金有限元方法 | |
| He et al. | Coupled analysis of 3D structural–acoustic problems using the edge-based smoothed finite element method/finite element method | |
| An et al. | Reduced-order extrapolation spectral-finite difference scheme based on POD method and error estimation for three-dimensional parabolic equation | |
| CN111859645B (zh) | 冲击波求解的改进musl格式物质点法 | |
| CN104778321B (zh) | 一种复杂形状声学空间的分析预报方法 | |
| Su | Accurate and robust adaptive mesh refinement for aerodynamic simulation with multi‐block structured curvilinear mesh | |
| Chen et al. | Localized method of fundamental solutions for acoustic analysis inside a car cavity with sound-absorbing material | |
| CN104950261B (zh) | 电池的硬件在环仿真测试方法及系统 | |
| Hejranfar et al. | Chebyshev collocation spectral lattice Boltzmann method for simulation of low-speed flows | |
| Tang et al. | A novel four-node quadrilateral element with continuous nodal stress | |
| CN109101752B (zh) | 一种复杂水工建筑物局部结构自振频率计算方法 | |
| Li et al. | Research on the iterative method for model updating based on the frequency response function | |
| CN112231863B (zh) | 太阳翼电池阵基板建模方法、装置、设备及存储介质 | |
| Ezzatneshan et al. | Simulation of three‐dimensional incompressible flows in generalized curvilinear coordinates using a high‐order compact finite‐difference lattice Boltzmann method | |
| CN109948253A (zh) | 薄板无网格Galerkin结构模态分析的GPU加速方法 | |
| Ngoc Nguyen et al. | A consecutive‐interpolation polyhedral finite element method for solid structures | |
| CN102938017B (zh) | 基于无网格模型计算周期结构板声学散射系数的方法 | |
| Chen et al. | A domain decomposition method for analyzing a coupling between multiple acoustical spaces (L) | |
| Okaue et al. | Study on analysis method for ferrofluid | |
| Shi et al. | Analysis of acoustic characteristics of arbitrary triangular prism and quadrangular prism acoustic cavities | |
| CN107169207A (zh) | 基于有限元理论的反应堆中子噪声频谱计算方法 | |
| Luo | A finite volume method based on weno reconstruction for compressible flows on hybrid grids | |
| Xie et al. | A fast direct solver for non-intrusive reduced order modeling of vibroacoustic problems | |
| Zhang et al. | An implicit algorithm for high-order DG/FV schemes for compressible flows on 2D arbitrary grids |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171117 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |