基于CFD-FEM的浮式结构水弹性响应数值模拟

开发并验证了一种基于CFD-FEM耦合的弹性浮体水弹性响应计算模拟方法。采用CFD方法建立黏性数值水池模拟非线性波浪，弹性浮板进行有限元离散，并在交界面进行数据交互实现耦合计算；通过与水池试验数据和三维板理论在各种波浪环境下的浮体垂向位移结果对比，证实CFD-FEM耦合方法的有效性。并进一步研究了浮板的厚度、入射波波幅和浮板的三维效应对浮板水弹性响应的影响。结论表明，波幅的增加会加剧弹性浮板的水弹性响应，浮板各点处的垂向位移随波幅的增加而增大；当浮板厚度改变时，不同厚度浮板自由端处的垂向位移差异较小，而在中部等位置处，厚度对浮板的水弹性响应有较大的影响。     

Hydroelastic response of a circular plate in waves on a two-layer fluid of finite depth

The hydroelastic response of a circular, very large floating structure （VLFS）, idealized as a floating circular elastic thin plate, is investigated for the case of time-harmonic incident waves of the surface and interfacial wave modes, of a given wave frequency, on a two-layer fluid of finite and constant depth. In linear potential-flow theory, with the aid of angular eigenfunction expansions, the diffraction potentials can be expressed by the Bessel functions. A system of simultaneous equations is derived by matching the velocity and the pressure between the open-water and the plate-covered regions, while incorporating the edge conditions of the plate. Then the complex nested series are simplified by utilizing the orthogonality of the vertical eigenfunctions in the open-water region. Numerical computations are presentedto investigate the effects of different physical quantities, such as the thickness of the plate, Young＇s modulus, the ratios ofthe densities and of the layer depths, on the dispersion relations of the flexural-gravity waves for the two-layer fluid.Rapid convergence of the method is observed, but is slower at higher wave frequency. At high frequency, it is found that there is some energy transferred from the interfacial mode to the surface mode.     

Hydroelastic interaction between obliquely incident waves and a semi-infinite elastic plate on a two-layer fluid

The hydroelastic response of a semi-infinite thin elastic plate floating on a two-layer fluid of finite depth due to obliquely incident waves is investigated. The upper and lower fluids with different densities separated by a sharp and stable interface are assumed to be inviscid and incompressible and the motion to be irrotational. Simply time-harmonic incident waves of the surface and interfacial wave modes with a given angular frequency are considered within the framework of linear potential flow theory. With the aid of the methods of matched eigenfunction expansion and the inner product of the two-layer fluid, a closed system of simultaneous linear equations is derived for the reflection and transmission coefficients of the series solutions. Based on the dispersion relations for the gravity waves and the flexural–gravity waves in a two-layer fluid and Snell’s law for refraction, we obtain a critical angle for the incident waves of the surface wave mode and three critical angles for the incident waves of the interfacial wave mode, which are related to the existence of the propagating waves. Graphical representations of the series solutions show the interaction between the water waves and the plate. The effects of several physical parameters, including the density and depth ratios of the fluid and the thickness of the plate, on the wave scattering and the hydroelastic response of the plate are studied. It is found that the variation of the thickness of the plate may change the wave numbers and the critical angles. The density ratio is the main factor to influence the wave numbers of the interfacial wave modes. Finally, the stress state is considered.     

Surface gravity wave interaction with elastic bottom

In the present paper, a hydroelastic model is developed to deal with surface gravity wave interaction with an elastic bed based on the small amplitude water wave theory and plate deflection in finite water depth. The elastic bottom bed is modelled as a thin elastic plate and is based on the Euler-Bernoulli beam equation. The wave characteristics in the presence of the elastic bed is analyzed in both the cases of deep and shallow water waves. Further, the linearized long wave equation is generalized to include bottom flexibility. A generalized expansion formula for the velocity potential is derived to deal with the boundary value problems associated with surface gravity waves having an elastic bed. The utility of the expansion formula is illustrated by demonstrating specific physical problems which will play significant role in the analysis of wave structure interaction problems. Behavior of the wave spectra are discussed in the case of closed basin having a free surface and an elastic bottom topography.     

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.     

Hydroelastic Response of VLFS with An Attached Submerged Horizontal Solid/Porous Plate Under Wave Action

