CN110110408A - 刚性圆柱体横流与顺流方向涡激振动耦合响应预测方法 - Google Patents
刚性圆柱体横流与顺流方向涡激振动耦合响应预测方法 Download PDFInfo
- Publication number
- CN110110408A CN110110408A CN201910337055.4A CN201910337055A CN110110408A CN 110110408 A CN110110408 A CN 110110408A CN 201910337055 A CN201910337055 A CN 201910337055A CN 110110408 A CN110110408 A CN 110110408A
- Authority
- CN
- China
- Prior art keywords
- formula
- follows
- vibration
- coefficient
- vortex
- 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
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000008878 coupling Effects 0.000 title claims abstract description 26
- 238000010168 coupling process Methods 0.000 title claims abstract description 26
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 26
- 230000004044 response Effects 0.000 title claims abstract description 26
- 238000006073 displacement reaction Methods 0.000 claims abstract description 19
- 239000012530 fluid Substances 0.000 claims description 27
- 238000013016 damping Methods 0.000 claims description 12
- 238000002474 experimental method Methods 0.000 claims description 11
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 6
- 230000010355 oscillation Effects 0.000 claims description 6
- 238000005381 potential energy Methods 0.000 claims description 6
- 239000004576 sand Substances 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 4
- 239000007788 liquid Substances 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 238000010200 validation analysis Methods 0.000 claims description 3
- 238000012795 verification Methods 0.000 claims description 2
- 238000010276 construction Methods 0.000 abstract description 3
- 238000011160 research Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 3
- 230000000737 periodic effect Effects 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 2
- 239000010931 gold Substances 0.000 description 2
- 229910052737 gold Inorganic materials 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种刚性圆柱体横流与顺流方向涡激振动耦合响应预测方法，主要包括以下步骤：1)建立横流方向与顺流方向互为耦合的振动控制方程；2)基于有限差分法对耦合的振动控制方程进行求解；3)基于分析数据，对实例进行计算分析。本发明针对刚性圆柱体横流与顺流方向耦合振动响应进行了研究。建立了完整的横流、顺流方向的耦合振动模型，用于分析同时考虑横流方向振动、顺流方向振动以及结构几何非线性的刚性圆柱体涡激振动响应预测问题。该模型以很好地模拟出结构的锁定以及位移突变现象。
Description
技术领域
本发明涉及圆柱体涡激振动响应领域，具体地说，特别涉及一种刚性圆柱体横流与顺流方向涡激振动耦合响应预测方法。
背景技术
圆柱体广泛存在于各种实际工程中，圆柱体在一定的来流中，会在结构两侧形成交替脱落的漩涡，周期性的漩涡脱落会在结构上产生周期性的拖曳力以及升力。若结构为弹性支撑，周期性的拖曳力和升力会引起结构在顺流以及横流方向产生周期性的振动，统称为“涡激振动”(Gao,Y,Fu,SX,Xiong,YM,Zhao,Y,and Liu,LM,Experimental study onresponse performance of vortex-induced vibration of a flexible cylinder,Shipsand Offshore Structures,2017,12(1):116-134)。当漩涡脱落频率出现在结构固有频率附近时，便会发生锁定现象。锁定区域内，圆柱体会发生大幅的、危险的涡激振动响应，这种响应会对结构带来很大的疲劳损伤。
现有的研究圆柱体涡激振动响应的方法主要可分为实验方法以及数值方法。与数值方法相比，实验方法得到的数据更为可靠、得到的现象更为直观，但实验方法通常研究成本较高。数值方法通常主要包括计算流体动力学(Computational Fluid Dynamics,CFD)方法(Gao,Y,Zong,Z,Zou,L,Takagi,S,Jiang,ZY,Numerical simulation of vortex-induced vibration of a circular cylinder with different surface roughnesses,Marine Structures,2018,57:165-179)以及经验模型方法。与经验模型方法相比，CFD方法的计算精度要高，但CFD方法的数值计算成本要比经验模型方法高很多。为了适应海洋工程中的多工况问题，有必要建立一种能够快速预报结构涡激振动响应主要特性的经验模型方法，而尾流振子模型方法则正是近些年广泛应用于涡激振动响应预报中的一种经验模型方法。
目前基于尾流振子模型方法对刚性圆柱体涡激振动响应的研究主要集中只考虑结构在横流方向的振动，而忽略了结构在顺流方向的振动。而实际工程中横流方向与顺流方向振动互为耦合、互相影响。在少量的同时考虑横流以及顺流方向振动响应的研究中，针对结构几何非线性对涡激振动响应的影响研究则更少。对于这种同时考虑结构几何非线性以及CF方向与IL方向互为耦合的圆柱体模型，其涡激振动响应特性有待进一步展开研究。
发明内容
为了解决现有技术的问题，本发明实施例提供了一种刚性圆柱体横流与顺流方向涡激振动耦合响应预测方法。所述技术方案如下：
本发明提供了一种刚性圆柱体横流与顺流方向涡激振动耦合响应预测方法，主要包括以下步骤：
1)建立横流方向与顺流方向互为耦合的振动控制方程；
2)基于有限差分法对耦合的振动控制方程进行求解；
3)基于分析数据，对实例进行计算分析。
进一步的，所述步骤1)具体包括：
a、建立均匀来流下刚性圆柱体二自由度涡激振动耦合模型；
b、在所述均匀来流下刚性圆柱体二自由度涡激振动耦合模型中，考虑一单位长度、直径为D的刚性圆柱体在均匀来流U作用下引起的CF方向以及IL方向互为耦合的涡激振动响应问题，其中X方向为IL方向，Y方向为CF方向。
其中，在分析过程中假设系统在X方向和Y方向呈对称分布，因此在X方向的质量、刚度、阻尼以及固有频率均与Y方向对应的参数相等；
所述均匀来流下刚性圆柱体二自由度涡激振动耦合模型中设有4弹簧刚性圆柱体系统。
进一步的，所述4弹簧刚性圆柱体系统在X和Y方向的振动方程描述如下：
式(1)d表示对变量求全导；表示对变量求偏导；T为时间；为拉格朗日函数，表示为：其中K为动能，P为势能；u1以及u2分别为X方向和Y方向位移；和分别为X方向和Y方向速度；F1和F2分别表示作用在结构上的X方向以及Y方向水动力载荷Fx以及Fy；r为阻尼系数，包括结构阻尼系数rs以及流体阻尼系数rf两项，流体阻尼系数rf表示为：rf＝ωfγρD2＝(2πStU/D)γρD2，其中ωf为漩涡脱落频率；St为斯脱哈尔数，这里取为0.2；ρ为流体密度；γ为粘滞力系数，取为0.8。
进一步的，假设与圆柱体相连的4根弹簧初始长度均为S，弹簧弹性系数均为k，并假设圆柱体T时刻在X方向以及Y方向的结构位移分别为X(T)以及Y(T)，便得到此时刻圆柱体的动能K以及势能P，表示如下：
式(2)中m为单位长度系统总质量，由结构质量ms以及流体附加质量mf两部分构成，流体附加质量表示为：mf＝CMρD2π/4，其中CM为附加质量系数，对于圆柱体取1.0；将式(2)和式(3)代入式(1)中，得到圆柱体在X方向和Y方向的振动方程为:
方程式(4)和(5)中的非线性项体现了结构振动的几何非线性特征。为了求解式(4)和式(5)，需要对二式进行一定的简化处理；函数f(X,Y)在(X0,Y0)的二元泰勒级数展开式表示如下：
将式(4)和式(5)中的非线性项按照式(6)在(0,0)处进行泰勒级数展开，取前三阶得到：
将式(7)代入式(4)、式(8)代入(5)中得到：
式(9)中h＝2k，为了方便表达，记：式(9)进一步写作：
式(10)中，h为系统刚度，已知系统总质量m以及系统刚度h，便可得到系统在X方向以及Y方向的固有频率，表示为：式(10)中，X3以及Y3前面的系数以及反应了结构在轴向变形的几何非线性特性；XY2以及YX2前面的系数以及反应了CF以及IL方向振动互为耦合的的几何非线性特性；本发明将这以及这4个参数统称为结构的几何非线性系数。
进一步的，对于静止圆柱体，作用在圆柱体上的Fx以及Fy的方向分别与拖曳力FD以及升力FL方向一致；但对于振荡圆柱体，由于结构的振动，导致结构与流体之间的相对速度方向不再沿X方向，使得：Fx将不再与FD方向一致，Fy也将不再与FL方向一致。
进一步的，假设结构的振动速度V方向与流速U方向的夹角为θ，通常情况下振动速度V与流速U比起来可以看作是小量，因此夹角θ也为小量；此时Fx以及Fy表示如下：
式(11)中为平均拖曳力，FD为振荡拖曳力，FL为升力，分别表示如下：
式(12)中，CD以及CL分别表示平均拖曳力系数、振荡拖曳力系数以及升力系数，这里平均拖曳力系数取为1.2，CD和CL分别写作：
式(13)中，CD0与CL0分别表示静止圆柱体上的振荡拖曳力系数以及升力系数，分别取为0.2以及0.3；p(T)以及q(T)分别表示在X以及Y方向的尾流振子的运动；
采用Van der Pol方程来描述尾流振子特性，表示如下：
式(14)中εx、εy、Ax以及Ay分别为经验参数，Ax和Ay取为12，εx取为0.3，而εy则参照数值与实验的对比结果，对预设数据进行拟合，得到εy关于m*的表示式如下：
将方程(10)和(14)转换为无量纲形式，引入IL以及CF方向的无量纲位移x、y以及无量纲时间t，如下：
x＝X/D,y＝Y/D,t＝T·ωn (16)
将式(16)带入方程式(10)和(14)中，整理得到X方向以及Y方向结构和尾流振子的无量纲方程为：
式(17)中，αx、βx、αy以及βy为结构的几何非线性无量纲系数，ξ为结构阻尼比，与结构阻尼系数rs之间的关系为：ξ＝rs/2mωn；式(17)中频率比Ω、质量比μ、以及系统无量纲参数MD以及ML分别表示如下：
式(17)以及式(18)中Ur为折合速度：Ur＝2πU/ωnD，在分析中使用的质量比为m*，即：结构质量与被结构排开的液体质量的比值，其与μ之间的关系为：m*＝4μ/π-CM。
进一步的，所述步骤2)具体如下：
在时间上采用二阶精度中心差分格式对方程式(17)进行先离散后迭代求解；将计算无量纲总时间ttotal划分为N段，得到每个时间步长Δt＝ttotal/N，被离散后的N+1个时刻依次记为：t＝tj(j＝0,1,2,..,N)；假设tn时刻对应的无量纲参数x、y、p以及q分别表示为xn、yn、pn以及qn；那么tn时刻式(17)中各导数项的二阶精度中心差分格式表示为：
将式(19)代入式(17)整理得到：
x和y的初始条件(t＝t0)设为位移以及速度均为0，即：x0＝y0＝0以及dx/dt＝dy/dt＝0；p和q的初始条件(t＝t0)为：p0以及q0均取一微小的振幅以及dp/dt＝dq/dt＝0；由dx/dt＝dy/dt＝dp/dt＝dq/dt＝0得到：
x-1＝x1,y-1＝y1,p-1＝p1,q-1＝q1 (21)
将式(21)代入式(20)得到t1时刻x、y、p以及q的值，表示如下：
依据式(22)得到y1、q1、x1以及p1；至此得到了t0以及t1时刻的x、y、p以及q的值；在此基础上，依据式(20)通过迭代依次求解得到t2,t3,…,tN时刻的x、y、p以及q的值。
进一步的，所述步骤3)具体如下：
对尾流振子数值模型进行数值验证；在验证过程中，选取的计算参数m*＝2.36,ζ＝0.006，计算总时间ttotal取为1000个无量纲时间，时间步长Δt取为0.1，折合速度Ur区间取为1.0-14.0，区间间隔ΔUr取为0.1。
本发明实施例提供的技术方案带来的有益效果是：
本发明针对刚性圆柱体横流与顺流方向耦合振动响应进行了研究。建立了完整的横流、顺流方向的耦合振动模型，用于分析同时考虑横流方向振动、顺流方向振动以及结构几何非线性的刚性圆柱体涡激振动响应预测问题。该模型可以很好地模拟出结构的锁定以及位移突变现象。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例的均匀来流下刚性圆柱体二自由度涡激振动耦合模型的示意图；
图2是本发明实施例的振动圆柱体受力示意图；
图3是本发明实施例的数值模型CF与IL方向振动幅值与Stappenbelt等人的实验结果对比示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
本发明提供了一种刚性圆柱体横流与顺流方向涡激振动耦合响应预测方法，主要包括以下步骤：
1)建立横流方向与顺流方向互为耦合的振动控制方程
2)基于有限差分法对耦合的振动控制方程进行求解；
3)基于分析数据，对实例进行计算分析。
步骤1：建立结构和流场互为耦合的振动控制方程，具体如下：
如图1所示，考虑一单位长度、直径为D的刚性圆柱体在均匀来流U作用下引起的CF方向以及IL方向互为耦合的涡激振动响应问题，其中X方向为IL方向，Y方向为CF方向。本发明在分析过程中假设系统在X方向和Y方向呈对称分布，因此在X方向的质量、刚度、阻尼以及固有频率均与Y方向对应的参数相等。
依据Hamilton原理，图1所示的4弹簧刚性圆柱体系统在X和Y方向的振动方程可描述如下：
式(1)中T为时间；为拉格朗日函数，可表示为：其中K为动能，P为势能；u1以及u2分别为X方向和Y方向位移；和分别为X方向和Y方向速度；F1和F2分别表示作用在结构上的X方向以及Y方向水动力载荷Fx以及Fy；r为阻尼系数，包括结构阻尼系数rs以及流体阻尼系数rf两项，流体阻尼系数rf可表示为：rf＝ωfγρD2＝(2πStU/D)γρD2，其中ωf为漩涡脱落频率；St为斯脱哈尔数，这里取为0.2；ρ为流体密度；γ为粘滞力系数，可取为0.8(Facchinetti,ML,de Langre,E,Biolley,F,Coupling of structure and wakeoscillators in vortex-induced vibrations,Journal of Fluids and Structures,2004,19:123-140)。如图1所示，假设与圆柱体相连的4根弹簧初始长度均为S，弹簧弹性系数均为k，并假设圆柱体T时刻在X方向以及Y方向的结构位移分别为X(T)以及Y(T)，便可得到此时刻圆柱体的动能K以及势能P，表示如下：
式(2)中m为单位长度系统总质量，由结构质量ms以及流体附加质量mf两部分构成，流体附加质量可表示为：mf＝CMρD2π/4，其中CM为附加质量系数，对于圆柱体通常取1.0(Govardhan,R,Williamson,CHK,Modes of vortex formation and frequency responseof a freely vibrating cylinder,Journal of Fluid Mechanics,2000,420:85-130)。将式(2)和式(3)代入式(1)中，可以得到圆柱体在X方向和Y方向的振动方程为:
方程式(4)和(5)中的非线性项体现了结构振动的几何非线性特征。为了求解式(4)和式(5)，需要对二式进行一定的简化处理。函数f(X,Y)在(X0,Y0)的二元泰勒级数展开式可表示如下：
将式(4)和式(5)中的非线性项按照式(6)在(0,0)处进行泰勒级数展开，取前三阶得到：
将式(7)代入式(4)、式(8)代入(5)中得到：
式(9)中h＝2k，为了方便表达，记：式(9)可进一步写作：
式(10)中，h为系统刚度，已知系统总质量m以及系统刚度h，便可得到系统在X方向以及Y方向的固有频率，可表示为：式(10)中，X3以及Y3前面的系数以及反应了结构在轴向变形的几何非线性特性；XY2以及YX2前面的系数以及反应了CF以及IL方向振动互为耦合的的几何非线性特性；本发明将这以及这4个参数统称为结构的几何非线性系数。
对于静止圆柱体，作用在圆柱体上的Fx以及Fy的方向分别与拖曳力FD以及升力FL方向一致。但对于振荡圆柱体，由于结构的振动，导致结构与流体之间的相对速度方向不再沿X方向，使得：Fx将不再与FD方向一致，Fy也将不再与FL方向一致。
如图2所示，假设结构的振动速度V方向与流速U方向的夹角为θ，通常情况下振动速度V与流速U比起来可以看作是小量，因此夹角θ也为小量。此时Fx以及Fy可表示如下(Wang,XQ,So,RMC,Chan,KT,A non-linear fluid force model for vortex-inducedvibration of an elastic cylinder,Journal of Sound and Vibration,2003,260:287-305)：
式(11)中为平均拖曳力，FD为振荡拖曳力，FL为升力，分别表示如下：
式(12)中，CD以及CL分别表示平均拖曳力系数、振荡拖曳力系数以及升力系数，这里平均拖曳力系数取为1.2(Ge,F,Long,X,Wang,L,Hong,Y,Flow-inducedvibration of a long circular cylinders modelled by coupled nonlinearoscillators,Science in China Series G:Physics,Mechanics and Astronomy,2009,52(7):1086-1093)，CD和CL分别写作(Facchinetti,ML,de Langre,E,Biolley,F,Coupling ofstructure and wake oscillators in vortex-induced vibrations,Journal of Fluidsand Structures,2004,19:123-140)：
式(13)中，CD0与CL0分别表示静止圆柱体上的振荡拖曳力系数以及升力系数，分别取为0.2以及0.3(Srinil,N.,Zanganeh,H.Modelling of coupled cross-flow/in-linevortex-induced vibrations using double Duffing and van der Poloscillators.Ocean Engineering,2012；53:83-97)。p(T)以及q(T)分别表示在X以及Y方向的尾流振子的运动。
采用Van der Pol方程来描述尾流振子特性，可表示如下(Nayfeh,N.H.Introduction to Perturbation Techniques.New York:Wiley,1993)：
式(14)中εx、εy、Ax以及Ay分别为经验参数，Ax和Ay可取为12，εx可取为0.3(Facchinetti,ML,de Langre,E,Biolley,F,Coupling of structure and wakeoscillators in vortex-induced vibrations,Journal of Fluids and Structures,2004,19:123-140)，而εy则参考文献中(Stappenbelt,B.,Laiji,F.,Tan,G.Low massratio vortex-induced motion.In:the 16th Australian Fluid Mechanics Conference,Gold,Coast,Australia,2007,pp.1491-1497)表1(b)中数值与实验的对比结果，对表中数据进行拟合，可得到εy关于m*的表示式如下：
将方程(10)和(14)转换为无量纲形式，引入IL以及CF方向的无量纲位移x、y以及无量纲时间t，如下：
x＝X/D,y＝Y/D,t＝T·ωn (16)
将式(16)带入方程式(10)和(14)中，整理得到X方向以及Y方向结构和尾流振子的无量纲方程为：
式(17)中，αx、βx、αy以及βy为结构的几何非线性无量纲系数，ξ为结构阻尼比，与结构阻尼系数rs之间的关系为：ξ＝rs/2mωn。式(17)中频率比Ω、质量比μ、以及系统无量纲参数MD以及ML分别表示如下：
式(17)以及式(18)中Ur为折合速度：Ur＝2πU/ωnD，这里值得注意的是通常在分析中大多时使用的质量比为m*，即：结构质量与被结构排开的液体质量的比值，其与μ之间的关系为：m*＝4μ/π-CM。
步骤2：基于有限差分法对耦合的振动控制方程进行求解，具体如下：
在时间上采用二阶精度中心差分格式对方程式(17)进行先离散后迭代求解。将计算无量纲总时间ttotal划分为N段，得到每个时间步长Δt＝ttotal/N，被离散后的N+1个时刻依次记为：t＝tj(j＝0,1,2,..,N)。假设tn时刻对应的无量纲参数x、y、p以及q分别表示为xn、yn、pn以及qn。那么tn时刻式(17)中各导数项的二阶精度中心差分格式可表示为：
将式(19)代入式(17)整理得到：
x和y的初始条件(t＝t0)设为位移以及速度均为0，即：x0＝y0＝0以及dx/dt＝dy/dt＝0；p和q的初始条件(t＝t0)为：p0以及q0均取一微小的振幅以及dp/dt＝dq/dt＝0。由dx/dt＝dy/dt＝dp/dt＝dq/dt＝0得到：
x-1＝x1,y-1＝y1,p-1＝p1,q-1＝q1 (21)
将式(21)代入式(20)可得到t1时刻x、y、p以及q的值，表示如下：
依据式(22)可以此得到y1、q1、x1以及p1。至此得到了t0以及t1时刻的x、y、p以及q的值。在此基础上，依据式(20)可通过迭代依次求解得到t2,t3,…,tN时刻的x、y、p以及q的值。
步骤3：基于分析数据，对实例进行计算分析。
为了验证本发明中尾流振子数值模型的可靠性，这里对其进行了数值验证。在验证过程中，选取的计算参数与Stappenbelt等(2007)的实验参数相同(m*＝2.36,ξ＝0.006)(Stappenbelt,B.,Laiji,F.,Tan,G.Low mass ratio vortex-induced motion.In:the16th Australasian Fluid Mechanics Conference,Gold,Coast,Australia,2007,pp.1491-1497)，几何非线性无量纲系数αx、βx、αy以及βy均取为0.7(Srinil,N.,Zanganeh,H.Modelling of coupled cross-flow/in-line vortex-induced vibrations usingdouble Duffing and van der Pol oscillators.Ocean Engineering,2012；53:83-97)，计算总时间ttotal取为1000个无量纲时间，时间步长Δt取为0.1，折合速度Ur区间取为1.0-14.0，区间间隔ΔUr取为0.1。
图3给出了Stappenbelt等人实验方法以及本发明中尾流振子模型方法得到的CF以及IL方向的振动幅值。这里值得注意的是：由于本发明中数值模型得到的Ax/D包括两部分，即平均拖曳力引起的平均位移以及振荡拖曳力FD引起的振荡位移。因此，在后续分析中需要将引起的平均位移加以消除，只留下FD引起的振荡位移。由图3可以看出：本发明中数值方法以及他人实验方法得到的Ay/D以及Ax/D随折合速度Ur的变化趋势一致：随着Ur的增加，振动幅值明显增大，响应区间从初始分支(initial branch)进入上分支(upperbranch)；随后Ur的进一步增加，振动幅值出现突降，进入低分支(lower branch)。实验方法以及本发明中数值方法得到的不同折合速度下Ay/D以及Ax/D的最大值均非常接近：Ay/D最大值均在1.5附近；Ax/D最大值均在0.25附近。且两种方法得到的振动幅值最大值均出现在Ur＝8.5附近。总体上来说，本发明中尾流振子模型方法的数值计算结果与Stappenbelt等人的实验结果吻合较好，从而验证了本发明中数值模型的可靠性。
本发明实施例提供的技术方案带来的有益效果是：
本发明针对刚性圆柱体横流与顺流方向耦合振动响应进行了研究。建立了完整的横流、顺流方向的耦合振动模型，用于分析同时考虑横流方向振动、顺流方向振动以及结构几何非线性的刚性圆柱体涡激振动响应预测问题。该模型可以很好地模拟出结构的锁定以及位移突变现象。
以上仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种刚性圆柱体横流与顺流方向涡激振动耦合响应预测方法，其特征在于，主要包括以下步骤：
1)建立横流方向与顺流方向互为耦合的振动控制方程；
2)基于有限差分法对耦合的振动控制方程进行求解；
3)基于分析数据，对实例进行计算分析。
2.如权利要求1所述的方法，其特征在于，所述步骤1)具体包括：
a、建立均匀来流下刚性圆柱体二自由度涡激振动耦合模型；
b、在所述均匀来流下刚性圆柱体二自由度涡激振动耦合模型中，考虑一单位长度、直径为D的刚性圆柱体在均匀来流U作用下引起的CF方向以及IL方向互为耦合的涡激振动响应问题，其中X方向为IL方向，Y方向为CF方向；
其中，在分析过程中假设系统在X方向和Y方向呈对称分布，因此在X方向的质量、刚度、阻尼以及固有频率均与Y方向对应的参数相等；
所述均匀来流下刚性圆柱体二自由度涡激振动耦合模型中设有4弹簧刚性圆柱体系统。
3.如权利要求2所述的方法，其特征在于，所述4弹簧刚性圆柱体系统在X和Y方向的振动方程描述如下：
式(1)d表示对变量求全导；表示对变量求偏导；T为时间；为拉格朗日函数，表示为：其中K为动能，P为势能；u1以及u2分别为X方向和Y方向位移；和分别为X方向和Y方向速度；F1和F2分别表示作用在结构上的X方向以及Y方向水动力载荷Fx以及Fy；r为阻尼系数，包括结构阻尼系数rs以及流体阻尼系数rf两项，流体阻尼系数rf表示为：rf＝ωfγρD2＝(2πStU/D)γρD2，其中ωf为漩涡脱落频率；St为斯脱哈尔数，这里取为0.2；ρ为流体密度；γ为粘滞力系数，取为0.8。
4.如权利要求3所述的方法，其特征在于，假设与圆柱体相连的4根弹簧初始长度均为S，弹簧弹性系数均为k，并假设圆柱体T时刻在X方向以及Y方向的结构位移分别为X(T)以及Y(T)，便得到此时刻圆柱体的动能K以及势能P，表示如下：
式(2)中m为单位长度系统总质量，由结构质量ms以及流体附加质量mf两部分构成，流体附加质量表示为：mf＝CMρD2π/4，其中CM为附加质量系数，对于圆柱体取1.0；将式(2)和式(3)代入式(1)中，得到圆柱体在X方向和Y方向的振动方程为:
方程式(4)和(5)中的非线性项体现了结构振动的几何非线性特征；为了求解式(4)和式(5)，需要对二式进行一定的简化处理；函数f(X,Y)在(X0,Y0)的二元泰勒级数展开式表示如下：
将式(4)和式(5)中的非线性项按照式(6)在(0,0)处进行泰勒级数展开，取前三阶得到：
将式(7)代入式(4)、式(8)代入(5)中得到：
式(9)中h＝2k，为了方便表达，记：式(9)进一步写作：
式(10)中，h为系统刚度，已知系统总质量m以及系统刚度h，便得到系统在X方向以及Y方向的固有频率，表示为：式(10)中，X3以及Y3前面的系数以及反应了结构在轴向变形的几何非线性特性；XY2以及YX2前面的系数以及反应了CF以及IL方向振动互为耦合的的几何非线性特性；本发明将这以及这4个参数统称为结构的几何非线性系数。
5.如权利要求4所述的方法，其特征在于，对于静止圆柱体，作用在圆柱体上的Fx以及Fy的方向分别与拖曳力FD以及升力FL方向一致；但对于振荡圆柱体，由于结构的振动，导致结构与流体之间的相对速度方向不再沿X方向，使得：Fx将不再与FD方向一致，Fy也将不再与FL方向一致。
6.如权利要求5所述的方法，其特征在于，假设结构的振动速度V方向与流速U方向的夹角为θ，通常情况下振动速度V与流速U比起来可以看作是小量，因此夹角θ也为小量；此时Fx以及Fy表示如下：
式(11)中为平均拖曳力，FD为振荡拖曳力，FL为升力，分别表示如下：
式(12)中，CD以及CL分别表示平均拖曳力系数、振荡拖曳力系数以及升力系数，这里平均拖曳力系数取为1.2，CD和CL分别写作：
式(13)中，CD0与CL0分别表示静止圆柱体上的振荡拖曳力系数以及升力系数，分别取为0.2以及0.3；p(T)以及q(T)分别表示在X以及Y方向的尾流振子的运动；
采用Van der Pol方程来描述尾流振子特性，表示如下：
式(14)中εx、εy、Ax以及Ay分别为经验参数，Ax和Ay取为12，εx取为0.3，而εy则参照数值与实验的对比结果，对预设数据进行拟合，得到εy关于m*的表示式如下：
将方程(10)和(14)转换为无量纲形式，引入IL以及CF方向的无量纲位移x、y以及无量纲时间t，如下：
x＝X/D,y＝Y/D,t＝T·ωn (16)
将式(16)带入方程式(10)和(14)中，整理得到X方向以及Y方向结构和尾流振子的无量纲方程为：
式(17)中，αx、βx、αy以及βy为结构的几何非线性无量纲系数，ξ为结构阻尼比，与结构阻尼系数rs之间的关系为：ξ＝rs/2mωn；式(17)中频率比Ω、质量比μ、以及系统无量纲参数MD以及ML分别表示如下：
式(17)以及式(18)中Ur为折合速度：Ur＝2πU/ωnD，在分析中使用的质量比为m*，即：结构质量与被结构排开的液体质量的比值，其与μ之间的关系为：m*＝4μ/π-CM。
7.如权利要求6所述的方法，其特征在于，所述步骤2)具体如下：
在时间上采用二阶精度中心差分格式对方程式(17)进行先离散后迭代求解；将计算无量纲总时间ttotal划分为N段，得到每个时间步长Δt＝ttotal/N，被离散后的N+1个时刻依次记为：t＝tj(j＝0,1,2,..,N)；假设tn时刻对应的无量纲参数x、y、p以及q分别表示为xn、yn、pn以及qn；那么tn时刻式(17)中各导数项的二阶精度中心差分格式表示为：
将式(19)代入式(17)整理得到：
x和y的初始条件(t＝t0)设为位移以及速度均为0，即：x0＝y0＝0以及dx/dt＝dy/dt＝0；p和q的初始条件(t＝t0)为：p0以及q0均取一微小的振幅以及dp/dt＝dq/dt＝0；由dx/dt＝dy/dt＝dp/dt＝dq/dt＝0得到：
x-1＝x1,y-1＝y1,p-1＝p1,q-1＝q1 (21)
将式(21)代入式(20)得到t1时刻x、y、p以及q的值，表示如下：
依据式(22)得到y1、q1、x1以及p1；至此得到了t0以及t1时刻的x、y、p以及q的值；在此基础上，依据式(20)通过迭代依次求解得到t2,t3,…,tN时刻的x、
y、p以及q的值。
8.如权利要求7所述的方法，其特征在于，所述步骤3)具体如下：
对尾流振子数值模型进行数值验证；在验证过程中，选取的计算参数m*＝2.36,ζ＝0.006，计算总时间ttotal取为1000个无量纲时间，时间步长Δt取为0.1，折合速度Ur区间取为1.0-14.0，区间间隔ΔUr取为0.1。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910337055.4A CN110110408B (zh) | 2019-04-25 | 2019-04-25 | 刚性圆柱体横流与顺流方向涡激振动耦合响应预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910337055.4A CN110110408B (zh) | 2019-04-25 | 2019-04-25 | 刚性圆柱体横流与顺流方向涡激振动耦合响应预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110110408A true CN110110408A (zh) | 2019-08-09 |
| CN110110408B CN110110408B (zh) | 2022-06-07 |
Family
ID=67486511
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910337055.4A Expired - Fee Related CN110110408B (zh) | 2019-04-25 | 2019-04-25 | 刚性圆柱体横流与顺流方向涡激振动耦合响应预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110110408B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110598337A (zh) * | 2019-09-17 | 2019-12-20 | 中国海洋大学 | 一种圆柱体涡激振动的流固耦合时域分析方法 |
| CN110795827A (zh) * | 2019-10-11 | 2020-02-14 | 扬州大学 | 非线性能量阱作用下弹性支撑刚性柱体结构的涡激振动仿真方法 |
| CN111177926A (zh) * | 2019-12-30 | 2020-05-19 | 清华大学深圳国际研究生院 | 铺管过程中的涡激振动预报方法 |
| CN111985138A (zh) * | 2020-08-21 | 2020-11-24 | 哈尔滨工业大学（威海） | 一种柔性结构横流与顺流方向涡激振动耦合响应预测方法 |
| CN112084560A (zh) * | 2020-09-09 | 2020-12-15 | 长沙理工大学 | 一种多幅桥面桥梁涡激振动数值模拟计算方法 |
| CN113111420A (zh) * | 2021-04-19 | 2021-07-13 | 哈尔滨工业大学（威海） | 一种边界激励细长张力梁不稳定区间的快速预测方法 |
| CN115391881A (zh) * | 2022-08-09 | 2022-11-25 | 哈尔滨工业大学 | 一种桥塔尾流区吊索风致振动数值预测方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101539477A (zh) * | 2009-05-08 | 2009-09-23 | 中国海洋大学 | 一种深水顶张式立管涡激振动与疲劳分析的方法 |
| CN102507084A (zh) * | 2011-09-30 | 2012-06-20 | 中国海洋大学 | 一种尾流立管的时域涡激升力确定方法 |
| CN105205280A (zh) * | 2015-10-10 | 2015-12-30 | 天津大学 | 一种单跨管道横流向涡激振动的预报方法 |
| CN105279373A (zh) * | 2015-10-10 | 2016-01-27 | 天津大学 | 一种海底自由悬跨管道纯顺流向涡激振动预报方法 |
| CN108414191A (zh) * | 2018-02-09 | 2018-08-17 | 山东交通学院 | 一种考虑长径比影响的浮式圆柱涡激运动分析系统及方法 |
-
2019
- 2019-04-25 CN CN201910337055.4A patent/CN110110408B/zh not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101539477A (zh) * | 2009-05-08 | 2009-09-23 | 中国海洋大学 | 一种深水顶张式立管涡激振动与疲劳分析的方法 |
| CN102507084A (zh) * | 2011-09-30 | 2012-06-20 | 中国海洋大学 | 一种尾流立管的时域涡激升力确定方法 |
| CN105205280A (zh) * | 2015-10-10 | 2015-12-30 | 天津大学 | 一种单跨管道横流向涡激振动的预报方法 |
| CN105279373A (zh) * | 2015-10-10 | 2016-01-27 | 天津大学 | 一种海底自由悬跨管道纯顺流向涡激振动预报方法 |
| CN108414191A (zh) * | 2018-02-09 | 2018-08-17 | 山东交通学院 | 一种考虑长径比影响的浮式圆柱涡激运动分析系统及方法 |
Non-Patent Citations (4)
| Title |
|---|
| 唐友刚 等: "高雷诺数下圆柱顺流向和横向涡激振动分析", 《振动与冲击》 * |
| 唐友刚 等: "高雷诺数下圆柱顺流向和横向涡激振动分析", 《振动与冲击》, vol. 32, no. 13, 31 July 2013 (2013-07-31), pages 88 - 92 * |
| 高云 等: "基于有限差分法刚性圆柱体涡激振动响应特性数值研究", 《计算物理》 * |
| 高云 等: "基于有限差分法刚性圆柱体涡激振动响应特性数值研究", 《计算物理》, vol. 36, no. 1, 30 January 2018 (2018-01-30), pages 53 - 59 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110598337A (zh) * | 2019-09-17 | 2019-12-20 | 中国海洋大学 | 一种圆柱体涡激振动的流固耦合时域分析方法 |
| CN110795827A (zh) * | 2019-10-11 | 2020-02-14 | 扬州大学 | 非线性能量阱作用下弹性支撑刚性柱体结构的涡激振动仿真方法 |
| CN110795827B (zh) * | 2019-10-11 | 2024-02-06 | 扬州大学 | 非线性能量阱作用下弹性支撑刚性柱体结构的涡激振动仿真方法 |
| CN111177926B (zh) * | 2019-12-30 | 2022-06-21 | 清华大学深圳国际研究生院 | 铺管过程中的涡激振动预报方法 |
| CN111177926A (zh) * | 2019-12-30 | 2020-05-19 | 清华大学深圳国际研究生院 | 铺管过程中的涡激振动预报方法 |
| CN111985138A (zh) * | 2020-08-21 | 2020-11-24 | 哈尔滨工业大学（威海） | 一种柔性结构横流与顺流方向涡激振动耦合响应预测方法 |
| CN111985138B (zh) * | 2020-08-21 | 2023-12-19 | 哈尔滨工业大学（威海） | 一种柔性结构横流与顺流方向涡激振动耦合响应预测方法 |
| CN112084560B (zh) * | 2020-09-09 | 2023-03-24 | 长沙理工大学 | 一种多幅桥面桥梁涡激振动数值模拟计算方法 |
| CN112084560A (zh) * | 2020-09-09 | 2020-12-15 | 长沙理工大学 | 一种多幅桥面桥梁涡激振动数值模拟计算方法 |
| CN113111420B (zh) * | 2021-04-19 | 2022-04-22 | 哈尔滨工业大学（威海） | 一种边界激励细长张力梁不稳定区间的快速预测方法 |
| CN113111420A (zh) * | 2021-04-19 | 2021-07-13 | 哈尔滨工业大学（威海） | 一种边界激励细长张力梁不稳定区间的快速预测方法 |
| CN115391881A (zh) * | 2022-08-09 | 2022-11-25 | 哈尔滨工业大学 | 一种桥塔尾流区吊索风致振动数值预测方法 |
| CN115391881B (zh) * | 2022-08-09 | 2023-04-18 | 哈尔滨工业大学 | 一种桥塔尾流区吊索风致振动数值预测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110110408B (zh) | 2022-06-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110110408A (zh) | 刚性圆柱体横流与顺流方向涡激振动耦合响应预测方法 | |
| Assi et al. | The role of wake stiffness on the wake-induced vibration of the downstream cylinder of a tandem pair | |
| Govardhan et al. | Resonance forever: existence of a critical mass and an infinite regime of resonance in vortex-induced vibration | |
| Wu et al. | Experimental and numerical investigation of hydroelastic response of a flexible hydrofoil in cavitating flow | |
| Mittal et al. | Finite element study of vortex‐induced cross‐flow and in‐line oscillations of a circular cylinder at low Reynolds numbers | |
| Raphaël et al. | Capillary gravity waves caused by a moving disturbance: wave resistance | |
| CN111985138B (zh) | 一种柔性结构横流与顺流方向涡激振动耦合响应预测方法 | |
| Toebes et al. | Hydroelastic vibrations of flat plates related to trailing edge geometry | |
| Wu et al. | Numerical investigation of vortex-induced vibration of a circular cylinder with a hinged flat plate | |
| Han et al. | Energy harvesting from flow-induced vibration of a low-mass square cylinder with different incidence angles | |
| CN110046451A (zh) | 变张力细长柔性圆柱体涡激振动响应预测方法 | |
| Luongo et al. | A continuous approach to the aeroelastic stability of suspended cables in 1: 2 internal resonance | |
| Fang et al. | Theoretical and experimental investigation on the advantages of auxetic nonlinear vortex-induced vibration energy harvesting | |
| Ding et al. | Large eddy simulation of unsteady aerodynamic behavior of long-span vaulted roofs | |
| Du et al. | Enhancement of the Piezoelectric Cantilever Beam Performance via Vortex‐Induced Vibration to Harvest Ocean Wave Energy | |
| Liu et al. | Experimental investigation on high-mode vortex-induced vibration of a flexible stay cable in smooth flow | |
| Verma et al. | Dynamics of vortex-induced-vibrations of a slit-offset circular cylinder for energy harvesting at low Reynolds number | |
| Hafizh et al. | A hybrid piezoelectric–electromagnetic nonlinear vibration energy harvester excited by fluid flow | |
| CN102353516A (zh) | 一种深水立管的时域脉动拖曳力确定方法 | |
| González Cornejo et al. | Modelling low Reynolds number vortex-induced vibration problems with a fixed mesh fluid-solid interaction formulation | |
| CN108090263A (zh) | 输电线脱冰振动缩尺试验气动阻尼的验证方法 | |
| Metzner | Behavior of suspended matter in rapidly accelerating viscoelastic fluids: The Uebler effect | |
| CN102507082B (zh) | 一种深水立管的时域涡激升力确定方法 | |
| Paré-Lambert et al. | A parametric study of energy extraction from vortex-induced vibrations | |
| Song et al. | A piezoelectric energy harvester with vortex induced vibration |
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 | ||
| 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: 20220607 |