删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

基于RTCP功能的五轴联动伺服匹配优化

本站小编 Free考研考试/2021-12-20

陈彦羽1, 关立文1, 常佳豪1, 胡蓝2, 王林泉3
1. 清华大学 机械工程系, 北京 100084;
2. 上海航天设备制造总厂有限公司, 上海 200245;
3. 上海拓璞数控股份有限公司, 上海 201100
收稿日期:2020-11-09
基金项目:国家重大科技专项(2019ZX04025001)
作者简介:陈彦羽(1997-), 男, 博士研究生
通讯作者:关立文, 研究员, E-mail: guanlw@tsinghua.edu.cn

摘要:五轴联动机床的动态性能不仅取决于各单轴伺服进给系统动态特性,各轴之间的动态特性是否匹配也对其有着显著影响。该文针对具有刀具绕中心点旋转(RTCP)功能的AC型双摆头五轴联动铣床,提出了一种通过圆度测试法进行多轴间伺服匹配以提高机床整体动态精度的优化方法。该方法在选取了合理的匹配对象后,可在有效降低用于搭建各单轴机电耦合系统成本的同时保证系统较高的整体动态精度。首先基于该装有AC双摆头的五轴联动铣床的整体结构进行机构运动学分析,基于机床整体结构对其各轴伺服进给系统进行机电耦合建模并在SIMULINK中搭建各轴模型;接着基于RTCP功能设计了一条多轴参与运动的特定轨迹,对各轴在该运动轨迹下的运动关系进行分析;通过合理选择匹配对象,基于圆度测试法对XYZ三轴两两轴间进行伺服匹配;而后基于SIMULINK模型,对优化前后机床在特定轨迹下运动的整体动态性能进行仿真对比。结果表明:在该优化方法下,几乎不需要提高伺服进给系统搭建成本,系统就可达到较高的整体动态性能。
关键词:五轴联动机床机电耦合伺服匹配
Optimization of servo matching for a five-axis machine tool based on the RTCP function
CHEN Yanyu1, GUAN Liwen1, CHANG Jiahao1, HU Lan2, WANG Linquan3
1. Department of Mechanical Engineering, Tsinghua University, Beijing 100084, China;
2. Shanghai Aerospace Equipment Manufacturing General Factory, Shanghai 200245, China;
3. Shanghai Topu CNC Technology Co., Ltd., Shanghai 201100, China

Abstract: The performance of a five-axis machine tool not only depends on the dynamic characteristics of each single-axis servo feed system, but also how well the dynamic characteristics of each axis are matched. An optimization method is developed for A/C axes included precision five-axis machine tools with the rotation tool center point (RTCP) function using multi-axis servo matching based on roundness testing to improve the overall machine accuracy. This method uses matched single-axis servo feed systems to reduce the cost of building each single-axis electromechanical coupling system while ensuring the dynamic accuracy. A kinematics analysis of the machine tool structure is used to build an electromechanical coupling model for each axis servo feed system that is implemented in SIMULINK. Then, a specific trajectory is designed for the RTCP function to analyze the motion relationships for each axis for the trajectory. The method is evaluated by selecting matching single-axis servo feed systems for servo matching between three linear axes of X, Y, and Z based on the roundness test method. The SIMULINK model is then used to compare the dynamic performance of the machine tool before and after optimization for the specific trajectory. The results show that this optimization method improves the machine dynamics with little cost for each single-axis servo feed system.
Key words: five-axis machine toolelectromechanical couplingservo matching
随着近年来与高速转轴相结合的五轴联动机床在高精密的复杂曲线加工中应用越来越广泛,其刀头的动态精度与整体的动态性能也成为国内外科研工作者重视的研究热点之一[1-3]
伺服进给系统性能的高低对于数控机床的动态性能与加工精度具有重要影响[4-5]。对伺服进给系统的机电耦合建模主要有两种方式:1) 在未知系统具体结构的条件下通过辨识来获得整体宏观上的系统模型[6-8];2) 在已知系统理论结构的条件下直接建立系统模型[9-11]
目前针对五轴联动机床伺服进给系统的研究主要集中在各单轴伺服进给系统导致跟踪误差的原因及其对最终加工精度的影响[12-14],而对于各单轴之间伺服匹配关系对最终加工精度的影响的研究却非常匮乏[15]。伺服匹配是指当参与插补曲线加工的伺服轴数量为2个或2个以上时,对参与轴基于伺服控制参数进行联动性能优化与调整。这种基于已有的机械结构与控制设备进行的优化与调整仅需改变内部参数,因此比其他优化方法成本大幅度降低,具有十分重要的工程意义。
在伺服匹配方面,对于两轴间的伺服匹配优化已有较为成熟的伺服插补圆度测试法,但对于多轴参与运动的机械系统如五轴联动机床等尚未有成熟、系统的伺服匹配优化方法。本文提出了一种基于圆度测试法的多轴间伺服匹配优化方法,从机床整体结构、各轴负载及运动关系出发,对各轴的伺服控制参数进行调整优化。本文提出的方法可行性强,且仿真结果验证了其优化的有效性。
1 五轴联动铣床三维模型1.1 机床整体结构描述由AC型双摆头五轴联动铣床实际整体结构简化所得的三维模型如图 1所示。X轴导轨与工作台均与床身固连,X方向两边导轨各有一驱动系统同时对Y轴横梁及其上负载进行驱动,Y轴、Z轴、A轴及C轴上负载均由单一驱动系统进行驱动。
图 1 (网络版彩图)机床实际结构简化三维模型
图选项





该装有AC双摆头的五轴联动铣床的三维原理图如图 2所示。首先,建立工件坐标系。以工作台中心点Op作为坐标系原点,以过该点竖直向上的工作平面法向量为z轴,以过该点的Y轴横梁平行向量为y轴,由右手定则确定x轴。考虑到工作台与床身固连,工件坐标系也作为该系统的惯性坐标系。
图 2 五轴联动铣床三维原理图
图选项





