CN105631122B - 一种机床大件的热变形仿真分析与建模方法 - Google Patents

一种机床大件的热变形仿真分析与建模方法 Download PDF

Info

Publication number
CN105631122B
CN105631122B CN201511000318.0A CN201511000318A CN105631122B CN 105631122 B CN105631122 B CN 105631122B CN 201511000318 A CN201511000318 A CN 201511000318A CN 105631122 B CN105631122 B CN 105631122B
Authority
CN
China
Prior art keywords
thermal deformation
machine tool
thermal
analysis
large part
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
Application number
CN201511000318.0A
Other languages
English (en)
Other versions
CN105631122A (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.)
Dingqi Tianjin Spindle Co ltd
Original Assignee
Dingqi Tianjin Spindle Co ltd
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 Dingqi Tianjin Spindle Co ltd filed Critical Dingqi Tianjin Spindle Co ltd
Priority to CN201511000318.0A priority Critical patent/CN105631122B/zh
Publication of CN105631122A publication Critical patent/CN105631122A/zh
Application granted granted Critical
Publication of CN105631122B publication Critical patent/CN105631122B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Numerical Control (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种适用于机床大件的热变形仿真分析与建模的方法，包括对机床大件的热变形仿真分析与热变形建模方法；所述机床大件的热变形仿真分析方法是基于通用有限元分析软件Ansys进行的，并在发明中详细阐述了运用“间接法”进行热‑结构耦合分析的具体步骤；所述机床大件的热变形建模方法，包括在有限元分析完成后，根据提取出的需要建立热变形模型位置处的热变形结果以及按照一定规则在机床大件预先选定位置处所提取的特定方向的热变形结果，利用多项式拟合及多元线性回归的方法，建立机床大件上确定位置处的热变形的数学模型。

Description

一种机床大件的热变形仿真分析与建模方法
技术领域
本发明涉及数控机床技术领域，特别是一种机床大件的热变形仿真分析与建模方法。
背景技术
随着数控机床整机及零部件设计、制造、装配和材料等相关技术的不断进步，几何误差、刀具磨损、伺服误差等在数控机床整体误差中所占的比例逐渐减少。在高速高精度极端加工条件下，热变形日益成为影响机床加工精度的重要因素。大量研究与加工实践表明，对于高速高精度机床，由热变形引起的加工制造误差所占的比例为40％～70％。因此，对于机床热变形建模技术的研究就显得十分必要。目前，常采用的热变形建模技术，一般是根据机床的物理结构，分析各部位的温度变化情况，以便寻找机床的热关键点，对机床的测温点进行优化布置，建立温度关键点与机床特定部位热变形之间的关系，从而得到其热误差模型。这种热变形建模方法，依赖于准确的选择测温关键点来建立热变形模型，一旦测温关键点的选择不合理，建立的热变形模型将会偏差很大。
发明内容
本发明为解决公知技术中存在的技术问题而提供一种机床大件的热变形仿真分析与建模方法。
本发明为解决公知技术中存在的技术问题所采取的技术方案是：一种适合于机床大件的热变形仿真分析与建模方法，包括对机床大件的热变形仿真分析与热变形建模方法。
所述机床大件的热变形仿真分析是基于通用有限元分析软件Ansys进行的，包括如下步骤：
步骤一，采用间接法进行热-结构耦合分析，选择适用于间接法进行热-结构耦合分析的单元；所述的间接法是指，Ansys提供的一种热-结构耦合分析方法，这种方法要求首先进行热分析，然后将求得的节点温度作为体载荷施加在结构应力分析中，本热变形仿真分析方法是基于间接法建立的；
步骤二，对机床几何模型进行简化，去掉机床大件的热变形有限元分析影响较小的结构，如小的倒角、圆角、小孔等细小结构；
步骤三，将机床几何模型导入Ansys中，添加材料属性；
步骤四，热载荷的施加，包括热源的施加、环境温度设定、机床表面传热系数的设定以及接触位置的热阻定义；
步骤五，按照上述步骤定义完成后首先进行热分析求解；
步骤六，热分析完成后，保存热分析结果文件；进行单元的转换，将热单元转换成对应的结构单元；删除热分析中定义的接触热阻；定义约束，以及设置接触面处的刚度、阻尼；载荷的施加，将热分析结果文件作为载荷施加到结构分析中；
步骤七，进行机床结构应力分析，得到有限元分析结果；
步骤八，提取有限元分析结果，包括在有限元分析完成后，确定机床大件上需要建立热变形数学模型的位置并提取热变形结果，以及按照一定规则在机床大件上预先选定位置并提取其特定方向的热变形结果。
所述步骤八，按照一定规则在机床大件上预先选定位置，包括预先选定的位置所占的范围要尽可能大，这样热变形结果的变化就会更明显；所述特定方向，最好与需要建立热变形模型的位置处的方向垂直。
所述机床大件的热变形建模方法，包括利用多项式拟合及多元线性回归的方法，建立机床大件确定位置热变形的数学模型。其中，首先利用多项式拟合的方法，根据提取出的机床大件上需要建立热变形模型位置处的热变形结果，建立相同阶次的不同加载情况下的热变形数学模型；其次，利用多元线性回归的方法，建立按照一定规则在机床大件预先选定位置处所提取的特定方向的热变形结果与第一次所建立的数学模型系数之间的函数关系。
附图说明
图1为机床大件热变形的仿真分析与建模方法的流程示意图；
图2为一台卧式加工中心的CAD模型示意图；
图中：1、预先选定位置1；2、预先选定位置2；3、立柱上导轨；4、预先选定位置3；5、预先选定位置4；
图3为图2所示的卧式加工中心的16种加载工况；
图中：溜板位置-左，主轴箱位置-上，工作台位置-前，与其各自中间位置的距离为0.2m；
图4为图3中两种热载荷所对应的加载情况；
图中：a)热载荷1所对应的加载情况；b)热载荷2所对应的加载情况；
图5为图2所示的卧式加工中心的有限元网格模型示意图；
图中：a)施加热载荷后的网格模型；b)定义约束后的网格模型。
图6为图2所示的卧式加工中心02工况下的温度场及变形云图示意图；
图中：a)温度场云图；b)变形云图。
图7为图2所示的卧式加工中心立柱上导轨位置处的Z向热变形数学模型的效果验证示意图。
具体实施方式
为能进一步了解本发明的发明内容、特点及功效，兹例举以下实施例，并配合附图详细说明如下：
请参阅图1-7，一种机床大件的热变形仿真分析与建模方法，包括确定机床大件上需要建立热变形数学模型的位置及方向；对机床大件进行有限元仿真分析；利用多项式拟合及多元线性回归的方法，建立关于机床大件确定位置及方向的热变形的数学模型。
参照图1所示流程，对图2所示的某一卧式加工中心的立柱上导轨的热变形进行仿真分析与建模。
首先，确定对立柱上导轨的Z向热变形建立数学模型，如图2中“黄线”所示；
其次，按照图3所示的16中工况进行有限元仿真分析，有限元仿真分析的具体步骤如下：
步骤一，本次分析选用间接法进行热-结构耦合分析，选用Solid87单元；
步骤二，对机床几何模型进行简化，去掉机床大件的热变形有限元分析影响较小的结构，如小的倒角、圆角、小孔等细小结构；
步骤三，将简化好的机床几何模型导入Ansys中添加材料属性，机床大件包括床身、立柱、溜板、转台以及主轴箱按照QT500-7的属性添加，包括弹性模量1730GPa，泊松比0.3，密度7300kg/m3，热膨胀系数9.1×10-61/K，热导率47W/(m·K)，X、Y及Z三向导轨按照45钢的材料属性添加，包括弹性模量2190GPa，泊松比0.3，密度7850kg/m3，热膨胀系数1.17×10-51/K，热导率48W/(m·K)；
步骤四，热载荷的施加，包括热源的施加、环境温度设定、机床表面传热系数的设定以及接触位置的热阻定义；其中，热源的施加位置包括电机座位置、轴承座位置、导轨位置、滑块位置、主轴位置，大小按照图4所示进行施加；环境温度设定为20℃；机床表面传热系数10W/(m2·K)；接触热阻位置的定义，包括滑块-大件接触面、导轨-滑块接触面，导轨-机床大件接触面；
步骤五，按照上述方式定义完成后首先进行热分析求解；
步骤六，热分析完成后，保存热分析结果文件；进行单元的转换，将热单元转换成Solid92结构单元；删除热分析中定义的接触热阻；定义约束，包括对机床底面的X、Y与Z方向的自由度的限制，以及接触面处的刚度、阻尼的设置；载荷的施加，将热分析结果文件作为载荷施加到结构分析中；
步骤七，进行机床结构应力分析，得到有限元分析结果，图6中列出了该加工中心在02工况下的有限元分析结果。
步骤八，有限元分析完成后，提取有限元分析结果文件。包括机床上导轨的Z向热变形，以及预先选定位置1、2、3及4的Y向热变形。
再者，根据不同加载工况下提取的立柱上导轨的Z向热变形结果，利用多项式拟合的方法建立相同阶次的立柱上导轨的Z向热变形数学模型，
δzi＝pi·x (1)
式中pi＝(pi1 pi2 pi3)，i＝1,2…,16，x＝(x2 x 1)T，δzi为i工况下立柱上导轨的Z向热变形，pi为表达式的系数，x为坐标值；再利用多元线性回归的方法，建立预先选定位置处所提取的热变形结果与式(1)系数之间的函数关系，
pi＝A·δyi+B (2)
其中，δyi＝(δyi1 δyi2 δyi3 δyi4)T，i＝1,2…,16，δyi为提取出的i工况下预先选定位置处的1、2、3及4的Y向热变形结果；联立式(1)与式(2)，最终可以得到立柱上导轨的Z向热变形数学模型，
δz＝(A·δy+B)·x (3)
式中，δz为立柱上导轨的Z向热变形，δy为提取出的预先选定位置处的1、2、3及4的Y向热变形结果，x＝(x2 x 1)T，x为坐标值，A是3×4的系数矩阵，
Figure BDA0000890024970000041
B是3×1的常数向量，
Figure BDA0000890024970000042
最后，将工况6、工况9提取出的预先选定位置1、2、3及4的Y向热变形代入得到的热变形数学模型(3)中，得到两种情况下立柱上导轨的Z向变形情况，并与直接从有限元分析软件中提取出Z向变形情况进行比对，对比结果如图7所示。
尽管上面结合附图对本发明的优选实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，并不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可以做出很多形式，这些均属于本发明的保护范围之内。

