WO2020048110A1 - Thermodynamic modeling and thermal design method for overall structure of high-speed processing machine - Google Patents

Thermodynamic modeling and thermal design method for overall structure of high-speed processing machine Download PDF

Info

Publication number
WO2020048110A1
WO2020048110A1 PCT/CN2019/078377 CN2019078377W WO2020048110A1 WO 2020048110 A1 WO2020048110 A1 WO 2020048110A1 CN 2019078377 W CN2019078377 W CN 2019078377W WO 2020048110 A1 WO2020048110 A1 WO 2020048110A1
Authority
WO
WIPO (PCT)
Prior art keywords
machine tool
processing machine
speed processing
thermal
calculation
Prior art date
Application number
PCT/CN2019/078377
Other languages
French (fr)
Chinese (zh)
Inventor
蒋书运
朱国振
Original Assignee
东南大学
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 东南大学 filed Critical 东南大学
Publication of WO2020048110A1 publication Critical patent/WO2020048110A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Definitions

  • the invention relates to a method for thermodynamic modeling and thermal design of a high-speed processing machine tool structure, and belongs to the field of numerical control machine tool design.
  • thermodynamic model of the high-speed processing machine tool conducting the thermodynamic analysis of the whole machine, and then completing the thermodynamic design of the machine tool are the first technical links in the research and development of high-speed machine tools.
  • thermodynamic model The establishment of a high-precision high-speed machining machine thermodynamic model is the theoretical basis for developing the high-speed machining machine thermodynamic design.
  • scholars at home and abroad have carried out a lot of research around machine tool thermodynamic modeling and calculation of thermal characteristics.
  • the main research work includes:
  • thermodynamics Considering the rolling bearings in the spindle as well as the built-in motor rotor and stator as the main heat source, the heat transfer coefficients of the surfaces of the spindle are calculated, and the thermodynamic model of the electro-spindle is established. The temperature field and thermal displacement of the electro-spindle are calculated by the finite element method field.
  • Thermal resistance at the joint of the machine tool Establish the main joint surface thermal resistance of the machine tool, take the machine tool bed as the research object, establish a local thermodynamic model of the machine tool, and calculate the temperature field and thermal displacement field of the machine tool bed considering the thermal resistance through the finite element method.
  • thermodynamic single-factor modeling and design methods of machine tools cannot meet the design requirements of high-speed processing machine tools.
  • the integration of high-speed processing machine tools that consider the thermal resistance of the main heat sources (high-speed electric spindle, workpiece spindle, secondary heat source of cutting fluid) and the joints (plane joints, rolling joints) should be established.
  • the thermodynamic characteristics analysis of high-speed processing machine tools can be carried out, and then the thermal design of the whole machine of high-speed processing machine tools is proposed to improve the thermal performance of the machine tools.
  • the invention patent proposes the main heat source of the high-speed processing machine tool (high-speed electric spindle, workpiece spindle, secondary heat source of cutting fluid) and the typical joint (plane (Joint, rolling joint) contact thermal resistance modeling method, and then established the whole machine thermodynamic model of high-speed processing machine tools; Based on the analysis results of the whole machine thermodynamic sensitivity of high-speed processing machine tools, a thermodynamic design method of high-speed processing machine tools was proposed.
  • the invention patent not only facilitates the forward design of high-speed processing machine tools, but also improves the thermal design accuracy of the machine tools and the design success rate at one time.
  • thermodynamic modeling and thermal design method of the present invention includes the following steps:
  • Step 1 Three-dimensional digital modeling of high-speed processing machine tools: The three-dimensional digital software is used to convert the preliminary structure of the proposed high-speed processing machine tools into three-dimensional CAD models;
  • Step 2 Calculation of the heating power and related heat transfer coefficient of the main heat source of high-speed processing machine tools: Calculate the heating power of rolling bearings using Palmgren's empirical formula, and calculate the convective heat transfer coefficient of relevant parts by using heat transfer theory; Basic data;
  • Step 3 Calculation of thermal resistance parameters of the plane joint of the machine tool: using fractal geometry theory, contact thermal resistance of the plane joint of the computer bed;
  • Step 4 Thermodynamic modeling and thermal characteristics calculation of the high-speed processing machine tool structure: Establish a complete machine thermodynamic model of the high-speed processing machine tool including the contact thermal resistance between the main heat source and the typical joint; develop the key parameters of the machine tool, and the physical parameters of the water-cooling system against Analysis and calculation of the thermodynamic sensitivity of the entire machine tool to find the weak links that affect the thermodynamic performance of the machine tool;
  • Step 5 The thermal design method of the high-speed processing machine tool structure: Aiming at the weak links affecting the thermodynamic performance of the machine tool found in step 4, aiming at improving the thermal performance of the machine tool, a dynamic design method of the high-speed processing machine tool machine structure is proposed.
  • Step 1 Three-dimensional digital modeling of the high-speed processing machine tool.
  • the three-dimensional digital software Solidworks or Pro-E is used to convert the preliminary structure of the proposed high-speed processing machine tool into a three-dimensional CAD model.
  • the step 2 the calculation of the heating power of the main heat source and the related heat transfer coefficient of the high-speed processing machine tool is divided into the following three steps:
  • Step 2a Calculation of heating power of electric spindle and workpiece spindle of high-speed processing machine
  • M 0 is the hydrodynamic loss of the lubricant
  • M 1 is the elastic hysteresis and the friction loss of the local differential sliding
  • n is the speed of the bearing, r / min
  • d m is the bearing diameter, mm
  • f 0 is the bearing Type and empirical constants related to the lubrication method
  • v is the kinematic viscosity of the lubricant at the working temperature
  • cst, f 1 are parameters related to the bearing type and the load
  • p 1 is the calculated load N to determine the bearing friction torque
  • Step 2b Calculation of convective heat transfer coefficient between high-speed machining machine electric spindle and workpiece spindle
  • is the thermal conductivity of the fluid
  • Nu f is the Nusser number
  • D is the fixed dimension of the geometric feature
  • is the thermal conductivity of the fluid
  • Nu f is the Nusser number
  • Pr f is the Prandtl number
  • Re f is the Reynolds number
  • L is the length of the cooling water channel
  • r 1 is the outer diameter of the rotor (m); ⁇ is fixed, the air gap between the rotors (m);
  • Convection heat transfer coefficient of the surface to which the cutting fluid splashes Multiply the convection coefficient of the cutting fluid channel by 0.6 to obtain the convection coefficient of the surface.
  • the qualitative temperature of formula (9) is the average temperature T m of the fluid and the wall. Length L of wall surface through which fluid flows.
  • Step 3 The calculation of the thermal resistance parameters of the main joint of the high-speed processing machine tool is divided into the following two steps:
  • the total contact thermal resistance TCR is NCR
  • R ci is the i-th micro-convex contact thermal resistance.
  • the dimensionless total contact thermal resistance is the dimensionless total contact thermal resistance
  • R ci is the i-th micro-convex contact thermal resistance
  • N (a ′ s ) is the number of contact points
  • Step 3b) Calculation of the contact thermal resistance between the bearing outer ring and the bearing seat
  • the bearing contact thermal conductivity is,
  • is the contact thermal conductivity (W / K); h ring and h gap are the average thickness of the bearing outer ring, the outer ring and the bearing seat, respectively; ⁇ ring and ⁇ gap are the air in the bearing outer ring and the gap, respectively.
  • Thermal conductivity W / (m 2 ⁇ K), A is the cylindrical outer surface area of the bearing ring;
  • the average gap h gap can be calculated by the following formula
  • h gap h gap0- (T ring -T housing ) a′r housing (13)
  • h gap is the original gap between the bearing outer ring and the bearing seat, m;
  • a ′ is the linear expansion coefficient, which can be 11.7 ⁇ 10 -6 / K for steel.
  • Step 4 Thermodynamic modeling and thermal characteristic calculation of the high-speed processing machine tool structure, establishing a complete machine thermodynamic model of the high-speed processing machine tool including the thermal resistance of the main heat source and the typical joint contact thermal resistance; developing the heating power of the main heat source of the high-speed processing machine tool, Sensitivity analysis and calculation of key machine parts and water cooling system physical parameters to the thermodynamic performance of high-speed processing machine tools, looking for weak links and main factors that affect the thermal performance of machine tools.
  • the main heat source is a high-speed electric main shaft, a workpiece main shaft, and a secondary heat source for cutting fluid; typical joints are a planar joint and a rolling joint.
  • thermodynamic performance of the high-speed processing machine tool includes a temperature field and a thermal displacement field.
  • Step 5 The thermal state design method of the high-speed processing machine tool structure, surrounding the weak links and main factors affecting the thermal performance of the machine tool found in step 4, with the goal of reducing thermal deformation and improving processing accuracy during the process of modifying the machine tool Structural parameters of main components, redesign related physical parameters, and complete thermal design of high-speed processing machine tools.
  • thermodynamic modeling and thermal design method of the high-speed processing machine tool structure provided by the present invention is based on the thermodynamic analysis and calculation results of a complete high-speed processing machine tool thermodynamic model and system. Therefore, by adopting the "thermodynamic modeling and thermal design method for the whole structure of a high-speed processing machine tool" of the present invention, the thermal design accuracy of the machine tool and the design success rate at one time can be improved.
  • Figure 2 Three-dimensional digital model of high-speed processing machine
  • a is a 3D-CAD model of the internal compound grinding machine
  • d is the relationship between the dimensionless contact area and the contact thermal resistance under different contact deformation situations
  • Figure 4 Cloud diagram of temperature field distribution under rated conditions of the whole machine
  • a is the cloud diagram of the temperature range of the entire machine under rated working conditions at an isometric perspective
  • b is the distribution cloud diagram of the rated operating temperature field of the whole machine from the rear view perspective
  • a is the cloud map of the thermal displacement field of the entire machine under the rated operating conditions at an isometric perspective
  • Figure 6 Cloud diagram of the whole machine temperature field after thermal design.
  • a is the cloud map of the whole machine temperature field after thermal design from the perspective of the top view
  • b is the cloud map of the whole machine's temperature field after thermal design from the perspective of the back view.
  • thermodynamic design method of the entire structure of the high-speed processing machine tool of the present invention is used to further describe the thermodynamic design method of the entire structure of the high-speed processing machine tool of the present invention.
  • FIG. 1 shows the specific content of the thermodynamic design method of the high-speed processing machine tool structure of the present invention, including the following steps:
  • Step 1 3D digital modeling of high-speed processing machine
  • 3D modeling software (such as Solidworks) is used to establish a 3D digital (CAD, CAE) model of a precision CNC internal cylindrical composite grinder, as shown in Figure 2.
  • Step 2 Calculation of heating power of main heat source of high-speed processing machine
  • Step 2a Use formulas (1)-(3) to calculate the heating power Q of the high-speed machining machine grinding electric spindle rolling bearing. The calculation results are shown in Table 1 at different temperatures.
  • the electric motor's built-in motor has a rated power of 16kW, a power factor of 0.9 under normal working conditions, and a power loss of 1.5kW. Assuming that all the losses are converted into heat, of which 2/3 (1kW) is emitted by the stator and 1/3 (0.5kW) is emitted by the rotor, the motor heating and rotor heat generation rates are 254250W / m 3 and 255056W / m 3 respectively .
  • the efficiency of the belt drive is 96%.
  • the belt drive input power P beltin is the friction loss power of the tilting pad dynamic pressure bearing and thrust sliding bearing, which is 5.5kW.
  • the frictional heat loss of the driven pulley on the workpiece spindle (q belt ) It is 43.86W.
  • Step 2b Calculation of convective heat transfer coefficient between high-speed machining machine electric spindle and workpiece spindle
  • the convection heat transfer coefficients corresponding to the radial surfaces of the pulley diameters of 245mm, 200mm, and 144mm are 81W / (m 2 ⁇ K), 70W / (m 2 ⁇ K), and 56W / (m 2 ⁇ K);
  • the convective heat transfer coefficients of the radial surfaces of the grinding wheel with a diameter of 95mm and 56mm are 43W / (m 2 ⁇ K) and 31W / (m 2 ⁇ K)
  • Step 2c Use formula (9) to calculate the relevant convective heat transfer coefficient of the cutting fluid of the high-speed processing machine tool as a "secondary heat source".
  • the flow rate of the cutting fluid is 100 L / min under rated conditions.
  • the convection heat transfer coefficients of the main surfaces are shown in Table 3:
  • Step 3 Calculation of thermal resistance parameters of the plane joint of the machine tool
  • Step 3a Use formulas (10) and (11) to calculate the total contact thermal resistance R and the dimensionless contact thermal resistance R * . The calculation results are shown in FIG.
  • Step 3b Use formula (12) to calculate the contact thermal resistance between the outer ring of the bearing and the bearing block.
  • h ring 5 ⁇ 10 -6 m
  • ⁇ ring 24W / ( m ⁇ K)
  • ⁇ gap 2.84 ⁇ 10 -2 W / (m ⁇ K)
  • Step 4 Thermodynamic Modeling and Thermal Characteristic Calculation of High Speed Machining Machine Structure
  • the internal grinding hole diameter is 150 mm, and the grinding wheel diameter is 80 mm. All internal heat sources such as electric spindle, bedside box and cutting fluid are applied to the thermodynamic analysis model of the machine tool, and the temperature field distribution cloud diagram of the whole machine under rated working conditions is calculated, as shown in Figure 4.
  • the cloud diagram of the thermal displacement field of the whole machine under rated conditions is shown in Figure 5.
  • the black grid in the figure indicates the outline when no thermal deformation has occurred.
  • the internal hole and tapered surface grinding error caused by the thermal deformation of the whole machine is 43.30 ⁇ m.
  • Step 5 Thermal design method of high-speed processing machine tool structure
  • this embodiment improves the design of the machine tool structure.
  • the flat plate with the same size as the flow path is covered on the flow path, so that the cutting fluid does not directly contact the bed, and the heat transmitted to the bed is reduced, thereby reducing the bed.
  • Thermal deformation, the thermal displacement of the main parts are listed in Table 6.
  • the grinding error of the inner hole and the tapered surface caused by the thermal deformation of the whole machine is 1.15 ⁇ m.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Automatic Control Of Machine Tools (AREA)
  • Numerical Control (AREA)

