CN105488246B - 一种静压支承转台速度对热油携带影响的数值计算方法 - Google Patents

一种静压支承转台速度对热油携带影响的数值计算方法 Download PDF

Info

Publication number
CN105488246B
CN105488246B CN201510773768.7A CN201510773768A CN105488246B CN 105488246 B CN105488246 B CN 105488246B CN 201510773768 A CN201510773768 A CN 201510773768A CN 105488246 B CN105488246 B CN 105488246B
Authority
CN
China
Prior art keywords
deep fat
hydrostatic support
carries
heavy wool
turntable
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
CN201510773768.7A
Other languages
English (en)
Other versions
CN105488246A (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.)
Harbin University of Science and Technology
Original Assignee
Harbin University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Priority to CN201510773768.7A priority Critical patent/CN105488246B/zh
Publication of CN105488246A publication Critical patent/CN105488246A/zh
Application granted granted Critical
Publication of CN105488246B publication Critical patent/CN105488246B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Bearings And Hydrostatic Bearings (AREA)

Abstract

多油垫式静压支承在进行切削加工过程中，随着转台速度的增加，支承油膜的发热量加剧，主要来源于热油携带的影响，因此，本发明提出一种静压支承转台速度对热油携带影响的数值计算方法，目的是提高静压支承的加工效率及加工精度。步骤A、建立多油垫静压支承单个油垫模型及热油携带面。步骤B、将本发明涉及的转速工况参数代入推导出的热油携带方程得出公式计算值。步骤C、确定多油垫式静压支承油膜数值模拟条件。步骤D、进行不同转速下多油垫静压支承热油携带数值计算。步骤E、将计算收敛后的数据结果进行后处理，并与公式计算值对比得出对比曲线图验证此发明方法的有效性。本发明应用于一种静压支承转台速度对热油携带影响的数值计算方法。

Description

