CN106990756B - 一种数控机床几何精度在线监测方法 - Google Patents

一种数控机床几何精度在线监测方法 Download PDF

Info

Publication number
CN106990756B
CN106990756B CN201710195633.6A CN201710195633A CN106990756B CN 106990756 B CN106990756 B CN 106990756B CN 201710195633 A CN201710195633 A CN 201710195633A CN 106990756 B CN106990756 B CN 106990756B
Authority
CN
China
Prior art keywords
machine tool
strain
guide rail
along guide
numerically
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
Application number
CN201710195633.6A
Other languages
English (en)
Other versions
CN106990756A (zh
Inventor
王永青
吴嘉锟
刘志松
刘海波
况康
刘阔
盛贤君
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN201710195633.6A priority Critical patent/CN106990756B/zh
Publication of CN106990756A publication Critical patent/CN106990756A/zh
Application granted granted Critical
Publication of CN106990756B publication Critical patent/CN106990756B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical 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/406Numerical 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/4065Monitoring tool breakage, life or condition
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/37Measurements
    • G05B2219/37616Use same monitoring tools to monitor tool and workpiece

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)
  • Arrangements For Transmission Of Measured Signals (AREA)

Abstract

本发明一种数控机床几何精度在线监测方法属于数控机床技术领域,涉及一种数控机床几何精度在线监测方法。该方法中,首先在机床基础件表面粘贴应变片;通过搭建的无线应变监测***,采集并远程传输各离散测点应变数据;将离散应变数据作为边界条件,求解满足双调和函数的应力方程,预测机床基础件表面应变场分布,通过积分得到机床基础件变形状态;再通过端点连线法求出机床基础件直线度。该方法具有测量精度高、适应性好的特点,实现了测量—加工无干扰同步进行。应变数据通过无线检测***进行数据传输,具有占用空间小,为机床维护、故障预紧与自主维修提供了重要数据支持。

Description

