一种微小间隙内油膜网格划分方法及系统的制作方法
【专利摘要】本发明公开了一种微小间隙内油膜网格划分方法及系统，具体的，对微小间隙内的油膜进行几何变换，将油膜厚度方向上的尺寸放大，使其与轴承直径尺度达到同一量级；然后对变换后的油膜进行网格划分；最后再对划分好的网格坐标进行逆变换。通过坐标变换的方法，将滑动轴承微小间隙内油膜尺度放大，使其在径向上的尺度和轴向、周向上的尺度相当，对放大后的油膜进行网格划分，再通过坐标逆变换的方法，将划分好的网格坐标转换回去，得到真实的油膜网格坐标。本发明解决了滑动轴承微小间隙内因不同方向上尺度相差巨大所带来的网格划分困难问题，可以应用于基于计算流体力学技术的滑动轴承性能分析。
【专利说明】
-种微小间隙内油膜网格划分方法及系统
技术领域
[0001] 本发明设及轴承检测技术领域，更具体的，本发明设及一种微小间隙内油膜网格 划分方法及系统。
【背景技术】
[0002] 滑动轴承被广泛应用于汽轮机、发电机、压缩机和累等各类旋转机械，起着支撑转 子的重要作用，对转子动力特性和机组安全运行的影响很大，分析和掌握滑动轴承性能具 有重要意义。随着计算技术的发展，计算流体力学技术在滑动轴承性能分析中的应用越来 越普遍，目前开始在原始纳维艾一斯托克斯方程基础上开展滑动轴承性能分析，运种方法 可W考虑任意复杂轴承结构型式，能更全面、准确地反映轴承内流体特性。
[0003] 基于计算流体力学技术开展滑动轴承性能计算分析前，首先需要对滑动轴承间隙 内的油膜进行空间离散，划分油膜计算网格，根据计算流体力学理论可知，网格划分质量对 于轴承性能计算准确性至关重要。但是由于滑动轴承轴颈和轴瓦之间的油膜很薄，其厚度 只有轴承直径和长度的1.5%〇~3%〇,轴承间隙内油膜径向尺寸相对于周向和轴向尺寸很 小，不是同一尺度，网格划分时多条轮廓线混为一条曲线，采用传统网格划分软件时，如不 对轴承间隙进行放大，油膜轮廓线很难精确选取，甚至不能选取;如果对轴承间隙充分放 大，油膜轮廓线又有可能出现线型失真，甚至产生不正确的几何结构，因此，由于滑动轴承 微小间隙的特点，使得油膜计算网格划分困难，而且所划分的网格质量较差，影响计算结 果。
【发明内容】
[0004] 本发明的目的在于提出一种微小间隙内油膜网格划分方法及系统，W解决现有技 术中由于滑动轴承微小间隙的特点，使得油膜计算网格划分困难，而且所划分的网格质量 较差，影响计算结果的问题。
[0005] 为达到上述目的，本发明提供了 W下技术方案：
[0006] 本发明公开了 一种微小间隙内油膜网格划分方法，包括：
[0007] 对滑动轴承微小间隙内的油膜进行几何变换，并选取油膜轮廓线；
[000引将所述油膜的厚度方向上的尺寸放大，使其与所述滑动轴承的直径尺度达到同一 量级；
[0009] 对变换后的所述油膜进行网格划分，将划分好的网格进行坐标逆转换，得到所述 微小间隙内油膜网格坐标。
[0010] 其中，所述对滑动轴承微小间隙内的油膜进行几何变换，并选取油膜轮廓线具体 为：
[0011] 对所述滑动轴承间隙进行放大，其轴径位置与尺寸不变，选取缩放因子AW所述滑 动轴承的圆屯、为基点，将所述滑动轴承的半径放大；
[0012] 选取所述滑动轴承上的所述油膜轮廓线。
[0013] 其中，将所述油膜的厚度方向上的尺寸放大，使其与所述滑动轴承的直径尺度达 到同一量级具体为：
[0014] 根据所述缩放因子A将所述滑动轴承的半径放大，放大后的所述油膜的厚度为h = (l+ecos0)(AR-r)/2；
[0015] 根据选取的所述所述缩放因子A将所述所述油膜的厚度与所述滑动轴承的直径尺 度达到同一量级。
[0016] 其中，所述对变换后的所述油膜进行网格划分，将划分好的网格进行坐标逆转换， 得到所述微小间隙内油膜网格坐标具体为：
[0017] 将变换后的所述油膜进行网格划分，得到网格节点；
[0018] 对所述网格节点进行编号，得到所述网格节点的坐标；
[0019] 根据网格坐标逆变换将所述网格节点的坐标进行坐标逆转换，得到所述微小间隙 内油膜网格坐标。
[0020] 本发明还公开了一种微小间隙内油膜网格划分系统，包括：
[0021] 变换及选取单元，用于对滑动轴承微小间隙内的油膜进行几何变换，并选取油膜 轮廓线；
[0022] 尺寸放大单元，用于将所述油膜的厚度方向上的尺寸放大，使其与所述滑动轴承 的直径尺度达到同一量级；
[0023] 逆变换单元，用于对变换后的所述油膜进行网格划分，将划分好的网格进行坐标 逆变换，得到所述微小间隙内油膜网格坐标。
[0024] 其中，所述变换及选取单元包括：
[0025] 放大单元，用于对所述滑动轴承间隙进行放大，其轴径位置与尺寸不变，选取缩放 因子AW所述滑动轴承的圆屯、为基点，将所述滑动轴承的半径放大；
[0026] 选取单元，用于选取所述滑动轴承上的所述油膜轮廓线。
[0027] 其中，所述放大单元包括：
[0028] 半径放大单元，用于根据所述缩放因子A将所述滑动轴承的半径放大，放大后的所 述油膜的厚度为h=(l+ecos0)(AR-r)/2;
[0029] 变换单元，用于根据选取的所述缩放因子A将所述油膜的厚度与所述滑动轴承的 直径尺度变换达到同一量级。
[0030] 其中，所述逆变换单元包括：
[0031] 网格划分单元，用于将变换后的所述油膜进行网格划分，得到网格节点；
[0032] 编号单元，用于对所述网格节点进行编号，得到所述网格节点的坐标；
[0033] 坐标逆变换单元，用于根据网格坐标逆变换将所述网格节点的坐标进行坐标逆转 换，得到所述微小间隙内油膜网格坐标。
[0034] 经由上述的技术方案可知，与现有技术相比，本发明公开了一种微小间隙内油膜 网格划分方法及系统，具体的，对微小间隙内的油膜进行几何变换，将油膜厚度方向上的尺 寸放大，使其与轴承直径尺度达到同一量级;然后对变换后的油膜进行网格划分;最后再对 划分好的网格坐标进行逆变换。通过坐标变换的方法，将滑动轴承微小间隙内油膜尺度放 大，使其在径向上的尺度和轴向、周向上的尺度相当，对放大后的油膜进行网格划分，再通 过坐标逆变换的方法，将划分好的网格坐标转换回去，得到真实的油膜网格坐标。本发明解 决了滑动轴承微小间隙内因不同方向上尺度相差巨大所带来的网格划分困难问题，可W应 用于基于计算流体力学技术的滑动轴承性能分析。
【附图说明】
[0035] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可W根据 提供的附图获得其他的附图。
[0036] 图1为滑动轴承油膜的几何示意图；
[0037] 图2为本发明公开的一种微小间隙内油膜网格划分方法的流程示意图；
[0038] 图3为应用本发明提供的方法后的轴承油膜几何示意图；
[0039] 图4为应用本发明提供的方法后的离散的油膜划分网格示意图；
[0040] 图5为图2中步骤S203的具体方法的流程图；
[0041] 图6为本发明公开的一种微小间隙内油膜网格划分系统的结构示意图；
[0042] 图7为图6中逆变换单元的具体结构示意图。
【具体实施方式】
[0043] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0044] 将滑动轴承油膜看作为由油膜主体a、油槽b、进油口 C和轴瓦内表面d组成，轴颈外 表面和轴瓦内表面之间的间隙构成了油膜。附图1给出了滑动轴承油膜几何形状示意图。因 为油膜厚度非常微小，在图中无法将油膜和轴瓦分辨开来。将局部视图放大百余倍后，才可 W看出油膜形状。
[0045] 请参阅附图2,图2为本发明公开的一种微小间隙内油膜网格划分方法的流程示意 图。本发明公开了 一种微小间隙内油膜网格划分方法，具体包括如下步骤：
[0046] S201、对滑动轴承微小间隙内的油膜进行几何变换，并选取油膜轮廓线。
[0047] 具体的，几何变换时，轴颈位置和尺寸不变。选取缩放因子A，W轴承圆屯、为基点， 将轴承半径R放大M音，由R变为AR,经放大变换处理后的滑动轴承几何形状示意图如图3所 示，图3为应用本发明提供的方法后的轴承油膜几何示意图。轴承半径放大后，进油口和油 槽位置将发生偏移，关键点a、b和C点的纵坐标保持不变，横坐标变化如表1所示。
[004引表1 a、b和C点的坐标变换
[0049]
[0050] S202、将油膜的厚度方向上的尺寸放大，使其与滑动轴承的直径尺度达到同一量 级。
[0051] 根据上述变换W后，其油膜的厚度发生了变化：
[0052]变换前的油膜厚度为:h=(l+ecos目）（R-r)/2 [0053] 变换后的油膜厚度为:h=(l+ecos目）（AR-r)/2
[0054] 上述式子中，r为轴颈半径，R为轴承半径，A为缩放大因子，选定足够大的放大因子 入，就能够使油膜厚度h与轴承半径R达到同一量级。
[0055] S203、对变换后的油膜进行网格划分，将划分好的网格进行坐标逆转换，得到微小 间隙内油膜网格坐标。
[0056] 通过上述变换后，油膜厚度h与轴承半径R达到同一量级，可W采用常规网格划分 方法对变换放大后的油膜进行网格划分，划分后的网格如图4所示。将油膜网格标记为=个 区：油膜主体A区、油槽B区和和进油口 C区。
[0057] 如图5所示，图5为图2中步骤S203的具体流程图。对变换后的油膜进行网格划分， 将划分好的网格进行坐标逆转换，得到微小间隙内油膜网格坐标具体为：
[005引S501、将变换后的油膜进行网格划分，得到网格节点；
[0059] S502、对网格节点进行编号，得到网格节点的坐标；
[0060] S503、根据网格坐标逆变换将网格节点的坐标进行坐标逆转换，得到微小间隙内 油膜网格坐标。
[0061] 具体的，通过常规网格划分方法对变换放大后的油膜进行网格划分，划分后的网 格如图4所示。将油膜网格标记为立个区：油膜主体A区、油槽B区和和进油口C区，并对网格 节点进行编号。WNi,J表示任一网格节点，下标i表示径向层数，轴颈最内侧为1，径向逐层递 增;下标j表示周向节点号，逆时针递增。
[0062] 上述对所划分好的网格坐标进行逆变换运算，得到得到微小间隙内油膜网格坐 标，具体的网格坐标逆变换计算方法如下：
[00创计算节点Ni,j与x轴的夹角：目 = 1:an-l(yi,j-yl,j，Xi,j-Xl,j)，-3T<白<JT [0064] 对于A区节点：
[0067] 对于B区和C区节点：
[00 化]
[0066]
[006引
[0069]
[0070] 上述式子中I, J，/ I, J为实际油膜网格坐标。
[0071 ]本发明公开了一种微小间隙内油膜网格划分方法具体的，对微小间隙内的油膜进 行几何变换，将油膜厚度方向上的尺寸放大，使其与轴承直径尺度达到同一量级;然后对变 换后的油膜进行网格划分;最后再对划分好的网格坐标进行逆变换。通过坐标变换的方法， 将滑动轴承微小间隙内油膜尺度放大，使其在径向上的尺度和轴向、周向上的尺度相当，对 放大后的油膜进行网格划分，再通过坐标逆变换的方法，将划分好的网格坐标转换回去，得 到真实的油膜网格坐标。本发明解决了滑动轴承微小间隙内因不同方向上尺度相差巨大所 带来的网格划分困难问题，可W应用于基于计算流体力学技术的滑动轴承性能分析。
[0072] 请参阅附图6,图6为本发明公开的一种微小间隙内油膜网格划分系统的结构示意 图。在上述公开的方法的基础上，本发明还公开了一种微小间隙内油膜网格划分系统，该系 统具体包括：
[0073] 变换及选取单元U61，用于对滑动轴承微小间隙内的油膜进行几何变换，并选取油 膜轮廓线;尺寸放大单元U62,用于将油膜的厚度方向上的尺寸放大，使其与滑动轴承的直 径尺度达到同一量级;逆变换单元U63,用于对变换后的油膜进行网格划分，将划分好的网 格进行坐标逆变换，得到微小间隙内油膜网格坐标。
[0074] 该系统对应执行图2中的方法流程，具体其中的描述，请参阅附图2对应的相关内 容。
[0075] 具体的，所述变换及选取单元包括：
[0076] 放大单元，用于对所述滑动轴承间隙进行放大，其轴径位置与尺寸不变，选取缩放 因子AW所述滑动轴承的圆屯、为基点，将所述滑动轴承的半径放大；
[0077] 选取单元，用于选取所述滑动轴承上的所述油膜轮廓线。
[007引具体的，所述放大单元包括：
[0079] 半径放大单元，用于根据所述缩放因子A将所述滑动轴承的半径放大，放大后的所 述油膜的厚度为h=(l+ecos0)(AR-r)/2;
[0080] 变换单元，用于根据选取的所述缩放因子A将所述油膜的厚度与所述滑动轴承的 直径尺度变换达到同一量级。
[0081] 具体的，请参阅附图7,图7为图6中逆变换单元的具体结构示意图。如图所示，逆变 换单元U63具体包括：
[0082] 网格划分单元U71，用于将变换后的油膜进行网格划分，得到网格节点；编号单元 U72,用于对网格节点进行编号，得到网格节点的坐标;坐标逆变换单元U73,用于根据网格 坐标逆变换将网格节点的坐标进行坐标逆转换，得到微小间隙内油膜网格坐标。
[0083] 该系统对应执行图5中的方法流程，具体其中的描述，请参阅附图5对应的相关内 容。
[0084] 本发明公开了 一种微小间隙内油膜网格划分系统，具体的，对微小间隙内的油膜 进行几何变换，将油膜厚度方向上的尺寸放大，使其与轴承直径尺度达到同一量级;然后对 变换后的油膜进行网格划分;最后再对划分好的网格坐标进行逆变换。通过坐标变换的方 法，将滑动轴承微小间隙内油膜尺度放大，使其在径向上的尺度和轴向、周向上的尺度相 当，对放大后的油膜进行网格划分，再通过坐标逆变换的方法，将划分好的网格坐标转换回 去，得到真实的油膜网格坐标。本发明解决了滑动轴承微小间隙内因不同方向上尺度相差 巨大所带来的网格划分困难问题，可W应用于基于计算流体力学技术的滑动轴承性能分 析。
[0085] 综上所述，本发明公开了一种微小间隙内油膜网格划分方法及系统，具体的，对微 小间隙内的油膜进行几何变换，将油膜厚度方向上的尺寸放大，使其与轴承直径尺度达到 同一量级;然后对变换后的油膜进行网格划分;最后再对划分好的网格坐标进行逆变换。通 过坐标变换的方法，将滑动轴承微小间隙内油膜尺度放大，使其在径向上的尺度和轴向、周 向上的尺度相当，对放大后的油膜进行网格划分，再通过坐标逆变换的方法，将划分好的网 格坐标转换回去，得到真实的油膜网格坐标。本发明解决了滑动轴承微小间隙内因不同方 向上尺度相差巨大所带来的网格划分困难问题，可W应用于基于计算流体力学技术的滑动 轴承性能分析。
[0086] 需要说明的是，本说明书中的各个实施例均采用递进的方式描述，每个实施例重 点说明的都是与其他实施例的不同之处，各个实施例之间相同相似的部分互相参见即可。
[0087] 对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。 对运些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的 一般原理可W在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明 将不会被限制于本文所示的运些实施例，而是要符合与本文所公开的原理和新颖特点相一 致的最宽的范围。
【主权项】
1. 一种微小间隙内油膜网格划分方法，其特征在于，包括： 对滑动轴承微小间隙内的油膜进行几何变换，并选取油膜轮廓线； 将所述油膜的厚度方向上的尺寸放大，使其与所述滑动轴承的直径尺度达到同一量 级； 对变换后的所述油膜进行网格划分，将划分好的网格进行坐标逆转换，得到所述微小 间隙内油膜网格坐标。2. 根据权利要求1所述的方法，其特征在于，所述对滑动轴承微小间隙内的油膜进行几 何变换，并选取油膜轮廓线具体为： 对所述滑动轴承间隙进行放大，其轴径位置与尺寸不变，选取缩放因子λ以所述滑动轴 承的圆心为基点，将所述滑动轴承的半径放大； 选取所述滑动轴承上的所述油膜轮廓线。3. 根据权利要求1所述的方法，其特征在于，将所述油膜的厚度方向上的尺寸放大，使 其与所述滑动轴承的直径尺度达到同一量级具体为： 根据所述缩放因子λ将所述滑动轴承的半径放大，放大后的所述油膜的厚度为h=(l+e cos0)(AR-r)/2； 根据选取的所述所述缩放因子λ将所述所述油膜的厚度与所述滑动轴承的直径尺度达 到同一量级。4. 根据权利要求1所述的方法，其特征在于，所述对变换后的所述油膜进行网格划分， 将划分好的网格进行坐标逆转换，得到所述微小间隙内油膜网格坐标具体为： 将变换后的所述油膜进行网格划分，得到网格节点； 对所述网格节点进行编号，得到所述网格节点的坐标； 根据网格坐标逆变换将所述网格节点的坐标进行坐标逆转换，得到所述微小间隙内油 膜网格坐标。5. -种微小间隙内油膜网格划分系统，其特征在于，包括： 变换及选取单元，用于对滑动轴承微小间隙内的油膜进行几何变换，并选取油膜轮廓 线； 尺寸放大单元，用于将所述油膜的厚度方向上的尺寸放大，使其与所述滑动轴承的直 径尺度达到同一量级； 逆变换单元，用于对变换后的所述油膜进行网格划分，将划分好的网格进行坐标逆变 换，得到所述微小间隙内油膜网格坐标。6. 根据权利要求5所述的系统，其特征在于，所述变换及选取单元包括： 放大单元，用于对所述滑动轴承间隙进行放大，其轴径位置与尺寸不变，选取缩放因子 λ以所述滑动轴承的圆心为基点，将所述滑动轴承的半径放大； 选取单元，用于选取所述滑动轴承上的所述油膜轮廓线。7. 根据权利要求5所述的系统，其特征在于，所述放大单元包括： 半径放大单元，用于根据所述缩放因子λ将所述滑动轴承的半径放大，放大后的所述油 膜的厚度为h=(l+ecos9)(>R-；r)/2; 变换单元，用于根据选取的所述缩放因子λ将所述油膜的厚度与所述滑动轴承的直径 尺度变换达到同一量级。8.根据权利要求5所述的系统，其特征在于，所述逆变换单元包括： 网格划分单元，用于将变换后的所述油膜进行网格划分，得到网格节点； 编号单元，用于对所述网格节点进行编号，得到所述网格节点的坐标； 坐标逆变换单元，用于根据网格坐标逆变换将所述网格节点的坐标进行坐标逆转换， 得到所述微小间隙内油膜网格坐标。
【文档编号】G06F17/50GK106021775SQ201610374049
【公开日】2016年10月12日
【申请日】2016年5月30日
【发明人】张楚, 高庆水, 杨毅, 金格, 杜胜磊, 谭金, 邓小文
【申请人】广东电网有限责任公司电力科学研究院