一种静压支承转台速度对热油携带影响的数值计算方法
技术领域
本发明涉及了一种静压支承转台速度对热油携带影响的数值计算方法。
背景技术
多油垫式静压支承圆导轨由于其结构尺寸很大，多用于重型数控装备制造业中，与普通静压支承比较，在具有普通静压支承的一系列特点外，还具有承载能力强，线速度高，发热量大等特点。随着现代机械工业中的重型装备逐步地走向重载荷、高速度、高效率、高精度、更智能化的发展趋势，液体静压滑动支承定将成为一项有生命力的技术，也成为了学者们研究的热点问题。本发明所涉及的多油垫式静压支承在进行切削加工过程中，随着转台速度的增加，支承油膜的发热量加剧，主要来源于热油携带的影响，因此，本发明提出一种静压支承转台速度对热油携带影响的数值计算方法，获取不同转台速度下的静压支承热油携带情况，目的是为了揭示出静压支承的转台速度对热油携带的影响规律，提高静压支承的加工效率及加工精度。
发明内容
本发明是要解决现有多油垫静压支承加工过程中由于油膜温升过高而限制了转台速度提升的问题，为此而提供了一种静压支承转台速度对热油携带影响的数值计算方法。
一种静压支承转台速度对热油携带影响的数值计算方法按以下步骤实现：
步骤A、建立多油垫静压支承单个油垫模型及热油携带面如下：
本发明涉及到的多油垫静压支承用在十米双柱立式数控车床上，该车床的最大加工工件直径为10米，转台直径为7.5米，转台重量62吨。采用三维建模软件建立圆导轨中的单个油垫模型如图1所示，并在油膜单个油垫上创建了截面3，命名为face3，此热油携带面位于单个油垫上右侧的封油边与回油槽交接位置处；在油垫上封油边周向左侧又创建了截面2，命名为face2。
步骤B、利用前期推导得出的热油携带计算公式，将本发明涉及的多油垫静压支承结构参数及转速工况参数代入，列出了热油携带情况的理论计算值，包括face3和face2的面流量值及热油携带因子K的值。
步骤C、确定多油垫式静压支承油膜数值模拟条件：
根据静压支承实际运行情况，分别在转速为0.5r/min、1r/min、2r/min、4r/min、8r/min、16r/min、32r/min时，在承载40t工况下进行多油垫静压支承热油流量数值计算。
步骤D、进行不同转速下多油垫静压支承热油携带数值计算。
步骤E、将计算收敛后的数据结果进行后处理，得出了不同转速下的face3和face2的热油流量变化值数值模拟结果曲线，并与推导公式的计算值进行对比得出对比曲线图，通过对比验证此发明的多油垫静压支承热油携带数值计算方法的有效性。
发明效果：
本发明通过提出一种静压支承转台速度对热油携带影响的数值计算方法，获取不同转台速度下的静压支承热油携带情况，并判定出不同转速下的静压支承热油携带因子变化情况，且热油携带系数K模拟值较理论公式计算值偏差在5.15%范围内，误差主要来源于理论计算时对推导公式中模型尺寸的部分简化。本发明的计算方法可以推广到不同结构尺寸的系列静压支承研究中，不仅提高了计算结果的准确度，且计算周期也可提高两倍以上，大大改善了计算效率，实现了多油垫静压支承运行过程中发热问题的有效预测，提高了静压支承的加工效率及加工精度。
附图说明
图1是本发明单个油垫模型及封face3与face2位置图；
图2到图15是具体实施方式中不同转速工况下的face3和face2的流量模拟值，其中：图2为转速0.5r/min工况截面face2流量，图3为转速0.5r/min工况截面face3流量，图4为转速1r/min工况截面face2流量，图5为转速1r/min工况截面face3流量，图6 为转速2r/min工况截面face2流量，图7为转速2r/min工况截面face3流量，图8为转速4r/min工况截面face2流量，图9为转速4r/min工况截面face3流量，图10为转速8r/min工况截面face2流量，图11为转速8r/min工况截面face3流量，图12为转速16r/min工况截面face2流量，图13为转速16r/min工况截面face3流量，图14为转速32r/min工况截面face2流量，图15为转速32r/min工况截面face3流量；
图16是具体实施方式中热油携带因子随转台速度变化曲线。
具体实施方式
一种静压支承转台速度对热油携带影响的数值计算方法按以下步骤实现：
步骤A、建立多油垫静压支承单个油垫模型及热油携带面如下：
本发明涉及到的多油垫静压支承用在十米双柱立式数控车床上，该车床的最大加工工件直径为10米，转台直径为7.5米，转台重量62吨。采用三维建模软件建立圆导轨中的单个油垫模型如图1所示，并在油膜单个油垫上创建了截面3，命名为face3，此热油携带面位于单个油垫上右侧的封油边与回油槽交接位置处；在油垫上封油边周向左侧又创建了截面2，命名为face2。
步骤B、利用前期推导得出的热油携带计算公式，将本发明涉及的多油垫静压支承结构参数及转速工况参数代入，列出了表B中的热油携带情况的理论计算值，包括face3和face2的面流量值及热油携带因子K的值。
表B 40t时不同转台速度工况下的Q3、Q2、K理论公式值
R(r/min) 0.5 1 2 4 8 16 32
Q3 (10-6m3/s) 8.25 4.85 1.95 15.53 42.71 97.07 205.7
Q2 (10-6m3/s) 15.05 18.45 25.25 38.83 66.01 120.3 229.0
K 0 0 7.72% 40.0% 64.7% 80.6% 89.8%
步骤C、确定多油垫式静压支承油膜数值模拟条件：
根据静压支承实际运行情况，分别在转速为0.5r/min、1r/min、2r/min、4r/min、8r/min、16r/min、32r/min时，在承载40t工况下进行多油垫静压支承热油流量数值计算。
步骤D、进行不同转速下多油垫静压支承热油携带数值计算：
步骤D1、打开流体力学计算软件，通过File中的Import导入网格模型文件jyzc-youmo. mesh，并通过Grid进行网格质量检查；
步骤D2、通过Define进行求解器Solver的参数设置，采用三维定常流体流动设置；
步骤D3、通过Define进行材料属性Materials设置，主要给出所研究静压支承工作介质的粘度和密度参数等；
步骤D4、进行Define-Boundary Conditions边界条件设定，分别给出不同载荷下的入口速度，设定静压支承工作台的旋转速度，出口处选择压力出口边界；
步骤D5、进行监视器平面设置，将face3和face2作为要监视的平面，选择质量流Mass-flow参数监视，获取质量流曲线；
步骤D6、初始化流场，进行迭代计算，迭代800步数后提示流场数值计算收敛，获取计算结果。
步骤E、将计算收敛后的数据结果进行后处理，得出了如图2所示不同转速下的face3和face2的热油流量变化值数值模拟结果曲线，并与推导公式的计算值进行对比，图3为对比曲线图，通过对比验证此发明的多油垫静压支承热油携带数值计算方法的有效性。
本实施方式效果：
通过本发明的数值计算结果分析可知，多油垫静压支承转台转速小于或等于1r/min时，不会发生热油携带现象；当转台转速大于1r/min，静压支承有热油携带现象发生，且随着转速的增加，热油携带系数呈逐渐增大的趋势。且热油携带系数K模拟值较理论公式计算值偏差在5.15%范围内，误差主要来源于理论计算时对推导公式中模型尺寸的部分简化。本发明的一种静压支承转台速度对热油携带影响的数值计算方法可以推广到不同结构尺寸的系列静压支承研究中，不仅提高了计算结果的准确度，且计算周期也可提高两倍以上，大大改善了计算效率，实现了多油垫静压支承运行过程中发热问题的有效预测，可避免静压支承的烧瓦现象发生。