The hydroelastic responses of a submerged horizontal solid/porous plate attached at the front of a very large rectangular floating structure(VLFS) under wave action has been investigated in the context of linear water wave theory. Darcy's law is adopted to represent energy dissipation in pores. It is assumed that the porous plates are made of material with very fine pores so that the normal velocity across the perforated porous is linearly associated with the pressure drop. In the analytic method, the eigenfunction expansion-matching method(EEMM) for multiple domains is applied to solve the hydrodynamic problem and the elastic equation of motion is solved by the modal expansion method. The performance of the proposed submerged horizontal solid/porous plate can be significantly enhanced by selecting optimal design parameters, such as plate length, horizontal position, submerged depth and porosity. It is concluded that good damping effect can be achieved through installation of solid and porous plate.Porous plate has better damping effect at low frequencies, while solid plate has better damping effect at high frequencies. The optimal ratio of plate length to water depth is 0.25-0.375, and the optimal ratio of submerged depth to water depth is 0.09-0.181.     

Hydroelastic analysis of flexible floating interconnected structures

Three-dimensional hydroelasticity theory is used to predict the hydroelastic response of flexible floating interconnected structures. The theory is extended to take into account hinge rigid modes, which are calculated from a numerical analysis of the structure based on the finite element method. The modules and connectors are all considered to be flexible, with variable translational and rotational connector stiffness. As a special case, the response of a two-module interconnected structure with very high connector stiffness is found to compare well to experimental results for an otherwise equivalent continuous structure. This model is used to study the general characteristics of hydroelastic response in flexible floating interconnected structures, including their displacement and bending moments under various conditions. The effects of connector and module stiffness on the hydroelastic response are also studied, to provide information regarding the optimal design of such structures.     

Hydroelastic analysis of the floating plate optimized for maximum radiation damping

In previous work, the problem of optimizing the shape of a thin floating plate to maximize radiation damping was investigated. The plate was modelled with zero draft and floated on the surface of an irrotational, incompressible ocean of infinite extent. For simplicity, only rigid heave motions were considered and the damping coefficient at one wave number was maximized. In the present work, the hydroelastic properties of the optimized plate are determined and compared with those of circular and square plates. The added mass, damping, and diffraction force coefficients in each mode are determined as a function of wave number. The amplitude responses of the plate deflection and bending moments are also presented. The finite element method is used to determine the vibration mode shapes and the flow problem is analysed using the Chen and Mei variational principle wherein the potential field inside a hemisphere surrounding the plate is represented using a spherical harmonic expansion and matched on the hemisphere to an outer field described by distributing sources on the hemisphere.     

Radiation and diffraction of linear water waves by an infinitely long submerged rectangular structure parallel to a vertical wall

The radiation and the diffraction of linear water waves by an infinitely long floating rectangular structure submerged in water of finite depth with leeward boundary being a vertical wall are analyzed in this paper by using the method of separation of variables. Analytical expressions for the radiated and diffracted potentials are derived as infinite series with unknown coefficients determined by the eigenfunction expansion matching method. The expressions for wave forces and hydrodynamic coefficients are given. A comparison is made between the results obtained by the present analytical solution and those obtained by the boundary element method. By using the present analytical solution, the hydrodynamic influences of the submergence, the width, the thickness of the structure, and the distance between the structure and the wall on the wave forces and hydrodynamic coefficients are discussed in detail.     

The vertical mode method in the problems of flexural-gravity waves diffracted by a vertical cylinder

The linear three-dimensional problem of ice loads acting on a vertical circular cylinder frozen in an ice cover of infinite extent is studied. The loads are caused by an uni-directional hydroelastic wave propagating in the ice cover towards the cylinder mounted to the see bottom in water of constant depth. There are no open water surfaces in this problem. The deflection of the ice cover is described by the Bernoulli–Euler equation of a thin elastic plate of constant thickness. At the contact line between the ice cover and the surface of the cylinder, some edge conditions are imposed. In this study, the edge of the ice plate is either clamped to the cylinder or has no contact with the cylinder surface, with the plate edge being free of stresses and shear forces. The water is of finite constant depth, inviscid and incompressible. The problem is solved by both the vertical mode method and using the Weber integral transform in the radial coordinate. Each vertical mode corresponds to a root of the dispersion relation for flexural-gravity waves. It is proved that these two solutions are identical for the clamped edge conditions. This result is non-trivial because the vertical modes are non-orthogonal in a standard sense, they are linearly dependent, the roots of the dispersion relation can be double and even triple, and the set of the modes could be incomplete. A general solution of the wave-cylinder interaction problem is derived by the method of vertical modes and applied to different edge conditions on the contact line. There are three conditions of solvability in this problem. It is shown that these conditions are satisfied for any parameters of the problem.     