然后,分别建立XYZ 3个方向上的导轨-滑块坐标系。分别取XYZ方向导轨中心轴线一端端点为坐标系原点OxOyOz,3坐标系下的xyz轴正方向均与工件坐标系相同。
最后,建立AC两旋转轴坐标系。取摆头顶端平面中心点为坐标原点Ocxyz轴正方向与工件坐标系相同;取C轴回转轴线与A轴回转轴线交点为A轴坐标系原点Oax轴正方向沿刀主轴轴线方向,y轴正方向沿A轴回转轴线方向,z轴正方向由右手定则确定。
1.2 理想运动变换矩阵设a坐标系下一点坐标转换到b坐标系下的齐次变换矩阵为 ${ }_{a}^{b} \boldsymbol{T} $。根据图 2中系统拓扑结构,建立该机床的正运动学变换矩阵,如式(1)所示。
${ }_{A}^{P} \boldsymbol{T}={ }_{X}^{P} \boldsymbol{T}_{Y}^{X} \boldsymbol{T}_{Z}^{Y} \boldsymbol{T}_{C}^{Z} \boldsymbol{T}_{A}^{C} \boldsymbol{T} .$ (1)
其中各变换矩阵表达式为:
${ }_{Z}^{Y} \boldsymbol{T}=\left[\begin{array}{cccc}1 & 0 & 0 & \Delta x_{z y} \\0 & 1 & 0 & y+\Delta y_{z y} \\0 & 0 & 1 & \Delta z_{z y} \\0 & 0 & 0 & 1\end{array}\right],$
${ }_{Y}^{X} \boldsymbol{T}=\left[\begin{array}{cccc}1 & 0 & 0 & x+\Delta x_{y x} \\0 & 1 & 0 & \Delta y_{y x} \\0 & 0 & 1 & \Delta z_{y x} \\0 & 0 & 0 & 1\end{array}\right],{ }_{X}^{P} \boldsymbol{T}=\left[\begin{array}{cccc}1 & 0 & 0 & \Delta x_{x p} \\0 & 1 & 0 & \Delta y_{x p} \\0 & 0 & 1 & \Delta z_{x p} \\0 & 0 & 0 & 1\end{array}\right],$
${ }_{C}^{Z} \boldsymbol{T}=\left[\begin{array}{cccc}1 & 0 & 0 & \Delta x_{c z} \\0 & 1 & 0 & 0 \\0 & 0 & 1 & z+\Delta z_{c z} \\0 & 0 & 0 & 1\end{array}\right]\left[\begin{array}{cccc}\cos \theta_{c} & -\sin \theta_{c} & 0 & 0 \\\sin \theta_{c} & \cos \theta_{c} & 0 & 0 \\0 & 0 & 1 & 0 \\0 & 0 & 0 & 1\end{array}\right],$
${ }_{A}^{C} \boldsymbol{T}=\left[\begin{array}{cccc}1 & 0 & 0 & 0 \\0 & 1 & 0 & 0 \\0 & 0 & 1 & \Delta z_{a c} \\0 & 0 & 0 & 1\end{array}\right]\left[\begin{array}{cccc}\sin \theta_{a} & 0 & \cos \theta_{a} & 0 \\0 & 1 & 0 & 0 \\-\cos \theta_{a} & 0 & \sin \theta_{a} & 0 \\0 & 0 & 0 & 1\end{array}\right].$
在上述坐标变换矩阵中,xyz分别为XYZ 3个方向上的伺服控制系统的进给量;Δx、Δy、Δz为各坐标系间转换时因整体机械结构导致的XYZ 3个方向上的位移增量常数;θaθc为AC双摆头绕AC两轴分别旋转的角度,即A轴与C轴的伺服控制系统的进给量。通过式(1)即可得到一组工件坐标系下的位姿(xyzθaθc)以及该位姿下五轴铣床各轴进给量的相互转换关系。
2 伺服进给系统建模2.1 平动轴模型本文选取研究的AC型双摆头五轴联动铣床中XYZ方向平动轴使用的是传统滚珠丝杠机械传动结构,具体结构经简化后的平动轴模型如图 3所示。
图 3 平动轴简化模型
图选项





