CN109376377A - 高速加工机床整机结构热力学建模与热设计方法 - Google Patents
高速加工机床整机结构热力学建模与热设计方法 Download PDFInfo
- Publication number
- CN109376377A CN109376377A CN201811041397.3A CN201811041397A CN109376377A CN 109376377 A CN109376377 A CN 109376377A CN 201811041397 A CN201811041397 A CN 201811041397A CN 109376377 A CN109376377 A CN 109376377A
- Authority
- CN
- China
- Prior art keywords
- machine tool
- speed processing
- processing machine
- thermal
- complete machine
- 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
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automatic Control Of Machine Tools (AREA)
- Numerical Control (AREA)
Abstract
本发明提供了一种高速加工机床整机结构热力学建模与热设计方法,其包括以下步骤:步骤1:高速加工机床三维数字化建模;步骤2:高速加工机床主要热源发热功率和相关换热系数计算计算;步骤3:机床平面结合部热阻参数计算;步骤4:高速加工机床整机结构热力学建模与热特性计算;步骤5:高速加工机床整机结构热态设计方法。采用本发明提供的高速加工机床整机结构热力学设计方法,能够大幅提高高速加工机床整机结构热力学建模精度,缩短设计周期。不仅便于高速加工机床的正向设计,而且提高一次设计成功率。
Description
技术领域
本发明涉及一种高速加工机床整机结构热力学建模与热设计方法,属于数控机床设计领域。
背景技术
数控机床正朝着高速与高精度方向发展,产业界对高速加工机床的热态性能要求越来越高。建立高速加工机床整机热力学模型,开展整机热力学分析,继而完成机床的热力学设计,是高速机床研发的首要技术环节。
建立高精度的高速加工机床热力学模型,是开展高速加工机床热力学设计的理论基础。近年来,国内外学者围绕着机床热力学建模与热特性计算,开展了大量的研究,概括起来,其主要研究工作包括:
①电主轴热力学:考虑主轴中滚动轴承和内置电机转子、定子作为主要热源,计算主轴各表面的换热系数,建立电主轴的热力学模型,通过有限元法,计算出电主轴的温度场和热位移场。
②机床结合部热阻:建立机床主要结合面热阻,以机床床身为研究对象,建立机床局部热力学模型,通过有限元法,计算考虑热阻的机床床身温度场和热位移场。
应该看出,以上两种机床热力学单一因素建模与设计方法,均不能满足高速加工机床的设计要求。从机床的正向热设计角度出发,应该建立考虑主要热源(高速电主轴、工件主轴、切削液二次热源)与结合部(平面结合部、滚动结合部)接触热阻的高速加工机床的整机热力学模型,基于该模型,可以开展高速加工机床热力学特性分析,继而提出高速加工机床整机热设计,提升机床的热态性能。
发明内容
技术问题:针对高速加工机床整机热力学建模与热设计中存在的问题,本发明专利提出了高速加工机床主要热源(高速电主轴、工件主轴、切削液二次热源)与典型结合部(平面结合部、滚动结合部)接触热阻的建模方法,继而建立了高速加工机床的整机热力学模型;基于高速加工机床整机热力学灵敏度分析结果,提出了高速加工机床整机热力学设计方法。本发明专利不仅便于高速加工机床正向设计,而且提高机床热态设计准确性和一次设计成功率。
技术方案:本发明的一种高速加工机床整机结构热力学建模与热设计方法包括以下步骤:
步骤1:高速加工机床三维数字化建模:采用三维数字化软件,将拟设计的高速加工机床的初步结构转化为三维CAD模型;
步骤2:高速加工机床主要热源发热功率和相关换热系数计算:采用Palmgren经验公式计算滚动轴承发热功率,采用传热学理论计算相关部位对流换热系数;为高速加工机床整机热态建模提供基础数据;
步骤3:机床平面结合部热阻参数计算:采用分形几何理论,计算机床平面结合部接触热阻;
步骤4:高速加工机床整机结构热力学建模与热特性计算:建立含主要热源与典型结合部接触热阻的高速加工机床的整机热力学模型;开展机床关键零部件、水冷系统物理参数对高速加工机床整机热力学灵敏度分析与计算,寻找影响机床热力学性能的薄弱环节;
步骤5:高速加工机床整机结构热态设计方法:针对步骤4寻找到的影响机床热力学性能的薄弱环节,以提高机床热态性能为目标,提出高速加工机床整机结构动态设计方法。
所述步骤1:高速加工机床三维数字化建模,采用三维数字化软件Solidworks或者Pro-E,将拟设计的高速加工机床的初步结构转化为三维CAD模型。
所述所述步骤2:高速加工机床主要热源发热功率和相关换热系数计算,分为以下三步:
步骤2a:高速加工机床电主轴与工件主轴发热功率计算
1)滚动轴承发热功率Q1计算,
Q1=1.047×10-4(M0+M1)·n (1)
M1=f1p1dm (3)
式中:M0为润滑剂的流体动力损耗,M1为弹性滞后和局部差动滑动的摩擦损耗,n为轴承转速,r/min;dm为轴承中径,mm;f0为与轴承类型和润滑方式有关的经验常数,v为在工作温度下润滑剂的运动粘度,cst,f1为与轴承类型和所受负荷有关的参数;p1为确定轴承摩擦力矩的计算负荷N;
2)电主轴电机转子、定子发热功率计算
P0为额定功率,电机功率系数α,效率η,则热损耗Q=P0α(1-η),其中热损耗2/3分布在电机转子,热损耗1/3分布在电机定子;
3)工件主轴带传动发热功率qbelt计算
qbelt=Pbeltin(1-η) (3)
式中:Pbeltin带传动输入功率W,η为带传动效率,带传动的热损耗一半传递到皮带上,
一半传递到带轮上,且主动轮和从动轮上的摩擦生热按两带轮的皮带包角等比例分配;
步骤2b:高速加工机床电主轴与工件主轴对流换热系数计算
换热系数α和努赛尔数Nu的关系:
式中:λ为流体导热系数,Nuf为努赛尔数,D为几何特征的定型尺寸;
1)电主轴内置电机定子冷却水套强制对流换热系数h1计算
式中:λ为流体导热系数,Nuf为努赛尔数;Prf为普朗特数,Ref为雷诺数,L为冷却水流道长度;
2)电主轴电机定子内表面与转子外表面的传热系数h2计算
式中:r1为转子外径(m);σ为定、转子间气隙(m);
3)工件主轴旋转表面与空气间的对流换热系数h3计算
4)工件主轴带轮对流换热系数h4计算
h4=9.7+5.33×U0.8 (8)
式中:h4—对流换热系数W/(m2·K);U—带轮表面线速度m/s;
5)主轴外表面与空气间的对流换热系数,取复合传热系数为h5=9.7W/(m2·K)
步骤2c)切削液“二次热源”的建模方法
切削液流道部分的对流换热系数Nu计算:
切削液溅到的表面的对流换热系数:将切削液流道对流系数乘以0.6,得到该表面的对流系数;式(9)的定性温度取流体和壁面的平均温度Tm,定型尺寸取流体流过的壁面长度L。
所述步骤3:高速加工机床主要结合部热阻参数计算,分为以下2步:
步骤3a)固定结合面接触热阻Rc的计算
总的接触热阻TCR为
式中:Rci为第i个微凸体接触热阻。
无量纲总接触热阻为
式中:Rci为第i个微凸体接触热阻,为相对面积接触比,为第i个微凸体的截面积,N(a′s)为接触点数目;
步骤3b)轴承外圈与轴承座间的接触热阻计算
轴承接触热导为,
式中:Π为接触热导(W/K);hring、hgap分别为轴承外圈、外圈与轴承座平均间隙厚度m;λring、λgap分别为轴承外圈和间隙中空气的导热系数W/(m2·K),A为轴承套圈的圆柱形外表面面积;
平均间隙hgap可由下式计算
hgap=hgap0-(Tring-Thousing)a′·rhousing (13)
式中:hgap0为轴承外圈与轴承座之间的原始间隙,m;a′为线性膨胀系数,对于钢材可取其值为11.7×10-6/K。
所述步骤4:高速加工机床整机结构热力学建模与热特性计算,建立含主要热源与典型结合部接触热阻的高速加工机床的整机热力学模型;开展高速加工机床主要热源的发热功率、机床关键零部件、水冷系统物理参数对高速加工机床整机热力学性能的灵敏度分析计算,寻找影响机床热性能的薄弱环节与主要因素。
所述主要热源为高速电主轴、工件主轴、切削液二次热源;典型结合部为平面结合部、滚动结合部。
所述高速加工机床整机热力学性能包括温度场与热位移场。
所述步骤5:高速加工机床整机结构热态设计方法,围绕步骤4发现的影响机床热性能的薄弱环节与主要因素,以降低加工过程中的热变形、提高加工精度为设计目标,修改机床主要零部件的结构参数,重新设计相关物理参数,完成高速加工机床热设计。
有益效果:本发明提供的高速加工机床整机结构热力学建模与热设计方法,是基于完善的高速加工机床热力学模型与系统的热力学分析计算结果而提出的。因此采用本发明“高速加工机床整机结构热力学建模与热设计方法”,可以提高机床热态设计准确性和一次设计成功率。
附图说明
图1本发明专利的设计步骤与内容,
图2高速加工机床三维数字化模型,
a为内圆复合磨床3D-CAD模型,
b为整机有限元分析模型,
图3无量纲接触热阻R*在不同参数影响下计算结果,
a为无量纲接触热阻与接触载荷在不同G*下的关系(E/σ=106,D=1.3),
b为无量纲接触热阻与接触载荷在不同D下的关系(E/σ=106,G=10-7),
c为无量纲接触热阻与接触载荷在不同E/σ下的关系(D=1.3,G=10-7),
d为无量纲接触面积与接触热阻在不同接触变形情形下的关系,
图4整机额定工况下的温度场分布云图,
a为等距测视角的整机额定工况温度场分布云图,
b为背视图视角的整机额定工况温度场分布云图,
图5整机额定工况下的热位移场分布云图,
a为等距测视角的整机额定工况热位移场分布云图,
b为背视图视角的整机额定工况热位移场分布云图,
图6热设计后的整机温度场分布云图,
a为俯视图视角的热设计后的整机温度场分布云图,
b为背视图视角的热设计后的整机温度场分布云图。
具体实施方式
下面结合一个实施例(精密数控内圆复合磨床),对本发明的高速加工机床整机结构热力学设计方法作进一步详细说明。
图1给出了本发明高速加工机床整机结构热力学设计方法的具体内容,包括如下步骤:步骤1:高速加工机床三维数字化建模
采用三维建模软件(如Solidworks)建立精密数控内圆复合磨床的三维数字化(CAD、CAE)模型,如图2所示。
步骤2:高速加工机床主要热源发热功率计算
步骤2a:采用公式(1)-(3),计算高速加工机床磨削电主轴滚动轴承发热功率Q,在不同温度情况下,计算结果如表1所示
表1内圆磨削电主轴轴承生热率q与温度的对应关系(q=Q/V)
采用公式(1)-(3),计算高速加工机床工件主轴滚动轴承发热功率Q,在不同温度情况下,得到不同温度下轴承生热率,见表2所示。
表2轴承生热率与温度的对应关系(q=Q/V)
电主轴内置电机的额定功率为16kW,正常工作状态下功率系数为0.9,功率损耗为1.5kW。假设该损耗全部转化成热量,其中2/3(1kW)由定子发出,1/3(0.5kW)由转子发出,则电机定、转子生热率分别为254250W/m3、255056W/m3。
带传动的效率为96%,带传动输入功率Pbeltin为可倾瓦动压轴承和推力滑动轴承等的摩擦损耗功率,为5.5kW,工件主轴上从动带轮的摩擦热损(qbelt)为43.86W。
步骤2b:高速加工机床电主轴与工件主轴对流换热系数计算
采用公式(5)计算强制冷却水对流换热系数h1=13719W/(m2·K)
采用公式(6)电机内置转子、定子对流换热系数h2=116.2W/(m2·K)
采用公式(7)工件主轴旋转表面与空气间的对流换热系数计算h3=100W/(m2·K)
根据公式(8),皮带轮直径为245mm、200mm和144mm的径向表面所对应的对流换热系数分别为81W/(m2·K)、70W/(m2·K)和56W/(m2·K);砂轮主轴直径为95mm和56mm的径向表面所对应的对流换热系数分别为43W/(m2·K)和31W/(m2·K)
主轴外表面与空气间的对流换热系数,取复合传热系数为h5=9.7W/(m2·K)
步骤2c:采用公式(9)计算高速加工机床切削液作为“二次热源”相关对流换热系数计算,额定工况下切削液流量100L/min。主要表面的对流换热换热系数见表3所示:
表3床身主要表面的对流换热系数
步骤3:机床平面结合部热阻参数计算
步骤3a:采用公式(10)(11)计算总的接触热阻R和无量纲接触热阻R*,计算结果如图3所示。
步骤3b:采用公式(12)计算轴承外圈与轴承座间的接触热阻,对于普通主轴轴承(取4级加工精度),取hring=5×10-6m,λring=24W/(m·K),λgap=2.84×10-2W/(m·K),则主轴轴承与配合面的接触系数见表4所示。
表4主轴轴承与配合面间的接触系数h(单位:W/m2·K)
步骤4:高速加工机床整机结构热力学建模与热特性计算
本实施例,内圆磨削孔径为150mm,砂轮直径为80mm。将内圆磨削电主轴、床头箱和切削液等热源,全部施加于机床热力学分析模型,计算得到整机在额定工况下的温度场分布云图,见图4。
整机在额定工况下的热位移场分布云图,见图5,图中黑色网格表示未发生热变形时的轮廓。
额定工况下,整机主要部位的热位移,列于表5。
表5整机在额定工况下各构件平均热位移
经换算,整机热变形引起的内孔和锥面磨削误差为43.30μm。
步骤5:高速加工机床整机结构热态设计方法
针对步骤4分析结果,本实施例对机床结构进行改进设计,将与流道尺寸相同的平板覆盖在流道上,使切削液不直接与床身接触,减少传递给床身热量,进而减少床身热变形,其主要部位的热位移列于表6。
表6热设计之后各构件的平均热位移
经热设计之后,整机热变形引起的内孔和锥面磨削误差为1.15μm。
由设计结果看,在磨削孔径为150mm,内圆磨削砂轮直径为80mm状态下,对磨床进行热设计前、后的加工热误差对比,见表7。可以看出,采用本发明专利“高速加工机床整机结构热力学建模与热设计方法”,可以大幅度提高机床的磨削精度。
表7热设计前后加工精度对比
Claims (8)
1.一种高速加工机床整机结构热力学建模与热设计方法,其特征在于该方法包括以下步骤:
步骤1:高速加工机床三维数字化建模:采用三维数字化软件,将拟设计的高速加工机床的初步结构转化为三维CAD模型;
步骤2:高速加工机床主要热源发热功率和相关换热系数计算:采用Palmgren经验公式计算滚动轴承发热功率,采用传热学理论计算相关部位对流换热系数;为高速加工机床整机热态建模提供基础数据;
步骤3:机床平面结合部热阻参数计算:采用分形几何理论,计算机床平面结合部接触热阻;
步骤4:高速加工机床整机结构热力学建模与热特性计算:建立含主要热源与典型结合部接触热阻的高速加工机床的整机热力学模型;开展机床关键零部件、水冷系统物理参数对高速加工机床整机热力学灵敏度分析与计算,寻找影响机床热力学性能的薄弱环节;
步骤5:高速加工机床整机结构热态设计方法:针对步骤4寻找到的影响机床热力学性能的薄弱环节,以提高机床热态性能为目标,提出高速加工机床整机结构动态设计方法。
2.根据权利要求1所述的高速加工机床整机结构热力学建模与热设计方法,其特征在于,所述步骤1:高速加工机床三维数字化建模,采用三维数字化软件Solidworks或者Pro-E,将拟设计的高速加工机床的初步结构转化为三维CAD模型。
3.根据权利要求1所述的高速加工机床整机结构热力学建模与热设计方法,其特征在于,所述所述步骤2:高速加工机床主要热源发热功率和相关换热系数计算,分为以下三步:
步骤2a:高速加工机床电主轴与工件主轴发热功率计算
1)滚动轴承发热功率Q1计算,
Q1=1.047×10-4(M0+M1)·n (1)
M1=f1p1dm (3)
式中:M0为润滑剂的流体动力损耗,M1为弹性滞后和局部差动滑动的摩擦损耗,n为轴承转速,r/min;dm为轴承中径,mm;f0为与轴承类型和润滑方式有关的经验常数,v为在工作温度下润滑剂的运动粘度,cst,f1为与轴承类型和所受负荷有关的参数;p1为确定轴承摩擦力矩的计算负荷N;
2)电主轴电机转子、定子发热功率计算
P0为额定功率,电机功率系数α,效率η,则热损耗Q=P0α(1-η),其中热损耗2/3分布在电机转子,热损耗1/3分布在电机定子;
3)工件主轴带传动发热功率qbelt计算
qbelt=Pbeltin(1-η) (3)
式中:Pbeltin带传动输入功率W,η为带传动效率,带传动的热损耗一半传递到皮带上,一半传递到带轮上,且主动轮和从动轮上的摩擦生热按两带轮的皮带包角等比例分配;
步骤2b:高速加工机床电主轴与工件主轴对流换热系数计算
换热系数α和努赛尔数Nu的关系:
式中:λ为流体导热系数,Nuf为努赛尔数,D为几何特征的定型尺寸;
1)电主轴内置电机定子冷却水套强制对流换热系数h1计算
式中:λ为流体导热系数,Nuf为努赛尔数;Prf为普朗特数,Ref为雷诺数,L为冷却水流道长度;
2)电主轴电机定子内表面与转子外表面的传热系数h2计算
式中:r1为转子外径(m);σ为定、转子间气隙(m);
3)工件主轴旋转表面与空气间的对流换热系数h3计算
4)工件主轴带轮对流换热系数h4计算
h4=9.7+5.33×U0.8 (8)
式中:h4-对流换热系数W/(m2·K);U-带轮表面线速度m/s;
5)主轴外表面与空气间的对流换热系数,取复合传热系数为h5=9.7W/(m2·K)
步骤2c)切削液“二次热源”的建模方法
切削液流道部分的对流换热系数Nu计算:
切削液溅到的表面的对流换热系数:将切削液流道对流系数乘以0.6,得到该表面的对流系数;式(9)的定性温度取流体和壁面的平均温度Tm,定型尺寸取流体流过的壁面长度L。
4.根据权利要求1所述的高速加工机床整机结构热力学建模与热设计方法,其特征在于,所述步骤3:高速加工机床主要结合部热阻参数计算,分为以下2步:
步骤3a)固定结合面接触热阻Rc的计算
总的接触热阻TCR为
式中:Rci为第i个微凸体接触热阻。
无量纲总接触热阻为
式中:Rci为第i个微凸体接触热阻,为相对面积接触比,为第i个微凸体的截面积,N(a′s)为接触点数目;
步骤3b)轴承外圈与轴承座间的接触热阻计算
轴承接触热导为,
式中:Π为接触热导(W/K);hring、hgap分别为轴承外圈、外圈与轴承座平均间隙厚度m;λring、λgap分别为轴承外圈和间隙中空气的导热系数W/(m2·K),A为轴承套圈的圆柱形外表面面积;
平均间隙hgap可由下式计算
hgap=hgap0-(Tring-Thousing)a′·rhousing (13)
式中:hqap0为轴承外圈与轴承座之间的原始间隙,m;a′为线性膨胀系数,对于钢材可取其值为11.7×10-6/K。
5.根据权利要求1所述的高速加工机床整机结构热力学建模与热设计方法,其特征在于,所述步骤4:高速加工机床整机结构热力学建模与热特性计算,建立含主要热源与典型结合部接触热阻的高速加工机床的整机热力学模型;开展高速加工机床主要热源的发热功率、机床关键零部件、水冷系统物理参数对高速加工机床整机热力学性能的灵敏度分析计算,寻找影响机床热性能的薄弱环节与主要因素。
6.根据权利要求5所述的高速加工机床整机结构热力学建模与热设计方法,其特征在于,所述主要热源为高速电主轴、工件主轴、切削液二次热源;典型结合部为平面结合部、滚动结合部。
7.根据权利要求5所述的高速加工机床整机结构热力学建模与热设计方法,其特征在于,所述高速加工机床整机热力学性能包括温度场与热位移场。
8.根据权利要求1所述的高速加工机床整机结构热力学建模与热设计方法,其特征在于,所述步骤5:高速加工机床整机结构热态设计方法,围绕步骤4发现的影响机床热性能的薄弱环节与主要因素,以降低加工过程中的热变形、提高加工精度为设计目标,修改机床主要零部件的结构参数,重新设计相关物理参数,完成高速加工机床热设计。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811041397.3A CN109376377B (zh) | 2018-09-07 | 2018-09-07 | 高速加工机床整机结构热力学建模与热设计方法 |
PCT/CN2018/104779 WO2020047880A1 (zh) | 2018-09-07 | 2018-09-10 | 高速加工机床整机结构热力学建模与热设计方法 |
PCT/CN2019/078377 WO2020048110A1 (zh) | 2018-09-07 | 2019-03-15 | 高速加工机床整机结构热力学建模与热设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811041397.3A CN109376377B (zh) | 2018-09-07 | 2018-09-07 | 高速加工机床整机结构热力学建模与热设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109376377A true CN109376377A (zh) | 2019-02-22 |
CN109376377B CN109376377B (zh) | 2019-10-11 |
Family
ID=65404726
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811041397.3A Active CN109376377B (zh) | 2018-09-07 | 2018-09-07 | 高速加工机床整机结构热力学建模与热设计方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN109376377B (zh) |
WO (2) | WO2020047880A1 (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110076974A (zh) * | 2019-04-18 | 2019-08-02 | 甘肃机电职业技术学院 | 基于增材制造技术的注塑模具随形冷却水道的设计方法 |
CN110705147A (zh) * | 2019-09-18 | 2020-01-17 | 北京工业大学 | 一种数控机床主轴热态特性综合性理论建模与分析方法 |
WO2020048110A1 (zh) * | 2018-09-07 | 2020-03-12 | 东南大学 | 高速加工机床整机结构热力学建模与热设计方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112364444A (zh) * | 2020-09-23 | 2021-02-12 | 江苏赛洋机电科技有限公司 | 一种基于有限元模型数值仿真的数控机床温度测点优化方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102009370A (zh) * | 2010-11-03 | 2011-04-13 | 北京航空航天大学 | 一种高速切削数控机床的反馈补偿方法 |
CN102452020A (zh) * | 2010-10-22 | 2012-05-16 | 西安交通大学 | 一种数控机床刀具温度场和热变形定量分析方法 |
CN102867088A (zh) * | 2012-09-12 | 2013-01-09 | 西安交通大学 | 一种应用于电主轴稳态温度场的热网络建模方法 |
CN103034169A (zh) * | 2012-12-26 | 2013-04-10 | 西安理工大学 | 一种数控机床热误差建模与补偿的方法 |
US20140081607A1 (en) * | 2009-09-15 | 2014-03-20 | Invensys Systems, Inc. | Thermodynamic Phase Equilibrium Analysis Based On a Reduced Composition Domain |
CN104200019A (zh) * | 2014-08-27 | 2014-12-10 | 电子科技大学 | 一种重型数控落地铣镗床热特性数值仿真方法 |
CN105022344A (zh) * | 2014-04-22 | 2015-11-04 | 上海开通数控有限公司 | 数控机床热误差补偿方法 |
CN105631122A (zh) * | 2015-12-25 | 2016-06-01 | 鼎奇(天津)主轴科技有限公司 | 一种机床大件的热变形仿真分析与建模方法 |
CN106407528A (zh) * | 2016-09-05 | 2017-02-15 | 西安交通大学 | 一种基于动力学模型的机床主轴轴承产热计算方法 |
CN107498391A (zh) * | 2017-08-03 | 2017-12-22 | 南京航空航天大学 | 机床切削加工过程的热力载荷模拟和误差检测系统及方法 |
CN108241787A (zh) * | 2018-01-12 | 2018-07-03 | 哈尔滨理工大学 | 极端工况下静压回转工作台热特性研究方法 |
CN108304657A (zh) * | 2018-02-02 | 2018-07-20 | 重庆大学 | 基于有限元的机床关重件残余应力多工序连续建模仿真方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3405965B2 (ja) * | 2000-07-06 | 2003-05-12 | ファナック株式会社 | 工作機械の熱変位補正方法 |
CN102063549B (zh) * | 2011-01-07 | 2012-07-04 | 西安交通大学 | 一种机床主轴的流固热耦合数值仿真方法 |
CN103366050B (zh) * | 2013-06-25 | 2016-01-13 | 华中科技大学 | 一种高速电主轴力-热耦合建模方法 |
CN109376377B (zh) * | 2018-09-07 | 2019-10-11 | 东南大学 | 高速加工机床整机结构热力学建模与热设计方法 |
-
2018
- 2018-09-07 CN CN201811041397.3A patent/CN109376377B/zh active Active
- 2018-09-10 WO PCT/CN2018/104779 patent/WO2020047880A1/zh active Application Filing
-
2019
- 2019-03-15 WO PCT/CN2019/078377 patent/WO2020048110A1/zh active Application Filing
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140081607A1 (en) * | 2009-09-15 | 2014-03-20 | Invensys Systems, Inc. | Thermodynamic Phase Equilibrium Analysis Based On a Reduced Composition Domain |
CN102452020A (zh) * | 2010-10-22 | 2012-05-16 | 西安交通大学 | 一种数控机床刀具温度场和热变形定量分析方法 |
CN102009370A (zh) * | 2010-11-03 | 2011-04-13 | 北京航空航天大学 | 一种高速切削数控机床的反馈补偿方法 |
CN102867088A (zh) * | 2012-09-12 | 2013-01-09 | 西安交通大学 | 一种应用于电主轴稳态温度场的热网络建模方法 |
CN103034169A (zh) * | 2012-12-26 | 2013-04-10 | 西安理工大学 | 一种数控机床热误差建模与补偿的方法 |
CN105022344A (zh) * | 2014-04-22 | 2015-11-04 | 上海开通数控有限公司 | 数控机床热误差补偿方法 |
CN104200019A (zh) * | 2014-08-27 | 2014-12-10 | 电子科技大学 | 一种重型数控落地铣镗床热特性数值仿真方法 |
CN105631122A (zh) * | 2015-12-25 | 2016-06-01 | 鼎奇(天津)主轴科技有限公司 | 一种机床大件的热变形仿真分析与建模方法 |
CN106407528A (zh) * | 2016-09-05 | 2017-02-15 | 西安交通大学 | 一种基于动力学模型的机床主轴轴承产热计算方法 |
CN107498391A (zh) * | 2017-08-03 | 2017-12-22 | 南京航空航天大学 | 机床切削加工过程的热力载荷模拟和误差检测系统及方法 |
CN108241787A (zh) * | 2018-01-12 | 2018-07-03 | 哈尔滨理工大学 | 极端工况下静压回转工作台热特性研究方法 |
CN108304657A (zh) * | 2018-02-02 | 2018-07-20 | 重庆大学 | 基于有限元的机床关重件残余应力多工序连续建模仿真方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020048110A1 (zh) * | 2018-09-07 | 2020-03-12 | 东南大学 | 高速加工机床整机结构热力学建模与热设计方法 |
CN110076974A (zh) * | 2019-04-18 | 2019-08-02 | 甘肃机电职业技术学院 | 基于增材制造技术的注塑模具随形冷却水道的设计方法 |
CN110705147A (zh) * | 2019-09-18 | 2020-01-17 | 北京工业大学 | 一种数控机床主轴热态特性综合性理论建模与分析方法 |
CN110705147B (zh) * | 2019-09-18 | 2023-08-04 | 北京工业大学 | 一种数控机床主轴热态特性综合性理论建模与分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109376377B (zh) | 2019-10-11 |
WO2020047880A1 (zh) | 2020-03-12 |
WO2020048110A1 (zh) | 2020-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109376377B (zh) | 高速加工机床整机结构热力学建模与热设计方法 | |
Min et al. | An improved thermal model for machine tool bearings | |
Ma et al. | Simulation and experimental study on the thermally induced deformations of high-speed spindle system | |
Ma et al. | Thermal characteristics analysis and experimental study on the high-speed spindle system | |
Uhlmann et al. | Thermal modelling of a high speed motor spindle | |
CN103366050B (zh) | 一种高速电主轴力-热耦合建模方法 | |
Meng et al. | Research on thermal resistance network modeling of motorized spindle based on the influence of various fractal parameters | |
Liu et al. | Improved thermal resistance network model of motorized spindle system considering temperature variation of cooling system | |
Tan et al. | An optimal convective heat transfer coefficient calculation method in thermal analysis of spindle system | |
Li et al. | Thermal-mechanical coupling calculation method for deformation error of motorized spindle of machine tool | |
Oezkaya et al. | A new reverse engineering method to combine FEM and CFD simulation three-dimensional insight into the chipping zone during the drilling of Inconel 718 with internal cooling | |
Chen et al. | Simulation on thermal characteristics of high-speed motorized spindle | |
Li et al. | Calculation method of convective heat transfer coefficients for thermal simulation of a spindle system based on RBF neural network | |
Fang et al. | An improved thermal performance modeling for high-speed spindle of machine tool based on thermal contact resistance analysis | |
Li et al. | Method for modifying convective heat transfer coefficients used in the thermal simulation of a feed drive system based on the response surface methodology | |
Zhao et al. | Thermal analysis and experimental study on the spindle of the high-speed machining center | |
Lyu et al. | Texture optimization and verification for the thrust bearing used in rotary compressors based on a transient tribo-dynamics model | |
Fan et al. | Thermoelectric-based cooling system for high-speed motorized spindle I: design and control mechanism | |
Fan et al. | An optimization method for thermal behavior of high-speed spindle of gear form grinding machine | |
Deng et al. | A predictive model for temperature rise of spindle–bearing integrated system | |
Gao et al. | Thermal behavior study and experimental testing of deep groove ball bearings | |
Deng et al. | A model for investigating the temperature of trochoidal machining | |
Kamenev | Method for automated building of spindle thermal model with use of CAE system | |
Chen et al. | A thermal analyzing model for high speed motorized spindles—thermal properties | |
Zhang et al. | Thermal performance analysis for the machine tool's spindle |
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 |