Time-domain,shallow-water hydroelastic analysis of VLFS elastically connected to the seabed

In order to ensure the safe operation of a VLFS, a combination of mooring, breakwater and other motion reducing systems is employed. In the present work, the transient hydroelastic response of a floating, thin elastic plate, elastically connected to the seabed, is examined. The plate is modelled as an Euler-Bernoulli strip, while the linearized shallow water equations are used for the hydrodynamic modelling. Elastic connectors are approximated by a series of simple spring-dashpot systems positioned along the strip. A higher order finite element scheme is employed for the calculation of the hydroelastic response of the strip-connector configuration over the shallow bathymetry. After the definition of the initial-boundary value problem, its variational form is derived and discussed. Next, on the basis of the aforementioned formulation, an energy balance expression is obtained. The effect of variable bathymetry on the response of a two connector-strip system is examined by means of three seabed profiles featuring a flat bottom, an upslope and a downslope environment. For the flat bottom case, the strip response mitigating effect exerted by the employment of two and three elastic connectors is considered. Finally, by means of the derived energy balance equation, the energy exchange is monitored, providing a valuable insight into the transient phenomena that take place in the studied configurations.     

Fluid-structure-material coupling analysis for a floating laminated structure consisting of high-stiffness panels and a soft core

This paper presents a fluid-structure-material coupling analysis for the interaction between water waves and a very large floating laminated structure (VLFLS), which is consisted of two enhanced ultrahigh-performance concrete (UHPC) panels and a middle lightweight foamed rubber core. The representative volume element (RVE) method is used to design the mechanical properties of enhanced UHPC and foamed rubber, and the parameterized formulas are presented to reveal the dependency between macroscale mechanical properties and mesoscale hierarchical characteristics. By idealizing the rubber core as a uniformly distributed spring layer, an eighth-order differential equation of motion of the laminated structure is derived. In the context of linear potential flow theory, a hydroelastic analytical model is developed for the floating laminated structure with finite length under wave action. In the process of solving velocity potentials, a complicated dispersion equation for the wave motion below the laminated structure is derived, and this equation contains two pairs of conjugate complex roots with positive real parts. The various hydrodynamic quantities, including reflection coefficient, transmission coefficient, deflection, shear force, and bending moment, are calculated. The hydroelastic model is confirmed by considering the convergence of calculation results and the energy conservation of wave propagation. The coupled effects of wave action, material characteristics, structural parameters, and edge conditions on the hydroelastic and mechanical response of the floating laminated structure are clarified to provide important information regarding the optimal design of such structures.     

Hydroelastic analysis of oblique irregular waves with a pontoon-type VLFS edged with dual inclined perforated plates

The fluid-structure interaction of oblique irregular waves with a pontoon-type very large floating structure (VLFS) edged with dual horizontal/inclined perforated plates has been investigated in the context of the direct time domain modal expansion theory. For the hydroelastic analysis, the boundary element method (BEM) based on time domain Kelvin sources is implemented to establish water wave model including the viscous effect of the perforated plates through the Darcy’s law, and the finite element method (FEM) is adopted for solving the deflections of the VLFS modeled as an equivalent Mindlin thick plate. In order to enhance the computing efficiency, the interpolation-tabulation scheme is applied to assess rapidly and accurately the free-surface Green function and its partial derivatives in finite water depth, and the boundary integral equation of a half or quarter VLFS model is further established taking advantage of symmetry of flow field and structure. Also, the numerical solutions are validated against a series of experimental tests. In the comparison, the empirical relationship between the actual porosity and porous parameter is successfully applied. Numerical solutions and model tests are executed to determine the hydroelastic response characteristics of VLFS with an attached anti-motion device. This study examines the effects of porosity, submerged depth, inclined angle and gap distance of such dual perforated anti-motion plates on the hydroelastic response to provide information regarding the optimal design. The effects of oblique wave angle on the performance of anti-motion and hydroelastic behavior of VLFS are also emphatically examined.     

超大型浮体结构考虑非线性项的波浪力响应计算

为求解超大型浮体结构波浪响应问题,将其简化成弹性梁板模型。以有限元软件为基础来建立超大型浮体结构的弹性梁板模型,在势流理论的基础上,将速度势分解为入射势、绕射势和辐射势三部分,并利用格林函数法推导出一阶、二阶波浪速度势的表达式。以此为理论基础,编制了相应的计算程序来计算浮体结构所受到的二阶波浪力,最后建立了梁板模型的水弹性响应方程,对非线性波浪力作用下的浮体结构的水弹性响应进行求解分析,并分析了波向对响应振幅的影响。     

Time domain springing analysis on a floating barge under oblique wave