图 3中:J1为主动轴转动惯量;J2为从动轴转动惯量;Tm为电机输出扭矩;θm为电机输出转角;T2为从动轴驱动扭矩;θ2为从动轴转角;Th为从动轴传递至滚珠丝杠的输出扭矩;θh为从动轴驱动滚珠丝杠的输出转角;K2为从动轴的扭转刚度系数;D2为从动轴阻尼系数;Mt为动平台质量;Dt为负载阻尼系数;Xt为动平台输出位移;ft为动平台负载阻力;i为主动轴到从动轴的第1级减速装置减速比;ih为滚珠丝杠减速比,对于单线螺纹有ih=2π/Ph,其中Ph为滚珠丝杠螺距。
由从动轴到动平台间能量守恒可知从动轴转动1周有:
$2 {\rm{ \mathsf{ π} }} T_{\mathrm{h}\_\mathrm{M}}=M_{\mathrm{t}} \ddot{X}_{\mathrm{t}} P_{\mathrm{h}}, 2 {\rm{ \mathsf{ π} }} T_{\mathrm{h}\_\mathrm{D}}=D_{\mathrm{t}} \dot{X}_{\mathrm{t}} P_{\mathrm{h}},$
$2 {\rm{ \mathsf{ π} }} T_{\mathrm{h}\_\mathrm{f}}=f_{\mathrm{t}} P_{\mathrm{h}}, i_{\mathrm{h}}=2 {\rm{ \mathsf{ π} }} / P_{\mathrm{h}}, X_{\mathrm{t}}=\theta_{\mathrm{h}} / i_{\mathrm{h}}.$
其中:Th_MTh_DTh_f分别为从动轴输出转矩中驱动平台移动的转矩、克服阻尼的转矩以及克服负载阻力的转矩。联立上面两式可得
$T_{\mathrm{h}}=T_{\mathrm{h}\_\mathrm{M}}+T_{\mathrm{h}\_\mathrm{D}}+T_{\mathrm{h}\_ \mathrm{f}}=\frac{M_{\mathrm{t}} \ddot{\theta}_{\mathrm{h}}+D_{\mathrm{t}} \dot{\theta}_{\mathrm{h}}+f_{\mathrm{t}} i_{\mathrm{h}}}{i_{\mathrm{h}}^{2}}.$ (2)
由主动轴到从动轴的力矩平衡以及从动轴到动平台的力矩平衡可得:
$T_{\mathrm{m}}=J_{1} \ddot{\theta}_{\mathrm{m}}+\frac{T_{2}}{i_{0}},$ (3)
$\begin{gathered}T_{2}=K_{2}\left(\theta_{2}-\theta_{\mathrm{h}}\right)+D_{2}\left(\dot{\theta}_{2}-\dot{\theta}_{\mathrm{h}}\right)=J_{2} \ddot{\theta}_{\mathrm{h}}+T_{\mathrm{h}}= \\J_{2} \ddot{\theta}_{\mathrm{h}}+\frac{M_{\mathrm{t}} \ddot{\theta}_{\mathrm{h}}+D_{\mathrm{t}} \dot{\theta}_{\mathrm{h}}+f_{\mathrm{t}} i_{\mathrm{h}}}{i_{\mathrm{h}}^{2}}.\end{gathered}$ (4)
联立式(2)—(4),并进行Laplace变换,可得:
$T_{\mathrm{m}}(s)=J_{1} s^{2} \varTheta_{\mathrm{m}}(s)+\frac{\left(K_{2}+D_{2} s\right)\left(\varTheta_{2}(s)-\varTheta_{\mathrm{h}}(s)\right)}{i_{0}},$ (5)
$\begin{gathered}\left(K_{2}+D_{2} s\right)\left(\varTheta_{2}(s)-\varTheta_{\mathrm{h}}(s)\right)=J_{2} s^{2} \varTheta_{\mathrm{h}}(s)+ \\\frac{M_{\mathrm{t}} s^{2} \varTheta_{\mathrm{h}}(s)+D_{\mathrm{t}} s \varTheta_{\mathrm{h}}(s)+F_{\mathrm{t}}(s) i_{\mathrm{h}}}{i_{\mathrm{h}}^{2}}.\end{gathered}$ (6)
联立式(5)和(6),可得平动轴传动系统传递函数为
$\begin{gathered}\varTheta_{\mathrm{h}}(s)=\frac{i_{\mathrm{h}}^{2}\left(K_{2}+D_{2} s\right) \varTheta_{2}(s)}{\left(J_{2} i_{\mathrm{h}}^{2}+M_{\mathrm{t}}\right) s^{2}+\left(D_{\mathrm{t}}+i_{\mathrm{h}}^{2} D_{2}\right) s+i_{\mathrm{h}}^{2} K_{2}}- \\\frac{i_{\mathrm{h}} F_{\mathrm{t}}(s)}{\left(J_{2} i_{\mathrm{h}}^{2}+M_{\mathrm{t}}\right) s^{2}+\left(D_{\mathrm{t}}+i_{\mathrm{h}}^{2} D_{2}\right) s+i_{\mathrm{h}}^{2} K_{2}}.\end{gathered}$ (7)
式(7)可化为系统输入电机转角与动平台位移的传递函数,且由典型二阶系统结构不难验证该系统稳定可控。
2.2 旋转轴模型本文选取研究的AC型双摆头五轴联动铣床中AC方向旋转轴采用的是传统蜗轮蜗杆机械传动结构,经简化后的转动轴模型如图 4所示。
图 4 旋转轴简化模型
图选项





