CN108320049A - 数控车床多工位回转刀架自动换刀能耗准确预测方法 - Google Patents
数控车床多工位回转刀架自动换刀能耗准确预测方法 Download PDFInfo
- Publication number
- CN108320049A CN108320049A CN201810024999.1A CN201810024999A CN108320049A CN 108320049 A CN108320049 A CN 108320049A CN 201810024999 A CN201810024999 A CN 201810024999A CN 108320049 A CN108320049 A CN 108320049A
- Authority
- CN
- China
- Prior art keywords
- energy consumption
- tool changer
- automatic tool
- lathe
- indicate
- 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
- 238000005265 energy consumption Methods 0.000 title claims abstract description 134
- 238000000034 method Methods 0.000 title claims abstract description 87
- 230000001052 transient effect Effects 0.000 claims abstract description 27
- 238000004364 calculation method Methods 0.000 claims description 20
- 238000005259 measurement Methods 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 5
- 230000001186 cumulative effect Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 4
- 238000003754 machining Methods 0.000 description 4
- 239000000571 coke Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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/19—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 positioning or contouring control systems, e.g. to control position from one programmed point to another or to control movement along a programmed continuous path
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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/408—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 data handling or data format, e.g. reading, buffering or conversion of data
- G05B19/4086—Coordinate conversions; Other special calculations
-
- 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/406—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 monitoring or safety
- G05B19/4063—Monitoring general control system
-
- 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
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0218—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
- G05B23/0224—Process history based detection method, e.g. whereby history implies the availability of large amounts of data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/04—Manufacturing
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q3/00—Devices holding, supporting, or positioning work or tools, of a kind normally removable from the machine
- B23Q3/155—Arrangements for automatic insertion or removal of tools, e.g. combined with manual handling
- B23Q3/15503—Processes characterized by special sequencing of operations or the like, e.g. for optimizing tool changing time or capacity in tool storage
-
- 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/50—Machine tool, machine tool null till machine tool work handling
- G05B2219/50249—Tool, probe, pen changer
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Manufacturing & Machinery (AREA)
- Automation & Control Theory (AREA)
- Tourism & Hospitality (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Human Computer Interaction (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Numerical Control (AREA)
Abstract
本发明公开了一种数控车床多工位回转刀架自动换刀过程能耗准确预测方法。该方法首先测量多组旋转刀位数下自动换刀持续时间,拟合获得自动换刀持续时间的计算模型。采集并运算得到机床基本模块功率,并基于机床基本模块功率和自动换刀持续时间计算得到机床基本模块能耗。采集并运算得到换刀装置稳态功率,进一步计算得到换刀装置稳态能耗。通过获得数控车床执行自动换刀过程中所有功率峰值引起的能耗,累加获得换刀装置瞬态能耗。基于已获得的机床基本模块能耗、换刀装置稳态能耗、换刀装置瞬态能耗,得到自动换刀过程能耗预测模型,实现数控车床自动换刀过程能耗准确预测,为机床节能优化提供基础模型支持。
Description
技术领域
本发明涉及数控车床自动换刀过程能耗预测领域,尤其是一种数控车床多工位回转刀架自动换刀过程能耗准确预测方法。
背景技术
数控车床在执行机械加工工艺过程中由于加工特征(圆柱面、螺纹、切槽)的不同,需要多次进行自动换刀,在外圆车刀、内孔车刀、切槽刀,螺纹刀等刀具间进行切换。数控车床普遍采用多工位回转刀架进行自动换刀,换刀过程仍然消耗能量且换刀过程频繁发生,因此,数控车床多工位回转刀架自动换刀过程能耗是数控车床加工过程能耗的重要组成该部分。数控车床自动换刀过程能耗预测模型的建立对于提高数控车床机械加工全过程能耗预测精度以及数控车床的能量优化具有重要意义。
数控车床多工位回转刀架自动换刀过程中会产生多次功率峰值,因此,换刀全过程能耗既包括功率稳定部分的能耗,也包括功率峰值引起的能耗。当前仍然缺乏一种精确的预测数控车床多工位回转刀架全过程能耗且综合考虑稳态及瞬态能耗的有效方法。本发明方法通过建立数控车床多工位回转刀架自动换刀过程中机床基本模块能耗、换刀装置稳态能耗、换刀装置瞬态能耗的计算模型,进一步建立自动换刀全过程能耗预测模型。本发明方法不仅考虑数控车床自动换刀过程中的稳态能耗还综合考虑了瞬态能耗,与实际情况更加吻合,能耗预测精度较高,是一种实用的数控车床多工位回转刀架自动换刀过程能耗准确预测方法。
发明内容
本发明的目的在于提供一种在数控车床自动换刀之前能够准确预测不同旋转刀位数下的自动换刀过程能耗的方法。该方法预测准确度高,对于数控车床加工过程的能耗精确预测及优化提供重要的模型和方法支持。
一种数控车床多工位回转刀架自动换刀过程能耗准确预测方法,包括如下步骤:
步骤1,根据数控车床多工位回转刀架的当前刀位号和目标刀位号,确定旋转刀位数。由于数控车床多工位回转刀架普遍采用单向顺序旋转方式,因此旋转刀位数计算公式可表示为:
其中:Δ表示回转刀架旋转刀位数,Tt表示回转刀架目标刀位号,Ti表示回转刀架当前刀位号,Tp表示回转刀架总刀位数。
步骤2,测量获取数控车床多个不同旋转刀位数(Δ1,Δ2,…ΔN)下的自动换刀持续时间值(ttc1,ttc2,…ttcN),其中ΔN表示第N个旋转刀位数,ttcN表示在第N个旋转刀位数下自动换刀持续时间值。
步骤3,根据获取得到多组不同旋转刀位数下的自动换刀持续时间,以旋转刀位数Δ为自变量,以自动换刀持续时间ttc为响应变量进行数据拟合,建立自动换刀持续时间的计算模型,其模型表示为:
其中:ttc表示自动换刀持续时间,Δ表示回转刀架旋转刀位数,α,β表示公式中的系数。
步骤4,通过采集多个机床基本模块功率并求平均值得到机床基本模块功率,其功率计算公式表示为:
其中:Pbasic表示机床基本模块功率,Pbasic_i机床基本模块功率第i个测量值;Nb表示机床基本模块功率测量值总个数。
步骤5,基于已获得的机床基本模块功率和自动换刀持续时间计算获得数控车床执行自动换刀过程中的机床基本模块能耗,其计算公式如下:
Ebasic=Pbasic×ttc
其中:Ebasic表示数控车床执行自动换刀过程中的机床基本模块能耗,Pbasic表示机床基本模块功率,ttc表示自动换刀持续时间。
步骤6,通过采集多个换刀装置稳态功率并求平均值得到数控车床执行自动换刀过程中的换刀装置稳态功率,其计算公式如下:
其中:Pms表示换刀装置稳态功率,Pms,Δ_j表示回转刀架旋转刀位数为Δ时换刀装置稳态功率第j个测量值;NsΔ表示回转刀架旋转刀位数为Δ时换刀装置稳态功率测量值的总个数。
步骤7,基于已获得的换刀装置稳态功率和自动换刀持续时间计算获得数控车床执行自动换刀过程中的换刀装置稳态模块能耗,其计算公式如下:
Ems=Pms×ttc
其中:Ems表示数控车床执行自动换刀过程中的换刀装置稳态能耗,Pms表示换刀装置稳态功率,ttc表示自动换刀持续时间。
步骤8,通过获得数控车床执行自动换刀过程中所有功率峰值引起的能耗,累加获得换刀装置瞬态能耗,其计算公式如下:
其中:Emt表示数控车床执行自动换刀过程中的换刀装置瞬态能耗,Δ表示回转刀架旋转刀位数,Emt,Δ_u表示回转刀架旋转刀位数为Δ时第u个功率峰值引起的能耗,NΔ表示回转刀架旋转刀位数为Δ时功率峰值的总个数。
步骤9,基于已获得的数控车床自动换刀过程中的机床基本模块能耗、换刀装置稳态能耗和换刀装置瞬态能耗,计算得到自动换刀过程能耗,其计算公式如下:
Etc=Ebasic+Ems+Emt
其中:Etc表示数控车床多工位回转刀架自动换刀过程能耗,Ebasic表示数控车床执行自动换刀过程中的机床基本模块能耗,Ems表示数控车床执行自动换刀过程中的换刀装置稳态能耗,Emt表示数控车床执行自动换刀过程中的换刀装置瞬态能耗。
在步骤3中,利用数据分析软件中的“Modified Hyperbola I”函数进行数据拟合。
在步骤6中,回转刀架旋转刀位数为Δ时换刀装置稳态功率第j个测量值的采集方法为:采集旋转刀位数为Δ时自动换刀过程中的第j个机床稳态功率Pmst,Δ_j,进一步减去机床基本模块功率Pbasic得到该转速下主轴旋转功率,表示为:Pms,Δ_j=Pmst,Δ_j-Pbasic。
与现有技术相比,本发明具有如下有益效果:
本发明方法通过建立数控车床自动换刀持续时间计算模型,进一步得到数控车床自动换刀过程中机床基本模块能耗、换刀装置稳态能耗、换刀装置瞬态能耗,上述三部分能耗共同构成数控车床自动换刀过程能耗预测模型。该方法不仅考虑数控车床自动换刀过程中的稳态能耗还综合考虑了瞬态能耗,与实际情况更加吻合,能耗预测精度较高,为建立数控车床加工全过程的能耗精确预测模型,提供重要的模型支持。
本方法对数控车床多工位回转刀架自动换刀过程能耗预测精度高,并且可以方便地推广至其他类型机床(如数控加工中心)的自动换刀过程能耗预测。
附图说明
图1为本发明方法的流程示意图;
图2为自动换刀持续时间拟合图;
图3为实现本发明实施例所采用的功率-能耗采集装置示意图;
图4为数控车床自动换刀过程功率曲线图;
图5为换刀装置瞬态能耗构成示意图。
具体实施方式
现结合实施例及附图对本发明进行详细解释。
本发明提出一种数控车床多工位回转刀架自动换刀过程能耗准确预测方法。本发明方法的流程示意图如图1所示,首先,根据回转刀架当前刀位号和目标刀位号计算得到旋转刀位数,测量多组旋转刀位数下自动换刀持续时间,进一步拟合获得自动换刀持续时间的计算模型。采集并运算得到机床基本模块功率,根据已获得的机床基本模块功率和自动换刀持续时间计算得到机床基本模块能耗。采集并运算得到换刀装置稳态功率,根据已获得的换刀装置稳态功率和自动换刀持续时间计算得到换刀装置稳态能耗。通过获得数控车床执行自动换刀过程中所有功率峰值引起的能耗,累加获得换刀装置瞬态能耗。基于已获得的机床基本模块能耗、换刀装置稳态能耗、换刀装置瞬态能耗,得到自动换刀过程能耗预测模型,实现数控车床自动换刀过程能耗准确预测,为机床节能优化提供基础模型支持。
本发明实施例以数控车床CK6153i自动换刀过程为例,该数控车床使用的回转刀架为四工位回转刀架,旋转方式为单向顺序旋转。采用本发明方法其自动换刀过程能耗预测模型,事前精确预测不同旋转刀位数下的自动换刀能耗。
1.确定旋转刀位数
根据数控车床多工位回转刀架的当前刀位号和目标刀位号,确定旋转刀位数。数控车CK6153i采用的是四工位回转刀架,因此,回转刀架总刀位数Tp=4。旋转刀位数计算公式为其中:Δ表示回转刀架旋转刀位数,Tt表示回转刀架目标刀位号,Ti表示回转刀架当前刀位号,Tp表示回转刀架总刀位数。假设当前刀位号Ti=1,目标刀位号Tt=4时,满足Tt≥Ti,因此,旋转刀位数Δ=Tt-Ti=4-1=3。假设当前刀位号Ti=4,目标刀位号Tt=2时,满足Tt<Ti,因此,旋转刀位数Δ=Tp-|Tt-Ti|=4-|2-4|=2。根据上述方法可以获得任意当前和目标刀位号下的旋转刀位数。
2.获取多个不同旋转刀位数下的自动换刀持续时间
当旋转刀位数Δ=0,Δ=1,Δ=2和Δ=3时,用秒表分别测量自动换刀持续时间ttc。所测量得到的不同旋转刀位数下的自动换刀持续时间如表1所示。
表1
3.拟合建立自动换刀持续时间计算模型
上述获得的多组旋转刀位数下的自动换刀持续时间满足如下公式:
以旋转刀位数Δ为自变量,以自动换刀持续时间ttc为响应变量进数据拟合,数据拟合使用Sigmaplot 13软件实现。可以得到系数值α=7.3766、β=1.5127。拟合结果如图2所示,其中相关系数R-sqr为0.9958,非常接近1,表明该拟合公式能够较好地表示不同旋转刀位数下的自动换刀持续时间值。取系数值小数点后三位,可得到数控车床CK6153i自动换刀持续时间计算公式为:其中,ttc为自动换刀持续时间,单位为秒(s);Δ表示旋转刀位数。
4.获取机床基本模块功率
机床基本模块功率Pbasic通过采集多个机床基本模块功率求平均值得到。功率采集所使用的功率-能耗采集装置如图3所示。三个电压夹分别与数控车床空气开关输出端三根相线连接,用于测量数控车床自动换刀过程的三相电压信号;三个电流钳分别套在空气开关输出端三根相线上,用于测量数控车床自动换刀过程的三相电流信号;电压及电流信号传输至功率分析仪中进行分析处理获得数控车床自动换刀过程的功率及能耗信息,并存储于计算机中的SQL数据库里。该装置与数控车床CK6153i进行连接。启动数控车床CK6153i,不进行任何操作,仅仅使机床基本模块处于运行状态,测量50组机床基本模块功率值,如表2所示。
表2
机床基本模块功率计算公式为其中,Pbasic表示机床基本模块功率,单位为瓦(W);Pbasic_i表示所采集的第i个机床基本模块功率值,单位为瓦(W);Nb为采集的机床基本模块功率值的总数。上述共采集50组机床基本模块功率值,因此Nb=50,将上述采集得到的50组数据,代入公式可以得到数控车床CK6153i的机床基本模块功率为
5.计算机床基本模块能耗
根据已获得的机床基本模块功率和自动换刀持续时间,计算得到数控车床执行自动换刀过程中的机床基本模块能耗,其计算公式为Ebasic=Pbasic×ttc,其中,Ebasic表示数控车床执行自动换刀过程中的机床基本模块能耗,单位为焦(J);Pbasic表示机床基本模块功率,单位为瓦(W);ttc表示自动换刀持续时间,单位为秒(s)。对于数控车床CK6153i,已得到机床基本模块功率为Pbasic=332.1 (W),其自动换刀持续时间公式为当旋转刀位数Δ=1时,自动换刀持续时间为进一步计算得到数控车床执行自动换刀过程中的机床基本模块能耗为Ebasic=332.1×2.9=963.1 (J)。
6.获取换刀装置稳态功率
通过采集多个换刀装置稳态功率并求平均值得到数控车床执行自动换刀过程中的换刀装置稳态功率。如图4所示,换刀装置稳态功率采集方法为:当旋转刀位数为Δ时,采集自动换刀过程中的机床稳态功率Pmst,换刀装置稳态功率值Pms通过采集自动换刀过程中的机床稳态功率Pmst减去机床基本模块功率Pbasic得到。例如:对于数控车床CK6153i,旋转刀位数Δ=1时,采集得到自动换刀过程中的第1个机床稳态功率值Pmst,1_1=395.6 (W),且已得到数控车床CK6153i的机床基本模块功率Pbasic=332.1 (W)。因此,可得到旋转刀位数Δ=1时,换刀装置稳态功率第1个测量值为Pms,1_1=Pmst,1_1-Pbasic=395.6-332.1=63.5(W)。利用同样的方法,可以采集得到旋转刀位数Δ=1时,20组换刀装置稳态功率值,如表3所示。
表3
换刀装置稳态功率计算公式为其中,Pms表示换刀装置稳态功率,单位为瓦(W);Pms,Δ_j表示回转刀架旋转刀位数为Δ时换刀装置稳态功率第j个测量值,单位为瓦(W);NsΔ表示回转刀架旋转刀位数为Δ时换刀装置稳态功率测量值的总个数。上述共采集旋转刀位数Δ=1时,20组换刀装置稳态功率值,因此,Ns1=20,将上述采集得到的20组数据,代入公式可计算得到旋转刀位数Δ=1时,数控车床CK6153i执行自动换刀过程中的换刀装置稳态功率利用同样的方法可以得到数控车床CK6153i,旋转刀位数Δ=2时,换刀装置稳态功率为66.5W,旋转刀位数Δ=3时,换刀装置稳态功率为62.8W。
7.计算换刀装置稳态能耗
根据已获得的换刀装置稳态功率和自动换刀持续时间,计算得到数控车床执行自动换刀过程中的换刀装置稳态能耗,其计算公式为Ems=Pms×ttc,其中,Ems表示数控车床执行自动换刀过程中的换刀装置稳态能耗,单位为焦(J);Pms表示换刀装置稳态功率,单位为瓦(W);ttc表示自动换刀持续时间,单位为秒(s)。对于数控车床CK6153i,当旋转刀位数Δ=1时,已得到换刀装置稳态功率为Pms=63.1 (W),其自动换刀持续时间ttc=2.9 (s)。进一步计算得到旋转刀位数Δ=1时,数控车床CK6153i执行自动换刀过程中的换刀装置稳态能耗为Ems=63.1×2.9=183.0 (J)。
8.计算换刀装置瞬态能耗
数控车床执行自动换刀过程中会出现多次功率峰值,所有功率峰值引起的能耗总和为换刀装置瞬态能耗,其计算公式为其中,Emt表示数控车床执行自动换刀过程中的换刀装置瞬态能耗,单位为焦(J);Δ表示回转刀架旋转刀位数;Emt,Δ_u表示回转刀架旋转刀位数为Δ时第u个功率峰值引起的能耗,单位为焦(J);NΔ表示回转刀架旋转刀位数为Δ时功率峰值的总个数。承接前例,对于数控车床CK6153i,当旋转刀位数Δ=1时,其自动换刀过程功率曲线如图5所示,自动换刀过程中会出3次功率峰值,因此,N1=3,此时换刀装置瞬态能耗值为三次功率峰值引起能耗的总和。根据能耗采集数据可得到旋转刀位数Δ=1时,三次功率峰值引起的能耗分别为Emt,1_1=4.6 (J),Emt,1_2=23.9 (J),Emt,1_2=88.3 (J)。根据公式可计算得到数控车床CK6153i旋转刀位数为1时,换刀装置瞬态能耗为:
9.计算自动换刀过程能耗
如图4所示,数控车床自动换刀过程能耗由该过程中机床基本模块能耗、换刀装置稳态能耗和换刀装置瞬态能耗构成。因此,基于已获得的数控车床自动换刀过程中的机床基本模块能耗、换刀装置稳态能耗和换刀装置瞬态能耗,可计算得到自动换刀过程能耗,其计算公式为:Etc=Ebasic+Ems+Emt。其中,Etc表示数控车床多工位回转刀架自动换刀过程能耗,单位为焦(J);Ebasic表示数控车床执行自动换刀过程中的机床基本模块能耗,单位为焦(J);Ems表示数控车床执行自动换刀过程中的换刀装置稳态能耗,单位为焦(J);Emt表示数控车床执行自动换刀过程中的换刀装置瞬态能耗,单位为焦(J)。承接前例,已得到数控车床CK6153i旋转刀位数Δ=1时,机床基本模块能耗为Ebasic=963.1 (J);换刀装置稳态能耗为Ems=183.0 (J);换刀装置瞬态能耗为Emt=116.8 (J)。因此,计算得到数控车床自动换刀过程能耗为Etc=963.1+183.0+116.8=1262.9 (J)。根据上述同样的方法,可以得到数控车床CK6153i旋转刀位数Δ=2,Δ=3时的自动换刀过程能耗预测值如表4所示。同时,通过图3所示的功率-能耗采集装置采集数控车床CK6153i旋转刀位数为1,2,3时的自动换刀过程能耗测量值,并通过精度计算公式:精度=(1-|Etc预测-Etc测量|/Etc测量)×100%,验证本发明的有效性。精度计算结果如表示所示。
表4
通过上述对比结果发现,应用本发明方法对数控车床多工位回转刀架自动换刀过程能耗的预测精度都在95%以上,说明本发明方法能够较为精确地预测数控车床多工位回转刀架自动换刀过程的能耗。
本发明方法可以用于数控车床多工位回转刀架自动换刀过程的能耗事前准确预测,计算结果可应用于数控车床机械加工全过程的能量预测,提高其预测精度且为数控车床能量优化提供基础。
最后说明的是,以上实施案例仅用以说明本发明的技术方案而非限制,对本发明的技术方案进行修改或者等同替换,而不脱离本发明方法的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (3)
1.一种数控车床多工位回转刀架自动换刀过程能耗准确预测方法,其特征在于,包括如下步骤:
步骤1,根据数控车床多工位回转刀架的当前刀位号和目标刀位号,确定旋转刀位数。由于数控车床多工位回转刀架普遍采用单向顺序旋转方式,因此旋转刀位数计算公式可表示为:
其中:Δ表示回转刀架旋转刀位数,Tt表示回转刀架目标刀位号,Ti表示回转刀架当前刀位号,Tp表示回转刀架总刀位数。
步骤2,测量获取数控车床多个不同旋转刀位数(Δ1,Δ2,…ΔN)下的自动换刀持续时间值(ttc1,ttc2,…ttcN),其中ΔN表示第N个旋转刀位数,ttcN表示在第N个旋转刀位数下自动换刀持续时间值。
步骤3,根据获取得到多组不同旋转刀位数下的自动换刀持续时间,以旋转刀位数Δ为自变量,以自动换刀持续时间ttc为响应变量进行数据拟合,建立自动换刀持续时间的计算模型,其模型表示为:
其中:ttc表示自动换刀持续时间,Δ表示回转刀架旋转刀位数,α,β表示公式中的系数。
步骤4,通过采集多个机床基本模块功率并求平均值得到机床基本模块功率,其功率计算公式表示为:
其中:Pbasic表示机床基本模块功率,Pbasic_i机床基本模块功率第i个测量值;Nb表示机床基本模块功率测量值总个数。
步骤5,基于已获得的机床基本模块功率和自动换刀持续时间计算获得数控车床执行自动换刀过程中的机床基本模块能耗,其计算公式如下:
Ebasic=Pbasic×ttc
其中:Ebasic表示数控车床执行自动换刀过程中的机床基本模块能耗,Pbasic表示机床基本模块功率,ttc表示自动换刀持续时间。
步骤6,通过采集多个换刀装置稳态功率并求平均值得到数控车床执行自动换刀过程中的换刀装置稳态功率,其计算公式如下:
其中:Pms表示换刀装置稳态功率,Pms,Δ_j表示回转刀架旋转刀位数为Δ时换刀装置稳态功率第j个测量值;NsΔ表示回转刀架旋转刀位数为Δ时换刀装置稳态功率测量值的总个数。
步骤7,基于已获得的换刀装置稳态功率和自动换刀持续时间计算获得数控车床执行自动换刀过程中的换刀装置稳态模块能耗,其计算公式如下:
Ems=Pms×ttc
其中:Ems表示数控车床执行自动换刀过程中的换刀装置稳态能耗,Pms表示换刀装置稳态功率,ttc表示自动换刀持续时间。
步骤8,通过获得数控车床执行自动换刀过程中所有功率峰值引起的能耗,累加获得换刀装置瞬态能耗,其计算公式如下:
其中:Emt表示数控车床执行自动换刀过程中的换刀装置瞬态能耗,Δ表示回转刀架旋转刀位数,Emt,Δ_u表示回转刀架旋转刀位数为Δ时第u个功率峰值引起的能耗,NΔ表示回转刀架旋转刀位数为Δ时功率峰值的总个数。
步骤9,基于已获得的数控车床自动换刀过程中的机床基本模块能耗、换刀装置稳态能耗和换刀装置瞬态能耗,计算得到自动换刀过程能耗,其计算公式如下:
Etc=Ebasic+Ems+Emt
其中:Etc表示数控车床多工位回转刀架自动换刀过程能耗,Ebasic表示数控车床执行自动换刀过程中的机床基本模块能耗,Ems表示数控车床执行自动换刀过程中的换刀装置稳态能耗,Emt表示数控车床执行自动换刀过程中的换刀装置瞬态能耗。
2.如权利要求1所述一种数控车床多工位回转刀架自动换刀过程能耗准确预测方法,其特征在于,在步骤3中,利用数据分析软件中的“Modified Hyperbola I”函数进行数据拟合。
3.如权利要求1所述一种数控车床多工位回转刀架自动换刀过程能耗准确预测方法,其特征在于,在步骤6中,回转刀架旋转刀位数为Δ时换刀装置稳态功率第j个测量值的采集方法为:采集旋转刀位数为Δ时自动换刀过程中的第j个机床稳态功率Pmst,Δ_j,进一步减去机床基本模块功率Pbasic得到该转速下主轴旋转功率,表示为:Pms,Δ_j=Pmst,Δ_j-Pbasic。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810024999.1A CN108320049B (zh) | 2018-01-11 | 2018-01-11 | 数控车床多工位回转刀架自动换刀能耗准确预测方法 |
PCT/CN2018/088876 WO2019136906A1 (zh) | 2018-01-11 | 2018-05-29 | 数控车床多工位回转刀架自动换刀能耗准确预测方法 |
US16/483,995 US10691098B2 (en) | 2018-01-11 | 2018-05-29 | Method of accurately predicting energy consumption of automatic tool change for multi-position rotary tool holder of numerical control machine |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810024999.1A CN108320049B (zh) | 2018-01-11 | 2018-01-11 | 数控车床多工位回转刀架自动换刀能耗准确预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108320049A true CN108320049A (zh) | 2018-07-24 |
CN108320049B CN108320049B (zh) | 2019-06-11 |
Family
ID=62893122
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810024999.1A Active CN108320049B (zh) | 2018-01-11 | 2018-01-11 | 数控车床多工位回转刀架自动换刀能耗准确预测方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US10691098B2 (zh) |
CN (1) | CN108320049B (zh) |
WO (1) | WO2019136906A1 (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109933002A (zh) * | 2019-03-28 | 2019-06-25 | 河海大学常州校区 | 一种面向节能的机械加工过程数控机床能耗建模方法 |
CN110560920A (zh) * | 2019-08-22 | 2019-12-13 | 浙江科技学院 | 一种回转件激光切割加工的能耗预测方法 |
CN112230601A (zh) * | 2020-10-21 | 2021-01-15 | 山东科技大学 | 数控机床主轴停转节能临界时间确定及节能方法 |
CN114415595A (zh) * | 2021-11-05 | 2022-04-29 | 山东科技大学 | 一种车削优化方法、系统、计算机设备和存储介质 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113050541B (zh) * | 2021-03-24 | 2023-04-28 | 武汉科技大学 | 一种数控机床能耗与加工状态在线监测系统 |
CN113110288B (zh) * | 2021-04-23 | 2022-09-27 | 重庆大学 | 一种滚齿机床机械加工系统集成优化设计方法 |
CN114895630B (zh) * | 2022-05-31 | 2024-06-21 | 安徽工业大学 | 一种机床热误差预测模型的温度敏感点的选择方法 |
CN115167279B (zh) * | 2022-09-07 | 2022-11-29 | 中科航迈数控软件(深圳)有限公司 | 针对数控机床的能耗预测方法、系统及相关设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104391479A (zh) * | 2014-11-26 | 2015-03-04 | 上海交通大学 | 基于主轴功率监控的机床工作状态的自动统计方法及系统 |
CN105652791A (zh) * | 2015-12-04 | 2016-06-08 | 西安交通大学 | 订单驱动的离散制造过程能耗优化方法 |
CN105929689A (zh) * | 2016-04-22 | 2016-09-07 | 江南大学 | 基于粒子群算法的机床制造系统加工节能优化方法 |
CN107193258A (zh) * | 2017-06-22 | 2017-09-22 | 重庆大学 | 面向能耗的数控加工工艺路线与切削参数优化模型与方法 |
CN107220734A (zh) * | 2017-06-26 | 2017-09-29 | 江南大学 | 基于决策树的数控车床车削过程能耗预测系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011138463A (ja) * | 2010-01-04 | 2011-07-14 | Mitsubishi Electric Corp | 数値制御装置 |
-
2018
- 2018-01-11 CN CN201810024999.1A patent/CN108320049B/zh active Active
- 2018-05-29 US US16/483,995 patent/US10691098B2/en active Active
- 2018-05-29 WO PCT/CN2018/088876 patent/WO2019136906A1/zh active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104391479A (zh) * | 2014-11-26 | 2015-03-04 | 上海交通大学 | 基于主轴功率监控的机床工作状态的自动统计方法及系统 |
CN105652791A (zh) * | 2015-12-04 | 2016-06-08 | 西安交通大学 | 订单驱动的离散制造过程能耗优化方法 |
CN105929689A (zh) * | 2016-04-22 | 2016-09-07 | 江南大学 | 基于粒子群算法的机床制造系统加工节能优化方法 |
CN107193258A (zh) * | 2017-06-22 | 2017-09-22 | 重庆大学 | 面向能耗的数控加工工艺路线与切削参数优化模型与方法 |
CN107220734A (zh) * | 2017-06-26 | 2017-09-29 | 江南大学 | 基于决策树的数控车床车削过程能耗预测系统 |
Non-Patent Citations (1)
Title |
---|
吕景祥: "面向低碳制造的数控机床能量供给建模研究", 《中国博士学位论文全文数据库 工程科技II辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109933002A (zh) * | 2019-03-28 | 2019-06-25 | 河海大学常州校区 | 一种面向节能的机械加工过程数控机床能耗建模方法 |
CN110560920A (zh) * | 2019-08-22 | 2019-12-13 | 浙江科技学院 | 一种回转件激光切割加工的能耗预测方法 |
CN112230601A (zh) * | 2020-10-21 | 2021-01-15 | 山东科技大学 | 数控机床主轴停转节能临界时间确定及节能方法 |
CN114415595A (zh) * | 2021-11-05 | 2022-04-29 | 山东科技大学 | 一种车削优化方法、系统、计算机设备和存储介质 |
CN114415595B (zh) * | 2021-11-05 | 2024-05-10 | 山东科技大学 | 一种车削优化方法、系统、计算机设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN108320049B (zh) | 2019-06-11 |
US20200019135A1 (en) | 2020-01-16 |
WO2019136906A1 (zh) | 2019-07-18 |
US10691098B2 (en) | 2020-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108320049B (zh) | 数控车床多工位回转刀架自动换刀能耗准确预测方法 | |
CN105759719B (zh) | 一种基于无偏估计拆分模型的数控机床热误差预测方法及系统 | |
CN102122146B (zh) | 用于高速精密加工的热误差实时补偿系统及其补偿方法 | |
CN104808584B (zh) | 机床多源能耗系统多信息在线检测系统 | |
CN102621932A (zh) | 一种数控机床服役过程的能量消耗预测方法 | |
CN105607579B (zh) | 一种机床加工智能节能控制方法及系统 | |
CN103235554B (zh) | 一种基于nc代码的数控车床加工工件能耗获取方法 | |
CN108427841A (zh) | 一种数控机床刀具剩余寿命实时预测方法 | |
CN102637014B (zh) | 数控机床加工过程机电主传动系统能量效率获取的方法 | |
CN108133091B (zh) | 一种基于刀具状态建立机床碳排放优化模型的方法 | |
CN107193259A (zh) | 一种基于实时电量的数控机床主传动系统能耗和时间信息获取方法 | |
CN102880771B (zh) | 高速切削加工中工件表面粗糙度的预测方法 | |
CN103231279A (zh) | 一种数控机床切削状态下的机床主轴动态测试装置 | |
CN103941644A (zh) | 一种基于时间参数的数控铣床能耗预测方法 | |
CN112147950B (zh) | 数控侧铣加工系统及其阈值与待监控信号的精确同步方法 | |
CN103676782B (zh) | 数控铣床加工过程中能量效率在线检测方法 | |
CN109753017B (zh) | 数控机床的固有能效要素函数获取系统与获取方法 | |
CN108037734A (zh) | 数控机床钻削过程功率及能耗获取与节能控制方法 | |
CN105955198B (zh) | 基于最小二乘迭代算法的机床工步能耗监测方法 | |
CN103170403A (zh) | 球磨机群监测系统的数据采集监测模块及其监测方法 | |
CN110262392A (zh) | 一种基于虚拟样件的机床固有能效属性测量系统与方法 | |
CN107024907B (zh) | 一种嵌入式全生命周期机床热误差补偿系统及方法 | |
CN208945752U (zh) | 一种数控机床电主轴的热补偿装置 | |
CN109342812A (zh) | 一种机械加工设备能量消耗测量装置 | |
CN111007801B (zh) | 一种基于多维属性状态判断的刀具寿命实时采集方法 |
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 |