Single point mooring system modal parameter identification based on empirical mode decomposition and time-varying autoregressive model

A method based on empirical mode decomposition (EMD) and time-varying autoregressive (TVAR) model is proposed here to identify the modal parameters of time-varying systems, such as the Floating Production Storage and Offloading (FPSO) single point mooring system. For the EMD–TVAR method, the original signal is decomposed into a finite number of ‘intrinsic mode functions’ (IMFs) by the EMD. Each IMF can be represented as a TVAR model. Then, the time-varying modal parameters i.e., instantaneous frequency (IF) and modal dumping, can be obtained by the basis functions expansion method. The proposed EMD–TVAR method has good results in two experiments compared with the Huang–Hilbert transformation and Short Time Fourier Transform method, and it has been used to analysis the modal parameters of FPSO single point mooring system successfully. The system's time-varying characteristic and its frequency distribution can be known from the modal analysis results.     

Comparative study of broadband damage localization methods applied to test data

This application-oriented paper presents comparison of various broadband frequency based structural damage localization indices using experimental data from a full-scale structure known as the I-40 Bridge. First, three ‘damage-sensitive’ response parameters determined in the time and/or spectral domain, modal domain and wavelet domain are presented for damage localization in the context of a ‘non-model-based’ damage identification approach. Secondly, experimental modal data (namely natural frequencies, mode shapes and modal damping) obtained from this full-scale bridge subjected to various damage conditions is employed for assessment of the robustness of these methods on ‘real-world’ applications. Consequently, the results obtained are compared with those obtained from existing resonance frequency based damage identification methods. The results obtained demonstrate the improved capability of the broadband methods to localize damage in a full-scale structure despite sparse modal information and limited measurement grid points.     

Hydroelastic analysis of pontoon-type circular VLFS with an attached submerged plate

This paper is concerned with the hydroelastic analysis of a pontoon-type, circular, very large floating structure (VLFS) with a horizontal submerged annular plate attached around its perimeter. The coupled fluid–structure interaction problem may be solved by using the modal expansion method in the frequency domain. It involves, firstly, the decomposition of the deflection of a circular Mindlin plate with free edges into vibration modes that are obtained analytically. Then the hydrodynamic diffraction and radiation forces are evaluated by using the eigenfunction expansion matching method which can also be done in an exact manner. The hydroelastic equation of motion is solved by the Rayleigh–Ritz method for the modal amplitudes, and then the modal responses are summed up to obtain the total response. The effectiveness of the attached submerged annular plate in reducing the motion of VLFS has been confirmed by the analysis.     

基于OrcaFlex的钢悬链立管疲劳损伤预报

涡激振动是钢悬链立管疲劳损伤的重要诱导因素之一,因此涡激振动分析在钢悬链立管设计中占据重要地位。OrcaFlex在海洋工程领域功能强大、应用广泛,能够提供两种涡激振动疲劳损伤预报方法,分别是基于模态叠加法的频域预报方法和基于Iwan-Blevins尾流振子模型的时域预报方法。OrcaFlex的涡激振动时域预报模型只考虑了横流向涡激振动造成的疲劳损伤,基于OrcaFlex的时域预报方法开发了双自由度涡激振动尾流振子模型的接口,考虑了浮式平台运动对钢悬链立管的涡激振动激励作用,实现了钢悬链立管双自由度涡激振动疲劳损伤预报。     

Efficient estimation of extreme response of drag-dominated offshore structures by Monte Carlo simulation

The paper describes a novel approach to the problem of estimating the extreme response statistics of a drag-dominated offshore structure exhibiting a pronounced dynamic behaviour when subjected to harsh weather conditions. It is shown that the key quantity for extreme response prediction is the mean upcrossing rate function, which can be simply extracted from simulated response time histories. A commonly adopted practice for obtaining adequate extremes for design purposes requires the execution of 20 or more 3-h time domain analyses for several extreme sea states. For early phase considerations, it would be convenient if extremes of a reasonable accuracy could be obtained based on shorter and fewer simulations. The aim of the work reported in the present paper has therefore been to develop specific methods which make it possible to extract the necessary information about the extreme response from relatively short time histories.The method proposed in this paper opens up the possibility to predict simply and efficiently both short-term and long-term extreme response statistics. The results presented are based on extensive simulation results for the Kvitebjørn jacket structure, in operation on the Norwegian Continental Shelf. Specifically, deck response time histories for different sea states simulated from an MDOF model were used as the basis for our analyses.     