一种数控机床几何精度在线监测方法
技术领域
本发明属于数控机床技术领域,特别涉及一种数控机床几何精度在线监测方法。
背景技术
数控机床基础件(如床身、立柱等)的几何精度是机床各项精度基础,例如床身导轨安装面的直线度。在数控机床长期服役过程中,由于受振动、热等因素影响,机床基础件的内部残余内应力和装配应力释放、重平衡显著,进而造成基础件产生复杂的应力变形,导致机床几何精度不可避免地恶化、衰退。机床几何精度变化规律的获取,不仅对于机床精度保持能力评价,而且对于机床在使用中的精度维护和保障措施的施加具有重要意义。实际生产中,机床厂家和用户单位大多采用几何精度精密测量仪器(例如激光干涉仪)定期、间歇式地实施测量。这种检测方法比较复杂、耗时,且需要停机测量,势必给正常生产带来很大压力。特别对于生产线,需要坚决杜绝频繁停机。为此,寻求一种数控机床几何精度在线监测方法,迫在眉睫。
研究表明,多因素影响下的数控机床几何精度在线监测,需满足数据处理能力强、预估精度高、***抗干扰能力强、不影响机床正常运行等基本条件,具有极大的工程挑战性。通过检测基础件上特征点的应变状态,反算结构变形,完成几何精度预估,为数控机床几何精度在线监测提供了可行性。
四川普什宁江机床有限公司在专利CN104483896A中公开了《数控机床丝杠热补偿的实现方法》,结合热补偿数据库数据,利用热补偿控制器,进行丝杠热变形补偿。北京工研精机股份有限公司在专利CN104400563A中公开了《一种机床加工过程变形的在线补偿方法》,通过布置在大型结构内侧的低膨胀系数的测量框架测量机床特征点位置变化,按照机床的结构模型得知整个机床的变形情况以及主轴刀具中心点变化,最后通过修调刀具位置完成变形补偿。
然而,上述方法均未提及一种数控机床几何精度在线监测方法。
发明内容
本发明目的在于克服现有方法不足,针对服役过程的数控机床几何精度长期、在线监测的问题,发明了一种数控机床几何精度在线监测方法。通过建立参数自适应化几何精度反算算法,覆盖任意复杂测点分布形式应力场函数拟合,实用性强;建立应力函数满足双调和方程,防止拟合场微元出现“裂缝”、“嵌入”现象,提高应力场预测精度,实现直线度精确反算;监测***采用多应变片分布测量形式,占用空间小,便于组装,实现测量—加工无干扰同步进行,工况适应性强;应变数据进行无线、实时传输与分析,为机床远程智能维护、故障预警与自主维修提供重要数据依据。
本发明采用的技术方案是一种数控机床几何精度在线监测方法,该方法首先在机床基础件表面粘贴应变片;通过搭建的无线应变监测***,采集并远程传输各离散测点应变数据;将离散应变数据作为边界条件,求解满足双调和函数的应力方程,预测机床基础件表面应变场分布,通过积分得到机床基础件变形状态;再通过端点连线法求出机床基础件直线度。方法的具体步骤如下:
第一步布置应变测点
在机床基础件1的第一侧表面S1、第二侧表面S2的上、下边缘,上表面S3的导轨内侧处的边缘位置粘贴三轴应变片,各表面应变片数量均为k,
其中,[]为取整符号,即[x]为不超过x的最大整数。l为导轨长度。各表面沿导轨方向s相邻应变片间距为d,
第二步采集应变数据
利用数控机床几何精度无线监测***进行应变数据采集。该***是一个分布式的数据采集***,符合“控制中心-分布式采集节点”的构架模式,主要包括应变片组2、无线应变节点3、无线网关4、数据中心服务器5、上位机6组成。无线应变节点3内置电桥电阻与放大调理电路,无线应变节点3采集应变片组2的应变值,并将其无线传输至无线网关4;通过通用无线分组业务GPRS将无线网关4实时监测的测点应变数据传送到数据中心服务器5;上位机6利用数据中心服务器5中的应变数据进行几何精度计算。
第三步反算机床基础件变形
首先,设定平面Si(i=1,2,3)表示为第i个表面,即第一侧表面S1、第二侧表面S2或上表面S3。在平面Si建立的直角坐标系为xSi-oSi-ySi。坐标轴xS1、xS2、xS3正方向沿导轨方向s,坐标轴yS1、yS2正方向沿导轨安装底面法向n,坐标轴yS3正方向沿导轨L1安装侧面法向t。那么,各应变片的坐标集合为 为等效平面Si上编号w的应变片坐标。平面Si沿坐标轴ySi方向变形曲线为:
其中,p按如下公式取值,
其中,E为基础件材料的杨氏模量;ν为泊松比;为坐标集合中最小与最大的纵坐标;xSi为沿着坐标轴xSi方向自变量;为m/n幂系数,组成幂系数数组为零的m/n幂系数组成0元素数组为,
非零的m/n幂系数通过联立求解(5-1),(5-2),(5-3),(5-4)计算:
其中,分别为平面Si上编号w应变片采集的沿xSi方向、ySi方向正应变;为等效平面Si上编号w应变片采集的切应变。
第四步直线度评估
采用端点连线法计算第一侧表面S1、第二侧表面S2与上表面S3变形曲线的平面直线度分别为因此,沿导轨安装底面法向n基础件直线度为:
沿导轨L1安装侧面法向t基础件直线度为:
通过监测***第N(N=1,2,…,+∞)次采集各测点应变值并通过上述方法计算得到的沿导轨安装底面法向n与沿导轨L1安装侧面法向t的基础件直线度为fn(N)、ft(N)。
第五步直线度归一化处理
采用直线度测量仪器测量基础件沿导轨安装底面法向n、沿导轨L1安装侧面法向t直线度测量过程中使监测***中所有应变片平衡清零。一个月后,采用直线度测量仪器再次测量直线度同时,通过监测***首次采集测点应变值并通过计算得到基础件的直线度fn(1)、ft(1)。因此,归一化处理下通过监测***第N次监测的基础件沿导轨安装底面法向n与沿导轨L1安装侧面法向t的修正直线度Fn、Ft为:
通过以上步骤完成数控机床几何精度在线监测。
本发明的有益效果为采用了基于分布式应变片测量方法,实现了不影响正常加工情况下,对数控机床基础件几何精度的实时测量。通过无线信号实现了数据传输,方便对数控机床床身几何精度的远程、长期监控。具有测量精度高、适应性好的特点,实现了测量—加工无干扰同步进行。应变数据通过无线检测***进行数据传输,具有占用空间小,为机床维护、故障预紧与自主维修提供了重要数据支持。
附图说明
附图1-数控机床床身示意图,其中,1-数控机床床身,S1为数控机床床身1的第一侧表面,S2为数控机床床身1的第二侧表面,S3为数控机床床身1的上表面,L1为第一侧表面S1侧导轨,L2为二侧表面S2侧导轨,n方向为导轨L1、L2底部安装面的法向,t方向为导轨L1侧向安装面法向,s为导轨方向,2-应变片组,xS1、yS1、oS1分别为第一侧表面S1的xS1坐标轴、yS2坐标轴与坐标原点,xS2、yS2与oS2分别为第二侧表面S2的xS2坐标轴、yS2坐标轴与坐标原点,xS3、yS3与oS3分别为上表面S3的xS3坐标轴、yS3坐标轴与坐标原点。
附图2-数控机床床身几何精度的监测***结构示意图,其中,1-数控机床床身,2-应变片组,3-无线应变节点,4-无线网关,5-数据中心服务器,6-PC端上位机。
具体实施方式
下面将结合附图以及具体实施例来详细说明本发明。
被测数控机床为立式加工中心床身,导轨L1、L2长1060mm,导轨L1、L2内侧间距770mm,床身高550mm。
第一步 应变测点布置
在数控机床床身1的第一侧表面S1、第二侧表面S2的上、下边缘,上表面S3的两导轨内侧处的边缘均匀粘贴应变片,各表面应变片数量,k=[(1060/250)+3]=3。第一侧表面S1、第二侧表面S2与上表面S3沿导轨方向s相邻应变片间距为,d=(2×1060)/(6-2)=530。第一侧表面S1、第二侧表面S2上粘贴的应变片中部敏感栅沿导轨安装底面法向n粘贴,其中沿导轨安装底面法向n相邻应变片间距500mm;上表面S3上粘贴的导轨L1、L2侧应变片中部敏感栅沿导轨L1安装侧面法向t粘贴,其中沿导轨L1安装侧面法向t相邻应变片间距770mm。
第二步 应变数据采集
应变片组2与无线应变节点3连接,每节点8通道,应变片桥路为1/4桥,架设数据中心服务器5,通过PC端上位机6的应变采集软件完成无线网关4数据目标端口设置与网关绑定,连接IP地址和端口完成PC端上位机6对各无线应变节点3数据的采集。
第三步 机床基础件变形反算
建立第一侧表面S1、第二侧表面S2与上表面S3上坐标系xS1-oS1-yS1、xS2-oS2-yS2、xS3-oS3-yS3。其中坐标轴xS1、xS2、xS3沿导轨方向s,坐标轴yS1、yS2沿导轨安装底面法向n,坐标轴yS3沿导轨L1安装侧面法向t。使各面端部的应变片坐标为(0.01,0.01),保证其余应变片坐标为正。采用激光干涉仪测量床身沿导轨安装底面法向n、沿导轨L1安装侧面法向t直线度测量过程中使监测***中所有应变片平衡清零。机床服役一个月后,首次第一侧表面S1、第二侧表面S2与上表面S3采集应变片坐标及对应采集到的应变值分别为表1、表2、表3:
表1
表2
表3
床身为灰铸铁,弹性模量E为1.01×1011Pa,泊松比v为0.28。k=6,由公式3得p=6。第一侧表面S1、第二侧表面S2与上表面S3变形曲线为:
其中,求得第一侧表面S1各m/n幂系数
第二侧表面S2各m/n幂系数
上表面S3各m/n幂系数
第四步 直线度评估
采用端点连线法得到变形曲线的平面直线度分别为:因此,沿导轨安装底面法向n与沿导轨L1安装侧面法向t的测量初始直线度分别:ft(1)=112.42μm/m。采用激光干涉仪测量再次测量床身直线度为
第五步 直线度归一化处理
对数控机床床身几何精度变化进行连续监测,数控机床服役3个月后,第二测量数控机床床身初始直线度为ft(2)=41.91μm/m、fn(2)=62.26μm/m。因此,归一化处理下的基础件沿导轨安装底面法向n与沿导轨L1安装侧面法向t的修正直线度:
本发明所述的数控机床床身几何精度在线监测方法,它基于数控机床床身上粘贴的应变片组采集实时应变数据,通过应变场预测并积分实现数控机床床身变形反算,再通过端点连线法求出机床基础件直线度,具有测量精度高、适应性好的优点,实现了测量—加工无干扰同步进行。应变数据通过无线检测***进行数据传输,具有占用空间小,为机床只能维护、故障预紧与自主维修提供了重要的数据支持。