Nonlinear hydroelastic analysis of a moored floating body

By integration of the second-order fluid pressure over the instantaneous wetted surface, the generalized first- and second-order fluid forces used in nonlinear hydroelastic analysis are obtained. The expressions for coefficients of the generalized first- and second-order hydrodynamic forces in irregular waves are also given. The coefficients of the restoring forces of a mooring system acting on a flexible floating body are presented. The linear and nonlinear three-dimensional hydroelastic equations of motion of a moored floating body in frequency domain are established. These equations include the second-order forces, induced by the rigid body rotations of large amplitudes in high waves, the variation of the instantaneous wetted surface and the coupling of the first order wave potentials. The first-order and second-order principal coordinates of the hydrelastic vibration of a moored floating body are calculated. The frequency characteristics of the principal coordinates are discussed. The numerical results indicate that the rigid resonance and the coupling resonance of a moored floating body can occur in low frequency domain while the flexible resonance can occur in high frequency domain. The hydroelastic responses of a moored box-type barge are also given in this paper. The effects of the second-order forces on the modes are investigated in detail.     

Gravity wave interaction with floating membrane due to abrupt change in water depth

Using the recently developed expansion formulae for wave structure interaction problems, the scattering of surface water waves by a semi-infinite floating membrane due to abrupt change in bottom topography is analyzed. Both the cases of finite and infinite steps are analyzed. In the present paper, the analysis is based on the linearized theory of water waves and small amplitude membrane response. Combining the linearized kinematic and dynamic surface conditions on the water surface with the dynamic pressure condition on the membrane, a third order differential equation is derived to describe the membrane covered free surface condition. General wave energy relation for wave scattering by floating horizontal membrane is derived by the application of law of conservation of energy flux and alternately by the direct application of Green's second identity. In the floating membrane covered region, the wave energy density is a combination of the kinetic and potential energy density due to the surface gravity waves, and the surface energy density which is due to the existence of the floating membrane on the free surface. Gravity wave transformations due to an abrupt change in bottom topography in the presence of a floating membrane in finite water depth are analyzed based on shallow water approximation. Numerical results are computed and analyzed to understand the wave transformation due to the floating membrane when there is an abrupt change in topography in different cases.     

Transient Hydroelastic Response of VLFS by FEM with Impedance Boundary Conditions in Time Domain

A time-dependent finite element method (FEM) is developed to analyze the transient hydroelastie responses of very large floating structures (VLFS) subjected to dynamic loads. The hydrodynamic problem is formulated based on the linear theory of fluid and the structural response is analyzed based on the thin plate theory. The FEM truncates the unbounded fluid domain by introducing an artificial boundary surface, thus defining a finite computational domain. At this boundary surface an impedance boundary conditions are applied so that no wave reflections occur. In the proposed scheme, all of the procedures are processed directly in time domain, which is efficient for nonlinear analyses of structure floating on unbounded fluid. Numerical results indicate acceptable accuracy of the proposed method.     

Hydroelastic Analysis of Very Large Floating Structures Using Plate Green Functions

1 .IntroductionRecentlygreatinteresthasbeenshowninthedevelopmentofverylargefloatingstructuressuchasMegaFloatofJapan (Isobe ,1 999)andMOBofUSA (Remmers ,1 999) .Owingtotheirextremelargesizeandgreatflexibility ,thecouplingbetweenthestructuraldeformationandfluidmotionissignifi cant.Thisisatypicalproblemofhydroelasticity .Efficientandaccurateestimationofthehydroelasticresponseofverylargefloatingstructuresinwavesisveryimportantfordesign .Manymethodshavebeenproposedinliteratureforthepredictiono…     

Hydroelasticity of a floating plate in multidirectional waves

The membrane forces are included in the hydroelastic analysis of a floating plate undergoing large vertical deflections in regular monochromatic multidirectional waves. The first-order vertical displacements induced by the linear wave exciting forces are calculated by the mode expansion method in the frequency domain. The second-order vertical displacements induced by the membrane forces are calculated by the von Karman plate theory. The results show that the membrane contribution both in terms of the axial stresses and the effect on the bending stresses can be important.     

