发明名称 线性调频双光束激光外差测量金属线膨胀系数的装置及方法
摘要 线性调频双光束激光外差测量金属线膨胀系数的装置及方法,涉及工程领域。本发明是为了解决现有的测量方法的测量误差大、测量精度低的问题。本发明包括线性调频激光器、第一平面反射镜、会聚透镜、薄玻璃板、第二平面反射镜、光电探测器、信号处理系统和电热炉,第二平面与薄玻璃板相互平行并等高,采用经线性调频后的激光以入射角θ<sub>0</sub>斜入射至薄玻璃板上,光电探测器输出光电流经傅里叶变换后获得外差信号频率,获得薄玻璃板和平面反射镜之间的距离d,距离的变化量Δd等于待测金属杆的伸长量Δl,数显温控仪读取温度显示值,通过金属线膨胀系数公式获得金属线膨胀系数α。它可用于超精密测量、检测、加工设备和激光雷达系统中。
申请公布号 CN103954645B 申请公布日期 2016.03.02
申请号 CN201410205956.5 申请日期 2014.05.15
申请人 黑龙江大学 发明人 李彦超;杨九如;冉玲苓;高扬;柳春郁;杨瑞海;杜军;丁群;王春晖;马立峰
分类号 G01N25/16(2006.01)I 主分类号 G01N25/16(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 岳泉清
主权项 线性调频双光束激光外差测量金属线膨胀系数的方法,该方法是基于线性调频双光束激光外差测量金属线膨胀系数的装置实现的,所述线性调频双光束激光外差测量金属线膨胀系数的装置包括线性调频激光器(1)、第一平面反射镜(2)、会聚透镜(5)、薄玻璃板(3)、第二平面反射镜(4)、光电探测器(6)、信号处理系统(7)和电热炉(8),线性调频激光器(1)发出的线偏振光经过第一平面反射镜(2)反射之后斜入射到薄玻璃板(3)上,经薄玻璃板(3)透射后入射至第二平面反射镜(4),该入射光经该薄玻璃板(3)与第二平面反射镜(4)多次反射后获得多束反射光,该多束反射光经薄玻璃板(3)的透射之后汇聚至光电探测器(6)的光敏面上;光电探测器(6)的电信号输出端与信号处理系统(7)的电信号输入端连接,薄玻璃板(3)与第二平面反射镜(4)之间的距离为正实数d,第二平面反射镜(4)的非反射面中心与待测金属杆(9)的一端固定连接,所述电热炉(8)用于给待测金属杆(9)加热,所述方法包括:调节电热炉(8)的位置,使得与待测金属杆(9)固定连接的第二平面反射镜(4)的反射面与薄玻璃板(3)相互平行并等高,并使第二平面反射镜(4)的反射面与薄玻璃板(3)之间的距离d为20mm~30mm,采用电热炉(8)对待测金属杆(9)均匀加热,利用数显温控仪(10)监测电热炉(8)内部的温度,读取并记录温度显示值,获得温度变化量为ΔT,同时,线性调频激光器(1)连续采集光电探测器(6)输出的电信号,并对电信号进行处理,获得薄玻璃板和第二平面反射镜之间距离的变化量Δd,该变化量Δd等于待测金属杆的伸长量Δl,根据金属线膨胀系数公式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>&alpha;</mi><mo>=</mo><mfrac><mrow><mi>&Delta;</mi><mi>l</mi></mrow><mrow><msub><mi>l</mi><mn>0</mn></msub><mi>&Delta;</mi><mi>T</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000801192280000011.GIF" wi="206" he="143" /></maths>        (公式11)获得金属线膨胀系数α,式中,l<sub>0</sub>为待测金属杆(9)的原长;其特征在于,线性调频激光器(1)连续采集光电探测器(6)输出的电信号,并对电信号进行处理,获得薄玻璃板和第二平面反射镜之间距离的变化量Δd的过程为:经线性调频后的激光以入射角θ<sub>0</sub>斜入射至薄玻璃板上时,入射光场为:E(t)=E<sub>0</sub>exp{i(ω<sub>0</sub>t+kt<sup>2</sup>)}            (公式1)其中,<img file="FDA0000801192280000012.GIF" wi="173" he="126" />为调频带宽的变化率,T为调频周期,△F为调频带宽,E<sub>0</sub>为入射光场振幅,t为时间,ω<sub>0</sub>为光场角频率,i表示虚数,到达薄玻璃板前表面的光程为L,则t‑L/c时刻到达薄玻璃板前表面的反射光场为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>E</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&alpha;</mi><mn>1</mn></msub><msub><mi>E</mi><mn>0</mn></msub><mi>exp</mi><mo>{</mo><mi>i</mi><mo>&lsqb;</mo><msub><mi>&omega;</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mfrac><mi>L</mi><mi>c</mi></mfrac><mo>)</mo></mrow><mo>+</mo><mi>k</mi><msup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mfrac><mi>L</mi><mi>c</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo><mo>}</mo></mrow>]]></math><img file="FDA0000801192280000021.GIF" wi="837" he="126" /></maths>       (公式2)其中,α<sub>1</sub>=r,r为薄玻璃板的反射率,c为光速,同一时刻经过薄玻璃板后表面反射后透射出前表面的光场为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>E</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&alpha;</mi><mn>2</mn></msub><msub><mi>E</mi><mn>0</mn></msub><mi>exp</mi><mo>{</mo><mi>i</mi><mo>&lsqb;</mo><msub><mi>&omega;</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mfrac><mrow><mi>L</mi><mo>+</mo><mn>2</mn><mi>n</mi><mi>d</mi><mi>cos</mi><mi>&theta;</mi></mrow><mi>c</mi></mfrac><mo>)</mo></mrow><mo>+</mo><mi>k</mi><msup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mfrac><mrow><mi>L</mi><mo>+</mo><mn>2</mn><mi>n</mi><mi>d</mi><mi>cos</mi><mi>&theta;</mi></mrow><mi>c</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mfrac><mrow><mn>2</mn><msub><mi>&omega;</mi><mn>0</mn></msub><mi>n</mi><mi>d</mi><mi>cos</mi><mi>&theta;</mi></mrow><mi>c</mi></mfrac><mo>&rsqb;</mo><mo>}</mo></mrow>]]></math><img file="FDA0000801192280000022.GIF" wi="1597" he="128" /></maths>  (公式3)其中,α<sub>2</sub>=β<sup>2</sup>r′,β为薄玻璃板的透射率,r′为平面反射镜的反射率,d为薄玻璃板和第二平面反射镜之间的距离,θ为入射光折射角,n为平板玻璃折射率,光电探测器接收到的总光场为:E′(t)=E<sub>1</sub>(t)+E<sub>2</sub>(t)          (公式4)则光电探测器输出的光电流为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>I</mi><mo>=</mo><mfrac><mrow><mi>&eta;</mi><mi>e</mi></mrow><mrow><mi>h</mi><mi>&nu;</mi></mrow></mfrac><mfrac><mn>1</mn><mi>Z</mi></mfrac><munder><mrow><mo>&Integral;</mo><mo>&Integral;</mo></mrow><mi>D</mi></munder><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>&lsqb;</mo><msub><mi>E</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>E</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo><msup><mrow><mo>&lsqb;</mo><msub><mi>E</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>E</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mo>*</mo></msup><mi>d</mi><mi>s</mi></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mfrac><mrow><mi>&eta;</mi><mi>e</mi></mrow><mrow><mn>2</mn><mi>h</mi><mi>&nu;</mi></mrow></mfrac><mfrac><mn>1</mn><mi>Z</mi></mfrac><munder><mrow><mo>&Integral;</mo><mo>&Integral;</mo></mrow><mi>D</mi></munder><mo>&lsqb;</mo><msubsup><mi>E</mi><mn>1</mn><mn>2</mn></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>E</mi><mn>2</mn><mn>2</mn></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><msub><mi>E</mi><mn>1</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><msup><msub><mi>E</mi><mn>2</mn></msub><mo>*</mo></msup><mo>(</mo><mi>t</mi><mo>)</mo><mo>+</mo><msup><msub><mi>E</mi><mn>1</mn></msub><mo>*</mo></msup><mo>(</mo><mi>t</mi><mo>)</mo><msub><mi>E</mi><mn>2</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow><mo>&rsqb;</mo><mi>d</mi><mi>s</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000801192280000023.GIF" wi="1235" he="286" /></maths>      (公式5)其中,e为电子电量,Z为探测器表面介质的本征阻抗,η为量子效率,D为探测器光敏面的面积,h为普朗克常数,v为激光频率,差频信号直流项经过低通滤波器滤除后,交流项为中频电流,整理后得到中频电流为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>I</mi><mrow><mi>I</mi><mi>F</mi></mrow></msub><mo>=</mo><mfrac><mrow><mi>&eta;</mi><mi>e</mi></mrow><mrow><mn>2</mn><mi>h</mi><mi>&nu;</mi></mrow></mfrac><mfrac><mn>1</mn><mi>Z</mi></mfrac><munder><mrow><mo>&Integral;</mo><mo>&Integral;</mo></mrow><mi>D</mi></munder><mo>&lsqb;</mo><mrow><mo>(</mo><msub><mi>E</mi><mn>1</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msup><msub><mi>E</mi><mn>2</mn></msub><mo>*</mo></msup><mo>(</mo><mi>t</mi><mo>)</mo><mo>+</mo><msup><msub><mi>E</mi><mn>1</mn></msub><mo>*</mo></msup><mo>(</mo><mi>t</mi><mo>)</mo><msub><mi>E</mi><mn>2</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo><mo>&rsqb;</mo><mi>d</mi><mi>s</mi></mrow>]]></math><img file="FDA0000801192280000024.GIF" wi="942" he="142" /></maths>         (公式6)将(公式2)和(公式3)带入(公式6),得到:<maths num="0006" id="cmaths0006"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><msub><mi>I</mi><mrow><mi>I</mi><mi>F</mi></mrow></msub><mo>=</mo><mfrac><mrow><mi>&eta;</mi><mi>e</mi></mrow><mrow><mi>h</mi><mi>&nu;</mi></mrow></mfrac><mfrac><mi>&pi;</mi><mi>Z</mi></mfrac><msubsup><mi>E</mi><mn>0</mn><mn>2</mn></msubsup><msub><mi>&alpha;</mi><mn>1</mn></msub><msub><mi>&alpha;</mi><mn>2</mn></msub><mi>cos</mi><mo>{</mo><mo>&lsqb;</mo><msub><mi>&omega;</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mfrac><mi>L</mi><mi>c</mi></mfrac><mo>)</mo></mrow><mo>+</mo><mi>k</mi><msup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mfrac><mi>L</mi><mi>c</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo><mo>-</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>&lsqb;</mo><msub><mi>&omega;</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mfrac><mrow><mi>L</mi><mo>+</mo><mn>2</mn><mi>n</mi><mi>d</mi><mi> </mi><mi>cos</mi><mi>&theta;</mi></mrow><mi>c</mi></mfrac><mo>)</mo></mrow><mo>+</mo><mi>k</mi><msup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mfrac><mrow><mi>L</mi><mo>+</mo><mn>2</mn><mi>n</mi><mi>d</mi><mi> </mi><mi>cos</mi><mi>&theta;</mi></mrow><mi>c</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mfrac><mrow><mn>2</mn><msub><mi>&omega;</mi><mn>0</mn></msub><mi>n</mi><mi>d</mi><mi> </mi><mi>cos</mi><mi>&theta;</mi></mrow><mi>c</mi></mfrac><mo>&rsqb;</mo><mo>}</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mfrac><mrow><mi>&eta;</mi><mi>e</mi></mrow><mrow><mi>h</mi><mi>&nu;</mi></mrow></mfrac><mfrac><mi>&pi;</mi><mi>Z</mi></mfrac><msubsup><mi>E</mi><mn>0</mn><mn>2</mn></msubsup><msub><mi>&alpha;</mi><mn>1</mn></msub><msub><mi>&alpha;</mi><mn>2</mn></msub><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mn>4</mn><mi>k</mi><mi>n</mi><mi>d</mi><mi> </mi><mi>cos</mi><mi>&theta;</mi></mrow><mi>c</mi></mfrac><mi>t</mi><mo>-</mo><mfrac><mrow><mn>4</mn><mi>k</mi><mi>n</mi><mi>d</mi><mi> </mi><mi>cos</mi><mi>&theta;</mi><mrow><mo>(</mo><mi>L</mi><mo>+</mo><mi>n</mi><mi>d</mi><mi> </mi><mi>cos</mi><mi>&theta;</mi><mo>)</mo></mrow></mrow><msup><mi>c</mi><mn>2</mn></msup></mfrac><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000801192280000025.GIF" wi="1691" he="390" /></maths>  (公式7)中频电流经过傅里叶变换后获得外差信号的频率为:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>f</mi><mrow><mi>I</mi><mi>F</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn>2</mn><mi>k</mi><mi>n</mi><mi>d</mi><mi> </mi><mi>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi></mrow><mrow><mi>&pi;</mi><mi>c</mi></mrow></mfrac><mo>=</mo><mi>K</mi><mi>d</mi></mrow>]]></math><img file="FDA0000801192280000026.GIF" wi="464" he="135" /></maths>            (公式8)光电探测器(6)输出的光电流表达式经过傅里叶变换之后获得线性调频双光束激光外差信号频谱,从频谱中得到斜入射时线性调频双光束激光外差信号频谱第一个主峰的中心频率和正入射时理论曲线的中心频率数值,根据两个中心频率的比值:ζ=cosθ      (公式10)获得激光经薄玻璃板后折射角θ,根据折射定律获得入射角θ<sub>0</sub>,再根据比例系数公式:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>K</mi><mo>=</mo><mfrac><mrow><mn>2</mn><mi>k</mi><mi>n</mi><mi> </mi><mi>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi></mrow><mrow><mi>&pi;</mi><mi>c</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000801192280000031.GIF" wi="293" he="127" /></maths>     (公式9)获得K的值,当d改变时,根据公式8测出对应d的薄玻璃板和第二平面反射镜之间距离的变化量Δd。
地址 150080 黑龙江省哈尔滨市南岗区学府路74号