CN107153773A - 一种随机不规则波浪特征参数提取方法 - Google Patents

一种随机不规则波浪特征参数提取方法 Download PDF

Info

Publication number
CN107153773A
CN107153773A CN201710355031.2A CN201710355031A CN107153773A CN 107153773 A CN107153773 A CN 107153773A CN 201710355031 A CN201710355031 A CN 201710355031A CN 107153773 A CN107153773 A CN 107153773A
Authority
CN
China
Prior art keywords
wave
mrow
msub
irregular water
random irregular
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
Application number
CN201710355031.2A
Other languages
English (en)
Other versions
CN107153773B (zh
Inventor
金永平
万步炎
刘德顺
彭佑多
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hunan University of Science and Technology
Original Assignee
Hunan University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hunan University of Science and Technology filed Critical Hunan University of Science and Technology
Priority to CN201710355031.2A priority Critical patent/CN107153773B/zh
Publication of CN107153773A publication Critical patent/CN107153773A/zh
Application granted granted Critical
Publication of CN107153773B publication Critical patent/CN107153773B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00—Subject matter not provided for in other main groups of this subclass

Landscapes

  • Revetment (AREA)

Abstract

本发明公开了一种随机不规则波浪特征参数提取方法，首先输入随机不规则波浪所期望达到某等级海况所对应的有义波高及波能谱；然后，计算对应频率等分的波幅、频率和波数，并输入相位角；接着，利用随机不规则波浪运动的波面方程，生成波面坐标，通过波面坐标提取随机不规则波浪的波高，并按从大到小依次排列，再从波高大的部分开始提取全部波高的三分之一求平均值；最后，记录连续十个波高所对应波高三分之一部分的平均值，并通过判断所得平均值与给定海况的有义波高之间的稳定偏差是否均小于给定偏差来得到所需等级海况的波面坐标。本发明实现了随机不规则波浪特征参数提取，生成了精确的对应不同等级海况的随机不规则波浪的坐标参数。

Description

一种随机不规则波浪特征参数提取方法
技术领域
本发明涉及海洋工程领域，特别涉及一种随机不规则波浪特征参数提取方法。
背景技术
目前，在有关海上吊装作业的研究文献中，大多采用规则波浪作为外部激励进行分析，这种分析方式比较简单，但是这与实际海上吊装作业环境存在很大差异，实际海面上的波浪是随机不规则的，即随机不规则波浪的特征参数——有义波高随着时间的变化将呈现无规则的变化。人们通过长期的大量观测和分析，发现随机不规则波浪的特征参数存在一定的统计规律，并以随机不规则波浪的特征参数(有义波高)为参考来划分不同的海况等级，以便相关部门进行预报及分析计算。
在利用随机不规则波浪作为外部激励进行海上吊装作业装备的性能分析过程中，如何准确获取对应不同等级海况的随机不规则波浪的坐标参数将变得至关重要，也是现有分析方法必须突破的技术难点。
发明内容
为了解决上述技术问题，本发明提供一种实施方便、提取精度高的随机不规则波浪特征参数提取方法。
本发明解决上述问题的技术方案是：一种随机不规则波浪特征参数提取方法，包括以下步骤：
1)输入随机不规则波浪所期望达到某等级海况所对应的有义波高H1/3及波能谱S(ωj)；
2)对波能谱S(ωj)的频率进行离散化处理，将波能谱S(ωj)的频率范围ωmin～ωmax分成N等分，同时给定随机不规则波浪的浪向角χ；
3)计算对应频率等分的波幅、频率和波数；
4)输入在[0-2π]之间均匀分布的随机不规则波浪的相位角εj；
5)利用随机不规则波浪运动的波面方程，生成波面坐标；
6)提取随机不规则波浪的波高；并将从第1个随机不规则波浪的波高到第i个随机不规则波浪的波高按从大到小依次排列，再从波高大的部分开始提取全部波高的三分之一求平均值
7)给定偏差Δ；
8)记录从当前第i个随机不规则波浪的波高开始连续十个随机不规则波浪的波高所对应波高三分之一部分的平均值
9)判断步骤8)所得平均值与给定海况的有义波高之间的稳定偏差是否均小于给定偏差，即：
若是，则整个计算终止，所生成的随机不规则波浪的波面坐标即代表该等级海况的波面坐标；若不是，则令波面方程的计算时间t＝t+Δt，Δt为计算时间步长，并返回步骤3)继续产生随机不规则波浪的波面坐标，直至满足计算终止条件为止。
上述随机不规则波浪特征参数提取方法，所述步骤1)中，波能谱S(ωj)的公式为：
上述随机不规则波浪特征参数提取方法，所述步骤3)中，第j个频率等分对应的波幅ζj、频率ωj和波数kj的计算公式为：
和
式中：Tj为第j个频率等分对应的周期，g为重力加速度g＝9.8；当j＝1时，ω0＝ωmin。
上述随机不规则波浪特征参数提取方法，所述步骤5)中，波面方程为：
上述随机不规则波浪特征参数提取方法，所述步骤6)中，提取随机不规则波浪的波高的过程为：利用计算机分别提取步骤5)所生成的随机不规则波浪的波面坐标的波峰和波谷，将第i个波峰与相邻第i个波谷之间的垂直距离定义为对应第i个随机不规则波浪的波高并提取该波高
本发明的有益效果在于：本发明首先输入随机不规则波浪所期望达到某等级海况所对应的有义波高及波能谱；然后，计算随机不规则波浪的波幅、频率和波数，并输入相位角；接着，利用随机不规则波浪运动的波面方程，生成波面坐标，通过波面坐标提取随机不规则波浪的波高，并按从大到小依次排列，再从波高大的部分开始提取全部波高的三分之一求平均值；最后，记录连续十个波高所对应波高三分之一部分的平均值，并通过判断所得平均值与给定海况的有义波高之间的稳定偏差是否均小于给定偏差来得到所需等级海况的波面坐标，从而实现了随机不规则波浪特征参数提取，生成了精确的对应不同等级海况的随机不规则波浪的坐标参数，利用这些随机不规则波浪的坐标参数，开展海上吊装作业装备性能分析研究，将较以往利用规则波浪进行研究更加符合实际。
附图说明
图1为本发明的流程图。
图2为随机不规则波浪的有义波高逼近四级海况(有义波高2.1米)的过程图。
图3为四级海况对应的随机不规则波浪波形图。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明。
如图1所示，一种随机不规则波浪特征参数提取方法，包括以下步骤：
1)输入随机不规则波浪所期望达到某等级海况所对应的有义波高H1/3(四级海况对应的有义波高H1/3＝2.1米)及波能谱S(ωj)；波能谱S(ωj)的公式为：
2)对波能谱S(ωj)的频率进行离散化处理，将波能谱S(ωj)的频率范围ωmin～ωmax分成N等分，本例中N＝30，同时给定随机不规则波浪的浪向角χ，本例中χ＝30°。
3)计算对应频率等分的波幅、频率和波数；第j个随机不规则波浪的波幅ζj、频率ωj和波数kj的计算公式为：
和
式中：Tj为第j个频率等分对应的周期，g为重力加速度g＝9.8；当j＝1时，ω0＝ωmin。
4)输入在[0-2π]之间均匀分布的随机不规则波浪的相位角εj。
5)将步骤2)、步骤3)和步骤4)的计算结果代入到随机不规则波浪运动的波面方程：
利用随机不规则波浪运动的波面方程，生成波面坐标。
6)利用计算机分别提取步骤5)所生成的随机不规则波浪的波面坐标的波峰和波谷，将第i个波峰与相邻第i个波谷之间的垂直距离定义为对应第i个随机不规则波浪的波高并提取该波高然后，将从第1个随机不规则波浪的波高到第i个随机不规则波浪的波高按从大到小依次排列，再从波高大的部分开始提取全部波高的三分之一求平均值
7)给定偏差Δ，本例中Δ＝0.001m；
8)按照步骤6)的方式，计算并记录从当前第i个随机不规则波浪的波高开始连续十个随机不规则波浪的波高所对应波高三分之一部分的平均值
9)判断步骤8)所得平均值与给定海况的有义波高之间的稳定偏差是否均小于给定偏差，即：
若是，则整个计算终止，所生成的随机不规则波浪的波面坐标即代表该等级海况的波面坐标；若不是，则令波面方程的计算时间t＝t+Δt，Δt为计算时间步长，并返回步骤3)继续产生随机不规则波浪的波面坐标，直至满足计算终止条件为止。
随机不规则波浪的有义波高逼近四级海况(有义波高2.1米)过程如图2所示，四级海况对应的随机不规则波浪波形图如图3所示。

