波浪作用下海洋石油井架模态参数识别与承载力评价

提出基于试验模态参数识别和优化算法修改动力有限元模型的海洋石油井架安全承载能力评价方法.推导基于波浪脉动的试验模态参数识别公式,介绍有限元模型动力修正的一阶优化方法.以海洋随机波浪脉动作为激励,对勘探三号海洋石油井架进行了现场模态试验,识别出前二阶固有频率.以此作为有限元模型修正的基础,应用一阶搜索优化算法,对动力有限元模型进行了修正,结果表明:该修正模型比较准确反映了井架目前的结构状态,能够用于进一步的静、动力分析和安全承载能力评价.     

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

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

Numerical Simulation on the Resistance Performance of Ice-Going Container Ship Under Brash Ice Conditions

Ice resistance prediction is a critical issue in the preliminary design of ships navigating brash ice conditions, which is closely related to the safety of a ship to navigate encounter brash ice, and has significant effects on the kinds of propellers and motor power needed. In research on this topic, model tests and full-scale tests on ships have thus far been the primary approaches. In recent years, the application of the finite element method (FEM) has also attracted interest. Some researchers have conducted numerical simulations on ship–ice interactions using the fluid–structure interaction (FSI) method. This study used this method to predict and analyze the resistance of an ice-going ship, and compared the results with those of model ship tests conducted in a towing tank with synthetic ice to discuss the feasibility of the FEM. A numerical simulation and experimental methods were used to predict the brash ice resistance of an ice-going container ship model in a condition with three concentrations of brash ice (60%, 80%, and 90%). A comparison of the results yielded satisfactory agreement between the numerical simulation and the experiments in terms of both observed phenomena and resistance values, indicating that the proposed numerical simulation has significant potential for use in related studies in the future.     

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.     

Parametric study for underwater blast-induced pipeline response embedded in marine sediments

Abstract

Blast response of submerged pipelines has been a research focus in recent years. In this article, a three-dimensional numerical model is established to investigate dynamic response of pipelines due to underwater explosion. The u–p approximation is integrated into finite element method (FEM) to simulate pore water effect in the seabed. Numerical continuity between hydraulic pressure in the flow field and pore pressure in the marine sediment is guaranteed to realize the blast response of submerged pipelines in ocean environment. Both fluid–structure interaction (FSI) and pipeline–seabed interaction (PSI) have been considered in the proposed model simultaneously. A comprehensive parametric study is carried out after validation of the present model with test data from underground explosion and underwater explosion, respectively. The effect of embedment depth, TNT equivalent, stand-off distance, pipeline diameter, and pipeline thickness to blast response of the submerged pipelines is investigated based on numerical results. Variation of deformation patterns and stress distribution of the pipeline with various installation and structure parameters has been illustrated and discussed to facilitate engineering practice.     

Numerical modeling on the interaction of internal solitary wave with slope-shelf and modal analysis

On the basis of a nonhydrostatic numerical model,the interaction of internal solitary wave with slope-shelf was studied.The breaking and polarity transformation were analyzed.A "kink" structure,due to shoaling topography and higher nonlinear effect,was found to be generated by the leading wave before breaking.Coherent vortex shedding behind the leading wave was presented.The evolution characteristics of the modal structure were analyzed based on the empirical orthogonal function method.The modal structure was complicated due to the effect of the variable topography,especially when breaking occurred.In the performed experiments,the contributions to the total variance from higher mode jumped from no more than 20% to over 40%.     

一种基于TIN模型的LOD简化算法

为了实现不规则三角网（TIN）地形模型的动态模型的动态细节简化模型,本文提出了一种有效的方法,即在TIN的凸壳构网算法和累进格网算法的基础上,设计一种混合型数据结构,提出边折叠的限制条件,解决了边折叠时易产生的边自交现象,实现了TIN模型的LOD快速简化。     

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.     

Development of structural modeling techniques for evaluating HDPE plastic net pens used in marine aquaculture

Finite-element modeling (FEM) techniques are developed to determine the structural capabilities of net pen flotation structures made of high-density polyethylene (HDPE). The modeling approach uses shell elements and localized failure criteria to predict critical loading conditions. Finite element modeling simulations were performed using values for the modulus of elasticity for weathered HDPE determined from a series of tensile tests. Poisson's ratio and yield stress used with the approach were obtained from standard manufacturers values. To investigate the method, a series of experiments were performed in the laboratory by testing circular sections of HDPE pipe to localized failure (“kinking”). The same test was replicated with the FEM using localized failure techniques. Laboratory and model results were within 16%. A FEM was then built to represent the complex geometry of a net pen flotation structure deployed at an operational fish farm located in Eastport, ME, USA. Simulations were performed using attachment line tension values measured at the site. The goal was to assess flotation pipe stress levels for typical operational conditions. Simulations were also performed to investigate the maximum capabilities of the net pen structure with different attachment line configurations. Quantifying the operational limits will become more important as these systems are considered for more exposed, energetic environments.     