Claims (2)

1.一种机床大件的热变形仿真分析与建模方法，其特征在于，包括如下步骤：
步骤一，采用间接法进行热-结构耦合分析，选择适用于间接法进行热-结构耦合分析的单元；所述的间接法是指，Ansys提供的一种热-结构耦合分析方法，这种方法要求首先进行热分析，然后将求得的节点温度作为体载荷施加在结构应力分析中，本热变形仿真分析方法是基于间接法建立的；
步骤二，对机床几何模型进行简化，去掉倒角、圆角、小孔；
步骤三，将机床几何模型导入Ansys中，添加材料属性；
步骤四，热载荷的施加，包括热源的施加、环境温度设定、机床表面传热系数的设定以及接触位置的热阻定义；
步骤五，按照上述步骤定义完成后首先进行热分析求解；
步骤六，热分析完成后，保存热分析结果文件；进行单元的转换，将热单元转换成对应的结构单元；删除热分析中定义的接触热阻；定义约束，以及设置接触面处的刚度、阻尼；载荷的施加，将热分析结果文件作为载荷施加到结构分析中；
步骤七，进行机床结构应力分析，得到有限元分析结果；
步骤八，提取有限元分析结果，包括在有限元分析完成后，确定机床大件上需要建立热变形数学模型的位置并提取热变形结果，以及按照机床大件上预先选定位置提取其特定方向的热变形结果；
所述机床大件的热变形建模方法，包括利用多项式拟合及多元线性回归的方法，建立关于机床大件确定位置的热变形的数学模型，其中，首先利用多项式拟合的方法，根据提取出的机床大件上需要建立热变形模型位置处的热变形结果，建立相同阶次的不同加载情况下的热变形数学模型；其次，利用多元线性回归的方法，建立机床大件预先选定位置处所提取的特定方向的热变形结果与第一次所建立的数学模型系数之间的函数关系，具体方式为：有限元分析完成后，提取有限元分析结果文件，包括机床上导轨的Z向热变形，以及预先选定位置1、2、3及4的Y向热变形，再者，根据不同加载工况下提取的立柱上导轨的Z向热变形结果，利用多项式拟合的方法建立相同阶次的立柱上导轨的Z向热变形数学模型，
δzi＝pi·x，式中pi＝(pi1 pi2 pi3)，i＝1,2…,16，x＝(x2x1)T，δzi 为i工况下立柱上导轨的Z向热变形，pi为表达式的系数，x为坐标值；再利用多元线性回归的方法，建立预先选定位置处所提取的热变形结果与δzi＝pi·x系数之间的函数关系，
pi＝A·δyi+B，其中，δyi＝(δyi1 δyi2 δyi3 δyi4)T，i＝1,2…,16，δyi为提取出的i工况下预先选定位置处的1、2、3及4的Y向热变形结果；联立δzi＝pi·x和pi＝A·δyi+B得到立柱上导轨的Z向热变形数学模型，
δz＝(A·δy+B)·x
式中，δz为立柱上导轨的Z向热变形，δy为提取出的预先选定位置处的1、2、3及4的Y向热变形结果，x＝(x2x1)T，x为坐标值，A是3×4的系数矩阵，B是3×1的常数向量。
2.根据权利要求1所述的一种机床大件的热变形仿真分析与建模方法，其特征在于，所述步骤八，所述特定方向与建立热变形模型的位置处的方向垂直。
CN201511000318.0A 2015-12-25 2015-12-25 一种机床大件的热变形仿真分析与建模方法 Expired - Fee Related CN105631122B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201511000318.0A CN105631122B (zh) 2015-12-25 2015-12-25 一种机床大件的热变形仿真分析与建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201511000318.0A CN105631122B (zh) 2015-12-25 2015-12-25 一种机床大件的热变形仿真分析与建模方法