Abstract

Provided is a thermodynamic modeling and thermal design method for the overall structure of a high-speed processing machine, comprising the following steps: step 1: performing three-dimensional digital modeling of a high-speed processing machine; step 2: calculating a heating power of a main heat source of the high-speed processing machine and a related heat transfer coefficient of the high-speed processing machine; step 3: calculating a thermal resistance parameter of a surface joint portion of the high-speed processing machine; step 4: performing thermodynamic modeling of the overall structure of the high-speed processing machine and calculating a thermal property thereof; and step 5: generating a thermal design method for the overall structure of the high-speed processing machine. The method can be adopted to greatly improve the accuracy of thermodynamic modeling of the overall structure of a high-speed processing machine and shortens a design cycle. The present invention not only facilitates forward design of a high-speed processing machine, but also increases the success rate of a first-time design.

Description

高速加工机床整机结构热力学建模与热设计方法Thermodynamic Modeling and Thermal Design Method for High Speed Machining Machine Structure 技术领域Technical field
本发明涉及一种高速加工机床整机结构热力学建模与热设计方法,属于数控机床设计领域。The invention relates to a method for thermodynamic modeling and thermal design of a high-speed processing machine tool structure, and belongs to the field of numerical control machine tool design.
背景技术Background technique
数控机床正朝着高速与高精度方向发展,产业界对高速加工机床的热态性能要求越来越高。建立高速加工机床整机热力学模型,开展整机热力学分析,继而完成机床的热力学设计,是高速机床研发的首要技术环节。CNC machine tools are developing in the direction of high speed and high precision, and the industry has higher and higher requirements for the thermal performance of high-speed processing machine tools. Establishing the thermodynamic model of the high-speed processing machine tool, conducting the thermodynamic analysis of the whole machine, and then completing the thermodynamic design of the machine tool are the first technical links in the research and development of high-speed machine tools.
建立高精度的高速加工机床热力学模型,是开展高速加工机床热力学设计的理论基础。近年来,国内外学者围绕着机床热力学建模与热特性计算,开展了大量的研究,概括起来,其主要研究工作包括:The establishment of a high-precision high-speed machining machine thermodynamic model is the theoretical basis for developing the high-speed machining machine thermodynamic design. In recent years, scholars at home and abroad have carried out a lot of research around machine tool thermodynamic modeling and calculation of thermal characteristics. In summary, the main research work includes:
①电主轴热力学:考虑主轴中滚动轴承和内置电机转子、定子作为主要热源,计算主轴各表面的换热系数,建立电主轴的热力学模型,通过有限元法,计算出电主轴的温度场和热位移场。①Electro-spindle thermodynamics: Considering the rolling bearings in the spindle as well as the built-in motor rotor and stator as the main heat source, the heat transfer coefficients of the surfaces of the spindle are calculated, and the thermodynamic model of the electro-spindle is established. The temperature field and thermal displacement of the electro-spindle are calculated by the finite element method field.
②机床结合部热阻:建立机床主要结合面热阻,以机床床身为研究对象,建立机床局部热力学模型,通过有限元法,计算考虑热阻的机床床身温度场和热位移场。② Thermal resistance at the joint of the machine tool: Establish the main joint surface thermal resistance of the machine tool, take the machine tool bed as the research object, establish a local thermodynamic model of the machine tool, and calculate the temperature field and thermal displacement field of the machine tool bed considering the thermal resistance through the finite element method.
应该看出,以上两种机床热力学单一因素建模与设计方法,均不能满足高速加工机床的设计要求。从机床的正向热设计角度出发,应该建立考虑主要热源(高速电主轴、工件主轴、切削液二次热源)与结合部(平面结合部、滚动结合部)接触热阻的高速加工机床的整机热力学模型,基于该模型,可以开展高速加工机床热力学特性分析,继而提出高速加工机床整机热设计,提升机床的热态性能。It should be noted that the above two types of thermodynamic single-factor modeling and design methods of machine tools cannot meet the design requirements of high-speed processing machine tools. From the perspective of the forward thermal design of the machine tool, the integration of high-speed processing machine tools that consider the thermal resistance of the main heat sources (high-speed electric spindle, workpiece spindle, secondary heat source of cutting fluid) and the joints (plane joints, rolling joints) should be established. Based on this model, the thermodynamic characteristics analysis of high-speed processing machine tools can be carried out, and then the thermal design of the whole machine of high-speed processing machine tools is proposed to improve the thermal performance of the machine tools.
发明内容Summary of the Invention
技术问题:针对高速加工机床整机热力学建模与热设计中存在的问题,本发明专利提出了高速加工机床主要热源(高速电主轴、工件主轴、切削液二次热源)与典型结合部(平面结合部、滚动结合部)接触热阻的建模方法,继而建立了高速加工机床的整机热力学模型;基于高速加工机床整机热力学灵敏度分析结果,提出了高速加工机床整机热力学设计方法。本发明专利不仅便于高速加工机床正向设计,而且提高机床热态设计准确性和一次设计成功率。Technical problem: In view of the problems existing in the thermodynamic modeling and thermal design of the high-speed processing machine tool, the invention patent proposes the main heat source of the high-speed processing machine tool (high-speed electric spindle, workpiece spindle, secondary heat source of cutting fluid) and the typical joint (plane (Joint, rolling joint) contact thermal resistance modeling method, and then established the whole machine thermodynamic model of high-speed processing machine tools; Based on the analysis results of the whole machine thermodynamic sensitivity of high-speed processing machine tools, a thermodynamic design method of high-speed processing machine tools was proposed. The invention patent not only facilitates the forward design of high-speed processing machine tools, but also improves the thermal design accuracy of the machine tools and the design success rate at one time.
技术方案:本发明的一种高速加工机床整机结构热力学建模与热设计方法包括以下步骤:Technical solution: A high-speed processing machine tool structure thermodynamic modeling and thermal design method of the present invention includes the following steps:
步骤1:高速加工机床三维数字化建模:采用三维数字化软件,将拟设计的高速加工机床的初步结构转化为三维CAD模型;Step 1: Three-dimensional digital modeling of high-speed processing machine tools: The three-dimensional digital software is used to convert the preliminary structure of the proposed high-speed processing machine tools into three-dimensional CAD models;
步骤2:高速加工机床主要热源发热功率和相关换热系数计算:采用Palmgren经验公式计算滚动轴承发热功率,采用传热学理论计算相关部位对流换热系数;为高速加工机床整机热态建模提供基础数据;Step 2: Calculation of the heating power and related heat transfer coefficient of the main heat source of high-speed processing machine tools: Calculate the heating power of rolling bearings using Palmgren's empirical formula, and calculate the convective heat transfer coefficient of relevant parts by using heat transfer theory; Basic data;
步骤3:机床平面结合部热阻参数计算:采用分形几何理论,计算机床平面结合部接触热阻;Step 3: Calculation of thermal resistance parameters of the plane joint of the machine tool: using fractal geometry theory, contact thermal resistance of the plane joint of the computer bed;
步骤4:高速加工机床整机结构热力学建模与热特性计算:建立含主要热源与典型结合部接触热阻的高速加工机床的整机热力学模型;开展机床关键零部件、水冷系统物理参数对高速加工机床整机热力学灵敏度分析与计算,寻找影响机床热力学性能的薄弱环节;Step 4: Thermodynamic modeling and thermal characteristics calculation of the high-speed processing machine tool structure: Establish a complete machine thermodynamic model of the high-speed processing machine tool including the contact thermal resistance between the main heat source and the typical joint; develop the key parameters of the machine tool, and the physical parameters of the water-cooling system against Analysis and calculation of the thermodynamic sensitivity of the entire machine tool to find the weak links that affect the thermodynamic performance of the machine tool;
步骤5:高速加工机床整机结构热态设计方法:针对步骤4寻找到的影响机床热力学性能的薄弱环节,以提高机床热态性能为目标,提出高速加工机床整机结构动态设计方法。Step 5: The thermal design method of the high-speed processing machine tool structure: Aiming at the weak links affecting the thermodynamic performance of the machine tool found in step 4, aiming at improving the thermal performance of the machine tool, a dynamic design method of the high-speed processing machine tool machine structure is proposed.
所述步骤1:高速加工机床三维数字化建模,采用三维数字化软件Solidworks或者Pro-E,将拟设计的高速加工机床的初步结构转化为三维CAD模型。Step 1: Three-dimensional digital modeling of the high-speed processing machine tool. The three-dimensional digital software Solidworks or Pro-E is used to convert the preliminary structure of the proposed high-speed processing machine tool into a three-dimensional CAD model.
所述所述步骤2:高速加工机床主要热源发热功率和相关换热系数计算,分为以下三步:The step 2: the calculation of the heating power of the main heat source and the related heat transfer coefficient of the high-speed processing machine tool is divided into the following three steps:
步骤2a:高速加工机床电主轴与工件主轴发热功率计算Step 2a: Calculation of heating power of electric spindle and workpiece spindle of high-speed processing machine
1)滚动轴承发热功率Q 1计算, 1) Calculation of heating power Q 1 for rolling bearings,
Q 1=1.047×10 -4(M 0+M 1)·n         (1) Q 1 = 1.047 × 10 -4 (M 0 + M 1 ) · n (1)
Figure PCTCN2019078377-appb-000001
Figure PCTCN2019078377-appb-000001
M 1=f 1p 1d m            (3) M 1 = f 1 p 1 d m (3)
式中:M 0为润滑剂的流体动力损耗,M 1为弹性滞后和局部差动滑动的摩擦损耗,n为轴承转速,r/min;d m为轴承中径,mm;f 0为与轴承类型和润滑方式有关的经验常数,v为在工作温度下润滑剂的运动粘度,cst,f 1为与轴承类型和所受负荷有关的参数;p 1为确定轴承摩擦力矩的计算负荷N; In the formula: M 0 is the hydrodynamic loss of the lubricant, M 1 is the elastic hysteresis and the friction loss of the local differential sliding, n is the speed of the bearing, r / min; d m is the bearing diameter, mm; f 0 is the bearing Type and empirical constants related to the lubrication method, v is the kinematic viscosity of the lubricant at the working temperature, cst, f 1 are parameters related to the bearing type and the load; p 1 is the calculated load N to determine the bearing friction torque;
2)电主轴电机转子、定子发热功率计算2) Calculation of heating power of rotor and stator of electric spindle motor
P 0为额定功率,电机功率系数α,效率η,则热损耗Q=P 0α(1-η),其中热损耗2/3分布在电机转子,热损耗1/3分布在电机定子; P 0 is the rated power, the motor power coefficient α, and the efficiency η, then the heat loss Q = P 0 α (1-η), where 2/3 of the heat loss is distributed in the motor rotor, and 1/3 of the heat loss is distributed in the motor stator;
3)工件主轴带传动发热功率q belt计算 3) Calculation of the heating power q belt of the workpiece spindle belt
q belt=P beltin(1-η)           (3) q belt = P beltin (1-η) (3)
式中:P beltin带传动输入功率W,η为带传动效率,带传动的热损耗一半传递到皮带上, 一半传递到带轮上,且主动轮和从动轮上的摩擦生热按两带轮的皮带包角等比例分配; In the formula: P beltin belt transmission input power W, η is the belt transmission efficiency, half of the heat loss of the belt transmission is transmitted to the belt, and half is transmitted to the belt wheel, and the frictional heat generated on the driving wheel and the driven wheel is based on the two belt wheels. Angle distribution of belts
步骤2b:高速加工机床电主轴与工件主轴对流换热系数计算Step 2b: Calculation of convective heat transfer coefficient between high-speed machining machine electric spindle and workpiece spindle
换热系数α和努赛尔数Nu的关系:Relationship between heat transfer coefficient α and Nusser number Nu:
Figure PCTCN2019078377-appb-000002
Figure PCTCN2019078377-appb-000002
式中:λ为流体导热系数,Nu f为努赛尔数,D为几何特征的定型尺寸; Where: λ is the thermal conductivity of the fluid, Nu f is the Nusser number, and D is the fixed dimension of the geometric feature;
1)电主轴内置电机定子冷却水套强制对流换热系数h 1计算 1) Calculation of forced convection heat transfer coefficient h 1 of the stator cooling water jacket of the built-in motor of the electric spindle
Figure PCTCN2019078377-appb-000003
Figure PCTCN2019078377-appb-000003
式中:λ为流体导热系数,Nu f为努赛尔数;Pr f为普朗特数,Re f为雷诺数,L为冷却水流道长度; In the formula: λ is the thermal conductivity of the fluid, Nu f is the Nusser number; Pr f is the Prandtl number, Re f is the Reynolds number, and L is the length of the cooling water channel;
2)电主轴电机定子内表面与转子外表面的传热系数h 2计算 2) Calculation of the heat transfer coefficient h 2 of the inner surface of the stator and the outer surface of the rotor of the electric spindle motor
Figure PCTCN2019078377-appb-000004
Figure PCTCN2019078377-appb-000004
式中:r 1为转子外径(m);σ为定、转子间气隙(m); In the formula: r 1 is the outer diameter of the rotor (m); σ is fixed, the air gap between the rotors (m);
3)工件主轴旋转表面与空气间的对流换热系数h 3计算 3) Calculation of the convective heat transfer coefficient h 3 between the rotating surface of the workpiece spindle and the air
Figure PCTCN2019078377-appb-000005
Figure PCTCN2019078377-appb-000005
4)工件主轴带轮对流换热系数h 4计算 4) Convection heat transfer coefficient h 4 of workpiece spindle
h 4=9.7+5.33×U 0.8            (8) h 4 = 9.7 + 5.33 × U 0.8 (8)
式中:h 4—对流换热系数W/(m 2·K);U—带轮表面线速度m/s; In the formula: h 4 —Convection heat transfer coefficient W / (m 2 · K); U—Linear surface speed m / s;
5)主轴外表面与空气间的对流换热系数,取复合传热系数为h 5=9.7W/(m 2·K) 5) The convective heat transfer coefficient between the outer surface of the main shaft and the air, taking the composite heat transfer coefficient as h 5 = 9.7W / (m 2 · K)
步骤2c)切削液“二次热源”的建模方法Step 2c) Modeling method of "secondary heat source" of cutting fluid
切削液流道部分的对流换热系数Nu计算:Calculation of convection heat transfer coefficient Nu of cutting fluid flow channel part:
Figure PCTCN2019078377-appb-000006
Figure PCTCN2019078377-appb-000006
切削液溅到的表面的对流换热系数:将切削液流道对流系数乘以0.6,得到该表面的对流系数;式(9)的定性温度取流体和壁面的平均温度T m,定型尺寸取流体流过的壁面长度L。 Convection heat transfer coefficient of the surface to which the cutting fluid splashes: Multiply the convection coefficient of the cutting fluid channel by 0.6 to obtain the convection coefficient of the surface. The qualitative temperature of formula (9) is the average temperature T m of the fluid and the wall. Length L of wall surface through which fluid flows.
所述步骤3:高速加工机床主要结合部热阻参数计算,分为以下2步:Step 3: The calculation of the thermal resistance parameters of the main joint of the high-speed processing machine tool is divided into the following two steps:
步骤3a)固定结合面接触热阻R c的计算 Step 3a) Calculation of contact thermal resistance R c of the fixed bonding surface
总的接触热阻TCR为The total contact thermal resistance TCR is
Figure PCTCN2019078377-appb-000007
Figure PCTCN2019078377-appb-000007
式中:R ci为第i个微凸体接触热阻。 In the formula: R ci is the i-th micro-convex contact thermal resistance.
无量纲总接触热阻为The dimensionless total contact thermal resistance is
Figure PCTCN2019078377-appb-000008
Figure PCTCN2019078377-appb-000008
式中:R ci为第i个微凸体接触热阻,
Figure PCTCN2019078377-appb-000009
为相对面积接触比,
Figure PCTCN2019078377-appb-000010
为第i个微凸体的截面积,N(a′ s)为接触点数目;
Where R ci is the i-th micro-convex contact thermal resistance,
Figure PCTCN2019078377-appb-000009
Is the relative area contact ratio,
Figure PCTCN2019078377-appb-000010
Is the cross-sectional area of the i-th microconvex, and N (a ′ s ) is the number of contact points;
步骤3b)轴承外圈与轴承座间的接触热阻计算Step 3b) Calculation of the contact thermal resistance between the bearing outer ring and the bearing seat
轴承接触热导为,The bearing contact thermal conductivity is,
Figure PCTCN2019078377-appb-000011
Figure PCTCN2019078377-appb-000011
式中:Π为接触热导(W/K);h ring、h gap分别为轴承外圈、外圈与轴承座平均间隙厚度m;λ ring、λ gap分别为轴承外圈和间隙中空气的导热系数W/(m 2·K),A为轴承套圈的圆柱形外表面面积; Where: Π is the contact thermal conductivity (W / K); h ring and h gap are the average thickness of the bearing outer ring, the outer ring and the bearing seat, respectively; λ ring and λ gap are the air in the bearing outer ring and the gap, respectively. Thermal conductivity W / (m 2 · K), A is the cylindrical outer surface area of the bearing ring;
平均间隙h gap可由下式计算 The average gap h gap can be calculated by the following formula
h gap=h gap0-(T ring-T housing)a′·r housing        (13) h gap = h gap0- (T ring -T housing ) a′r housing (13)
式中:h gap为轴承外圈与轴承座之间的原始间隙,m;a′为线性膨胀系数,对于钢材可取其值为11.7×10 -6/K。 In the formula: h gap is the original gap between the bearing outer ring and the bearing seat, m; a ′ is the linear expansion coefficient, which can be 11.7 × 10 -6 / K for steel.
所述步骤4:高速加工机床整机结构热力学建模与热特性计算,建立含主要热源与典型结合部接触热阻的高速加工机床的整机热力学模型;开展高速加工机床主要热源的发热功率、机床关键零部件、水冷系统物理参数对高速加工机床整机热力学性能的灵敏度分析计算,寻找影响机床热性能的薄弱环节与主要因素。Step 4: Thermodynamic modeling and thermal characteristic calculation of the high-speed processing machine tool structure, establishing a complete machine thermodynamic model of the high-speed processing machine tool including the thermal resistance of the main heat source and the typical joint contact thermal resistance; developing the heating power of the main heat source of the high-speed processing machine tool, Sensitivity analysis and calculation of key machine parts and water cooling system physical parameters to the thermodynamic performance of high-speed processing machine tools, looking for weak links and main factors that affect the thermal performance of machine tools.
所述主要热源为高速电主轴、工件主轴、切削液二次热源;典型结合部为平面结合部、滚动结合部。The main heat source is a high-speed electric main shaft, a workpiece main shaft, and a secondary heat source for cutting fluid; typical joints are a planar joint and a rolling joint.
所述高速加工机床整机热力学性能包括温度场与热位移场。The thermodynamic performance of the high-speed processing machine tool includes a temperature field and a thermal displacement field.
所述步骤5:高速加工机床整机结构热态设计方法,围绕步骤4发现的影响机床热性能的薄弱环节与主要因素,以降低加工过程中的热变形、提高加工精度为设计目标,修改机床主要零部件的结构参数,重新设计相关物理参数,完成高速加工机床热设计。Step 5: The thermal state design method of the high-speed processing machine tool structure, surrounding the weak links and main factors affecting the thermal performance of the machine tool found in step 4, with the goal of reducing thermal deformation and improving processing accuracy during the process of modifying the machine tool Structural parameters of main components, redesign related physical parameters, and complete thermal design of high-speed processing machine tools.
有益效果:本发明提供的高速加工机床整机结构热力学建模与热设计方法,是基于完善的高速加工机床热力学模型与系统的热力学分析计算结果而提出的。因此采用本发明“高速加工机床整机结构热力学建模与热设计方法”,可以提高机床热态设计准确性和一次设计成功率。Beneficial effect: The thermodynamic modeling and thermal design method of the high-speed processing machine tool structure provided by the present invention is based on the thermodynamic analysis and calculation results of a complete high-speed processing machine tool thermodynamic model and system. Therefore, by adopting the "thermodynamic modeling and thermal design method for the whole structure of a high-speed processing machine tool" of the present invention, the thermal design accuracy of the machine tool and the design success rate at one time can be improved.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
图1本发明专利的设计步骤与内容,Figure 1 Design steps and content of the invention patent,
图2高速加工机床三维数字化模型,Figure 2 Three-dimensional digital model of high-speed processing machine,
a为内圆复合磨床3D-CAD模型,a is a 3D-CAD model of the internal compound grinding machine,
b为整机有限元分析模型,b is the finite element analysis model of the whole machine,
图3无量纲接触热阻R *在不同参数影响下计算结果, Fig. 3 Dimensionless contact thermal resistance R * Calculation results under the influence of different parameters,
a为无量纲接触热阻与接触载荷在不同G *下的关系(E/σ=106,D=1.3), a is the relationship between the dimensionless contact thermal resistance and the contact load under different G * (E / σ = 106, D = 1.3),
b为无量纲接触热阻与接触载荷在不同D下的关系(E/σ=106,G=10 -7), b is the relationship between the dimensionless contact thermal resistance and the contact load under different D (E / σ = 106, G = 10 -7 ),
c为无量纲接触热阻与接触载荷在不同E/σ下的关系(D=1.3,G=10 -7), c is the relationship between the dimensionless contact thermal resistance and the contact load at different E / σ (D = 1.3, G = 10 -7 ),
d为无量纲接触面积与接触热阻在不同接触变形情形下的关系,d is the relationship between the dimensionless contact area and the contact thermal resistance under different contact deformation situations,
图4整机额定工况下的温度场分布云图,Figure 4 Cloud diagram of temperature field distribution under rated conditions of the whole machine,
a为等距测视角的整机额定工况温度场分布云图,a is the cloud diagram of the temperature range of the entire machine under rated working conditions at an isometric perspective,
b为背视图视角的整机额定工况温度场分布云图,b is the distribution cloud diagram of the rated operating temperature field of the whole machine from the rear view perspective,
图5整机额定工况下的热位移场分布云图,Fig. 5 Distribution diagram of thermal displacement field under rated conditions of the whole machine,
a为等距测视角的整机额定工况热位移场分布云图,a is the cloud map of the thermal displacement field of the entire machine under the rated operating conditions at an isometric perspective,
b为背视图视角的整机额定工况热位移场分布云图,b is the cloud diagram of the thermal displacement field distribution of the complete machine under rated operating conditions,
图6热设计后的整机温度场分布云图,Figure 6 Cloud diagram of the whole machine temperature field after thermal design.
a为俯视图视角的热设计后的整机温度场分布云图,a is the cloud map of the whole machine temperature field after thermal design from the perspective of the top view,
b为背视图视角的热设计后的整机温度场分布云图。b is the cloud map of the whole machine's temperature field after thermal design from the perspective of the back view.
具体实施方式detailed description
下面结合一个实施例(精密数控内圆复合磨床),对本发明的高速加工机床整机结构热力学设计方法作进一步详细说明。In the following, an embodiment (precision numerical control inner circle composite grinder) is used to further describe the thermodynamic design method of the entire structure of the high-speed processing machine tool of the present invention.
图1给出了本发明高速加工机床整机结构热力学设计方法的具体内容,包括如下步骤:FIG. 1 shows the specific content of the thermodynamic design method of the high-speed processing machine tool structure of the present invention, including the following steps:
步骤1:高速加工机床三维数字化建模Step 1: 3D digital modeling of high-speed processing machine
采用三维建模软件(如Solidworks)建立精密数控内圆复合磨床的三维数字化(CAD、CAE)模型,如图2所示。3D modeling software (such as Solidworks) is used to establish a 3D digital (CAD, CAE) model of a precision CNC internal cylindrical composite grinder, as shown in Figure 2.
步骤2:高速加工机床主要热源发热功率计算Step 2: Calculation of heating power of main heat source of high-speed processing machine
步骤2a:采用公式(1)-(3),计算高速加工机床磨削电主轴滚动轴承发热功率Q,在不同温度情况下,计算结果如表1所示Step 2a: Use formulas (1)-(3) to calculate the heating power Q of the high-speed machining machine grinding electric spindle rolling bearing. The calculation results are shown in Table 1 at different temperatures.
表1 内圆磨削电主轴轴承生热率q与温度的对应关系(q=Q/V)Table 1 Correspondence between the heat generation rate q and temperature of the electric spindle bearing for internal grinding (q = Q / V)
Figure PCTCN2019078377-appb-000012
Figure PCTCN2019078377-appb-000012
Figure PCTCN2019078377-appb-000013
Figure PCTCN2019078377-appb-000013
采用公式(1)-(3),计算高速加工机床工件主轴滚动轴承发热功率Q,在不同温度情况下,得到不同温度下轴承生热率,见表2所示。Using formulas (1)-(3), calculate the heating power Q of the rolling bearing of the spindle of the high-speed machining machine tool. Under different temperature conditions, the bearing heat generation rate at different temperatures is obtained, as shown in Table 2.
表2 轴承生热率与温度的对应关系(q=Q/V)Table 2 Correspondence between bearing heat generation rate and temperature (q = Q / V)
Figure PCTCN2019078377-appb-000014
Figure PCTCN2019078377-appb-000014
电主轴内置电机的额定功率为16kW,正常工作状态下功率系数为0.9,功率损耗为1.5kW。假设该损耗全部转化成热量,其中2/3(1kW)由定子发出,1/3(0.5kW)由转子发出,则电机定、转子生热率分别为254250W/m 3、255056W/m 3The electric motor's built-in motor has a rated power of 16kW, a power factor of 0.9 under normal working conditions, and a power loss of 1.5kW. Assuming that all the losses are converted into heat, of which 2/3 (1kW) is emitted by the stator and 1/3 (0.5kW) is emitted by the rotor, the motor heating and rotor heat generation rates are 254250W / m 3 and 255056W / m 3 respectively .
带传动的效率为96%,带传动输入功率P beltin为可倾瓦动压轴承和推力滑动轴承等的摩擦损耗功率,为5.5kW,工件主轴上从动带轮的摩擦热损(q belt)为43.86W。 The efficiency of the belt drive is 96%. The belt drive input power P beltin is the friction loss power of the tilting pad dynamic pressure bearing and thrust sliding bearing, which is 5.5kW. The frictional heat loss of the driven pulley on the workpiece spindle (q belt ) It is 43.86W.
步骤2b:高速加工机床电主轴与工件主轴对流换热系数计算Step 2b: Calculation of convective heat transfer coefficient between high-speed machining machine electric spindle and workpiece spindle
采用公式(5)计算强制冷却水对流换热系数h 1=13719W/(m 2·K) Use formula (5) to calculate the convective heat transfer coefficient of forced cooling water h 1 = 13719W / (m 2 · K)
采用公式(6)电机内置转子、定子对流换热系数h 2=116.2W/(m 2·K) Adopt formula (6) Convection heat transfer coefficient of motor with built-in rotor and stator h 2 = 116.2W / (m 2 · K)
采用公式(7)工件主轴旋转表面与空气间的对流换热系数计算h 3=100W/(m 2·K) Use formula (7) to calculate the convective heat transfer coefficient between the rotating surface of the workpiece spindle and the air h 3 = 100W / (m 2 · K)
根据公式(8),皮带轮直径为245mm、200mm和144mm的径向表面所对应的对流换热系数分别为81W/(m 2·K)、70W/(m 2·K)和56W/(m 2·K);砂轮主轴直径为95mm和56mm的径向表面所对应的对流换热系数分别为43W/(m 2·K)和31W/(m 2·K) According to formula (8), the convection heat transfer coefficients corresponding to the radial surfaces of the pulley diameters of 245mm, 200mm, and 144mm are 81W / (m 2 · K), 70W / (m 2 · K), and 56W / (m 2 · K); the convective heat transfer coefficients of the radial surfaces of the grinding wheel with a diameter of 95mm and 56mm are 43W / (m 2 · K) and 31W / (m 2 · K)
主轴外表面与空气间的对流换热系数,取复合传热系数为h 5=9.7W/(m 2·K) Convective heat transfer coefficient between the outer surface of the main shaft and the air, taking the composite heat transfer coefficient as h 5 = 9.7W / (m 2 · K)
步骤2c:采用公式(9)计算高速加工机床切削液作为“二次热源”相关对流换热系数计算,额定工况下切削液流量100L/min。主要表面的对流换热换热系数见表3所示:Step 2c: Use formula (9) to calculate the relevant convective heat transfer coefficient of the cutting fluid of the high-speed processing machine tool as a "secondary heat source". The flow rate of the cutting fluid is 100 L / min under rated conditions. The convection heat transfer coefficients of the main surfaces are shown in Table 3:
表3 床身主要表面的对流换热系数Table 3 Convection heat transfer coefficients of the main surfaces of the bed
Figure PCTCN2019078377-appb-000015
Figure PCTCN2019078377-appb-000015
Figure PCTCN2019078377-appb-000016
Figure PCTCN2019078377-appb-000016
步骤3:机床平面结合部热阻参数计算Step 3: Calculation of thermal resistance parameters of the plane joint of the machine tool
步骤3a:采用公式(10)(11)计算总的接触热阻R和无量纲接触热阻R *,计算结果如图3所示。 Step 3a: Use formulas (10) and (11) to calculate the total contact thermal resistance R and the dimensionless contact thermal resistance R * . The calculation results are shown in FIG.
步骤3b:采用公式(12)计算轴承外圈与轴承座间的接触热阻,对于普通主轴轴承(取4级加工精度),取h ring=5×10 -6m,λ ring=24W/(m·K),λ gap=2.84×10 -2W/(m·K),则主轴轴承与配合面的接触系数见表4所示。 Step 3b: Use formula (12) to calculate the contact thermal resistance between the outer ring of the bearing and the bearing block. For ordinary spindle bearings (take 4 levels of machining accuracy), take h ring = 5 × 10 -6 m, λ ring = 24W / ( m · K), λ gap = 2.84 × 10 -2 W / (m · K), then the contact coefficient between the spindle bearing and the mating surface is shown in Table 4.
表4 主轴轴承与配合面间的接触系数h(单位:W/m 2·K) Table 4 Contact coefficient h between the main shaft bearing and the mating surface (unit: W / m 2 · K)
Figure PCTCN2019078377-appb-000017
Figure PCTCN2019078377-appb-000017
步骤4:高速加工机床整机结构热力学建模与热特性计算Step 4: Thermodynamic Modeling and Thermal Characteristic Calculation of High Speed Machining Machine Structure
本实施例,内圆磨削孔径为150mm,砂轮直径为80mm。将内圆磨削电主轴、床头箱和切削液等热源,全部施加于机床热力学分析模型,计算得到整机在额定工况下的温度场分布云图,见图4。In this embodiment, the internal grinding hole diameter is 150 mm, and the grinding wheel diameter is 80 mm. All internal heat sources such as electric spindle, bedside box and cutting fluid are applied to the thermodynamic analysis model of the machine tool, and the temperature field distribution cloud diagram of the whole machine under rated working conditions is calculated, as shown in Figure 4.
整机在额定工况下的热位移场分布云图,见图5,图中黑色网格表示未发生热变形时的轮廓。The cloud diagram of the thermal displacement field of the whole machine under rated conditions is shown in Figure 5. The black grid in the figure indicates the outline when no thermal deformation has occurred.
额定工况下,整机主要部位的热位移,列于表5。The thermal displacements of the main parts of the machine under rated conditions are listed in Table 5.
表5 整机在额定工况下各构件平均热位移Table 5 The average thermal displacement of each component of the machine under rated conditions
Figure PCTCN2019078377-appb-000018
Figure PCTCN2019078377-appb-000018
经换算,整机热变形引起的内孔和锥面磨削误差为43.30μm。After conversion, the internal hole and tapered surface grinding error caused by the thermal deformation of the whole machine is 43.30 μm.
步骤5:高速加工机床整机结构热态设计方法Step 5: Thermal design method of high-speed processing machine tool structure
针对步骤4分析结果,本实施例对机床结构进行改进设计,将与流道尺寸相同的平板覆盖在流道上,使切削液不直接与床身接触,减少传递给床身热量,进而减少床身热变形,其主要部位的热位移列于表6。According to the analysis result of step 4, this embodiment improves the design of the machine tool structure. The flat plate with the same size as the flow path is covered on the flow path, so that the cutting fluid does not directly contact the bed, and the heat transmitted to the bed is reduced, thereby reducing the bed. Thermal deformation, the thermal displacement of the main parts are listed in Table 6.
表6 热设计之后各构件的平均热位移Table 6 Average thermal displacement of each component after thermal design
Figure PCTCN2019078377-appb-000019
Figure PCTCN2019078377-appb-000019
经热设计之后,整机热变形引起的内孔和锥面磨削误差为1.15μm。After thermal design, the grinding error of the inner hole and the tapered surface caused by the thermal deformation of the whole machine is 1.15 μm.
由设计结果看,在磨削孔径为150mm,内圆磨削砂轮直径为80mm状态下,对磨床进行热设计前、后的加工热误差对比,见表7。可以看出,采用本发明专利“高速加工机床整机结构热力学建模与热设计方法”,可以大幅度提高机床的磨削精度。From the design results, under the condition that the grinding hole diameter is 150mm and the internal grinding wheel diameter is 80mm, the comparison of the processing thermal errors before and after the thermal design of the grinding machine is shown in Table 7. It can be seen that the use of the invention patent "the method of thermodynamic modeling and thermal design of the whole structure of a high-speed processing machine tool" can greatly improve the grinding accuracy of the machine tool.
表7 热设计前后加工精度对比Table 7 Comparison of processing accuracy before and after thermal design
Figure PCTCN2019078377-appb-000020
Figure PCTCN2019078377-appb-000020