Claims (2)

1.一种静压支承转台速度对热油携带影响的数值计算方法， 其特征在于提出一种静压支承转台速度对热油携带影响的数值计算方法，获取不同转台速度下的静压支承热油携带情况；
一种静压支承转台速度对热油携带影响的数值计算方法按以下步骤实现：
步骤A、建立多油垫静压支承单个油垫模型及热油携带面如下：
多油垫静压支承用在十米双柱立式数控车床上，该车床的最大加工工件直径为10米，转台直径为7.5米，转台重量62吨， 采用三维建模软件建立圆导轨中的单个油垫模型，并在单个油垫油膜模型上创建了两个热油携带面；
步骤B、利用前期推导得出的热油携带计算公式，将涉及的多油垫静压支承结构参数及转速工况参数代入，列出了热油携带情况的理论计算值；
步骤C、确定多油垫式静压支承油膜数值模拟条件：
根据静压支承实际运行情况，分别在转速为0.5r/min、1r/min、2r/min、4r/min、8r/min、16r/min、32r/min时，在承载40t工况下进行多油垫静压支承热油流量数值计算；
步骤D、进行不同转速下多油垫静压支承热油携带数值计算；
步骤E、将计算收敛后的数据结果进行后处理，得出了不同转速下的热油携带面的热油流量变化值数值模拟结果曲线，并与推导公式的计算值进行对比得出对比曲线图，通过对比验证多油垫静压支承热油携带数值计算方法的有效性。
2.根据权利要求1所述的一种静压支承转台速度对热油携带影响的数值计算方法，其特征在于当多油垫静压支承转台转速小于或等于1r/min时，不会发生热油携带现象，当转台转速大于1r/min，静压支承有热油携带现象发生，且随着转速的增加，热油携带系数呈逐渐增大的趋势； 且热油携带系数K模拟值较理论公式计算值偏差在5.15%范围内，误差主要来源于理论计算时对推导公式中模型尺寸的部分简化。
CN201510773768.7A 2015-11-13 2015-11-13 一种静压支承转台速度对热油携带影响的数值计算方法 Expired - Fee Related CN105488246B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510773768.7A CN105488246B (zh) 2015-11-13 2015-11-13 一种静压支承转台速度对热油携带影响的数值计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510773768.7A CN105488246B (zh) 2015-11-13 2015-11-13 一种静压支承转台速度对热油携带影响的数值计算方法

Publications (2)