图 4中:Jt为负载转动惯量,Dt为负载阻尼系数,θt为蜗轮转轴输出转角,Tt为蜗轮转轴负载扭矩。
由做功相等可得:
$T_{\mathrm{h}} \theta_{\mathrm{h}}=\left(J_{\mathrm{t}} \ddot{\theta}_{\mathrm{t}}+D_{\mathrm{t}} \dot{\theta}_{\mathrm{t}}+T_{\mathrm{t}}\right) \theta_{\mathrm{t}},$ (8)
$\begin{gathered} T_{\mathrm{h}}=\left(J_{\mathrm{t}} \ddot{\theta}_{\mathrm{t}}+D_{\mathrm{t}} \dot{\theta}_{\mathrm{t}}+T_{\mathrm{t}}\right) \frac{\theta_{\mathrm{t}}}{\theta_{\mathrm{h}}}=\frac{J_{\mathrm{t}} \ddot{\theta}_{\mathrm{t}}+D_{\mathrm{t}} \dot{\theta}_{\mathrm{t}}+T_{\mathrm{t}}}{i_{\mathrm{h}}}= \\\frac{J_{\mathrm{t}} \ddot{\theta}_{\mathrm{h}}+D_{\mathrm{t}} \dot{\theta}_{\mathrm{h}}+T_{\mathrm{t}} i_{\mathrm{h}}}{i_{\mathrm{h}}^{2}} .\end{gathered}$ (9)
由主动轴与蜗杆之间力矩平衡方程可得:
$T_{\mathrm{m}}=J_{1} \ddot{\theta}_{\mathrm{m}}+\frac{T_{2}}{i},$ (10)
$\begin{gathered} T_{2}=K_{2}\left(\theta_{2}-\theta_{\mathrm{h}}\right)+D_{2}\left(\dot{\theta}_{2}-\dot{\theta}_{\mathrm{h}}\right)=J_{2} \ddot{\theta}_{\mathrm{h}}+T_{\mathrm{h}}= \\J_{2} \ddot{\theta}_{\mathrm{h}}+\frac{J_{\mathrm{t}} \ddot{\theta}_{\mathrm{h}}+D_{\mathrm{t}} \dot{\theta}_{\mathrm{h}}+T_{\mathrm{t}} i_{\mathrm{h}}}{i_{\mathrm{h}}^{2}}.\end{gathered}$ (11)
联立式(8)—(11)并进行Laplace变换,与平动轴同理可得系统传递函数为
$\begin{gathered} \varTheta_{\mathrm{h}}(s)=\frac{i_{\mathrm{h}}^{2}\left(K_{2}+D_{2} s\right) \varTheta_{2}(s)}{\left(J_{2} i_{\mathrm{h}}^{2}+J_{\mathrm{t}}\right) s^{2}+\left(D_{\mathrm{t}}+i_{\mathrm{h}}^{2} D_{2}\right) s+i_{\mathrm{h}}^{2} K_{2}}-\\\frac{i_{\mathrm{h}} T_{\mathrm{t}}(s)}{\left(J_{2} i_{\mathrm{h}}^{2}+J_{\mathrm{t}}\right) s^{2}+\left(D_{\mathrm{t}}+i_{\mathrm{h}}^{2} D_{2}\right) s+i_{\mathrm{h}}^{2} K_{2}}.\end{gathered}$ (12)
基于式(12)对传递函数进行分析,根据Routh判据可知该系统亦为稳定可控。
2.3 伺服电机控制系统模型比例-积分-微分(proportion-integration-differ-entiation, PID)控制器在工业控制领域是一种常见、简单并且十分有效的控制器。在许多应用中,可以将PID控制器简化为P、PI或PD控制器进行控制。在通常的伺服电机控制方案中,采用PI控制的电流环和速度环以及P控制的位置环,就可以保证控制系统具有较好的稳定性、较高的精度以及较快的响应。
从电流环开始,以输出角速度作为输出、输入电流作为输入,从内向外解算该控制系统的传递函数。解算各环的推导过程在文[15]中已给出,在此不再赘述,仅给出最后得到的系统整体控制框图,如图 5所示。ΘiΘo分别为输入转角与输出转角。
图 5 伺服电机控制系统框图[15]
图选项





2.4 SIMULINK伺服进给系统机电耦合仿真模型以平动轴为例讨论机电的耦合方式,考虑到两方面问题:1) 若电机系统与传动系统通过直接拼接的方式搭建伺服进给系统模型,将会出现“角加速度 $ \rightarrow J_{1} \rightarrow T_{\mathrm{m}}\left(M_{\mathrm{C}}\right) \rightarrow 1 / J_{\mathrm{e}} \rightarrow$角加速度”的代数环,导致仿真无法进行;2) 实际伺服进给系统中,最终接受位置精度检测并反馈到输入的输出为动平台位移,故位置环反馈信号须由电机轴转角改为动平台位移。
基于直接拼接方式的伺服进给系统模型,将J1折合进入Je中并赋J1零值,即Je=Je+J1J1=0,从而在不改变系统整体性质的条件下消除了代数环;对伺服进给系统位置环进行修改,反馈信号由最终的平动轴位移反馈到输入端。修正后的模型如图 6所示。图中:XinputX轴输入位置指令;ActT_XX轴电机实际输出转矩;Actaccelaration_XX轴电机实际输出角加速度;Real_XX轴上动平台实际位移;Ft_XX轴所受干扰力矩。
图 6 伺服进给系统的修正SIMULINK模型
图选项





旋转轴伺服进给系统建模与平动轴建模同理。仿真过程中涉及的参数初始值如表 1所示。
表 1 五轴机床伺服进给系统初始参数
X Y Z A C
位置环比例增益
Kpp/s-1
20.855 27.452 13.036 5.851 1.575
速度环比例增益
Kpv/(A·s·rad-1)
5.517 4.242 11.843 3.257 3.203
电流环比例增益
Kpi/(V·A-1)
13.433 24.767 7.136 36.521 36.521
速度环积分时间常数
Tiv/s
0.01 0.006 0.006 0.006 0.014
电流环积分时间常数
Tii/s
0.002 0.002 0.002 0.002 0.002
电机力矩常数
Kt/(N·m·A-1)
2.68 1.82 1.92 2.22 2.21
电机反电动势常数
Ke/(V·s·rad-1)
10.23 6.93 7.41 8.52 8.46
电枢绕组电感
La/H
0.003 2 0.005 9 0.001 7 0.008 7 0.008 7
电枢绕组电阻
Ra/Ω
0.157 0.350 0.084 0.520 0.530
电机及丝杠等效转动惯量
Je/(kg·m2)
0.017 8 0.004 5 0.020 5 0.004 5 0.004 5
一级减速比
i0
38.5 38.5 3 1 1
二级减速比
ih
50π 50π 62.5π 45.454 45.454
主动轴惯量
J1/(kg·m2)
0.017 8 0.004 5 0.020 5 0.004 5 0.004 5
从动轴转动惯量
J2/(kg·m2)
0.110 0.028 0.041 0.030 0.030
从动轴阻尼系数
D2/ (kg·m2·s-1)
0.2 0.3 0.5 0.4 0.4
负载阻尼系数
Dt/(kg·s-1)
0.4 0.18 0.12 0.4 0.4
负载(动平台)质量
Mt/kg
1 240 1 003 777 207 557


表选项