Claims (1)

1.一种数控机床几何精度在线监测方法,其特征在于,该方法中,首先在数控机床床身表面粘贴应变片;通过搭建的无线应变监测***,采集并远程传输各离散测点应变数据;将离散应变数据作为边界条件,求解满足双调和函数的应力方程,预测机床基础件表面应变场分布,通过积分得到机床基础件变形状态;再通过端点连线法求出机床基础件直线度;方法的具体步骤如下:
第一步布置应变测点
在数控机床床身(1)的第一侧表面S1、第二侧表面S2的上、下边缘,上表面S3的导轨内侧处的边缘位置粘贴三轴应变片,各表面应变片数量均为k,
其中,[]为取整符号,即[x]为不超过x的最大整数;l为导轨长度;各表面沿导轨方向s相邻应变片间距为d,
第二步采集应变数据
利用数控机床几何精度无线应变监测***进行应变数据采集;该***是一个分布式的数据采集***,符合“控制中心-分布式采集节点”的构架模式,包括应变片组(2)、无线应变节点(3)、无线网关(4)、数据中心服务器(5)、上位机(6);无线应变节点(3)内置电桥电阻与放大调理电路,无线应变节点(3)采集应变片组(2)的应变值,并将其无线传输至无线网关(4);通过通用无线分组业务GPRS将无线网关(4)实时监测的测点应变数据传送到数据中心服务器(5);上位机(6)利用数据中心服务器(5)中的应变数据进行几何精度计算;
第三步反算机床基础件变形
首先,设定平面Si(i=1,2,3)表示为第i个表面,即第一侧表面S1、第二侧表面S2或上表面S3;在平面Si建立的直角坐标系为xSi-oSi-ySi;坐标轴xS1、xS2、xS3正方向沿导轨方向s,坐标轴yS1、yS2正方向沿导轨安装底面法向n,坐标轴yS3正方向沿导轨L1安装侧面法向t;那么,各应变片的坐标集合为 为等效平面Si上编号w的应变片坐标;平面Si沿坐标轴ySi方向变形曲线为:
其中,p按如下公式取值,
其中,E为基础件材料的杨氏模量;ν为泊松比;为坐标集合中最小与最大的纵坐标;xSi为沿着坐标轴xSi方向自变量; 为m/n幂系数,组成幂系数数组为零的m/n幂系数组成0元素数组为,
非零的m/n幂系数通过联立求解(5-1),(5-2),(5-3),(5-4)计算:
其中,分别为平面Si上编号w应变片采集的沿xSi方向、ySi方向正应变;为等效平面Si上编号w应变片采集的切应变;
第四步评估直线度
采用端点连线法计算第一侧表面S1、第二侧表面S2与上表面S3变形曲线的平面直线度分别为因此,沿导轨安装底面法向n基础件直线度为:
沿导轨L1安装侧面法向t基础件直线度为:
通过监测***第N(N=1,2,…,+∞)次采集各测点应变值并通过上述方法计算得到的沿导轨安装底面法向n与沿导轨L1安装侧面法向t的基础件直线度为fn(N)、ft(N);
第五步直线度归一化处理
采用直线度测量仪器测量基础件沿导轨安装底面法向n、沿导轨L1安装侧面法向t直线度F1 n、F1 t;测量过程中使监测***中所有应变片平衡清零;一个月后,采用直线度测量仪器再次测量直线度同时,通过监测***首次采集测点应变值并通过计算得到基础件的直线度fn(1)、ft(1);因此,归一化处理下通过监测***第N次监测的基础件沿导轨安装底面法向n与沿导轨L1安装侧面法向t的修正直线度Fn、Ft为:
通过以上步骤完成数控机床几何精度在线监测。
CN201710195633.6A 2017-03-29 2017-03-29 一种数控机床几何精度在线监测方法 Active CN106990756B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710195633.6A CN106990756B (zh) 2017-03-29 2017-03-29 一种数控机床几何精度在线监测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710195633.6A CN106990756B (zh) 2017-03-29 2017-03-29 一种数控机床几何精度在线监测方法

