CN112817274A - Machine tool acceleration and deceleration time optimization method and system based on load inertia - Google Patents
Machine tool acceleration and deceleration time optimization method and system based on load inertia Download PDFInfo
- Publication number
- CN112817274A CN112817274A CN202011592795.1A CN202011592795A CN112817274A CN 112817274 A CN112817274 A CN 112817274A CN 202011592795 A CN202011592795 A CN 202011592795A CN 112817274 A CN112817274 A CN 112817274A
- Authority
- CN
- China
- Prior art keywords
- acceleration
- load inertia
- deceleration time
- machine tool
- deceleration
- 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.)
- Granted
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/416—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 of velocity, acceleration or deceleration
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
本发明公开了一种基于负载惯量的机床加减速时间优化方法及系统,属于数控加工技术领域,方法包括:通过解析加工的G代码,识别出G代码中非切削区间和切削区间;利用初始加减速时间参数,对首个工件进行加工,并获取加工过程中的状态数据,以辨识得到每行G代码对应的初始负载惯量;基于所述初始负载惯量和机床运动控制策略,调整每行G代码对应的加减速时间参数,得到第一次优化后的加减速时间参数表;利用所述第一次优化后的加减速时间参数表,对下一个工件进行加工,得到第二次优化后的加减速时间参数表;直至加工完最后一个工件。如此,本发明能对数控系统加减速时间参数快速优化,最大程度实现数控系统加工质量与效率的提升。
The invention discloses a method and system for optimizing acceleration and deceleration time of a machine tool based on load inertia, belonging to the technical field of numerical control machining. Deceleration time parameter, process the first workpiece, and obtain state data during processing to identify the initial load inertia corresponding to each line of G code; based on the initial load inertia and machine tool motion control strategy, adjust each line of G code Corresponding acceleration and deceleration time parameters, obtain the acceleration and deceleration time parameter table after the first optimization; use the acceleration and deceleration time parameter table after the first optimization to process the next workpiece, and obtain the acceleration and deceleration time after the second optimization. Deceleration time parameter table; until the last workpiece is machined. In this way, the present invention can rapidly optimize the acceleration and deceleration time parameters of the numerical control system, thereby realizing the improvement of the processing quality and efficiency of the numerical control system to the greatest extent.
Description
技术领域technical field
本发明属于数控加工技术领域,更具体地,涉及一种基于负载惯量的机床加减速时间优化方法及系统。The invention belongs to the technical field of numerical control machining, and more particularly relates to a method and system for optimizing acceleration and deceleration time of a machine tool based on load inertia.
背景技术Background technique
制造业是发展国民经济的根本,而数控机床的发展是制造业发展的基础。随着制造业对高档数控机床的性能要求不断提高,尤其是要求机床在保证零件加工质量的前提下,提高生产效率、降低生产成本,机床进给系统正朝着高速、精密方向发展。高速主轴也必须配合高速进给系统,才能充分发挥高速切削的众多优点。The manufacturing industry is the foundation of the development of the national economy, and the development of CNC machine tools is the foundation of the development of the manufacturing industry. With the continuous improvement of the performance requirements of high-end CNC machine tools in the manufacturing industry, especially the requirements for machine tools to improve production efficiency and reduce production costs on the premise of ensuring the quality of parts processing, the machine tool feed system is developing in the direction of high speed and precision. The high-speed spindle must also cooperate with the high-speed feed system in order to give full play to the many advantages of high-speed cutting.
加减速时间参数决定了机床轴移动的加速度(加加速度)大小,加减速参数越大,加速度(加加速度)就越小,机床轴到达指定速度的时间就越长。为了实现高速高精加工的要求,合理的设置不同加工进给状态下的加减速时间常数变的很有必要。合理的加减速时间常数,能够防止主轴的高速旋转时的过冲、缩短低速旋转时的加工时间。同时,加减速时间常数也会受到机床进给系统负载惯量的影响。当负载惯量与加减速时间常数不匹配时,机床在高速度和加速度状态下进行频繁的加减速会对机床造成冲击,激发机床进给系统的振动模态,影响工件加工质量。The acceleration/deceleration time parameter determines the acceleration (jerk) of the machine axis movement. The larger the acceleration/deceleration parameter, the smaller the acceleration (jerk), and the longer the time it takes for the machine axis to reach the specified speed. In order to achieve the requirements of high-speed and high-precision machining, it is necessary to reasonably set the acceleration and deceleration time constants under different machining feed states. Reasonable acceleration and deceleration time constants can prevent overshoot during high-speed rotation of the spindle and shorten the machining time during low-speed rotation. At the same time, the acceleration and deceleration time constants are also affected by the load inertia of the machine tool feed system. When the load inertia does not match the acceleration and deceleration time constants, the frequent acceleration and deceleration of the machine tool in the high speed and acceleration state will impact the machine tool, stimulate the vibration mode of the machine tool feed system, and affect the machining quality of the workpiece.
为了提高生产效率、缩短不必要的加工时间、提高产品质量,加工工艺参数优化得到广泛的关注和研究,但加工工艺参数优化的优化变量一般为加工中的主轴转速、进给速度、背吃刀量、铣削宽度等工艺参数。加减速时间常数往往不在其考虑范围之内。同时,对机床进给系统进行惯量辨识,在此基础上进行的参数自整定,也仅仅聚焦在对伺服参数的调整,忽视了加减速时间常数的优化。In order to improve production efficiency, shorten unnecessary processing time, and improve product quality, the optimization of machining process parameters has received extensive attention and research. Quantity, milling width and other process parameters. Acceleration and deceleration time constants are often out of consideration. At the same time, the inertia identification of the machine tool feed system, and the parameter self-tuning based on this, only focuses on the adjustment of the servo parameters, ignoring the optimization of the acceleration and deceleration time constants.
机床加工时,在不同的进给速度、背吃刀量、加工材料等因素的影响下,其进给系统的负载惯量时不同的,因此要通过对机床进行惯量辨识,识别不同情况下的系统惯量,根据负载惯量调整系统的加减速时间常数。目前而言,关于根据负载惯量优化加减速时间常数方面的研究和成果还非常薄弱。During machine tool processing, under the influence of factors such as different feed speeds, back-feeding amount, processing materials, etc., the load inertia of the feed system is different. Therefore, it is necessary to identify the inertia of the machine tool to identify the system under different conditions. Inertia, adjust the acceleration and deceleration time constant of the system according to the load inertia. At present, the research and achievements on optimizing the acceleration and deceleration time constants according to the load inertia are still very weak.
发明内容SUMMARY OF THE INVENTION
针对现有技术的以上缺陷或改进需求,本发明提供了一种基于负载惯量的机床加减速时间优化方法及系统,通过实时辨识各加工状态下的系统负载惯量,并根据负载惯量匹配合适的加减速时间参数,实现对机床轴加减速时间参数的优化。In view of the above defects or improvement needs of the prior art, the present invention provides a method and system for optimizing the acceleration and deceleration time of a machine tool based on the load inertia. Deceleration time parameters, realize the optimization of the acceleration and deceleration time parameters of the machine axis.
为实现上述目的,一方面,本发明提供了一种基于负载惯量的机床加减速时间优化方法,包括以下步骤:In order to achieve the above object, on the one hand, the present invention provides a method for optimizing the acceleration and deceleration time of a machine tool based on the load inertia, including the following steps:
S1、通过解析加工的G代码,识别出G代码中非切削区间和切削区间;S1. Identify the non-cutting interval and the cutting interval in the G code by analyzing the processed G code;
S2、利用初始加减速时间参数,对首个工件进行加工,并获取加工过程中的状态数据,以辨识得到每行G代码对应的初始负载惯量;S2. Use the initial acceleration and deceleration time parameters to process the first workpiece, and obtain the status data during the processing to identify and obtain the initial load inertia corresponding to each line of G code;
S3、基于所述初始负载惯量和机床运动控制策略,调整每行G代码对应的加减速时间参数,得到第一次优化后的加减速时间参数表;其中,所述非切削区间对应的加减速时间参数不进行调整;S3. Based on the initial load inertia and the machine tool motion control strategy, adjust the acceleration and deceleration time parameters corresponding to each line of G code, and obtain the first optimized acceleration and deceleration time parameter table; wherein, the acceleration and deceleration corresponding to the non-cutting interval The time parameter is not adjusted;
S4、利用所述第一次优化后的加减速时间参数表,对下一个工件进行加工,重复步骤S2和S3,得到第二次优化后的加减速时间参数表;直至加工完最后一个工件。S4. Use the first optimized acceleration and deceleration time parameter table to process the next workpiece, and repeat steps S2 and S3 to obtain the second optimized acceleration and deceleration time parameter table until the last workpiece is processed.
进一步地,辨识负载惯量具体包括:Further, identifying the load inertia specifically includes:
(1)计算k时刻的输出估计值:其中,为k时刻的机床速度,为k时刻的电磁转矩,T为转置符号,为k-1时刻的负载惯量估计值;(1) Calculate the estimated output value at time k: in, is the machine tool speed at time k, is the electromagnetic torque at time k, T is the transpose sign, is the estimated value of the load inertia at time k-1;
(2)计算误差其中,y(k)为系统实际输出;将所述误差e(k)反馈到辨识算法中,通过修正项对上一时刻的负载惯量估计值进行修正,从而得到k时刻的负载惯量估计值 (2) Calculation error Among them, y(k) is the actual output of the system; the error e(k) is fed back to the identification algorithm, and the estimated value of the load inertia at the previous moment is estimated by the correction term Correction is made to obtain the estimated value of load inertia at time k
(3)利用k时刻的负载惯量估计值估算出k时刻的负载转矩TL(k),然后利用TL(k)、估算k+1时刻的负载惯量估计值 (3) Use the estimated value of load inertia at time k Estimate the load torque T L (k) at time k, and then use T L (k), Estimate the estimated value of the load inertia at time k+1
(4)循环迭代,直至对应的准则函数取最小值。(4) Loop iteration until the corresponding criterion function takes the minimum value.
进一步地,调整后的加减速时间参数表示为:Further, the adjusted acceleration and deceleration time parameters are expressed as:
其中,t为基准的初始加减速时间参数;J为机床非加工状态的负载惯量;J′是机床运行过程中辨识的负载惯量;m为调节比例系数,与机床运动控制策略、加工指令有关。Among them, t is the initial acceleration and deceleration time parameter of the benchmark; J is the load inertia of the machine tool in the non-processing state; J ' is the load inertia identified during the operation of the machine tool; m is the adjustment proportional coefficient, which is related to the machine tool motion control strategy and machining instructions.
进一步地,所述加工过程中的状态数据包括机床的电流、速度。Further, the state data during the processing includes the current and speed of the machine tool.
进一步地,根据加工要求和加工条件选择机床运动控制策略,所述机床运动控制策略包括直线加减速、指数加减速、S形曲线加减速。Further, a machine tool motion control strategy is selected according to processing requirements and processing conditions, and the machine tool motion control strategy includes linear acceleration and deceleration, exponential acceleration and deceleration, and S-curve acceleration and deceleration.
另一方面,本发明提供了一种基于负载惯量的机床加减速时间优化系统,包括:On the other hand, the present invention provides a system for optimizing the acceleration and deceleration time of a machine tool based on load inertia, including:
解析模块,用于通过解析加工的G代码,识别出G代码中非切削区间和切削区间;The analysis module is used to identify the non-cutting interval and the cutting interval in the G code by analyzing the processed G code;
辨识模块,用于利用初始加减速时间参数,对首个工件进行加工,并获取加工过程中的状态数据,以辨识得到每行G代码对应的初始负载惯量;The identification module is used to process the first workpiece by using the initial acceleration and deceleration time parameters, and obtain the status data during the processing, so as to identify the initial load inertia corresponding to each line of G code;
优化模块,用于基于所述初始负载惯量和机床运动控制策略,调整每行G代码对应的加减速时间参数,得到第一次优化后的加减速时间参数表;其中,所述非切削区间对应的加减速时间参数不进行调整;The optimization module is used to adjust the acceleration and deceleration time parameters corresponding to each line of G codes based on the initial load inertia and the machine tool motion control strategy, and obtain the acceleration and deceleration time parameter table after the first optimization; wherein, the non-cutting interval corresponds to The parameter of acceleration and deceleration time is not adjusted;
迭代模块,用于利用所述第一次优化后的加减速时间参数表,对下一个工件进行加工,重复执行所述辨识模块和优化模块的操作,得到第二次优化后的加减速时间参数表;直至加工完最后一个工件。The iteration module is used to process the next workpiece by using the acceleration and deceleration time parameter table after the first optimization, and repeat the operations of the identification module and the optimization module to obtain the acceleration and deceleration time parameters after the second optimization. table; until the last workpiece is machined.
进一步地,所述辨识模块具体用于,Further, the identification module is specifically used for,
(1)计算k时刻的输出估计值:其中,为k时刻的机床速度,为k时刻的电磁转矩,T为转置符号,为k-1时刻的负载惯量估计值;(1) Calculate the estimated output value at time k: in, is the machine tool speed at time k, is the electromagnetic torque at time k, T is the transpose sign, is the estimated value of the load inertia at time k-1;
(2)计算误差其中,y(k)为系统实际输出;将所述误差e(k)反馈到辨识算法中,通过修正项对上一时刻的负载惯量估计值进行修正,从而得到k时刻的负载惯量估计值 (2) Calculation error Among them, y(k) is the actual output of the system; the error e(k) is fed back to the identification algorithm, and the estimated value of the load inertia at the previous moment is estimated by the correction term Correction is made to obtain the estimated value of load inertia at time k
(3)利用k时刻的负载惯量估计值估算出k时刻的负载转矩TL(k),然后利用TL(k)、估算k+1时刻的负载惯量估计值 (3) Use the estimated value of load inertia at time k Estimate the load torque T L (k) at time k, and then use T L (k), Estimate the estimated value of the load inertia at time k+1
(4)循环迭代,直至对应的准则函数取最小值。(4) Loop iteration until the corresponding criterion function takes the minimum value.
进一步地,调整后的加减速时间参数表示为:Further, the adjusted acceleration and deceleration time parameters are expressed as:
其中,t为基准的初始加减速时间参数;J为机床非加工状态的负载惯量;J′是机床运行过程中辨识的负载惯量;m为调节比例系数,与机床运动控制策略、加工指令有关。Among them, t is the initial acceleration and deceleration time parameter of the benchmark; J is the load inertia of the machine tool in the non-processing state; J ' is the load inertia identified during the operation of the machine tool; m is the adjustment proportional coefficient, which is related to the machine tool motion control strategy and machining instructions.
进一步地,所述加工过程中的状态数据包括机床的电流、速度。Further, the state data during the processing includes the current and speed of the machine tool.
进一步地,根据加工要求和加工条件选择机床运动控制策略,所述机床运动控制策略包括直线加减速、指数加减速、S形曲线加减速。Further, a machine tool motion control strategy is selected according to processing requirements and processing conditions, and the machine tool motion control strategy includes linear acceleration and deceleration, exponential acceleration and deceleration, and S-curve acceleration and deceleration.
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:In general, compared with the prior art, the above technical solutions conceived by the present invention can achieve the following beneficial effects:
(1)本发明以指令域为基础,通过实时辨识各加工状态下的系统负载惯量,并根据负载惯量匹配合适的加减速时间参数,实现对机床当前G代码行的加减速时间参数的优化;在下一轮加工时,当G代码运行到该行时,自动调取优化后的参数值。如此,本发明能对数控系统加减速时间参数快速优化,最大程度实现数控系统加工质量与效率的提升。(1) The present invention is based on the command domain, realizes the optimization of the acceleration and deceleration time parameters of the current G code line of the machine tool by identifying the system load inertia in real time under each processing state, and matching the appropriate acceleration and deceleration time parameters according to the load inertia; In the next round of processing, when the G code runs to this line, the optimized parameter values are automatically retrieved. In this way, the present invention can rapidly optimize the acceleration and deceleration time parameters of the numerical control system, thereby realizing the improvement of the processing quality and efficiency of the numerical control system to the greatest extent.
(2)本发明提供了集成于数控系统中的操作人机交互界面,该界面提供基于指令序列的加工状态信息曲线图,加工过程中实时采集当前状态数据,根据加工数据在线辨识机床的负载惯量,同时绘制指令序列曲线图包括波状和柱状两种形式。其中波状图是负载惯量辨识结果,负载惯量的辨识结果随着辨识周期的增加趋于稳定,用波状图便于观察负载惯量辨识的结果;柱状图是优化后的加减速时间参数,其显示时序略低于波状图,柱状图以加工G代码程序的行为观察轴,对每行的加减速时间参数进行优化,把优化后的结果时时显示,机床轴运动G代码行与其他G代码行的柱状图进行区分,对于轴运动行的加减速时间才需要进行优化调整。如此,本发明与数控系统高度无缝集成,所有数据的交互均在数据系统内部数据进行,参数优化与数控系统周期同步,避免了繁琐、不变的手工调整,从而极大的缩短加减速时间参数优化周期,保证参数优化的及时性、有效性。(2) The present invention provides an operating man-machine interaction interface integrated in the numerical control system, which provides a processing state information graph based on an instruction sequence, collects the current state data in real time during the processing, and identifies the load inertia of the machine tool online according to the processing data. , and draws the instruction sequence curve graph at the same time, including two forms of wave and column. The waveform diagram is the identification result of the load inertia, and the identification result of the load inertia tends to be stable with the increase of the identification period. The waveform diagram is convenient to observe the result of the identification of the load inertia. Below the waveform graph, the bar graph observes the axis with the behavior of processing the G code program, optimizes the acceleration and deceleration time parameters of each line, and displays the optimized result from time to time. The bar graph of the machine axis motion G code line and other G code lines To distinguish, the acceleration and deceleration time of the axis motion line only needs to be optimized and adjusted. In this way, the present invention is highly seamlessly integrated with the numerical control system, all data interactions are carried out in the internal data of the data system, and the parameter optimization is synchronized with the numerical control system cycle, which avoids tedious and constant manual adjustment, thereby greatly shortening the acceleration and deceleration time. The parameter optimization cycle ensures the timeliness and effectiveness of parameter optimization.
附图说明Description of drawings
图1为本发明提供的一种基于负载惯量的机床加减速时间优化方法流程图;1 is a flowchart of a method for optimizing acceleration and deceleration time of a machine tool based on load inertia provided by the present invention;
图2为三种加减速控制策略速度变化曲线图;其中,(a)、(b)、(c)分别为直线型加减速、指数型加减速、S形(铃型)加减速;T为从0加速到一指定速度或从一指定速度减速到0的时间;S型加减速中,加速过程先是变加速,加速度从0开始增加,增加到以稳定值,开始匀速加速,然后再变减速,T2为加速度变化(增加或减小)时间,T1是从开始加速到匀加速阶段结束时间,由设定的速度和T1推算得到;Figure 2 shows the speed change curves of three acceleration and deceleration control strategies; among them, (a), (b), (c) are linear acceleration and deceleration, exponential acceleration and deceleration, and S-shaped (bell-shaped) acceleration and deceleration; T is Time to accelerate from 0 to a specified speed or decelerate from a specified speed to 0; in S-type acceleration and deceleration, the acceleration process first changes to acceleration, the acceleration starts to increase from 0, increases to a stable value, starts to accelerate at a constant speed, and then changes to deceleration , T2 is the acceleration change (increase or decrease) time, T1 is the time from the start of acceleration to the end of the uniform acceleration stage, calculated from the set speed and T1;
图3为本发明的负载惯量辨识结果曲线波形显示图;其中,N1、N2、N3……Nn表示G代码程序行号;纵坐标Y轴表示辨识出的惯量结果,在每一轴移动G代码行进行惯量辨识,随着采集的数据,惯量辨识的结果逐渐收敛,辨识结果稳定;在非轴移动行,不对机床进给系统进行惯量辨识,辨识结果曲线为一水平线;Fig. 3 is the curve waveform display diagram of the load inertia identification result of the present invention; wherein, N1, N2, N3... Inertia identification is carried out on the line, and with the collected data, the inertia identification results gradually converge, and the identification results are stable; in non-axis movement, inertia identification is not performed on the machine tool feed system, and the identification result curve is a horizontal line;
图4为本发明的基于负载惯量加减速时间参数优化曲线柱状显示图;其中,N1、N2、N3……Nn表示G代码程序行号,a、b、b……表示对应行优化后的加减速时间参数,浅色的柱状图表示该行程序是轴移动程序需要修改加减速时间参数,深色的柱状图表示是非轴移动行程序,不需要对加减速时间参数进行修改;纵坐标Y轴代表当前加减速时间参数值,每个柱状图的高度由该G代码行辨识出的负载惯量和选择的加减速控制策略决定;Fig. 4 is a columnar display diagram of the optimization curve of acceleration and deceleration time parameters based on the load inertia of the present invention; wherein, N1, N2, N3... The deceleration time parameter, the light bar chart indicates that the line program is an axis movement program and needs to modify the acceleration and deceleration time parameters, and the dark bar chart indicates that it is a non-axis movement line program, and the acceleration and deceleration time parameters do not need to be modified; the ordinate Y axis Represents the current acceleration/deceleration time parameter value, and the height of each histogram is determined by the load inertia identified by the G code line and the selected acceleration/deceleration control strategy;
图5为本发明提供的在线修改加减速时间参数过程示意图。FIG. 5 is a schematic diagram of a process of online modification of acceleration and deceleration time parameters provided by the present invention.
具体实施方式Detailed ways
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。In order to make the objectives, technical solutions and advantages of the present invention clearer, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention, but not to limit the present invention. In addition, the technical features involved in the various embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.
参阅图1,本发明提供了一种基于负载惯量的机床加减速时间优化方法,包括以下步骤:Referring to FIG. 1, the present invention provides a method for optimizing the acceleration and deceleration time of a machine tool based on the load inertia, including the following steps:
步骤1:首先通过解析加工的G代码,识别出G代码中非切削区间和切削区间,对于非切削过程如G00,各运动轴主要是完成快移定位的任务,对于相同工装、夹具及零件的情况下静态负载惯量不会发生变化;对于切削过程,由于刀具的切削力的作用,机床轴的动态负载惯量会变发生变化;因此这两种情况需要分别辨识负载惯量,根据不同的负载惯量,优化不同轴的加减速时间常数,例如fanuc系统1620/1621参数。Step 1: First, identify the non-cutting interval and cutting interval in the G code by analyzing the processed G code. For the non-cutting process such as G00, each motion axis mainly completes the task of rapid traverse positioning. In this case, the static load inertia will not change; for the cutting process, due to the cutting force of the tool, the dynamic load inertia of the machine tool axis will change; therefore, the load inertia needs to be identified separately in these two cases, according to different load inertia, Optimize the acceleration and deceleration time constants of different axes, such as fanuc system 1620/1621 parameters.
步骤2:机床运动控制策略即加减速控制策略(如直线加减速(速度连续),指数加减速,S形曲线加减速(加速度连续)等),通常数控系统会根据不同加工需要手工调整加减速策略,例如fanuc系统直线型/铃型加减速。不同的加减速控制策略在加工过程中会对机床造成不同的冲击、振动,影响工件的加工质量。对于加工精度要求不高的场合如粗加工中,或者用于启停、进退刀等辅助运动中采用直线加减速控制,主要保证速度的连续,在速度加减速阶段的起点、终点处加速度会有突变,另外速度的过渡不够平滑,运动精度低。指数加减速相对于直线加减速平滑性好,运动进度高,而且加减速的起点与终点仍然存在加速度突变,具有一定的柔性冲击,对于机床的跟踪响应要求较高,可以用在切削进给或手动进给运动中。S形加减速在任意一点的加速度都是连续变化的,从而避免了柔性冲击,速度的平滑性很好,运动精度高,对机床和系统的要求比较高,适用于高速、高精加工场合。具体的加减速策略的选择,需要依据加工人员的经验根据实际的加工要求、加工条件进行选择。不同的加减速策略对于加减速时间的调整不同,如图2所示,根据系统设定好的加减速控制策略,实现根据惯量变化优化加减速时间参数。Step 2: The machine tool motion control strategy is the acceleration and deceleration control strategy (such as linear acceleration and deceleration (continuous speed), exponential acceleration and deceleration, S-curve acceleration and deceleration (continuous acceleration), etc.), usually the CNC system will manually adjust the acceleration and deceleration according to different processing needs. Strategies such as linear/bell acceleration and deceleration for fanuc systems. Different acceleration and deceleration control strategies will cause different shocks and vibrations to the machine tool during the machining process, which will affect the machining quality of the workpiece. For occasions with low requirements on machining accuracy, such as rough machining, or for auxiliary motions such as start-stop, advance and retreat, etc., linear acceleration and deceleration control is adopted, which mainly ensures the continuity of speed. In addition, the transition of speed is not smooth enough, and the motion accuracy is low. Compared with linear acceleration and deceleration, exponential acceleration and deceleration have better smoothness and high movement progress, and there is still a sudden acceleration of acceleration at the start and end points of acceleration and deceleration, which has a certain flexible impact. During manual feed movement. The acceleration of S-shaped acceleration and deceleration is continuously changed at any point, thus avoiding the flexible impact, the smoothness of the speed is very good, the motion precision is high, the requirements for the machine tool and the system are relatively high, and it is suitable for high-speed and high-precision machining occasions. The choice of the specific acceleration and deceleration strategy needs to be selected according to the experience of the processing personnel and the actual processing requirements and processing conditions. Different acceleration and deceleration strategies adjust the acceleration and deceleration time differently. As shown in Figure 2, according to the acceleration and deceleration control strategy set by the system, the acceleration and deceleration time parameters are optimized according to the inertia change.
步骤3:通过首个工件的加工获取加工过程中的数据,此步骤不修改系统的加减速时间参数。本发明需要先对惯量进行识别,根据惯量对加减速时间参数进行优化。惯量识别首选需要获取机床在切削加工中的状态数据,因此需要试加工一工件,通过此次工件获得初始的惯量和优化后的加减速时间参数。在此工件以后加工的工件可以进行迭代辨识和优化。惯量辨识的结果如图3,在人机界面中实时显示。通过获取的加工过程中的电流数据,根据转矩常数(如果电机不提供转矩常数也可根据最大电流、最大力矩或者堵转电流、堵转力矩推算)计算出轴运动加速过程中的电磁转矩。Step 3: Obtain the data in the processing process through the processing of the first workpiece. This step does not modify the acceleration and deceleration time parameters of the system. The present invention needs to identify the inertia first, and optimize the acceleration and deceleration time parameters according to the inertia. The first choice of inertia identification needs to obtain the state data of the machine tool during cutting, so it is necessary to test a workpiece, and obtain the initial inertia and the optimized acceleration and deceleration time parameters through this workpiece. Parts machined after this part can be iteratively identified and optimized. The result of inertia identification is shown in Figure 3, which is displayed in real time on the man-machine interface. Through the acquired current data during the machining process, according to the torque constant (if the motor does not provide the torque constant, it can also be calculated based on the maximum current, maximum torque or stall current and stall torque) to calculate the electromagnetic rotation during the acceleration of the shaft motion. moment.
伺服系统机械部分动态模型通过微分方程描述如下:The dynamic model of the mechanical part of the servo system is described by differential equations as follows:
其中:Te—为电机的电磁转矩;ωr—电机转速(rad/s);Jm—伺服电机转动惯量;JL—电机转子等效转动惯量;TL—伺服电机受到负载转矩。Among them: T e - the electromagnetic torque of the motor; ω r - the motor speed (rad/s); J m - the moment of inertia of the servo motor; J L - the equivalent moment of inertia of the motor rotor; T L - the load torque of the servo motor .
惯量辨识主要是过程是利用系统可测量的输入输出数据,由上式可以看出惯量辨识中涉及到的可测量的数据为电流、速度等数据,根据某种准则,估测出负载惯量。The main process of inertia identification is to use the measurable input and output data of the system. From the above formula, it can be seen that the measurable data involved in inertia identification are current, speed and other data. According to a certain criterion, the load inertia is estimated.
在离线辨识中,可以通过设计周期性激励信号消除掉微分方程中的TL,但是在加工过程中的在线辨识,机床是在按照实际加工程序运行,TL项无法被消除,只能通过进行迭代的方法求得。首先初始化负载惯量、负载转矩,通过迭代方法逐渐逼近精确值。In offline identification, TL in the differential equation can be eliminated by designing periodic excitation signals, but in online identification during machining, the machine tool is running according to the actual machining program, and the TL term cannot be eliminated, only by performing obtained by an iterative method. First, the load inertia and load torque are initialized, and the precise value is gradually approached by an iterative method.
为了获得转动惯量的估计值采用逐渐逼近的方法,在k时刻,根据k-1时刻对象参数估计值转矩估计值TL(k-1)以及输入的数据(根据电流计算的电磁转矩),来计算出当前时刻的输出估计值(机床速度): To obtain an estimate of the moment of inertia Using the method of gradual approximation, at time k, according to the estimated value of the object parameters at time k-1 Torque estimate TL (k-1) and input data (Electromagnetic torque calculated according to the current), to calculate the estimated output value (machine tool speed) at the current moment:
计算出估计误差y(k)为系统实际输出。然后将输出估计误差反馈到辨识算法中去,通过修正项对上一时刻的估计值进行修正,计算出k时刻的估计值。利用k时刻负载惯量的估计值估算出k时刻的负载转矩TL(k),然后利用TL(k)、去估算k+1时刻的值。Calculate the estimation error y(k) is the actual output of the system. Then the output estimation error is fed back to the identification algorithm, and the estimated value at the previous moment is estimated by the correction term. Correction is made to calculate the estimated value at time k. Using the estimated value of the load inertia at time k Estimate the load torque T L (k) at time k, and then use T L (k), to estimate the value at time k+1.
如此循环迭代下去,直至对应的准则函数取最小值,这时估计出的输出值也已经逼近系统的实际输出值y(k),转动惯量的估计值也基本接近实际的参数值。This loop iterates until the corresponding criterion function takes the minimum value, then the estimated output value The actual output value y(k) of the system has also been approximated, the estimated value of the moment of inertia It is also basically close to the actual parameter value.
关于负载惯量和负载转矩的初始值,可以按照经验估计出一个范围,一般,系统等效转动惯量最小为电机参数手册上转子惯量的出厂值,负载惯量等效值不超过转子惯量的10倍;电机负载等效转矩在空载(非切削加工阶段)时为0,最高不超过电机额定转矩的3倍,因此可以得到一个负载惯量和负载转矩的经验范围[Jmin,Jmax]与[TLmin,TLmax],初始值可在这个是范围内确定。Regarding the initial values of load inertia and load torque, a range can be estimated based on experience. Generally, the minimum equivalent moment of inertia of the system is the factory value of rotor inertia in the motor parameter manual, and the equivalent value of load inertia does not exceed 10 times of rotor inertia. ; The equivalent torque of the motor load is 0 at no-load (non-cutting stage), and the maximum does not exceed 3 times the rated torque of the motor, so an empirical range of load inertia and load torque can be obtained [J min , J max ] and [T Lmin ,T Lmax ], the initial value can be determined within this range.
步骤4:根据辨识的惯量结果和选定的加减速控制策略匹配适合的加减速时间常数。根据机床的伺服参数等设置的初始加减速时间参数为t,根据辨识的惯量调整后的加减速时间参数有:Step 4: Match the appropriate acceleration/deceleration time constant according to the identified inertia result and the selected acceleration/deceleration control strategy. The initial acceleration and deceleration time parameters set according to the servo parameters of the machine tool are t, and the acceleration and deceleration time parameters adjusted according to the identified inertia are as follows:
其中,t是作为基准的初始加减速时间参数;J为机床非加工状态的惯量;J′是机床运行过程中辨识的负载惯量;m是与加减速控制策略、加工指令速度等有关的调节比例系数,根据加工经验得出。Among them, t is the initial acceleration and deceleration time parameter as the reference; J is the inertia of the machine tool in the non-processing state; J ′ is the load inertia identified during the operation of the machine tool; m is the adjustment ratio related to the acceleration and deceleration control strategy, processing command speed, etc. The coefficient is obtained according to the processing experience.
生成每一行G代码对应的优化后的加减速时间参数和机床轴的信息。加减速时间常数的优化结果的产生时序落后于惯量辨识结果一行G代码。因为需要先获得该行的负载惯量辨识结果才能进一步的匹配优化加减速时间常数。优化结果如图4所示,优化结果的显示在时序上落后于辨识结果一个G代码行。Generate the optimized acceleration/deceleration time parameters and machine axis information corresponding to each line of G code. The generation sequence of the optimization result of the acceleration and deceleration time constant lags behind the inertia identification result by one line of G code. Because the load inertia identification result of this row needs to be obtained first, the acceleration and deceleration time constants can be further matched and optimized. The optimization results are shown in Figure 4. The display of the optimization results is one G code line behind the identification results in timing.
步骤5:在线修改加减速时间参数。在上一个工件完整的加工过程中生成的优化结果,会保存在一个参数表中,在下一个工件加工时,会根据加工的进程,对加减速时间参数进行修改,如图5。Step 5: Modify the acceleration and deceleration time parameters online. The optimization results generated during the complete machining process of the previous workpiece will be saved in a parameter table. When the next workpiece is machined, the acceleration and deceleration time parameters will be modified according to the machining process, as shown in Figure 5.
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。Those skilled in the art can easily understand that the above are only preferred embodiments of the present invention, and are not intended to limit the present invention. Any modifications, equivalent replacements and improvements made within the spirit and principles of the present invention, etc., All should be included within the protection scope of the present invention.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011592795.1A CN112817274B (en) | 2020-12-29 | 2020-12-29 | A method and system for optimizing acceleration and deceleration time of machine tool based on load inertia |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011592795.1A CN112817274B (en) | 2020-12-29 | 2020-12-29 | A method and system for optimizing acceleration and deceleration time of machine tool based on load inertia |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112817274A true CN112817274A (en) | 2021-05-18 |
CN112817274B CN112817274B (en) | 2022-04-01 |
Family
ID=75855764
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011592795.1A Active CN112817274B (en) | 2020-12-29 | 2020-12-29 | A method and system for optimizing acceleration and deceleration time of machine tool based on load inertia |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112817274B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113359619A (en) * | 2021-06-09 | 2021-09-07 | 意特利(上海)科技有限公司 | Adaptive control method for rotary table inertia of vertical-horizontal switching five-axis machine tool |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPWO2017077607A1 (en) * | 2015-11-04 | 2017-11-02 | 三菱電機株式会社 | Numerical controller |
CN109085801A (en) * | 2017-06-14 | 2018-12-25 | 发那科株式会社 | The control device of motor |
CN109085802A (en) * | 2017-06-14 | 2018-12-25 | 发那科株式会社 | The control device of motor |
-
2020
- 2020-12-29 CN CN202011592795.1A patent/CN112817274B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPWO2017077607A1 (en) * | 2015-11-04 | 2017-11-02 | 三菱電機株式会社 | Numerical controller |
CN109085801A (en) * | 2017-06-14 | 2018-12-25 | 发那科株式会社 | The control device of motor |
CN109085802A (en) * | 2017-06-14 | 2018-12-25 | 发那科株式会社 | The control device of motor |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113359619A (en) * | 2021-06-09 | 2021-09-07 | 意特利(上海)科技有限公司 | Adaptive control method for rotary table inertia of vertical-horizontal switching five-axis machine tool |
CN113359619B (en) * | 2021-06-09 | 2022-09-16 | 意特利(上海)科技有限公司 | Adaptive control method for rotary table inertia of vertical-horizontal switching five-axis machine tool |
Also Published As
Publication number | Publication date |
---|---|
CN112817274B (en) | 2022-04-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102821910B (en) | Tool path generating apparatus, tool path computational methods and tool path generator | |
JP6868161B1 (en) | Information processing equipment, information processing methods and information processing programs | |
CN109901518B (en) | An acceleration and deceleration speed planning method for CNC machine tools under constant force constraints | |
US20100138018A1 (en) | Method for producing a parts program | |
CN105814503A (en) | Numerical control device and numerical control method | |
CN105229543A (en) | Numerical control device | |
JP6740199B2 (en) | Numerical control device, CNC machine tool, numerical control method, and numerical control program | |
CN102081376A (en) | Machining load control system based on instruction sequence optimization | |
CN112817274B (en) | A method and system for optimizing acceleration and deceleration time of machine tool based on load inertia | |
CN112859754A (en) | Machine tool machining control method, machine tool machining control device, storage medium, edge device and server | |
CN112916883B (en) | A Method for Predicting the Deformation of the Allowed Tool in Axle Turning | |
JP2007179314A (en) | Machine tool and program conversion method therefor | |
US11003161B2 (en) | Numerical controller | |
US10444721B2 (en) | Simulation apparatus, program generating device, controller, and display method for computer | |
CN113874799A (en) | NC program generation system and NC program generation method | |
CN109648387B (en) | Control device | |
Kim et al. | Integrated design methodology of ball-screw driven servomechanisms with discrete controllers. Part II: Formulation and synthesis of the integrated design | |
JP7636549B2 (en) | Travel route determination device and computer program | |
CN116257020A (en) | Control method and device of turning and milling composite control system | |
WO2023157244A1 (en) | Machining time prediction device and machining time prediction method | |
CN116449770B (en) | Machining method, device and equipment of numerical control machine tool and computer storage medium | |
JP6587763B1 (en) | Computing device | |
JP2022101268A (en) | Information processing device | |
CN115291529B (en) | Numerical control batch machining cutting parameter optimization method responding to cutter wear time-varying characteristic | |
Wang et al. | Research and implementation of high-precision biaxial tracking control system based on NURBS interpolator |
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 |