CN116805028A - 一种基于浮体运动响应的波面反演方法及系统 - Google Patents
一种基于浮体运动响应的波面反演方法及系统 Download PDFInfo
- Publication number
- CN116805028A CN116805028A CN202311036372.5A CN202311036372A CN116805028A CN 116805028 A CN116805028 A CN 116805028A CN 202311036372 A CN202311036372 A CN 202311036372A CN 116805028 A CN116805028 A CN 116805028A
- Authority
- CN
- China
- Prior art keywords
- floating body
- motion
- heave
- wave surface
- freedom
- 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
- 238000007667 floating Methods 0.000 title claims abstract description 156
- 230000004044 response Effects 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000006073 displacement reaction Methods 0.000 claims abstract description 40
- 238000005316 response function Methods 0.000 claims description 30
- 239000011159 matrix material Substances 0.000 claims description 28
- 239000013598 vector Substances 0.000 claims description 27
- 238000012512 characterization method Methods 0.000 claims description 18
- 230000001364 causal effect Effects 0.000 claims description 17
- 230000014509 gene expression Effects 0.000 claims description 12
- 230000009471 action Effects 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 6
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 6
- RFHAOTPXVQNOHP-UHFFFAOYSA-N fluconazole Chemical compound C1=NC=NN1CC(C=1C(=CC(F)=CC=1)F)(O)CN1C=NC=N1 RFHAOTPXVQNOHP-UHFFFAOYSA-N 0.000 claims description 3
- 238000011084 recovery Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 2
- 238000007665 sagging Methods 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 5
- 230000001052 transient effect Effects 0.000 abstract description 4
- 230000007547 defect Effects 0.000 abstract description 3
- 230000001788 irregular Effects 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 2
- 230000002706 hydrostatic effect Effects 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Complex Calculations (AREA)
Abstract
本发明属于海洋工程技术领域,提供一种基于浮体运动响应的波面反演方法及系统。基于浮体运动响应信号初始值,确定由非零初始条件引发的浮体垂荡自由度运动信号的自振响应分量,将自振响应分量从浮体垂荡自由度运动信号中剔除,得到零初始条件下浮体垂荡位移序列;基于零初始条件下浮体垂荡位移序列,以及浮体的极点‑留数表征模型,计算得到输入波面信号的傅里叶系数;基于输入波面信号的傅里叶系数,进行傅里叶逆变换,得到入射波面时程。本发明提供的波面反演方法和系统以拉普拉斯域极留数理论框架为基础,充分考虑了浮体实测运动响应的固有瞬态特征,突破了传统技术将浮体响应全部假设为稳态响应的缺陷,可提高波面反演的精度。
Description
技术领域
本发明涉及海洋工程技术领域,尤其涉及一种基于浮体运动响应的波面反演方法及系统。
背景技术
海洋动力环境,特别是海浪环境要素,对于海洋工程的设计和运维至关重要,波面运动是海洋动力环境的直接反应,通过对波面运动的测量,可以有效反应海洋运动规律,指导海洋工程建设和海上设施运维。波面运动的测量通常需要借助海上目标的运动进行反演。
浮标是目前国际上最为主流和通用的海浪环境观测手段,观测浮标通常具有良好的随波性,通过记录浮标自身垂荡运动实现波面观测。近年来,人们逐渐将浮标观测海浪技术推广延伸至大型油气平台和船舶等一般浮体,利用浮体自身响应进行波浪反演,为平台结构实时健康监测、船舶航行辅助决策、海上施工作业等提供环境数据支撑。
然而,目前这一技术在实际应用中尚处于起步探索阶段,特别是在不规则波面时历反演方面仍存在较大不足。主要表现在,现有技术建立在浮体运动响应稳态假定的频域理论框架基础上,忽略了浮体动力响应本身的固有瞬态特性,严重影响了波面的识别精度,成为目前制约浮体响应反演波面的一大重要技术瓶颈。
现有技术中,尚缺少可精确进行波面反演的算法。浮体漂浮在海面上,其运动是对波面运动的直接反应,浮体的运动参数采集技术相对成熟。利用浮体运动参数进行波面运动的反演,将能获得较好的反演结果。
发明内容
本发明的目的在于解决以上技术问题,提供一种波面识别精度高的基于浮体运动响应的波面反演方法及系统,该方法旨在借助海面浮体目标的运动响应,充分利用其瞬态特性,对波面的运动进行更加精确的反演,该系统用于实现该波面反演方法。
为实现上述目的,本发明采用的技术方案是:
一种基于浮体运动响应的波面反演方法,包括以下步骤:
S1:基于浮体运动响应信号初始值,确定由非零初始条件引发的浮体垂荡自由度运动信号的自振响应分量,将所述自振响应分量从浮体垂荡自由度运动信号/>中剔除,得到零初始条件下浮体垂荡位移序列/>;
S2:基于零初始条件下浮体垂荡位移序列,以及浮体的极点-留数表征模型,计算得到输入波面信号的傅里叶系数/>;
S3:基于输入波面信号的傅里叶系数,进行傅里叶逆变换,得到入射波面时程。
本发明一些实施例中,步骤S1的执行步骤进一步包括:
S11:基于浮体六自由度运动位移信号和浮体六自由度速度信号,确定浮体在第个自由度的运动位移向量和浮体在第/>个自由度的运动速度向量,其中:
浮体六自由度运动位移信号:;
浮体六自由度速度信号:;
浮体在第个自由度的运动位移向量:/>=/>;
浮体在第个自由度的运动速度向量:/>=/>;
其中:下标1, 2, 3, 4, 5, 6分别代表浮体纵荡、横荡、垂荡、横摇、纵摇和首摇运动自由度,1, 2, 3, 4, 5, 6,/>,/>=0,1,……,/>-1,/>是时间,/>是时间间隔,/>为时间点序号,/>是时间点个数,上标/>代表向量转置;
令=0,确定浮体六自由度运动初始位移和浮体六自由度运动初始速度,其中:
浮体六自由度运动初始位移:
浮体六自由度运动初始速度:
S12:基于浮体六自由度运动初始位移和浮体六自由度运动初始速度/>,采用四阶-龙格库塔算法求解浮体自振运动控制方程:
;
基于浮体自振运动控制方程,得到非零初始条件引发的浮体垂荡自由度运动信号的自振响应分量;
其中,为浮体质量矩阵,/>为浮体附加质量矩阵,/>为浮体静水回复刚度矩阵,为浮体时延函数,/>为时间哑变量;
S13:将所述自振响应分量从浮体垂荡自由度信号/>中剔除,得到零初始条件下浮体垂荡位移序列/>:
。
本发明一些实施例中,步骤S2的执行步骤进一步包括:
S21:确定待求入射波面时程的傅里叶级数表征模型:
;
其中,是欧拉数,/>为指数函数,/>为虚数,/>为离散频率,/>为频率间隔,/>为傅里叶系数,/>为频率点序号,/>是频率点个数;
S22:对上述傅里叶级数表征模型函数进行拉普拉斯变换,得到波面在拉普拉斯域的极点-留数表征函数:
;
其中,为波面时程函数/>的拉普拉斯变换,/>为拉普拉斯域变量;
S23:计算波浪作用下浮体垂荡频响函数,在拉普拉斯域内构建浮体因果化波浪-垂荡运动系统传递函数的极点-留数表征模型:
对浮体垂荡频响函数进行傅里叶逆变换,获得浮体垂荡脉冲响应函数/>:
;
其中,是圆周率,/>是圆频率,/>是时间;
基于浮体垂荡脉冲响应函数,计算得到因果化浮体垂荡脉冲响应函数/>:
;
其中,为因果化时间,满足当/>时,/>取值为零;
对进行复指数分解,将/>表示成一系列有限数量复指数函数之和:
;
对复指数分解后的所述因果化浮体垂荡脉冲响应函数进行拉普拉斯变换,得到所述浮体垂荡脉冲响应函数在拉普拉斯域内的极点-留数表征模型;
;
其中,为/>的拉普拉斯变换,/>为极点,/>为留数,/>为极点和留数的序号,为极点和留数的个数;
S24:构建零初始条件下入射波面作用引起的浮体垂荡位移序列计算模型,包括:
将上述得到的和/>相乘,得到拉普拉斯域内垂荡响应/>:
;
将写成极点-留数形式:
;
其中:为对应极点/>的留数,/>为对应极点/>的留数,其表达式分别为:
;
;
对进行拉普拉斯逆变换,得到时域函数/>表达式:
;
其中,函数:
;
令,可得零初始条件下浮体垂荡响应序列:
;
S25:构建基于零初始条件下浮体垂荡自由度运动信号自振响应分量的波面傅里叶系数/>的求解模型,包括:
使用的后/>个数据,将/>表示成矩阵形式:
;
其中,为垂荡响应列向量;
;
为傅里叶系数列向量:
;
为系数矩阵:
;
其中,为列向量,表达式为:
;
与/>满足如下关系:
;
;
其中,符号代表复共轭运算;
将矩阵表达写成如下显式表达:
;
令,/>,将/>的表达式改写成:
;
其中,和/>为/>的实部和虚部,/>和/>为/>的实部和虚部,/>表示取实部运算,/>表示取虚部运算;
将的表达式进一步写成矩阵形式:
;
其中,为系数矩阵:
;
为待求系数列向量:
;
求解关于的矩阵方程,可计算得到系数矩阵/>;
基于系数矩阵的计算结果,可以得到波面的傅里叶系数/>。
本发明一些实施例中,步骤S3的执行步骤包括:基于计算得到的波面傅里叶系数,计算入射波面时程:
。
本发明一些实施例中,步骤S23中,通过水动力软件,分析获得波浪作用下浮体垂荡频响函数。
本发明一些实施例中,浮体为浮标或船舶。
本发明一些实施例进一步提供一种波面反演系统,用于执行以上所述的波面反演方法。
本发明提供的基于浮体运动响应的波面反演方法及系统,其有益效果在于:
本发明提出一种高精度的利用浮体运动响应识别不规则波面时历的新方法。该方法建立在全新的极留数理论框架之上,充分考虑了浮体实测运动响应的固有瞬态特征,突破了传统技术将浮体响应全部假设为稳态响应的缺陷与不足,极大地提高了波面反演精度。
本发明提出的波面反演系统,可被搭载设置在海面浮体或陆地端,配置灵活,可实时进行数据反演计算,查看计算结果。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为发明基于浮体运动响应的波面反演方法流程图;
图2为本发明实施例公开的不规则入射波面时程图;
图3为本发明的实施例公开的FPSO垂荡运动响应时程;
图4为本发明的实施例公开的非零垂荡初始位移条件下FPSO垂荡自振响应分量时程;
图5为本发明的实施例公开的零初始条件下FPSO垂荡响应分量时程;
图6为本发明的实施例公开的FPSO垂荡脉冲响应函数以及因果化垂荡脉冲响应函数;图7a为本发明方法得到的入射波面傅里叶系数的实部结果;
图7b为本发明方法得到的入射波面傅里叶系数的虚部结果;
图8为本发明方法得到的入射波面时程结果与真实值的对比图。
具体实施方式
为了使本发明所要解决的技术问题、技术方案及有益效果更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明第一实施例提供一种基于浮体运动响应的波面反演方法,其中,浮体被设置在海面上,可以为浮标、船舶等。
在本发明一些示意性实施例中,波面反演方法具体包括以下步骤。
S1:基于浮体运动响应信号初始值,确定由非零初始条件引发的浮体垂荡自由度运动信号的自振响应分量,将所述自振响应分量从浮体垂荡自由度运动信号/>中剔除,得到零初始条件下浮体垂荡位移序列/>。
本发明一些实施例中,步骤S1的执行步骤进一步包括:
S11:基于浮体六自由度运动位移信号和浮体六自由度速度信号,确定浮体在第个自由度的运动位移向量和浮体在第/>个自由度的运动速度向量,其中:
浮体六自由度运动位移信号:;
浮体六自由度速度信号:;
浮体在第个自由度的运动位移向量:/>=/>;
浮体在第个自由度的运动速度向量:/>=/>;
其中:下标1, 2, 3, 4, 5, 6分别代表浮体纵荡、横荡、垂荡、横摇、纵摇和首摇运动自由度,1, 2, 3, 4, 5, 6,/>,/>=0,1,……,/>-1,/>是时间,/>是时间间隔,/>为时间点序号,/>是时间点个数,上标/>代表向量转置;
令 =0,确定浮体六自由度运动初始位移和浮体六自由度运动初始速度,其中:
浮体六自由度运动初始位移:;
浮体六自由度运动初始速度:;
S12:基于浮体六自由度运动初始位移和浮体六自由度运动初始速度/>,采用四阶-龙格库塔算法求解浮体自振运动控制方程:
;
基于浮体自振运动控制方程,得到非零初始条件引发的浮体垂荡自由度运动信号的自振响应分量;
其中,为浮体质量矩阵,/>为浮体附加质量矩阵,/>为浮体静水回复刚度矩阵,为浮体时延函数,/>为时间哑变量;
S13:将所述自振响应分量从浮体垂荡自由度运动信号/>中剔除,得到零初始条件下浮体垂荡位移序列/>:
。
S2:基于零初始条件下浮体垂荡位移序列,以及浮体的极点-留数表征模型,计算得到输入波面信号的傅里叶系数/>。
本发明一些实施例中,步骤S2的执行步骤进一步包括:
S21:确定待求波面时程的傅里叶级数表征模型:
;
其中,是欧拉数,/>为指数函数,/>为虚数,/>为离散频率,/>为频率间隔,/>为傅里叶系数,/>为频率点序号,/>是频率点个数;
S22:对上述傅里叶级数表征模型函数进行拉普拉斯变换,得到波面在拉普拉斯域的极点-留数表征函数:
;
其中,为波面时程/>的拉普拉斯变换,/>为拉普拉斯域变量;
S23:计算波浪作用下浮体垂荡频响函数,在拉普拉斯域内构建浮体因果化波浪-垂荡运动系统传递函数的极点-留数表征模型:
对浮体垂荡频响函数进行傅里叶逆变换,获得浮体垂荡脉冲响应函数/>:
;
其中,是圆周率,/>是圆频率,/>是时间;
基于浮体垂荡脉冲响应函数,计算得到因果化浮体垂荡脉冲响应函数/>:
;
其中,为因果化时间,满足当/>时,/>取值为零;
对进行复指数分解,将/>表示成一系列有限数量复指数函数之和:
;
对复指数分解后的所述因果化浮体垂荡脉冲响应函数进行拉普拉斯变换,得到所述浮体垂荡脉冲响应函数在拉普拉斯域内的极点-留数表征模型;
;
其中,为/>的拉普拉斯变换,/>为极点,/>为留数,/>为极点和留数的序号,为极点和留数的个数;
S24:构建零初始条件下入射波面作用引起的浮体垂荡位移序列计算模型,包括:
将上述得到的和/>相乘,得到拉普拉斯域内垂荡响应/>:
;
将写成极点-留数形式:
;
其中:为对应极点/>的留数,/>为对应极点/>的留数,其表达式分别为:
;
;
对进行拉普拉斯逆变换,得到时域函数/>表达式:
;
其中,函数:
;
令,可得零初始条件下浮体垂荡响应序列:
;
S25:构建基于零初始条件下浮体垂荡自由度运动信号自振响应分量的波面傅里叶系数/>的求解模型,包括:
使用的后/>个数据,将/>表示成矩阵形式:
;
其中,为垂荡响应列向量;
;
为傅里叶系数列向量:
;
为系数矩阵:
;
其中,为列向量,表达式为:
;
与/>满足如下关系:
;
;
其中,符号代表复共轭运算;
将矩阵表达写成如下显式表达:
;
令,/>,将/>的表达式改写成:
;
其中,和/>为/>的实部和虚部,/>和/>为/>的实部和虚部,/>表示取实部运算,/>表示取虚部运算;
将的表达式进一步写成矩阵形式:
;
其中,为系数矩阵:
;
为待求系数列向量:
;
求解关于的矩阵方程,可计算得到系数矩阵/>;
基于系数矩阵的计算结果,可以得到波面的傅里叶系数/>。
本发明一些实施例中,步骤S23中,通过水动力软件,分析获得波浪作用下浮体垂荡频响函数。
S3:基于计算得到的波面傅里叶系数,计算入射波面时程。
具体的,基于计算得到的波面傅里叶系数,由下式可得入射波面时程:
。
本发明采用如下实际应用来验证本发明方法的有效性。
本发明的实施例选取DNV船级社商业水动力软件SESAM的一个六自由度FPSO(浮式生产储油船)数值算例,该FPSO长224.12 m,宽29.7 m,吃水16 m;重1.257×108Kg,横摇、纵摇和艏摇惯性半径分别为16.24 m,54.14 m和54.90 m。重心坐标(−4.18 m,0,−2.51 m),其中坐标系平面位于静水面,Z轴正向向上。
图2为时长51.2 s、时间间隔0.0256 s的不规则入射波面时程。考虑垂荡初始非零位移m,通过数值仿真得到不规则波零度角入射下FPSO的垂荡响应时程,见图3。本实施例将利用零度入射波浪作用下FPSO的垂荡响应识别不规则波面时程。
基于给定的FPSO垂荡响应非零初始位移条件m,通过龙格-库塔方法求解自振响应控制方程,得到垂荡非零初始位移引发的FPSO垂荡自振响应时程,见图4。进一步地,将自振响应分量从垂荡响应中剔除,得到零初始条件下FPSO垂荡响应分量,见图5。
下一步需要计算FPSO系统的极点和留数。对SESAM软件计算得到的该FPSO系统垂荡响应离散频响函数(频率范围0.03-3.12 rad/s,频率间隔0.03 rad/s)进行傅里叶逆变换,得到垂荡脉冲响应函数,见图6,观察图6可确定因果化时间/>s,从而得到因果化垂荡脉冲响应函数/>,/>曲线也见图6。对信号/>进行复指数分解,得到9对复共轭的因果化FPSO波浪-垂荡系统的极点与留数,见表1。
表1 FPSO波浪-垂荡系统的极点与留数值
基于图5得到零初始条件下FPSO的垂荡响应时程,以及表1得到的FPSO波浪-垂荡系统的极点与留数结果,利用本发明进行入射波面识别。
图7a和图7b分别为本发明得到的入射波面傅里叶系数的实部和虚部的识别结果,基于图7a和图7b得到的傅里叶系数,通过傅里叶逆变换得到波面时程。
图8展示了识别得到的入射波面时程与给定真实入射波面时程的对比结果,可以看出二者吻合较好,证明了本发明方法的有效性和高识别精度。
本发明第二实施例进一步提供一种波面反演系统,该系统用于实现实施例1所述的浮体运动响应的波面反演方法。该波面反演系统基于数据处理器而实现,可被搭载在浮体上,或设置在陆地端。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种基于浮体运动响应的波面反演方法,其特征在于,包括以下步骤:
S1:基于浮体运动响应信号初始值,确定由非零初始条件引发的浮体垂荡自由度运动信号的自振响应分量,将所述自振响应分量从浮体垂荡自由度运动信号/>中剔除,得到零初始条件下浮体垂荡位移序列/>;
S2:基于零初始条件下浮体垂荡位移序列,以及浮体的极点-留数表征模型,计算得到输入波面信号的傅里叶系数/>;
S3:基于输入波面信号的傅里叶系数,进行傅里叶逆变换,得到入射波面时程/>。
2.如权利要求1所述的基于浮体运动响应的波面反演方法,其特征在于,步骤S1的执行步骤进一步包括:
S11:基于浮体六自由度运动位移信号和浮体六自由度速度信号,确定浮体在第个自由度的运动位移向量和浮体在第/>个自由度的运动速度向量,其中:
浮体六自由度运动位移信号:;
浮体六自由度速度信号:;
浮体在第个自由度的运动位移向量:/>=/>;
浮体在第个自由度的运动速度向量:/>=/>;
其中:下标1, 2, 3, 4, 5, 6分别代表浮体纵荡、横荡、垂荡、横摇、纵摇和首摇运动自由度,1, 2, 3, 4, 5, 6, />,/>=0,1,……,/>-1,/>是时间,/>是时间间隔,/>为时间点序号,/>是时间点个数,上标/>代表向量转置;
令=0,确定浮体六自由度运动初始位移和浮体六自由度运动初始速度,其中:
浮体六自由度运动初始位移:
浮体六自由度运动初始速度:
S12:基于浮体六自由度运动初始位移和浮体六自由度运动初始速度/>,采用四阶-龙格库塔算法求解浮体自振运动控制方程:
;
基于浮体自振运动控制方程,得到非零初始条件引发的浮体垂荡自由度运动信号的自振响应分量;
其中,为浮体质量矩阵,/>为浮体附加质量矩阵,/>为浮体静水回复刚度矩阵,/>为浮体时延函数,/>为时间哑变量;
S13:将所述自振响应分量从浮体垂荡自由度运动信号/>中剔除,得到零初始条件下浮体垂荡位移序列/>:
。
3.如权利要求1所述的基于浮体运动响应的波面反演方法,其特征在于,步骤S2的执行步骤进一步包括:
S21:确定待求入射波面时程的傅里叶级数表征模型:
;
其中,是欧拉数,/>为指数函数,/>为虚数,/>为离散频率,/>为频率间隔,/>为傅里叶系数,/>为频率点序号,/>是频率点个数;
S22:对上述傅里叶级数表征模型函数进行拉普拉斯变换,得到波面在拉普拉斯域的极点-留数表征函数:
;
其中,为波面时程/>的拉普拉斯变换,/>为拉普拉斯域变量;
S23:计算波浪作用下浮体垂荡频响函数,在拉普拉斯域内构建浮体因果化波浪-垂荡运动系统传递函数的极点-留数表征模型:
对浮体垂荡频响函数进行傅里叶逆变换,获得浮体垂荡脉冲响应函数/>;
基于浮体垂荡脉冲响应函数,计算得到因果化浮体垂荡脉冲响应函数/>;
对进行复指数分解,将/>表示成一系列有限数量复指数函数之和:
;
对复指数分解后的所述因果化浮体垂荡脉冲响应函数进行拉普拉斯变换,得到所述浮体垂荡脉冲响应函数在拉普拉斯域内的极点-留数表征模型:
;
其中,为/>的拉普拉斯变换,/>为极点,/>为留数,/>为极点和留数的序号,/>为极点和留数的个数;
S24:基于和/>,构建零初始条件下入射波面作用引起的浮体垂荡位移序列计算模型;
S25:构建基于零初始条件下浮体垂荡自由度运动信号自振响应分量的波面傅里叶系数/>的求解模型,包括:
使用的后/>个数据,将/>表示成矩阵形式:
;
为垂荡响应列向量;求解关于/>的矩阵方程,可计算得到系数矩阵/>;基于系数矩阵/>的计算结果,得到波面的傅里叶系数/>。
4.如权利要求3所述的基于浮体运动响应的波面反演方法,其特征在于,步骤S24的执行步骤包括:
将和/>相乘,得到拉普拉斯域内垂荡响应/>:
;
将写成极点-留数形式:
;
其中:为对应极点/>的留数,/>为对应极点/>的留数,其表达式分别为:
;
;
对进行拉普拉斯逆变换,得到时域函数/>表达式:
;
其中,函数:
;
令,可得零初始条件下浮体垂荡位移序列:
。
5.如权利要求3所述的基于浮体运动响应的波面反演方法,其特征在于,步骤S25中:
;
为傅里叶系数列向量:
;
为系数矩阵:
;
其中,为列向量,表达式为:
;
与/>满足如下关系:
;
;
其中,符号代表复共轭运算;
将矩阵表达写成如下显式表达:
;
令,/>,将/>的表达式改写成:
;
其中,和/>为/>的实部和虚部,/>和/>为/>的实部和虚部,/>表示取实部运算,表示取虚部运算;
将的表达式进一步写成矩阵形式:
;
其中,为系数矩阵:
;
为待求系数列向量:
;
求解关于的矩阵方程,可计算得到系数矩阵/>;
基于系数矩阵的计算结果,可以得到波面的傅里叶系数/>。
6.如权利要求5所述的基于浮体运动响应的波面反演方法,其特征在于,步骤S23中:
浮体垂荡脉冲响应函数:
;
其中,是圆周率,/>是圆频率,/>是时间;
因果化浮体垂荡脉冲响应函数:
;
其中,为因果化时间,满足当/>时,/>取值为零。
7.如权利要求3或4或5或6所述的基于浮体运动响应的波面反演方法,其特征在于,步骤S3的执行步骤包括:基于计算得到的波面傅里叶系数,计算入射波面时程:
。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311036372.5A CN116805028B (zh) | 2023-08-17 | 2023-08-17 | 一种基于浮体运动响应的波面反演方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311036372.5A CN116805028B (zh) | 2023-08-17 | 2023-08-17 | 一种基于浮体运动响应的波面反演方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116805028A true CN116805028A (zh) | 2023-09-26 |
CN116805028B CN116805028B (zh) | 2023-12-08 |
Family
ID=88080814
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311036372.5A Active CN116805028B (zh) | 2023-08-17 | 2023-08-17 | 一种基于浮体运动响应的波面反演方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116805028B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117648829A (zh) * | 2024-01-29 | 2024-03-05 | 中国海洋大学 | 非线性海工结构动力响应的拉普拉斯域计算方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002013923A (ja) * | 2000-06-29 | 2002-01-18 | Kawasaki Heavy Ind Ltd | 大型浮体による波浪情報計測方法と装置 |
JP2008197018A (ja) * | 2007-02-14 | 2008-08-28 | National Maritime Research Institute | 複数の隣接浮体模型の波浪中動揺応答特性推定方法 |
US20100046599A1 (en) * | 2006-11-22 | 2010-02-25 | Dong-Kyoo Kim | Apparatus and method for acquiring initial coefficient of decision feedback equalizer using fast fourier transform |
US8618961B1 (en) * | 2012-03-12 | 2013-12-31 | Audience, Inc. | Sample rate conversion using infinite impulse response filters |
US20180372494A1 (en) * | 2015-11-20 | 2018-12-27 | Fluid Techno Co., Ltd. | Sea state estimation device and sea state estimation method |
CN112394332A (zh) * | 2020-11-27 | 2021-02-23 | 中国舰船研究设计中心 | 一种基于空间位移的船舶姿态运动成分去除方法 |
CN113673076A (zh) * | 2021-07-09 | 2021-11-19 | 华南理工大学 | 一种适用于海洋浮式结构的动力响应求解方法 |
CN114880619A (zh) * | 2022-05-10 | 2022-08-09 | 中国海洋大学 | 浮式海洋结构随机动力响应解析计算方法 |
CN115292930A (zh) * | 2022-08-04 | 2022-11-04 | 中国海洋大学 | 浮式海洋结构物瞬态动力响应的频域计算方法 |
CN116227242A (zh) * | 2023-05-09 | 2023-06-06 | 中国海洋大学 | 基于白噪声滤波的随机不规则波浪模拟方法 |
CN116226728A (zh) * | 2023-05-09 | 2023-06-06 | 中国海洋大学 | 基于单次规则波激励的浮式结构频响函数识别方法 |
WO2023102960A1 (zh) * | 2021-12-06 | 2023-06-15 | 中国人民解放军国防科技大学 | 一种实验室波浪智能测量系统 |
CN116304647A (zh) * | 2023-05-19 | 2023-06-23 | 中国海洋大学 | 一种基于瞬态响应的浮式结构频响函数识别方法 |
-
2023
- 2023-08-17 CN CN202311036372.5A patent/CN116805028B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002013923A (ja) * | 2000-06-29 | 2002-01-18 | Kawasaki Heavy Ind Ltd | 大型浮体による波浪情報計測方法と装置 |
US20100046599A1 (en) * | 2006-11-22 | 2010-02-25 | Dong-Kyoo Kim | Apparatus and method for acquiring initial coefficient of decision feedback equalizer using fast fourier transform |
JP2008197018A (ja) * | 2007-02-14 | 2008-08-28 | National Maritime Research Institute | 複数の隣接浮体模型の波浪中動揺応答特性推定方法 |
US8618961B1 (en) * | 2012-03-12 | 2013-12-31 | Audience, Inc. | Sample rate conversion using infinite impulse response filters |
US20180372494A1 (en) * | 2015-11-20 | 2018-12-27 | Fluid Techno Co., Ltd. | Sea state estimation device and sea state estimation method |
CN112394332A (zh) * | 2020-11-27 | 2021-02-23 | 中国舰船研究设计中心 | 一种基于空间位移的船舶姿态运动成分去除方法 |
CN113673076A (zh) * | 2021-07-09 | 2021-11-19 | 华南理工大学 | 一种适用于海洋浮式结构的动力响应求解方法 |
WO2023102960A1 (zh) * | 2021-12-06 | 2023-06-15 | 中国人民解放军国防科技大学 | 一种实验室波浪智能测量系统 |
CN114880619A (zh) * | 2022-05-10 | 2022-08-09 | 中国海洋大学 | 浮式海洋结构随机动力响应解析计算方法 |
CN115292930A (zh) * | 2022-08-04 | 2022-11-04 | 中国海洋大学 | 浮式海洋结构物瞬态动力响应的频域计算方法 |
CN116227242A (zh) * | 2023-05-09 | 2023-06-06 | 中国海洋大学 | 基于白噪声滤波的随机不规则波浪模拟方法 |
CN116226728A (zh) * | 2023-05-09 | 2023-06-06 | 中国海洋大学 | 基于单次规则波激励的浮式结构频响函数识别方法 |
CN116304647A (zh) * | 2023-05-19 | 2023-06-23 | 中国海洋大学 | 一种基于瞬态响应的浮式结构频响函数识别方法 |
Non-Patent Citations (4)
Title |
---|
SUN, JINWEI等: "Laplace domain method for computing the transient response of floating structures to irregular waves", OCEAN ENGINEERING * |
曾东;叶恒奎;孙江龙;: "箱式浮体在波浪中的运动分析", 中国舰船研究, no. 1 * |
王华洁;孙金伟;邵萌;范秀涛;: "频域内海洋资料浮标垂荡运动特性研究", 山东科学, no. 06 * |
范秀涛 等: "大型海洋资料浮标的运动性能研究", 海洋技术学报 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117648829A (zh) * | 2024-01-29 | 2024-03-05 | 中国海洋大学 | 非线性海工结构动力响应的拉普拉斯域计算方法 |
CN117648829B (zh) * | 2024-01-29 | 2024-05-03 | 中国海洋大学 | 非线性海工结构动力响应的拉普拉斯域计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN116805028B (zh) | 2023-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116805028B (zh) | 一种基于浮体运动响应的波面反演方法及系统 | |
CN109345875A (zh) | 一种提高船舶自动识别系统测量精度的估计方法 | |
CN101707016A (zh) | 航海模拟器用船舶六自由度运动数学模型的建立方法 | |
CN107330164B (zh) | 一种三体船纵向运动模型辨识方法 | |
CN107179693A (zh) | 基于Huber估计的鲁棒自适应滤波和状态估计方法 | |
WO2023216946A1 (zh) | 漂移速度预测方法、装置、计算设备及存储介质 | |
CN112883488A (zh) | 一种基于船波相对运动的砰击载荷预报方法 | |
Jiao et al. | Predictions of wave induced ship motions and loads by large-scale model measurement at sea and numerical analysis | |
CN115421483A (zh) | 一种无人船操纵运动预报方法 | |
CN113468658B (zh) | 一种风浪中船舶失速系数的直接计算方法 | |
Zhou et al. | Nonparametric modeling of ship maneuvering motions in calm water and regular waves based on R-LSTM hybrid method | |
CN116202492A (zh) | 一种基于高低频分离的随船波浪测量方法 | |
CN110937082A (zh) | 一种基于随机风场和海浪的船舶倾覆风险测试方法 | |
Mu et al. | Time-varying encounter angle trajectory tracking control of unmanned surface vehicle based on wave modeling | |
CN116738124B (zh) | 浮式结构运动响应信号端点瞬态效应消除方法 | |
CN115409160B (zh) | 一种基于定深数据的全海深温度剖面反演方法及系统 | |
CN118052170B (zh) | 一种基于航行体水动力计算内波传播速度的方法、系统、设备及介质 | |
Islam et al. | Estimation of linear hydrodynamic derivatives of a VLCC using static drift simulation | |
Yu et al. | Research on ship coupled motion modeling in waves based on KT maneuvering response model | |
Dejhalla et al. | A Review of the state-of-the art in marine hydrodynamics | |
ZHANG et al. | Control system of ROV for underwater safety inspection and operation | |
Zhang | Stability and chaos analysis of nonlinear roll motion of trimaran ship with variable lay-out under wind and waves | |
Bai et al. | V-shaped wing design and hydrodynamic analysis based on moving base for recovery AUV | |
Cheng et al. | Depth control analysis of Smartfloat using sea trail data and parameter identification | |
Duan et al. | Maneuvering Prediction of a Twin-Propeller Twin-Rudder High-Speed Craft Based on Integrated CFD and Empirical Methods |
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 |