基于无限长直同心圆柱的挤压膜力计算方法、系统及介质与流程

1.本发明涉及流体润滑技术领域，具体涉及基于无限长直同心圆柱的挤压膜力计算方法、系统及介质。


背景技术：

2.蒸汽发生器是核反应堆一回路系统中关键设备之一，传热管是其核心组成部分。传热管在运行过程中传热管会发生流致振动现象；带有小间隙支撑的传热管振动一直是结构动力学的研究热点及难点，也是影响传热管振动磨蚀的关键输入。
3.目前对小间隙支撑反力的描述多采用唯象的非线性弹簧模型，有学者使用润滑理论计算同心圆柱挤压膜力，但模型从多角度简化，所得的公式有不同的使用范围。缺乏统一表达的挤压膜力公式限制了对蒸汽发生器传热管振动的全局稳定性研究，制约了对蒸汽发生器传热管流致振动机制的认识，可能导致工程上未能预见的损坏失效。


技术实现要素：

4.本发明所要解决的技术问题是：传统的小间隙支撑反力挤压膜力计算模型较为局限，缺乏统一表示，限制了对传热管振动的全局稳定性研究；本发明目的在于提供基于无限长直同心圆柱的挤压膜力计算方法、系统及介质，将小间隙结构简化为无限长直同心圆柱，通过对无限长直同心圆柱进行流体控制方程量级分析，获得简化数学模型；对简化数学模型进行近似求解获得挤压膜力模型；获得一般的广泛使用的挤压膜力模型，为小间隙支承管束振动分析奠定重要理论基础。
5.本发明通过下述技术方案实现：
6.本发明提供基于无限长直同心圆柱的挤压膜力计算方法，包括步骤：
7.步骤一：将小间隙支撑结构简化为无限长直同心圆柱；
8.步骤二：基于无限长直同心圆柱进行流体控制方程量级分析，获得简化数学模型；对简化数学模型进行近似求解获得挤压膜力模型；
9.步骤三：基于挤压膜力模型计算出小间隙结构的挤压膜力。
10.进一步优化方案为，所述小间隙结构包括：振动管件、支撑件和防振件；
11.所述防振件用于限制振动管件的振动幅度，支撑件用以提供振动管件的横向刚度；
12.振动管件与支撑件或防振件之间具有间隙，且间隙内具有间隙流体。
13.本方案工作原理：传统的小间隙支撑反力挤压膜力计算模型较为局限，缺乏统一表示，限制了对传热管振动的全局稳定性研究；本发明目的在于提供基于无限长直同心圆柱的挤压膜力计算方法及系统，将小间隙结构简化为无限长直同心圆柱，通过对无限长直同心圆柱进行流体控制方程量级分析，获得简化数学模型；对简化数学模型进行近似求解获得挤压膜力模型；获得一般的广泛使用的挤压膜力模型，为小间隙支承管束振动分析奠定重要理论基础。
14.进一步优化方案为，所述振动管件的直径远大于间隙厚度，且间隙流体为单一均匀流体，流动绝热不可压缩，雷诺数re在103量级以下，周向流动速度沿径向呈抛物线分布。
15.进一步优化方案为，所述振动管件为核反应堆一回路系统中蒸汽发生器的传热管，所述传热管为圆柱形，所述间隙流体为蒸汽发生器二次侧的水。
16.进一步优化方案为，步骤二包括以下子步骤：
17.步骤2.1，基于无量纲化方法对流体控制方程进行量级分析，获得简化数学模型；
18.步骤2.2，基于积分动量方程和边界条件对简化数学模型进行近似求解；
19.步骤2.3，以简化数学模型的近似求解作为挤压膜力模型。
20.进一步优化方案为，所述步骤2.1包括过程：
21.基于无量纲化方法，在极坐标下有t＝ωt，ur(r,θ,t)＝cωu(α,θ,t)，u
θ
(r,θ,t)＝aωv(α,θ,t)，p(r,θ,t)＝ρa2ω2p(α,θ,t)；
22.在时；将navier-stokes方程简化获得简化数学模型：
[0023][0024]
其中，re≡c2ων；r表示极坐标中的极径，θ表示极坐标中的角度，ur表示间隙流体的径向速度，u
θ
表示间隙流体的周向速度，b表示支撑件的半径，ω表示振动管件的振动特征频率，c表示振动管件平衡状态下与支撑件或防振件之间的间隙；t表示时间；ρ是间隙流体的密度，ν是间隙流体的运动学粘度；p(r,θ,t)表示压强
[0025]
无量纲边界条件为：
[0026][0027][0028]
v(α,0,t)＝v(α,π,t)＝0
[0029]
其中e(t)≡e(t)/c，e(t)为振动管件的偏心度e关于时间t的函数，|e
max
＜c；＜c；
[0030]
进一步优化方案为，所述步骤2.2包括过程：
[0031]
将积分动量方程表示为：
[0032][0033]
将所述积分动量方程与无量纲边界条件联立后解得
[0034][0035]
因此
[0036][0037][0038]
进一步优化方案为，挤压膜力模型获取过程包括：
[0039]
积分平均的压力梯度得到：
[0040][0041][0042]
积分后，得到挤压膜力模型：
[0043][0044]
其中f0表示间隙流体的总流体力。
[0045]
本方案还提供基于无限长直同心圆柱的挤压膜力计算系统，用于实现上述方案所述的基于无限长直同心圆柱的挤压膜力计算方法，包括：
[0046]
预处理模块，用于将小间隙支撑结构简化为无限长直同心圆柱；
[0047]
建模模块，用于基于无限长直同心圆柱进行流体控制方程量级分析，获得简化数学模型；对简化数学模型进行近似求解获得挤压膜力模型；
[0048]
计算模块，基于挤压膜力模型计算出小间隙结构的挤压膜力。
[0049]
本方案还提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时能实现上述方案所述的基于无限长直同心圆柱的挤压膜力计算方法。
[0050]
本发明与现有技术相比，具有如下的优点和有益效果：
[0051]
本发明提供的基于无限长直同心圆柱的挤压膜力计算方法及系统及介质；将小间隙结构简化为无限长直同心圆柱，通过对无限长直同心圆柱进行流体控制方程量级分析，获得简化数学模型；对简化数学模型进行近似求解获得挤压膜力模型；获得一般的广泛使用的挤压膜力模型，为小间隙支承管束振动分析奠定重要理论基础。
附图说明
[0052]
为了更清楚地说明本发明示例性实施方式的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被
看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。在附图中：
[0053]
图1为基于无限长直同心圆柱的挤压膜力计算方法流程示意图；
[0054]
图2为挤压膜力模型获取方法流程示意图；
[0055]
图3为小间隙支撑结构简化原理图。
具体实施方式
[0056]
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。
[0057]
实施例1
[0058]
本实施例提供基于无限长直同心圆柱的挤压膜力计算方法，如图1所示，包括步骤：
[0059]
步骤一：将小间隙支撑结构简化为无限长直同心圆柱；
[0060]
步骤二：基于无限长直同心圆柱进行流体控制方程量级分析，获得简化数学模型；对简化数学模型进行近似求解获得挤压膜力模型；
[0061]
步骤三：基于挤压膜力模型计算出小间隙结构的挤压膜力。
[0062]
所述小间隙结构包括：振动管件、支撑件和防振件；
[0063]
所述防振件用于限制振动管件的振动幅度，支撑件用以提供振动管件的横向刚度；
[0064]
振动管件与支撑件或防振件之间具有间隙，且间隙内具有间隙流体。
[0065]
所述振动管件的直径远大于间隙厚度，且间隙流体为单一均匀流体，流动绝热不可压缩，雷诺数re在103量级以下，周向流动速度沿径向呈抛物线分布。
[0066]
所述振动管件为核反应堆一回路系统中蒸汽发生器的传热管，所述传热管为圆柱形，所述间隙流体为蒸汽发生器二次侧的水。
[0067]
本方案将小间隙支撑结构简化为无限长直同心圆柱，如图3所示，图中大圆表示支撑件，小实线圆表示振动管件在偏心度为e的位置，虚线的小圆表示振动管件在平衡位置，点o表示平衡点的圆心，o’表示平偏心度为e的圆心，以o点为极坐标的原点，r为极径。本发明针对无限长直同心圆柱绝热不可压缩挤压膜力，使用积分平均方法，通过严格的理论推导，获得了小间隙情形下具有一阶精度的挤压膜力模型；该挤压膜力模型是研究小间隙支承管束振动的重要理论基础，也是诸如轴承润滑、挤压膜阻尼器等微流体机械模型问题的理论解。
[0068]
如图2所示，步骤二包括以下子步骤：
[0069]
步骤2.1，基于无量纲化方法对流体控制方程进行量级分析，获得简化数学模型；
[0070]
步骤2.2，基于积分动量方程和边界条件对简化数学模型进行近似求解；
[0071]
步骤2.3，以简化数学模型的近似求解作为挤压膜力模型。
[0072]
所述步骤2.1包括过程：
[0073]
基于无量纲化方法，在极坐标下有t＝ωt，ur(r,θ,t)＝cωu(α,θ,t)，u
θ
(r,θ,t)＝aωv(α,θ,t)，p(r,θ,t)＝ρa2ω2p(α,θ,t)；
[0074]
在时；将navier-stokes方程简化获得简化数学模型：
[0075][0076]
其中，re≡c2ων；r表示极坐标中的极径，θ表示极坐标中的角度，ur表示间隙流体的径向速度，u
θ
表示间隙流体的周向速度，b表示支撑件的半径，ω表示振动管件的振动特征频率，c表示振动管件平衡状态下与支撑件或防振件之间的间隙；t表示时间；ρ是间隙流体的密度，ν是间隙流体的运动学粘度；p(r,θ,t)表示压强；
[0077]
无量纲边界条件为：
[0078][0079][0080]
v(α,0,t)＝v(α,π,t)＝0
[0081]
其中e(t)≡e(t)/c，e(t)为振动管件的偏心度e关于时间t的函数，|e
max
＜c；＜c；
[0082]
所述步骤2.2包括过程：
[0083]
将积分动量方程表示为：
[0084][0085]
将所述积分动量方程与无量纲边界条件联立后解得
[0086][0087]
因此
[0088]
[0089][0090]
挤压膜力模型获取过程包括：
[0091]
积分平均的压力梯度得到：
[0092][0093][0094]
积分后，得到挤压膜力模型：
[0095][0096]
其中f0表示间隙流体的总流体力。
[0097]
实施例2
[0098]
本实施例提供基于无限长直同心圆柱的挤压膜力计算系统，用于实现实施例1所述的基于无限长直同心圆柱的挤压膜力计算方法，包括：
[0099]
预处理模块，用于将小间隙支撑结构简化为无限长直同心圆柱；
[0100]
建模模块，用于基于无限长直同心圆柱进行流体控制方程量级分析，获得简化数学模型；对简化数学模型进行近似求解获得挤压膜力模型；
[0101]
计算模块，基于挤压膜力模型计算出小间隙结构的挤压膜力。
[0102]
实施例3
[0103]
本实施例提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时能实现根据实施例1所述的基于无限长直同心圆柱的挤压膜力计算方法。
[0104]
本发明方法得出的挤压膜力模型，用于无限长直同心圆柱间充满单一均匀流体，流体受内圆柱运动诱导产生挤压膜反力的计算，要求圆柱间间隙相对于直径为小量，流动绝热不可压缩，雷诺数re在103量级以下，周向流动速度沿径向接近于抛物线分布。
[0105]
该模型是研究小间隙支承管束振动的重要理论基础，也是诸如轴承润滑、挤压膜阻尼器等微流体机械模型问题的近似解。
[0106]
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。