3 基于RTCP功能的多轴伺服匹配优化3.1 RTCP功能简述将刀具绕中心点旋转(rotation tool center point, RTCP)功能应用于数控机床运动误差检测,其主要的优势在于:1) 可大幅减小旋转轴运动导致的离散化误差与非线性误差,提高刀尖点实际运动轨迹与理想轨迹的重合度;2) 简化数控机床运动时刀尖点空间位置测量工作,刀尖点在运动过程中理想状态应为相对工件坐标系静止,因而实际运动中刀尖点的运动即为出现的运动误差[12]。RTCP功能关闭与开启时摆头与刀尖点运动状态如图 7所示。
图 7 RTCP功能示意图
图选项





在RTCP功能开启后,当输入的工件坐标系下后置指令为(x, y, z, θa, θc)时,由式(1)可得各轴进给运动指令(x, y, z, θa, θc)的关系表达式为
$\left\{\begin{array}{l}x+\Delta x_{0}=x^{\prime}-L_{0} \cdot \sin \theta_{a} \cdot \cos \theta_{c}, \\y+\Delta y_{0}=y^{\prime}-L_{0} \cdot \sin \theta_{a} \cdot \sin \theta_{c}, \\z+\Delta z_{0}=z^{\prime}+L_{0} \cdot \cos \theta_{a}.\end{array}\right.$ (13)
其中:L0为刀尖点到A轴回转轴线的距离;Δx0、Δy0、Δz0分别为因整体几何结构而存在的XYZ 3个方向上的位移增量常数总和。
3.2 基于转动轴的RTCP轨迹规划由式(13)可得,对AC两轴的转动进行规划即可控制整个数控机床的位姿变化。考虑到AC型双摆头五轴联动铣床中A轴的转动范围为±90°,为避免转动轴运动中出现速度或加速度突变,采用正弦变化的方式规划转动轴的运动。同时,作为后续多轴运动伺服匹配的运动轨迹,在沿该轨迹运动时各轴必须同时参与实际运动。综上,对AC两轴转动进行“8”字轨迹规划,其表达式为
$\left\{\begin{array}{l}\theta_{a}=\frac{{\rm{ \mathsf{ π} }}}{6} \cdot \sin (\omega \cdot t) ,\\\theta_{c}=\frac{{\rm{ \mathsf{ π} }}}{6} \cdot \sin (2 \omega \cdot t).\end{array}\right.$ (14)
Oa在空间中一个完整周期内运动形成的理想轨迹如图 8所示。
图 8 点Oa在空间中运动形成的“8”字轨迹
图选项





机床的数控系统对于机床各轴的运动性能有限制,如最大进给速度等。若在轨迹规划中存在运动性能超限,数控系统调整则会导致机床的整体运动过程与理想轨迹出现偏差,此误差将会影响到对数控机床本身运动误差的测量与判断。因此,须在轨迹规划后对轨迹特性进行分析,使得机床按该轨迹运动时各轴运动状态处于数控系统允许的范围内。
对式(13)求导可得
$\left\{\begin{array}{l}v_{x}=-L_{0} \cdot\left(\omega_{a} \cdot \cos \theta_{a} \cdot \cos \theta_{c}-\omega_{c} \cdot \sin \theta_{a} \cdot \sin \theta_{c}\right) ,\\v_{y}=-L_{0} \cdot\left(\omega_{a} \cdot \cos \theta_{a} \cdot \sin \theta_{c}+\omega_{c} \cdot \sin \theta_{a} \cdot \cos \theta_{c}\right) ,\\v_{z}=-L_{0} \cdot \omega_{a} \cdot \sin \theta_{a}.\end{array}\right.$ (15)
分析式(15)可得各平动轴速度与AC两轴角速度约束关系为
$|v| \leqslant L_{0} \cdot \max \left\{\left|\omega_{a \max }\right|,\left|\omega_{c \max }\right|\right\}.$ (16)
对式(15)求导可得
$\left\{\begin{aligned}a_{x}=&-L_{0} \cdot\left(\alpha_{a} \cdot \cos \theta_{a} \cdot \cos \theta_{c}-\alpha_{c} \cdot \sin \theta_{a} \cdot \sin \theta_{c}-\right.\\&\left(\omega_{a}^{2}+\omega_{c}^{2}\right) \cdot \sin \theta_{a} \cdot \cos \theta_{c}-\\&\left.2 \omega_{a} \omega_{c} \cdot \cos \theta_{a} \cdot \sin \theta_{c}\right), \\a_{y}=&-L_{0} \cdot\left(\alpha_{a} \cdot \cos \theta_{a} \cdot \sin \theta_{c}-\alpha_{c} \cdot \sin \theta_{a} \cdot \cos \theta_{c}-\right.\\&\left(\omega_{a}^{2}+\omega_{c}^{2}\right) \cdot \sin \theta_{a} \cdot \sin \theta_{c}+\\&\left.2 \omega_{a} \omega_{c} \cdot \cos \theta_{a} \cdot \cos \theta_{c}\right), \\a_{z}=&-L_{0} \cdot\left(\alpha_{a} \cdot \sin \theta_{a}+\omega_{a}^{2} \cos \theta_{a}\right).\end{aligned}\right.$ (17)
分析式(17)可得各平动轴加速度与AC两轴角速度、角加速度的约束关系为
$\begin{array}{c}|a| \leqslant L_{0} \cdot\left[\max \left\{\left|\alpha_{a \max }\right|,\left|\alpha_{c \max }\right|\right\}+\right. \\\left.\omega_{a \max }^{2}+\omega_{c \max }^{2}\right].\end{array}$ (18)
对于本节中前述的“8”字轨迹,对式(14)求导可得:
$\begin{array}{c}\left|\omega_{c \max }\right|=2\left|\omega_{a \max }\right|=\frac{{\rm{ \mathsf{ π} }}|\omega|}{3}, \\\left|\alpha_{c \max }\right|=4\left|\alpha_{a \max }\right|=\frac{2 {\rm{ \mathsf{ π} }}|\omega|^{2}}{3}.\end{array}$ (19)
将式(19)代入式(16)和(18)中即可得到在该轨迹下各单轴最大速度和最大加速度与转动轴运动状态间的约束关系。通过调整轨迹定义公式(19)中参数ω就可以限制平动轴的速度、加速度在运行过程中的最大值。通过这种约束关系,可以方便地验证所设计的运动轨迹下的运动状态是否在数控机床允许的范围内,从而保证编程轨迹与数控机床实际驱动轨迹一致。
3.3 基于圆度测试的多轴间伺服匹配本文选取的伺服匹配优化实例为一台AC型双摆头五轴联动铣床,其整体结构与优化前伺服进给系统相关参数已分别在节1.1与节2.4中给出。
圆度测试法是比较成熟的在两不同进给轴伺服系统间进行伺服匹配的优化方法。然而,对于同时有三轴及三轴以上伺服系统参与运动的情况,圆度测试法无法直接应用。本节提出一种基于圆度测试的多轴间伺服匹配方法:在多轴中选取一轴作为匹配对象,将其他各轴分别与该轴进行圆度测试伺服匹配,匹配过程中保持该轴的伺服系统参数不变,仅改变其他各轴的伺服系统参数。
在选取匹配对象时考虑误差累积以及电机性能与伺服性能相匹配两方面问题:1) YZAC 4轴系统结构均安装于X轴导轨上,机床实际运行过程中X轴动态性能如高速运动时的微振等将会影响到其余4轴的定位精度;2) X方向负载最大,行程也较大,XYZAC 5个轴若要达到相同的最大速度及运动精度,X轴所需要的成本更大。
综合考虑以上两方面问题,在保证机床运行平稳的条件下,为了以最小的成本消耗实现相对较好的机床性能,选择X轴伺服系统为匹配对象。
文[6]中给出了包含5个单轴跟随误差及几何误差的AC型五轴联动机床刀尖点位姿表达式,即机床各单轴运动过程中速度越快,各单轴跟随误差通过刀尖点位姿体现得越明显,为避免因跟随误差远小于几何误差而被忽略,测试应在尽量高速的运动状态下进行。
考虑到本文中涉及的五轴联动机床平动轴采用的是西门子1FT7某型电机,其额定转速为3 000 r/min;旋转轴采用的是西门子1FT7某型电机,其额定转速为4 500 r/min,且为与后续在特定轨迹下多轴联动时运动速度相匹配,因此选取运动周期T=30 s、半径r=250 mm的圆轨迹分别对XY两轴及XZ两轴进行圆度测试。
首先通过式(20)计算三平动轴在该圆轨迹下的频率响应特性参数[16]
$f_{\text {circle }} / \mathrm{Hz}=\frac{v /\left(\mathrm{mm} \cdot \mathrm{min}^{-1}\right)}{60 \times 2 {\rm{ \mathsf{ π} }} \ \mathrm{r} / \mathrm{mm}}.$ (20)
计算得到的三平动轴频率响应特性参数均小于0.5 Hz,为较低动态响应状态,轴间可进行圆度测试匹配;A轴与C轴的频率响应特性参数大于0.5 Hz,处于较高动态响应状态,无法与XYZ进行匹配[16]。考虑到处于较高动态响应水平的两旋转轴在节3.2中给出的基于RTCP功能的“8”字轨迹下运动时跟随误差较小,因此本文仅对XYZ三平动轴进行轴间伺服匹配。
选取实际运动过程中Oa点到圆心距离与理论圆半径的方差 $ \sigma_{\text {round }}^{2}$作为机床整体轮廓加工精度的评价指标;同时,为维持各单轴跟随性能在一定水平,选取各单轴在圆度测试中的跟随误差最大值εmax作为各单轴动态特性的评价指标。根据以上两项指标,优化目标函数为
$\rho=\sigma_{\text {round }}^{2}+\lambda \cdot \varepsilon_{\max }.$
其中λ为权重系数,根据对机床的实际性能要求进行选定,为常数。本节优化中取1。
由文[16]可知,位置环比例增益Kpp以及速度环比例增益Kpv之间的差异是圆度测试结果的主要影响因素。采用遗传算法分别对XY轴间控制参数KppKpv以及XZ轴间控制参数KppKpv进行匹配。
YKppKpv优化值搜索范围以及ZKppKpv优化值搜索范围分别为:
$y: K_{\mathrm{pp}} \in[20.855,27.452];$
$K_{\mathrm{pv}} \in[4.242,5.517].$
$z: K_{\mathrm{pp}} \in[13.036,20.855] ;$
$K_{\mathrm{pv}} \in[5.517,11.843].$
适应度函数为
$f=\frac{1}{c \cdot \rho} .$
其中c为大于1的常数。本节优化中取100。
种群大小取50,迭代次数为100,交叉概率为0.5,变异因子为0.000 1。
在SIMULINK中进行圆度测试仿真,XY两轴间匹配及XZ两轴间匹配结果分别如图 910所示。
图 9 (网络版彩图)XY两轴间圆度测试前后精度对比
图选项





图 10 (网络版彩图)XZ两轴间圆度测试前后精度对比
图选项





优化后Y轴伺服系统参数Kpp=22.582,Kpv=5.142,其他参数不变;Z轴伺服系统参数Kpp=17.056,Kpv=10.043,其他参数不变;X轴与两旋转轴伺服系统参数不变。
3.4 基于特定轨迹的多轴伺服匹配优化效果验证根据节3.2中的轨迹规划及特性分析,令ω=π/15;根据机床刀具实际尺寸,令L0=331.445 mm。通过式(13)、(14)及式(16)、(18)、(19),得到具体的空间“8”字运动轨迹,且运动过程中速度始终不超过电机额定转速。为使本文优化方法能对提高实际机床加工精度更具指导意义,选取实际运动过程中Oa点到刀尖点距离与L0的方差 $\sigma_{\text {test }}^{2} $作为精度评价指标。在SIMULINK中进行仿真,结果如图 11所示。
图 11 (网络版彩图)五轴联动伺服匹配优化前后精度对比
图选项





图 11可知,经过以X轴伺服系统为匹配对象的XY两轴伺服匹配以及XZ两轴伺服匹配后,机床多轴联动的动态精度明显提高。
4 结论本文针对具有RTCP功能的AC型双摆头五轴联动铣床,提出了一种精度易于测量、优化成本较低且工程实践性强的基于匹配对象的多轴间伺服匹配精度优化方法。首先通过多体拓扑机构学对铣床进行了运动学分析;而后建立了各轴的传动模型以及电机模型,并在SIMULINK中分别搭建了平动轴与旋转轴的机电耦合伺服进给系统模型;基于RTCP功能规划了一条多轴同时参与的运动轨迹并分析了该轨迹下各轴的运动特性,解决了五轴联动铣床刀尖点空间位置精度测量问题;考虑到误差累积以及电机性能与伺服性能相匹配两方面问题,选取了合理的匹配对象,基于圆度测试法以及遗传算法进行两两轴间匹配,实现了多轴间伺服匹配;最后在完成了轴间匹配后,通过对比匹配优化前后的“8”字轨迹运动精度验证了本文优化方法的有效性。

参考文献
[1] 李杰, 谢福贵, 刘辛军, 等. 五轴数控机床空间定位精度改善方法研究现状[J]. 机械工程学报, 2017, 53(7): 113-128.
LI J, XIE F G, LIU X J, et al. Analysis on the research status of volumetric positioning accuracy improvement methods for five-axis NC machine tools[J]. Journal of Mechanical Engineering, 2017, 53(7): 113-128. (in Chinese)
[2] 王伟, 陶文坚, 李晴朝. 五轴数控机床动态精度检验试件特性研究[J]. 机械工程学报, 2017, 53(1): 101-109.
WANG W, TAO W J, LI Q Z. Research on characteristic of test specimen for five-axis CNC machine tools[J]. Journal of Mechanical Engineering, 2017, 53(1): 101-109. (in Chinese)
[3] 吕盾, 李润泽, 刘辉, 等. 数控机床高速高加速进给下的跟随误差控制策略[J]. 西安交通大学学报, 2018, 52(12): 25-31.
LYU D, LI R Z, LIU H, et al. A control strategy of tracking errors for numerical control machine tools at high speed and acceleration feeding[J]. Journal of Xi'an Jiaotong University, 2018, 52(12): 25-31. (in Chinese)
[4] 徐少华, 曹文成. 数控机床伺服系统中伺服电机的应用[J]. 泰州职业技术学院学报, 2009, 9(3): 48-50.
XU S H, CAO W C. Application of the servo motor in the NC machine tool servo feed system[J]. Journal of Taizhou Polytechnical Institute, 2009, 9(3): 48-50. DOI:10.3969/j.issn.1671-0142.2009.03.018 (in Chinese)
[5] 吴南星, 余冬玲. 数控机床伺服进给系统综合模型的探讨[J]. 机床与液压, 2008, 36(8): 66-68.
WU N X, YU D L. Discussion on the synthesis model of servo feed system for NC machine tool[J]. Machine Tool & Hydraulics, 2008, 36(8): 66-68. (in Chinese)
[6] 吕盾, 刘青, 刘沛, 等. 五轴机床加工零件轮廓误差预测方法[J]. 西安交通大学学报, 2020, 54(2): 9-15.
LYU D, LIU Q, LIU P, et al. Contour error prediction method of part for five-axis machining[J]. Journal of Xi'an Jiaotong University, 2020, 54(2): 9-15. (in Chinese)
[7] CHEN Y Y, HUANG P Y, YEN J Y. Frequency-domain identification algorithms for servo systems with friction[J]. IEEE Transactions on Control Systems Technology, 2002, 10(5): 654-665. DOI:10.1109/TCST.2002.801804
[8] PISLARU C, FREEMAN J M, FORD D G. Modal parameter identification for CNC machine tools using wavelet transform[J]. International Journal of Machine Tools and Manufacture, 2003, 43(10): 987-993. DOI:10.1016/S0890-6955(03)00104-4
[9] FENG B, YANG J, REN J R, et al. Control parameters optimization for servo feed system using an improved genetic algorithm[C]//Proceeding of the 11th World Congress on Intelligent Control and Automation. Shenyang, 2014: 4865-4870.
[10] LIN R K, YAO B, CHEN M H, et al. Application analysis of servo-control system based on PMAC with feed-forward control[C]//Fourth International Conference on Machine Vision (ICMV 2011): Computer Vision and Image Analysis; Pattern Recognition and Basic Technologies. Paris, France, 2012, 8350: 835020.
[11] XIE D, ZHU J Q, WANG F. Fuzzy PID control to feed servo system of CNC machine tools[J]. Procedia Engineering, 2012, 29: 2853-2858. DOI:10.1016/j.proeng.2012.01.403
[12] 邓梦. 基于RTCP的五轴联动数控机床动态误差溯源方法研究[D]. 成都: 电子科技大学, 2014.
DENG M. Research on tracing the source of dynamic error for five-axis CNC machine tool based on RTCP[D]. Chengdu: University of Electronic Science and Technology of China, 2014. (in Chinese)
[13] XI Z Y, LU C H, ZHANG J C. Modeling of computer-controlled AC servo feed system[J]. Materials Science Forum, 2004, 471-472: 216-220. DOI:10.4028/www.scientific.net/MSF.471-472.216
[14] NIE X J. Modeling and simulation for feed servo system in CNC machine tool[J]. Advanced Materials Research, 2012, 510: 345-349. DOI:10.4028/www.scientific.net/AMR.510.345
[15] 闻邦椿. 机械设计手册: 第6卷[M]. 5版. 北京: 机械工业出版社, 2010.
WEN B C. Mechanical design handbook: Vol[M]. 5th ed. Beijing: Machinery Industry Press, 2010. (in Chinese)
[16] 张义红, 张立群. 西门子840D SL伺服插补轴圆度测试的参数优化与研究[J]. 中国科技纵横, 2017(19): 75-76, 78.
ZHANG Y H, ZHANG L Q. Parameter optimization and research on roundness test of Siemens 840D SL servo interpolation axis[J]. China Science & Technology Overview, 2017(19): 75-76, 78. (in Chinese)

相关话题/运动 系统 优化 结构 测试

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 采用力矩电机直驱的数控机床进给系统伺服刚度优化
    陈彦羽1,付萌1,关立文1,常佳豪1,刁磊2,胡蓝21.清华大学机械工程系,北京100084;2.上海航天设备制造总厂有限公司,上海200245收稿日期:2020-12-07基金项目:国家重大科技专项(2019ZX04025001);国家自然科学基金项目(51675301)作者简介:陈彦羽(1997 ...
    本站小编 Free考研考试 2021-12-20
  • PLC指令表程序的一种语法分析方法及其在代码静态测试上的应用
    李学良,赵千川,杨文,SyedNaeemHAIDER清华大学自动化系,北京信息科学与技术国家研究中心,智能与网络化系统研究中心,北京100084收稿日期:2020-12-16基金项目:国家重点研发计划项目(2017YFC0704100,2016YFB0901900);国家自然科学基金项目(61425 ...
    本站小编 Free考研考试 2021-12-20
  • 微型涡喷发动机燃烧室优化设计
    李东杰1,周伯豪1,梁骞2,兰旭东11.清华大学航天航空学院,航空发动机研究中心,北京100084;2.北京航天长征飞行器研究所,北京100076收稿日期:2020-07-27作者简介:李东杰(1994-),男,博士研究生通讯作者:兰旭东,副研究员,E-mail:lanxd@mail.tsinghu ...
    本站小编 Free考研考试 2021-12-20
  • 三维激光断面仪调查系统工作原理及其应用
    徐粒1,呙润华2,彭慧婷1,施鹏程21.新疆大学建筑工程学院,乌鲁木齐830047;2.清华大学土木工程系,北京100084收稿日期:2020-07-07基金项目:国家自然科学基金-地区基金项目(51568063)作者简介:徐粒(1994-),男,硕士研究生通讯作者:呙润华,副教授,E-mail:g ...
    本站小编 Free考研考试 2021-12-20
  • 双余度线控转向机系统的电流均衡余度控制
    米峻男,王通,蔡智凯,连小珉清华大学汽车安全与节能国家重点实验室,北京100084收稿日期:2020-11-06作者简介:米峻男(1992-),男,博士研究生通讯作者:连小珉,教授,lianxm@tsinghua.edu.cn摘要:线控转向系统的安全可靠性不足是制约其应用的瓶颈,其中转向机系统的安全 ...
    本站小编 Free考研考试 2021-12-20
  • 面向中文AMR标注体系的兼语语料库构建及兼语结构识别
    侯文惠1,曲维光1,2,魏庭新2,3,李斌2,顾彦慧1,周俊生11.南京师范大学计算机与电子信息学院,南京210023;2.南京师范大学文学院,南京210097;3.南京师范大学国际文化教育学院,南京210097收稿日期:2020-11-30基金项目:国家自然科学基金面上项目(61772278);江 ...
    本站小编 Free考研考试 2021-12-20
  • 钻铆机器人静刚度建模及优化
    关立文1,陈志雄2,刘春3,薛俊41.清华大学机械工程系,北京100084;2.电子科技大学机械与电气工程学院,成都611173;3.成都飞机工业(集团)有限责任公司,成都610092;4.中国航空制造技术研究院,北京100025收稿日期:2020-10-29基金项目:国家重点研发计划(2017YF ...
    本站小编 Free考研考试 2021-12-20
  • 基于粒子群算法的工件三维膨胀变形下转站参数优化
    刘华森1,陈恳1,王国磊21.清华大学机械工程系,北京100084;2.清华大学摩擦学国家重点实验室,北京100084收稿日期:2020-06-14基金项目:国家自然科学基金资助项目(51975308);工信部民机专项(MJZ-2018-G-75)作者简介:刘华森(1988-),男,博士后通讯作者: ...
    本站小编 Free考研考试 2021-12-20
  • 地震作用下岩体结构及岩性对高陡岩质边坡动力响应特征的影响
    宋丹青,黄进,刘晓丽,王恩志清华大学水沙科学与水利水电工程国家重点实验室,北京100084收稿日期:2021-02-17基金项目:国家重点研发计划项目(2018YFC1504902);中国博士后科学基金(2020M680583);国家自然科学基金项目(52079068);博士后创新人才支持计划(BX ...
    本站小编 Free考研考试 2021-12-20
  • 基于海量光纤测温数据的混凝土坝三维温度场分析系统
    周华维1,2,赵春菊1,2,陈文夫3,周宜红1,2,谭尧升3,刘全4,潘志国5,游皓6,梁志鹏1,2,王放1,2,龚攀31.三峡大学水电工程施工与管理湖北省重点实验室,宜昌443002;2.三峡大学水利与环境学院,宜昌443002;3.中国三峡建设管理有限公司,成都610041;4.武汉大学水利水电 ...
    本站小编 Free考研考试 2021-12-20