发明名称 一种测量水声材料透声性能角谱的方法
摘要 本发明涉及一种测量水声材料透声性能角谱的方法。本发明在自由场或近自由场水域中,指向性声源发射宽带脉冲信号,通过测量水声材料试样表面及试样背面的复声压分布,利用时频傅立叶变换及二维空间傅立叶变换,分离出不同频率、不同入射角度样品表面的入射波分量和透射波分量,从而得到样品不同频率时声压透射系数角谱。本发明有益的效果是:一是实现水声材料构件透声性能的角谱测量;二是无需反复发射不同单频信号,无需进行样品多角度地旋转,就能够测试声波斜入射时,水声材料不同频率的透声性能分布,极大地提高了工作效率。
申请公布号 CN103675100A 申请公布日期 2014.03.26
申请号 CN201310711013.5 申请日期 2013.12.20
申请人 中国船舶重工集团公司第七一五研究所 发明人 易燕;李水
分类号 G01N29/04(2006.01)I 主分类号 G01N29/04(2006.01)I
代理机构 杭州九洲专利事务所有限公司 33101 代理人 陈健
主权项 一种测量水声材料透声性能角谱的方法,其特征在于:该方法通过三个测试面上复声压测试,分离出不同频率、不同入射角度样品表面的入射波分量和透射波分量,从而得到样品的声压透射系数角谱,具体包括如下步骤:步骤一:假定声源位于待测样品上方z=zs处,在声源与待测材料间及待测材料背面,z=z1、z=z2和z=z3设定为三个测量面,测量面与待测材料互相平行,并假定z=0平面与待测材料重合,且原点与材料几何中心重合;步骤二:声源发射宽带信号,利用二维傅立叶变换则z=z1、z=z2和z=z3三个测量面上的复声压分解为各个方向上的平面波分量: <mrow> <mi>P</mi> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mi>x</mi> </msub> <mo>,</mo> <msub> <mi>k</mi> <mi>y</mi> </msub> <mo>,</mo> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mo>&Integral;</mo> <munderover> <mo>&Integral;</mo> <mrow> <mo>-</mo> <mo>&infin;</mo> </mrow> <mo>&infin;</mo> </munderover> <mi>p</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mi>x</mi> </msub> <mi>x</mi> <mo>+</mo> <msub> <mi>k</mi> <mi>y</mi> </msub> <mi>y</mi> <mo>)</mo> </mrow> </mrow> </msup> <mi>dxdy</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>P</mi> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mi>x</mi> </msub> <mo>,</mo> <msub> <mi>k</mi> <mi>y</mi> </msub> <mo>,</mo> <msub> <mi>z</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mo>&Integral;</mo> <munderover> <mo>&Integral;</mo> <mrow> <mo>-</mo> <mo>&infin;</mo> </mrow> <mo>&infin;</mo> </munderover> <mi>p</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <msub> <mi>z</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mi>x</mi> </msub> <mi>x</mi> <mo>+</mo> <msub> <mi>k</mi> <mi>y</mi> </msub> <mi>y</mi> <mo>)</mo> </mrow> </mrow> </msup> <mi>dxdy</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>P</mi> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mi>x</mi> </msub> <mo>,</mo> <msub> <mi>k</mi> <mi>y</mi> </msub> <mo>,</mo> <msub> <mi>z</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mo>&Integral;</mo> <munderover> <mo>&Integral;</mo> <mrow> <mo>-</mo> <mo>&infin;</mo> </mrow> <mo>&infin;</mo> </munderover> <mi>p</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <msub> <mi>z</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mi>x</mi> </msub> <mi>x</mi> <mo>+</mo> <msub> <mi>k</mi> <mi>y</mi> </msub> <mi>y</mi> <mo>)</mo> </mrow> </mrow> </msup> <mi>dxdy</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>步骤三:z1、z2面上的各平面波分量又可以用被测样品平面上的不同频率入射波和反射波平面分量表示:P(kx,ky,z1)=Pi(kx,ky,0)exp(jkzz1)+Pr(kx,ky,0)exp(‑jkzz1)   (4)P(kx,ky,z2)=Pi(kx,ky,0)exp(jkzz2)+Pr(kx,ky,0)exp(‑jkzz2)   (5)通过式(4)和(5),可以分离出z=0面上入射平面波: <mrow> <msub> <mi>P</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mi>x</mi> </msub> <mo>,</mo> <msub> <mi>k</mi> <mi>y</mi> </msub> <mo>,</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>P</mi> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mi>x</mi> </msub> <mo>,</mo> <msub> <mi>k</mi> <mi>y</mi> </msub> <mo>,</mo> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <msub> <mi>k</mi> <mi>z</mi> </msub> <msub> <mi>z</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>P</mi> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mi>x</mi> </msub> <mo>,</mo> <msub> <mi>k</mi> <mi>y</mi> </msub> <mo>,</mo> <msub> <mi>z</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <msub> <mi>k</mi> <mi>z</mi> </msub> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mi>j</mi> <msub> <mi>k</mi> <mi>z</mi> </msub> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>-</mo> <mi>j</mi> <msub> <mi>k</mi> <mi>z</mi> </msub> <msub> <mi>z</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mi>j</mi> <msub> <mi>k</mi> <mi>z</mi> </msub> <msub> <mi>z</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>j</mi> <msub> <mi>k</mi> <mi>z</mi> </msub> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>步骤三:根据平面波传播原理,z=z3面上无样品时入射平面波为:Pi(kx,ky,z3)=Pi(kx,ky,0)exp(‑jkzz3)       (7)z=z3面上有样品时入射平面波为:Pt(kx,ky,z3)=P(kx,ky,z3)       (8)步骤四:由式(9)可以计算得到声波任意入射方向的声压透射系数: <mrow> <mfenced open='' close=''> <mtable> <mtr> <mtd> <msub> <mi>t</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mi>x</mi> </msub> <mo>,</mo> <msub> <mi>k</mi> <mi>y</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>P</mi> <mi>t</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mi>x</mi> </msub> <mo>,</mo> <msub> <mi>k</mi> <mi>y</mi> </msub> <mo>,</mo> <msub> <mi>z</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>P</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mi>x</mi> </msub> <mo>,</mo> <msub> <mi>k</mi> <mi>y</mi> </msub> <mo>,</mo> <msub> <mi>z</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mo>=</mo> <mfrac> <mrow> <mi>P</mi> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mi>x</mi> </msub> <mo>,</mo> <msub> <mi>k</mi> <mi>y</mi> </msub> <mo>,</mo> <msub> <mi>z</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <mrow> <mo>(</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mi>j</mi> <msub> <mi>k</mi> <mi>z</mi> </msub> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>-</mo> <mi>j</mi> <msub> <mi>k</mi> <mi>z</mi> </msub> <msub> <mi>z</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mi>j</mi> <msub> <mi>k</mi> <mi>z</mi> </msub> <msub> <mi>z</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>j</mi> <msub> <mi>k</mi> <mi>z</mi> </msub> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mrow> <mrow> <mi>P</mi> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mi>x</mi> </msub> <mo>,</mo> <msub> <mi>k</mi> <mi>y</mi> </msub> <mo>,</mo> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <msub> <mi>k</mi> <mi>z</mi> </msub> <msub> <mi>z</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>j</mi> <msub> <mi>k</mi> <mi>z</mi> </msub> <msub> <mi>z</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>P</mi> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mi>x</mi> </msub> <mo>,</mo> <msub> <mi>k</mi> <mi>y</mi> </msub> <mo>,</mo> <msub> <mi>z</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <msub> <mi>k</mi> <mi>z</mi> </msub> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>-</mo> <mi>j</mi> <msub> <mi>k</mi> <mi>z</mi> </msub> <msub> <mi>z</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow> 入射波的入射角度由下式确定: <mrow> <mi>&theta;</mi> <mo>=</mo> <msup> <mi>sin</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>[</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mi>k</mi> <mi>x</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>k</mi> <mi>y</mi> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> <mo>/</mo> <msub> <mi>k</mi> <mn>0</mn> </msub> <mo>]</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>式中 <mrow> <msub> <mi>k</mi> <mn>0</mn> </msub> <mo>=</mo> <msqrt> <msubsup> <mi>k</mi> <mi>x</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>k</mi> <mi>y</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>k</mi> <mi>z</mi> <mn>2</mn> </msubsup> </msqrt> <mo>=</mo> <mn>2</mn> <mi>&pi;f</mi> <mo>/</mo> <mi>c</mi> <mo>,</mo> </mrow>f为声波频率。
地址 310012 浙江省杭州市西湖区华星路96号