发明名称 一种基于刚度匹配的重载三坐标数控定位器设计方法
摘要 本发明公开了一种基于刚度匹配的重载三坐标数控定位器设计方法。方法的步骤为:1)建立重载三坐标数控定位器的运动特征矩阵;2)采用有限元法计算重载三坐标数控定位器各部件的刚度并建立其误差模型;3)建立重载三坐标数控定位器的空间定位误差模型,并根据误差模型匹配各部件的刚度并优化其结构。本发明的优点在于:1)根据重载三坐标数控定位器的空间定位误差模型可以获取其空间定位误差数据;2)通过对重载三坐标数控定位器的结构优化,可以适度控制三坐标数控定位器的外形尺寸,降低其制造成本;3)能够提高飞机的调姿精度。
申请公布号 CN102736554A 申请公布日期 2012.10.17
申请号 CN201210232844.X 申请日期 2012.07.05
申请人 浙江大学;西安飞机工业(集团)有限责任公司 发明人 柯映林;盖宇春;梁青霄;何胜强;杨国荣
分类号 G05B19/18(2006.01)I;B64F5/00(2006.01)I 主分类号 G05B19/18(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 张法高
主权项 1.一种基于刚度匹配的重载三坐标数控定位器设计方法,其特征在于它的步骤如下:1)重载三坐标数控定位器能够沿X、Y、Z三个方向进给,通过多个重载三坐标数控定位器的协同运动能够实现飞机或部件的位姿调整,重载三坐标定位器划分为如下四个部件:(1)底座部件:包括底座7、Y向传动系统9、Y轴拖板6;(2)立柱部件:包括立柱5、Z向传动系统4;(3)Z轴滑台部件:包括Z轴滑台3和X向传动系统8;(4)X轴滑台部件:包括X轴滑台2和入位夹紧机构1,球托10;2)根据重载三坐标定位器的结构与运动方式建立坐标系系统,坐标系系统包括底座坐标系{base},立柱坐标系{Oy},Z轴滑台坐标系{Oz},X轴滑台坐标系{Ox},装配坐标系{Oa},并建立重载三坐标数控定位器的运动特征矩阵;底座部件到立柱部件之间的运动特征矩阵为:<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>T</mi><mi>y</mi><mi>base</mi></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><msubsup><mi>x</mi><mi>y</mi><mi>base</mi></msubsup></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><msubsup><mi>y</mi><mi>y</mi><mi>base</mi></msubsup><mo>+</mo><msub><mi>p</mi><mi>y</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><msubsup><mi>z</mi><mi>y</mi><mi>base</mi></msubsup></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>立柱部件到Z轴滑台部件之间的运动特征矩阵为:<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>T</mi><mi>z</mi><mi>y</mi></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><msubsup><mi>x</mi><mi>z</mi><mi>y</mi></msubsup></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><msubsup><mi>y</mi><mi>z</mi><mi>y</mi></msubsup></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><msubsup><mi>z</mi><mi>z</mi><mi>y</mi></msubsup><mo>+</mo><msub><mi>p</mi><mi>z</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>Z轴滑台部件到X轴滑台部件之间的运动特征矩阵为:<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>T</mi><mi>x</mi><mi>z</mi></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><msubsup><mi>x</mi><mi>x</mi><mi>z</mi></msubsup><mo>+</mo><msub><mi>p</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><msubsup><mi>y</mi><mi>x</mi><mi>z</mi></msubsup></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><msubsup><mi>z</mi><mi>x</mi><mi>z</mi></msubsup></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths><img file="FDA00001855977100014.GIF" wi="255" he="51" />为重载三坐标定位器在调姿之前坐标系{n}的原点在坐标系{m}中的坐标值,p<sub>i</sub>为重载三坐标数控定位器各轴的进给量;球托球心在X轴滑台坐标系{Ox}中的坐标为<img file="FDA00001855977100015.GIF" wi="523" he="62" />在装配坐标系{Oa}中的坐标为<img file="FDA00001855977100016.GIF" wi="523" he="62" />底座坐标系{base}的原点在装配坐标系{Oa}中的坐标为<img file="FDA00001855977100021.GIF" wi="489" he="56" />则得到球托球心在装配坐标系下的理论位置:<maths num="0004"><![CDATA[<math><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>P</mi><mi>joint</mi><mi>a</mi></msubsup></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>=</mo><msubsup><mi>T</mi><mi>x</mi><mi>z</mi></msubsup><msubsup><mi>T</mi><mi>z</mi><mi>y</mi></msubsup><msubsup><mi>T</mi><mi>y</mi><mi>base</mi></msubsup><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>P</mi><mi>joint</mi><mi>x</mi></msubsup></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>P</mi><mi>base</mi><mi>a</mi></msubsup></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>3)根据重载三坐标定位器的承载情况建立各部件误差模型;(1)底座部件的误差模型:<maths num="0005"><![CDATA[<math><mrow><mi>&Delta;</mi><msubsup><mi>T</mi><mi>y</mi><mi>base</mi></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mi>&Delta;</mi><msub><mi>y</mi><mi>base</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mi>&Delta;</mi><msub><mi>z</mi><mi>base</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中,Δz<sub>base</sub>为底座部件在Z轴方向上产生的变形,Δy<sub>base</sub>为Y向传动系统中丝杆产生的轴向变形,底座沿Z轴方向的变形为:<maths num="0006"><![CDATA[<math><mrow><mi>&Delta;</mi><msub><mi>z</mi><mi>base</mi></msub><mo>=</mo><mfrac><msubsup><mi>F</mi><mi>z</mi><mi>base</mi></msubsup><msubsup><mi>k</mi><mi>z</mi><mi>base</mi></msubsup></mfrac></mrow></math>]]></maths>其中,<img file="FDA00001855977100025.GIF" wi="92" he="57" />为底座部件所受的Z向载荷,<img file="FDA00001855977100026.GIF" wi="81" he="57" />为底座部件的Z向刚度,Y向传动系统中丝杆变形为:<maths num="0007"><![CDATA[<math><mrow><mi>&Delta;</mi><msubsup><mi>y</mi><mi>y</mi><mi>base</mi></msubsup><mo>=</mo><mfrac><mrow><mi>&mu;</mi><mo>&CenterDot;</mo><msubsup><mi>F</mi><mi>z</mi><mi>base</mi></msubsup><mo>+</mo><msubsup><mi>F</mi><mi>y</mi><mi>base</mi></msubsup></mrow><msubsup><mi>k</mi><mi>y</mi><mi>base</mi></msubsup></mfrac></mrow></math>]]></maths>其中,μ为Y向传动系统中导轨、滑块之间的摩擦系数,<img file="FDA00001855977100028.GIF" wi="92" he="63" />为底座部件所受的Y向载荷,丝杆的轴向刚度为:<maths num="0008"><![CDATA[<math><mrow><msubsup><mi>k</mi><mi>y</mi><mi>base</mi></msubsup><mo>=</mo><mfrac><mrow><mi>A</mi><mo>&CenterDot;</mo><msub><mi>E</mi><mi>S</mi></msub><mo>&CenterDot;</mo><mi>L</mi></mrow><mrow><mi>a</mi><mo>&CenterDot;</mo><mi>b</mi></mrow></mfrac></mrow></math>]]></maths>其中,A为丝杆轴断面面积,E<sub>S</sub>为丝杆弹性模量,L为丝杆安装间距,a和b为丝杆螺母距丝杆两端支撑位置距离;(2)立柱部件误差模型:<maths num="0009"><![CDATA[<math><mrow><mi>&Delta;</mi><msubsup><mi>T</mi><mi>z</mi><mi>y</mi></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mi>&Delta;&beta;</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mo>-</mo><mi>&Delta;&alpha;</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>-</mo><mi>&Delta;&beta;</mi></mtd><mtd><mi>&Delta;&alpha;</mi></mtd><mtd><mn>1</mn></mtd><mtd><mi>&Delta;</mi><msub><mi>z</mi><mi>y</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中Δα为Z向传动系统中导轨绕立柱坐标系{Oy}下X轴的转角,Δβ为Z向传动系统中导轨绕立柱坐标系{Oy}下Y轴的转角,Δz<sub>y</sub>为立柱的Z向变形与Z向传动系统中丝杆的轴向变形之和,导轨的偏转角度为:<maths num="0010"><![CDATA[<math><mrow><mi>&Delta;&theta;</mi><mo>=</mo><mfrac><msubsup><mi>F</mi><mi>i</mi><mi>y</mi></msubsup><msubsup><mi>w</mi><mrow><mi>i</mi><mo>.</mo><mi>j</mi></mrow><mi>y</mi></msubsup></mfrac></mrow></math>]]></maths>其中Δθ为偏转角,<img file="FDA00001855977100032.GIF" wi="55" he="51" />为立柱所受的沿i向的载荷,<img file="FDA00001855977100033.GIF" wi="50" he="56" />为立柱的在i向载荷作用下沿j轴偏转的抗弯刚度,Z向传动系统的丝杆刚度计算公式与Y向传动系统相同;(3)Z轴滑台部件误差模型:<maths num="0011"><![CDATA[<math><mrow><mi>&Delta;</mi><msubsup><mi>T</mi><mi>x</mi><mi>z</mi></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mi>&Delta;</mi><msub><mi>x</mi><mi>z</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mi>&Delta;</mi><msub><mi>z</mi><mi>z</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中Δz<sub>z</sub>为Z轴滑台沿Z轴的压缩变形,Δx<sub>z</sub>为Z轴滑台沿X轴的压缩变形与X向传动系统中的丝杆轴向变形之和,X向传动系统的丝杆刚度计算公式与Y向传动系统相同;(4)X轴滑台部件误差模型:<maths num="0012"><![CDATA[<math><mrow><mi>&Delta;</mi><msubsup><mi>T</mi><mi>x</mi><mi>x</mi></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mi>&Delta;</mi><msub><mi>x</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mi>&Delta;</mi><msub><mi>y</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mi>&Delta;</mi><msub><mi>z</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中Δi<sub>x</sub>轴滑台组件沿i轴的压缩变形,4).建立重载三坐标数控定位器的空间定位误差模型,当考虑重载三坐标定位器的变形情况时,重载三坐标定位器的运动特征矩阵变为:底座部件到立柱部件之间的运动特征矩阵<img file="FDA00001855977100036.GIF" wi="427" he="69" />立柱部件到Z轴滑台之间的运动特征矩阵<img file="FDA00001855977100037.GIF" wi="297" he="64" />Z轴滑台到X轴滑台之间的运动特征矩阵<img file="FDA00001855977100038.GIF" wi="378" he="61" />则得到球窝球心在装配坐标系中的实际位置方程:<maths num="0013"><![CDATA[<math><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mmultiscripts><mi>P</mi><mi>joint</mi><mi>a</mi><mprescripts/><none/><mi>R</mi></mmultiscripts></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>=</mo><mmultiscripts><mi>T</mi><mi>x</mi><mi>z</mi><mprescripts/><none/><mi>R</mi></mmultiscripts><mmultiscripts><mi>T</mi><mi>z</mi><mi>y</mi><mprescripts/><none/><mi>R</mi></mmultiscripts><mmultiscripts><mi>T</mi><mi>y</mi><mi>base</mi><mprescripts/><none/><mi>R</mi></mmultiscripts><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>P</mi><mi>joint</mi><mi>x</mi></msubsup></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>P</mi><mi>base</mi><mi>a</mi></msubsup></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>重载三坐标定位器的空间位置误差模型为:<maths num="0014"><![CDATA[<math><mrow><mo>[</mo><mmultiscripts><mi>P</mi><mi>joint</mi><mi>a</mi><mprescripts/><none/><mi>R</mi></mmultiscripts><mo>]</mo><mo>-</mo><mo>[</mo><msubsup><mi>P</mi><mi>joint</mi><mi>a</mi></msubsup><mo>]</mo></mrow></math>]]></maths>重载三坐标定位器各部件为整体铸造件,各部件的制造成本f(k)随部件刚度k增加而增加,<img file="FDA00001855977100043.GIF" wi="113" he="56" />为重载三坐标定位器的空间定位精度要求,重载三坐标数控定位器的设计优化模型为:<maths num="0015"><![CDATA[<math><mrow><mi>min</mi><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>f</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>k</mi><mn>1</mn></msub><mo>,</mo><msub><mi>k</mi><mn>2</mn></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><msub><mi>k</mi><mi>m</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0016"><![CDATA[<math><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>[</mo><mi>&delta;</mi><msubsup><mi>P</mi><mi>joint</mi><mi>a</mi></msubsup><mo>]</mo><mo>&lt;</mo><mo>[</mo><mmultiscripts><mi>P</mi><mi>joint</mi><mi>a</mi><mprescripts/><none/><mi>R</mi></mmultiscripts><mo>]</mo><mo>-</mo><mo>[</mo><msubsup><mi>P</mi><mi>joint</mi><mi>a</mi></msubsup><mo>]</mo><mo>.</mo></mrow></math>]]></maths>
地址 310027 浙江省杭州市西湖区浙大路38号