发明名称 一种介观尺度零件加工误差预测以及控制方法
摘要 本发明公开了一种介观尺度零件加工误差预测及控制方法,该方法包括以下步骤:首先对待预测的介观尺度零件进行特征分析,提取加工特征的加工精度要求,并根据该加工精度要求确定关键零件特征KPC;然后获取所述待预测的介观尺度零件加工的多工位加工制造过程中每一道工序中的刀具信息、机床信息、装夹信息以及切削参数信息;对于第k道工序,计算因铣削力作用引起的刀具变形量<img file="DDA0000576382950000011.GIF" wi="70" he="79" />和工件变形量<img file="DDA0000576382950000012.GIF" wi="112" he="102" />获得则刀具变形误差<img file="DDA0000576382950000013.GIF" wi="75" he="80" />和工件变形误差<img file="DDA0000576382950000014.GIF" wi="106" he="80" />根据待预测的介观尺度零件的多工位加工制造过程,建立动态工序流误差传递状态空间模型,将动态工序流误差传递状态空间模型进行线性化,获得第k道工序的KPC的误差预测结果q<sub>(k)</sub>,并使用该结果对加工过程进行控制。
申请公布号 CN104298171B 申请公布日期 2017.02.15
申请号 CN201410494301.4 申请日期 2014.09.24
申请人 北京理工大学 发明人 焦黎;王西彬;余璐云;谭方浩;高守峰;刘志兵;梁志强;解丽静
分类号 G05B19/18(2006.01)I 主分类号 G05B19/18(2006.01)I
代理机构 北京理工大学专利中心 11120 代理人 高燕燕;仇蕾安
主权项 一种介观尺度零件加工误差预测方法,其特征在于,该方法包括以下步骤:步骤1、针对待预测的介观尺度零件进行特征分析,提取加工特征的加工精度要求,并根据该加工精度要求确定关键零件特征KPC;步骤2、获取所述待预测的介观尺度零件加工的多工位加工制造过程中每一道工序中的刀具信息、机床信息、装夹信息以及切削参数信息;步骤3、对于第k道工序,其包括l个特征面的加工过程;根据装夹信息预测夹具误差<img file="FDA0001117988030000011.GIF" wi="102" he="78" />根据机床信息确定切削时的刀具跳动、颤振状态,结合刀具信息及切削参数信息预测铣削力,并计算因铣削力作用引起的刀具变形量<img file="FDA0001117988030000012.GIF" wi="72" he="63" />和工件变形量<img file="FDA0001117988030000013.GIF" wi="102" he="70" />同时计算得到第k道工序中刀具误差传递矩阵<img file="FDA0001117988030000014.GIF" wi="83" he="79" />以及第k道工序中工件误差传递矩阵<img file="FDA0001117988030000015.GIF" wi="107" he="71" />则刀具变形误差<img file="FDA0001117988030000016.GIF" wi="283" he="70" />工件变形误差<img file="FDA0001117988030000017.GIF" wi="283" he="71" /><img file="FDA0001117988030000018.GIF" wi="826" he="87" />对于i∈[1,l],i为正整数,N<sub>(k),i</sub>为第k道工序中加工第i个特征面时的刀具误差传递矩阵,<img file="FDA0001117988030000019.GIF" wi="654" he="149" /><img file="FDA00011179880300000110.GIF" wi="83" he="70" />和<img file="FDA00011179880300000111.GIF" wi="73" he="71" />分别为刀具坐标系到第i个特征面坐标系的齐次变换矩阵中的旋转变换矩阵和平移变换矩阵;<img file="FDA00011179880300000112.GIF" wi="846" he="86" />M<sub>(k),i</sub>为第k道工序中加工第i个特征面时的工件误差传递矩阵,其中对于i∈[1,l],i为正整数,<img file="FDA00011179880300000113.GIF" wi="1294" he="150" />I为单位矩阵;步骤4、根据待预测的介观尺度零件的多工位加工制造过程,建立动态工序流误差传递状态空间模型:<img file="FDA0001117988030000021.GIF" wi="894" he="162" />其中x<sub>(*)</sub>为第*道工序中产生的工件表面加工特征制造误差,<img file="FDA0001117988030000022.GIF" wi="82" he="79" />为第k道工序中重定位误差传递矩阵,<img file="FDA0001117988030000023.GIF" wi="84" he="79" />为第k道工序中基准误差传递矩阵,<img file="FDA0001117988030000024.GIF" wi="86" he="78" />为第k道工序中夹具误差传递矩阵,<img file="FDA0001117988030000025.GIF" wi="86" he="85" />为第k道工序中计算误差传递矩阵,<img file="FDA0001117988030000026.GIF" wi="82" he="77" />为第k道工序中计算所有新生成特征的误差传递矩阵,<img file="FDA0001117988030000027.GIF" wi="529" he="86" />和<img file="FDA0001117988030000028.GIF" wi="77" he="71" />均在多工位加工制造过程中计算获得;x<sub>(k‑1)</sub>表示第k‑1道工序中产生的工件表面加工特征制造误差,A<sub>(k)</sub>为x<sub>(k‑1)</sub>在第k道工序中的动态工序流传递矩阵,<img file="FDA0001117988030000029.GIF" wi="581" he="79" /><img file="FDA00011179880300000210.GIF" wi="74" he="76" />表示第k道工序中的夹具误差,<img file="FDA00011179880300000211.GIF" wi="74" he="77" />为第k道工序中<img file="FDA00011179880300000212.GIF" wi="70" he="78" />的动态工序流传递矩阵<img file="FDA00011179880300000213.GIF" wi="325" he="71" /><img file="FDA00011179880300000214.GIF" wi="78" he="70" />表示与切削过程中产生的刀具变形误差和工件变形误差,<img file="FDA00011179880300000215.GIF" wi="356" he="71" /><img file="FDA00011179880300000216.GIF" wi="83" he="70" />为第k道工序中<img file="FDA00011179880300000217.GIF" wi="77" he="71" />的动态工序流传递矩阵,<img file="FDA00011179880300000218.GIF" wi="405" he="70" />w<sub>(k)</sub>为第k道工序的制造过程中的噪声,w<sub>(k)</sub>~N(0,W<sub>k</sub>),其中W<sub>k</sub>为w<sub>(k)</sub>的方差;C<sub>(k)</sub>x<sub>(k)</sub>表示在不考虑测量误差时测量得到的第k道工序中关键零件特征KPC的误差值,C<sub>(k)</sub>为系数;μ<sub>(k)</sub>为测量噪声,μ<sub>(k)</sub>~N(0,U<sub>k</sub>),其中U<sub>k</sub>为μ<sub>(k)</sub>的方差;步骤5、将所述动态工序流误差传递状态空间模型进行线性化,获得第k道工序的KPC的误差预测结果q<sub>(k)</sub>:<maths num="0001"><math><![CDATA[<mrow><msub><mi>q</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msub><mo>=</mo><msub><mi>C</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msub><msubsup><mi>&Phi;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mn>1</mn></msubsup><msub><mi>x</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msub><mo>+</mo><msub><mi>C</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msub><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msubsup><mi>&Phi;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msubsup><msubsup><mi>B</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mi>f</mi></msubsup><msubsup><mi>u</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mi>f</mi></msubsup><mo>+</mo><msub><mi>C</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msub><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msubsup><mi>&Phi;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msubsup><msubsup><mi>B</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mi>m</mi></msubsup><msubsup><mi>u</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mi>m</mi></msubsup><mo>+</mo><msub><mi>C</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msub><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msubsup><mi>&Phi;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msubsup><msub><mi>w</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msub><mo>+</mo><msub><mi>&mu;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msub></mrow>]]></math><img file="FDA00011179880300000219.GIF" wi="1518" he="127" /></maths>式中,<img file="FDA00011179880300000220.GIF" wi="62" he="62" />为状态转移矩阵,且<img file="FDA00011179880300000221.GIF" wi="462" he="70" />其中j=1,2,…,k‑1,当j=k时,<img file="FDA00011179880300000222.GIF" wi="182" he="63" />
地址 100081 北京市海淀区中关村南大街5号