CN111400976B - 波浪作用下圆柱阵列中j柱对k柱激振力影响的计算方法 - Google Patents
波浪作用下圆柱阵列中j柱对k柱激振力影响的计算方法 Download PDFInfo
- Publication number
- CN111400976B CN111400976B CN202010184471.8A CN202010184471A CN111400976B CN 111400976 B CN111400976 B CN 111400976B CN 202010184471 A CN202010184471 A CN 202010184471A CN 111400976 B CN111400976 B CN 111400976B
- Authority
- CN
- China
- Prior art keywords
- column
- columns
- exciting force
- representing
- cylindrical array
- 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.)
- Active
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
技术领域
本发明涉及海洋工程领域,特别涉及任意排布的穿透水面圆柱阵列在波浪作用下,基于线性势流理论计算任意j柱对k柱激振力影响的方法。
背景技术
圆柱阵列是一种常用的海上结构物,广泛应用于海洋石油平台、跨海大桥、波浪发电单元等。布置在海洋中的坐底圆柱阵列所承受的激振力是决定圆柱阵列设计方案,保证结构安全性的一个重要因素,为此需要掌握各柱激振力的具体组成成分。目前对该问题的研究主要采用势流理论和线性波理论假设下的特征函数展开法。该方法由Spring&Monkeyer(1974)首先提出,后来Linton&Evans(1990)又对该方法进行了简化。圆柱阵列中某一柱受到的激振力可以划分为两个部分:一部分是入射波对该柱的作用,另一部分是其他圆柱造成的波浪绕射对该柱的作用。运用现有的计算方法能够计算出阵列中各柱的激振力,但无法明确其组成成分,不利于对圆柱阵列进行优化设计。
发明内容
本发明为解决公知技术中存在的技术问题而提供一种,波浪作用下圆柱阵列中j柱对k柱激振力影响的计算方法,该方法能够明确各柱所受激振力的具体组成成分。
本发明为解决公知技术中存在的技术问题所采取的技术方案是:一种波浪作用下圆柱阵列中j柱对k柱激振力影响的计算方法,采用以下公式:
式中:Fx(j→k)表示第j柱的存在对第k柱产生的x方向激振力;ρ表示流体密度,d表示水深,A表示入射波波幅,g表示重力加速度,κ是色散关系κtanh(κd)=ω2/g的正实根,Hn+1、Hn-1分别表示阶数为n+1和n-1的第一类整数阶汉开尔函数,H-1′、H1′分别表示阶数为-1和1的第一类整数阶汉开尔函数的一阶导数;Rjk表示第j柱圆心和第k柱圆心的距离,αjk表示第j柱圆心和第k柱圆心连线与x轴夹角;由等式定义,a表示圆柱半径;表示未知系数,可由下式求得:
本发明具有的优点和积极效果是:明确了各柱所受激振力的具体组成成分,有利于对圆柱阵列排布的优化设计。例如,为了降低阵列中某一柱的激振力,可以首先计算出其激振力具体来源,通过改变阵列排布方式或者增减圆柱数量来达到目的。
附图说明
图1为本发明的圆柱阵列一般形式的坐标定义示意图。
具体实施方式
为能进一步了解本发明的发明内容、特点及功效,兹例举以下实施例,并配合附图详细说明如下:
一种波浪作用下圆柱阵列中j柱对k柱激振力影响的计算方法,采用以下步骤:
一)建立坐标系
请参见图1,圆柱阵列一般形式的坐标定义示意图。假设圆柱阵列共有N个圆柱,每个圆柱半径均为a。以迎浪侧首柱横截面圆心o为原点,建立总体平面直角坐标系xoy,该平面位于静止的自由表面上,原点o位于该面上x轴和y轴的交点处,z轴垂直向上。单个柱体的圆心oj在总体直角坐标系下的坐标为(xj,yj)。建立N个原点在圆心oj的局部极坐标系ojrjθj。这样,第k个柱体中心Ok在第j柱局部极坐标系中的坐标可记为(Rjk,αjk)j,k=1…N,其中Rjk表示第j柱圆心和第k柱圆心的距离,αjk表示第j柱圆心和第k柱圆心连线与x轴夹角。图中β表示波浪的入射角度。
二)根据入射波频率ω,获得特征值κ
在线性波浪理论的假设下,流场内速度势可表示为如下形式:
Φ(x,y,z,t)=Re[φ(x,y)f(z)e-iωt]
式中,ω表示波浪运动的角频率,Re表示取复数的实部。f(z)有如下表达式:
上式中,d表示水深,A表示入射波波幅,g表示重力加速度,κ是色散关系中的正实根,色散关系表示如下:
κtanh(κd)=ω2/g
根据入射波频率ω求解上述超越方程,得到特征值κ。
三)分解速度势
在圆柱阵列第k柱附近的速度势可进行如下分解:
φ(rk,θk)=φI+φDk+φIk2
=φI+φDk1+φDk2+φIk2
上式中:φ(rk,θk)表示第k柱附近的总速度势;φDk表示第k柱附近的总绕射速度势;φDk1表示第k柱上由入射波作用引起的绕射速度势;φDk2表示第k柱上由其他柱子的绕射速度势作为二次入射波作用到第k柱上引起的二次绕射势;φIk2表示其他所有圆柱在第k柱附近产生的绕射势。φI+φDk1表示入射波对孤立k柱的影响;φDk2+φIk2表示其他柱子的存在对第k柱的影响。
φI表示环境入射势,在局部极坐标系ojrjθj(j=1,...,N)下,点(rj,θj)的入射速度势可表示为:
四)φDk1的求解
φDk1具有以下形式:
上式中,Hm表示第一类整数阶汉开尔函数,m表示阶数,cm0表示未知系数。
最终得到φDk1的表达式如下:
五)φIk2及φDk2的求解
在第k柱极坐标系下,φIk2有如下形式:
(k=1,...,N,-∞<m<∞)
φDk2具有以下形式:
故得到φDk2的表达式如下:
综上可得到其他柱子的存在对第k柱的影响φDk2+φIk2,其表达式如下:
六)第j柱的存在对第k柱激振力的影响
第j柱的存在对第k柱产生的速度势可表示为如下形式:
上式中,Fx(j→k)表示第j柱的存在对第k柱产生的x方向激振力。
七)其余所有圆柱对第k柱的影响
除了第k柱,其余所有圆柱的存在对第k柱x方向激振力的影响可表示为下式:
Claims (1)
1.一种波浪作用下圆柱阵列中j柱对k柱激振力影响的计算方法,其特征在于,采用以下公式:
式中:Fx(j→k)表示第j柱的存在对第k柱产生的x方向激振力;a表示圆柱半径;ρ表示流体密度,d表示水深,A表示入射波波幅,g表示重力加速度,κ是色散关系κtanh(κd)=ω2/g的正实根,Hn+1、Hn-1分别表示阶数为n+1和n-1的第一类整数阶汉开尔函数,H-1′、H1′分别表示阶数为-1和1的第一类整数阶汉开尔函数的一阶导数;Rjk表示第j柱圆心和第k柱圆心的距离,αjk表示第j柱圆心和第k柱圆心连线与x轴夹角;由等式定义,Jn′表示第一类整数阶贝塞尔函数的一阶导数;表示未知系数,可由下式求得:
k=1,...,N,-∞<m<∞
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010184471.8A CN111400976B (zh) | 2020-03-16 | 2020-03-16 | 波浪作用下圆柱阵列中j柱对k柱激振力影响的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010184471.8A CN111400976B (zh) | 2020-03-16 | 2020-03-16 | 波浪作用下圆柱阵列中j柱对k柱激振力影响的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111400976A CN111400976A (zh) | 2020-07-10 |
CN111400976B true CN111400976B (zh) | 2022-05-06 |
Family
ID=71434264
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010184471.8A Active CN111400976B (zh) | 2020-03-16 | 2020-03-16 | 波浪作用下圆柱阵列中j柱对k柱激振力影响的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111400976B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2039296A1 (en) * | 1990-04-13 | 1991-10-14 | Kenneth J. Gibson | Methods and apparatus for optimizing prefetch caching |
CN109726423A (zh) * | 2018-07-17 | 2019-05-07 | 中国科学院力学研究所 | 圆柱阵列波浪力幅值波动特性、步长和包络线的获取方法 |
CN109726418A (zh) * | 2018-07-17 | 2019-05-07 | 中国科学院力学研究所 | 圆柱阵列的波浪力幅值在非捕获区波动间距的描述模型 |
-
2020
- 2020-03-16 CN CN202010184471.8A patent/CN111400976B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2039296A1 (en) * | 1990-04-13 | 1991-10-14 | Kenneth J. Gibson | Methods and apparatus for optimizing prefetch caching |
CN109726423A (zh) * | 2018-07-17 | 2019-05-07 | 中国科学院力学研究所 | 圆柱阵列波浪力幅值波动特性、步长和包络线的获取方法 |
CN109726418A (zh) * | 2018-07-17 | 2019-05-07 | 中国科学院力学研究所 | 圆柱阵列的波浪力幅值在非捕获区波动间距的描述模型 |
Non-Patent Citations (2)
Title |
---|
"The interaction of waves with arrays of vertical circular cylinders";C.M.LINTON.etc;《Fluid Mech》;19901231;第215卷;全文 * |
"波浪作用下的直立圆柱阵列水动力波动特性研究";张华庆等;《水道港口》;20191031;第40卷(第5期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111400976A (zh) | 2020-07-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kristiansen et al. | Higher harmonic wave loads on a vertical cylinder in finite water depth | |
CN110794459B (zh) | 一种海底近断层地震动的拟合方法 | |
Achenbach | Reciprocity in elastodynamics | |
Lighthill | Fundamentals concerning wave loading on offshore structures | |
Jiao et al. | Time-domain numerical and segmented model experimental study on ship hydroelastic responses and whipping loads in harsh irregular seaways | |
Hayatdavoodi et al. | Diffraction of cnoidal waves by vertical cylinders in shallow water | |
Jian et al. | Short crested wave–current forces around a large vertical circular cylinder | |
Wang et al. | An analytical method for the solution of earthquake-induced hydrodynamic forces and wave forces on multiple vertical cylinders with arbitrary smooth sections | |
CN111400976B (zh) | 波浪作用下圆柱阵列中j柱对k柱激振力影响的计算方法 | |
Isaacson et al. | Transverse forces on vertical cylinders in waves | |
CN104483702B (zh) | 一种适用于非均匀运动水体的地震正演模拟方法 | |
Olaya et al. | Hydrodynamic coefficients and wave loads for a WEC device in heaving mode | |
Tian et al. | Theoretical study on hydroelastic responses of very large floating structures near islands and reefs | |
LIN et al. | Observation and Theoretical Analysis for the Sand‐Waves Migration in the North Gulf of the South China Sea | |
Gryazin et al. | A new method to determine directional spectrum of sea waves and its application to wave buoys | |
Teng et al. | Motion response of hinged multiple floating bodies on local seabed | |
Milgram | Waves and wave forces | |
Liu et al. | Using a time-domain higher-order boundary element method to simulate wave and current diffraction from a 3-D body | |
Sen | Time-domain simulation of motions of large structures in nonlinear waves | |
Qiao et al. | Numerical simulation for truncated model tests of deepwater semi-submersible platform with viscous damper compensated system in mooring lines | |
Harun et al. | Analytical Solution For a Long Waves Propagation in Two-Layer Fluid Over a Submerged Hump | |
Ali et al. | On hydrodynamic interaction between several freely floating vertical cylinders in waves | |
Ma et al. | An Experimental Investigation on Concomitant Vortex-Induced Vibration and Cross Flow Top-Motion Excitation with Flexible Risers | |
Nτίντας | Wind wave analysis of coastal waters using OCEANLYZ toolbox | |
Kim et al. | Hydrodynamic characteristics of prismatic barges |
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 | ||
CB02 | Change of applicant information |
Address after: No.2618, No.2 Road, Tanggu, Binhai New Area, Tianjin Applicant after: Tianjin water transport engineering survey and Design Institute Co.,Ltd. Address before: No.2618, No.2 Road, Tanggu, Binhai New Area, Tianjin Applicant before: TIANJIN SURVEY AND DESIGN INSTITUTE OF WATER TRANSPORT ENGINEERING |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |