发明名称 弹性固体中声速测量方法
摘要 本发明的目的在于提供弹性固体中声速测量方法。将矩形弹性板通过棉线与旋转装置相连接,收发换能器通过连接杆与固定装置相连接,矩形弹性板和收发换能器均位于水中;以声波垂直入射弹性板正横方向为起始位置,收发换能器发射脉冲信号,接收并测量弹性板的回波信号,通过测量的所有旋转角度的回波信号,得到弹性散射波出现的角度,从而得到弹性板中Lamb波传播的相速度,建立水中弹性板Lamb波传播模型,得到频散方程,将相速度分别代入水中弹性板的对称和反对称频散方程,得到弹性固体中的纵波和横波波速。本发明可以同时测量得到固体中的横波和纵波波速,不需分别测量,提高了测量效率。
申请公布号 CN101915602B 申请公布日期 2011.09.14
申请号 CN201010226607.3 申请日期 2010.07.15
申请人 哈尔滨工程大学 发明人 孙辉;朱广平;张明辉;陈文剑;孙筱逸
分类号 G01H5/00(2006.01)I;G01N29/07(2006.01)I 主分类号 G01H5/00(2006.01)I
代理机构 代理人
主权项 1.弹性固体中声速测量方法,其特征是:(1)将矩形弹性板通过棉线与旋转装置相连接,收发换能器通过连接杆与固定装置相连接,矩形弹性板和收发换能器均位于水中,矩形弹性板上边缘与水面平行,且弹性板中心与收发换能器中心处在相同深度;(2)确定测量的起始位置:以声波垂直入射弹性板正横方向为起始位置;(3)收发换能器发射CW脉冲信号,接收并测量弹性板的回波信号,以0.1度为步长旋转弹性板,待旋转弹性板在水中静止后,重复所述的测量,直至完成旋转弹性板旋转90度结束测量;(4)通过测量的所有旋转角度的回波信号,得到弹性散射波出现的角度θ,通过公式sinθ=c<sub>1</sub>/c<sub>p</sub>计算得到弹性板中Lamb波传播的相速度c<sub>p</sub>,其中c<sub>1</sub>为水中声速;(5)建立水中弹性板Lamb波传播模型,得到频散方程为:对称模态<maths num="0001"><![CDATA[<math><mrow><mo>{</mo><msup><mrow><mo>(</mo><msup><mi>&beta;</mi><mn>2</mn></msup><mo>-</mo><msup><mi>&xi;</mi><mn>2</mn></msup><mo>)</mo></mrow><mn>2</mn></msup><mi>sin</mi><mrow><mo>(</mo><mi>&beta;h</mi><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mi>&alpha;h</mi><mo>)</mo></mrow><mo>+</mo><mn>4</mn><msup><mi>&xi;</mi><mn>2</mn></msup><mi>&alpha;&beta;</mi><mi>sin</mi><mrow><mo>(</mo><mi>&alpha;h</mi><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mi>&beta;h</mi><mo>)</mo></mrow><mo>}</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mo>-</mo><mi>i</mi><mfrac><mrow><msub><mi>&rho;</mi><mn>1</mn></msub><msup><mi>&omega;</mi><mn>2</mn></msup><mi>&alpha;</mi></mrow><mrow><mi>&rho;&eta;</mi><msubsup><mi>c</mi><mi>T</mi><mn>2</mn></msubsup></mrow></mfrac><mrow><mo>{</mo><msup><mrow><mo>(</mo><msup><mi>&beta;</mi><mn>2</mn></msup><mo>+</mo><msup><mi>&xi;</mi><mn>2</mn></msup><mo>)</mo></mrow><mn>2</mn></msup><mi></mi><mi>sin</mi><mrow><mo>(</mo><mi>&alpha;h</mi><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><mi>&beta;h</mi><mo>)</mo></mrow><mo>}</mo></mrow><mo>=</mo><mn>0</mn><mo>,</mo></mrow></math>]]></maths>反对称模态<maths num="0003"><![CDATA[<math><mrow><mo>{</mo><msup><mrow><mo>(</mo><msup><mi>&beta;</mi><mn>2</mn></msup><mo>-</mo><msup><mi>&xi;</mi><mn>2</mn></msup><mo>)</mo></mrow><mn>2</mn></msup><mi>sin</mi><mrow><mo>(</mo><mi>&alpha;h</mi><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mi>&beta;h</mi><mo>)</mo></mrow><mo>+</mo><mn>4</mn><msup><mi>&xi;</mi><mn>2</mn></msup><mi>&alpha;&beta;</mi><mi>cos</mi><mrow><mo>(</mo><mi>&alpha;h</mi><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><mi>&beta;h</mi><mo>)</mo></mrow><mo>}</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mo>+</mo><mi>i</mi><mfrac><mrow><msub><mi>&rho;</mi><mn>1</mn></msub><msup><mi>&omega;</mi><mn>2</mn></msup><mi>&alpha;</mi></mrow><mrow><mi>&rho;&eta;</mi><msubsup><mi>c</mi><mi>T</mi><mn>2</mn></msubsup></mrow></mfrac><mrow><mo>{</mo><msup><mrow><mo>(</mo><msup><mi>&beta;</mi><mn>2</mn></msup><mo>+</mo><msup><mi>&xi;</mi><mn>2</mn></msup><mo>)</mo></mrow><mn>2</mn></msup><mi>cos</mi><mrow><mo>(</mo><mi>&alpha;h</mi><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mi>&beta;h</mi><mo>)</mo></mrow><mo>}</mo></mrow><mo>=</mo><mn>0</mn><mo>,</mo></mrow></math>]]></maths>其中<img file="FSA00000190067900015.GIF" wi="296" he="129" /><img file="FSA00000190067900016.GIF" wi="297" he="128" /><img file="FSA00000190067900017.GIF" wi="288" he="129" /><img file="FSA00000190067900018.GIF" wi="166" he="115" />ω=2πf,ω为角频率,f为声波频率,ρ<sub>1</sub>为水的密度,c<sub>1</sub>为水中声速,ρ为弹性板密度,2h为板的厚度,c<sub>L</sub>和c<sub>T</sub>分别表示固体中的纵波和横波速度,c<sub>p</sub>为弹性板中Lamb波传播的相速度,<img file="FSA00000190067900019.GIF" wi="156" he="55" />为虚数单位;将Lamb波相速度c<sub>p</sub>分别代入水中弹性板的对称和反对称频散方程,得到一组非线性方程组,采用优化算法中的遗传算法对方程组进行求解,定义代价函数F为:<maths num="0005"><![CDATA[<math><mrow><mi>F</mi><mo>=</mo><mfrac><mn>1</mn><mi>M</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mfrac><mrow><msup><mrow><mo>(</mo><msubsup><mi>c</mi><mi>Ln</mi><mo>&prime;</mo></msubsup><mo>-</mo><mi>E</mi><mrow><mo>(</mo><msubsup><mi>c</mi><mi>L</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msubsup><mi>c</mi><mi>Tn</mi><mo>&prime;</mo></msubsup><mo>-</mo><mi>E</mi><mrow><mo>(</mo><msubsup><mi>c</mi><mi>T</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><mi>E</mi><msup><mrow><mo>(</mo><msubsup><mi>c</mi><mi>L</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>E</mi><msup><mrow><mo>(</mo><msubsup><mi>c</mi><mi>T</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow></math>]]></maths>其中M表示重复实验的次数,<img file="FSA00000190067900022.GIF" wi="55" he="58" />和<img file="FSA00000190067900023.GIF" wi="51" he="57" />表示进行第n次实验得到的纵波和横波波速,<img file="FSA00000190067900024.GIF" wi="114" he="74" />和<img file="FSA00000190067900025.GIF" wi="117" he="74" />表示M次实验结果的平均值,M的值大于5,当F接近0时,优化计算得到的<img file="FSA00000190067900026.GIF" wi="116" he="74" />和<img file="FSA00000190067900027.GIF" wi="115" he="75" />即分别为弹性固体中的纵波和横波波速。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室