This paper studies the springing response of a flexible floating barge under an oblique wave. A time domain Rankine panel method was used to represent fluid motion surrounding a flexible seagoing vessel, and a finite element method was used for structural response. For accurate prediction of the structural response under an oblique wave, special attention was given to the structural model, such as the effect of warping distortion and bending-torsion coupling. The Vlasov assumption was followed for a deformable beam element to take into account the effect of warping distortion so that the cross section of the beam deforms out of its original plane without changing its cross-section contour. The coupled equation for both the fluid and structural domain was solved by using the implicit iterative method. A fixed point iteration with a relaxation scheme was employed in this study with the aid of the Aitken δ² process seeking acceleration of solution convergence. Accuracy of a developed computer program was verified through the comparison with experimental work done by Remy et al. (Experimental and numerical study of the wave response of a flexible barge, Hydroelasticity in Marine Technology, Wuxi, 2006) resulting in good correspondence between the two results.     

Study on sloshing in cargo tanks including hydroelastic effects

The sloshing problem in cargo tanks is studied through experiments and numerical analysis. The fluid motion is described using a higher-order boundary element method and the structural response by a thin plate theory. It was found that hydraulic jumps are formed when the excitation frequency is close to the resonance frequency in the case of low filling depth. In the case of high filling depth, the flow resonates and hits the top of the tank, thus inducing a large impact pressure. The pressure on the flexible plate shows amplified initial peaks followed by oscillatory components, the frequency of which coincides with the natural frequencies of the plate in water as a result of hydroelastic effects. Received for publication on Nov. 18, 1998; accepted on May 14, 1999     

Hydroelastic analysis of a nonlinearly connected floating bridge subjected to moving loads

In this paper, the motion equations for the nonlinearly connected floating bridge, considering the nonlinear properties of connectors and vehicles’ inertia effects, are proposed. The super-element method is used to condense the whole calculation scale, and the direct integration and Newton–Raphson iteration method are applied to solve the reduced equations. Based on the modal and static analyses, the dynamic displacement and connection forces characteristics of a floating bridge with nonlinear connectors subjected to moving loads are investigated. It is found that nonlinearity and initial gap of the connectors are important for the hydroelastic response of a nonlinearly connected floating bridge.     

基于ALE算法的无转角和有转角铝制加筋板楔形体水弹性砰击的数值模拟研究

文章采用基于任意拉格朗日—欧拉(ALE)算法的显式有限元技术研究水弹性砰击现象,针对已开展的铝制加筋板楔形体结构入水砰击模型实验,开展了数值模拟比较工作。该楔形体底部斜升角为20度,底部两侧是包含三根纵骨和两根横梁的加筋板结构,两侧结构刚度不同。预报了模型无转角和有转角典型工况的砰击入水过程,得到的入水加速度、底部加筋板结构纵骨应力和横梁响应与模型实验结果吻合较好。研究表明该ALE算法具备模拟船舶局部结构的水弹性砰击流固耦合问题的能力。     

三维水弹性力学中的结构广义质量计算

该文给出了三维水弹性力学中弹性模态位移振型归一化时广义质量的计算方法,并应用三维有限元方法计算了箱式大型浮体结构的固有频率和广质量,经与均匀梁解析解的比较表明,两种方法差异不大,因此,采用均匀染解析解来计算这种结构的固有频率和广义质量是完全可行的,简化了计算工作。     

支撑条件对浮动式厚弹性板水弹性行为的影响分析（英文）

The hydroelastic response of very large floating structures(VLFS) under the action of ocean waves is analysed considering the small amplitude wave theory. The very large floating structure is modelled as a floating thick elastic plate based on TimoshenkoMindlin plate theory, and the analysis for the hydroelastic response is performed considering different edge boundary conditions.The numerical study is performed to analyse the wave reflection and transmission characteristics of the floating plate under the influence of different support conditions using eigenfunction expansion method along with the orthogonal mode-coupling relation in the case of finite water depth. Further, the analysis is extended for shallow water depth, and the continuity of energy and mass flux is applied along the edges of the plate to obtain the solution for the problem. The hydroelastic behaviour in terms of reflection and transmission coefficients, plate deflection, strain, bending moment and shear force of the floating thick elastic plate with support conditions is analysed and compared for finite and shallow water depth. The study reveals an interesting aspect in the analysis of thick floating elastic plate with support condition due to the presence of the rotary inertia and transverse shear deformation. The present study will be helpful for the design and analysis of the VLFS in the case of finite and shallow water depth.     

基于一个二维修正水弹性方法的集装箱船体梁动态前极限强度评估研究

