专利名称：一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法
技术领域：
本发明涉及润滑油膜模型的方法。
背景技术：
在静压技术推广应用的同时国内外高校和科研单位进行了大量的实验及理论研究，包括静压轴承的轴承结构的改进、节流器的控制方式、设计计算方法的完善等。而在重型静压推力轴承性能研究方面，都侧重与采用定压供油方式的静压推力轴承的研究，且局限于其实验装置、节流器型式 的选择、结构优化、不同类型节流器的比较、腔型结构优化、轴承特性和阻尼特性等方面的研究。而在对静压推力轴承承受力和热产生变形后的润滑油膜参数的获取研究几乎空白，在以往工程计算和数值模拟中一般忽略其变形对油膜态的影响，而实际上，这对工件效率与精度有很大影响。静压推力轴承在实际工作中由于承载及长时间运行发热将会产生局部变形，而变形后间隙油膜的模型参数的获取通过理论计算是不能实现。
发明内容
本发明是要解决现有建立油膜模型的方法对工件效率与精度有很大影响的问题，而提供了一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法。一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法按以下步骤实现:步骤A、依据试验机床的尺寸结构与给定初始油膜厚度值，使用三维建模软件建立工作台、底座、扇形油腔和初始油膜的模型，并将工作台、底座与扇形油腔装配，导出zhuangpe1.x_t 和 youm0.x_t 文件；步骤B、用GAMBITZ软件对初始油膜进行网格划分，得到网格，指定流体域边界条件并输出 youm0.fluent5/6mesh 文件，将 youm0.fluent5/6mesh 文件导入 FLUENT 中，输入或选择油膜入口温度Tin初值、油膜固定壁面温度初值、油膜入口流量Q、出口压力、旋转壁面角速度和固定壁面边界条件，然后进行模拟计算进行稳态分析最终得到youm0.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、依据试验机床的尺寸结构与给定初始油膜厚度值，使用三维建模软件建立工作台、底座、扇形油腔和初始油膜的模型，并将工作台、底座与油腔装配，导出zhuangpe1.x_t和youm0.x_t文件(其中工作台和底座为1/12模型);·
步骤B、用GAMBITZ软件对初始油膜进行网格划分，得到高质量的网格，指定流体域边界条件并输出youm0.fluent5/6mesh文件，将youm0.fluent5/6mesh文件导入FLUENT中输入或选择油膜入口温度Tin初值、油膜固定壁面温度初值、油膜入口流量Q、出口压力、旋转壁面角速度和固定壁面边界条件，然后进行模拟计算进行稳态分析最终得到youm0.cas文件；步骤C、在ANSYS WOERKBECH软件中建立工作流程图，对静压推力轴承进行热力耦合计算分析:步骤D、封闭工作台的油腔之间的间隙，在UG中导入bxh.stp文件，将油腔四周拉伸，高度保证嵌入与油膜接触的扇形工作台面，厚度保证与油膜接触的扇形工作台面的扇形边界嵌入其中，然后导出fb.x_t文件；步骤E、得到间隙油膜的模型，实用ANSYS WOERKBECH软件重新建立工作流程图，添加Geometry几何模块将fb.x_t文件导入，经布尔运算和填充操作得到间隙油膜模型，即完成了一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法。本实施方式效果:本实施方式方法利用计算机模拟现场实际工况下静压推力轴承承受热和力产生变形后的变形情况，对静压推力轴承热力耦合变形后的润滑油膜参数的获取研究，获取润滑油膜参数，揭示静压推力轴承的变形对其支撑特性、润滑特性的影响规律。这为静压推力轴承实现高效、高精奠定了技术基础，为静压推力轴承的油腔设计提供重要参考依据，最终实现静压推力轴承的结构优化设计和润滑系统设计提供重要的参考，避免出现静压轴承润滑失效。计算机数值模拟过程符合实际工况，获取变形后间隙油膜的模型参数，使设计人员可快速地获取静压推力轴承承受热和力产生变形后对其间隙油膜厚度影响的结果，为静压轴承设计更合理的结构提供了更有价值的理论依据数值模拟结果具有更重要的实用价值。在实际应用中，使用该方法验证的大型机床并对其进行结构优化以及对润滑油流量的合理控制，以保证其正常、安全运转并避免发生干摩擦和边界润滑，在保证一定精度的前提下，可以使机床转速普遍提高了 30% 60%。
具体实施方式
二:本实施方式与具体实施方式
一不同的是:步骤A中工作台与底座采用1/12模型分析。其它步骤及参数与具体实施方式
一相同。
具体实施方式
三:本实施方式与具体实施方式
一或二不同的是:步骤C在ANSYSffOERKBECH软件中建立工作流程图，对静压推力轴承进行热力耦合计算分析具体为:步骤Cl、导入流体分析结果,添加FLUENT模块,将步骤B中youm0.cas文件导入到FLUENT中重新迭代计算；步骤C2、导入装配模型,添加Geometry几何模块,将步骤A中zhuangpe1.x_t文件导入；步骤C 3、热分析,插入Steady-State Thermal模块,设置工作台、底座和油腔的材料属性，由于只需要进行热分析和线性结构分析，故只需在材料属性窗口填上对应的密度、杨氏模量、泊松比和导热系数；将流体分析的结果导入到热分析中以温度边界条件作用在工作台和油腔相应接触面上，对装配体分别划分网格，设置工作台、油腔、底座的网格尺寸大小分别为40mm、20mm、45mm ;设置边界条件对空气的热对流系数，添加热载荷工作台和对坐两侧对称面的完全绝热，然后进行求解计算；步骤C4、线性静力结构分析,插入Static Structural模块,将流体分析结果导入到线性静力结构分析中以压力载荷作用在工作台和油腔相应接触面上，插入重力加速度载荷，设置底座与地面接触面的固定约束、工作台和对坐两侧对称面的无摩擦约束、工作台的圆柱面约束，求解计算；步骤C5、更新变形后的实体模型，复制Static Structural模块，在StaticStructural 模块中右键点击 geometry,再点 update geometry from result file (要先打开beta选项),找结果文件,就在保存的文件夹***_files (或临时文件夹)里面,dpO— >SYS(-n)- > MECH- > file, rst，更新实体模型，为了方便计算由于底座对下面进行的油膜抽取没有影响故将其抑制(即不显示不参与运算)，对工作台和油腔重新画网格网格尺寸大小分别为40mm和20mm ；步骤C6、工作台和油腔有限元模型转化为实体模型，插入FE模块，将步骤C5中重新画好的网格导入到FE中，调节容差15，工作台和底座的缝补公差分别为0.001和
0.0051，得到实体模型；步骤C7、导出工作台和油腔实体,插入Geometry几何模块,在Geometry中将实体导出bxh.Stp以便在三维建模软件UG中进行操作。其它步骤及参数与具体实施方式
一或二相同。
具体实施方式
四:本实施方式与具体实施方式
一至三之一不同的是:步骤C3中所述工作台和底座为灰口铸铁，油腔铝，在材料属性窗口填上对应的密度、杨氏模量、泊松比和导热系数。其 它步骤及参数与具体实施方式
一至三之一相同。
权利要求
1.一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法，其特征在于建立静压推力轴承热力耦合变形后润滑油膜模型的方法按以下步骤实现: 步骤A、依据试验机床的尺寸结构与给定初始油膜厚度值，使用三维建模软件建立工作台、底座、扇形油腔和初始油膜的模型，并将工作台、底座与扇形油腔装配，导出zhuangpe1.x_t 和 youm0.x_t 文件； 步骤B、用GAMBITZ软件对初始油膜进行网格划分，得到网格，指定流体域边界条件并输出 youm0.fluent5/6mesh 文件，将 youm0.fluent5/6mesh 文件导入 FLUENT 中，输入或选择油膜入口温度Tin初值、油膜固定壁面温度初值、油膜入口流量Q、出口压力、旋转壁面角速度和固定壁面边界条件，然后进行模拟计算进行稳态分析最终得到youm0.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软件中建立工作流程图，对静压推力轴承进行热力耦合计算分析具体为: 步骤Cl、导入流体分析结果,添加FLUENT模块,将步骤B中youm0.cas文件导入到FLUENT中重新迭代计算； 步骤C2、导入装配模型,添加Geometry几何模块,将步骤A中zhuangpe1.x_t文件导A ； 步骤C3、热分析,插入Steady-State Thermal模块,设置工作台、底座和扇形油腔的材料属性，将流体分析结果导入到热分析中，对装配体分别划分网格，设置热载荷和热边界条件，然后进行求解计算； 步骤C4、线性静力结构分析,插入Static Structural模块,将流体分析结果导入到线性静力结构分析中，插入载荷和约束，求解计算； 步骤C5、更新变形后的实体模型，对工作台和油腔重新画网格； 步骤C6、工作台和扇形油腔有限元模型转化为实体模型，插入FE模块，将步骤C5中重新画好的网格导入到FE中，调节容差和缝补公差得到实体模型； 步骤C7、导出工作台和油腔实体,插入Geometry几何模块,在Geometry中将实体导出bxh.stp以便在三维建模软件UG中进行操作。
4.根据权利要求3所述的一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法，其特征在于所述步骤C3中所述工作台和底座为灰口铸铁，油腔铝，在材料属性窗口填上对应的密度、杨氏模量、泊松比和导热系数。
全文摘要
一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法，本发明涉及润滑油膜模型的方法。本发明是要解决现有建立油膜模型的方法对工件效率与精度有很大影响的问题，而提供了一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法。步骤A建立工作台、底座、扇形油腔和初始油膜的模型，并将工作台、底座与油腔装配；步骤B前处理初始油膜的模型；步骤C在workbench中建立流程图，对静压推力轴承进行热力耦合计算分析；步骤D封闭工作台和油腔之间的间隙；步骤E获取间隙润滑油膜的模型。本发明应用于润滑油膜模型领域。
文档编号G06F17/50GK103246786SQ20131019864
公开日2013年8月14日 申请日期2013年5月24日 优先权日2013年5月24日
发明者张艳芹, 范立国, 于晓东, 李锐, 杨晓冬, 陈瑶 申请人:哈尔滨理工大学