发明名称 一种采用多波长剪切干涉测量物体表面形貌的方法
摘要 一种采用多波长剪切干涉测量物体表面形貌的方法,涉及被测物体表面形貌测量、光学实验技术领域。该方法利用多波长剪切干涉测量被测物体表面形貌,同时得到多个波长的激光干涉条纹图,利用这些干涉条纹图像得到被测物体表面的斜率、曲率分布和形貌。本发明装置简单,易于实现,较普通单个波长剪切干涉测量被测物体表面形貌精度大大提高,实现全场实时在线高精度测量被测物体表面形貌。
申请公布号 CN103234478B 申请公布日期 2015.08.12
申请号 CN201310123513.7 申请日期 2013.04.10
申请人 清华大学 发明人 冯雪;张长兴;董雪林
分类号 G01B11/24(2006.01)I 主分类号 G01B11/24(2006.01)I
代理机构 北京鸿元知识产权代理有限公司 11327 代理人 邸更岩
主权项 一种采用多波长剪切干涉测量物体表面形貌的方法,其特征在于,所述方法采用多波长剪切干涉物体表面形貌测量系统,该测量系统包括k个波长不同的激光器、合束装置、扩束装置、分光镜、第一透射光栅、第二透射光栅、透镜、光阑、反射光栅、k个CCD相机,其中k为大于等于2的整数;不同波长激光器发出的激光经过合束装置合成一束激光,经过扩束装置扩束后透过分光镜照射到表面经过处理能反光的被测物体表面,从被测物体表面反射后的激光由分光镜反射后依次经过第一透射光栅和第二透射光栅剪切干涉形成干涉条纹,再透过透镜和光阑仅使+1级次干涉条纹到达反射光栅,通过反射光栅将不同波长的干涉条纹分开反射,分别用CCD相机记录各个不同波长的激光干涉条纹,针对每个不同波长的激光干涉条纹,利用傅里叶光学分析得到+1级次干涉条纹的条纹级数、激光波长和被测物体表面斜率的关系:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>f</mi><mrow><mo>,</mo><mi>x</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mrow><msub><mi>n</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></msub><mo>,</mo><msub><mi>y</mi><mrow><msub><mi>n</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msubsup><mi>n</mi><mi>i</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></msubsup><mi>p</mi></mrow><mrow><mn>2</mn><mi>&Delta;</mi></mrow></mfrac><mo>-</mo><mfrac><msub><mi>&lambda;</mi><mi>i</mi></msub><mrow><mn>4</mn><mi>p</mi></mrow></mfrac><mo>,</mo></mtd><mtd><msubsup><mi>n</mi><mi>i</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mn>0</mn><mo>,</mo><mo>&PlusMinus;</mo><mn>1</mn><mo>,</mo><mo>&PlusMinus;</mo><mn>2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>f</mi><mrow><mo>,</mo><mi>y</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mrow><msub><mi>n</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow></msub><mtext>,</mtext><msub><mi>y</mi><mrow><msub><mi>n</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msubsup><mi>n</mi><mi>i</mi><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></msubsup><mi>p</mi></mrow><mrow><mn>2</mn><mi>&Delta;</mi></mrow></mfrac><mo>-</mo><mfrac><msub><mi>&lambda;</mi><mi>i</mi></msub><mrow><mn>4</mn><mi>p</mi></mrow></mfrac><mo>,</mo></mtd><mtd><msubsup><mi>n</mi><mi>i</mi><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mn>0</mn><mo>,</mo><mo>&PlusMinus;</mo><mn>1</mn><mo>,</mo><mo>&PlusMinus;</mo><mn>2</mn><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr></mtable></mfenced><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>k</mi></mrow>]]></math><img file="FDA00003032569100011.GIF" wi="1954" he="476" /></maths>其中:x为水平方向,y为竖直方向,f<sub>,x</sub>、f<sub>,y</sub>分别为表面形貌函数z=f(x,y)对x、y方向的偏导数,λ<sub>i</sub>分别为所用不同波长的激光波长,<img file="FDA00003032569100012.GIF" wi="375" he="98" /><img file="FDA00003032569100013.GIF" wi="364" he="101" />为波长λ<sub>i</sub>激光干涉条纹对应的物体表面点坐标,k为所用不同波长激光器个数,<img file="FDA00003032569100014.GIF" wi="98" he="83" />为第一透射光栅和第二透射光栅沿水平方向时波长λ<sub>i</sub>激光剪切干涉+1级次干涉条纹的条纹级数,<img file="FDA00003032569100015.GIF" wi="92" he="88" />为第一透射光栅和第二透射光栅沿竖直方向时波长λ<sub>i</sub>激光剪切干涉+1级次干涉条纹的条纹级数,Δ为第一透射光栅和第二透射光栅之间距离,p为透射光栅的栅线间距;这样由多个波长激光剪切干涉得到被测物体表面离散条纹对应的一系列离散点的表面斜率值,利用分段样条插值得到被测物体表面全场的斜率分布S<sub>x</sub>(x,y)和S<sub>y</sub>(x,y),在离散点<img file="FDA00003032569100016.GIF" wi="349" he="99" />和<img file="FDA00003032569100017.GIF" wi="357" he="106" />上有:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>S</mi><mi>x</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mrow><msub><mi>n</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></msub><mo>,</mo><msub><mi>y</mi><mrow><msub><mi>n</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>f</mi><mrow><mo>,</mo><mi>x</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mrow><msub><mi>n</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></msub><mo>,</mo><msub><mi>y</mi><mrow><msub><mi>n</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>S</mi><mi>y</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mrow><msub><mi>n</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow></msub><mo>,</mo><msub><mi>y</mi><mrow><msub><mi>n</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow></msub><mo>)</mo></mrow><mtext>=</mtext><msub><mi>f</mi><mrow><mo>,</mo><mi>y</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mrow><msub><mi>n</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow></msub><mtext>,</mtext><msub><mi>y</mi><mrow><msub><mi>n</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA00003032569100021.GIF" wi="997" he="257" /></maths>利用分段样条插值得到的被测物体表面全场斜率分布S<sub>x</sub>(x,y)和S<sub>y</sub>(x,y)计算被测物体表面曲率κ分布:<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>&kappa;</mi><mi>xx</mi></msub><mo>=</mo><mfrac><mrow><msub><mrow><mo>&PartialD;</mo><mi>S</mi></mrow><mi>x</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><mi>x</mi></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mi>&kappa;</mi><mi>yy</mi></msub><mo>=</mo><mfrac><mrow><msub><mrow><mo>&PartialD;</mo><mi>S</mi></mrow><mi>y</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><mi>y</mi></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mi>&kappa;</mi><mi>xy</mi></msub><mo>=</mo><mfrac><mrow><msub><mrow><mo>&PartialD;</mo><mi>S</mi></mrow><mi>x</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><mi>y</mi></mrow></mfrac></mtd></mtr></mtable></mfenced>]]></math><img file="FDA00003032569100022.GIF" wi="355" he="456" /></maths>其中:κ<sub>xx</sub>、κ<sub>yy</sub>分别为x、y方向物体表面曲率,κ<sub>xy</sub>为物体表面扭率;<img file="FDA00003032569100023.GIF" wi="202" he="126" />为S<sub>x</sub>(x,y)对x的偏导数;<img file="FDA00003032569100024.GIF" wi="205" he="146" />为S<sub>y</sub>(x,y)对y的偏导数;<img file="FDA00003032569100025.GIF" wi="194" he="134" />为S<sub>x</sub>(x,y)对y的偏导数;被测物体表面形貌z=f(x,y)由下式表示:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>z</mi><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mo>&Integral;</mo><msub><mi>S</mi><mi>x</mi></msub><mi>dx</mi><mo>+</mo><mo>&Integral;</mo><mrow><mo>(</mo><msub><mi>S</mi><mi>y</mi></msub><mo>-</mo><mfrac><mrow><mi>d&Integral;</mi><msub><mi>S</mi><mi>x</mi></msub><mi>dx</mi></mrow><mi>dy</mi></mfrac><mo>)</mo></mrow><mi>dy</mi></mrow>]]></math><img file="FDA00003032569100026.GIF" wi="916" he="177" /></maths>其中z=f(x,y)为被测物体表面形貌函数,∫S<sub>x</sub>dx为S<sub>x</sub>(x,y)对x的积分原函数,<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mo>&Integral;</mo><mrow><mo>(</mo><msub><mi>S</mi><mi>y</mi></msub><mo>-</mo><mfrac><mrow><mi>d&Integral;</mi><msub><mi>S</mi><mi>x</mi></msub><mi>dx</mi></mrow><mi>dy</mi></mfrac><mo>)</mo></mrow><mi>dy</mi></mrow>]]></math><img file="FDA00003032569100027.GIF" wi="383" he="162" /></maths>为<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mo>(</mo><msub><mi>S</mi><mi>y</mi></msub><mo>-</mo><mfrac><mrow><mi>d&Integral;</mi><msub><mi>S</mi><mi>x</mi></msub><mi>dx</mi></mrow><mi>dy</mi></mfrac><mo>)</mo></mrow>]]></math><img file="FDA00003032569100028.GIF" wi="310" he="166" /></maths>对y的积分原函数。
地址 100084 北京市海淀区100084信箱82分箱清华大学专利办公室