[发明专利]一种工件台微动部分机械参数估计方法有效
申请号: | 201310295948.X | 申请日: | 2013-07-15 |
公开(公告)号: | CN103366063A | 公开(公告)日: | 2013-10-23 |
发明(设计)人: | 陈兴林;王一光;李晓杰 | 申请(专利权)人: | 哈尔滨工业大学 |
主分类号: | G06F17/50 | 分类号: | G06F17/50 |
代理公司: | 哈尔滨市松花江专利商标事务所 23109 | 代理人: | 杨立超 |
地址: | 150001 黑龙*** | 国省代码: | 黑龙江;23 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 一种工件台微动部分机械参数估计方法,主要涉及工件台微动部分机械参数估计方法。本发明是要解决光刻机工件台微动部分三个水平音圈电机的安装角度和质心位置难于通过传统精密仪器精确测量的问题。一、建立转换关系;二、建立位置关系;三、将位置关系二次求导;四、得出相对驱动力及转矩与曝光区域中心加速度的关系式;五、得出曝光中心位置和音圈电机驱动力之间的关系式;六、设定实际质心位置和安装角度;七、由平面几何关系得出L1-L10的实际值;八、利用步骤三中的方法计算实际Q;九、利用步骤一中的方法计算实际L;十、建立误差模型;十一、列写位移差方程;十二、通过试验测得[x′t y′t r′zt z′t r′xt r′y t]T;十三、仿真得[xt yt rzt zt rxt ryt]T;十四、解方程组得△x0,△y0,△z0,△θ1,△θ2,△θ3。本发明属于超精密制造领域。 | ||
搜索关键词: | 一种 工件 微动 部分 机械 参数估计 方法 | ||
【主权项】:
1.一种工件台微动部分机械参数估计方法,其特征在于工件台微动部分机械参数估计方法按以下步骤实现:一、根据音圈电机的安装位置和角度建立工件台微动部分质心驱动力及转矩与音圈电机驱动力之间的转换关系;二、建立曝光中心与工件台微动部分质心之间的位置关系;三、将步骤二得到的位置关系二次求导并忽略一次微分项,得到工件台微动部分质心与曝光中心的加速度关系;四、根据牛顿第二定律得出作用在工件台微动部分质心上的驱动力及转矩与曝光区域中心加速度的关系式F cog = M cog · Q · x . . y . . r . . z z . . r . . x r . . y ]]> 其中M cog = m 0 0 0 0 0 0 m 0 0 0 0 0 0 J z 0 0 0 0 0 0 m 0 0 0 0 0 0 J x 0 0 0 0 0 0 J y ]]> m为工件台微动部分质量,Jx,Jy,Jz为工件台微动部分相对于x,y,z轴的转动惯量;五、利用步骤一和步骤四所得结果列写关于质心驱动力及转矩的等式,得到曝光中心位置和音圈电机驱动力之间的关系式;六、计算实际含有测量误差的工件台微动部分质心位置和水平方向音圈电机安装角度为x ‾ 0 = x 0 + Δ x 0 y ‾ 0 = y 0 + Δ y 0 z ‾ 0 = z 0 + Δ z 0 θ ‾ 1 = θ 1 + Δ θ 1 θ ‾ 2 = θ 2 + Δ θ 2 θ ‾ 3 = θ 3 + Δ θ 3 , ]]> 其中Δx0,Δy0,Δz0为工件台微动部分质心位置误差,Δθ1,Δθ2,Δθ3为工件台微动部分水平方向音圈电机安装角度误差,为实际含有测量位置误差的工件台微动部分质心位置,为实际含有测量角度误差的工件台微动部分水平方向音圈电机安装角度;七、在工件台微动部分包含质心位置误差Δx0,Δy0,Δz0,并且包含水平方向音圈电机安装角度误差Δθ1,Δθ2,Δθ3的情况下,由平面几何计算出含有测量误差的L1-L10的实际值;L ‾ 1 = L 1 - Δ x 0 2 + Δ y 0 2 · cos [ arctan Δ y 0 Δ x 0 - ( π 2 - θ 1 - Δ θ 1 ) ] L ‾ 2 = L 2 + Δ x 0 2 + Δ y 0 2 · sin ( θ 2 + Δ θ 2 - arctan Δ y 0 Δ x 0 ) L ‾ 3 = ( L 3 + Δ y 0 ) cos ( θ 3 + Δ θ 3 ) L ‾ 4 = L 4 + Δ z 0 L ‾ 5 = L 5 - Δ y 0 L ‾ 6 = L 6 + Δ y 0 L ‾ 7 = L 7 + Δ z 0 L ‾ 8 = L 8 + Δ x 0 L ‾ 9 = L 9 + Δ x 0 L ‾ 10 = L 10 - Δ x 0 ]]> 其中为含有位置误差的质心到f1-f3方向上的距离;为含有位置误差的质心到f1-f3所在平面的竖直距离;为含有位置误差的质心到f4,f5(f6)方向上的距离;为含有位置误差的质心到f4-f6方向上的距离;八、在工件台微动部分包含质心位置误差Δx0,Δy0,Δz0的情况下,利用步骤三中的方法计算含有测量误差的Q可得Q Δ = 1 0 - ( y + y ‾ 0 ) 0 0 z + z ‾ 0 0 1 x + x ‾ 0 0 - ( z + z ‾ 0 ) 0 0 0 1 0 0 0 0 0 0 1 y + y ‾ 0 - ( x + x ‾ 0 ) 0 0 0 0 1 0 0 0 0 0 0 1 ]]> 其中QΔ为工件台微动部分含有质心位置误差Δx0,Δy0,Δz0条件下的转换矩阵Q;九、在工件台微动部分包含质心位置误差Δx0,Δy0,Δz0,并且包含水平方向音圈电机安装角度误差Δθ1,Δθ2,Δθ3的情况下,利用步骤一中的方法计算对应的含有测量误差的转换关系矩阵L可得L Δ = - cos θ ‾ 1 - cos θ ‾ 2 cos θ ‾ 3 0 0 0 sin θ ‾ 1 - sin θ ‾ 2 sin θ ‾ 3 0 0 0 L ‾ 1 L ‾ 2 L ‾ 3 0 0 0 0 0 0 1 1 1 - L ‾ 4 sin θ ‾ 1 L ‾ 4 sin θ ‾ 2 L ‾ 4 sin θ ‾ 3 - L ‾ 5 L ‾ 6 L ‾ 6 - L ‾ 7 cos θ ‾ 1 - L ‾ 7 cos θ ‾ 2 L ‾ 7 cos θ ‾ 3 - L ‾ 8 - L ‾ 9 L ‾ 10 ]]> 其中LΔ为工件台微动部分包含质心位置误差Δx0,Δy0,Δz0、水平方向音圈电机安装角度误差Δθ1,Δθ2,Δθ3条件下的转换矩阵;十、利用步骤八和步骤九所得结果建立误差模型,列写在工件台微动部分包含质心位置误差Δx0,Δy0,Δz0,并且包含水平方向音圈电机安装角度误差Δθ1,Δθ2,Δθ3的情况下,曝光中心位置和音圈电机驱动力之间的关系式:S ′ = D ′ · f 1 f 2 f 3 f 4 f 5 f 6 ]]> 其中D ′ = 1 s 2 · Q Δ - 1 · M - 1 · L Δ = Q Δ - 1 · 1 m s 2 0 0 0 0 0 0 1 m s 2 0 0 0 0 0 0 1 J z s 2 0 0 0 0 0 0 1 m s 2 0 0 0 0 0 0 1 J x s 2 0 0 0 0 0 0 1 J y s 2 · L Δ ]]> S′为在工件台微动部分包含质心位置误差Δx0,Δy0,Δz0,并且包含水平方向音圈电机安装角度误差Δθ1,Δθ2,Δθ3的情况下,工件台微动部分的六自由度输出;十一、列写位移差方程:S ′ - S = 1 s 2 · ( Q Δ - 1 · M - 1 · L Δ - Q - 1 · M - 1 · L ) · f 1 f 2 f 3 f 4 f 5 f 6 ]]> 十二、对工件台微动部分六个音圈电机施加常值恒定力f1-f6,施力时间t,通过试验测得[x′t y′t r′zt z′t r′xt r′yt]T,其中x′t,y′t,r′zt,z′t,r′xt,r′yt为工件台微动部分六个自由度在时刻t的实际输出;十三、利用步骤五的结果,用计算机仿真在给定常值恒力f1-f6作用下,施力时间t,仿真得[xt yt rzt zt rxt ryt]T,其中xt,yt,rzt,zt,rxt,ryt为工件台在无质心位置误差和音圈电机安装角度误差的情况下微动部分六个自由度的仿真输出;十四、将步骤十二和步骤十三的结果代入步骤十一中的位移差方程,令S′=[x′t y′tr′zt z′t r′xt r′yt]T,S=[xt yt rzt zt rxt ryt]T,通过解方程组得Δx0,Δy0,Δz0,Δθ1,Δθ2,Δθ3。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于哈尔滨工业大学,未经哈尔滨工业大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201310295948.X/,转载请声明来源钻瓜专利网。