Claims (8)

  1. 一种高速加工机床整机结构热力学建模与热设计方法,其特征在于该方法包括以下步骤:A high-speed processing machine tool structure thermodynamic modeling and thermal design method is characterized in that the method includes the following steps:
    步骤1:高速加工机床三维数字化建模:采用三维数字化软件,将拟设计的高速加工机床的初步结构转化为三维CAD模型;Step 1: Three-dimensional digital modeling of high-speed processing machine tools: The three-dimensional digital software is used to convert the preliminary structure of the proposed high-speed processing machine tools into three-dimensional CAD models;
    步骤2:高速加工机床主要热源发热功率和相关换热系数计算:采用Palmgren经验公式计算滚动轴承发热功率,采用传热学理论计算相关部位对流换热系数;为高速加工机床整机热态建模提供基础数据;Step 2: Calculation of the heating power and related heat transfer coefficient of the main heat source of high-speed processing machine tools: Calculate the heating power of rolling bearings using Palmgren's empirical formula, and calculate the convective heat transfer coefficient of relevant parts by using heat transfer theory; Basic data;
    步骤3:机床平面结合部热阻参数计算:采用分形几何理论,计算机床平面结合部接触热阻;Step 3: Calculation of thermal resistance parameters of the plane joint of the machine tool: using fractal geometry theory, contact thermal resistance of the plane joint of the computer bed;
    步骤4:高速加工机床整机结构热力学建模与热特性计算:建立含主要热源与典型结合部接触热阻的高速加工机床的整机热力学模型;开展机床关键零部件、水冷系统物理参数对高速加工机床整机热力学灵敏度分析与计算,寻找影响机床热力学性能的薄弱环节;Step 4: Thermodynamic modeling and thermal characteristics calculation of the high-speed processing machine tool structure: Establish a complete machine thermodynamic model of the high-speed processing machine tool including the contact thermal resistance between the main heat source and the typical joint; develop the key parameters of the machine tool, and the physical parameters of the water-cooling system against high-speed Analysis and calculation of the thermodynamic sensitivity of the entire machine tool to find the weak links that affect the thermodynamic performance of the machine tool;
    步骤5:高速加工机床整机结构热态设计方法:针对步骤4寻找到的影响机床热力学性能的薄弱环节,以提高机床热态性能为目标,提出高速加工机床整机结构动态设计方法。Step 5: The thermal design method of the high-speed processing machine tool structure: Aiming at the weak links affecting the thermodynamic performance of the machine tool found in step 4, aiming at improving the thermal performance of the machine tool, a dynamic design method of the high-speed processing machine tool machine structure is proposed.
  2. 根据权利要求1所述的高速加工机床整机结构热力学建模与热设计方法,其特征在于,所述步骤1:高速加工机床三维数字化建模,采用三维数字化软件Solidworks或者Pro-E,将拟设计的高速加工机床的初步结构转化为三维CAD模型。The method for thermodynamic modeling and thermal design of a high-speed processing machine tool structure according to claim 1, wherein the step 1: three-dimensional digital modeling of the high-speed processing machine tool, using three-dimensional digital software Solidworks or Pro-E, The preliminary structure of the designed high-speed machining machine was transformed into a three-dimensional CAD model.
  3. 根据权利要求1所述的高速加工机床整机结构热力学建模与热设计方法,其特征在于,所述所述步骤2:高速加工机床主要热源发热功率和相关换热系数计算,分为以下三步:The method for thermodynamic modeling and thermal design of a high-speed processing machine tool structure according to claim 1, wherein said step 2: the calculation of the heating power and the related heat transfer coefficient of the main heat source of the high-speed processing machine tool are divided into the following three step:
    步骤2a:高速加工机床电主轴与工件主轴发热功率计算Step 2a: Calculation of heating power of electric spindle and workpiece spindle of high-speed processing machine
    1)滚动轴承发热功率Q 1计算, 1) Calculation of heating power Q 1 for rolling bearings,
    Q 1=1.047×10 -4(M 0+M 1)·n          (1) Q 1 = 1.047 × 10 -4 (M 0 + M 1 ) · n (1)
    Figure PCTCN2019078377-appb-100001
    Figure PCTCN2019078377-appb-100001
    M 1=f 1p 1d m       (3) M 1 = f 1 p 1 d m (3)
    式中:M 0为润滑剂的流体动力损耗,M 1为弹性滞后和局部差动滑动的摩擦损耗,n为轴承转速,r/min;d m为轴承中径,mm;f 0为与轴承类型和润滑方式有关的经验常数,v为在工作温度下润滑剂的运动粘度,cst,f 1为与轴承类型和所受负荷有关的参数;p 1为确定轴承摩擦力矩的计算负荷N; In the formula: M 0 is the hydrodynamic loss of the lubricant, M 1 is the elastic hysteresis and the friction loss of the local differential sliding, n is the speed of the bearing, r / min; d m is the bearing diameter, mm; f 0 is the bearing Type and empirical constants related to the lubrication method, v is the kinematic viscosity of the lubricant at the working temperature, cst, f 1 are parameters related to the bearing type and the load; p 1 is the calculated load N to determine the bearing friction torque;
    2)电主轴电机转子、定子发热功率计算2) Calculation of heating power of rotor and stator of electric spindle motor
    P 0为额定功率,电机功率系数α,效率η,则热损耗Q=P 0α(1-η),其中热损耗2/3分布在电机转子,热损耗1/3分布在电机定子; P 0 is the rated power, the motor power coefficient α, and the efficiency η, then the heat loss Q = P 0 α (1-η), where 2/3 of the heat loss is distributed in the motor rotor, and 1/3 of the heat loss is distributed in the motor stator;
    3)工件主轴带传动发热功率q belt计算 3) Calculation of the heating power q belt of the workpiece spindle belt
    q belt=P beltin(1-η)          (3) q belt = P beltin (1-η) (3)
    式中:P beltin带传动输入功率W,η为带传动效率,带传动的热损耗一半传递到皮带上,一半传递到带轮上,且主动轮和从动轮上的摩擦生热按两带轮的皮带包角等比例分配; In the formula: P beltin belt transmission input power W, η is the belt transmission efficiency, half of the heat loss of the belt transmission is transmitted to the belt, and half is transmitted to the pulley, and the frictional heat generation on the driving and driven wheels is based on the two pulleys. Angle distribution of belts
    步骤2b:高速加工机床电主轴与工件主轴对流换热系数计算Step 2b: Calculation of convective heat transfer coefficient between high-speed machining machine electric spindle and workpiece spindle
    换热系数α和努赛尔数Nu的关系:Relationship between heat transfer coefficient α and Nusser number Nu:
    Figure PCTCN2019078377-appb-100002
    Figure PCTCN2019078377-appb-100002
    式中:λ为流体导热系数,Nu f为努赛尔数,D为几何特征的定型尺寸; Where: λ is the thermal conductivity of the fluid, Nu f is the Nusser number, and D is the fixed dimension of the geometric feature;
    1)电主轴内置电机定子冷却水套强制对流换热系数h 1计算 1) Calculation of forced convection heat transfer coefficient h 1 of the stator cooling water jacket of the built-in motor of the electric spindle
    Figure PCTCN2019078377-appb-100003
    Figure PCTCN2019078377-appb-100003
    式中:λ为流体导热系数,Nu f为努赛尔数;Pr f为普朗特数,Re f为雷诺数,L为冷却水流道长度; In the formula: λ is the thermal conductivity of the fluid, Nu f is the Nusser number; Pr f is the Prandtl number, Re f is the Reynolds number, and L is the length of the cooling water channel;
    2)电主轴电机定子内表面与转子外表面的传热系数h 2计算 2) Calculation of the heat transfer coefficient h 2 of the inner surface of the stator and the outer surface of the rotor of the electric spindle motor
    Figure PCTCN2019078377-appb-100004
    Figure PCTCN2019078377-appb-100004
    式中:r 1为转子外径(m);σ为定、转子间气隙(m); In the formula: r 1 is the outer diameter of the rotor (m); σ is fixed, the air gap between the rotors (m);
    3)工件主轴旋转表面与空气间的对流换热系数h 3计算 3) Calculation of the convective heat transfer coefficient h 3 between the rotating surface of the workpiece spindle and the air
    Figure PCTCN2019078377-appb-100005
    Figure PCTCN2019078377-appb-100005
    4)工件主轴带轮对流换热系数h 4计算 4) Convection heat transfer coefficient h 4 of workpiece spindle
    h 4=9.7+5.33×U 0.8          (8) h 4 = 9.7 + 5.33 × U 0.8 (8)
    式中:h 4—对流换热系数W/(m 2·K);U—带轮表面线速度m/s; In the formula: h 4 —Convection heat transfer coefficient W / (m 2 · K); U—Linear surface speed m / s;
    5)主轴外表面与空气间的对流换热系数,取复合传热系数为h 5=9.7W/(m 2·K) 5) The convective heat transfer coefficient between the outer surface of the main shaft and the air, taking the composite heat transfer coefficient as h 5 = 9.7W / (m 2 · K)
    步骤2c)切削液“二次热源”的建模方法Step 2c) Modeling method of "secondary heat source" of cutting fluid
    切削液流道部分的对流换热系数Nu计算:Calculation of convection heat transfer coefficient Nu of cutting fluid flow channel part:
    Figure PCTCN2019078377-appb-100006
    Figure PCTCN2019078377-appb-100006
    切削液溅到的表面的对流换热系数:将切削液流道对流系数乘以0.6,得到该表面的对流系数;式(9)的定性温度取流体和壁面的平均温度T m,定型尺寸取流体流过的壁面长度L。 Convection heat transfer coefficient of the surface to which the cutting fluid splashes: Multiply the convection coefficient of the cutting fluid channel by 0.6 to obtain the convection coefficient of the surface. The qualitative temperature of formula (9) is the average temperature T m of the fluid and the wall. Length L of wall surface through which fluid flows.
  4. 根据权利要求1所述的高速加工机床整机结构热力学建模与热设计方法,其特征在于,所述步骤3:高速加工机床主要结合部热阻参数计算,分为以下2步:The method for thermodynamic modeling and thermal design of a high-speed processing machine tool structure according to claim 1, wherein the step 3: calculating the thermal resistance parameters of the main joint of the high-speed processing machine tool is divided into the following two steps:
    步骤3a)固定结合面接触热阻R c的计算 Step 3a) Calculation of contact thermal resistance R c of the fixed bonding surface
    总的接触热阻TCR为The total contact thermal resistance TCR is
    Figure PCTCN2019078377-appb-100007
    Figure PCTCN2019078377-appb-100007
    式中:R ci为第i个微凸体接触热阻。 In the formula: R ci is the i-th micro-convex contact thermal resistance.
    无量纲总接触热阻为The dimensionless total contact thermal resistance is
    Figure PCTCN2019078377-appb-100008
    Figure PCTCN2019078377-appb-100008
    式中:R ci为第i个微凸体接触热阻,
    Figure PCTCN2019078377-appb-100009
    为相对面积接触比,
    Figure PCTCN2019078377-appb-100010
    为第i个微凸体的截面积,N(a′ s)为接触点数目;
    Where R ci is the i-th micro-convex contact thermal resistance,
    Figure PCTCN2019078377-appb-100009
    Is the relative area contact ratio,
    Figure PCTCN2019078377-appb-100010
    Is the cross-sectional area of the i-th microconvex, and N (a ′ s ) is the number of contact points;
    步骤3b)轴承外圈与轴承座间的接触热阻计算Step 3b) Calculation of the contact thermal resistance between the bearing outer ring and the bearing seat
    轴承接触热导为,The bearing contact thermal conductivity is,
    Figure PCTCN2019078377-appb-100011
    Figure PCTCN2019078377-appb-100011
    式中:Π为接触热导(W/K);h ring、h gap分别为轴承外圈、外圈与轴承座平均间隙厚度m;λ ring、λ gap分别为轴承外圈和间隙中空气的导热系数W/(m 2·K),A为轴承套圈的圆柱形外表面面积; Where: Π is the contact thermal conductivity (W / K); h ring and h gap are the average thickness of the bearing outer ring, the outer ring and the bearing seat, respectively; λ ring and λ gap are the air in the bearing outer ring and the gap, respectively. Thermal conductivity W / (m 2 · K), A is the cylindrical outer surface area of the bearing ring;
    平均间隙h gap可由下式计算 The average gap h gap can be calculated by the following formula
    h gap=h gap0-(T ring-T housing)a′·r housing     (13) h gap = h gap0- (T ring -T housing ) a′r housing (13)
    式中:h gap为轴承外圈与轴承座之间的原始间隙,m;a′为线性膨胀系数,对于钢材可取其值为11.7×10 -6/K。 In the formula: h gap is the original gap between the bearing outer ring and the bearing seat, m; a ′ is the linear expansion coefficient, which can be 11.7 × 10 -6 / K for steel.
  5. 根据权利要求1所述的高速加工机床整机结构热力学建模与热设计方法,其特征在于,所述步骤4:高速加工机床整机结构热力学建模与热特性计算,建立含主要热源与典型结合部接触热阻的高速加工机床的整机热力学模型;开展高速加工机床主要热源的发热功率、机床关键零部件、水冷系统物理参数对高速加工机床整机热力学性能的灵敏度分析计算,寻找影响机床热性能的薄弱环节与主要因素。The method for thermodynamic modeling and thermal design of a high-speed processing machine tool structure according to claim 1, characterized in that, said step 4: thermodynamic modeling and thermal characteristic calculation of the high-speed processing machine tool structure, establishing the main heat source and typical The whole machine thermodynamic model of the high-speed processing machine tool with contact thermal resistance at the joint; the sensitivity analysis and calculation of the thermal power of the main heat source of the high-speed machine tool, the key components of the machine tool, and the physical parameters of the water-cooling system on the thermodynamic performance of the high-speed machine tool are sought to find the effect on the machine tool Weak links and main factors of thermal performance.
  6. 根据权利要求5所述的高速加工机床整机结构热力学建模与热设计方法,其特征在于,所述主要热源为高速电主轴、工件主轴、切削液二次热源;典型结合部为平面结合部、滚动结合部。The method for thermodynamic modeling and thermal design of a high-speed processing machine tool structure according to claim 5, wherein the main heat source is a high-speed electric spindle, a workpiece spindle, and a secondary heat source for cutting fluid; a typical joint is a planar joint Scrolling joint.
  7. 根据权利要求5所述的高速加工机床整机结构热力学建模与热设计方法,其特征在于, 所述高速加工机床整机热力学性能包括温度场与热位移场。The method for thermodynamic modeling and thermal design of a high-speed processing machine tool structure according to claim 5, wherein the thermodynamic performance of the high-speed processing machine tool includes a temperature field and a thermal displacement field.
  8. 根据权利要求1所述的高速加工机床整机结构热力学建模与热设计方法,其特征在于,所述步骤5:高速加工机床整机结构热态设计方法,围绕步骤4发现的影响机床热性能的薄弱环节与主要因素,以降低加工过程中的热变形、提高加工精度为设计目标,修改机床主要零部件的结构参数,重新设计相关物理参数,完成高速加工机床热设计。The method for thermodynamic modeling and thermal design of a high-speed processing machine tool structure according to claim 1, characterized in that said step 5: a method for thermal state design of the high-speed processing machine tool structure, affecting the thermal performance of the machine tool found around step 4 Weak links and main factors are designed to reduce thermal deformation during processing and improve processing accuracy, modify the structural parameters of the main parts of the machine tool, redesign the relevant physical parameters, and complete the thermal design of the high-speed processing machine tool.