光子晶体光纤带隙及模场的一种综合模拟分析方法

提出直接由MATLAB读入光纤端面图实现复杂结构光子晶体光纤(PCF)模型的快速建立,并综合平面波法(PWM)和频域有限差分(FDFD)法,模拟分析带隙型光子晶体光纤(PBG-PCF)的带隙和模场分布.利用PWM计算得到了在PBG-PCF中传输的光波频率及模式有效折射率范围;基于FDFD法在给定波长及模式有效折射率范围情况下,模拟得到了PBG-PCF中可能存在的模场分布及其他特性.以市售的PBG-PCF为例,验证了数值模拟的正确性,随后系统地分析了结构参数(晶格结构、"原子"占有率、背景材料折射率及"原子     

Prediction of the dynamic response of a slender marine structure under an irregular ocean wave using the NARX-based quadratic Volterra series

The time series of the dynamic response of a slender marine structure was predicted in approximate sense using a truncated quadratic Volterra series. The wave-structure interaction system was identified using the NARX (Nonlinear Autoregressive with Exogenous Input) technique, and the network parameters were determined through supervised training using prepared datasets. The dataset used for network training was obtained by nonlinear finite element analysis of the slender marine structure under random ocean waves of white noise. The nonlinearities involved in the analysis were both large deformation of the structure under consideration and the quadratic term of the relative velocity between the water particle and structure in the Morison formula. The linear and quadratic frequency response functions of the given system were extracted using the multi-tone harmonic probing method and the time series of the response of the structure was predicted using the quadratic Volterra series. To check the applicability of the method, the response of a slender marine structure under a realistic ocean wave environment with a given significant wave height and modal period was predicted and compared with the nonlinear time domain simulation results. The predicted time series of the response of structure with quadratic Volterra series successfully captured the slowly varying response with reasonably good accuracy. This method can be used to predict the response of the slender offshore structure exposed to a Morison type load without relying on the computationally expensive time domain analysis, especially for screening purposes.     

Wavelet transform based coherence analysis of freak wave and its impact

The present paper presents the results of a wavelet transform-based coherence analysis of freak wave and its impact. Wavelet transform has been used as a tool in analyzing signals in the time domain as well as in the frequency domain. The analysis was applied to laboratory-generated freak waves. The wavelet transform of the time history of the freak wave and its impact force revealed that a wide range of frequency components were contained in them. The coherence analysis was conducted on the wave and its impact force time histories. The coherence analysis revealed that some high-frequency components were highly correlated with the impact forces. The present study demonstrates that the wavelet transform can be an alternative tool in the analysis of strongly nonlinear freak wave and its impact.     

Curvature of Flexibility Matrix for Damage Identification in Legs of Jacket Platforms

In this paper a new nondestructive damage identification method is introduced. The method based on flexibility matrix can be used to detect and locate structm＇al damage and evaluate the severity of damage in legs of jacket platforms by modal parameters of a structure. With the modal data for only the few lower modes in both the intact and damaged states, the one-dimensional and two-dimensional distributed curvatures can be used to analyze damage location and the severity. Instead of directly comparing the curvatures before and ＇after damage, the method here uses modal parameters only in the damaged structure to detect the damage and it consists of three parts. First, ilexibility matrix is obtained by use of the absolute maximum in each column. Second, because the legs of jacket platforms are the pipe-like structure, the circumferential flexibility curvature matrix is obtained by use of the circular curvature. At last, equivalent curvature ratio is defined and the curve meaning equivalent curvature ratio and the severity of damage relationship for one element is given through the data of damage severity from ten percent to ninety percent by numerical simulation. Many existing damage detection methods need two steps, locate the damage firstly and evaluate the severity of the damage. However, the method present- ed! in this paper can locate and then evaluate the severity of damage at the same time. The numerical analysis results in- dicate that the present method is effective, useful and only need the first and the second mode data of the structure.     

Method for estimating parameters of practical ship manoeuvring models based on the combination of RANSE computations and System Identification

In this work a method for estimating parameters of practical ship manoeuvring models based on the combination of RANSE computations and System Identification procedure is investigated, considering as test case a rather slender twin screw and two rudders ship. The approach consists in the estimation of the hydrodynamic coefficients applying System Identification to a set of free running manoeuvres obtained from an in-house unsteady RANS equations solver, which substitute the usually adopted experimental tests at model or full scale. In this alternative procedure the numerical quasi-trials (in terms of kinematic parameters time histories and, if needed, forces time histories) are used as input for the System Identification procedure; the aim of this approach is to reduce external disturbances that, if not properly considered in the mathematical model, may compromise the identification results, or at least amplify the well-known “cancellation effects”. Furthermore, the CFD results provide information both in terms of flow field variables and hydrodynamic forces on the manoeuvring ship. These data may be adopted for a better understanding of the complex flow during manoeuvres, especially at stern, providing also additional information about the interaction between the various appendages (including rudders) and the hull. The identification procedure is based on an off-line genetic algorithm used for minimizing the discrepancy between the reference manoeuvres from CFD and those simulated with the system based modular model. The discrepancy was measured considering different metric functions and simplified formulations which consider only the main macroscopic parameters of the manoeuvre; the metrics have been analyzed in terms of their capability in reproducing the time histories and in limiting the cancellation effect of the hydrodynamic derivatives.     

Localization of a broadband source using a matched-mode procedurein the time-frequency domain

In this paper, we present a nonconventional matched-mode procedure for localizing a broadband source in the time-frequency domain. This hybrid coherent and incoherent approach exploits both the temporal and spatial characteristics of the multimode arrival structure at a receiving sensor array. In the previous work, a time-domain technique was developed to deal with narrowband signals coherently. It consists of the following three steps. The first step employs a receiving sensor array to separate the modes by the conventional modal filtering approach. The second step is to estimate the energy and relative arrival times of the various modes which arrive at the receiver. The last step uses the differences of modal travel times to estimate the source range, and uses the ratios of modal energies to estimate the source depth. Here, we employ bandpass filters to divide the received broadband signal into several subfrequency bands, and apply the first and second steps of the previously developed coherent narrowband technique to the subfrequency bands in the time domain. The results obtained from subfrequency bands are then combined incoherently in the frequency domain to produce an estimate of the source position. Numerical simulation of an experiment with explosive sources at the shallow water site of the Yellow Sea is presented     

New Method for Predicting the Motion Responses of A Flexible Joint Multi-Body Floating System to Irregular Waves

A new hybrid method of frequency domain and time domain is developed in this paper to predict the motion responses of a flexibly joint multi-body floating system to irregular waves. The main idea of the method is that the three-dimensional frequency method is used to obtain the hydrodynamic coefficients and the response equations are solved in time domain step by step. All the forces can be obtained at the same time. The motions and nonlinear mooring forces of a box type six-body floating system are predicted. A comparison of the theoretical method-based solutions with experimental results has shown good agreement.     

Two phase modal analysis of nonlinear sloshing in a rectangular container

Sloshing, or liquid free surface oscillation, in containers has many important applications in a variety of engineering fields. The modal method can be used to solve linear sloshing problems and is the most efficient reduced order method that has been used during the previous decade. In the present article, the modal method is used to solve a nonlinear sloshing problem. The method is based on a potential flow solution that implements a two-phase analysis on sloshing in a rectangular container. According to this method, the solution to the mass conservation equation, with a nonpenetration condition at the tank walls, results in velocity potential expansion; this is similar to the mode shapes used in modal method. The kinematic and dynamic boundary conditions create a set of two-space-dimensional differential equations with respect to time. The numerical solution of this set of differential equations, in the time domain, predicts the time response of interfacial oscillations. Modal method solutions for the time response of container sloshing due to lateral harmonic oscillations show a good agreement with experimental and numerical results reported in the literature.     

Model Simplification and Parameter Modification for Offshore Platforms Using Experiment Data

- In this paper, a computation method has been developed so as to compare the finite element method (FEM) with the test results directly. The structure is divided into the "master" and "slave" degrees of freedom. The simplified model can be obtained with modal reduction. Then the design sensitivity analysis of the eigenvalues and eigenvectors has been carried out using the modal frequency and modal shape of the test. A two-story frame structure and a jacket model structure have been calculated. Meanwhile, the modified coefficient, the FEM computational and experimental values have been given. It has been shown that the FEM model modified using the test modal value is efficient.     

Wave grouping characteristics in nearshore Great Lakes

The recently advanced approach of wavelet transform is applied to the analysis of wave data measured in the nearshore areas of the Great Lakes. The conventional spectrum analysis of wave time series in the frequency domain can be readily generalized to the frequency and time domain using the wavelet transform. The traditional Fourier transform approach has not been able to directly assess the time localized nature of wave groups. With the application of wavelet transformation, the relatively unexplored wave grouping characteristics come to light as the predominant feature of wave processes.     

X波段雷达海面回波谱宽特征研究

海面电磁回波频谱宽度与海浪波高密切相关,可应用频谱宽度进行海浪有效波高反演。本文应用线性滤波法仿真出了海表散射面元在雷达视向上的投影速度,建立了回波谱宽模型,分析了雷达空间分辨率、回波时间序列长度及海洋环境参数等因素对频谱宽度的影响,同时还针对如何在实际观测过程中选择回波时间序列长度、观测方位角等参数进行了讨论。最后还将理论结果与CSIR-X波段雷达实测数据谱宽估计结果进行了比较。结果表明,剔除雷达噪声以及频率泄露的影响后,基于高斯分布标准偏差的谱宽估计方法所得结果与理论结果吻合很好,这从而证明了理论结果的可靠性。本文所得结果对海浪有效波高反演具有一定参考价值。     

Using incomplete modal data for damage detection in offshore jacket structures

The development of robust techniques for early damage detection for offshore structures is crucial to avoid the possible catastrophe caused by structural failures. This article applies the cross-model cross-mode (CMCM) method for damage detection that is capable of identifying the damage to individual members of offshore jacket platforms, when limited, spatially incomplete modal data is available. Basically, the CMCM method is classified as a direct, physical property adjustment model updating method. Implementing this method requires only a few modes measured from the damaged structure. In dealing with spatial incompleteness, this paper investigates both model reduction and modal expansion techniques. Specifically, either Guyan (static condensation) or SEREP (system equivalent reduction expansion process) transformation matrix, between the master and slave degrees-of-freedom, is employed in the model reduction or modal expansion process. One theoretical development is an iterative procedure to compute the transformation matrix associated with the (unknown) damaged structure. Numerical studies are conducted for a jacket platform with multiple damaged members based on synthetic data generated from finite-element models. The results suggest that (i) Guyan scheme always outperforms SEREP, (ii) model reduction is always better than modal expansion, and (iii) the CMCM method in conjunction with iterative Guyan reduction approach yields the best damage location and severity estimate.     

Time Domain Simulation of Transient Responses of Very Large Floating Structures Under Unsteady External Loads

A time domain finite element method （FEM） for the analysis of transient elastic response of a very large floating structure （VLFS） subjected to arbitrary time-dependent external loads is presented. This method is developed directly in time domain and the hydrodynamic problem is formulated based on linear, inviscid and slightly compressible fluid theory and the structural response is analyzed on the thin plate assumption. The time domain finite element procedure herein is validated by comparing numerical results with available experimental data. Finally, the transient elastic response of a pontoon-type VLFS under the landing of an airplane is computed by the proposed time domain FEM. The time histories of the applied force and the position and velocity of an airplane during landing are modeled with data from a Boeing 747-400 jumbo jet.     

基于VMD对SAR海洋内波参数的自动反演

经验模态分解(empirical mode decomposition,简称EMD)是反演海洋内波参数的有效方法之一,但由于EMD存在模态混叠等问题,对海洋内波进行参数反演时会产生一定误差。相比较于EMD,变分模态分解(variational mode decomposition,简称VMD)能够有效地抑制模态混叠现象。为了更好地对海洋内波进行参数反演,提出了一种基于VMD对合成孔径雷达(synthetic aperture radar,简称SAR)遥感图像中的内波参数进行自动反演的方法。该方法先对SAR图像进行Canny处理,获取图像中的内波条纹信息,再根据内波传播方向自动选取灰度剖面;然后利用集合经验模态分解(ensemble empirical mode decomposition,简称EEMD)信号特征自适应分解模态函数的特点,再将分解得到的有效模态数作为VMD中参数K的参考值;最后利用VMD分解后的数据进行内波参数反演。试验结果表明:通过对Canny预处理后的条纹信息进行灰度剖面自动选取,解决了人为选取剖面所可能导致的误差;通过对剖面信号进行VMD处理不仅解决了EMD模态混叠的问题,成功地反演出内波的前导波振幅,而且所反演的结果与EEMD反演参数以及实测资料数据吻合得很好。     

基于谱分析的GPS浮标测量参数提取方法

介绍采用谱分析原理提取海浪信息的方法。根据GPS浮标测量得到的海浪高度数据,先通过滑动平均的方法分离海浪信息和潮流潮位信息,然后对海浪信息采用AR模型法进行功率谱估计,并对海浪信息进行连续小波变换,对AR模型法功率谱估计结果和小波变换结果进行综合研究,提取海浪特征参数———周期。最后通过一个GPS浮标试验对上述方法进行了验证。