Publications (2)

Publication Number Publication Date
CN105631122A CN105631122A (zh) 2016-06-01
CN105631122B true CN105631122B (zh) 2021-08-24

Family

ID=56046051

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201511000318.0A Expired - Fee Related CN105631122B (zh) 2015-12-25 2015-12-25 一种机床大件的热变形仿真分析与建模方法

Country Status (1)

Country Link
CN (1) CN105631122B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107784142B (zh) * 2016-08-31 2021-10-19 株洲中车时代半导体有限公司 半导体功率组件的热仿真方法
CN107498391B (zh) * 2017-08-03 2019-07-02 南京航空航天大学 机床切削加工过程的热力载荷模拟和误差检测系统及方法
CN109376377B (zh) * 2018-09-07 2019-10-11 东南大学 高速加工机床整机结构热力学建模与热设计方法
CN110232206B (zh) * 2019-04-28 2023-07-25 广东工业大学 基于ansys的脉冲功率器件的热变形仿真分析方法
CN112364444A (zh) * 2020-09-23 2021-02-12 江苏赛洋机电科技有限公司 一种基于有限元模型数值仿真的数控机床温度测点优化方法
CN113051686B (zh) * 2021-04-01 2023-03-10 重庆大学 倾斜工作条件下的主轴系统热-结构耦合特性模型创建方法及热误差模型建模方法
CN115688522B (zh) * 2022-11-03 2023-11-03 四川大学 一种悬伸式主轴机床结构优化方法及系统

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06288864A (ja) * 1993-04-06 1994-10-18 Hitachi Ltd 光及び構造の複合解析システム
JP2009036669A (ja) * 2007-08-02 2009-02-19 Toshiba Corp 溶接残留応力解析方法および溶接残留応力解析システム
CN101629801B (zh) * 2009-08-18 2012-01-04 上海理工大学 数控磨床导轨热误差确定方法
CN102063090B (zh) * 2010-12-30 2012-05-16 大连理工大学 数控机床滑枕悬伸变形实时补偿方法
CN104155491A (zh) * 2013-05-13 2014-11-19 木本军生 探针
CN103473421B (zh) * 2013-09-17 2016-06-29 西安电子科技大学 地基天线温度补偿方法
CN103984287A (zh) * 2014-03-12 2014-08-13 江苏齐航数控机床有限责任公司 一种数控机床热误差补偿灰色神经网络建模方法
CN104502625A (zh) * 2014-11-24 2015-04-08 北京交通大学 一种三维浑水脉动流速仪