Publications (2)

Publication Number Publication Date
CN106990756A CN106990756A (zh) 2017-07-28
CN106990756B true CN106990756B (zh) 2019-03-05

Family

ID=59412452

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710195633.6A Active CN106990756B (zh) 2017-03-29 2017-03-29 一种数控机床几何精度在线监测方法

Country Status (1)

Country Link
CN (1) CN106990756B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108415368B (zh) * 2018-01-23 2019-08-13 大连理工大学 一种数控机床直线运动轴几何精度快速检测方法
CN108919688B (zh) * 2018-08-14 2020-09-22 南京工大数控科技有限公司 一种离散式测控在线执行方法
CN111360581A (zh) * 2020-04-22 2020-07-03 机械科学研究总院海西(福建)分院有限公司 检验液体静压导轨副实际工况下导向精度的装置及方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1209855A2 (en) * 2000-11-27 2002-05-29 Newport Electronics Inc. Remote internet control of instruments
CN1430721A (zh) * 2001-12-14 2003-07-16 香港中文大学 在线监控冲压过程的方法及***
CN105511402A (zh) * 2015-12-01 2016-04-20 西安交通大学 一种重型机床床身变形自动监测与报警装置
DE102015100655A1 (de) * 2015-01-19 2016-07-21 Cord Winkelmann Linearführungseinrichtung für eine Vorschubachse einer Werkzeugmaschine
CN106404537A (zh) * 2016-11-09 2017-02-15 北京工业大学 一种监测重型机床地基基础变形的实验方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1209855A2 (en) * 2000-11-27 2002-05-29 Newport Electronics Inc. Remote internet control of instruments
CN1430721A (zh) * 2001-12-14 2003-07-16 香港中文大学 在线监控冲压过程的方法及***
DE102015100655A1 (de) * 2015-01-19 2016-07-21 Cord Winkelmann Linearführungseinrichtung für eine Vorschubachse einer Werkzeugmaschine
CN105511402A (zh) * 2015-12-01 2016-04-20 西安交通大学 一种重型机床床身变形自动监测与报警装置
CN106404537A (zh) * 2016-11-09 2017-02-15 北京工业大学 一种监测重型机床地基基础变形的实验方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
《地脚螺栓蠕变松弛对大型数控机床几何精度衰退的影响》;刘海波 等;《西安交通大学学报》;20150930;第49卷(第9期);全文
《数控机床智能化状态监测与故障诊断***》;赵敏;《中国优秀硕士学位论文全文数据库 工程科技I辑》;20120415;全文
《数控机床运行可靠性控制技术研究》;唐贤进;《中国优秀硕士学位论文全文数据库 工程科技I辑》;20110115;全文
《螺栓预紧序列与机床精度保持的相关性研究》;吴嘉锟;《中国优秀硕士学位论文全文数据库 工程科技I辑》;20160315;全文