On the radiation and diffraction of linear water waves by an infinitely long rectangular structure submerged in oblique seas

The radiation and diffraction of linear water waves by an infinitely long rectangular structure submerged in oblique seas of finite depth is investigated. The analytical expressions for the radiated and diffracted potentials are derived as infinite series by use of the method of separation of variables. The unknown coefficients in the series are determined by the eigenfunction expansion matching method. The expressions for wave forces, hydrodynamic coefficients and reflection and transmission coefficients are given and verified by the boundary element method. Using the present analytical solution, the hydrodynamic influences of the angle of incidence, the submergence, the width and the thickness of the structure on the wave forces, hydrodynamic coefficients, and reflection and transmission coefficients are discussed in detail.     

The hydrodynamics of thin floating plates

The radiation and diffraction problems are considered in the frequency domain for a thin elastic plate of rectangular planform floating in an irrotational, incompressible ocean of infinite depth. The inner potential field inside a hemisphere surrounding the plate is represented using a spherical harmonic expansion which suits the geometry and zero-draft nature of the plate. Problems associated with distributing sources in the free surface are avoided. The Chen and Mei variational principle is used to weakly match this inner solution and its normal derivative to an outer field described by distributing sources on the exterior of the hemisphere. The validity of the procedure is first illustrated by considering a heaving circular disk. Numerous hydrodynamic coefficients are presented as benchmark data for floating flexible structures. The transient motion of the plate is simulated using rational approximations (in the frequency domain) to the radiation impedance and diffraction mapping which are implemented as ODE's in the time domain.     

Performance of pile-restrained flexible floating breakwaters

The overall performance of pile-restrained flexible floating breakwaters is investigated under the action of linear monochromatic incident waves in the frequency domain. The aforementioned floating breakwaters undergo only vertical structural deflections along their length and are held in place by means of vertical piles. The total number of degrees of freedom equals the six conventional body modes, when the breakwater moves as a rigid body, plus the extra bending modes. These bending modes are introduced to represent the structural deflections of the floating breakwater and are described by the Bernoulli–Euler flexible beam equation. The number of bending modes introduced is determined through an appropriate iterative procedure. The hydrostatic coefficients corresponding to the bending modes are also derived. The numerical analysis of the flexible floating breakwaters is based on a three-dimensional hydrodynamic formulation of the floating body. A parametric study is carried out for a wide range of structural stiffness parameters and wave headings, to investigate their effect on the performance of flexible floating breakwaters. Moreover, this performance is compared with that of the corresponding pile-restrained rigid floating breakwater. Results indicated that the degree of structural stiffness and the wave heading strongly affect the performance of flexible floating breakwaters. The existence of an “optimum” value of structural stiffness is demonstrated for the entire wave frequency range.     

Multipole expansions for wave diffraction and radiation in deep water

A multipole expansion of the velocity potential is described for two- and three-dimensional wave diffraction and radiation problems. The velocity potential is expressed in terms of a series of multipole potentials. The wave terms and the local disturbance terms are represented by separated multipole potentials. Floating bodies and submerged bodies are treated in the same way. This approach differs from that of some other authors, who considered floating bodies and submerged bodies separately and derived entirely different multipoles. Semi-analytical solutions for a circular cylinder in two-dimensional motions are given. It is found that the local disturbance decays rapidly and steadily. The general application of the multipole expansion to arbitrary geometries is also presented, based on a method coupling multipoles to a boundary integral expression. Numerical results for several floating and submerged cylinders are presented.     

On the Interaction of Surface Waves with An Elastic Plate of Finite Length in Head Seas

An eigen-function expansion method based on a new orthogonal inner product is proposed by Sahoo et al. (2000) for the study of the hydroelastic response of mat-type VLFS in head seas. However, their main emphasis is on the effect of edge conditions and they assume that the plate is of a semi-infinite length. In reality, the plate is of finile length. For consider-alion of the finite length effect, the reflection and transmission from the other end must be considered. The effect of this reflection and transmission on the hydroelaslic response of VLFS is of interest for praclical application. Furthermore, the physi-cal meaning of the new inner producl was not given in their paper. In this paper, it is shown that the new inner product can be derived from the governing equation and the bottom boundary conditions. Then the same eigen-function expansion method is adopted for the study of the hydroelastic response of an elastic plate of finite length in surface waves. Delailed comparisons are made between the