一种基于Fluent的径向油膜轴承的动网格更新方法技术

一种基于Fluent的径向油膜轴承的动网格更新方法，包括以下步骤：（1）创建轴承油膜网格文件；（2）Fluent软件中UDF控制宏程序的编写；（4）利用Excel工具对（3）中得到的数据文件进行处理；（5）启动流体计算软件Fluent，导入网格文件，挂载UDF程序并设置UDF控制条件；（6）利用UDF动网格程序控制流场计算；（7）判断轴颈中心轨迹是否收敛。本发明专利技术解决了轴承油膜网格长宽比过大容易在动网格运动中发生网格畸变的问题，保证了移动前后径向方向的网格线沿圆周方向上等均分布。利用该方法计算动网格节点坐标的过程中，使网格线始终垂直于轴颈表面，网格不发生扭曲变形，避免了网格计算的累计误差，加快了轴承性能的计算速度。

A dynamic mesh updating method for radial oil film bearing based on Fluent

A radial oil film bearing Fluent dynamic update method based on grid, which comprises the following steps: (1) to create the bearing oil film mesh file; (2) the preparation of UDF macro program control software in Fluent; (4) by the tool of Excel (3) in the data file for processing; (5) start the fluid calculation software Fluent, import the grid file, mount UDF program and set UDF control conditions; (6) using UDF dynamic grid program control flow field calculation; (7) to determine whether the convergence trajectory of Journal center. The invention solves the problem that the mesh length of the oil film is too large and the grid distortion is easy to occur in the moving mesh movement, so that the grid lines in the radial direction of the moving front and the back are evenly distributed along the circumferential direction. Using this method to calculate the dynamic grid coordinates of the nodes in the process of the grid lines perpendicular to the shaft neck surface, mesh distortion does not occur, the cumulative error of grid computing is avoided, speed up the calculation of bearing performance.