PCT/CN2019/078377 2018-09-07 2019-03-15 Thermodynamic modeling and thermal design method for overall structure of high-speed processing machine WO2020048110A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201811041397.3A CN109376377B (en) 2018-09-07 2018-09-07 High-speed processing machine tool complete machine structure Thermodynamics modeling and Thermal design
CN201811041397.3 2018-09-07

Publications (1)

Publication Number Publication Date
WO2020048110A1 true WO2020048110A1 (en) 2020-03-12

Family

ID=65404726

Family Applications (2)

Application Number Title Priority Date Filing Date
PCT/CN2018/104779 WO2020047880A1 (en) 2018-09-07 2018-09-10 Thermodynamic modeling and thermal design method for overall structure of high-speed machine tool
PCT/CN2019/078377 WO2020048110A1 (en) 2018-09-07 2019-03-15 Thermodynamic modeling and thermal design method for overall structure of high-speed processing machine

Family Applications Before (1)

Application Number Title Priority Date Filing Date
PCT/CN2018/104779 WO2020047880A1 (en) 2018-09-07 2018-09-10 Thermodynamic modeling and thermal design method for overall structure of high-speed machine tool

Country Status (2)

Country Link
CN (1) CN109376377B (en)
WO (2) WO2020047880A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112364444A (en) * 2020-09-23 2021-02-12 江苏赛洋机电科技有限公司 Numerical control machine tool temperature measuring point optimization method based on finite element model numerical simulation

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109376377B (en) * 2018-09-07 2019-10-11 东南大学 High-speed processing machine tool complete machine structure Thermodynamics modeling and Thermal design
CN110076974B (en) * 2019-04-18 2021-04-27 甘肃机电职业技术学院 Design method of injection mold conformal cooling water channel based on additive manufacturing technology
CN110705147B (en) * 2019-09-18 2023-08-04 北京工业大学 Numerical control machine tool spindle thermal state characteristic comprehensive theoretical modeling and analysis method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102063549A (en) * 2011-01-07 2011-05-18 西安交通大学 Method for stimulating fluid-solid-heat coupling value of machine tool spindle
CN103366050A (en) * 2013-06-25 2013-10-23 华中科技大学 High-speed electric spindle power-heat coupling modeling method
CN108241787A (en) * 2018-01-12 2018-07-03 哈尔滨理工大学 Static pressure rotary worktable thermal characteristics research method under extreme operating condition
CN109376377A (en) * 2018-09-07 2019-02-22 东南大学 High-speed processing machine tool complete machine structure Thermodynamics modeling and Thermal design

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3405965B2 (en) * 2000-07-06 2003-05-12 ファナック株式会社 Thermal displacement compensation method for machine tools
US9026413B2 (en) * 2009-09-15 2015-05-05 Invensys Systems, Inc. Thermodynamic phase equilibrium analysis based on a reduced composition domain
CN102452020B (en) * 2010-10-22 2016-08-10 西安交通大学 A kind of cutting tool for CNC machine temperature field and thermal deformation quantitative analysis method
CN102009370B (en) * 2010-11-03 2012-07-11 北京航空航天大学 Feedback compensation method of high-speed cutting CNC (computer numerical control) machine tool
CN102867088B (en) * 2012-09-12 2015-04-29 西安交通大学 Thermal network modeling method applied to electric spindle steady temperature field
CN103034169B (en) * 2012-12-26 2014-11-26 西安理工大学 Modeling and compensation method of heat error of numerical control machine tool
CN105022344B (en) * 2014-04-22 2017-11-28 上海开通数控有限公司 Numerical control machine heat error compensation method
CN104200019B (en) * 2014-08-27 2018-03-30 电子科技大学 A kind of heavy type numerical control floor type boring and milling machine thermal characteristics numerical value emulation method
CN105631122B (en) * 2015-12-25 2021-08-24 鼎奇(天津)主轴科技有限公司 Thermal deformation simulation analysis and modeling method for machine tool large part
CN106407528B (en) * 2016-09-05 2019-05-24 西安交通大学 A kind of machine tool mainshaft bearing heat production calculation method based on kinetic model
CN107498391B (en) * 2017-08-03 2019-07-02 南京航空航天大学 The heating power of machine cut process is load simulated and error detecting system and method
CN108304657B (en) * 2018-02-02 2020-11-10 重庆大学 Finite element-based multi-process continuous modeling simulation method for residual stress of machine tool weight parts

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102063549A (en) * 2011-01-07 2011-05-18 西安交通大学 Method for stimulating fluid-solid-heat coupling value of machine tool spindle
CN103366050A (en) * 2013-06-25 2013-10-23 华中科技大学 High-speed electric spindle power-heat coupling modeling method
CN108241787A (en) * 2018-01-12 2018-07-03 哈尔滨理工大学 Static pressure rotary worktable thermal characteristics research method under extreme operating condition
CN109376377A (en) * 2018-09-07 2019-02-22 东南大学 High-speed processing machine tool complete machine structure Thermodynamics modeling and Thermal design

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
HU, JINXIANG ET AL.: "Thermal Characteristics Analysis of the Work-piece Axle System in High Speed Internal Grinder", PRECISE MANUFACTURING & AUTOMATION, no. 1, 25 February 2007 (2007-02-25), pages 10 - 13, ISSN: 1009-962X *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112364444A (en) * 2020-09-23 2021-02-12 江苏赛洋机电科技有限公司 Numerical control machine tool temperature measuring point optimization method based on finite element model numerical simulation

