发明名称 利用脉冲激光加热产生的瞬态光热信号测量半透明材料热物性参数的方法
摘要 利用脉冲激光加热产生的瞬态光热信号测量半透明材料热物性参数的方法,本发明涉及一种测量半透明材料热物性参数的方法。其步骤为:激光照射半透明材料的一侧表面,利用热电偶测温仪测量并记录材料两表面温度随时间的变化,同步使用激光功率计分别测量试件激光入射侧的半球反射辐射信号和激光出射侧的半球透射辐射信号。根据激光透射辐射信号、反射辐射信号以及两表面随时间变化的温度,通过逆问题求解技术获得半透明材料的吸收系数、散射系数和导热系数。本发明通过建立测量半透明材料吸收系数、散射系数和导热系数的正、逆问题模型,能简单、快速、准确的利用逆问题求解技术同时测量半透明材料吸收系数、散射系数和导热系数。
申请公布号 CN103528978A 申请公布日期 2014.01.22
申请号 CN201310533554.3 申请日期 2013.11.01
申请人 哈尔滨工业大学 发明人 齐宏;牛春洋;姚禹辰;阮立明
分类号 G01N21/31(2006.01)I;G01N21/49(2006.01)I;G01N25/20(2006.01)I 主分类号 G01N21/31(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 一种利用脉冲激光加热产生的瞬态光热信号测量半透明材料热物性参数的方法,其特征在于,该方法的具体步骤为:步骤一、将待测半透明材料制作成厚度为δ的扁平试件,将激光器对准试件左侧中心,并使激光方向垂直于试件表面,并且分别将两只连接到热电偶温度测量仪上的热电偶探头固定在试件左右两侧;步骤二、打开脉冲激光器对试件左侧表面进行一个脉冲的加热,与此同时使用热电偶测温仪同步测量得到左侧激光入射表面温度随时间的变化为Tw1(t),右侧表面温度随时间的变化为Tw2(t),并且同步使用激光功率计分别测量试件激光入射侧的半球反射辐射信号QR(t)和激光出射侧的半球透射辐射信号QT(t);其中测量时间间隔为dt,总测量时间为N倍激光脉冲宽度tp;步骤三、根据逆问题算法假设待测材料的吸收系数κa、散射系数κs和导热系数λ;通过对稳态辐射传输方程(1)和瞬态导热微分方程(2)的耦合求解,得到在假设吸收系数κa、散射系数κs和导热系数λ条件下的材料两侧面边界温度响应估计值T′w1(t)和T′w2(t),以及激光入射侧的半球反射辐射信号估计值Q′R(t)和激光出射侧的半球透射辐射信号估计值Q′T(t); <mrow> <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> <msub> <mo>&Integral;</mo> <mrow> <mn>4</mn> <mi>&pi;</mi> </mrow> </msub> <mi>I</mi> <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> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>式中,I为辐射强度;x为距离;Ω为立体角;β为衰减系数;κa为吸收系数;κs为散射系数;Φ(Ω′,Ω)为散射相函数; <mrow> <msub> <mi>&rho;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> <msup> <mrow> <mo>&PartialD;</mo> <mi>x</mi> </mrow> <mn>2</mn> </msup> </mfrac> <mo>-</mo> <mfrac> <msub> <mrow> <mo>&PartialD;</mo> <mi>q</mi> </mrow> <mi>r</mi> </msub> <mrow> <mo>&PartialD;</mo> <mi>x</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>式中,T为温度;ρ为材料密度;cp为材料比定压热容;λ为导热系数;qr为辐射项热流;步骤四、将步骤二得到的左边界温度测量值Tw1(t),右边界温度测量值Tw2(t),半球反射信号测量值QR(t)以及半球透射信号测量值QT(t)和步骤三得到的左边界温度估计值T′w1(t),右边界温度估计值T′w2(t),半球反射信号估计值Q′R(t)以及半球透射信号估计值 Q′T(t)代入如下目标函数计算公式,计算得到目标函数值Fobj; <mrow> <msub> <mi>F</mi> <mi>obj</mi> </msub> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mfenced open='{' close='}'> <mtable> <mtr> <mtd> <msup> <mrow> <mo>[</mo> <mfrac> <mrow> <msubsup> <mi>T</mi> <mrow> <mi>w</mi> <mn>1</mn> </mrow> <mo>&prime;</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>&kappa;</mi> <mi>a</mi> </msub> <mo>,</mo> <msub> <mi>&kappa;</mi> <mi>s</mi> </msub> <mo>,</mo> <mi>&lambda;</mi> <msub> <mrow> <mo>,</mo> <mi>t</mi> </mrow> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>T</mi> <mrow> <mi>w</mi> <mn>1</mn> </mrow> </msub> <mo>(</mo> <msub> <mi>&kappa;</mi> <mi>a</mi> </msub> <mo>,</mo> <msub> <mi>&kappa;</mi> <mi>s</mi> </msub> <mo>,</mo> <mi>&lambda;</mi> <msub> <mrow> <mo>,</mo> <mi>t</mi> </mrow> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mrow> <msub> <mi>T</mi> <mrow> <mi>w</mi> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>&kappa;</mi> <mi>a</mi> </msub> <mo>,</mo> <msub> <mi>&kappa;</mi> <mi>s</mi> </msub> <mo>,</mo> <mi>&lambda;</mi> <msub> <mrow> <mo>,</mo> <mi>t</mi> </mrow> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>]</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>[</mo> <mfrac> <mrow> <msubsup> <mi>T</mi> <mrow> <mi>w</mi> <mn>2</mn> </mrow> <mo>&prime;</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>&kappa;</mi> <mi>a</mi> </msub> <mo>,</mo> <msub> <mi>&kappa;</mi> <mi>s</mi> </msub> <mo>,</mo> <mi>&lambda;</mi> <msub> <mrow> <mo>,</mo> <mi>t</mi> </mrow> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>T</mi> <mrow> <mi>w</mi> <mn>2</mn> </mrow> </msub> <mo>(</mo> <msub> <mi>&kappa;</mi> <mi>a</mi> </msub> <mo>,</mo> <msub> <mi>&kappa;</mi> <mi>s</mi> </msub> <mo>,</mo> <mi>&lambda;</mi> <msub> <mrow> <mo>,</mo> <mi>t</mi> </mrow> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mrow> <msub> <mi>T</mi> <mrow> <mi>w</mi> <mn>2</mn> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>&kappa;</mi> <mi>a</mi> </msub> <mo>,</mo> <msub> <mi>&kappa;</mi> <mi>s</mi> </msub> <mo>,</mo> <mi>&lambda;</mi> <msub> <mrow> <mo>,</mo> <mi>t</mi> </mrow> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>]</mo> </mrow> <mn>2</mn> </msup> </mtd> </mtr> <mtr> <mtd> <msup> <mrow> <mo>[</mo> <mfrac> <mrow> <msubsup> <mi>Q</mi> <mi>R</mi> <mo>&prime;</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>&kappa;</mi> <mi>a</mi> </msub> <mo>,</mo> <msub> <mi>&kappa;</mi> <mi>s</mi> </msub> <mo>,</mo> <mi>&lambda;</mi> <msub> <mrow> <mo>,</mo> <mi>t</mi> </mrow> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>Q</mi> <mi>R</mi> </msub> <mo>(</mo> <msub> <mi>&kappa;</mi> <mi>a</mi> </msub> <mo>,</mo> <msub> <mi>&kappa;</mi> <mi>s</mi> </msub> <mo>,</mo> <mi>&lambda;</mi> <msub> <mrow> <mo>,</mo> <mi>t</mi> </mrow> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mrow> <msub> <mi>Q</mi> <mi>R</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&kappa;</mi> <mi>a</mi> </msub> <mo>,</mo> <msub> <mi>&kappa;</mi> <mi>s</mi> </msub> <mo>,</mo> <mi>&lambda;</mi> <msub> <mrow> <mo>,</mo> <mi>t</mi> </mrow> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>]</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>[</mo> <mfrac> <mrow> <msubsup> <mi>Q</mi> <mi>T</mi> <mo>&prime;</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>&kappa;</mi> <mi>a</mi> </msub> <mo>,</mo> <msub> <mi>&kappa;</mi> <mi>s</mi> </msub> <mo>,</mo> <mi>&lambda;</mi> <msub> <mrow> <mo>,</mo> <mi>t</mi> </mrow> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>Q</mi> <mi>T</mi> </msub> <mo>(</mo> <msub> <mi>&kappa;</mi> <mi>a</mi> </msub> <mo>,</mo> <msub> <mi>&kappa;</mi> <mi>s</mi> </msub> <mo>,</mo> <mi>&lambda;</mi> <msub> <mrow> <mo>,</mo> <mi>t</mi> </mrow> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mrow> <msub> <mi>Q</mi> <mi>T</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&kappa;</mi> <mi>a</mi> </msub> <mo>,</mo> <msub> <mi>&kappa;</mi> <mi>s</mi> </msub> <mo>,</mo> <mi>&lambda;</mi> <msub> <mrow> <mo>,</mo> <mi>t</mi> </mrow> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>]</mo> </mrow> <mn>2</mn> </msup> </mtd> </mtr> </mtable> </mfenced> </mrow>步骤五、判断步骤四中的目标函数值Fobj是否小于设定阈值ξ,若小于,则步骤三中所假设的半透明材料的吸收系数κa、散射系数κs和导热系数λ即为该待测材料的真实吸收系数、散射系数和导热系数;若大于,则返回步骤三,重新根据逆问题算法假设材料的吸收系数κa、散射系数κs和导热系数λ的设定值,并重新计算,直至步骤四中的目标函数值Fobj小于设定阈值ξ,完成利用脉冲激光加热产生的瞬态光热信号的半透明材料吸收系数、散射系数以及导热系数的测量。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号