Claims (5)

1.一种随机不规则波浪特征参数提取方法，包括以下步骤：
1)输入随机不规则波浪所期望达到某等级海况所对应的有义波高H1/3及波能谱S(ωj)；
2)对波能谱S(ωj)的频率进行离散化处理，将波能谱S(ωj)的频率范围ωmin～ωmax分成N等分，同时给定随机不规则波浪的浪向角χ；
3)计算对应频率等分的波幅、频率和波数；
4)输入在[0-2π]之间均匀分布的随机不规则波浪的相位角εj；
5)利用随机不规则波浪运动的波面方程，生成波面坐标；
6)提取随机不规则波浪的波高；并将从第1个随机不规则波浪的波高到第i个随机不规则波浪的波高按从大到小依次排列，再从波高大的部分开始提取全部波高的三分之一求平均值
7)给定偏差Δ；
8)记录从当前第i个随机不规则波浪的波高开始连续十个随机不规则波浪的波高所对应波高三分之一部分的平均值
9)判断步骤8)所得平均值与给定海况的有义波高之间的稳定偏差是否均小于给定偏差，即：
<mrow> <mo>|</mo> <msub> <mover> <mover> <mi>h</mi> <mo>~</mo> </mover> <mo>&amp;OverBar;</mo> </mover> <msub> <mi>i</mi> <mrow> <mn>1</mn> <mo>/</mo> <mn>3</mn> </mrow> </msub> </msub> <mo>-</mo> <msub> <mi>H</mi> <mrow> <mn>1</mn> <mo>/</mo> <mn>3</mn> </mrow> </msub> <mo>|</mo> <mo>&lt;</mo> <mi>&amp;Delta;</mi> <mo>,</mo> <mo>|</mo> <msub> <mover> <mover> <mi>h</mi> <mo>~</mo> </mover> <mo>&amp;OverBar;</mo> </mover> <mrow> <mi>i</mi> <mo>+</mo> <msub> <mn>1</mn> <mrow> <mn>1</mn> <mo>/</mo> <mn>3</mn> </mrow> </msub> </mrow> </msub> <mo>-</mo> <msub> <mi>H</mi> <mrow> <mn>1</mn> <mo>/</mo> <mn>3</mn> </mrow> </msub> <mo>|</mo> <mo>&lt;</mo> <mi>&amp;Delta;</mi> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mo>|</mo> <msub> <mover> <mover> <mi>h</mi> <mo>~</mo> </mover> <mo>&amp;OverBar;</mo> </mover> <mrow> <mi>i</mi> <mo>+</mo> <msub> <mn>9</mn> <mrow> <mn>1</mn> <mo>/</mo> <mn>3</mn> </mrow> </msub> </mrow> </msub> <mo>-</mo> <msub> <mi>H</mi> <mrow> <mn>1</mn> <mo>/</mo> <mn>3</mn> </mrow> </msub> <mo>|</mo> <mo>&lt;</mo> <mi>&amp;Delta;</mi> </mrow>
若是，则整个计算终止，所生成的随机不规则波浪的波面坐标即代表该等级海况的波面坐标；若不是，则令波面方程的计算时间t＝t+Δt，Δt为计算时间步长，并返回步骤3)继续产生随机不规则波浪的波面坐标，直至满足计算终止条件为止。
2.根据权利要求1所述的随机不规则波浪特征参数提取方法，其特征在于：所述步骤1)中，波能谱S(ωj)的公式为：
<mrow> <mi>S</mi> <mrow> <mo>(</mo> <msub> <mi>&amp;omega;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>0.78</mn> <msubsup> <mi>&amp;omega;</mi> <mi>j</mi> <mn>5</mn> </msubsup> </mfrac> <mi>exp</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>-</mo> <mn>3.11</mn> </mrow> <mrow> <msubsup> <mi>&amp;omega;</mi> <mi>j</mi> <mn>4</mn> </msubsup> <msubsup> <mi>H</mi> <mrow> <mn>1</mn> <mo>/</mo> <mn>3</mn> </mrow> <mn>2</mn> </msubsup> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
3.根据权利要求1所述的随机不规则波浪特征参数提取方法，其特征在于：所述步骤3)中，第j个频率等分所对应的波幅ζj、频率ωj和波数kj的计算公式为：
和
式中：Tj为第j个频率等分对应的周期，g为重力加速度g＝9.8；当j＝1时，ω0＝ωmin。
4.根据权利要求3所述的随机不规则波浪特征参数提取方法，其特征在于：所述步骤5)中，波面方程为：
<mrow> <mi>&amp;zeta;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msqrt> <mrow> <mn>2</mn> <mi>S</mi> <mrow> <mo>(</mo> <msub> <mi>&amp;omega;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>&amp;Delta;&amp;omega;</mi> <mi>j</mi> </msub> </mrow> </msqrt> <mi>cos</mi> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mi>j</mi> </msub> <mi>x</mi> <mi> </mi> <mi>cos</mi> <mi>&amp;chi;</mi> <mo>-</mo> <msub> <mi>k</mi> <mi>j</mi> </msub> <mi>y</mi> <mi> </mi> <mi>sin</mi> <mi>&amp;chi;</mi> <mo>-</mo> <msub> <mi>&amp;omega;</mi> <mi>j</mi> </msub> <mi>t</mi> <mo>+</mo> <msub> <mi>&amp;epsiv;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
5.根据权利要求1所述的随机不规则波浪特征参数提取方法，其特征在于：所述步骤6)中，提取随机不规则波浪的波高的过程为：利用计算机分别提取步骤5)所生成的随机不规则波浪的波面坐标的波峰和波谷，将第i个波峰与相邻第i个波谷之间的垂直距离定义为对应第i个随机不规则波浪的波高并提取该波高
CN201710355031.2A 2017-05-19 2017-05-19 一种随机不规则波浪特征参数提取方法 Active CN107153773B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710355031.2A CN107153773B (zh) 2017-05-19 2017-05-19 一种随机不规则波浪特征参数提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710355031.2A CN107153773B (zh) 2017-05-19 2017-05-19 一种随机不规则波浪特征参数提取方法

