CN101629801B - 数控磨床导轨热误差确定方法 - Google Patents
数控磨床导轨热误差确定方法 Download PDFInfo
- Publication number
- CN101629801B CN101629801B CN200910056599XA CN200910056599A CN101629801B CN 101629801 B CN101629801 B CN 101629801B CN 200910056599X A CN200910056599X A CN 200910056599XA CN 200910056599 A CN200910056599 A CN 200910056599A CN 101629801 B CN101629801 B CN 101629801B
- Authority
- CN
- China
- Prior art keywords
- guide rail
- control grinder
- error
- thermal
- model
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Automatic Control Of Machine Tools (AREA)
Abstract
本发明涉及一种数控磨床导轨热误差确定方法,具体步骤是:1.测量数控磨床有限位置上导轨热误差;2.建立数控磨床导轨热特性分析响应面模型;3.数控磨床导轨有限元模型的优化及热误差计算。在利用数控机床导轨热误差的测量基础上,以降低热误差数值模拟值的误差为优化目标,以有限元分析的边界条件为设计变量,采用响应面法建立设计变量和目标函数的近似模型。在有限元分析参数的优化过程中,利用有限位置导轨热变形的测量数据,修整有限元计算边界条件,从而得到接近实际测量值的导轨热变形误差值。通过本发明能解决数控磨床导轨热变形的测量问题,并给出了导轨有限元模型的优化方法,从而可以通过有限元计算得到接近测量值的导轨热变形误差。
Description
技术领域
本发明涉及一种机床导轨误差确定方法,尤其是一种磨床导轨热误差的确定方法。
背景技术
在精密机床的切削加工中,热源对加工精度的影响极大,提高工件的加工精度必须对机床的热变形作定量研究,并在加工过程中作合理控制与补偿。导轨部件是机床的主要热源之一,机床导轨与工作台的相对运动会产生大量的摩擦热,造成导轨与工作台的热变形,影响刀具与工件之间的相对位置关系,最终影响工件的加工精度。尤其对于数控主轴磨床,由于工作台承载着很重的磨头部件,而且工作台的行程很大,由摩擦引起的导轨热变形是导向误差的主要原因之一,导轨热变形对于零件的加工精度影响很大。因此,为提高其加工精度,需要获得导轨整个行程范围内的热变形动态变化量,以便进行实时误差补偿。但是,对导轨全行程范围内进行热误差实时测量目前还没有有效的方法。
随着有限元理论的不断完善和数值模拟技术的日趋成熟,使得数值模拟技术在热特性分析领域得到了广泛应用,成为热变形分析的重要手段。有限元数值模拟技术可以定量地计算出温度分布状态和由此产生的热位移、应力和应变等数据。但由于有限元分析边界条件的不确定性(例如热源和对流系数),无法对具体的机床导轨给出精确的计算结果,而只能给出定性的分析结果。
发明内容
本发明是要提供一种数控磨床导轨热误差确定方法,用于解决数控磨床导轨热误差测量与补偿问题。
为解决上述技术问题,本发明采用的技术方案是:一种数控磨床导轨热误差确定方法,包括以下具体步骤:
(1)测量数控磨床有限位置上导轨热误差
在Y方向沿导轨运动Z轴方向安装若干个位移传感器,其中二个位移传感器应放置在两个目标位置,其他位移传感器均匀分布在目标位置之间:
设传感器1、2、……r的读数分别为δ1(i)、δ2(i)、……δr(i),i=0,1,2,……N为测量序号,设测量开始时主轴热变形为0,第i次测量Y方向1、2、……r传感器位置处的热变形量分别为:
δ1i=δ1(i)-δ1(0)δ2i=δ2(i)-δ2(0)……δri=δr(i)-δr(0)
(2)建立数控磨床导轨热特性分析响应面模型
利用有限元分析软件建立数控磨床导轨有限元模型,把有限元模型的边界条件包括导轨表面对流系数(ξ1、ξ2……ξM)和热源的发热量(Q1、Q2……QM)作为设计变量,用X表示,取导轨上r个传感器安装位置的热变形计算误差平方和作为特征量,用F(X)表示,作为设计变量X的优化计算指标;
对r个测点的热变形误差模拟计算值Δ1i、Δ2i、……Δri与实际测量值δ1i、δ2i、……δri进行比较,取其误差平方和构造设计变量X的优化计算指标
采用上述响应面方法建立优化计算指标F(X)和设计变量X之间的近似模型,得到响应面函数表达式:
F(X)=f(ξ1,ξ2,…ξM,Q1,Q2,…QM)
(3)计算数控磨床导轨有限元模型的优化及热误差
导轨有限元分析模型优化的数学模型为:
其中Xmin、Xmax分别为设计变量的取值范围;采用上述优化数学模型并选择优化算法,对导轨有限元模型边界条件进行优化设计,得到一组最优设计变量X=[ξ1,ξ2,…ξM,Q1,Q2,…QM];根据优化后的边界条件,利用有限元分析软件计算导轨在整个行程范围内的热变形量。
本发明具有的有益效果:
通过本发明能解决数控磨床导轨热变形的测量问题,并给出了导轨有限元模型的优化方法,从而可以通过有限元计算得到接近测量值的导轨热变形误差。本发明将数值模拟技术和实际试验相结合,利用试验测试数据修正有限元分析边界条件,从而得到准确的导轨热变形计算结果,使其能应用于机床导轨的热变形实时补偿,提高机床的加工精度。
附图说明
图1是数控磨床沿Z轴移动引起的Y方向热变形测量的传感器安装示意图;其中:左目标测量位置1,中间测量位置2,右目标测量位置3,工件安装台4,传感器安装架5,位移传感器6。
具体实施方式
本发明要解决数控磨床热误差补偿建模中导轨热变形误差测量的技术问题,在有限位置测量导轨热变形的基础上,用数值模拟方法进行数控机床导轨的热变形误差计算的一种方法。
对于数控主轴磨床,轴向变形(Z向)和径向(X方向)对工件加工精度影响不大,所以主要考虑径向(Y方向)热变形误差沿导轨运动方向的分布,见图1。
在利用数控机床导轨热误差的基础上,以降低热误差数值模拟值的误差为优化目标,以有限元分析的边界条件为设计变量,采用响应面法建立设计变量和目标函数的近似模型。在有限元分析参数的优化过程中,利用有限位置导轨热变形的测量数据,修整有限元计算边界条件,从而得到接近实际测量值的导轨热变形误差值。
1.响应面近似模型方法概述
响应面法是以试验设计为基础的用于处理多变量问题建模和分析的一套统计处理技术,其基本思想是通过近似构造一个具有明确表达形式的多项式y=f(x1,x2,...,xk),以显式的响应面近似模型逼近特征量与设计变量间复杂的隐式关系,并在此基础上对原有模型进行修正。响应面方法的数学表达式是多元线性回归方程。在本发明中,该方法被用于根据导轨热变形的实际测量数据,修正有限元模型的边界条件(对流传热系数和热源)。
工程中最广泛采用的响应面近似函数为二阶模型:
其中:n为设计变量个数;xi为设计变量;α0、αi、αii、αij(i<j)为多项式的待定系数。
例如:对于4个设计变量的情况,二阶响应面函数展开为:
令x5=x1 2,x6=x2 2,x7=x3 2,x8=x4 2,x9=x1x2,x10=x1x3,x11=x1x4,x12=x2x3,x13=x2x4,x14=x3x4,β0=α0,β1=α1,β2=α2,β3=α3,β4=α4,β5=α11,β6=α22,β7=α33,β8=α44,β9=α12,β10=α13,β11=α14,β12=α23,β13=α24,β14=α34,则上式可转换成多元线型回归模型:
可得统一的简单形式:
其中k是待定系数βi的个数。为了确定βi,需要做m≥k次的独立试验,通过求解得到相应的系数βi。响应面可用矩阵形式表示为:
其中,y为实际的试验值,为响应面近似函数值,ε为服从正态分布N(0,σ2)的拟合误差。
为了估计二次多项式的系数β,可以用最小二乘法,使误差平方和最小,即:
根据微积分学极值定理,令:
很显然,上式是一个具有k个未知数和k个方程的线性方程组,写成矩阵形式如下所示:
(xβ-y)Tx=0 (8)
解上式可得β,即得到响应面函数表达式。
以上是响应面法的基本原理及响应面系数求解的具体方法。
2为解决上述技术问题,本发明的具体方法是:
(1)数控磨床有限位置上导轨热误差的测量
数控磨床的导轨结构示意图以及热误差测量的传感器安装方式如图1所示。在Y方向沿导轨运动Z轴方向安装r个位移传感器,其中2个位移传感器应放置在两个目标位置,其他位移传感器均匀分布在目标位置之间。目标位置的选择应考虑加工行程,尽可能靠近行程最大点。对于导轨热变形的测量要求定时(例如每5分钟)读取位移传感器的测量值,并根据测量值计算Y方向的变形量误差。考虑工作台的运动问题,传感器可以选用电容传感器或电涡流传感器等非接触位移传感器。
设传感器1、2、……r的读数分别为δ1(i)、δ2(i)、……δr(i),i=0,1,2,……N为测量序号。设测量开始时主轴热变形为0,第i次测量Y方向1、2、……r传感器位置处的热变形量分别为
δ1i=δ1(i)-δ1(0)δ2i=δ2(i)-δ2(0)……δri=δr(i)-δr(0)(9)
(2)数控磨床导轨热特性分析响应面模型的建立
利用有限元分析软件建立数控磨床导轨有限元模型。
在对研究导轨温度场和热变形误差有限元分析的过程中,由于有限元模型的边界条件(对流传热系数和热源)很难理论精确确定,使得有限元计算结果与试验结果之间往往存在明显的误差,为了缩小这种误差,可借助响应面方法进行模型修正。响应面方法的优势在于可以通过较少的试验(有限元计算)获得设计变量和特征量之间足够准确的相互关系,并且可以只用简单代数表达式展现出来。
因此把有限元模型的边界条件包括导轨表面对流系数(ξ1、ξ2......ξM)和热源的发热量(Q1、Q2......QM)作为设计变量,用X表示。在理论计算的基础上结合经验选取设计变量的因子水平范围,即:ξ1min≤ξ1≤ξ1max、ξ2min≤ξ2≤ξ2max、......ξMmin≤ξM≤ξMmax、Q1min≤Q1≤Q1max、Q2min≤Q2≤Q2max......QMmin≤QM≤QMmax。为了使数值模拟结果逼近实验结果,取导轨上r个传感器安装位置的热变形计算误差平方和作为特征量,用F(X)表示,作为设计变量X的优化计算指标。
对r个测点的热变形误差模拟计算值Δ1i、Δ2i、......Δri与实际测量值δ1i、δ2i、......δri进行比较,取其误差平方和构造设计变量X的优化计算指标
采用实验设计方法进行数据采样(实验设计方法有析因实验设计、中心复合设计(CCD)、正交实验设计以及D-最优设计等),利用有限元模型计算出一组设计变量X取值范围内的优化计算指标F(X)计算值,采用上述响应面方法(公式1-8)建立优化计算指标F(X)和设计变量X之间的近似模型。得到响应面函数表达式:
F(X)=f(ξ1,ξ2,…ξM,Q1,Q2,…QM)
(3)数控磨床导轨有限元模型的优化及热误差计算
对于数控机床导轨有限元模型,为了使数值模拟结果逼近实际测量结果,其优化目标应尽可能减少由于有限元模型边界条件难于精确确定所导致的有限元计算误差,所以导轨有限元分析模型优化的数学模型为:
其中Xmin、Xmax分别为设计变量的取值范围。
这样,有限元模型优化问题就转化为对式(11)的求解。
采用上述优化数学模型并选择优化算法,对导轨有限元模型边界条件进行优化设计。得到一组最优设计变量X=[ξ1,ξ2,…ξM,Q1,Q2,…QM]
根据优化后的边界条件,利用有限元分析软件计算导轨在整个行程范围内的热变形量。
Claims (1)
1.一种数控磨床导轨热误差确定方法,其特征在于:包括以下具体步骤:
(1)测量数控磨床有限位置上导轨热误差
沿导轨运动Z轴方向安装若干个位移传感器,其中二个位移传感器放置在两个目标位置,其他位移传感器均匀分布在目标位置之间;
设传感器1、2、……r的读数分别为δ1(i)、δ2(i)、……δr(i),i=0,1,2,……N,为测量序号,设测量开始时主轴热变形为0,第i次测量Y方向传感器1、2、……r位置处的热变形量分别为:
δ1i=δ1(i)-δ1(0)δ2i=δ2(i)-δ2(0)……δri=δr(i)-δr(0)
(2)建立数控磨床导轨热特性分析响应面模型
利用有限元分析软件建立数控磨床导轨有限元模型,把有限元模型的边界条件包括导轨表面对流系数(ξ1、ξ2……ξM)和热源的发热量(Q1、Q2……QM)作为设计变量,用X表示,取导轨上r个传感器安装位置的热变形计算误差平方和作为特征量,用F(X)表示,作为设计变量X的优化计算指标;
对r个测点的热变形误差模拟计算值Δ1i、Δ2i、……Δri与实际测量值δ1i、δ2i、……δri进行比较,取其误差平方和构造设计变量X的优化计算指标:
采用响应面方法建立优化计算指标F(X)和设计变量X之间的近似模型,得到响应面函数表达式:
F(X)=f(ξ1,ξ2,…ξM,Q1,Q2,…QM)
(3)数控磨床导轨有限元模型的优化及热误差计算
导轨有限元分析模型优化的数学模型为:
其中Xmin、Xmax分别为设计变量的取值范围;采用上述优化数学模型并选择优化算法,对导轨有限元模型边界条件进行优化设计,得到一组设计变量X=[ξ1,ξ2,…ξM,Q1,Q2,…QM];根据优化后的边界条件,利用有限元分析软件计算导轨在整个行程范围内的热变形量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910056599XA CN101629801B (zh) | 2009-08-18 | 2009-08-18 | 数控磨床导轨热误差确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910056599XA CN101629801B (zh) | 2009-08-18 | 2009-08-18 | 数控磨床导轨热误差确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101629801A CN101629801A (zh) | 2010-01-20 |
CN101629801B true CN101629801B (zh) | 2012-01-04 |
Family
ID=41575000
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910056599XA Expired - Fee Related CN101629801B (zh) | 2009-08-18 | 2009-08-18 | 数控磨床导轨热误差确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101629801B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101804581A (zh) * | 2010-03-23 | 2010-08-18 | 四川普什宁江机床有限公司 | 一种机床热变形自动补偿的实现方法 |
DE102010003303A1 (de) * | 2010-03-25 | 2011-09-29 | Deckel Maho Seebach Gmbh | Verfahren und Vorrichtung zum Kompensieren einer temperaturabhängigen Lageveränderung an einer Werkzeugmaschine |
CN103273380B (zh) * | 2013-06-18 | 2015-08-26 | 上海理工大学 | 基于仿真模型的数控机床热变形误差在线监测方法 |
CN105022344B (zh) * | 2014-04-22 | 2017-11-28 | 上海开通数控有限公司 | 数控机床热误差补偿方法 |
CN105631112B (zh) * | 2015-12-25 | 2021-08-24 | 鼎奇(天津)主轴科技有限公司 | 一种厚板类零件的热变形仿真分析与建模方法 |
CN105631122B (zh) * | 2015-12-25 | 2021-08-24 | 鼎奇(天津)主轴科技有限公司 | 一种机床大件的热变形仿真分析与建模方法 |
CN108188821B (zh) * | 2018-01-08 | 2019-04-26 | 东北大学 | 一种数控机床滚珠丝杠进给***热误差预测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006053644A1 (de) * | 2004-11-17 | 2006-05-26 | Tecosim Technische Simulation Gmbh | Verfahren zur erzeugung eines berechnungsmodells einer mechanischen struktur |
CN101290266A (zh) * | 2008-06-17 | 2008-10-22 | 上海理工大学 | 数控机床热误差补偿温度测点位置的确定方法 |
-
2009
- 2009-08-18 CN CN200910056599XA patent/CN101629801B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006053644A1 (de) * | 2004-11-17 | 2006-05-26 | Tecosim Technische Simulation Gmbh | Verfahren zur erzeugung eines berechnungsmodells einer mechanischen struktur |
CN101290266A (zh) * | 2008-06-17 | 2008-10-22 | 上海理工大学 | 数控机床热误差补偿温度测点位置的确定方法 |
Non-Patent Citations (5)
Title |
---|
JP特开2006-65716A 2006.03.09 |
胡金祥,蒋书运,黄国庆.高速内圆磨床工件主轴***热特性分析.《精密制造与自动化》.2007,(第1期,总第169期),全文. * |
郭学祥,胡友民,夏军勇,吴波.基于有限元分析的机床导轨热变形研究.《设计与研究》.2007,(第3期),全文. * |
闫占辉,于骏一.自准直导轨磨床床身有限元优化设计.《精密制造与自动化》.2003,(第4期,总第156期),全文. * |
龚伟,汪中厚,黄裕华,汪雪英.基于Pro/MECHANICA的磨床床身有限元热变形分析.《精密制造与自动化》.2007,(第3期,总第171期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN101629801A (zh) | 2010-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101629801B (zh) | 数控磨床导轨热误差确定方法 | |
CN101797704B (zh) | 数控滚齿机热变形误差补偿方法 | |
US6178389B1 (en) | Method of determining the measuring uncertainty of a coordinate measuring device | |
CN105759719B (zh) | 一种基于无偏估计拆分模型的数控机床热误差预测方法及*** | |
CN102176988B (zh) | 金属冲压成形品的裂缝判定方法、装置、程序以及记录介质 | |
CN103823991B (zh) | 一种考虑环境温度的重型机床热误差预测方法 | |
CN101446994A (zh) | 数控机床的热误差最小二乘支持向量机建模方法 | |
CN102179725B (zh) | 数控机床热特征监控测点的布置方法 | |
CN103273380B (zh) | 基于仿真模型的数控机床热变形误差在线监测方法 | |
CN105404237A (zh) | 一种基于空间网格补偿方式的数控机床空间误差建模方法 | |
CN105397560A (zh) | 一种干切数控滚齿机床及工件热变形误差补偿方法 | |
CN102540884A (zh) | 工具机热误差智能调适装置及其方法 | |
CN103885386A (zh) | 一种基于卡尔曼滤波的灰色模型热误差数据处理方法 | |
CN104915468A (zh) | 一种垂直载荷下滚柱直线导轨副的力变形误差计算方法 | |
CN111650881A (zh) | 数控机床定位误差预测与补偿方法、***及存储介质 | |
CN102658503A (zh) | 一种基于内置传感器的数控机床进给***的模态测试方法 | |
He et al. | Identification and adjustment of guide rail geometric errors based on BP neural network | |
Valentini et al. | Recent advancements in the hole-drilling strain-gage method for determining residual stresses | |
CN101498572B (zh) | 凸轮二维曲线测量方法 | |
Landkammer et al. | Investigations on residual stress generation in full-forward-extrusion | |
CN116611314A (zh) | 一种基于物理信息神经网络的机床加工过程热误差在线评估方法 | |
CN114895627B (zh) | 一种分区域机床热误差补偿方法 | |
Blaser | Adaptive learning control for thermal error compensation | |
Zhuang et al. | Autonomous calibration of hexapod machine tools | |
CN101908084B (zh) | 1000mw超临界汽轮机转子轮槽面温度场重构的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120104 Termination date: 20140818 |
|
EXPY | Termination of patent right or utility model |