CN108717287B - 半闭环控制方式下数控机床进给系统摩擦误差峰值预测方法 - Google Patents
半闭环控制方式下数控机床进给系统摩擦误差峰值预测方法 Download PDFInfo
- Publication number
- CN108717287B CN108717287B CN201810533269.4A CN201810533269A CN108717287B CN 108717287 B CN108717287 B CN 108717287B CN 201810533269 A CN201810533269 A CN 201810533269A CN 108717287 B CN108717287 B CN 108717287B
- Authority
- CN
- China
- Prior art keywords
- semi
- closed loop
- servo motor
- peak value
- feeding system
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
- G05B19/404—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by control arrangements for compensation, e.g. for backlash, overshoot, tool offset, tool wear, temperature, machine construction errors, load, inertia
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/34—Director, elements to supervisory
- G05B2219/34117—Machine workpiece along, parallel to largest dimension
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Engineering & Computer Science (AREA)
- Human Computer Interaction (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Numerical Control (AREA)
Abstract
本发明公开了一种半闭环控制方式下数控机床进给系统摩擦误差峰值预测方法,通过对半闭环控制方式下的数控机床的进给系统参数进行采集,并根据采集的数据,能够实现对摩擦误差峰值的有效预测,从而可有效评估摩擦误差及其对运行精度的影响,为半闭环控制方式下数控机床进给系统摩擦误差的有效分析与抑制奠定了基础,对提高半闭环控制方式下的数控机床运行精度具有重要意义与应用价值。
Description
技术领域
本发明属于数控机床技术领域,具体涉及一种半闭环控制方式下数控机床进给系统摩擦误差峰值预测方法。
背景技术
数控机床以采用的控制方式分为半闭环和全闭环。与全闭环控制方式下的数控机床相比,半闭环控制方式下的数控机床虽然其运行精度低于全闭环控制方式下的数控机床,但其未安装反馈工作台实际位置的高精度光栅尺,控制系统也未包含机械传动机构,因而具有成本低,使用稳定可靠,在中、低端市场,年均采购量与市场保有量占有绝对优势。
摩擦是存在于数控机床伺服进给系统中的一种复杂非线性外界干扰,摩擦误差由其而引起,呈现凸起尖峰形状,严重影响数控机床进给系统运动精度,非常不利于运动精度的控制。摩擦误差峰值大小直接决定了摩擦误差对运动精度的影响程度。由于半闭环控制方式下数控机床未安装反馈工作台实际位置的光栅尺,因而无法实时获得工作台的实际位置,导致难以有效实时监测数控机床换向过程中产生的摩擦误差,给半闭环控制方式下摩擦误差的评估、分析与控制带来了巨大困难,迫切需要一种半闭环控制方式下摩擦误差峰值预测方法。
发明内容
本发明的目的在于克服上述不足,提供一种半闭环控制方式下数控机床进给系统摩擦误差峰值预测方法,
为了达到上述目的,本发明包括以下步骤:
步骤一,将数控机床设置为在半闭环控制方式下工作,数控机床的进给系统依据位置指令Xr进行精密运动;
步骤二,采集运动过程中位置指令Xr、伺服电机编码器输出转角θm及伺服电机力矩指令电压um;
步骤三,依据数控机床产品与部件信息,获取数控机床进给系统的传动系数rg、力矩常数Kt与系统刚度Ke;
步骤四,通过伺服电机编码器输出转角θm和传动系数rg,得到伺服电机等效位置输出Xm,Xm(t)=θm(t)·rg;
步骤五,通过伺服电机力矩指令电压um和力矩常数Kt,得到伺服电机指令力矩输出Tm,Tm(t)=um(t)·Kt;
步骤六,通过采集获取的位置指令Xr及伺服电机等效位置输出Xm,得到半闭环控制方式下跟随误差Err,Err(t)=Xr(t)-Xm(t);
步骤七,通过半闭环控制方式下跟随误差Err和伺服电机指令力矩输出Tm,得到数控机床换向过程中由摩擦导致的跟随误差峰值Err_sp及对应的伺服电机指令力矩峰值Tm_sp;
步骤八,将跟随误差峰值Err_sp与伺服电机指令力矩峰值Tm_sp代入到半闭环摩擦误差峰值预测式中,求得预测的摩擦误差峰值Err_wp。
步骤一中,半闭环控制方式是将伺服电机编码器输出转角信号θm作为进给系统的位置反馈信号。
步骤七中,跟随误差峰值Err_sp即取半闭环控制方式下换向过程中跟随误差Err的绝对最大值。
步骤七中,伺服电机指令力矩峰值Tm_sp取半闭环控制方式下换向过程中伺服电机指令力矩输出Tm的绝对最大值。
步骤八中,得到半闭环摩擦误差峰值预测式的具体推导过程如下:
数控机床运动过程中,进给系统工作台等效位置输出Xl表示为:
Xl(t)=θl(t)·rg
式中:θl为进给系统工作台等效转角;由于进给系统各传动部件并非刚体,因而会产生变形量Δθ,结合进给系统动力学模型,变形量Δθ表示为:
式中:Tt为负载驱动力矩;Jm为伺服电机转子惯量;
推出进给系统工作台等效转角θl表示为:θl=θm-Δθ,从而半闭环控制方式下数控机床进给系统工作台运动误差Err_l表示为:
根据变形量Δθ和半闭环控制方式下数控机床进给系统工作台运动误差Err_l得到:
数控机床进给系统处于换向过程停驻阶段时,由于伺服电机编码器输出转角θm保持不变,因而d2θm/d2t=0,结合上式得到:
该式为推导出的半闭环摩擦误差峰值预测式,从而求得预测的tslip时刻摩擦误差Err_l,即预测的摩擦误差峰值Err_wp,表示为:
与现有技术相比,本发明通过对半闭环控制方式下的数控机床的进给系统参数进行采集,并根据采集的数据,能够实现对摩擦误差峰值的有效预测,从而可有效评估摩擦误差及其对运行精度的影响,为半闭环控制方式下数控机床进给系统摩擦误差的有效分析与抑制奠定了基础,对提高半闭环控制方式下的数控机床运行精度具有重要意义与应用价值。
附图说明
图1为本发明的控制框图;
图2为本发明中位置指令Xr图;
图3为本发明中伺服电机等效位置Xm与工作台等效位置Xl图;
图4为本发明中跟随误差Err与进给系统工作台运动误差Err_l图;
图5为本发明中换向过程伺服电机指令力矩Tm图。
具体实施方式
下面结合附图对本发明作进一步说明。
本发明包括以下步骤:
1)首先设置数控机床进给系统控制方式为半闭环,此时伺服电机编码器输出转角信号(θm)作为反馈位置信号;
2)数控机床进给系统按照位置指令(Xr)进行运动,在运动过程中采集位置指令(Xr)、伺服电机编码器输出转角(θm)及伺服电机力矩指令电压(um),一般取采样周期为0.001s;
3)依据数控机床产品与部件信息,获取数控机床进给系统的传动系数(rg)、力矩常数(Kt)与系统刚度(Ke);
4)由伺服电机编码器输出转角(θm)、传动系数(rg),可得伺服电机等效位置输出(Xm)可表示为:
Xm(t)=θm(t)·rg (1)
5)由运动过程中采集的位置指令(Xr)及由式(1)计算得到的伺服电机等效位置输出(Xm),获得半闭环控制方式下跟随误差(Err),可表示为:
Err(t)=Xr(t)-Xm(t) (2)
6)由采集获取的伺服电机力矩指令电压(um)、力矩常数(Kt),可得伺服电机指令力矩输出(Tm)可表示为:
Tm(t)=um(t)·Kt (3)
7)数控机床运动过程中,进给系统工作台等效位置输出(Xl)可表示为:
Xl(t)=θl(t)·rg (4)
式中:θl为进给系统工作台等效转角;由于进给系统各传动部件并非刚体,因而会产生变形量Δθ,结合进给系统动力学模型,变形量Δθ可表示为:
式中:Tt为负载驱动力矩;Jm为伺服电机转子惯量;
推出进给系统工作台等效转角θl可表示为:θl=θm-Δθ,从而半闭环控制方式下数控机床进给系统工作台运动误差(Err_l)可表示为:
将式(5)代入式(6)可得:
数控机床进给系统处于换向过程停驻阶段时,由于伺服电机编码器输出转角(θm)保持不变,因而d2θm/d2t=0,结合式(7)得到:
8)由获得的半闭环控制方式下跟随误差(Err)、伺服电机指令力矩输出(Tm),可得数控机床在该工况下某一换向过程中tslip时刻由摩擦导致的跟随误差峰值(Err_sp)及对应的伺服电机指令力矩峰值(Tm_sp);
9)将获得的跟随误差峰值(Err_sp)与伺服电机指令力矩峰值(Tm_sp),代入到式(8)中,该式即为推导出的半闭环摩擦误差峰值预测式,从而求得预测的tslip时刻摩擦误差(Err_l)即预测的摩擦误差峰值(Err_wp),可表示为:
本发明所述的半闭环控制方式下进给系统控制框图,如图1所示,该框图中的主要参数有:位置指令Xr、伺服电机等效位置输出Xm、跟随误差Err、位置环控制增益Kp、速度环比例增益Kvp、速度环积分增益Kvi、速度前馈增益Kfv、加速度前馈增益Kfa、力矩常数Kt、伺服电机力矩指令电压um、伺服电机指令力矩Tm输出、负载驱动力矩Tl、伺服电机转子惯量Jm、伺服电机转子角速度ωm、伺服电机编码器输出转角θm、变形量Δθ、系统刚度Ke、工作台等效惯量Jl、工作台等效角速度ωl、进给系统工作台等效转角θl、传动系数rg、摩擦力矩Tf、工作台等效位置输出Xl。
以下给出本发明的一种具体实施例的实验测试结果:
控制对象为采用三菱交流伺服系统的X-Y轴精密伺服工作台,X轴伺服电机的主要参数,如表1所示。
表1X轴进给系统主要参数
X-Y轴精密伺服工作台以半闭环控制方式工作,在F=3000mm/min,R=20mm工况下做圆轨迹运动,其中X轴位置指令Xr如图2所示,可见有A,B两处换向过程,下面以X轴换向过程B为例进行说明。实验过程中所采集的伺服电机等效位置Xm与工作台等效位置Xl,如图3所示。通过采集获得伺服电机等效位置Xm与工作台等效位置Xl,可得换向过程B处的跟随误差Err与进给系统工作台运动误差Err_l,如图4所示。由图4可得到跟随误差峰值Err_sp=0.0109mm,实际摩擦误差峰值Err_l_sp=0.0129mm,由实验过程中采集的伺服电机力矩指令电压um,可得换向过程B处伺服电机指令力矩Tm如图5所示,由图5可得到伺服电机指令力矩峰值Tm_sp=0.6371N·m。结合表1及上述获得的参数值,代入到半闭环摩擦误差峰值预测式(10)中,可得预测的摩擦误差峰值Err_wp=0.0109mm+0.637×2.5465/612.03=0.0136mm,预测误差=0.0136-0.0129=0.0007mm,可见预测误差值比较小,因此本发明所提出的半闭环控制方式下数控机床进给系统摩擦误差峰值预测方法,可实现对摩擦误差峰值的有效预测。
本发明所述的进给系统位置指令Xr,如图2所示,从图可以看出,运动过程中有A,B两处换向过程。
本发明所述的实验过程中所采集伺服电机等效位置输出Xm与工作台等效位置输出Xl,Xm与Xl的曲线基本重合,说明两者之间误差比较小,如图3所示。
本发明所述的实验过程中所采集换向过程B处的跟随误差Err与进给系统工作台运动误差Err_l,如图4所示。
本发明所述的实验过程中所采集的换向过程B处的伺服电机指令力矩Tm,如图5所示。
综上所述,本发明针对半闭环控制方式下的数控机床难以有效直接获取摩擦误差的缺点,为了实现对半闭环控制方式下数控机床摩擦误差峰值预测,提出了一种半闭环控制方式下数控机床进给系统摩擦误差峰值预测方法。由实验结果可知,由于实验与数据处理过程中,存在外界噪声、系统刚度存在偏差等不利影响,导致预测误差的产生。但预测误差较小,预测的摩擦误差峰值比较接近真实摩擦误差峰值,因而本发明所提出的方法可用于半闭环控制方式下数控机床进给系统摩擦误差峰值的有效预测。
Claims (4)
1.半闭环控制方式下数控机床进给系统摩擦误差峰值预测方法,其特征在于,包括以下步骤:
步骤一,将数控机床设置为在半闭环控制方式下工作,数控机床的进给系统依据位置指令Xr进行精密运动;
步骤二,采集运动过程中位置指令Xr、伺服电机编码器输出转角θm及伺服电机力矩指令电压um;
步骤三,依据数控机床产品与部件信息,获取数控机床进给系统的传动系数rg、力矩常数Kt与系统刚度Ke;
步骤四,通过伺服电机编码器输出转角θm和传动系数rg,得到伺服电机等效位置输出Xm,Xm(t)=θm(t)·rg;
步骤五,通过伺服电机力矩指令电压um和力矩常数Kt,得到伺服电机指令力矩输出Tm,Tm(t)=um(t)·Kt;
步骤六,通过采集获取的位置指令Xr及伺服电机等效位置输出Xm,得到半闭环控制方式下跟随误差Err,Err(t)=Xr(t)-Xm(t);
步骤七,通过半闭环控制方式下跟随误差Err和伺服电机指令力矩输出Tm,得到数控机床换向过程中由摩擦导致的跟随误差峰值Err_sp及对应的伺服电机指令力矩峰值Tm_sp;
步骤八,将跟随误差峰值Err_sp与伺服电机指令力矩峰值Tm_sp代入到半闭环摩擦误差峰值预测式中,求得预测的摩擦误差峰值Err_wp;得到半闭环摩擦误差峰值预测式的具体推导过程如下:
数控机床运动过程中,进给系统工作台等效位置输出Xl表示为:
Xl(t)=θl(t)·rg
式中:θl为进给系统工作台等效转角;由于进给系统各传动部件并非刚体,因而会产生变形量Δθ,结合进给系统动力学模型,变形量Δθ表示为:
式中:Tt为负载驱动力矩;Jm为伺服电机转子惯量;
推出进给系统工作台等效转角θl表示为:θl=θm-Δθ,从而半闭环控制方式下数控机床进给系统工作台运动误差Err_l表示为:
根据变形量Δθ和半闭环控制方式下数控机床进给系统工作台运动误差Err_l得到:
数控机床进给系统处于换向过程停驻阶段时,由于伺服电机编码器输出转角θm保持不变,因而d2θm/d2t=0,结合上式得到:
该式为推导出的半闭环摩擦误差峰值预测式,从而求得预测的tslip时刻摩擦误差Err_l,即预测的摩擦误差峰值Err_wp,表示为:
2.根据权利要求1所述的一种半闭环控制方式下数控机床进给系统摩擦误差峰值预测方法,其特征在于,半闭环控制方式是将伺服电机编码器输出转角信号θm作为进给系统的位置反馈信号。
3.根据权利要求1所述的一种半闭环控制方式下数控机床进给系统摩擦误差峰值预测方法,其特征在于,步骤七中,跟随误差峰值Err_sp即取半闭环控制方式下换向过程中跟随误差Err的绝对最大值。
4.根据权利要求1所述的一种半闭环控制方式下数控机床进给系统摩擦误差峰值预测方法,其特征在于,步骤七中,伺服电机指令力矩峰值Tm_sp取半闭环控制方式下换向过程中伺服电机指令力矩输出Tm的绝对最大值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810533269.4A CN108717287B (zh) | 2018-05-29 | 2018-05-29 | 半闭环控制方式下数控机床进给系统摩擦误差峰值预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810533269.4A CN108717287B (zh) | 2018-05-29 | 2018-05-29 | 半闭环控制方式下数控机床进给系统摩擦误差峰值预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108717287A CN108717287A (zh) | 2018-10-30 |
CN108717287B true CN108717287B (zh) | 2023-07-14 |
Family
ID=63911622
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810533269.4A Active CN108717287B (zh) | 2018-05-29 | 2018-05-29 | 半闭环控制方式下数控机床进给系统摩擦误差峰值预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108717287B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109290843B (zh) * | 2018-11-16 | 2021-02-05 | 西安科技大学 | 精密数控机床倾斜进给系统换向误差峰值预测方法 |
CN114281018A (zh) * | 2021-11-08 | 2022-04-05 | 陕西科控技术产业研究院有限公司 | 一种数控机床倾斜进给系统反向跃冲误差峰值离线预测方法 |
CN115576268B (zh) * | 2022-12-07 | 2023-03-10 | 深圳先进技术研究院 | 一种数控加工设备的跟随误差预测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6029098A (en) * | 1996-07-24 | 2000-02-22 | Toshiba Kikai Kabushiki Kaisha | Numerical controller for automatic axial feed |
CN103048949A (zh) * | 2012-12-22 | 2013-04-17 | 西安交通大学 | 一种数控机床伺服系统摩擦误差补偿方法 |
CN106774166A (zh) * | 2017-01-24 | 2017-05-31 | 西安科技大学 | 一种数控机床摩擦误差形态预测及形态特征参数确定方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001166805A (ja) * | 1999-12-13 | 2001-06-22 | Toshiba Mach Co Ltd | ハイブリッド制御方式の工作機械のロストモーション補正値設定方法およびその方法をコンピュータに実行させるプログラムを記録したコンピュータ読み取り可能な記録媒体および数値制御工作機械 |
CN1258431C (zh) * | 2004-03-31 | 2006-06-07 | 清华大学 | 数控机床误差补偿方法及其系统 |
CN104483900B (zh) * | 2014-10-28 | 2018-04-17 | 西安交通大学 | 半闭环控制数控机床滚珠丝杠进给系统定位误差建模方法 |
CN106323633A (zh) * | 2016-10-28 | 2017-01-11 | 华中科技大学 | 一种基于指令域分析的进给轴装配故障的诊断方法 |
CN107168242B (zh) * | 2017-05-05 | 2019-03-05 | 大连理工大学 | 一种数控机床半闭环进给轴的原点热漂移误差测试及建模方法 |
-
2018
- 2018-05-29 CN CN201810533269.4A patent/CN108717287B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6029098A (en) * | 1996-07-24 | 2000-02-22 | Toshiba Kikai Kabushiki Kaisha | Numerical controller for automatic axial feed |
CN103048949A (zh) * | 2012-12-22 | 2013-04-17 | 西安交通大学 | 一种数控机床伺服系统摩擦误差补偿方法 |
CN106774166A (zh) * | 2017-01-24 | 2017-05-31 | 西安科技大学 | 一种数控机床摩擦误差形态预测及形态特征参数确定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108717287A (zh) | 2018-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108907888B (zh) | 半闭环控制下数控机床进给系统换向误差峰值预测方法 | |
CN108717287B (zh) | 半闭环控制方式下数控机床进给系统摩擦误差峰值预测方法 | |
US5895181A (en) | Dynamic error correcting system in a numerically controlled machine tool | |
JP5143005B2 (ja) | 機械における加工プロセスの最適化のための方法 | |
CN103123477B (zh) | 一种基于电机和机床位置双反馈的轴运动控制方法 | |
US7366576B2 (en) | Position control device and position control method for machine tools | |
JP5657633B2 (ja) | 移動体が反転するときの位置誤差を補正するサーボ制御装置 | |
CN111687652A (zh) | 握持力调整装置以及握持力调整系统 | |
US20160039014A1 (en) | Machine tool | |
JP2010120150A (ja) | 工作機械の熱変形補正のための推定方法 | |
JP2010244256A (ja) | 干渉チェック機能を有する数値制御装置 | |
JP4796936B2 (ja) | 加工制御装置 | |
US9804583B2 (en) | Numerical control device | |
EP1288751A2 (en) | Control method and control apparatus for feed system | |
CN114274047B (zh) | 一种基于力感知测量的高效精密研抛轨迹优化方法 | |
CN112236729B (zh) | 数控装置 | |
KR100871456B1 (ko) | 5축 밀링가공 시간 예측방법 | |
JP7343319B2 (ja) | サーボ制御装置 | |
JP3997201B2 (ja) | 産業用ロボットおよびその異常判断方法 | |
CN109290843B (zh) | 精密数控机床倾斜进给系统换向误差峰值预测方法 | |
CN114029787B (zh) | 一种半闭环控制方式下数控机床进给系统弹性变形补偿方法 | |
CN114281018A (zh) | 一种数控机床倾斜进给系统反向跃冲误差峰值离线预测方法 | |
US11338406B2 (en) | Estimating device and estimating system | |
CN109884982B (zh) | 数值控制装置 | |
JP2021002194A (ja) | 数値制御装置、cnc工作機械、数値制御方法及び数値制御用プログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |