发明名称 超声波探测物体内部瞬态非均匀温度场的灵敏度方法
摘要 本发明公开了一种超声波探测物体内部瞬态非均匀温度场的灵敏度方法,根据标定不同温度条件下被测试件母材内超声波传播时间                                               <img file="dest_path_image002.GIF" wi="8" he="16" />与温度<img file="dest_path_image004.GIF" wi="16" he="18" />的相关关系以及由超声波脉冲回波法获得的测量全过程时间及对应时刻下的超声波传播时间,求解基于解灵敏度方程的热传导反问题,要求优化所得等效的温度边界条件使得数值模拟所得的超声在物体内传播时间与实测的超声波传播时间差值最小,再通过热传导正问题的求解即可获得试件内部不同时刻的温度场分布状态。实现快速准确地获得金属或金属合金物体表面温度和内部非均匀温度场的优点。实现对金属或合金物体表面温度和温度内部瞬态非均匀温度场的实时测量。
申请公布号 CN104596667A 申请公布日期 2015.05.06
申请号 CN201510003663.3 申请日期 2015.01.05
申请人 中国空气动力研究与发展中心计算空气动力研究所 发明人 魏东;石友安;杜雁霞;曾磊;桂业伟;肖光明;刘磊;杨肖锋
分类号 G01K11/22(2006.01)I 主分类号 G01K11/22(2006.01)I
代理机构 北京中恒高博知识产权代理有限公司 11249 代理人 宋敏
主权项 一种超声波探测物体内部瞬态非均匀温度场的灵敏度方法,其特征在于,包括以下步骤:步骤一:加工一个与被测物体母材相同的长方体试件,该长方体试件被测方向的长度大于一个超声波波长;并标定不同温度条件下长方体试件内超声波传播时间t与温度T的相关关系,由最小二乘法拟合出其线性或二次函数关系式T(t);步骤二:使用耐高温超声波探头,由超声波脉冲回波法,获得上述被测物体母材t<sub>i</sub>时刻的超声波传播时间t<sub>i,exp</sub>,并同步记录测量全过程时间及对应时刻下的超声波传播时间;步骤三:基于上述步骤一中的线性或二次函数关系式T(t),以及步骤二的超声波传播时间t<sub>i,exp</sub>和超声波传播时间,求解热传导的反问题:热传导反问题的目标函数为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>J</mi><mrow><mo>(</mo><msub><mi>T</mi><mi>equivalent</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><mi>&epsiv;</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>{</mo><msub><mi>t</mi><mrow><mi>i</mi><mo>,</mo><mi>cal</mi></mrow></msub><mo>-</mo><msub><mi>t</mi><mrow><mi>i</mi><mo>,</mo><mi>exp</mi></mrow></msub><mo>}</mo></mrow><mn>2</mn></msup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>{</mo><mn>2</mn><msubsup><mo>&Integral;</mo><mn>0</mn><mi>L</mi></msubsup><mfrac><mn>1</mn><mrow><mi>V</mi><mo>[</mo><mi>T</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow></mfrac><mi>dx</mi><mo>-</mo><msub><mi>t</mi><mrow><mi>i</mi><mo>,</mo><mi>exp</mi></mrow></msub><mo>}</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0000651806210000011.GIF" wi="1887" he="181" /></maths>式中,T<sub>equivalent</sub>为待辨识的被测物体等效温度边界条件;ε是实际测量值与计算值之间的误差;t<sub>i,cal</sub>为数值计算得到的t<sub>i</sub>时刻的超声波传播时间,下标i表示的测量时间序数,n表示总的测量时间点数;目标函数约束条件为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>S</mi><mo>.</mo><mi>t</mi><mo>.</mo><mfrac><mo>&PartialD;</mo><mrow><mo>&PartialD;</mo><mi>x</mi></mrow></mfrac><mo>[</mo><mi>k</mi><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><mfrac><mrow><mo>&PartialD;</mo><mi>T</mi></mrow><mrow><mo>&PartialD;</mo><mi>x</mi></mrow></mfrac><mo>]</mo><mo>=</mo><mi>&rho;</mi><msub><mi>C</mi><mi>p</mi></msub><mfrac><mrow><mo>&PartialD;</mo><mi>T</mi></mrow><mrow><mo>&PartialD;</mo><mi>t</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000651806210000012.GIF" wi="784" he="202" /></maths>T(x,t),t>0,x∈[0,L]式中,ρ为材料密度;C<sub>p</sub>为材料比热;k=k(T)为材料的热传导系数,是温度的函数;L为试件被测方向的长度;根据热传导反问题的目标函数和约束条件得出等效的温度边界条件T<sub>equivalent</sub>;步骤四:根据步骤三所得的等效的温度边界条件T<sub>equivalent</sub>,基于热传导的正问题求解获得试件内部不同时刻的温度场分布状态。
地址 621000 四川省绵阳市二环路南段6号