发明名称 采用多频调制激光加热与光热信息重建技术的半透明材料辐射特性测量方法
摘要 采用多频调制激光加热与光热信息重建技术的半透明材料辐射特性测量方法,涉及一种测量半透明材料辐射特性参数的方法。其步骤为:分别采用不同频率的激光照射半透明材料的一侧表面,使用光电探测器分别测量在该入射激光频率下材料激光入射侧的频域半球反射辐射信号和激光出射侧的频域半球透射辐射信号,同时利用热电偶测温仪测量并记录材料两表面温度随时间的变化。根据频域半球反射辐射信号、频域半球透射辐射信号以及两表面随时间变化的温度,通过逆问题求解技术获得半透明材料的衰减系数和反照率。本发明通过建立测量半透明材料衰减系数和反照率的正、逆问题模型,能简单、快速、准确的利用逆问题求解技术同时测量半透明材料衰减系数和反照率。
申请公布号 CN103528963A 申请公布日期 2014.01.22
申请号 CN201310533552.4 申请日期 2013.11.01
申请人 哈尔滨工业大学 发明人 齐宏;牛春洋;宫帅;阮立明
分类号 G01N21/17(2006.01)I 主分类号 G01N21/17(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 1.采用多频调制激光加热与光热信息重建技术的半透明材料辐射特性测量方法,其特征在于,该方法的具体步骤为:步骤一、将待测半透明材料制作成厚度为L的扁平试件,调节激光器位置,使其对准试件左侧中心位置,并且激光入射方向垂直于试件左侧表面,同时分别将两只热电偶探头固定在试件左右两侧表面上;步骤二、设定振幅可调的正弦激光光源的激光频率为ω<sub>1</sub>,打开激光光源对试件左侧表面进行加热,与此同时使用光电探测器分别测量试件激光入射侧的频域半球反射辐射信号<img file="FDA0000407021270000011.GIF" wi="214" he="90" />和激光出射侧的频域半球透射辐射信号<img file="FDA0000407021270000012.GIF" wi="232" he="82" />并且同步使用热电偶测温仪同步测量得到左侧激光入射表面温度随时间的变化为T<sub>w1</sub>(t),右侧表面温度随时间的变化为T<sub>w2</sub>(t);步骤三、依次改变激光光源频率为ω<sub>i</sub>(i=2,…,m),采用步骤二中方法分别测量得到相应入射激光频率下的试件激光入射侧的频域半球反射辐射信号<img file="FDA0000407021270000013.GIF" wi="204" he="90" />和激光出射侧的频域半球透射辐射信号<img file="FDA0000407021270000014.GIF" wi="221" he="82" />步骤四、根据逆问题算法假设待测材料的衰减系数β和反照率<img file="FDA0000407021270000015.GIF" wi="68" he="59" />通过分别对频域辐射传输方程(1)和瞬态导热辐射耦合方程组(2)的求解,得到在假设衰减系数β和反照率<img file="FDA0000407021270000016.GIF" wi="42" he="59" />条件下的材料激光入射侧的频域半球反射辐射信号估计值<img file="FDA0000407021270000017.GIF" wi="206" he="92" />和激光出射侧的频域半球透射辐射信号估计值<img file="FDA0000407021270000018.GIF" wi="226" he="84" />以及两侧面边界温度信号估计值T′<sub>w1</sub>(t)和T′<sub>w2</sub>(t);<maths num="0001"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mover><mi>I</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>&theta;</mi><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><mi>x</mi></mrow></mfrac><mo>=</mo><mo>-</mo><mi>&beta;</mi><msub><mover><mi>I</mi><mo>^</mo></mover><mi>c</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>&theta;</mi><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>+</mo><mfrac><msub><mi>&kappa;</mi><mi>s</mi></msub><mn>2</mn></mfrac><msubsup><mo>&Integral;</mo><mn>0</mn><mi>&pi;</mi></msubsup><mover><mi>I</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>x</mi><mo>,</mo><msup><mi>&theta;</mi><mo>&prime;</mo></msup><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow><mi>&Phi;</mi><mrow><mo>(</mo><msup><mi>&theta;</mi><mo>&prime;</mo></msup><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>sin</mi><msup><mi>&theta;</mi><mo>&prime;</mo></msup><mi>d</mi><msup><mi>&theta;</mi><mo>&prime;</mo></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,<img file="FDA00004070212700000110.GIF" wi="201" he="83" />为x处的频域辐射强度,<img file="FDA00004070212700000111.GIF" wi="176" he="85" />是入射的频域辐射强度,x为距离,θ辐射方向,β为衰减系数,ω为角频率,κ<sub>s</sub>为散射系数,Φ(θ′,θ)为散射相函数;<maths num="0002"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mrow><mo>&PartialD;</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>&Omega;</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><mi>x</mi></mrow></mfrac><mo>=</mo><mo>-</mo><mi>&beta;I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>&Omega;</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&kappa;</mi><mi>a</mi></msub><msub><mi>I</mi><mi>b</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>+</mo><mfrac><msub><mi>&kappa;</mi><mi>s</mi></msub><mrow><mn>4</mn><mi>&pi;</mi></mrow></mfrac><msup><msub><mo>&Integral;</mo><mrow><mn>4</mn><mi>&pi;</mi></mrow></msub><mi>I</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><msup><mi>&Omega;</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><mi>&Phi;</mi><mrow><mo>(</mo><msup><mi>&Omega;</mi><mo>&prime;</mo></msup><mo>,</mo><mi>&Omega;</mi><mo>)</mo></mrow><mi>d</mi><msup><mi>&Omega;</mi><mo>&prime;</mo></msup></mtd></mtr><mtr><mtd><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><mo>=</mo><mi>&lambda;</mi><mfrac><mrow><msup><mo>&PartialD;</mo><mn>2</mn></msup><mi>T</mi></mrow><mrow><mo>&PartialD;</mo><msup><mi>x</mi><mn>2</mn></msup></mrow></mfrac><mo>-</mo><mfrac><mrow><mo>&PartialD;</mo><msub><mi>q</mi><mi>r</mi></msub></mrow><mrow><mo>&PartialD;</mo><mi>x</mi></mrow></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,I为辐射强度;x为距离;Ω为立体角;β为衰减系数;κ<sub>a</sub>为吸收系数;κ<sub>s</sub>为散射系数;Φ(Ω′,Ω)为散射相函数,T为温度;ρ为材料密度;c<sub>p</sub>为材料比定压热容;λ为导热系数;q<sub>r</sub>为辐射项热流;步骤五、将步骤二和步骤三得到的试件激光入射侧的频域半球反射辐射信号测量值<img file="FDA0000407021270000021.GIF" wi="208" he="90" />和激光出射侧的频域半球透射辐射信号测量值<img file="FDA0000407021270000022.GIF" wi="229" he="82" />试件左边界温度测量值T<sub>w1</sub>(t),试件右边界温度测量值T<sub>w2</sub>(t)和步骤四得到的激光入射侧的频域半球反射辐射信号估计值<img file="FDA0000407021270000023.GIF" wi="210" he="92" />和激光出射侧的频域半球透射辐射信号估计值<img file="FDA0000407021270000024.GIF" wi="233" he="84" />以及左边界温度估计值T′<sub>w1</sub>(t),右边界温度估计值T′<sub>w2</sub>(t)代入如下目标函数计算公式,计算得到目标函数值F<sub>obj</sub>;<img file="FDA0000407021270000025.GIF" wi="1811" he="436" />步骤六、判断步骤五中的目标函数值F<sub>obj</sub>是否小于设定阈值ξ,若是,则步骤四中所假设的半透明材料的衰减系数β和反照率<img file="FDA0000407021270000026.GIF" wi="52" he="59" />即为该待测材料的真实衰减系数和反照率;若不是,则返回步骤四,重新根据逆问题算法假设材料的衰减系数β和反照率<img file="FDA0000407021270000027.GIF" wi="46" he="59" />的设定值,并重新计算,直至步骤五中的目标函数值F<sub>obj</sub>小于设定阈值ξ,完成采用多频调制激光加热与光热信息重建技术的半透明材料衰减系数和反照率的测量。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号
您可能感兴趣的专利