【技术实现步骤摘要】
一种基于Fluent的径向油膜轴承的动网格更新方法
本专利技术属于油膜轴承及其转子系统性能计算
，具体涉及一种基于Fluent的径向油膜轴承的动网格更新方法。
技术介绍
进行油膜轴承动态性能计算时，转子的轴颈中心在涡动，轴颈表面各点的绝对速度和加速度不一样，计算流场的形状和运动学参数时刻在改变，利用静网格模型不能进行计算。动网格模型可以用来模拟由于边界运动而引起的流域随时间的变化，在单相流和多相流下均适用。边界运动可以是规定动作（例如，用户可以定义刚体中心的速度和角速度），也可以是需要待求解的动作（例如，六自由度模型下以求解结果确定下一时刻的速度、加速度、角速度、角加速度）。对于动网格模型，用户需要确定初始网格和运动区域。常用的流体性能计算软件Fluent提供了三种动网格更新算法：Smoothing、Layering、Remeshing。根据运动边界的变化，Fluent软件对指定的网格自动更新。Smoothing方法即为光顺法，特点是对网格节点间的拓扑关系不改变，只进行网格形状的改变，当网格出现畸变时容易出现发散。Layering方法是随着动边界的移动，在边界处发生网格的增加或合并。Remeshing方法是将控制区内的所有网格重新划分，一般适用于非结构化网格。考虑到油膜径向尺寸很小的特点，从理论上来说，Smoothing动网格更新方式更适合用于轴承性能计算中，但会出现网格畸变。曾有人提出了基于弹性变形的动网格更新方法。该方法与Smoothing方法中的Laplace方法类似。网格节点位移按以下公式计算：（1）（2）式中，为膜厚方向网格的层数，为网格节点所在的层数，和是轴颈中心的位移。基于弹性变形的动网格调整法存在两个问题：首先，一旦偏心率过大，垂直于轴颈的网格线就会出现严重倾斜，这种网格更新算法容易出现发散、停止、闪退等现象，导致计算结果不准确；第二，这样的计算方式容易造成累积误差，长时间计算会有网格畸变甚至负体积的情况。纵观流场计算的动网格算法用于油膜轴承性能计算时，由于油膜轴承径向尺寸较小，网格发生变形时，极容易出现网格畸变和负体积网格，导致轴承性能无法计算或者累积误差较大。
技术实现思路
本专利技术为了解决现有技术中存在轴承油膜网格长宽比过大容易在动网格中发生网格畸变的问题，提供一种基于Fluent的径向油膜轴承的动网格更新方法。为解决上述技术问题，本专利技术采用如下技术方案：一种基于Fluent的径向油膜轴承的动网格更新方法，包括以下步骤：（1）创建轴承油膜网格文件；（2）Fluent软件中UDF控制宏程序的编写：包括三个宏：定义轴颈表面各点速度的宏DEFINE_PROFILE、自定义动网格更新方式的宏DEFINE_GRID_MOTION和用于创建网格节点数据文本文件的自定义宏DEFINE_GRID_MOTION；（3）启动流体计算软件Fluent，导入网格文件，挂载用于创建网格节点数据文本的宏DEFINE_GRID_MOTION，得到网格节点数据文件；（4）利用Excel工具对步骤（3）中得到的数据文件进行处理；（5）启动流体计算软件Fluent，导入网格文件，挂载UDF程序并设置UDF控制条件；（6）利用UDF动网格程序控制流场计算；（7）判断轴颈中心轨迹是否收敛。所述的步骤（1）具体为：根据研究对象选取有限单元的类型和网格划分方法，借助于网格划分软件Gambit将计算流场划分为混合网格，对于轴承油膜的进油口和进油槽部分采用非结构化网格，对油膜部分采用结构化网格。所述的步骤（2）中，由于轴颈中心存在变位速度，其坐标在时刻改变，所以轴颈表面各点的速度根据相对运动的公式进行计算，计算出各点的速度数值不一样，通过编写宏DEFINE_PROFILE程序控制轴颈表面各点的速度。所述的步骤（2）中，利用宏DEFINE_GRID_MOTION进行动网格程序的更新编写，控制指向轴颈中心的网格线方向始终不改变，从而避免网格倾斜，减少计算累计误差。所述的步骤（2）中，在Fluent软件中，利用自定义动网格技术计算轴承油膜特性时，需要建立4个文本文件以记录每个时间步的动力学数据，以便于用于下一时间步的计算；文件1记录轴颈受的油膜力；文件2记录每一时间步的轴颈中心的坐标；文件3记录每一时间步的轴颈中心的变位速度；文件4记录轴瓦的坐标。所述的步骤（3）中，在Fluent软件中，要精确控制油膜中的所有网格节点，需要获取网格节点的全局编号、所在层数和相对于轴颈中心的单位向量，并把这些数据保存在文件5中，以方便查找和控制；以上五个文本文件需要利用宏DEFINE_GRID_MOTION程序进行编写，在动网格预览中执行已编译的程序，得到所需的上述网格节点文本文件。所述的步骤（5）中，控制条件的设置分别为：ProblemSetup设置：选择mixture模型，定义润滑油物性参数，液相与气相之间的质量传递选择Singhal气化方式，设置进油口压力、出油口压力以及操作压力，转速选择UDF中已定义的宏DEFINE_PROFILE，动网格选择UDF中已定义的宏DEFINE_GRID_MOTION；Solution设置：压力-速度耦合项选择PISO，监视器中的连续性残差、速度残差、气相残差均设置，监视进油口和出油口的质量流量，其余均保持默认设置；初始化流场置零；时间步长设置为0.0001~0.0005，时间步数为10000；其余默认，点击Calculate计算。使网格不发生扭曲变形，就需要指向轴颈中心的网格线方向始终不改变，始终垂直于轴颈表面，同时要保证在移动前后径向方向的网格线沿圆周方向上等均分布，使径向网格线的夹角相等，即移动前的网格线与移动后的网格线共线或平行；所述的步骤（6）中的计算方法具体为：假设移动前和移动后轴承的中心点分别为和，移动前和移动后的轴颈的中心点为和，油膜网格在径向方向上分成层网格，假设为时刻第层网格上任一待更新节点，为节点更新后对应的第层网格节点；根据此刻的轴颈受的油膜力、轴颈坐标、轴颈中心变位速度、轴颈质量可以得出下一时刻的轴颈中心坐标，根据轴颈中心坐标和该节点相对轴颈中心的单位向量，做射线分别交轴颈圆和轴承圆于点、，根据移动前后网格线共线或平行，可知与平行或共线；又知待求点所在的层数为，即已知，为轴承半径，为轴颈半径；坐标求解算法步骤如下：（1）、求出；（2）、已知的单位法向量，由于与共线，所以得出的单位法向量即为；（3）、求解与的夹角余弦（4）、利用三角形关系式求解（5）、求出（6）、，已知坐标，即可求出坐标。所述的步骤（5）具体为：首先设定一个初始速度和初始位移，通过udf函数DEFINE_GRID_MOTION计算出油膜力和，再根据牛顿第二定律计算出该时间步的轴颈中心坐标、速度；如果需要，再计算出轴瓦的中心坐标以及速度；将计算出来的数据分别导入对应的文本文件：文件1、文件2、文件3和文件4；同时利用文件5精确控制网格，进行网格更新，在该时间步内迭代完后进入下一个时间步；下一个时间步用到文件2、文件3和文件4最后一次保存的数值；以此循环，直至轴颈中心趋近于一定点或者稳定的轨迹；当轴颈中心轨迹趋于一定点，或者轴颈中心轨迹稳定为某一封闭曲线时，判断为收敛，否则继续执行步骤（4）直到稳定为止。采用上述技术本文档来自技高网...