Hydroelastic impact of a wedge-shaped body

In this paper a hydroelastic method for the design of marine structures is used to study the impact of a wedge-shaped body onto a calm free-surface. The method combines computational fluid dynamics and dynamic finite-element techniques to predict the hydroelastic response of a structure in the time domain. The basic idea of the method is to predict the stress field on the fluid-structure interface due to rigid-body impact, and then to transfer the information to a wet modal model. The one-way transfer of information from the fluid simulation to the structure has two important advantages with respect to a fully coupled simulation: a single rigid-body motion computation can be reused for subsequent analyses after the structure has been redesigned, and cost of the fluid dynamic simulation is not increased due to iteration which is required in coupled simulations. The maximum displacement of an elastic wedge during the impact and penetration stages are compared with published results from a fully coupled theoretical model and a fully coupled numerical solver. The present method requires approximation of the flexural added mass, and the comparison with other results shows that the approximation is satisfactory for predicting maximum displacement.     

基于固有频率的海洋平台损伤检测方法的改进

针对海洋平台损伤检测方法中需要未损伤模型参数的不便,提出了一种相对频率变化量的方法。对不同类型构件的损伤及同一类型构件在不同位置损伤的情况,进行结构模态分析,研究平台结构的固有频率对损伤、损伤位置和损伤程度的敏感性,以及平台质量变化对频率的影响。利用两次测量的相对频率变化量,可以发现损伤的存在,并进行损伤的初步定位。     

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.     

Numerical simulation of a short flexible pipe subject to forced motion and vortex-induced vibration

A series of numerical simulations about a small scale(aspect ratio:63.2) flexible pipe undergoing forced harmonious oscillation and vortex-induced vibration(VIV) have been taken into account.The wake hydrodynamics and pipe deformation were accomplished by ANSYS MFX solution strategy designed for fluid-structure interaction(FSI) problem with well-performed LES model.The configuration of structured mesh,multi-domain design,different mesh stiffness admeasured by User Fortran ensured that the numerical task was competent to deal with large deformation related to this case.The introduction of instantaneous amplitude definition and modeless component decomposition method(Chen and Kim,2008) was helpful to reveal much more information from modal analysis.Most results from numerical simulation are generally consistent with those from model test(Choi and Hong,2000) via the comparison between them.As supplementary to model test,visualization of the vortex wake was also provided.It has been proved that the forced oscillation doesn't only excite a complicated dumbbell-like wake pattern around the outer thimble,but also results in inner flow inside the PVC pipe.The velocity of the inner flow increases with the frequency of forced oscillation.     

修正型缓坡方程的有限元模型

与缓坡方程相比,修正型缓坡方程增加了地形曲率项和坡度平方项,从而提高了数值求解的复杂性。本文将计算域划分为内域和外域,内域为水深变化区域,使用修正型缓坡方程,其中的地形曲率项和坡度平方项可用有限单元各节点的水深信息和单元插值函数表示,外域为水深恒定区,速度势满足Helmholtz方程,通过内外域的边界匹配建立有限元方程,并用高斯消去法求解。进而分别模拟了波浪传过Homma岛和圆形浅滩的变形,其结果与相关的解析解和实验数据吻合良好,证明了本文有限元模型的正确性。同时,通过与实验数据的对比也明显看出,在地形坡度较陡的情况下,修正型缓坡方程较缓坡方程具有更高的计算精度。     

On natural frequencies and modal shapes in two-dimensional asymmetric and symmetric moonpools in finite water depth

In this paper theoretical models are proposed for computing the natural frequencies and modal shapes of two-dimensional asymmetric and symmetric moonpools in the finite water depth. The boundary value problem is solved by using a domain decomposition approach. On the outer vertical boundary bounded by the beam of the two bodies, linearized velocity potential is assumed to be nil. Eigenvalue problem is formulated by matching the velocity potential and fluid flux on the common boundaries to obtain the natural frequencies and modal shapes of the free surface elevation. In the symmetric moonpool cases, so-called single mode approximations (SMA) have been derived and can be adopted for rapid estimation of the natural frequencies for both piston and sloshing modes. The present results have been extensively compared with the solutions using the two-dimensional infinite water depth model developed by Molin [1], the numerical solutions and experimental data by Faltinsen et al. [2]. It is found that the solutions have been improved from the infinite water depth model. It is demonstrated that the proposed models can well predict the resonance frequencies and modal shapes for the two-dimensional asymmetric and symmetric moonpools.     

基于实测数据的海上风电结构动力响应分析

地震是危害海上风电结构作业安全的重要环境因素,目前,国内尚未公开发表真实地震响应下,海上风电结构的实测动力响应数据。分析了某地震活动区海上风电结构的实测地震响应,采用随机子空间识别方法进行风机的模态识别,阐述了风机机舱偏航将引起前后、左右两个正交方向振动的耦合,并从理论上证明了利用耦合、解耦数据识别模态参数的差异。结果表明：1）耦合与解耦信号识别的频率、阻尼比完全相同,而耦合信号识别的模态振型与偏航角有关;2）地震作用会对结构产生巨大冲击;3）非地震作用下,风机塔筒前后、左右第一阶弯曲模态为主要模态,地震作用可以激发风机的高阶模态,使得塔筒中上部而不是顶部的振动响应最大。此分析对地震活动区海上风电结构的抗震设计具有一定的参考价值。     

A coupled-mode model for water wave scattering by horizontal, non-homogeneous current in general bottom topography

A coupled-mode model is developed for treating the wave–current–seabed interaction problem, with application to wave scattering by non-homogeneous, steady current over general bottom topography. The vertical distribution of the scattered wave potential is represented by a series of local vertical modes containing the propagating mode and all evanescent modes, plus additional terms accounting for the satisfaction of the free-surface and bottom boundary conditions. Using the above representation, in conjunction with unconstrained variational principle, an improved coupled system of differential equations on the horizontal plane, with respect to the modal amplitudes, is derived. In the case of small-amplitude waves, a linearised version of the above coupled-mode system is obtained, generalizing previous results by Athanassoulis and Belibassakis [J Fluid Mech 1999;389:275–301] for the propagation of small-amplitude water waves over variable bathymetry regions. Keeping only the propagating mode in the vertical expansion of the wave potential, the present system reduces to an one-equation model, that is shown to be compatible with mild-slope model concerning wave–current interaction over slowly varying topography, and in the case of no current it exactly reduces to the modified mild-slope equation. The present coupled-mode system is discretized on the horizontal plane by using second-order finite differences and numerically solved by iterations. Results are presented for various representative test cases demonstrating the usefulness of the model, as well as the importance of the first evanescent modes and the additional sloping-bottom mode when the bottom slope is not negligible. The analytical structure of the present model facilitates its extension to fully non-linear waves, and to wave scattering by currents with more general structure.     

Numerical Analysis of Factors Influencing Implosion

-By using gas-liquid two-phase flow theory, a modified mathematical model based on the computational fluid dynamics method SIMPLE (Semi-Implicit Method for Pressure-Linked Equations) is introduced to investigate implosion phenomena in high pressure chambers systematically. A theoretical simulation-prediction method, which is independent of experimental data, is developed in the paper and great improvement has been made on the topic. In the paper, various implosion situations have been simulated and predicted. Effects of a series of factors influencing implosion results and methods of reducing implosion danger have been analysed. The analysis results are of importance to underwater engineering practice.     

SIMULATION OF WIND FIELD AND THE NUMER-ICAL CALCULATION OF STORM SURGE OF THE BOHAI SEA

Up to now, available method of numerical forecasting and suitable wind field model for the Bohai Sea storm surge have been few. In this paper, through the analysis of the weather situation is presented a mathematical model for the wind fields involved mainly the deformation field of a high pressure matched with a low pressure, the temperate cyclone, the cold wave and the northword typhoon. Meanwhile, numerical computations combined with the nonlinear storm surge models are made by using "ADI" method. The computed results are generally coincident with the practical observations. It has showed a success in the simulated wind field and the feasibility of using "ADI" method to forecast the Bohai Sea storm surge.     

海洋平台基于隔振锥体装置的减振效应研究

针对导管架式海洋平台结构型式与特点,利用冲击隔振理论,将冰锥体与隔振技术结合起来,提出一种能有效降低冰荷载激励的"两级柔性隔振锥体"设计方案。建立了两极隔振锥体简化计算模型,确定了系统振动传递系数参数之间的关系。通过建立海洋平台有限元模型,采用有限元程序,进行了模拟冰荷载激励下的振动响应分析,并对安装隔振装置的小比例模型结构进行了冲击荷载下的减振试验。结果表明采用"两级柔性隔振锥体"能有效降低传递给海洋平台的冰力极值,减小海洋平台冰激励的响应。