海上极端波因其巨大的波高常常导致船体的极限破坏。该文提出了一个二维的修正水弹性方法来研究一个集装箱船船体梁在极端波中的动态前极限强度。传统的极限强度评估基于准静态方法,没有动态效应被考虑。而船体在波浪下的动态结构响应是基于水弹性方法,传统的水弹性方法并不能计算船体梁的动态非线性强度。该二维修正的水弹性方法考虑时域波浪和非线性船体梁之间的耦合,将水弹性方法和Smith方法结合,用Smith方法计算船体梁的刚度,而其刚度与船体梁的强度和变形曲率有关。所以该时域的非线性刚度被用于修改水弹性方法里的常数项的结构梁刚度。几组极端波模型被用以产生船体梁的大变形和非线性动态垂向弯矩。文中分别采用修正水弹性方法和普通水弹性方法,通过改变四个重要的极端波参数如极端波最大波高、规则波的波高、波速和波长等来研究其对船体梁船中处的大变形转角和非线性垂向弯矩的影响,通过采用修正的水弹性方法计算得来的结果与水弹性方法计算得来的结果相比较,得到了一些差异和结论。     

Efficient time-domain approach for hydroelastic-structural analysis including hydrodynamic pressure distribution on a moored SFT

This study investigates the hydroelastic analysis of a moored SFT (submerged floating tunnel) and the corresponding hydrodynamic pressure distribution under wave excitations. Time-domain discrete-module-beam (DMB) method, in which an elastic structure is modeled by multiple sub-bodies with beam elements, is employed to express the deformable tunnel with multiple mooring lines. Moreover, the top-down scheme is also adopted for detailed structure analyses with less computational cost, which applies the calculated hydrodynamic pressure distribution over SFT's surface to the three-dimensional finite element model. The hydrodynamic pressure includes both wave-induced diffraction pressure and motion-induced radiation pressure. For the validation of the developed numerical approach, comparisons are made with computationally intensive hydroelastic-structural direct-coupled method, two-dimensional wave flume experiment, and independently developed inhouse moored-SFT-simulation program. Furthermore, the influences of flexural motions with buoyancy-weight ratio (BWR) (or bending stiffness) and regular/irregular wave conditions on the dynamic pressure distribution and the resulting local stresses are investigated.     

Numerical analysis of the hydroelastic behavior of a vertical plate due to solitary waves

Interactions of a vertical elastic plate with fully nonlinear water waves were simulated. Utilizing the mixed Eulerian Lagrangian method for the free-surface flow and the finite element method for the deflection of an elastic plate, a fully coupled scheme for accurately determining fluid–plate motions was developed. Using this scheme, some modifications to the solvers for both fluid and plate were made. A hybrid wave-absorbing beach was installed to prevent wave reflection from the end of the wave tank. A fourth-order Runge–Kutta time-marching scheme with a uniform time step was applied to achieve numerical stability. The method was validated by simulating the wave generated by the initial deformation of a vertical plate and comparing the result with the corresponding analytical solution. For further validation, the hydroelastic behavior of a vertical plate induced by a pulse-type wave (where the initial pulse-type elevation of the free surface is specified) was computed, and the result was compared with another numerical result from a mode-expansion method. The interaction of a surface-piercing plate with nonzero initial free surface was then simulated, and the result was compared with the corresponding linear analytical solution. Finally, the hydroelastic response of a surface-piercing vertical plate due to a solitary wave (generated by actuating the vertical plate at the right end of the tank only at the beginning) was computed and investigated systematically.     

Seismic response of liquid-filled tank with baffles

In thispaper, the effects of a rigid baffle on the seismic response of liquid in a rigid cylindrical tank are evaluated. A baffle is an annular plate which supplies a kind of passive control on the effects of ground excitation. The contained liquid is assumed incompressible, inviscid and has irrotational motion. To estimate the seismic response, the method of superposition of modes has been applied. To analyze the rigid tank response, Laplace＆#39;s equation is considered as the governing equation of the fluid domain, in both time and frequency domains. The boundary element method （BEM） is employed to evaluate the natural modes of liquid in a cylindrical tank. To gain this goal, the fluid domain is divided into two upper and lower parts partitioned by the baffle. Linearized kinematic and dynamic boundary conditions of the free surface of the contained liquid have been considered.     

带有挡板的注液舱室地震响应分析（英文）

In thispaper, the effects of a rigid baffle on the seismic response of liquid in a rigid cylindrical tank are evaluated. A baffle is an annular plate which supplies a kind of passive control on the effects of ground excitation. The contained liquid is assumed incompressible, inviscid and has irrotational motion. To estimate the seismic response, the method of superposition of modes has been applied. To analyze the rigid tank response, Laplace's equation is considered as the governing equation of the fluid domain, in both time and frequency domains. The boundary element method(BEM) is employed to evaluate the natural modes of liquid in a cylindrical tank. To gain this goal, the fluid domain is divided into two upper and lower parts partitioned by the baffle. Linearized kinematic and dynamic boundary conditions of the free surface of the contained liquid have been considered.     

An advanced theory of thin-walled girders with application to ship vibrations