【技术保护点】
一种基于Fluent的径向油膜轴承的动网格更新方法，其特征在于：包括以下步骤：（1）创建轴承油膜网格文件；（2）Fluent软件中UDF控制宏程序的编写：包括三个宏：定义轴颈表面各点速度的宏DEFINE_PROFILE、自定义动网格更新方式的宏DEFINE_GRID_MOTION和用于创建网格节点数据文本文件的自定义宏DEFINE_GRID_MOTION；（3）启动流体计算软件Fluent，导入网格文件，挂载用于创建网格节点数据文本的宏DEFINE_GRID_MOTION，得到网格节点数据文件；（4）利用Excel工具对步骤（3）中得到的数据文件进行处理；（5）启动流体计算软件Fluent，导入网格文件，挂载UDF程序并设置UDF控制条件；（6）利用UDF动网格程序控制流场计算；（7）判断轴颈中心轨迹是否收敛。

【技术特征摘要】
1.一种基于Fluent的径向油膜轴承的动网格更新方法，其特征在于：包括以下步骤：（1）创建轴承油膜网格文件；（2）Fluent软件中UDF控制宏程序的编写：包括三个宏：定义轴颈表面各点速度的宏DEFINE_PROFILE、自定义动网格更新方式的宏DEFINE_GRID_MOTION和用于创建网格节点数据文本文件的自定义宏DEFINE_GRID_MOTION；（3）启动流体计算软件Fluent，导入网格文件，挂载用于创建网格节点数据文本的宏DEFINE_GRID_MOTION，得到网格节点数据文件；（4）利用Excel工具对步骤（3）中得到的数据文件进行处理；（5）启动流体计算软件Fluent，导入网格文件，挂载UDF程序并设置UDF控制条件；（6）利用UDF动网格程序控制流场计算；（7）判断轴颈中心轨迹是否收敛。2.根据权利要求1所述的一种基于Fluent的径向油膜轴承的动网格更新方法，其特征在于：所述的步骤（1）具体为：根据研究对象选取有限单元的类型和网格划分方法，借助于网格划分软件Gambit将计算流场划分为混合网格，对于轴承油膜的进油口和进油槽部分采用非结构化网格，对油膜部分采用结构化网格。3.根据权利要求1所述的一种基于Fluent的径向油膜轴承的动网格更新方法，其特征在于：所述的步骤（2）中，由于轴颈中心存在变位速度，其坐标在时刻改变，所以轴颈表面各点的速度根据相对运动的公式进行计算，计算出各点的速度数值不一样，通过编写宏DEFINE_PROFILE程序控制轴颈表面各点的速度。4.根据权利要求3所述的一种基于Fluent的径向油膜轴承的动网格更新方法，其特征在于：所述的步骤（2）中，利用宏DEFINE_GRID_MOTION进行动网格程序的更新编写，控制指向轴颈中心的网格线方向始终不改变，从而避免网格倾斜，减少计算累计误差。5.根据权利要求4所述的一种基于Fluent的径向油膜轴承的动网格更新方法，其特征在于：所述的步骤（2）中，在Fluent软件中，利用自定义动网格技术计算轴承油膜特性时，需要建立4个文本文件以记录每个时间步的动力学数据，以便于用于下一时间步的计算；文件1记录轴颈受的油膜力；文件2记录每一时间步的轴颈中心的坐标；文件3记录每一时间步的轴颈中心的变位速度；文件4记录轴瓦的坐标。6.根据权利要求5所述的一种基于Fluent的径向油膜轴承的动网格更新方法，其特征在于：所述的步骤（3）中，在Fluent软件中，要精确控制油膜中的所有网格节点，需要获取网格节点的全局编号、所在层数和相对于轴颈中心的单位向量，并把这些数据保存在文件5中，以方便查找和控制；以上五个文本文件需要利用宏DEFINE_GRID_MOTION程序进行编写，在动网格预览中执行已编译的程序，得到所需的上述网格节点文本文件。...

【专利技术属性】
技术研发人员：吴超，尹雪梅，李蒙蒙，朱茹敏，马明飞，李立伟，任菲，郭志强，杨春燕，巩晓赟，张德海，崔晓康，
申请(专利权)人：郑州轻工业学院，
类型：发明
国别省市：河南,41