Publication Number Publication Date
CN105488246A CN105488246A (zh) 2016-04-13
CN105488246B true CN105488246B (zh) 2018-08-10

Family

ID=55675220

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510773768.7A Expired - Fee Related CN105488246B (zh) 2015-11-13 2015-11-13 一种静压支承转台速度对热油携带影响的数值计算方法

Country Status (1)

Country Link
CN (1) CN105488246B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2195519A1 (fr) * 2007-10-12 2010-06-16 Renault S.A.S. Estimation de parametres d'etat d'un moteur par mesure de la pression interne d'un cylindre
CN102980755A (zh) * 2012-11-16 2013-03-20 北京工业大学 一种定量式静压转台动静态特性实验装置
CN103823940A (zh) * 2014-03-07 2014-05-28 北京工业大学 一种基于整体动力学模型的静压转台动态响应计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2195519A1 (fr) * 2007-10-12 2010-06-16 Renault S.A.S. Estimation de parametres d'etat d'un moteur par mesure de la pression interne d'un cylindre
CN102980755A (zh) * 2012-11-16 2013-03-20 北京工业大学 一种定量式静压转台动静态特性实验装置
CN103823940A (zh) * 2014-03-07 2014-05-28 北京工业大学 一种基于整体动力学模型的静压转台动态响应计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Numerical Calculation and Analysis of Dynamic Pressure Effect for Heavy Hydrostatic Bearing Disk with Multiple Oil Pads;Shao J 等;《Advanced Science Letters》;20111231;第540-545页 *
重型数控机床静压转台承载力分析;赵永胜 等;《高技术通讯》;20130531;第23卷(第5期);第2561-2565页 *

Also Published As

Publication number Publication date
CN105488246A (zh) 2016-04-13

Similar Documents

Publication Publication Date Title
CN104133417A (zh) 叶片式流体机械数控加工切削力的快速预测方法
CN103186698B (zh) 一种重型机床静压转台动、静态性能仿真优化方法
CN102063540A (zh) 一种机床床身结构优化设计方法
Williams et al. Experiments and computations on large tip clearance effects in a linear cascade
CN104537189B (zh) 一种静压转台运动误差建模及计算方法
CN105631122A (zh) 一种机床大件的热变形仿真分析与建模方法
CN102063549A (zh) 一种机床主轴的流固热耦合数值仿真方法
CN107505842A (zh) 一种数控机床广义空间切削稳定性预测与优化方法
CN106649986A (zh) 一种基于procast仿真平台对铜管水平连铸参数最优化匹配方法
CN111123833B (zh) 一种机床的热误差补偿计算方法、装置、设备和存储介质
CN104091021B (zh) 一种计算倾斜时油垫承载力的方法
CN105488246B (zh) 一种静压支承转台速度对热油携带影响的数值计算方法
CN105666244A (zh) 一种数控镗床热效应下镗杆热伸长误差温度测点约简的方法
CN114254460A (zh) 基于图卷积神经网络的透平机械气动鲁棒性优化方法
CN105095583B (zh) 一种微尺度下的静压主轴模态分析方法
CN109117522B (zh) 一种基于matlab-ansys软件的静压溜板流固耦合工作性能的计算方法
CN105095536B (zh) 一种考虑表面形貌的静压油垫流场特性数值仿真方法
CN102880766A (zh) 基于参数化设计的液体静压主轴制作方法
CN111274744B (zh) 一种粒子法模拟mcci中气泡夹带现象的边界处理方法
CN104598666A (zh) 大尺度静压支承临界载荷参数获取方法
CN104091004B (zh) 一种基于改进的粒子群算法的定量式静压转台优化设计方法
CN115994422A (zh) 基于拉削性能的拉刀刃形空间曲线参数化设计方法
US20120087786A1 (en) Stator Structure of Turbo Molecular Pump and Method for Manufacturing the Same
CN110245408A (zh) 一种汽轮机单圆弧压力面叶型设计方法
CN105317840B (zh) 一种重型静压轴承热油携带判定方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into 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: 20180810

Termination date: 20191113

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