工件载荷对机床定位误差的影响分析及补偿研究.pdf

侯家林，王宇林，刘璨,等.工件载荷对机床定位误差的影响分析及补偿研究J.机械科学与技术，2023，42（8）：1285-1292工件载荷对机床定位误差的影响分析及补偿研究侯家林，王宇林，刘璨，刘焕牢（广东海洋大学机械与动力工程学院，广东湛江524000）摘要：影响机床精度的因素众多，其中，工件载荷对机床产生的误差同样不可忽视。通过对机床的受载情况进行静力学分析，并利用激光干涉仪对机床的定位误差及其他各项几何误差进行测量，其结果显示，工件载荷对机床的几何误差存在很大影响，尤其在定位精度上，加载 100kg 与空载相差10%左右，并且在进行误差补偿之后，这一差距会提升到 40%左右，在此基础上提出了结合载荷的定位误差补偿方法。该方法解决了螺距误差补偿在实际工况下会产生欠补偿或过度补偿的问题。关键词：工件载荷；力学分析；定位误差；误差补偿中图分类号：TH17文献标志码：ADOI：10.13433/ki.1003-8728.20220051文章编号：1003-8728（2023）08-1285-08Analysis and Compensation Study on Influence of Workpiece Load onPositioning Error of Machine ToolHOUJialin，WANGYulin，LIUCan，LIUHuanlao(SchoolofMechanicalandPowerEngineering,GuangdongOceanUniversity,Zhanjiang524000,Guangdong,China)Abstract：Therearemanyfactorsaffectingtheaccuracyofmachinetools,amongwhichtheerrorcausedbyworkpieceloadcannotbe ignored.Through the macro statics analysis of the loading condition of the machine tool,the positioning error and othergeometricerrorsofthemachinetoolaremeasuredbyusingthelaserinterferometer.Theresultsshowthattheworkpieceloadhasthesignificantinfluenceonthegeometricerrorofthemachinetool,especiallyinthepositioningaccuracy,thedifferencebetweentheloadof100kgandtheno-loadisabout10%,andaftertheerrorcompensation,thedifferencewillbeincreasedtoabout40%.Therefore,thepositioningerrorcompensationmethodcombinedwiththeloadisproposed.Thismethodsolvestheproblemofundercompensationorovercompensationofpitcherrorcompensationinactualworkingconditions,andprovidestheoreticalguidanceforfurtherimprovingthepositioningaccuracyofmachinetoolinworkingconditions.Keywords：workpieceload；mechanicalanalysis；positioningerror；errorcompensation影响机床加工精度的因素很多，如几何误差、热误差、切削力误差和刀具磨损误差等1。除此之外，机床的工件和夹具的重量导致机床工作台受载所引起的误差也值得研究，这一因素也是影响定位精度的重要原因之一，因此科学合理地测量和减小载荷产生的误差是减小机床整体几何误差的必要条件。He 等2提出了一种基于 GRA-POT 的数控机床时域载荷外推方法，该方法可以获得数控机床寿命周期内可能发生的载荷，并且编制了数控机床在恒速条件下的动态切削力谱。章青等3通过分析机床结构和结合面的承载变形原理，分别用有限元法和接触变形理论计算载荷变形误差参数，并且在几何误差计算模型上，叠加载荷变形参数，使模型能适应多种误差的计算。陈振东4提出了将多体系统误差分析理论运用于机床载荷误差模型分析中，并且收稿日期：2021-06-30基金项目：国家自然科学基金项目（51375100）与广东省自然科学基金项目（2021A15150110591）作者简介：侯家林（1997），硕士研究生，研究方向为数控机床运动误差测量及误差补偿，通信作者：刘焕牢，教授，博士，2023 年8 月机械科学与技术August2023第 42 卷第 8 期MechanicalScienceandTechnologyforAerospaceEngineeringVol.42No.8http:/ 等6提出了一种求大跨度重载横梁荷载曲线及荷载误差的方法，并且基于有限元分析方法得到了载荷曲线，建立了基于载荷曲线和几何关系的载荷误差建模方法，为基于多体动力学的机床精度建模提供了理论依据。Zhong 等7通过对卧式机床装配平动轴的静力学分析，建立了机床基底、导轨和滑块以及工作台之间的误差传递模型，并且研究了机床主要零部件装备误差对机床空间几何误差的影响。除此之外，国内很多学者在力对机床几何误差存在影响方面也做了诸多研究，苏妙静8以机床单轴受力模拟的方法为基础，测算出了切削力对机床定位误差的影响并做了相应补偿。Liang 等9将切削力误差分解为10项分量，包括主轴变形、刀具变形等，建立了几何误差和切削力误差综合模型。这些研究为载荷对定位误差的影响试验的测量、分析提供了许多有益的启示。然而，它们主要集中在载荷误差模型建立和力学模型分析上，针对载荷对定位误差的具体影响以及补偿等研究较少。载荷对定位误差影响的实际情况也同样重要，应该引起足够地重视，以便为实际加工中的机床精度要求提供理论依据。在研究受载所引起的定位误差的同时，还需要选择科学合理地几何误差测量方法。几何误差的测量方法可分为直接测量法和间接测量法10。直接测量法是用相应的仪器直接测量单个几何误差。间接测量法测量几何误差时，应给出相应的辨识方法。激光干涉仪11-12和双球杆13-14作为成熟的商业仪器，广泛应用于机床几何误差的精密测量。针对直接测量机床特殊几何误差的方法效率低或无法实现的问题，人们广泛研究并提出了间接测量和识别方法，如 22 线、15 线、14 线、9 线、12 线、13 线法15及空间体对角线测量法16等，这些方法都是通过采集部分几何误差数据，再通过各项几何误差之间的几何关系，辨识出其他各项几何误差，但仍然无法避免了多轴联动调节和误差模型等带来的偏差。在本次试验中，为了直接运用采集到的几何误差数据，减少模型带来的误差，采用了一种基于直接测量的九线误差测量方法17-18它通过在机床不同的空间位置测量不同的误差类型来直接获得对应的几何误差。选择适合的测量方法能够快速获得可靠的几何误差数据，并且有利于后续的误差补偿试验及研究。本文以 VMC1000p 三轴数控加工中心为例，通过测量数控机床受载前后的 X 轴的各项几何误差数据，并结合定位误差理论值的计算方法，提出了结合工件载荷的定位误差补偿方法。1 机床受载情况及影响因素 1.1几何误差类型分析机床每根轴的几何误差包括 3 个平移误差和3 个角度误差。它们主要是由机床的几何误差引起的，如内装零件的不精确性和装配误差等19。几何误差主要来源于数控机床零、部件的制造及装配精度，是与机床零、部件的形位精度相关的一类误差元素。由于机床零、部件本身及其装配过程中存在形状和位置误差，当机床运动部件移动或转动时，这些形状和位置误差会反映到机床的运动部件上，从而产生几何误差。在此次研究中，当机床工作台受到工件载荷时，其装配关系的细微变化以及机床丝杠运动时产生的误差均会反映到机床的定位误差上，因此，在试验中，以测量 X 轴定位误差为主要研究点，并观测了工件载荷对垂直度误差、直线度误差等的影响情况。1.2机床受载参数辨识模型工作台与滑块的受载三维模型简图如图 1 所示，其中将滚动导轨简化成 4 个滑块和导轨的组合。为了方便设置描绘坐标系，现将工作台置于行程终点，它的体坐标系 O2-X2Y2Z2和体参考坐标系O1-X1Y1Z1均设置在工作台参考平面的右上角，方向与机床坐标系的方向一致，如图 2 所示。工件工作台导轨和滑块基座zxy图1工作台受载模型三维图Fig.1Tableloadmodel3DdiagramZ2Z1O2O1O2Y1Y2l5l1l6l3l2l4GxF1,F4X1,X2F2,F3XO2O1whP3P4P1P2X1,X2图2工作台受力简图Fig.2Tableforcediagram1286机械科学与技术第42卷http:/ X2Y2Z2-O2和体参考坐标系 X1Y1Z1-O1中的位置均为（）。当工作台沿 X 轴方向移动 x 时，即体坐标系 X2Y2Z2-O2相对于体参考坐标系 X1Y1Z1-O1沿 X 轴方向移动 x，故工件在工作台体参考坐标系中的位置应是（）。仅考虑工件重量和工作台自身重量的作用，以体参考坐标系 X1Y1Z1-O1为参考坐标系，以 O1点为原点建立静力学平衡方程：Fz=0My=0Mx=0（1）将图 2 中所示各个变量代入式（1），然后，得到平衡方程：F1+F2+F3+F4=Gx+w(F1+F2)l4+(F3+F4)(l2+l4)=Gxl6+wyw(F1+F4)l3+x)+(F2+F3)(l1+l3+x)=Gx(l5+x)+w(xw+x)（2）wGxl5l6式中：为工件的重量；为工作台的重量；和 分别为工作台重心点坐标的绝对值。式（2）用矩阵表示，令FT=F1F2F3F4TFW=FwMxwMywT则有FW=EWFT（3）式中：Fw=Gx+w;Mxw=Gx(l5+x)+w(xw+x);Myw=Gxl6+wyw;EW=1111l3+xl1+l3+xl1+l3+xl3+xl4l4l2+l4l2+l4。矩阵方程求解得FT=F0+kF?k F（4）式中：F0=Gxl1+l3l5l1+wl1+l3xwl1Gxl1l4+l2l5l2l3l1l6l1l2+wl1l4+l2lwl2l5l1ywl1l2Gxl6l4l2+wywl4l20;F=Gxl1+l3l5l1+wl1+l3xwl11Gxl1l4+l2l5l2l3l1l6l1l2+wl1l4+l2lwl2l5l1ywl1l2+1Gxl6l4l2+wywl4l211。F1F2F3F4WF1F2F3F41234P1P2P3P4由式（4）可知，的解并不唯一，一般情况下的标量值并不相同，而唯一的特殊情况（即将工件放置于工作台重心处）在现实情况中很难做到精准放置，故可不做考虑。现只考虑滚动导轨的受载变形，则可假设在，作用下，工作台 4 个滑块的变形量分别为，且变量数值也并不相同，再根据图 2 可知，工作台表面与滑块的接触点，在工作台体参考坐标系中的坐标如下：P1:(l3+x,l4,h1)P2:(l1+l3+x,l4,h2)P3:(l1+l3+x,l2+l4,h3)P4:(l3+x,l2+l4,h4)P1P2P3由，这 3 点互异，可以确定一个平面方程。计算平面上两矢量分别为：p12=P2P1=(l1,0,12)p13=P3P1=(l1,l2,13)根据这两矢量可求出平面的法向矢量 n 以及平面方程，即n=p12 p13=?ijkl1012l1l213?l2(12)(x1l3x)l1(23)(y1l4)+l1l2(z1+h+1)=0（5）P4:(l3+x,l2+l4,h4)由于不考虑工作台工件承载平面的变形，则可将代入式（5），可得滑块变第8期侯家林，等:工件载荷对机床定位误差的影响分析及补偿研究1287http:/ 3 所示。Z1L1X1O1OxhOxOxOxABx(x)zl图3工作台受载误差模型图Fig.3TableloaderrormodeldiagramOxx+xw,yw,0Oxx+xw,yw,hOxOxzx+xw,yw,h+zOx假设在工作台参考坐标系 X1Y1Z1-O1下工作台上 表 面 任 意 点在 该 坐 标 系 下 的 坐 标 为（），则其对应的工作台下表面的点为的坐标为（）。在工件和工作台的自重作用下，点下移到，距离为，则其在工作台参考坐标 系 X1Y1Z1-O1中 的 坐 标（）。显然，的坐标应满足式（5），将点坐标代入，解得z=xwl3l1(12)+ywl4l2(23)1（7）OxOxOxOxx（x）L1P1P2X1O1Z1P1P2为在载荷作用下，工作台体坐标系的上表面任意点产生方位变化后的新位置，与之间位置的水平变化为所求量；为线在平面的投影，空间直线的方程为x1(l3+x)l1=y1l40=z1+h+112（8）X1O1Z1L1在面上的投影方程，即直线的方程为x1(l3+x)l1=z1+h+112（9）从式（9）中可知此直线的斜率为1=12l1L1O1X1则直线与坐标系 X1Y1Z1-O1的轴的夹角为=arctan12l1（10）l1l2l3l4xwywzX(x)zX(x)由式（7）和式（10）可知，为恒定量，和为已知的载荷位置，则工作台出现的下移和微倾仅与滚动导轨的形变有关，因而工作台下移和微倾产生的下移量和倾斜角会引起 X 轴定位误差产生改变，现定义，向上为正，逆时针为正。其中，倾斜角存在的顺时针和逆时针方向的倾斜会使定位误差出现减小和增大两种情况，如图 4 和图 5 所示。ZX激光头线性干涉镜偏移方向测量反射镜x(x)受载前工作台上表面受载后工作台上表面图4定位误差偏移情况 1Fig.4Positioningerroroffsetcase1ZX激光头线性干涉镜偏移方向测量反射镜x(x)受载前工作台上表面受载后工作台上表面图5定位误差偏移情况 2Fig.5Positioningerroroffsetcase2工作台受载后，若倾斜情况如图 4 所示，测量反射镜会向接近激光头方向偏移，使定位误差减小；若倾斜情况如图 5 所示，测量反射镜会向远离激光头方向偏移，使定位误差增大。z1234z12zZ,1(x)Z,2(x)在实际运用中，还需要考虑工作台的厚度和工作台在体参考坐标系 X1Y1Z1-O1倾斜状态。由于与倾斜角和下移量相关的变量，为滚动导轨方便受力计算所简化出的假设量，不具备实际测量价值，因此，工作台的倾斜角和下移量直接由激光干涉仪和球杆仪获得，其中，角度为 X-Z轴 垂 直 度 误 差 的 变 化 量，由 无 载 荷 下 的 X-Z轴垂直度误差和加载荷下的 X-Z 轴垂直度误差获得；为 X 轴沿 Z 轴方向的直线度误差的变化量，由无载荷下的 X 轴沿 Z 轴方向的直线度误差和加载荷下的 X 轴沿 Z 轴方向的直线度误差获得。工作台误差关系如图 6 所示。x(x)OxOxOxOxABz图6工作台误差关系Fig.6Tableerrorrelationx(x)结合图 6 中几何关系可得到待求量的理论值，表达式为1288机械科学与技术第42卷http:/ OxAOxK=AOx/OxOx式中：系数 K=的相似比，即。AOxh OxOxlK=（l+h）/l由于工作台倾斜引起的垂直度误差与工作台的厚度和滑块的形变产生的垂直度误差相比都极小。因此，可近似于，可近似于，则。2 机床受载试验系统数控机床受载试验系统，通过对机床工作台施加载荷来模拟机床在实际加工过程中的工作状态，然后再借助激光干涉仪及相关的测量方法来获得机床的几何误差。2.1试验设计与仪器选用分别设置一组不施加载荷和一组施加 100kg的载荷，将其放置于工作台上偏向重心的同一位置，试验组中，选用 100kg 的载荷可以很好的模拟模具重量，也能够区分误差的产生是由工件载荷产生，而不是随机误差产生的。试验对象：VMC1000p 三轴立式加工中心，载荷选择：100kg 标准砝码，几何误差测量仪器：激光干涉仪（用于测量轴的移动误差和转角误差）、球杆仪（用于测量轴与轴之间的垂直度误差）。仪器安装采集示意图如图 7 所示。图7数控机床误差采集示意图Fig.7Numericalcontrolmachinetoolerroracquisitiondiagram 2.2测量方法与数据采集本次试验为了直接获得数控机床可靠的定位误差和直线度误差数据，采用了 1 种基于直接测量的九线误差测量方法，其测量方向为图 2 所示数控机床 X 轴运动方向。在安装仪器进行测量之前，需要对数控机床的周围环境和工作状态进行检测，首先要保证机床的周围环境没有振动源，其次要将机床运行 30min 后才能进行测量，以消除振动和环境温度对测量结果的影响，同时保证室温范围相对稳定。机床几何误差数据采集方法遵循 ISO230-2 标准，选择合适的间隔距离进行准静态测量，机床在数据采集点停留 10s，并采集 5 个周期。3 数据处理与理论验证Xt(x)以采集的数控机床的几何数据为基础，研究了工件载荷对机床垂直度误差和直线度误差及定位误差的定位精度的影响。通过对比实际测量值，验证理论计算值是有效的。3.1数据处理及分析Xt(x)借助 X-Z 轴的垂直度误差的变化，佐证了受载误差模型中的工作台在受载的情况下会发生微小的下移和倾斜，其中测量得到的下移量和倾斜角度可以为后续的计算提供条件。12选择已经形成商业化模式的球杆仪对机床垂直度误差数据进行采集。其结果如图 8 和图 9 所示，分别为工作台受载 0时的 X-Z 轴垂直度误差和工作台受载 100kg 的垂直度误差。因此可得到机床受载误差模型中的 和 分别为 102.5m/m 和 56.9m/m。从测量数据可知，倾斜角的值为负值，符合图 4 所示情况。运行 1运行 2拟合 1拟合 225%垂直度16%反向间隙 Z16%伺服不匹配14%反向跃冲 Z9%反向跃冲 X102.5 m/m6.4 m5.5 m3.5 m2.5 m21.0 m5.0 m/div圆度0.38 ms+Z+X图80载荷下的垂直度误差Fig.8Perpendicularityerrorunder0load运行 1运行 2拟合 1拟合 214%垂直度14%反向跃冲 Z14%横向间隙 X13%周期误差 Z11%反向间隙 Z圆度56.9 m/m5.6 m6.3 m4.5 m6.2 m3.6 m4.4 m5.0 m/div20.0 m+Z+X图9100kg 载荷下的垂直度误差Fig.9100kg 载荷下的垂直度误差借助激光干涉仪采集到的受载前后的直线度误差及定位误差如图 10 至图 13 所示，其定位精度受载前后对比如表 1 所示。第8期侯家林，等:工件载荷对机床定位误差的影响分析及补偿研究1289http:/ 200 300 400 500 600 700 800 900 1 00050050位置偏差/mX 轴测量目标位置的机床坐标值/mm200 400 600 800 1 000正向偏差平均值反向偏差平均值第 1 次正向移动第 1 次反向移动第 2 次正向移动第 2 次反向移动第 3 次正向移动第 3 次反向移动第 4 次正向移动第 4 次反向移动第 5 次正向移动第 5 次反向移动图100载荷下 X 轴沿 Z 轴方向的直线度误差Fig.10LinearityerrorofXaxisalongZaxisunder0load0100 200 300 400 500 600 700 800 900 1 000X 轴测量目标位置的机床坐标值/mm50050100150200250300350400位置偏差/m0200400600800 1 0005050150250350正向偏差平均值反向偏差平均值第 1 次正向移动第 1 次反向移动第 2 次正向移动第 2 次反向移动第 3 次正向移动第 3 次反向移动第 4 次正向移动第 4 次反向移动第 5 次正向移动第 5 次反向移动图11100kg 载荷下 X 轴沿 Z 轴方向的直线度误差Fig.11LinearityerrorofXaxisalongZaxisunder100kgload0100 200 300 400 500 600 700 800 900 1 000X 轴测量目标位置的机床坐标值/mm35302520151050位置偏差/m0200 400 600 8001 000304020100正向偏差平均值反向偏差平均值第 1 次正向移动第 1 次反向移动第 2 次正向移动第 2 次反向移动第 3 次正向移动第 3 次反向移动第 4 次正向移动第 4 次反向移动第 5 次正向移动第 5 次反向移动图120载荷下的 X 轴定位误差Fig.12X-axispositioningerrorunder0load0100 200 300 400 500 600 700 800 900 1 000X 轴测量目标位置的机床坐标值/mm3025201510505位置偏差/m0200400600800 1 00025351555正向偏差平均值反向偏差平均值第 1 次正向移动第 1 次反向移动第 2 次正向移动第 2 次反向移动第 3 次正向移动第 3 次反向移动第 4 次正向移动第 4 次反向移动第 5 次正向移动第 5 次反向移动图13100kg 载荷下的 X 轴定位误差Fig.13X-axispositioningerrorunder100kgload表 1 定位精度 A 受载前后对比Tab.1PositioningaccuracyAcomparisonbeforeandafterloading载荷重量/kgX轴定位误差/mX轴沿Z轴方向直线度误差034.26731.90010028.96628.490影响比例/%15.510.7上述数据表明机床工作台在有放置工件载荷的情况下，X-Z 轴的垂直度误差会发生明显的改变，并且 X轴沿 Z 轴方向是的直线度误差也会发生改变，从而产生几何误差之间的连锁效应，使 X 轴的定位误差产生改变，因此垂直度误差和直线度误差的改变是机床产生误差的两个重要参考指标。3.2机床受载误差模型理论验证3.2.1直线度误差处理Z1(x)Z2(x)l将在 0 载荷和 100kg 载荷下采集到的直线度误差和进行最小二乘法处理，即将图 10 和图 11 的直线度误差的单向平均偏差经最小二乘法处理，得到的结果如图 14 和图 15，则可以确定工作台的受载最大偏移量=4m。0100 200 300 400 500 600 700 800 900 1 000X 轴测量目标位置的机床坐标值/mm151050510位置偏差/m正向偏差平均值反向偏差平均值图140 载荷下的直线度误差Fig.14Linearityerrorunder0load0100 200 300 400 500 600 700 800 900 1 000X 轴测量目标位置的机床坐标值/mm151050510位置偏差/m正向偏差平均值反向偏差平均值图15100kg 载荷下的直线度误差Fig.15Linearityerrorunder100kgload3.2.2理论值与测量值Xt(x)Xm(x)Xt(x)Xm(x)为了方便观察定位误差的理论变化量和测量变化量，这里分别以图 16 和图 17 进行展示。其中，定位误差理论变化量利用式（11）结合已有的数据可以得到，定位误差测量变化量直接由测量出的 X 轴定位误差得到。1290机械科学与技术第42卷http:/ 200 300 400 500 600 700 800 900 1 000X 轴测量目标位置的机床坐标值/mm012345678单向平均位置偏差/m正向偏差平均值反向偏差平均值图16定位误差变化量的理论值Fig.16Theoreticalvalueofthechangeinpositioningerror0100 200 300 400 500 600 700 800 900 1 000X 轴测量目标位置的机床坐标值/mm0123456789单向平均位置偏差/m正向偏差平均值反向偏差平均值图17定位误差变化量的测量值Fig.17Measurementofthechangeinpositioningerror通过图 16 和图 17 的对比可知，定位误差的理论变化量比测量变化量更为稳定，并且误差范围相对较小，而测量变化量的波动更大，误差范围也更大。产生这种情况是因为定位误差理论变化量计算时，仅仅只考虑到了工作台受载的下移和微倾使定位误差发生的变化，而实际测量时，还存在工作台变形和环境温度改变等。4 定位误差补偿为了减小工件载荷对机床的精度影响，需要对受载机床进行螺距补偿。其考虑载荷和未考虑载荷的补偿结果如图 18，具体数据如表 2。0100 200 300 400 500 600 700 800 900 1 0001816141210864202位置偏差/m0 kg100 kg0200 400 600 800 1 000201510505正向偏差平均值反向偏差平均值X 轴测量目标位置的机床坐标值/mm第 1 次正向移动第 1 次反向移动第 2 次正向移动第 2 次反向移动第 3 次正向移动第 3 次反向移动第 4 次正向移动第 4 次反向移动第 5 次正向移动第 5 次反向移动图18补偿后的定位误差对比Fig.18Comparisonofpositioningerroraftercompensation表 2 定位精度 A 补偿前后对比Tab.2PositioningaccuracyAcomparisonbeforeandaftercompensation载荷重量/kgX轴定位误差/m补偿前补偿后034.26718.41410028.96610.968影响比例/%15.541.9上述数据表明，工件载荷对机床定位精度从未补偿前的 15.5%扩大到补偿后的 41.9%，即当机床定位精度越高时，工件载荷对其的影响也越明显。同样地，考虑到 100kg 工件载荷的螺距误差补偿与未考虑到工件载荷的螺距误差补偿相比，机床的定位精度在未考虑载荷的基础上提升了 41.9%。该成果弥补了仅仅只做空载补偿所造成的欠量或过量误差补偿，考虑到工件载荷对机床影响后，这一情况得到了解决，使得机床的定位精度得到进一步的提升。5 结论1）通过系统研究数控机床工作台的受载与定位误差之间的关系，分析出数控机床工作台在工件的载荷作用下受力是不均匀，这使得机床工作台会出现向左或向右两个方向上的倾斜，因此定位误差会出现增大或减少两种情况。其中，100kg 的载荷对机床定位精度的影响在 15%左右。2）结合九线几何误差辨识法计算出了定位误差理论变化量与实际测量出的定位误差测量变化量在图形上是吻合的，因此，数控机床工作台受载模型是有效的。3）利用数控机床误差补偿单元可减小工件载荷对定位误差的影响。其中，对 100kg 载荷的工件进行补偿后，定位精度比仅进行空载补偿提升了 41.9%。参考文献 余文利,邓小雷,谢长雄,等.结合POE旋量理论的三轴立式加工中心几何误差建模方法J.机械科学与技术,2022,41（1）:111-119.YUWL,DENGXL,XIECX,etal.Geometricerrormodelingmethodologyofthree-axisverticalmachiningcenter based on POE screw theoryJ.MechanicalScience and Technology for Aerospace Engineering,2022,41（1）:111-119.（inChinese）1HEJL,ZHAOXY,LIGF,etal.TimedomainloadextrapolationmethodforCNCmachinetoolsbasedonGRA-POT modelJ.The International Journal of2第8期侯家林，等:工件载荷对机床定位误差的影响分析及补偿研究1291http:/ Manufacturing Technology,2019,103（9）:3799-3812.章青,张海根,刘又午,等.数控机床误差补偿技术及应用载荷误差补偿技术J.制造技术与机床,1999（2）:8-9+20.ZHANGQ,ZHANGHG,LIUYW,etal.ApplicationoferrorcompensationtechniqueforNCmachinetoolsload error compensation techniqueJ.ManufacturingTechnology&MachineTools,1999（2）:8-9+20.（inChinese）3陈振东.数控机床载荷误差及其补偿的研究D.天津:天津大学,2007.CHEN Z D.The research on error of load in NCmachines and compensation methodD.Tianjin:TianjinUniversity,2007.(inChinese)4祝洪祥.考虑工况载荷影响的数控机床可靠性评估方法D.长春:吉林大学,2020.ZHU H X.Reliability evaluation method of NCmachine tools considering the influence of workingloadD.Changchun:Jilin University,2020(inChinese)5YANG W T,LUO R,LU H Q,et al.Load curveanalysis and load errors modeling of long span andheavyloadcrossbeamJ.AdvancedMaterialsResearch,2013,658:475-480.6ZHONGXM,LIUHQ,MAOXY,etal.Influenceand error transfer in assembly process of geometricerrors of a translational axis on volumetric error inmachinetoolsJ.Measurement,2019,140:450-461.7苏妙静.数控机床误差补偿及单轴受力对几何误差影响的研究D.湛江:广东海洋大学,2020.SU M J.Research on error compensation of CNCmachine tool and influence of uniaxial force ongeometric errorD.Zhanjiang:Guangdong OceanUniversity,2020.(inChinese)8LIANG J C,LI H F,YUAN J X,et al.Acomprehensiveerrorcompensationsystemforcorrecting geometric,thermal,and cutting force-induced errorsJ.The International Journal ofAdvanced Manufacturing Technology,1997,13（10）:708-712.9董泽圆,李杰,刘辛军,等.数控机床两种几何误差建模方法有效性试验研究J.机械工程学报,2019,55（5）:137-147.DONGZY,LIJ,LIUXJ,etal.Experimentalstudyonthe effectiveness of two different geometric errormodeling methods for machine toolsJ.Journal ofMechanical Engineering,2019,55（5）:137-147.（inChinese）10WUY,WANGDL,WANGZ,etal.Thekinematic11invariantsintestingerrormotionofmachinetoollinearaxesM/ZHANGXM,WANGNF,HUANGYJ.ProceedingsofASIANMMS2016&CCMMS2016.MechanismandMachineScience.Singapore:Springer,2017,408:1525-1540.VAHEBIM,AREZOOB.Accuracyimprovementofvolumetric error modeling in CNC machine toolsJ.TheInternationalJournalofAdvancedManufacturingTechnology,2018,95（5-8）:2243-2257.12JALALUDINAH,SHUKORMHA,MARDINA,etal.DevelopmentandevaluationofthemachiningperformanceofaCNCgantrydoublemotionmachinetoolindifferentmodesJ.TheInternationalJournalofAdvancedManufacturingTechnology,2017,93（1-4）:1347-1356.13LASEMI A,XUE D Y,GU P H.Accurateidentificationandcompensationofgeometricerrorsof5-axis CNC machine tools using double ball barJ.Measurement Science and Technology,2016,27（5）:055004.14LI J,XIE F G,LIU X J,et al.Geometric erroridentificationandcompensationoflinearaxesbasedonanovel13-linemethodJ.TheInternationalJournalofAdvancedManufacturingTechnology,2016,87（5-8）:2269-2283.15李金和.基于体对角线误差检测的数控机床几何误差辨识方法研究J.组合机床与自动化加工技术,2019（1）:79-85.LI J H.An identification strategy research of CNCmachinetoolgeometricerrorsbasedonbodydiagonalerrors measurementJ.Modular Machine Tool&AutomaticManufacturingTechnique,2019（1）:79-85.（inChinese）16LIU H L,RASHEED M,YOUNIS H H.A linemeasurementmethodforgeometricerrormeasurementof the vertical machining centerJ.SN AppliedSciences,2019,1（4）:324.17RASHEEDM.基于激光干涉仪和球杆仪的立式加工中心几何误差直线测量方法D.湛江:广东海洋大学,2019.RASHEED M.A line measurement method for thegeometricerroroftheverticalmachiningcenterbylaserinterferometerandball-barD.Zhanjiang:GuangdongOceanUnivers