Also Published As

Publication number Publication date
CN106990756A (zh) 2017-07-28

Similar Documents

Publication Publication Date Title
CN105945948B (zh) 一种应用于工业机器人的tcp在线快速标定方法及装置
CN106990756B (zh) 一种数控机床几何精度在线监测方法
CN104374453B (zh) 多通道数字化偏载误差智能修正称重传感器及其称量方法
CN112434359B (zh) 一种高铁桥墩沉降曲线预测方法及***
CN103335589B (zh) 无缝线路钢轨温度应力传感节点的标定方法
CN102353481A (zh) 基于二维正交函数的压力传感器温度和压力互补的方法及装置
CN1318155C (zh) 一种利用数据冗余提高轧机厚度控制精度的方法
CN102323095A (zh) 基于角度监测识别受损索松弛索支座广义位移的方法
CN108919063B (zh) 一种基于电容修正原理的电场遥测***及方法
CN105090084A (zh) 风机在线监测***及其监测方法
CN103091182A (zh) 基于双轴徐变仪的采集测量***
CN102323083A (zh) 基于应变监测识别受损索松弛索支座广义位移的方法
CN111339699A (zh) 基于宏观应力与残余应力拓度的箱型结构重变形量预测方法
CN208818234U (zh) 桥梁监测传感器在线校准用测量***
CN208606701U (zh) 一种热压机加热板平整度测试装置
CN102323077A (zh) 空间坐标监测识别受损索松弛索支座广义位移的方法
CN101545806B (zh) 一种大范围测量物体受力振动数据的修正方法
CN206514834U (zh) 一种新型专用高度测量仪
CN105571479A (zh) 数字疲劳传感器
CN105091821A (zh) 光伏压延玻璃在线自动测厚装置及方法
CN103631249A (zh) 一种检测工业控制***实时性的***及方法
CN112115540B (zh) 一种融合实测数据的海工支撑结构疲劳分析方法及***
CN203083517U (zh) 汽车玻璃电子检具
CN110487575B (zh) 基于支座反力和倾角对称斜率的梁结构损伤识别方法
CN208795186U (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