主权项 |
一种测量水声材料透声性能角谱的方法,其特征在于:该方法通过三个测试面上复声压测试,分离出不同频率、不同入射角度样品表面的入射波分量和透射波分量,从而得到样品的声压透射系数角谱,具体包括如下步骤:步骤一:假定声源位于待测样品上方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>∫</mo> <munderover> <mo>∫</mo> <mrow> <mo>-</mo> <mo>∞</mo> </mrow> <mo>∞</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>∫</mo> <munderover> <mo>∫</mo> <mrow> <mo>-</mo> <mo>∞</mo> </mrow> <mo>∞</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>∫</mo> <munderover> <mo>∫</mo> <mrow> <mo>-</mo> <mo>∞</mo> </mrow> <mo>∞</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>·</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>θ</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>πf</mi> <mo>/</mo> <mi>c</mi> <mo>,</mo> </mrow>f为声波频率。 |