Also Published As

Publication number Publication date
CN109376377A (en) 2019-02-22
CN109376377B (en) 2019-10-11
WO2020047880A1 (en) 2020-03-12

Similar Documents

Publication Publication Date Title
WO2020048110A1 (en) Thermodynamic modeling and thermal design method for overall structure of high-speed processing machine
Gorla et al. Theoretical and experimental analysis of a cycloidal speed reducer
CN112434452B (en) Method and system equipment for determining shaft diameter fit tolerance of high-speed motorized spindle bearing
CN102867088A (en) Thermal network modeling method applied to electric spindle steady temperature field
CN110705147B (en) Numerical control machine tool spindle thermal state characteristic comprehensive theoretical modeling and analysis method
Chen et al. Simulation on thermal characteristics of high-speed motorized spindle
CN102009370B (en) Feedback compensation method of high-speed cutting CNC (computer numerical control) machine tool
Dai et al. Application of an unstructured overset method for predicting the gear windage power losses
Fang et al. An improved thermal performance modeling for high-speed spindle of machine tool based on thermal contact resistance analysis
Zhao et al. Thermal analysis and experimental study on the spindle of the high-speed machining center
Deng et al. A predictive model for temperature rise of spindle–bearing integrated system
Fan et al. An optimization method for thermal behavior of high-speed spindle of gear form grinding machine
Fan et al. Thermoelectric-based cooling system for high-speed motorized spindle I: design and control mechanism
He et al. Investigation on heat transfer performance of heat pipe grinding wheel in dry grinding
Jiang et al. Research on thermal characteristics of spindle system with tilting pad bearing
Wu et al. Geometric design investigation of single screw compressor rotor grooves produced by cylindrical milling
Li et al. Analysis of fluid-solid-thermal coupling for ball screw in boring-milling machining center
Zhang et al. Hot oil carrying characteristic about hydrostatic bearing oil film of heavy vertical lathe in high speed
Zhang et al. Analysis for green grinding of Ti-6Al-4V titanium alloys with profile rotating heat pipe-grinding wheel
CN108959830B (en) High speed dynamic pressure lubricates accurate main shaft Thermodynamics modeling and Thermal design
Jan Analytical determination of rolling friction coefficient of angular bearings
Xue et al. An analytical prediction methodology of the correlation between transmission error and tooth pitch error for face-hobbed hypoid gears
Chen et al. A thermal analyzing model for high speed motorized spindles—thermal properties
Li et al. Analysis for Thermal Characteristics of High-Speed Angular Contact Ball Bearing under Different Preload
Chen et al. Model-based design of motorized spindles with different bearing configurations

Legal Events

Date Code Title Description
NENP Non-entry into the national phase

Ref country code: DE

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 24.01.2022)

122 Ep: pct application non-entry in european phase

Ref document number: 19857387

Country of ref document: EP

Kind code of ref document: A1