CN106294927B - 一种载流管道动力特性的计算方法 - Google Patents
一种载流管道动力特性的计算方法 Download PDFInfo
- Publication number
- CN106294927B CN106294927B CN201610575548.8A CN201610575548A CN106294927B CN 106294927 B CN106294927 B CN 106294927B CN 201610575548 A CN201610575548 A CN 201610575548A CN 106294927 B CN106294927 B CN 106294927B
- Authority
- CN
- China
- Prior art keywords
- formula
- pipeline
- liquid
- pipe
- filled pipe
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明涉及工程计算领域,尤其涉及一种载流管道动力特性的计算方法。本发明正确的描述了载流管道内流流速对管道横向弯曲振动固有频率的影响,提出了正确的计算方法,降低了管道在使用中存在的安全风险。
Description
技术领域
本发明涉及工程计算领域,尤其涉及一种载流管道动力特性的计算方法。
背景技术
管道作为物料输送的一种特种设备在现代化工业生产和人民生活中起着很重要的作用,管道事故时有发生,严重影响着人民的生命和财产安全。管道振动对安全生产造成很大的威胁。强烈的管道振动会使管路附件,特别是管道的连接部位和管道附件的连接部位等处发生松动和破裂,造成严重事故。
现阶段,对管道内的流体流动会影响管道的横向弯曲动力特性——固有频率进行计算时,得出的计算结果显示出,随着流速的增大管道的横向振动固有频率降低,但是与实际的情况不相符,在试验中,一根弯曲的软管会随着水流流速的增大而伸直,即管道的固有频率随流速的增大而提高。
发明内容
本发明要解决的技术问题是如何克服现有技术的不足,提供一种载流管道动力特性的计算方法。
本发明为实现上述目的采用的技术方案是:一种载流管道动力特性的计算方法,包括如下计算步骤:
(1)对管道微元段进行动力学分析,得出管道微元段的动力学方程,即公式1
式中:
m——管道单位长度的质量;
y——管道的横向变形(挠度);
x——管道长度方向的坐标;
ρ——管道内流体密度;
A——管道内壁横截面积;
v——管道内流流速;
EI——管道横截面抗弯模量;
t——时间;
(2)对公式1进行求解,由于公式1是建立在小变形假定基础上的,截面转动的角速度是高阶小量,因此,忽略科氏加速度后得公知振型方程,即公式2:
φⅣ(x)-αφ″(x)-βφ(x)=0
式中:
φ(x)——振型函数,
ω——载流管道横向弯曲振动固有频率;
(3)对公式2进行求解,设φ(x)=epx,代入公式2得关于p的4次方程:
p4-αp2-β=0
令p2=s,则上式为s的一元二次方程
s2-αs-β=0
由上式求得
则公式2的4个解为:
令
由此得
φ(x)=C1eδx+C2e-δx+C3eiεx+C4e-iεx
由欧拉公式
得到振型函数,即公式3:
φ(x)=D1chδx+D2shδx+D3cosεx+D4sinεx
式中:
D1、D2、D3、D4——待定系数;
δ、ε——管道两端约束条件有关的参数,
管道两端的约束条件包括:两端简支、两端固定、一端固定、一端固定一端简支和一端固定一端自由。
例如:两端简支条件为挠度为零和弯矩为零,即φ(0)=φ(l)=0,φ″(0)=φ″(l)=0,将它们代入上式即可求出δ,ε,而D1、D2、D3、D4这些均为结构动力学理论中的公知技术在上述求解过程中同时得到;
(4)根据管道两端实际约束条件,计算得出ε参数值;
(5)根据管道两端的约束条件有公式3求得关于ω的公式,即公式4:
(6)针对公式4中载流管道的EI、m、A进行必要的测算,并且根据实际情况测算出参数ρ、v的值,代入公式4中得出载流管道横向弯曲振动固有频率,即ω值。
本发明正确的描述了载流管道内流流速对管道横向弯曲振动固有频率的影响,提出了正确的计算方法,降低了管道在使用中存在的安全风险。
附图说明
图1本发明试验方法的试验结果图示。
图2本发明管道微元段的受力和运动示意图,
其中,T——张力,N、dN——剪力和剪力增量,M、dM——弯矩和弯矩增量,
dθ——管道微元段截面转角,——管道微元段出口和出口处的流体动量,——管道微元段横向加速度,dx——管道微元段长度
具体实施方式
建立正确的载流管道动力特性计算方法,能够正确的反应载流管道内的流体流动对管道横向弯曲动力特性(固有频率)的影响。
目前的载流管道动力特性是基于下列方程计算的:
式中:
m——管道单位长度的质量;
y——管道的横向变形(挠度);
x——管道长度方向的坐标;
ρ——管道内流体密度;
A——管道内壁横截面积;
v——管道内流流速;
由于管道的振动属于小变形,其截面转动的角速度非常小,因此,科氏加速度项此时,方程1简化为:
并计算得出载流管道的固有频率计算公式:
式中:ε——与管道两端约束条件有关的参数,管道两端的约束条件为
简支边界条件时,
从上述的固有频率计算公式可以看出,随着流速的增大管道的横向振动固有频率降低,但这与自然现象不符,也与试验现象不符;一根弯曲的软管会随着水流流速增大而伸直,如消防水龙带。图1中是本发明试验方法的试验结果,其表征了载流管道的固有频率随流速的增大而提高。
综上,目前的载流管道横向弯曲振动固有频率的计算方法是不正确的,原因在“方程1”不能正确的描述事实现象。
载流管道发生弯曲振动时,由于管道的弯曲变形,使得管内流体的动量发生变化,对于直管,尽管流速变化引起动量的变化,但由于动量的方向与管道轴线平行,因此,不会使管道产生横向作用力。但对于弯曲的管道或管道发生横向变形时,即使流速为常数(定长流速),由于流速方向的变化也将引起流体动量发生变化,而这部分变化的动量其方向垂直于流速(流速方向变化引起的动量增量垂直于流速),即垂直于管道,从而产生的横向作用力。
由于动量的增量等于力的冲量,因此,这部分动量的增量是管道的张力冲量在动量增量方向的投影。图2中示出了管道微元段的受力和运动示意图,对其进行分析,由此可得到管道微元段的动力学方程,即公式1
式中:
m——管道单位长度的质量;
y——管道的横向变形(挠度);
x——管道长度方向的坐标;
ρ——管道内流体密度;
A——管道内壁横截面积;
v——管道内流流速;
EI——管道横截面抗弯模量;
t——时间;
(2)对公式1进行求解,由于公式1是建立在小变形假定基础上的,截面转动的角速度是高阶小量,因此,忽略科氏加速度后得公知振型方程,即公式2:
φⅣ(x)-αφ″(x)-βφ(x)=0
式中:
φ(x)——振型函数,
ω——载流管道横向弯曲振动固有频率;
(3)对公式2进行求解,设φ(x)=epx,代入公式2得关于p的4次方程:
p4-αp2-β=0
令p2=s,则上式为s的一元二次方程
s2-αs-β=0
由上式求得
则公式2的4个解为:
令
由此得
φ(x)=C1eδx+C2e-δx+C3eiεx+C4e-iεx
由欧拉公式
得到振型函数,即公式3:
φ(x)=D1chδx+D2shδx+D3cosεx+D4sinεx
式中:
D1、D2、D3、D4——待定系数;
δ、ε——管道两端约束条件有关的参数,
管道两端的约束条件包括:两端简支、两端固定、一端固定、一端固定一端简支和一端固定一端自由。
其中一种约束条件:两端简支条件为挠度为零和弯矩为零,即φ(0)=φ(l)=0、φ″(0)=φ″(l)=0,将它们代入上式即可求出δ,ε,而D1、D2、D3、D4这些均为结构动力学理论中的公知技术在上述求解过程中同时得到;
(4)根据管道两端实际约束条件,计算得出ε参数值;
(5)根据管道两端的约束条件有公式3求得关于ω的公式,即公式4:
(6)针对公式4中载流管道的EI、m、A进行必要的测算,并且根据实际情况测算出参数ρ、v的值,代入公式4中得出载流管道横向弯曲振动固有频率,即ω值。
上述实施例只是为了说明本发明的技术构思及特点,其目的是在于让本领域内的普通技术人员能够了解本发明的内容并据以实施,并不能以此限制本发明的保护范围。凡是根据本发明内容的实质所作出的等效的变化或修饰,都应涵盖在本发明的保护范围内。
Claims (1)
1.一种载流管道动力特性的计算方法,其特征在于,包括如下计算步骤:
(1)对管道微元段进行动力学分析,得出管道微元段的动力学方程,即公式1
式中:
m——管道单位长度的质量;
y——管道的横向变形;
x——管道长度方向的坐标;
ρ——管道内流体密度;
A——管道内壁横截面积;
v——管道内流流速;
EI——管道横截面抗弯模量;
t——时间;
(2)对公式1进行求解,由于公式1是建立在小变形假定基础上的,截面转动的角速度是高阶小量,因此,忽略科氏加速度后得公知振型方程,即公式2:
φⅣ(x)-αφ″(x)-βφ(x)=0
式中:
φ(x)——振型函数,
ω——载流管道横向弯曲振动固有频率;
(3)对公式2进行求解,设φ(x)=epx,代入公式2得关于p的4次方程:
p4-αp2-β=0
令p2=s,则上式为s的一元二次方程
s2-αs-β=0
由上式求得
则公式2的4个解为:
令
由此得
φ(x)=C1eδx+C2e-δx+C3eiεx+C4e-iεx
由欧拉公式
得到振型函数,即公式3:
φ(x)=D1chδx+D2shδx+D3cosεx+D4sinεx
式中:
D1、D2、D3、D4——待定系数;
δ、ε——管道两端约束条件有关的参数,
(4)根据管道两端实际约束条件,计算得出ε参数值;
(5)根据管道两端的约束条件由公式3求得关于ω的公式,即公式4:
(6)针对公式4中载流管道的EI、m、A进行必要的测算,并且根据实际情况测算出参数ρ、v的值,代入公式4中得出载流管道横向弯曲振动固有频率,即ω值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610575548.8A CN106294927B (zh) | 2016-07-21 | 2016-07-21 | 一种载流管道动力特性的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610575548.8A CN106294927B (zh) | 2016-07-21 | 2016-07-21 | 一种载流管道动力特性的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106294927A CN106294927A (zh) | 2017-01-04 |
CN106294927B true CN106294927B (zh) | 2019-05-10 |
Family
ID=57651871
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610575548.8A Expired - Fee Related CN106294927B (zh) | 2016-07-21 | 2016-07-21 | 一种载流管道动力特性的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106294927B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110378060B (zh) * | 2019-07-26 | 2021-02-09 | 中国海洋大学 | 一种顶张式立管随机耦合振动的计算方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1876428A1 (de) * | 2006-07-04 | 2008-01-09 | Krohne AG | Verfahren zum Betreiben eines Coriolis-Massendurchflussmessgerätes |
CN105205280A (zh) * | 2015-10-10 | 2015-12-30 | 天津大学 | 一种单跨管道横流向涡激振动的预报方法 |
-
2016
- 2016-07-21 CN CN201610575548.8A patent/CN106294927B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1876428A1 (de) * | 2006-07-04 | 2008-01-09 | Krohne AG | Verfahren zum Betreiben eines Coriolis-Massendurchflussmessgerätes |
CN105205280A (zh) * | 2015-10-10 | 2015-12-30 | 天津大学 | 一种单跨管道横流向涡激振动的预报方法 |
Non-Patent Citations (2)
Title |
---|
充液管道动力学建模与振动特性分析;姚煜中;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20120115;C034-14 |
管道流固耦合振动的行波方法研究;任建亭等;《应用力学学报》;20050930;第22卷(第4期);530-536 |
Also Published As
Publication number | Publication date |
---|---|
CN106294927A (zh) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2386834B1 (en) | Multiphase flow measurement using a swirl inducer followed by a flow constriction | |
US7987733B2 (en) | Determination of density for metering a fluid flow | |
Liu et al. | Investigation on aerodynamic force nonlinear evolution for a central-slotted box girder under torsional vortex-induced vibration | |
CN101189508A (zh) | 用于确定管道粗糙度的方法和超声仪*** | |
Wang et al. | Pressure effects on Coriolis mass flowmeters | |
CN106934211A (zh) | 一种斜拉索索力与线形计算方法及应用 | |
CN105510627B (zh) | 用于测量由流体构成的流的参数的流动测量装置 | |
Chen et al. | Unsteady aerodynamics on a tapered prism under forced excitation | |
Djenidi et al. | Complete self-preservation on the axis of a turbulent round jet | |
CN106294927B (zh) | 一种载流管道动力特性的计算方法 | |
Rongmo et al. | Fluid-structure coupling analysis and simulation of viscosity effect on Coriolis mass flowmeter | |
Bi et al. | Numerical simulation and analysis of the effects of water-film morphological changes on the aerodynamic lift of stay cables | |
CN102353516B (zh) | 一种深水立管的时域脉动拖曳力确定方法 | |
CN102222168A (zh) | 一种深水钻井立管参激横向振动分析方法 | |
Chen et al. | A coupled model for vortex induced vibration of a circular cylinder with and without passive-jet flow control | |
US20090320552A1 (en) | Parallel flow meter device for measuring flow rate in pipes | |
CN102507082B (zh) | 一种深水立管的时域涡激升力确定方法 | |
US8919185B2 (en) | System and method for swirl generation | |
Ghozlani et al. | Numerical study of flow around an oscillating diamond prism and circular cylinder at low Keulegan-Carpenter number | |
CN102507083B (zh) | 一种尾流立管的脉动拖曳力确定方法 | |
Rongmo et al. | Experimental measurement and CFD simulation for characterization of Coriolis flowmeter performance | |
CN112834164A (zh) | 一种考虑约化速度和间距的尾流立管涡激升力确定方法 | |
Wen et al. | Supersonic separators for natural gas processing: Real gas effects | |
Miwa et al. | Two-phase flow induced force fluctuations on pipe bend | |
Okafor | Finite element analysis of vortex induced responses of multistory rectangular building |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190510 Termination date: 20200721 |