Publications (2)

Publication Number Publication Date
CN107153773A true CN107153773A (zh) 2017-09-12
CN107153773B CN107153773B (zh) 2018-07-13

Family

ID=59794272

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710355031.2A Active CN107153773B (zh) 2017-05-19 2017-05-19 一种随机不规则波浪特征参数提取方法

Country Status (1)

Country Link
CN (1) CN107153773B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108805100A (zh) * 2018-06-25 2018-11-13 大连理工大学 二维海况下基于波群特性的畸形波短期预报方法及应用
CN113703020A (zh) * 2021-01-30 2021-11-26 东南大学 一种海面小型目标升沉运动测量方法
CN115526052A (zh) * 2022-10-08 2022-12-27 南京航空航天大学 基于等分频率-能量-方向角策略不规则波浪谱解析方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080094940A1 (en) * 1998-08-04 2008-04-24 Teledyne Rd Instruments, Inc. System and method for measuring wave directional spectrum and wave height
CN104021288A (zh) * 2014-06-04 2014-09-03 中国石油化工集团公司 用于导管架平台频谱疲劳分析的基本波确定方法
CN106053016A (zh) * 2016-07-25 2016-10-26 中国海洋大学 一种推板造随机波的方法
CN106682391A (zh) * 2016-11-23 2017-05-17 大连理工大学 基于p‑m谱和jonswap谱的实测风浪情况筛选和理论频谱拟合方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080094940A1 (en) * 1998-08-04 2008-04-24 Teledyne Rd Instruments, Inc. System and method for measuring wave directional spectrum and wave height
CN104021288A (zh) * 2014-06-04 2014-09-03 中国石油化工集团公司 用于导管架平台频谱疲劳分析的基本波确定方法
CN106053016A (zh) * 2016-07-25 2016-10-26 中国海洋大学 一种推板造随机波的方法
CN106682391A (zh) * 2016-11-23 2017-05-17 大连理工大学 基于p‑m谱和jonswap谱的实测风浪情况筛选和理论频谱拟合方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108805100A (zh) * 2018-06-25 2018-11-13 大连理工大学 二维海况下基于波群特性的畸形波短期预报方法及应用
CN108805100B (zh) * 2018-06-25 2021-08-17 大连理工大学 二维海况下基于波群特性的畸形波短期预报方法及应用
CN113703020A (zh) * 2021-01-30 2021-11-26 东南大学 一种海面小型目标升沉运动测量方法
CN113703020B (zh) * 2021-01-30 2023-12-05 东南大学 一种海面小型目标升沉运动测量方法
CN115526052A (zh) * 2022-10-08 2022-12-27 南京航空航天大学 基于等分频率-能量-方向角策略不规则波浪谱解析方法

