CN115796079A - 基于无限长直同心圆柱的挤压膜力计算方法、系统及介质 - Google Patents
基于无限长直同心圆柱的挤压膜力计算方法、系统及介质 Download PDFInfo
- Publication number
- CN115796079A CN115796079A CN202211623167.4A CN202211623167A CN115796079A CN 115796079 A CN115796079 A CN 115796079A CN 202211623167 A CN202211623167 A CN 202211623167A CN 115796079 A CN115796079 A CN 115796079A
- Authority
- CN
- China
- Prior art keywords
- gap
- film force
- concentric cylinder
- extrusion film
- infinite
- 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
- 238000001125 extrusion Methods 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 32
- 239000012530 fluid Substances 0.000 claims abstract description 38
- 238000013178 mathematical model Methods 0.000 claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims abstract description 22
- 210000003722 extracellular fluid Anatomy 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 3
- 238000003860 storage Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000005461 lubrication Methods 0.000 abstract description 6
- 238000009795 derivation Methods 0.000 abstract description 2
- 238000005457 optimization Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Landscapes
- Extrusion Moulding Of Plastics Or The Like (AREA)
Abstract
本发明公开了基于无限长直同心圆柱的挤压膜力计算方法及系统,涉及流体润滑技术领域,将小间隙支撑结构简化为无限长直同心圆柱;基于无限长直同心圆柱进行流体控制方程量级分析,获得简化数学模型;对简化数学模型进行近似求解获得挤压膜力模型;最后基于挤压膜力模型计算出小间隙结构的挤压膜力;本方案针对无限长直同心圆柱绝热不可压缩挤压膜力,使用积分平均方法,通过严格的方法推导,获得了小间隙情形下具有一阶精度的挤压膜力模型;该挤压膜力模型是研究小间隙支承管束振动的重要理论基础,也是诸如轴承润滑、挤压膜阻尼器等微流体机械模型问题的理论解。
Description
技术领域
本发明涉及流体润滑技术领域,具体涉及基于无限长直同心圆柱的挤压膜力计算方法、系统及介质。
背景技术
蒸汽发生器是核反应堆一回路系统中关键设备之一,传热管是其核心组成部分。传热管在运行过程中传热管会发生流致振动现象;带有小间隙支撑的传热管振动一直是结构动力学的研究热点及难点,也是影响传热管振动磨蚀的关键输入。
目前对小间隙支撑反力的描述多采用唯象的非线性弹簧模型,有学者使用润滑理论计算同心圆柱挤压膜力,但模型从多角度简化,所得的公式有不同的使用范围。缺乏统一表达的挤压膜力公式限制了对蒸汽发生器传热管振动的全局稳定性研究,制约了对蒸汽发生器传热管流致振动机制的认识,可能导致工程上未能预见的损坏失效。
发明内容
本发明所要解决的技术问题是:传统的小间隙支撑反力挤压膜力计算模型较为局限,缺乏统一表示,限制了对传热管振动的全局稳定性研究;本发明目的在于提供基于无限长直同心圆柱的挤压膜力计算方法、系统及介质,将小间隙结构简化为无限长直同心圆柱,通过对无限长直同心圆柱进行流体控制方程量级分析,获得简化数学模型;对简化数学模型进行近似求解获得挤压膜力模型;获得一般的广泛使用的挤压膜力模型,为小间隙支承管束振动分析奠定重要理论基础。
本发明通过下述技术方案实现:
本发明提供基于无限长直同心圆柱的挤压膜力计算方法,包括步骤:
步骤一:将小间隙支撑结构简化为无限长直同心圆柱;
步骤二:基于无限长直同心圆柱进行流体控制方程量级分析,获得简化数学模型;对简化数学模型进行近似求解获得挤压膜力模型;
步骤三:基于挤压膜力模型计算出小间隙结构的挤压膜力。
进一步优化方案为,所述小间隙结构包括:振动管件、支撑件和防振件;
所述防振件用于限制振动管件的振动幅度,支撑件用以提供振动管件的横向刚度;
振动管件与支撑件或防振件之间具有间隙,且间隙内具有间隙流体。
本方案工作原理:传统的小间隙支撑反力挤压膜力计算模型较为局限,缺乏统一表示,限制了对传热管振动的全局稳定性研究;本发明目的在于提供基于无限长直同心圆柱的挤压膜力计算方法及系统,将小间隙结构简化为无限长直同心圆柱,通过对无限长直同心圆柱进行流体控制方程量级分析,获得简化数学模型;对简化数学模型进行近似求解获得挤压膜力模型;获得一般的广泛使用的挤压膜力模型,为小间隙支承管束振动分析奠定重要理论基础。
进一步优化方案为,所述振动管件的直径远大于间隙厚度,且间隙流体为单一均匀流体,流动绝热不可压缩,雷诺数Re在103量级以下,周向流动速度沿径向呈抛物线分布。
进一步优化方案为,所述振动管件为核反应堆一回路系统中蒸汽发生器的传热管,所述传热管为圆柱形,所述间隙流体为蒸汽发生器二次侧的水。
进一步优化方案为,步骤二包括以下子步骤:
步骤2.1,基于无量纲化方法对流体控制方程进行量级分析,获得简化数学模型;
步骤2.2,基于积分动量方程和边界条件对简化数学模型进行近似求解;
步骤2.3,以简化数学模型的近似求解作为挤压膜力模型。
进一步优化方案为,所述步骤2.1包括过程:
其中,Re≡c2ων;r表示极坐标中的极径,θ表示极坐标中的角度,ur表示间隙流体的径向速度,uθ表示间隙流体的周向速度,b表示支撑件的半径,ω表示振动管件的振动特征频率,c表示振动管件平衡状态下与支撑件或防振件之间的间隙;t表示时间;ρ是间隙流体的密度,ν是间隙流体的运动学粘度;p(r,θ,t)表示压强
无量纲边界条件为:
v(α,0,T)=v(α,π,T)=0
进一步优化方案为,所述步骤2.2包括过程:
将积分动量方程表示为:
将所述积分动量方程与无量纲边界条件联立后解得
因此
进一步优化方案为,挤压膜力模型获取过程包括:
积分平均的压力梯度得到:
积分后,得到挤压膜力模型:
本方案还提供基于无限长直同心圆柱的挤压膜力计算系统,用于实现上述方案所述的基于无限长直同心圆柱的挤压膜力计算方法,包括:
预处理模块,用于将小间隙支撑结构简化为无限长直同心圆柱;
建模模块,用于基于无限长直同心圆柱进行流体控制方程量级分析,获得简化数学模型;对简化数学模型进行近似求解获得挤压膜力模型;
计算模块,基于挤压膜力模型计算出小间隙结构的挤压膜力。
本方案还提供一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时能实现上述方案所述的基于无限长直同心圆柱的挤压膜力计算方法。
本发明与现有技术相比,具有如下的优点和有益效果:
本发明提供的基于无限长直同心圆柱的挤压膜力计算方法及系统及介质;将小间隙结构简化为无限长直同心圆柱,通过对无限长直同心圆柱进行流体控制方程量级分析,获得简化数学模型;对简化数学模型进行近似求解获得挤压膜力模型;获得一般的广泛使用的挤压膜力模型,为小间隙支承管束振动分析奠定重要理论基础。
附图说明
为了更清楚地说明本发明示例性实施方式的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。在附图中:
图1为基于无限长直同心圆柱的挤压膜力计算方法流程示意图;
图2为挤压膜力模型获取方法流程示意图;
图3为小间隙支撑结构简化原理图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
实施例1
本实施例提供基于无限长直同心圆柱的挤压膜力计算方法,如图1所示,包括步骤:
步骤一:将小间隙支撑结构简化为无限长直同心圆柱;
步骤二:基于无限长直同心圆柱进行流体控制方程量级分析,获得简化数学模型;对简化数学模型进行近似求解获得挤压膜力模型;
步骤三:基于挤压膜力模型计算出小间隙结构的挤压膜力。
所述小间隙结构包括:振动管件、支撑件和防振件;
所述防振件用于限制振动管件的振动幅度,支撑件用以提供振动管件的横向刚度;
振动管件与支撑件或防振件之间具有间隙,且间隙内具有间隙流体。
所述振动管件的直径远大于间隙厚度,且间隙流体为单一均匀流体,流动绝热不可压缩,雷诺数Re在103量级以下,周向流动速度沿径向呈抛物线分布。
所述振动管件为核反应堆一回路系统中蒸汽发生器的传热管,所述传热管为圆柱形,所述间隙流体为蒸汽发生器二次侧的水。
本方案将小间隙支撑结构简化为无限长直同心圆柱,如图3所示,图中大圆表示支撑件,小实线圆表示振动管件在偏心度为e的位置,虚线的小圆表示振动管件在平衡位置,点O表示平衡点的圆心,O’表示平偏心度为e的圆心,以O点为极坐标的原点,r为极径。本发明针对无限长直同心圆柱绝热不可压缩挤压膜力,使用积分平均方法,通过严格的理论推导,获得了小间隙情形下具有一阶精度的挤压膜力模型;该挤压膜力模型是研究小间隙支承管束振动的重要理论基础,也是诸如轴承润滑、挤压膜阻尼器等微流体机械模型问题的理论解。
如图2所示,步骤二包括以下子步骤:
步骤2.1,基于无量纲化方法对流体控制方程进行量级分析,获得简化数学模型;
步骤2.2,基于积分动量方程和边界条件对简化数学模型进行近似求解;
步骤2.3,以简化数学模型的近似求解作为挤压膜力模型。
所述步骤2.1包括过程:
其中,Re≡c2ων;r表示极坐标中的极径,θ表示极坐标中的角度,ur表示间隙流体的径向速度,uθ表示间隙流体的周向速度,b表示支撑件的半径,ω表示振动管件的振动特征频率,c表示振动管件平衡状态下与支撑件或防振件之间的间隙;t表示时间;ρ是间隙流体的密度,ν是间隙流体的运动学粘度;p(r,θ,t)表示压强;
无量纲边界条件为:
v(α,0,T)=v(α,π,T)=0
所述步骤2.2包括过程:
将积分动量方程表示为:
将所述积分动量方程与无量纲边界条件联立后解得
因此
挤压膜力模型获取过程包括:
积分平均的压力梯度得到:
积分后,得到挤压膜力模型:
实施例2
本实施例提供基于无限长直同心圆柱的挤压膜力计算系统,用于实现实施例1所述的基于无限长直同心圆柱的挤压膜力计算方法,包括:
预处理模块,用于将小间隙支撑结构简化为无限长直同心圆柱;
建模模块,用于基于无限长直同心圆柱进行流体控制方程量级分析,获得简化数学模型;对简化数学模型进行近似求解获得挤压膜力模型;
计算模块,基于挤压膜力模型计算出小间隙结构的挤压膜力。
实施例3
本实施例提供一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时能实现根据实施例1所述的基于无限长直同心圆柱的挤压膜力计算方法。
本发明方法得出的挤压膜力模型,用于无限长直同心圆柱间充满单一均匀流体,流体受内圆柱运动诱导产生挤压膜反力的计算,要求圆柱间间隙相对于直径为小量,流动绝热不可压缩,雷诺数Re在103量级以下,周向流动速度沿径向接近于抛物线分布。
该模型是研究小间隙支承管束振动的重要理论基础,也是诸如轴承润滑、挤压膜阻尼器等微流体机械模型问题的近似解。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.基于无限长直同心圆柱的挤压膜力计算方法,其特征在于,包括步骤:
步骤一:将小间隙支撑结构简化为无限长直同心圆柱;
步骤二:基于无限长直同心圆柱进行流体控制方程量级分析,获得简化数学模型;对简化数学模型进行近似求解获得挤压膜力模型;
步骤三:基于挤压膜力模型计算出小间隙结构的挤压膜力。
2.根据权利要求1所述的基于无限长直同心圆柱的挤压膜力计算方法,其特征在于,所述小间隙结构包括:振动管件、支撑件和防振件;
所述防振件用于限制振动管件的振动幅度,支撑件用以提供振动管件的横向刚度;
振动管件与支撑件或防振件之间具有间隙,且间隙内具有间隙流体。
3.根据权利要求2所述的基于无限长直同心圆柱的挤压膜力计算方法,其特征在于,所述振动管件的直径远大于间隙厚度,且间隙流体为单一均匀流体,流动绝热不可压缩,雷诺数Re在103量级以下,周向流动速度沿径向呈抛物线分布。
4.根据权利要求3所述的基于无限长直同心圆柱的挤压膜力计算方法,其特征在于,所述振动管件为核反应堆一回路系统中蒸汽发生器的传热管,所述传热管为圆柱形,所述间隙流体为蒸汽发生器二次侧的水。
5.根据权利要求1所述的基于无限长直同心圆柱的挤压膜力计算方法,其特征在于,步骤二包括以下子步骤:
步骤2.1,对无量纲流体控制方程进行量级分析,获得简化数学模型;
步骤2.2,基于积分动量方程和边界条件对简化数学模型进行近似求解;
步骤2.3,以简化数学模型的近似求解作为挤压膜力模型。
6.根据权利要求5所述的基于无限长直同心圆柱的挤压膜力计算方法,其特征在于,所述步骤2.1包括过程:
其中,Re≡c2ων;r表示极坐标中的极径,θ表示极坐标中的角度,ur表示间隙流体的径向速度,uθ表示间隙流体的周向速度,b表示支撑件的半径,ω表示振动管件的振动特征频率,c表示振动管件平衡状态下与支撑件或防振件之间的间隙;t表示时间;ρ是间隙流体的密度,ν是间隙流体的运动学粘度;p(r,θ,t)表示压强;
无量纲边界条件为:
v(α,0,T)=v(α,π,T)=0
9.基于无限长直同心圆柱的挤压膜力计算系统,其特征在于,用于实现权利要求1-8任一项所述的基于无限长直同心圆柱的挤压膜力计算方法,包括:
预处理模块,用于将小间隙支撑结构简化为无限长直同心圆柱;
建模模块,用于基于无限长直同心圆柱进行流体控制方程量级分析,获得简化数学模型;对简化数学模型进行近似求解获得挤压膜力模型;
计算模块,基于挤压膜力模型计算出小间隙结构的挤压膜力。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,
所述计算机程序被处理器执行时能实现根据权利要求1至8中任意一项所述的基于无限长直同心圆柱的挤压膜力计算方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211623167.4A CN115796079B (zh) | 2022-12-16 | 2022-12-16 | 基于无限长直同心圆柱的挤压膜力计算方法、系统及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211623167.4A CN115796079B (zh) | 2022-12-16 | 2022-12-16 | 基于无限长直同心圆柱的挤压膜力计算方法、系统及介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115796079A true CN115796079A (zh) | 2023-03-14 |
CN115796079B CN115796079B (zh) | 2023-11-17 |
Family
ID=85425382
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211623167.4A Active CN115796079B (zh) | 2022-12-16 | 2022-12-16 | 基于无限长直同心圆柱的挤压膜力计算方法、系统及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115796079B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111859752A (zh) * | 2020-07-18 | 2020-10-30 | 西安交通大学 | 一种核反应堆蒸汽发生器流致振动计算方法 |
WO2022101483A1 (en) * | 2020-11-15 | 2022-05-19 | Qu & Co R&D B.V. | Solving a set of (non)linear differential equations using a hybrid data processing system comprising a classical computer system and a quantum computer system |
CN114818199A (zh) * | 2022-05-17 | 2022-07-29 | 哈尔滨工程大学 | 动载荷作用下考虑轴颈惯性力轴承水膜压力分布计算方法 |
-
2022
- 2022-12-16 CN CN202211623167.4A patent/CN115796079B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111859752A (zh) * | 2020-07-18 | 2020-10-30 | 西安交通大学 | 一种核反应堆蒸汽发生器流致振动计算方法 |
WO2022101483A1 (en) * | 2020-11-15 | 2022-05-19 | Qu & Co R&D B.V. | Solving a set of (non)linear differential equations using a hybrid data processing system comprising a classical computer system and a quantum computer system |
CN114818199A (zh) * | 2022-05-17 | 2022-07-29 | 哈尔滨工程大学 | 动载荷作用下考虑轴颈惯性力轴承水膜压力分布计算方法 |
Non-Patent Citations (3)
Title |
---|
EIRINI VELALI: "Model-Based Performance Analysis of Filtration Devices and Membrane Adsorbers", 《AVT》 * |
王晓娟;韩知非;丁千;张微;: "流体惯性对弹性环式挤压油膜阻尼器-转子系统动力特性的影响", 航空动力学报, no. 12 * |
费建国;罗会信;左兵权;魏志鹏;杨洋;徐海洋;李忠祥;: "雷诺方程的数值计算方法概述", 润滑与密封, no. 04 * |
Also Published As
Publication number | Publication date |
---|---|
CN115796079B (zh) | 2023-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Paı | The canonical problem of the fluid-conveying pipe and radiation of the knowledge gained to other dynamics problems across applied mechanics | |
Hassan et al. | Modeling of streamwise and transverse fluidelastic instability in tube arrays | |
Hassan et al. | Pitch and mass ratio effects on transverse and streamwise fluidelastic instability in parallel triangular tube arrays | |
Lu et al. | Study on the effect of reciprocating pump pipeline system vibration on oil transportation stations | |
Ma et al. | Numerical investigation of trajectory and attitude robustness of an underwater vehicle considering the uncertainty of platform velocity and yaw angle | |
Wang et al. | An integrated thermal and mechanical performance analysis of effect of cold molten salt temperature for thermocline tank | |
Olala et al. | Prediction of streamwise fluidelastic instability of a tube array in two-phase flow and effect of frequency detuning | |
Ji et al. | Research on vibration and heat transfer in heat exchanger with vortex generator | |
CN115796079A (zh) | 基于无限长直同心圆柱的挤压膜力计算方法、系统及介质 | |
Dazhi et al. | Numerical simulation analysis of liquid transportation in capsule-type vane tank under microgravity | |
López-Queija et al. | A simplified modeling approach of floating offshore wind turbines for dynamic simulations | |
Habib et al. | Internal flow effects in OTEC cold water pipe: Finite element modelling in frequency and time domain approaches | |
Fadaei et al. | A survey of numerical simulation tools for offshore wind turbine systems | |
CN113673130B (zh) | 蒸发器管束两相流功率谱密度相关长度的获取方法及系统 | |
CN115795715A (zh) | 一种用于高温气冷堆换热装置热工水力的仿真方法及系统 | |
Kang et al. | Analytical solution for a vibrating simply-supported cylinder subjected to 2-D concentric annular flow, considering friction | |
CN116127629B (zh) | 基于量级分析的小间隙结构挤压膜力的计算方法及系统 | |
Xu et al. | Performance of the compliant foil gas seal with surface micro-textured top foil | |
Lin et al. | Non-Linear coupled dynamics of a flexible propeller-shaft system supported by water film bearings | |
Li et al. | Modal analysis of two typical fluid-filled pipes in aircraft | |
Gao et al. | Free vibration of rubber matrix cord-reinforced combined shells of revolution under hydrostatic pressure | |
Dong et al. | Modeling and control of fluid flow networks with application to a nuclear-solar hybrid plant | |
Diao et al. | Simulation analysis of cavitation effect of hydraulic mechanism and its influence research | |
Weisensee et al. | Springboard droplet bouncing on flexible superhydrophobic substrates | |
Lin et al. | Analysis of nonlinear time-domain lubrication characteristics of the hydrodynamic journal bearing system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |