发明名称 快速噪声诊断的声场重构方法
摘要 快速噪声诊断的声场重构方法是在传统HELS方法和无限元方法的基础上,根据声源形状已知的某声场中,已知若干个点(称为测量点)的声压,把声场划分为若干个区域,把每个区域的声场表示成区域节点声压的线性叠加。根据测量点的声压,通过最小二乘法和奇异值分解,可以得到区域节点的声压。然后通过区域节点的声压,可以重构出整个声场的声压。该方法比HELS方法更准确、更高效地描述声场,可以近似求解各种形状的声源、各个频段声场重构问题,并且避免了使用椭球声无限单元的繁琐的方程积分,因而大大减少了计算用时。本方法既保留了传统的HELS方法的高效率,又大大地拓宽了该方法的计算适用领域。该方法同样能方便地推广到磁场、波浪冲击等无界领域。
申请公布号 CN1472514A 申请公布日期 2004.02.04
申请号 CN03129381.6 申请日期 2003.06.19
申请人 上海交通大学 发明人 杨瑞梁;汪鸿振
分类号 G01H17/00 主分类号 G01H17/00
代理机构 上海交达专利事务所 代理人 周文娟
主权项 1、一种快速噪声诊断的声场重构方法,其特征在于:该方法是根据声源形状已知的某声场中,已知若干个点(称为测量点)的声压,把声场划分为若干个区域,把每个区域的声场表示成区域节点声压的线性叠加,根据测量点的声压,通过最小二乘法和奇异值分解,可以得到区域节点的声压,然后通过区域节点的声压,可以重构出整个声场的声压,具体方法是:(1)建立椭球坐标系在声源表面附近画一个与声源形状比较接近的椭球,以此为基础建立一个椭球坐标系;(2)划分区域和节点根据声场中测量点的数目和分布情况,沿椭球坐标系的坐标轴方向把声场分为若干个区域,如果是外部声场,就沿椭球坐标系的半径方向划分为若干个无限区域,每个无限区域沿半径方向向外无限延伸;如果是内部声场,就划分为有限区域;(3)计算每个测量点的插值函数利用无限元法区域的插值函数公式来计算:<math> <mrow> <mi>&psi;</mi> <mo>=</mo> <msubsup> <mi>&psi;</mi> <mi>&mu;</mi> <mrow> <mo>(</mo> <mi>r</mi> <mo>)</mo> </mrow> </msubsup> <msubsup> <mi>&psi;</mi> <mi>v</mi> <mrow> <mo>(</mo> <mi>&phi;</mi> <mo>)</mo> </mrow> </msubsup> <msubsup> <mi>&psi;</mi> <mi>w</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mo>)</mo> </mrow> </msubsup> <mo>;</mo> </mrow> </math> (4)计算出节点声压根据测量点的声压值以及测量点的插值函数,通过最小二乘法和奇异值分解来计算节点声压: {C}=[ζ]u{ξ};(5)计算出重构点的插值函数重构点插值函数的计算方法和3一样;(6)重构出声场中任意一点的声压根据重构点的插值函数和节点的声压,通过下述声压方程来计算:<math> <mrow> <mi>P</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>C</mi> <mi>i</mi> </msub> <msub> <mi>&psi;</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>.</mo> </mrow> </math>
地址 200030上海市华山路1954号