Also Published As

Publication number Publication date
CN105631122A (zh) 2016-06-01

Similar Documents

Publication Publication Date Title
CN105631122B (zh) 一种机床大件的热变形仿真分析与建模方法
Soori et al. Dimensional and geometrical errors of three-axis CNC milling machines in a virtual machining system
Majda Modeling of geometric errors of linear guideway and their influence on joint kinematic error in machine tools
Ibaraki et al. Machining tests to identify kinematic errors on five-axis machine tools
CN107553220A (zh) 一种数控机床综合误差实时补偿方法
Peng et al. Total differential methods based universal post processing algorithm considering geometric error for multi-axis NC machine tool
CN109765848B (zh) 五轴机床运动轴几何误差贡献值计算及影响评估方法
CN108445839B (zh) 一种基于误差增量的机床精度敏感度分析方法
CN105243218A (zh) 机床热误差精度转化和模型建立方法
CN107081638A (zh) 一种机床工作空间热误差辨识模型建立方法
CN110955979B (zh) 一种考虑几何误差偏相关性的机床加工精度可靠性灵敏度分析方法
Li et al. A volumetric positioning error compensation method for five-axis machine tools
Tian et al. Error Modeling and Sensitivity Analysis of a Five‐Axis Machine Tool
CN104915468A (zh) 一种垂直载荷下滚柱直线导轨副的力变形误差计算方法
Afkhamifar et al. Variational analysis for CNC milling process
Putilova et al. Loading of the manufacturing systems elements in the process of unsteady mode cutting and the models of their arrangement deviations
CN104200063B (zh) 机床空间加工误差的非确定性描述及预测方法
Liu et al. Thermal–mechanical coupling analysis and experimental study on CNC machine tool feed mechanism
Khodaygan Manufacturing error compensation based on cutting tool location correction in machining processes
Martin et al. Virtual manufacturing: prediction of work piece geometric quality by considering machine and set-up accuracy
CN116909209A (zh) 一种考虑动态热误差的数控机床误差建模和预测方法
CN110532667B (zh) 一种用于提高精密机床空间几何误差模型建模精度的方法
Terrier et al. VERNE-a five-axis parallel kinematics milling machine
Barari et al. NURBS representation of estimated surfaces resulting from machining errors
CN110110381B (zh) 一种卧式加工中心工艺信息的热误差感知方法

Legal Events

Date Code Title Description
C06 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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210824

Termination date: 20211225

CF01 Termination of patent right due to non-payment of annual fee