The paper presents an outline of the advanced theory of thin-walled girders. The improvement includes shear influence on torsion as an extension of shear influence on bending. The analogy between bending and torsion is recognized and pointed out throughout the paper. Complete differential equations of coupled flexural and torsional vibrations for a prismatic girder are derived. In addition, the 8 d.o.f. beam finite element, utilizing the energy approach, is constituted with stiffness and mass matrices, and load vectors. The paper describes determining of geometrical properties of multi-cell open cross-sections by employing the strip element method. Numerical procedures for vibration analyses are outlined. Furthermore, dry natural vibrations of a VLCS (Very Large Container Ship) are analysed by 1D FEM model as a prerogative for hydroelastic analyses of these relatively flexible vessels. Influence of transverse bulkheads is taken into account by increasing torsional stiffness of the ship hull proportionally to their deformation energies. Validation of 1D FEM model is checked by correlation analysis with the vibration response of the fine 3D FEM model.     

深海系缆动张力的超谐动力响应

文章考虑深水系缆的拉伸和弯曲变形,取正交坐标和切向量描述任意结点的6个广义位移,考虑缆变形的几何非线性和水动力影响,基于细长杆理论推导缆索的单元刚度矩阵,引入单元内张力一致假设,推导出12×12的非线性刚度矩阵,并与有限元软件对接。针对1800m水深平台的系缆,考虑系缆端点垂荡运动,以及端点垂荡、横荡和海流作用两种状态下系泊系统的动力响应。计算结果表明,端点运动和海流作用,引起系缆系统的二倍超谐动力响应,并且随着流速增大,系缆超谐动力响应增大,系缆的超谐动张力分量对于缆索疲劳累积损伤具有重要影响,在系缆疲劳分析时应该予以考虑。     

6750箱集装箱船非线性波激振动和颤振响应分析

本文采用三维时域非线性水弹性方法分析了一艘6750箱集装箱船的水弹性响应以及运动和垂向弯矩特征。通过考虑入射波力、静水恢复力、砰击效应的非线性,研究了在恶劣海况下船体的非线性运动和垂向弯矩响应,同时分析了波激振动及颤振对垂向弯矩的影响。数值计算结果表明:(1)非线性入射波力对运动的影响较小,但是对垂向弯矩的影响较大,使得其有明显的倍频成分,同时中垂弯矩显著大于中拱弯矩。另外,非线性入射波力也引起了明显的非线性波激振动;(2)非线性静水恢复力对运动和载荷的影响均较大,但是没有引起明显的非线性响应。非线性计算的垂荡响应小于线性结果,而纵摇和垂向弯矩响应大于线性结果;(3)砰击效应对运动的影响较小,但对垂向弯矩的影响较大,砰击效应引起了显著的船体弹性高频振动,增大了载荷幅值,但是其引起的合成中垂和中拱幅值相差不大;(4)非线性水动力的作用主要引起垂向弯矩的倍频响应,包括倍频可能引起的二节点垂向弯矩弹性共振,而砰击效应主要引起船体二节点垂向弯矩共振;(5)本文的非线性水弹性响应计算结果与Kim给出的数值计算结果吻合很好。     

Dynamic adaptive finite element analysis of acoustic wave propagation due to underwater explosion for fluid-structure interaction problems

In this paper, an investigation into the propagation of far field explosion waves in water and their effects on nearby structures are carried out. For the far field structure, the motion of the fluid surrounding the structure may be assumed small, allowing linearization of the governing fluid equations. A complete analysis of the problem must involve simultaneous solution of the dynamic response of the structure and the propagation of explosion wave in the surrounding fluid. In this study, a dynamic adaptive finite element procedure is proposed. Its application to the solution of a 2D fluid-structure interaction is investigated in the time domain. The research includes: a) calculation of the far-field scatter wave due to underwater explosion including solution of the time-depended acoustic wave equation, b) fluid-structure interaction analysis using coupled Euler-Lagrangian approach, and c) adaptive finite element procedures employing error estimates, and re-meshing. The temporal mesh adaptation is achieved by local regeneration of the grid using a time-dependent error indicator based on curvature of pressure function. As a result, the overall response is better predicted by a moving mesh than an equivalent uniform mesh. In addition, the cost of computation for large problems is reduced while the accuracy is improved.     

用三维非定常欧拉有限元方法进行波体干扰计算

本文利用理想流体的欧拉方程和用有限元方法建立了波体干扰计算模型.在自由表面上应用了线性自由面边界条件,非定常的欧拉方程用分步法求解,并用时间步进积分法分别处理压力泊松方程,速度方程和自由面上的波浪运动.为了证明本文提出的计算方法,在自由面上摇荡的半球用于作为数值计算例,把计算结果与用不同的理论方法得到的结果进行了比较.