Also Published As

Publication number Publication date
CN107153773B (zh) 2018-07-13

Similar Documents

Publication Publication Date Title
CN107153773A (zh) 一种随机不规则波浪特征参数提取方法
Tatum et al. Wave–current interaction effects on tidal stream turbine performance and loading characteristics
Lora et al. Design and performance of a nozzle-type rainfall simulator for landslide triggering experiments
Sergeeva et al. Rogue waves, rogue events and extreme wave kinematics in spatio-temporal fields of simulated sea states
Gopalan et al. A coupled mesoscale–microscale framework for wind resource estimation and farm aerodynamics
CN105069476B (zh) 基于两阶段集成学习的风电异常数据识别方法
CN105046705B (zh) 一种基于模糊理论的裂纹边缘检测方法
Tian et al. An experimental investigation on the aeromechanics and wake interferences of wind turbines sited over complex terrain
CN107944349A (zh) 基于粒子群‑支持向量机的gis局部放电类型模式识别方法
Hardy et al. A maximum capture width tracking controller for ocean wave energy converters in irregular waves
CN106202629A (zh) 一种跨海桥梁基础所受波浪力的计算方法
Kuo et al. Hydrodynamic characteristics of Oscillating Water Column caisson breakwaters
CN105375127B (zh) 一种变电站中打接地深井方式的降低接地电阻方法
CN104268662B (zh) 一种基于分步优化分位数回归的沉降预测方法
CN108984900B (zh) 一种基于深水区风浪条件的港内波高变化分析方法
CN104481511B (zh) 一种用载荷力测量抽油机井环空动液面及工况诊断方法
CN104849750B (zh) 基于目标波波形分析的核电楼层谱人工波拟合方法
Iungo et al. Reduced order model for optimization of power production from a wind farm
CN105093280A (zh) 表层模型对地震数据影响的低频与高频成分的分解方法
Eriksson et al. Analysis of long distance wakes of Horns Rev I using actuator disc approach
CN104729786B (zh) 一种孔内灌浆压力的动态测量方法
Castellani et al. Analyzing wind turbine flow interaction through vibration data
CA3139663C (en) Evaluating hydraulic fracturing breakdown effectiveness
CN104990660B (zh) 一种应用于浮式平台的内波流载荷获取方法
Su et al. Comparative study on collection performance of two back-end methods of double-row hydraulic sluicing structure in deep-sea mining

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