CN103246786A - 一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法 - Google Patents
一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法 Download PDFInfo
- Publication number
- CN103246786A CN103246786A CN2013101986475A CN201310198647A CN103246786A CN 103246786 A CN103246786 A CN 103246786A CN 2013101986475 A CN2013101986475 A CN 2013101986475A CN 201310198647 A CN201310198647 A CN 201310198647A CN 103246786 A CN103246786 A CN 103246786A
- Authority
- CN
- China
- Prior art keywords
- oil film
- model
- thrust bearing
- mechanical coupling
- worktable
- 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
Images
Abstract
一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法,本发明涉及润滑油膜模型的方法。本发明是要解决现有建立油膜模型的方法对工件效率与精度有很大影响的问题,而提供了一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法。步骤A:建立工作台、底座、扇形油腔和初始油膜的模型,并将工作台、底座与油腔装配;步骤B:前处理初始油膜的模型;步骤C:在workbench中建立流程图,对静压推力轴承进行热力耦合计算分析;步骤D:封闭工作台和油腔之间的间隙;步骤E:获取间隙润滑油膜的模型。本发明应用于润滑油膜模型领域。
Description
技术领域
本发明涉及润滑油膜模型的方法。
背景技术
在静压技术推广应用的同时国内外高校和科研单位进行了大量的实验及理论研究,包括静压轴承的轴承结构的改进、节流器的控制方式、设计计算方法的完善等。而在重型静压推力轴承性能研究方面,都侧重与采用定压供油方式的静压推力轴承的研究,且局限于其实验装置、节流器型式的选择、结构优化、不同类型节流器的比较、腔型结构优化、轴承特性和阻尼特性等方面的研究。而在对静压推力轴承承受力和热产生变形后的润滑油膜参数的获取研究几乎空白,在以往工程计算和数值模拟中一般忽略其变形对油膜态的影响,而实际上,这对工件效率与精度有很大影响。
静压推力轴承在实际工作中由于承载及长时间运行发热将会产生局部变形,而变形后间隙油膜的模型参数的获取通过理论计算是不能实现。
发明内容
本发明是要解决现有建立油膜模型的方法对工件效率与精度有很大影响的问题,而提供了一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法。
一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法按以下步骤实现:
步骤A、依据试验机床的尺寸结构与给定初始油膜厚度值,使用三维建模软件建立工作台、底座、扇形油腔和初始油膜的模型,并将工作台、底座与扇形油腔装配,导出zhuangpei.x_t和youmo.x_t文件;
步骤B、用GAMBITZ软件对初始油膜进行网格划分,得到网格,指定流体域边界条件并输出youmo.fluent5/6mesh文件,将youmo.fluent5/6mesh文件导入FLUENT中,输入或选择油膜入口温度Tin初值、油膜固定壁面温度初值、油膜入口流量Q、出口压力、旋转壁面角速度和固定壁面边界条件,然后进行模拟计算进行稳态分析最终得到youmo.cas文件;
步骤C、在ANSYS WOERKBECH软件中建立工作流程图,对静压推力轴承进行热力耦合计算分析:
步骤D、封闭工作台的油腔之间的间隙,在UG中导入bxh.stp文件,将油腔四周拉伸,然后导出fb.x_t文件;其中,所述油腔四周拉伸高度保证嵌入与油膜接触的扇形工作台面,厚度保证与油膜接触的扇形工作台面的扇形边界嵌入其中;
步骤E、得到间隙油膜的模型,实用ANSYS WOERKBECH软件重新建立工作流程图,添加Geometry几何模块将fb.x_t文件导入,经布尔运算和填充操作得到间隙油膜模型,即完成了一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法。
发明效果:
本发明方法利用计算机模拟现场实际工况下静压推力轴承承受热和力产生变形后的变形情况,对静压推力轴承热力耦合变形后的润滑油膜参数的获取研究,获取润滑油膜参数,揭示静压推力轴承的变形对其支撑特性、润滑特性的影响规律。这为静压推力轴承实现高效、高精奠定了技术基础,为静压推力轴承的油腔设计提供重要参考依据,最终实现静压推力轴承的结构优化设计和润滑系统设计提供重要的参考,避免出现静压轴承润滑失效。计算机数值模拟过程符合实际工况,获取变形后间隙油膜的模型参数,使设计人员可快速地获取静压推力轴承承受热和力产生变形后对其间隙油膜厚度影响的结果,为静压轴承设计更合理的结构提供了更有价值的理论依据数值模拟结果具有更重要的实用价值。在实际应用中,使用该方法验证的大型机床并对其进行结构优化以及对润滑油流量的合理控制,以保证其正常、安全运转并避免发生干摩擦和边界润滑,在保证一定精度的前提下,可以使机床转速普遍提高了30%~60%。
附图说明
图1是本发明流程图;
图2是本发明步骤C的具体流程图;
图3是具体实施方式一中最终获取的变形后间隙油膜模型图。
具体实施方式
具体实施方式一:本实施方式的一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法按以下步骤实现:
步骤A、依据试验机床的尺寸结构与给定初始油膜厚度值,使用三维建模软件建立工作台、底座、扇形油腔和初始油膜的模型,并将工作台、底座与油腔装配,导出zhuangpei.x_t和youmo.x_t文件(其中工作台和底座为1/12模型);
步骤B、用GAMBITZ软件对初始油膜进行网格划分,得到高质量的网格,指定流体域边界条件并输出youmo.fluent5/6mesh文件,将youmo.fluent5/6mesh文件导入FLUENT中输入或选择油膜入口温度Tin初值、油膜固定壁面温度初值、油膜入口流量Q、出口压力、旋转壁面角速度和固定壁面边界条件,然后进行模拟计算进行稳态分析最终得到youmo.cas文件;
步骤C、在ANSYS WOERKBECH软件中建立工作流程图,对静压推力轴承进行热力耦合计算分析:
步骤D、封闭工作台的油腔之间的间隙,在UG中导入bxh.stp文件,将油腔四周拉伸,高度保证嵌入与油膜接触的扇形工作台面,厚度保证与油膜接触的扇形工作台面的扇形边界嵌入其中,然后导出fb.x_t文件;
步骤E、得到间隙油膜的模型,实用ANSYS WOERKBECH软件重新建立工作流程图,添加Geometry几何模块将fb.x_t文件导入,经布尔运算和填充操作得到间隙油膜模型,即完成了一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法。
本实施方式效果:
本实施方式方法利用计算机模拟现场实际工况下静压推力轴承承受热和力产生变形后的变形情况,对静压推力轴承热力耦合变形后的润滑油膜参数的获取研究,获取润滑油膜参数,揭示静压推力轴承的变形对其支撑特性、润滑特性的影响规律。这为静压推力轴承实现高效、高精奠定了技术基础,为静压推力轴承的油腔设计提供重要参考依据,最终实现静压推力轴承的结构优化设计和润滑系统设计提供重要的参考,避免出现静压轴承润滑失效。计算机数值模拟过程符合实际工况,获取变形后间隙油膜的模型参数,使设计人员可快速地获取静压推力轴承承受热和力产生变形后对其间隙油膜厚度影响的结果,为静压轴承设计更合理的结构提供了更有价值的理论依据数值模拟结果具有更重要的实用价值。在实际应用中,使用该方法验证的大型机床并对其进行结构优化以及对润滑油流量的合理控制,以保证其正常、安全运转并避免发生干摩擦和边界润滑,在保证一定精度的前提下,可以使机床转速普遍提高了30%~60%。
具体实施方式二:本实施方式与具体实施方式一不同的是:步骤A中工作台与底座采用1/12模型分析。其它步骤及参数与具体实施方式一相同。
具体实施方式三:本实施方式与具体实施方式一或二不同的是:步骤C在ANSYSWOERKBECH软件中建立工作流程图,对静压推力轴承进行热力耦合计算分析具体为:
步骤C1、导入流体分析结果,添加FLUENT模块,将步骤B中youmo.cas文件导入到FLUENT中重新迭代计算;
步骤C2、导入装配模型,添加Geometry几何模块,将步骤A中zhuangpei.x_t文件导入;
步骤C3、热分析,插入Steady-State Thermal模块,设置工作台、底座和油腔的材料属性,由于只需要进行热分析和线性结构分析,故只需在材料属性窗口填上对应的密度、杨氏模量、泊松比和导热系数;将流体分析的结果导入到热分析中以温度边界条件作用在工作台和油腔相应接触面上,对装配体分别划分网格,设置工作台、油腔、底座的网格尺寸大小分别为40mm、20mm、45mm;设置边界条件对空气的热对流系数,添加热载荷工作台和对坐两侧对称面的完全绝热,然后进行求解计算;
步骤C4、线性静力结构分析,插入Static Structural模块,将流体分析结果导入到线性静力结构分析中以压力载荷作用在工作台和油腔相应接触面上,插入重力加速度载荷,设置底座与地面接触面的固定约束、工作台和对坐两侧对称面的无摩擦约束、工作台的圆柱面约束,求解计算;
步骤C5、更新变形后的实体模型,复制Static Structural模块,在Static Structural模块中右键点击geometry,再点update geometry from result file(要先打开beta选项),找结果文件,就在保存的文件夹***_files(或临时文件夹)里面,dp0-->SYS(-n)-->MECH-->file.rst,更新实体模型,为了方便计算由于底座对下面进行的油膜抽取没有影响故将其抑制(即不显示不参与运算),对工作台和油腔重新画网格网格尺寸大小分别为40mm和20mm;
步骤C6、工作台和油腔有限元模型转化为实体模型,插入FE模块,将步骤C5中重新画好的网格导入到FE中,调节容差15,工作台和底座的缝补公差分别为0.001和0.0051,得到实体模型;
步骤C7、导出工作台和油腔实体,插入Geometry几何模块,在Geometry中将实体导出bxh.stp以便在三维建模软件UG中进行操作。其它步骤及参数与具体实施方式一或二相同。
具体实施方式四:本实施方式与具体实施方式一至三之一不同的是:步骤C3中所述工作台和底座为灰口铸铁,油腔铝,在材料属性窗口填上对应的密度、杨氏模量、泊松比和导热系数。其它步骤及参数与具体实施方式一至三之一相同。
Claims (4)
1.一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法,其特征在于建立静压推力轴承热力耦合变形后润滑油膜模型的方法按以下步骤实现:
步骤A、依据试验机床的尺寸结构与给定初始油膜厚度值,使用三维建模软件建立工作台、底座、扇形油腔和初始油膜的模型,并将工作台、底座与扇形油腔装配,导出zhuangpei.x_t和youmo.x_t文件;
步骤B、用GAMBITZ软件对初始油膜进行网格划分,得到网格,指定流体域边界条件并输出youmo.fluent5/6mesh文件,将youmo.fluent5/6mesh文件导入FLUENT中,输入或选择油膜入口温度Tin初值、油膜固定壁面温度初值、油膜入口流量Q、出口压力、旋转壁面角速度和固定壁面边界条件,然后进行模拟计算进行稳态分析最终得到youmo.cas文件;
步骤C、在ANSYS WOERKBECH软件中建立工作流程图,对静压推力轴承进行热力耦合计算分析:
步骤D、封闭工作台的油腔之间的间隙,在UG中导入bxh.stp文件,将油腔四周拉伸,然后导出fb.x_t文件;其中,所述油腔四周拉伸高度保证嵌入与油膜接触的扇形工作台面,厚度保证与油膜接触的扇形工作台面的扇形边界嵌入其中;
步骤E、得到间隙油膜的模型,实用ANSYS WOERKBECH软件重新建立工作流程图,添加Geometry几何模块将fb.x_t文件导入,经布尔运算和填充操作得到间隙油膜模型,即完成了一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法。
2.根据权利要求1所述的一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法,其特征在于所述步骤A中工作台与底座采用1/12模型分析。
3.根据权利要求1所述的一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法,其特征在于所述步骤C在ANSYS WOERKBECH软件中建立工作流程图,对静压推力轴承进行热力耦合计算分析具体为:
步骤C1、导入流体分析结果,添加FLUENT模块,将步骤B中youmo.cas文件导入到FLUENT中重新迭代计算;
步骤C2、导入装配模型,添加Geometry几何模块,将步骤A中zhuangpei.x_t文件导入;
步骤C3、热分析,插入Steady-State Thermal模块,设置工作台、底座和扇形油腔的材料属性,将流体分析结果导入到热分析中,对装配体分别划分网格,设置热载荷和热边界条件,然后进行求解计算;
步骤C4、线性静力结构分析,插入Static Structural模块,将流体分析结果导入到线性静力结构分析中,插入载荷和约束,求解计算;
步骤C5、更新变形后的实体模型,对工作台和油腔重新画网格;
步骤C6、工作台和扇形油腔有限元模型转化为实体模型,插入FE模块,将步骤C5中重新画好的网格导入到FE中,调节容差和缝补公差得到实体模型;
步骤C7、导出工作台和油腔实体,插入Geometry几何模块,在Geometry中将实体导出bxh.stp以便在三维建模软件UG中进行操作。
4.根据权利要求3所述的一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法,其特征在于所述步骤C3中所述工作台和底座为灰口铸铁,油腔铝,在材料属性窗口填上对应的密度、杨氏模量、泊松比和导热系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310198647.5A CN103246786B (zh) | 2013-05-24 | 2013-05-24 | 一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310198647.5A CN103246786B (zh) | 2013-05-24 | 2013-05-24 | 一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103246786A true CN103246786A (zh) | 2013-08-14 |
CN103246786B CN103246786B (zh) | 2015-12-23 |
Family
ID=48926304
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310198647.5A Expired - Fee Related CN103246786B (zh) | 2013-05-24 | 2013-05-24 | 一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103246786B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103593537A (zh) * | 2013-11-27 | 2014-02-19 | 哈尔滨理工大学 | 一种建立静压中心架热力耦合变形后润滑油膜模型的方法 |
CN103632016A (zh) * | 2013-12-20 | 2014-03-12 | 哈尔滨理工大学 | 一种基于动网格技术的静压轴承压力场模拟方法 |
CN104537189A (zh) * | 2015-01-21 | 2015-04-22 | 北京工业大学 | 一种静压转台运动误差建模及计算方法 |
CN105069264A (zh) * | 2015-09-08 | 2015-11-18 | 哈尔滨理工大学 | 双矩形腔静压推力轴承旋转速度与承载合理匹配方法 |
CN108595150A (zh) * | 2018-04-27 | 2018-09-28 | 北京航空航天大学 | 人工智能程序员书写数字卫星力热耦合源程序的方法 |
CN109002569A (zh) * | 2018-01-12 | 2018-12-14 | 哈尔滨理工大学 | 一种建立静压油膜运动层边界条件的方法 |
CN110287576A (zh) * | 2019-06-20 | 2019-09-27 | 哈尔滨理工大学 | 一种基于Matlab的液体静压推力轴承润滑油膜三维温度场显示方法 |
CN113027921A (zh) * | 2021-02-09 | 2021-06-25 | 太原重工股份有限公司 | 获取静动压轴承油膜压力分布的方法和装置 |
CN113158369A (zh) * | 2021-04-19 | 2021-07-23 | 哈尔滨理工大学 | 静压推力轴承油垫封油边的油膜流量模拟监测方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109885881B (zh) * | 2019-01-18 | 2020-10-09 | 北京交通大学 | 一种造成轴承电腐蚀的临界轴电流的确定方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005106266A1 (de) * | 2004-05-04 | 2005-11-10 | Sms Demag Ag | Ölfilmlager für walzenzapfen mit hydrostatischer unterstützung |
CN102141084A (zh) * | 2011-03-30 | 2011-08-03 | 哈尔滨理工大学 | 基于膜厚可变的静压推力轴承间隙油膜温度与厚度关系数值模拟方法 |
CN102214260A (zh) * | 2011-07-01 | 2011-10-12 | 太原科技大学 | 用于油膜轴承润滑油膜的参数化计算分析系统及操作方法 |
-
2013
- 2013-05-24 CN CN201310198647.5A patent/CN103246786B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005106266A1 (de) * | 2004-05-04 | 2005-11-10 | Sms Demag Ag | Ölfilmlager für walzenzapfen mit hydrostatischer unterstützung |
CN102141084A (zh) * | 2011-03-30 | 2011-08-03 | 哈尔滨理工大学 | 基于膜厚可变的静压推力轴承间隙油膜温度与厚度关系数值模拟方法 |
CN102214260A (zh) * | 2011-07-01 | 2011-10-12 | 太原科技大学 | 用于油膜轴承润滑油膜的参数化计算分析系统及操作方法 |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103593537A (zh) * | 2013-11-27 | 2014-02-19 | 哈尔滨理工大学 | 一种建立静压中心架热力耦合变形后润滑油膜模型的方法 |
CN103632016A (zh) * | 2013-12-20 | 2014-03-12 | 哈尔滨理工大学 | 一种基于动网格技术的静压轴承压力场模拟方法 |
CN104537189A (zh) * | 2015-01-21 | 2015-04-22 | 北京工业大学 | 一种静压转台运动误差建模及计算方法 |
CN104537189B (zh) * | 2015-01-21 | 2018-02-02 | 北京工业大学 | 一种静压转台运动误差建模及计算方法 |
CN105069264A (zh) * | 2015-09-08 | 2015-11-18 | 哈尔滨理工大学 | 双矩形腔静压推力轴承旋转速度与承载合理匹配方法 |
CN105069264B (zh) * | 2015-09-08 | 2018-01-16 | 哈尔滨理工大学 | 双矩形腔静压推力轴承旋转速度与承载合理匹配方法 |
CN109002569A (zh) * | 2018-01-12 | 2018-12-14 | 哈尔滨理工大学 | 一种建立静压油膜运动层边界条件的方法 |
CN108595150A (zh) * | 2018-04-27 | 2018-09-28 | 北京航空航天大学 | 人工智能程序员书写数字卫星力热耦合源程序的方法 |
CN108595150B (zh) * | 2018-04-27 | 2020-09-25 | 北京航空航天大学 | 人工智能程序员书写数字卫星力热耦合源程序的方法 |
CN110287576A (zh) * | 2019-06-20 | 2019-09-27 | 哈尔滨理工大学 | 一种基于Matlab的液体静压推力轴承润滑油膜三维温度场显示方法 |
CN113027921A (zh) * | 2021-02-09 | 2021-06-25 | 太原重工股份有限公司 | 获取静动压轴承油膜压力分布的方法和装置 |
CN113027921B (zh) * | 2021-02-09 | 2022-10-11 | 太原重工股份有限公司 | 获取静动压轴承油膜压力分布的方法和装置 |
CN113158369A (zh) * | 2021-04-19 | 2021-07-23 | 哈尔滨理工大学 | 静压推力轴承油垫封油边的油膜流量模拟监测方法 |
CN113158369B (zh) * | 2021-04-19 | 2023-11-28 | 哈尔滨理工大学 | 静压推力轴承油垫封油边的油膜流量模拟监测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103246786B (zh) | 2015-12-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103246786A (zh) | 一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法 | |
CN103020358B (zh) | 针对机械装备的可适应性动态设计平台构建方法 | |
Parab et al. | Structural and thermal analysis of brake disc | |
Garaizar et al. | Integration of thermal effects into tolerancing using skin model shapes | |
CN102279126A (zh) | 测试与cae仿真相结合确定材料性能参数的方法 | |
CN103593537A (zh) | 一种建立静压中心架热力耦合变形后润滑油膜模型的方法 | |
Bressan et al. | Analysis of aluminium hot extrusion by finite volume method | |
CN105095536B (zh) | 一种考虑表面形貌的静压油垫流场特性数值仿真方法 | |
Li | Analyzing the flow field in the oil chamber of a hydrostatic guide rail used for ultra-precision machining: numerical simulation and performance optimization | |
CN103632016A (zh) | 一种基于动网格技术的静压轴承压力场模拟方法 | |
Meng et al. | Research on the virtual prototyping technology of the electric-hydraulic-mechanical system based on Simcape | |
Lu et al. | Motion simulation and finite element analysis of the manipulator based on SolidWorks | |
Tong et al. | The dynamic finite element analysis of shearer’s running gear based on LS-DYNA | |
Meng et al. | Effect of sink flow on dual-valve electro-hydraulic excitation system | |
Teramoto et al. | A framework of compositional machining simulation for versatile machining simulation | |
Cheng | Finite element analysis and structural optimization of the box on the ANSYS workbench | |
Wang et al. | An accelerated medial object transformation for whole engine optimisation | |
Li et al. | Dynamic Characteristics Analysis of a Micro Grinding Machine Tool | |
Zhou et al. | Lightweight Design of Valve Body Structure Based on Numerical Simulation | |
Tong et al. | Research on the Application of Structural Topology Optimisation in the High-Precision Design of a Press Machine Frame | |
Ha | Inventor in the Virtual Product Design | |
Li et al. | Definition of computational domain for flow simulation of wind turbine blades | |
Lan-tao | Inventor in the Virtual Product Design | |
Cheng et al. | The parametric design and automatic assembly of hydrostatic rotary table based on pro/engineer | |
Yong-xiang et al. | Parametric Design of Complex Elastic Components |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151223 Termination date: 20160524 |
|
CF01 | Termination of patent right due to non-payment of annual fee |