CN110598337A - Fluid-solid coupling time domain analysis method for vortex-induced vibration of cylinder - Google Patents
Fluid-solid coupling time domain analysis method for vortex-induced vibration of cylinder Download PDFInfo
- Publication number
- CN110598337A CN110598337A CN201910876543.2A CN201910876543A CN110598337A CN 110598337 A CN110598337 A CN 110598337A CN 201910876543 A CN201910876543 A CN 201910876543A CN 110598337 A CN110598337 A CN 110598337A
- Authority
- CN
- China
- Prior art keywords
- cylinder
- vortex
- fluid
- frequency
- formula
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Complex Calculations (AREA)
Abstract
The invention relates to a research method of an ocean deepwater riser, in particular to a fluid-solid coupling time domain analysis method of cylinder vortex-induced vibration. The key points of the invention comprise: calculating the fluid force and the vortex shedding frequency of the non-locking area by using the relative speed of the fluid and the cylinder; under the condition of fluid-solid coupling, the pulsating drag force and the vortex-induced lift force change along with the change of the relative speed of the fluid and the cylinder; the frequency of the fluid force is kept unchanged in the locking region, and the important characteristic of vortex-induced vibration, namely frequency locking, is fully embodied; the locking frequency depends not only on the wet modal frequency of the cylinder, but also on the mass ratio of the cylinder to the fluid. The invention has the following advantages: the influence of the cylinder motion on the magnitude and direction of the pulsating drag force and the vortex-induced lift force is correctly described; the property of the fluid force frequency of the locked zone and the unlocked zone, namely the relationship between the pulsating drag force frequency and the vortex induced lift force frequency is correctly described; the locking phenomenon of vortex shedding frequency is correctly described.
Description
Technical Field
The invention relates to a research method of an ocean deepwater riser, in particular to a fluid-solid coupling time domain analysis method of cylinder vortex-induced vibration.
Background
The time domain analysis method of the vortex-induced vibration of the cylinder mainly comprises a wake flow vibrator model and a force decomposition model, wherein the force decomposition model is used for decomposing fluid force borne by the cylinder which generates the vortex-induced vibration into two components in the flow velocity direction and the vertical flow velocity direction of a flow field and respectively calculating downstream vortex-induced vibration and transverse vortex-induced vibration of the cylinder.
Vortex-induced vibrations are generated when a cylinder is subjected to a fluid in a non-stationary flow field whose flow velocity is not parallel to its axis, the flow causing the vortex-induced vibrations being a component of the flow velocity perpendicular to the axis of the cylinder. Thus, the vortex-induced vibrations of the cylinder are all calculated using a flow velocity perpendicular to the axis of the cylinder, as shown in fig. 2.
For the purpose of describing the vortex-induced vibration of a cylinder, it is generally decomposed into two components parallel to the flow velocity and perpendicular to the flow velocity, which are called downstream vibration and transverse vibration, respectively, and the corresponding two fluid force components are called pulsating drag force and vortex-induced lift force, respectively. The pulsating drag force is defined as the component of the cylinder surface pressure integral in the direction of flow velocity; the vortex-induced lift is defined as the integral of the cylinder surface pressure in the direction perpendicular to the flow velocity. The above definition holds if the stagnation point of the fluid (see fig. 3 for the definition of the stagnation point) is located on the negative x-axis (see fig. 2) and the separation point is symmetrical to the x-axis (see fig. 3).
The stagnation point position of fig. 3 is the case when the cylinder is streaming (the cylinder is stationary in the flow field), when the cylinder is moving in the flow field non-parallel to the direction of the flow velocity, the direction of the velocity of the fluid relative to the cylinder is no longer the direction of the velocity of the flow field, so that the stagnation point and the separation point will also deviate from the position when the cylinder is streaming, as shown in fig. 4.
Therefore, when the cylinder generates vortex-induced vibration, the pulsating drag force applied to the cylinder is no longer along the x-axis direction, and the vortex-induced lift force is no longer along the y-axis direction. However, the current method defaults that even if the cylinder generates vortex-induced vibration, the pulsating drag force of the cylinder still follows the x-axis direction, and the vortex-induced lift force still follows the y-axis direction. Then, according to the theory of forced vibration, the downstream vortex-induced vibration frequency of the cylinder should be the same as the pulsating drag force frequency, and the transverse vortex-induced vibration frequency should be the same as the vortex-induced lift force frequency. This means that the components of the cylinder vortex-induced vibration in the x-axis and the y-axis respectively have only one frequency peak, and the result of the test is that each component has two frequency peaks, the main peak is the same as the corresponding fluid force frequency, and the secondary peak is the same as the fluid force frequency in the other direction, that is, the main peak of the downstream response frequency is the same as the pulsation drag force frequency, and the secondary peak is the same as the vortex-induced lift force frequency; the main peak of the transverse response frequency is the same as the vortex induced lift frequency and the secondary peak is the same as the pulsating drag frequency, as shown in figure 5.
In summary, the prior art has the following disadvantages:
1. the effect of the lateral motion of the cylinder on the fluid-solid coupling is not considered.
2. The effect of stagnation point and separation point shifts on fluid forces due to cylinder motion is not considered.
3. There is no reaction frequency locking phenomenon-in the locking area, the load frequency remains unchanged. Only the locked and unlocked regions are considered, the difference in the pulsating drag frequency and the vortex induced lift frequency.
4. The effect of the quality ratio on the locking frequency is not taken into account.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention provides a fluid-solid coupling time domain analysis method for vortex-induced vibration of a cylinder.
The technical scheme of the invention is as follows:
the invention provides a fluid-solid coupling time domain analysis method for vortex-induced vibration of a cylinder, which comprises the following steps:
1) establishing vortex-induced vibration equation of cylinder
Wherein:
locking zone (4)
Non-locking zone (5)
In the formula: f. ofD-a pulsating drag force;
fL-vortex induced lift;
ρ -fluid density;
d-cylinder diameter;
-vortex shedding frequency;
t is time;
St-Strouhal number;
m is cylinder mass;
ma-a fluid additional mass;
c-structural damping of the cylinder;
ca-additional damping of the fluid;
k-cylinder stiffness;
x,-displacement, velocity and acceleration of the cylinder in the forward direction;
y,-cylinder lateral displacement, velocity and acceleration;
C′D-a pulsating drag force coefficient;
C′L-vortex induced lift coefficient;
u-velocity vector of the flow field;
v-the velocity vector of the cylinder, the coordinate components of which are eachAndnamely:i and j are unit vectors of the x axis and the y axis respectively;
mu is the frequency coefficient of the pulsating drag force, the non-locking area mu is 1, and the locking area mu is 2;
beta is the included angle between the relative velocity vector U-v and the flow velocity vector U;
fn-the wet modal frequency of the cylinder;
m-mass ratio.
2) Establishing an incremental finite element equation of the formula (1), wherein the formula (1) is an implicit equation, and therefore, an iterative method is adopted for calculation:
in the formula: [ M ] -cylinder mass matrix;
[C] -a cylinder damping matrix;
[K] -a cylinder stiffness matrix;
-node increment additionThe speed of the motor is controlled by the speed of the motor,
-the node incremental speed is calculated,
{ Δ a } -node incremental displacement, { Δ a } - [ Δ x ═1,Δy1,…,Δxn,Δyn]T;
{ Δ F } -node incremental load, { Δ F } - [ Δ F } - ]x1,ΔFy1,…,ΔFxn,ΔFyn]T,
ΔFx=ΔfDcosβ-ΔfLsinβ,ΔFy=ΔfDsinβ+ΔfLcosβ;
n is the total number of nodes.
3) Let t be 0 and k be 0, and convert the initial values toAndcalculating by sequentially substituting the formulas (6), (4)/(5), (2) and (3)Then substituted into formula (7) and calculated by time course analysis methodAnd
4) will be provided withAndcalculating by sequentially substituting the formulas (6), (4)/(5), (2) and (3)Then substituted into formula (7) and calculated by time course analysis methodAndin this case, t is Δ t, and k is 1. And repeating the iteration until the convergence condition is:
and the precision condition that epsilon is set by a calculator is met. At this time, t ═ Δ t, k ═ s, and s is the number of iterations in the time step.
5) Let t be t +. DELTA.t and k be 0, and start the calculation of the next time step. I.e. the result of the previous time stepAnd(s is the number of iterations of the last time step) are sequentially substituted into the formulas (6), (4)/(5), (2) and (3) for calculationThen substituted into formula (7) and calculated by time course analysis methodAnd
6) will be provided withAndcalculating by sequentially substituting the formulas (6), (4)/(5), (2) and (3)Then substituted into formula (7) and calculated by time course analysis methodAndwhen t is t +. DELTA.t, k is 1; and repeating the iteration until the convergence condition is:
is satisfied; at this time, t ═ t +. DELTA.t, k ═ s, s is the number of iterations of the time step, and the two steps 5 and 6 are repeated until the time t reaches the calculation time length preset by the calculator.
The invention achieves the following beneficial effects:
(1) the invention correctly describes the influence of cylinder motion on the magnitude and direction of the pulsating drag force and the vortex induced lift force.
(2) The present invention correctly describes the property of the fluid force frequency of the lock zone and the unlock zone, namely the relationship between the pulsating drag force frequency and the vortex induced lift force frequency.
(3) The present invention correctly describes the locking phenomenon of vortex shedding frequency.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
FIG. 2 is a schematic of a cylinder and flow rate.
Fig. 3 is a schematic view of stagnation and separation points formed by the streaming.
FIG. 4 is a schematic diagram of the stagnation point and the separation point of the cylinder moving along the y-axis.
Fig. 5 is a graph of measured vortex excitation motion response spectra.
Where U is the flow rate and o is the cylinder axis.
Detailed Description
To facilitate an understanding of the present invention by those skilled in the art, specific embodiments thereof are described below with reference to the accompanying drawings.
As shown in fig. 1, the present invention provides a fluid-solid coupling time domain analysis method for vortex-induced vibration of a cylinder, including:
1) s1 establishing vortex-induced vibration equation of cylinder
Wherein:
locking zone (4)
Non-locking zone (5)
In the formula: f. ofD-a pulsating drag force;
fL-vortex induced lift;
ρ -fluid density;
d-cylinder diameter;
-vortex shedding frequency;
t is time;
St-Strouhal number;
m is cylinder mass;
ma-a fluid additional mass;
c-structural damping of the cylinder;
ca-additional damping of the fluid;
k-cylinder stiffness;
x,-displacement, velocity and acceleration of the cylinder in the forward direction;
y,-cylinder lateral displacement, velocity and acceleration;
C′D-a pulsating drag force coefficient;
C′L-vortex induced lift coefficient;
u-velocity vector of the flow field;
v-the velocity vector of the cylinder, the coordinate components of which are eachAndnamely:i and j are unit vectors of the x axis and the y axis respectively;
mu is the frequency coefficient of the pulsating drag force, the non-locking area mu is 1, and the locking area mu is 2;
beta is the included angle between the relative velocity vector U-v and the flow velocity vector U;
fn-the wet modal frequency of the cylinder;
m-mass ratio.
2) S2, establishing an incremental finite element equation of formula (1), and since formula (1) is an implicit equation, calculating by using an iterative method:
in the formula: [ M ] -cylinder mass matrix;
[C] -a cylinder damping matrix;
[K] -a cylinder stiffness matrix;
-the incremental acceleration of the node(s),
-the node incremental speed is calculated,
{ Δ a } -node incremental displacement, { Δ a } - [ Δ x ═1,Δy1,…,Δxn,Δyn]T;
{ Δ F } -node incremental load, { Δ F } - [ Δ F } - ]x1,ΔFy1,…,ΔFxn,ΔFyn]T,
ΔFx=ΔfDcosβ-ΔfLsinβ,ΔFy=ΔfDsinβ+ΔfLcosβ;
n is the total number of nodes.
3) S3 compares t with 0 and k with 0, and calculates the initial valueAndcalculating by sequentially substituting the formulas (6), (4)/(5), (2) and (3)Then substituted into formula (7) and calculated by time course analysis methodAnd
4) s4 is toAndcalculating by sequentially substituting the formulas (6), (4)/(5), (2) and (3)Then substituted into formula (7) and calculated by time course analysis methodAndin this case, t is Δ t, and k is 1. S5 iterates so many times until a convergence condition:
and the precision condition that epsilon is set by a calculator is met. At this time, t ═ Δ t, k ═ s, and s is the number of iterations in the time step.
5) S6 starting the next time step by making t equal to t +. DELTA.t and k equal to 0And (4) calculating. I.e. the result of the previous time stepAnd(s is the number of iterations of the last time step) are sequentially substituted into the formulas (6), (4)/(5), (2) and (3) for calculationThen substituted into formula (7) and calculated by time course analysis methodAnd
6) s7 is toAndcalculating by sequentially substituting the formulas (6), (4)/(5), (2) and (3)Then substituted into formula (7) and calculated by time course analysis methodAndin this case, t + Δ t and k are 1. S8 iterates so many times until a convergence condition:
is satisfied; at this time, t + Δ t, k is s, and s is the number of iterations in the time step. And S9, repeating the steps 5 and 6 until the time t reaches the calculation time preset by the calculator.
The above-described embodiments of the present invention do not limit the scope of the present invention. Any modification, equivalent replacement, and improvement made within the spirit and principle of the present invention should be included in the scope of the claims of the present invention.
Claims (1)
1. A fluid-solid coupling time domain analysis method for cylinder vortex-induced vibration is characterized by comprising the following steps:
1) establishing vortex-induced vibration equation of cylinder
Wherein:
in the formula: f. ofD-a pulsating drag force;
fL-vortex induced lift;
ρ -fluid density;
d-cylinder diameter;
-vortex shedding frequency;
t is time;
St-Strouhal number;
m is cylinder mass;
ma-a fluid additional mass;
c-structural damping of the cylinder;
ca-additional damping of the fluid;
k-cylinder stiffness;
x,-displacement, velocity and acceleration of the cylinder in the forward direction;
y,-cylinder lateral displacement, velocity and acceleration;
C′D-a pulsating drag force coefficient;
C′L-vortex induced lift coefficient;
u-velocity vector of the flow field;
v-the velocity vector of the cylinder, the coordinate components of which are eachAndnamely:
i and j are unit vectors of the x axis and the y axis respectively;
mu is the frequency coefficient of the pulsating drag force, the non-locking area mu is 1, and the locking area mu is 2;
beta is the included angle between the relative velocity vector U-v and the flow velocity vector U;
fn-the wet modal frequency of the cylinder;
m-mass ratio;
2) establishing an incremental finite element equation of the formula (1), wherein the formula (1) is an implicit equation, and therefore, an iterative method is adopted for calculation:
in the formula: [ M ] -cylinder mass matrix;
[C] -a cylinder damping matrix;
[K] -a cylinder stiffness matrix;
-the incremental acceleration of the node(s),
-the node incremental speed is calculated,
{ Δ a } -node incremental displacement, { Δ a } - [ Δ x ═1,Δy1,…,Δxn,Δyn]T;
{ Δ F } -node incremental load, { Δ F } - [ Δ F } - ]x1,ΔFy1,…,ΔFxn,ΔFyn]T,
ΔFx=ΔfDcosβ-ΔfLsinβ,ΔFy=ΔfDsinβ+ΔfLcosβ;
n is the total number of nodes;
3) let t be 0 and k be 0, and convert the initial values toAndcalculating by sequentially substituting the formulas (6), (4)/(5), (2) and (3)Then substituted into formula (7) and calculated by time course analysis methodAnd
4) will be provided withAndcalculating by sequentially substituting the formulas (6), (4)/(5), (2) and (3)Then substituted into formula (7) and calculated by time course analysis methodAndwhen t ═ Δ t, k ═ 1; and repeating the iteration until the convergence condition is:
the accuracy condition that epsilon is set by a calculator is met; at the moment, t is delta t, k is s, and s is the iteration number of the time step;
5) let t be t +. DELTA.t and k be 0, start the calculation of the next time step; i.e. the result of the previous time stepAnd(s is the number of iterations of the last time step) are sequentially substituted into the formulas (6), (4)/(5), (2) and (3) for calculationThen substituted into formula (7) and calculated by time course analysis methodAnd
6) will be provided withAndcalculating by sequentially substituting the formulas (6), (4)/(5), (2) and (3)Then substituted into formula (7) and calculated by time course analysis methodAndwhen t is t +. DELTA.t, k is 1; repeating the above stepsIterating until a convergence condition is met:
is satisfied; at this time, t ═ t +. DELTA.t, k ═ s, s is the number of iterations of the time step, and the two steps 5 and 6 are repeated until the time t reaches the calculation time length preset by the calculator.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910876543.2A CN110598337B (en) | 2019-09-17 | 2019-09-17 | Fluid-solid coupling time domain analysis method for vortex-induced vibration of cylinder |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910876543.2A CN110598337B (en) | 2019-09-17 | 2019-09-17 | Fluid-solid coupling time domain analysis method for vortex-induced vibration of cylinder |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110598337A true CN110598337A (en) | 2019-12-20 |
CN110598337B CN110598337B (en) | 2021-01-26 |
Family
ID=68860138
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910876543.2A Expired - Fee Related CN110598337B (en) | 2019-09-17 | 2019-09-17 | Fluid-solid coupling time domain analysis method for vortex-induced vibration of cylinder |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110598337B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111985138A (en) * | 2020-08-21 | 2020-11-24 | 哈尔滨工业大学(威海) | Flexible structure transverse flow and downstream direction vortex-induced vibration coupling response prediction method |
CN113758666A (en) * | 2021-08-25 | 2021-12-07 | 浙江大学 | Experiment table and method for testing flow-induced vibration of valve rod of valve |
CN115859748A (en) * | 2023-02-15 | 2023-03-28 | 山东科技大学 | Flexible cable vortex-induced vibration analysis method for dragging type thermohaline depth measuring instrument |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102353516A (en) * | 2011-09-29 | 2012-02-15 | 中国海洋大学 | Method for determining time-domain pulsating drag force of deepwater riser |
CN102507082A (en) * | 2011-09-30 | 2012-06-20 | 中国海洋大学 | Time domain vortex-induced lift force determination method of deep-water riser |
CN103353382A (en) * | 2013-07-10 | 2013-10-16 | 天津大学 | Analyzing method for parametric-excitation and vortex-induced vibratory fatigue of deep-ocean top tension riser |
KR20140046681A (en) * | 2012-10-10 | 2014-04-21 | 대우조선해양 주식회사 | Mooring system of tension leg platform |
CN105205280A (en) * | 2015-10-10 | 2015-12-30 | 天津大学 | Prediction method for single-span pipeline transverse flow direction vortex-induced vibration |
KR20160057153A (en) * | 2014-11-13 | 2016-05-23 | 대우조선해양 주식회사 | Rotating apparatus of reduced vortex induced vibration and offshore platform including the smae |
CN106679791A (en) * | 2016-12-15 | 2017-05-17 | 天津大学 | Simulation device for vortex-induced vibration of submarine pipeline and experimental method |
CN108205516A (en) * | 2017-12-26 | 2018-06-26 | 中国海洋大学 | A kind of ocean structure fluid structurecoupling dynamic response algorithm transmitted based on condition |
CN109740233A (en) * | 2018-12-27 | 2019-05-10 | 大连船舶重工集团有限公司 | Hard sail vortex-induced vibration based on fluid structurecoupling simulation responds preparation method |
CN109799049A (en) * | 2019-03-06 | 2019-05-24 | 北京理工大学 | A kind of elastic cylinder vortex-induced vibration rule and coupling mechanism measuring method |
CN109827734A (en) * | 2019-03-29 | 2019-05-31 | 西安建筑科技大学 | A kind of method that outflow acts on lower deep sea vertical pipe vortex-induced vibration in assessment |
CN109885864A (en) * | 2019-01-07 | 2019-06-14 | 长沙理工大学 | A kind of three-dimensional steel bridge tower vortex-induced vibration calculation method |
CN110110408A (en) * | 2019-04-25 | 2019-08-09 | 西南石油大学 | Rigid cylinder crossing current and downbeam vortex-induced vibration coupling response prediction technique |
-
2019
- 2019-09-17 CN CN201910876543.2A patent/CN110598337B/en not_active Expired - Fee Related
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102353516A (en) * | 2011-09-29 | 2012-02-15 | 中国海洋大学 | Method for determining time-domain pulsating drag force of deepwater riser |
CN102507082A (en) * | 2011-09-30 | 2012-06-20 | 中国海洋大学 | Time domain vortex-induced lift force determination method of deep-water riser |
KR20140046681A (en) * | 2012-10-10 | 2014-04-21 | 대우조선해양 주식회사 | Mooring system of tension leg platform |
CN103353382A (en) * | 2013-07-10 | 2013-10-16 | 天津大学 | Analyzing method for parametric-excitation and vortex-induced vibratory fatigue of deep-ocean top tension riser |
KR20160057153A (en) * | 2014-11-13 | 2016-05-23 | 대우조선해양 주식회사 | Rotating apparatus of reduced vortex induced vibration and offshore platform including the smae |
CN105205280A (en) * | 2015-10-10 | 2015-12-30 | 天津大学 | Prediction method for single-span pipeline transverse flow direction vortex-induced vibration |
CN106679791A (en) * | 2016-12-15 | 2017-05-17 | 天津大学 | Simulation device for vortex-induced vibration of submarine pipeline and experimental method |
CN108205516A (en) * | 2017-12-26 | 2018-06-26 | 中国海洋大学 | A kind of ocean structure fluid structurecoupling dynamic response algorithm transmitted based on condition |
CN109740233A (en) * | 2018-12-27 | 2019-05-10 | 大连船舶重工集团有限公司 | Hard sail vortex-induced vibration based on fluid structurecoupling simulation responds preparation method |
CN109885864A (en) * | 2019-01-07 | 2019-06-14 | 长沙理工大学 | A kind of three-dimensional steel bridge tower vortex-induced vibration calculation method |
CN109799049A (en) * | 2019-03-06 | 2019-05-24 | 北京理工大学 | A kind of elastic cylinder vortex-induced vibration rule and coupling mechanism measuring method |
CN109827734A (en) * | 2019-03-29 | 2019-05-31 | 西安建筑科技大学 | A kind of method that outflow acts on lower deep sea vertical pipe vortex-induced vibration in assessment |
CN110110408A (en) * | 2019-04-25 | 2019-08-09 | 西南石油大学 | Rigid cylinder crossing current and downbeam vortex-induced vibration coupling response prediction technique |
Non-Patent Citations (9)
Title |
---|
C.LE CUNFF,ET AL.: "Vortex-Induced Vibrations of Risers: Theoretical, Numerical and Experimental Investigation", 《OIL & GAS SCIENCE AND TECHNOLOGY》 * |
JIANHUA PANG,ET AL.: "A numerical simulation model for the vortex induced vibration of flexible risers using dynamic stiffness matrices", 《OCEAN ENGINEERING》 * |
吴学敏 等: "考虑大变形的大柔性立管涡激振动模型", 《振动与冲击》 * |
娄敏 等: "尾流干涉下考虑流固耦合作用的海洋立管涡激振动的数值模拟", 《船海工程》 * |
徐诚侃: "考虑大变形的含缺陷海洋管道非线性涡激振动研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
朱文波: "海洋热塑性增强管涡激振动流固耦合分析", 《中国优秀硕士学位论文全文数据库 工程科技I辑》 * |
王鹏: "深水钻井隔水管系统动力响应及力学耦合行为", 《中国博士学位论文全文数据库 工程科技I辑》 * |
赵鹏良 等: "海洋立管涡激振动的流固耦合模拟计算", 《海洋技术》 * |
黄维平 等: "大柔性圆柱体两自由度涡激振动试验研究", 《力学学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111985138A (en) * | 2020-08-21 | 2020-11-24 | 哈尔滨工业大学(威海) | Flexible structure transverse flow and downstream direction vortex-induced vibration coupling response prediction method |
CN111985138B (en) * | 2020-08-21 | 2023-12-19 | 哈尔滨工业大学(威海) | Flexible structure cross flow and forward flow direction vortex-induced vibration coupling response prediction method |
CN113758666A (en) * | 2021-08-25 | 2021-12-07 | 浙江大学 | Experiment table and method for testing flow-induced vibration of valve rod of valve |
CN115859748A (en) * | 2023-02-15 | 2023-03-28 | 山东科技大学 | Flexible cable vortex-induced vibration analysis method for dragging type thermohaline depth measuring instrument |
Also Published As
Publication number | Publication date |
---|---|
CN110598337B (en) | 2021-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110598337B (en) | Fluid-solid coupling time domain analysis method for vortex-induced vibration of cylinder | |
Blackburn et al. | Two-and three-dimensional simulations of vortex-induced vibration of a circular cylinder | |
Rannacher et al. | A posteriori error estimation and mesh adaptation for finite element models in elasto-plasticity | |
Aronsen | An experimental investigation of in-line and combined in-line and cross-flow vortex induced vibrations | |
Beyer et al. | Faraday instability: linear analysis for viscous fluids | |
CN111985138B (en) | Flexible structure cross flow and forward flow direction vortex-induced vibration coupling response prediction method | |
Galimov et al. | Parallel adaptive simulation of a plunging liquid jet | |
Cardona et al. | Numerical integration of second order differential—algebraic systems in flexible mechanism dynamics | |
CN114819347B (en) | Method for predicting vibration response of flow transmission pipeline under multiphase internal flow excitation | |
Whittaker et al. | The energetics of flow through a rapidly oscillating tube. Part 2. Application to an elliptical tube | |
Tan et al. | Flow past rectangular cylinders: receptivity to transverse forcing | |
Chalot et al. | A consistent finite element approach to large eddy simulation | |
Parent et al. | Validation of the wilcox k-omega model for flows characteristic to hypersonic airbreathing propulsion | |
CN102353516A (en) | Method for determining time-domain pulsating drag force of deepwater riser | |
CN102507082B (en) | Time domain vortex-induced lift force determination method of deep-water riser | |
JP2007286801A (en) | Computing apparatus for finite element method for discretely analyzing high order differential equation | |
Duan et al. | Lyapunov-based stability analysis for fluid conveying system with parallel nonlinear energy sinks | |
CN112765725B (en) | Analytic Riemann resolving method for multi-dimensional Euler equation | |
Kubomura | Component mode synthesis for damped structures | |
Meinke et al. | Simulation of unsteady flows | |
Vahala et al. | Entropic, LES and boundary conditions in lattice Boltzmann simulations of turbulence | |
Olsen et al. | Application of OVERFLOW to hypersonic perfect gas flowfields | |
Blackburn et al. | Two-and three-dimensional simulations of vortex-induced vibration or a circular cylinder | |
Bui et al. | Reduced-order modeling of parametrized finite element solutions by the POD-ISAT technique. Application to aircraft air control systems | |
Rixen | A Lanczos procedure for efficient mode superposition in dynamic analysis |
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 |
Granted publication date: 20210126 Termination date: 20210917 |
|
CF01 | Termination of patent right due to non-payment of annual fee |