发明名称 基于空间相移的面内变形的快速测量系统及测量方法
摘要 本发明公开了一种基于空间相移的面内变形的快速测量系统及测量方法,其特征是系统构成包括激光器、分光棱镜组、凸透镜、反射镜、扩束镜、载波光纤微调器件、镜头成像和CCD摄像机采集器件;本发明通过改变光纤照射到CCD阵列的倾斜角度和空间方向的不同,实现低频载波和高频散斑的分离,运用傅立叶变换将其低频项提取出来,对低频项进行傅立叶逆变换,得到表示物体面内变形的干涉波面的相位分布,通过加载前后相位图像的相减获得被测物体的真实面内变形。本发明在散斑干涉材料中只使用单幅图像实现快速相位检测的目的,简化了系统,降低了成本,实现了面内变形的动态高速无损测量。
申请公布号 CN103697829B 申请公布日期 2016.08.17
申请号 CN201310754607.4 申请日期 2013.12.31
申请人 合肥工业大学 发明人 王永红;孙建飞;杨连祥
分类号 G01B11/16(2006.01)I 主分类号 G01B11/16(2006.01)I
代理机构 安徽省合肥新安专利代理有限责任公司 34101 代理人 何梅生
主权项 基于空间相移的面内变形的快速测量系统,其特征是:设置激光器(1),其出射光经过分光棱镜(2)分成第一束光和第二束光;所述第一束光经过凸透镜(8)会聚到载波光纤(9)中,并由所述载波光纤(9)将第一束光引导照射到CCD摄像机(12)的靶面阵列上,构成参考光载波光路;所述第二束光经过分光棱镜组(3)分为两束物光,分别是第一束物光和第二束物光,所述两束物光分别依次经过扩束镜和反射镜以互补的角度叠加照射在被测物体(13)的表面形成激光散斑面内干涉场;所述激光散斑面内干涉场由所CCD摄像机(12)的靶面阵列进行采集,形成物光测量光路;所述参考光载波光路和物光测量光路的激光在所述CCD摄像机(12)的靶面阵列上形成干涉,对于所述CCD摄像机(12)的采集信息进行计算获得激光散斑面内干涉场的相位,比较被测物体的表面在变形前和变形后的激光散斑面内干涉相位,获得被测物体表面的面内变形信息;所述基于空间相移的面内变形的快速测量系统的测量方法是按如下步骤进行:a、构建物光与载波光的激光干涉模型H(ξ,η)定义物光测量光路中两束物光叠加照射到被测物体上时,由CCD摄像机采集的被测物体的反射光为U(ξ,η),其中ξ,η分别表示摄像机的CCD靶面阵列上的横坐标和纵坐标;所述两束物光在被测物体表面叠加照射由式(1)表征:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>U</mi><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>U</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>U</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>a</mi><mn>1</mn></msub><msup><mi>e</mi><mrow><msub><mi>i&phi;</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow></mrow></msup><mo>+</mo><msub><mi>a</mi><mn>2</mn></msub><msup><mi>e</mi><mrow><msub><mi>i&phi;</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001014174510000011.GIF" wi="1572" he="86" /></maths>式(1)中U<sub>1</sub>(ξ,η)表示第一束物光,U<sub>2</sub>(ξ,η)表示第二束物光,a<sub>1</sub>和a<sub>2</sub>分别为第一束物光和第二束物光的背景光强,φ<sub>1</sub>(ξ,η)和φ<sub>2</sub>(ξ,η)分别为第一束物光和第二束物光在CCD摄像机的CCD靶面阵列相应坐标中的相位;由于CCD摄像机所记录的是激光场的强度信息,则反射光U(ξ,η)的光强U<sub>T</sub>(ξ,η)由式(2)表征:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>U</mi><mi>T</mi></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>=</mo><mo>&lsqb;</mo><msub><mi>U</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>U</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>&rsqb;</mo><msup><mrow><mo>&lsqb;</mo><msub><mi>U</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>U</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mo>*</mo></msup></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msubsup><mi>a</mi><mn>1</mn><mn>2</mn></msubsup><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>a</mi><mn>2</mn><mn>2</mn></msubsup><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>+</mo><mn>2</mn><msub><mi>a</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><msub><mi>a</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mi>cos</mi><mo>&lsqb;</mo><msub><mi>&phi;</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&phi;</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001014174510000012.GIF" wi="1622" he="167" /></maths>定义参考光载波光路中由CCD摄像机采集的载波光为P(ξ,η),Q<sub>ξ</sub>,Q<sub>η</sub>表示载波光的激光在CCD摄像机的靶面阵列中横坐标ξ和纵坐标η方向上的分量,θ(ξ,η)为所要求解的激光散斑面内干涉相位,建立物光测量光路中的物光与参考光载波光路中的载波光的干涉模型H(ξ,η)如式(3):<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>b</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mi>exp</mi><mo>&lsqb;</mo><mn>2</mn><mi>&pi;</mi><mi>i</mi><mrow><mo>(</mo><msub><mi>Q</mi><mi>&xi;</mi></msub><mi>&xi;</mi><mo>+</mo><msub><mi>Q</mi><mi>&eta;</mi></msub><mi>&eta;</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>+</mo><msubsup><mi>c</mi><mi>k</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mi>exp</mi><mo>&lsqb;</mo><mo>-</mo><mn>2</mn><mi>&pi;</mi><mi>i</mi><mrow><mo>(</mo><msub><mi>Q</mi><mi>&xi;</mi></msub><mi>&xi;</mi><mo>+</mo><msub><mi>Q</mi><mi>&eta;</mi></msub><mi>&eta;</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001014174510000013.GIF" wi="1806" he="71" /></maths>式(3)中:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>b</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>U</mi><mi>T</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>+</mo><msup><mi>P</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001014174510000021.GIF" wi="606" he="70" /></maths>c<sub>k</sub>(ξ,η)=U<sub>T</sub>(ξ,η)P(ξ,η)exp{i[θ(ξ,η)]}θ(ξ,η)=φ<sub>1</sub>(ξ,η)‑φ<sub>2</sub>(ξ,η)b、采用傅立叶变换算法,对激光干涉模型按式(4)进行频域分离:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>F</mi><mi>T</mi><mo>&lsqb;</mo><mi>H</mi><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>=</mo><msub><mi>B</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>&xi;</mi></msub><mo>,</mo><msub><mi>f</mi><mi>&eta;</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>C</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>&xi;</mi></msub><mo>-</mo><msub><mi>Q</mi><mi>&xi;</mi></msub><mo>,</mo><msub><mi>f</mi><mi>&eta;</mi></msub><mo>-</mo><msub><mi>Q</mi><mi>&eta;</mi></msub><mo>)</mo></mrow><mo>+</mo><msubsup><mi>C</mi><mi>k</mi><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>f</mi><mi>&xi;</mi></msub><mo>+</mo><msub><mi>Q</mi><mi>&xi;</mi></msub><mo>,</mo><msub><mi>f</mi><mi>&eta;</mi></msub><mo>+</mo><msub><mi>Q</mi><mi>&eta;</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001014174510000022.GIF" wi="1661" he="71" /></maths>式(4)为物光与载波光的干涉模型的傅立叶变换,在傅立叶变换的频谱中B<sub>k</sub>(f<sub>ξ</sub>,f<sub>η</sub>)位于频谱的中心坐标(0,0),载波光的傅立叶变换频谱C<sub>k</sub>(f<sub>ξ</sub>‑Q<sub>ξ</sub>,f<sub>η</sub>‑Q<sub>η</sub>)和<img file="FDA0001014174510000023.GIF" wi="411" he="71" />是一对共轭的包含激光干涉相位的数组,其在频谱上分别分布在中心为(Q<sub>ξ</sub>,Q<sub>η</sub>)和(‑Q<sub>ξ</sub>,‑Q<sub>η</sub>)的周围邻域上,微调载波光纤的照射角度,等效于改变Q<sub>ξ</sub>,Q<sub>η</sub>的值,进而实现激光散斑干涉强度图像的分频的处理,提取频谱中的低频项,得到频域分离的干涉图像;c、进行相角计算,得到表示面内变形云图:对频域分离的干涉图像进行傅立叶逆变换和取相角,得到激光散斑面内干涉相位θ(ξ,η):θ(ξ,η)=arctan{Im[c<sub>k</sub>(ξ,η)]/Re[c<sub>k</sub>(ξ,η)]}   (5)式(5)中Im[c<sub>k</sub>(ξ,η)]为式(3)中c<sub>k</sub>(ξ,η)的虚部,Re[c<sub>k</sub>(ξ,η)]为c<sub>k</sub>(ξ,η)的实部,被测物体变形前和变形后的两个不同状态的激光散斑面内干涉的相位差相Δ(ξ,η)由式(6)获得:Δ(ξ,η)=θ<sub>1</sub>(ξ,η)‑θ<sub>2</sub>(ξ,η)   (6)式(6)中,θ<sub>2</sub>(ξ,η)为变形之前的激光散斑面内干涉相位,θ<sub>1</sub>(ξ,η)为变形之后的相位,由式(7)获得被测物体表面的面内变形D(ξ,η):<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>D</mi><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>&lambda;</mi><mi>&Delta;</mi><mrow><mo>(</mo><mi>&xi;</mi><mo>,</mo><mi>&eta;</mi><mo>)</mo></mrow></mrow><mrow><mn>4</mn><mi>&pi;</mi><mi>s</mi><mi>i</mi><mi>n</mi><mi>&alpha;</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001014174510000024.GIF" wi="710" he="119" /></maths>式(7)中,λ为激光器的激光波长,α为激光照射到被测物体表面的入射角。
地址 230009 安徽省合肥市包河区屯溪路193号