采用基于雷诺平均的N-S方程和势流方法进行散货船阻力和运动的预报研究（英文）

Resistance prediction of ships using computational fluid dynamics has gained popularity over the years because of its high accuracy and low cost. This paper conducts numerical estimations of the ship resistance and motion of a Japan bulk carrier model using SHIP_Motion, a Reynolds-averaged Navier–Stokes(RaNS)-based solver, and HydroSTAR, a commercial potential flow(PF)-based solver. The RaNS solver uses an overset-structured mesh and discretizes the flow field using the finite volume method, while the PF-based solver applies the three-dimensional panel method. In the calm water test, the total drag coefficient,sinkage, and trim were predicted using the RaNS solver following mesh dependency analysis, and the results were compared with the available experimental data. Next, calm water resistance was investigated for a range of Froude numbers. The added resistance in short-wave cases was simulated using both Ra NS and PF solvers, and the results were compared. The PF solver showed better agreement with the RaNS solver for predicting motion responses than for predicting added resistance. While the added resistance results could not be directly validated because of the absence of experimental data, considering the previous accuracy of the RaNS solver in added resistance prediction and general added resistance profile of similar hull forms(bulk carriers), the prediction results could be concluded to be reliable.     

Optimization of bow shape for a non ballast water ship

In this research,a commercial CFD code ＂Fluent＂ was applied to optimization of bulbous bow shape for a non ballast water ships（NBS）.The ship was developed at the Laboratory of the authors in Osaka Prefecture University,Japan.At first,accuracy of the CFD code was validated by comparing the CFD results with experimental results at towing tank of Osaka Prefecture University.In the optimizing process,the resistances acting on ships in calm water and in regular head waves were defined as the object function.Following features of bulbous bow shapes were considered as design parameters： volume of bulbous bow,height of its volume center,angle of bow bottom,and length of bulbous bow.When referring to the computed results given by the CFD like resistance,pressure and wave pattern made by ships in calm water and in waves,an optimal bow shape for ships was discovered by comparing the results in the series of bow shapes.In the computation on waves,the ship is in fully captured condition because shorter waves,λ/Lpp 0.6,are assumed.     

片体间距对高速三体船水动力性能影响的比较分析(英文)

This paper numerically investigates the influence of separation variation of the outriggers on the hydrodynamic performance of a high speed trimaran (HST) aiming at improving its applicability in diverse realistic disciplines. The present investigation was performed within the framework of the 2-D slender body method (SBM) by calculating the resistance of three symmetric trimaran series moving in a calm free surface of deep water. Each trimaran series comprises of 4681 configurations generated by considering 151 staggers (50%≤α≤+100%), and 31 separations (100%≤β≤400%) for 81 Froude numbers (0.20≤Fn≤1.0). In developing the three trimaran series, Wigley-st, AMECRC-09, and NPL-4a models were used separately for both the main and side hulls of each individual series models. A computer macro named Tri-PL was created using the Visual Basic for Applications. Tri-PL sequentially interfaced Maxsurf then Hullspeed to generate the models of the three trimaran series together with their detailed hydrostatic particulars, followed by their resistance components. The numerical results were partially validated against the available published numerical calculations and experimental results, to benchmark the Tri-PL macro and hence to rely on the analysis outcomes. A graph template was created within the framework of SigmaPlot to visualize the significant results of the Tri-PL properly.     

浮式海上风力机运动性能和锚泊系统(英文)

The development of offshore wind farms was originally carried out in shallow water areas with fixed(seabed mounted) structures.However,countries with limited shallow water areas require innovative floating platforms to deploy wind turbines offshore in order to harness wind energy to generate electricity in deep seas.The performances of motion and mooring system dynamics are vital to designing a cost effective and durable floating platform.This paper describes a numerical model to simulate dynamic behavior of a new semi-submersible type floating offshore wind turbine(FOWT) system.The wind turbine was modeled as a wind block with a certain thrust coefficient,and the hydrodynamics and mooring system dynamics of the platform were calculated by SESAM software.The effect of change in environmental conditions on the dynamic response of the system under wave and wind loading was examined.The results indicate that the semi-submersible concept has excellent performance and SESAM could be an effective tool for floating wind turbine design and analysis.     

Two-dimensional FEM model of Vortex-Induced Vibration of A Circular Cylinder

The finite element method was employed to solve the N-S equation.Incorporated with the vibration equa- tion,the vortex-induced vibration of the circular cylinder is studied by an Arbitrary-Lagrangian-Eulerian(ALE)al- gorithm.The Reynolds number based on the cylinder diameter and the undisturbed flow velocity ranges from 90 to 150. The motion of the cylinder was modeled by a spring-damper-mass system.The numerical model has been validated by comparison with the experimental data in literature.The"lock-in"and"beating"phenomena were successfully repro- duced in the numerical test.     

Experimental and numerical procedures of a sonar platform with a sound absorption wedge

Experiments involving a sonar platform with a sound absorption wedge were carried out for the purpose of obtaining the low frequency acoustic characteristics. Acoustic characteristics of a sonar platform model with a sound absorption wedge were measured, and the effects of different wedge laid areas on platform acoustic characteristic were tested. Vibration acceleration and self-noise caused by model vibration were measured in four conditions: 0%, 36%, 60%, and 100% of wedge laid area when the sonar platform was under a single frequency excitation force. An experiment was performed to validate a corresponding numerical calculation. The numerical vibration characteristics of platform area were calculated by the finite element method, and self-noise caused by the vibration in it was predicted by an experiential formula. The conclusions prove that the numerical calculation method can partially replace the experimental process for obtaining vibration and sound characteristics.     

Dynamic interaction of parallel moving ships in close proximity

Nowadays,there are many studies conducted in the field of marine hydrodynamics which focus on two vessels traveling and floating in sufficiently close proximity to experience significant interactions.The hydrodynamic behavior of parallel moving ships in waves is an interesting and important topic of late.A numerical investigation has been carried out for the prediction of wave exciting forces and motion responses of parallel moving ships in regular waves.The numerical solution was based on 3D distribution technique and using the linear wave theory to determine the exciting forces and ship＇s motion.The speed effects have been considered in the Green function for more realistic results.The numerical computations of wave exciting forces and motion responses were carried out for a Mariner and Series 60 for the purpose of discovering different Froude numbers and different separation distances in head sea conditions.Based on the numerical computations,it was revealed that the sway,roll and yaw have a significant effect due to hydrodynamic interaction.     

A method for numerical calculation of propeller hydrodynamics in unsteady inflow

The hydrodynamic performance of a propeller in unsteady inflow was calculated using the surface panel method. The surfaces of blades and hub were discreted by a number of hyperboloidal quadrilateral panels with constant source and doublet distribution. Each panel's corner coordinates were calculated by spline interpolation between the main parameter and the blade geometry of the propeller. The integral equation was derived using the Green Formula. The influence coefficient of the matrix was calculated by the Morino analytic formula. The tangential velocity distribution was calculated with the Yanagizawa method, and the pressure coefficient was calculated using the Bonuli equation. The pressure Kutta condition was satisfied at the trailing edge of the propeller blade using the Newton-Raphson iterative procedure, so as to make the pressure coefficients of the suction and pressure faces of the blade equal at the trailing edge. Calculated results for the propeller in steady inflow were taken as initialization values for the unsteady inflow calculation process. Calculations were carried out from the moment the propeller achieved steady rotation. At each time interval, a linear algebraic equation combined with Kutta condition was established on a key blade and solved numerically. Comparison between calculated results and experimental results indicates that this method is correct and effective.     

弹性圆柱壳结构的有源结构声辐射控制研究

A numerical and experimental study was presented on active control of structurally radiated sound from an elastic cylindrical shell.An analytical model was developed for the active structural acoustic control (ASAC) of the cylindrical shell.Both global and local control strategies were considered.The optimal control forces corresponding to each control strategy were obtained by using the linear quadratic optimal control theory.Numerical simulations were performed to examine and analyze the control performance under different control strategies.The results show that global sound attenuation of the cylindrical shell at resonance frequencies can be achieved by using point force as the control input of the ASAC system.Better control performance can be obtained under the control strategy of minimization of the radiated sound power.However,control spillover may occur at off-resonance frequencies with the control strategy of structural kinetic energy minimization in terms of the radiated sound power.Considerable levels of global sound attenuation can also be achieved in the on-resonance cases with the local control strategy,i.e.,minimization of the mean-square velocity of finite discrete locations.An ASAC experiment using an FXLMS algorithm was implemented,agreement was observed between the numerical and experimental results,and successful attenuation of structural vibration and radiated sound was achieved.     

基于改进的Logvinovich法对带横摇运动非对称楔形体入水砰击力的计算(英文)

The water entry problem of an asymmetric wedge with roll motion was analyzed by the method of a modified Logvinovich model (MLM). The MLM is a kind of analytical model based on the Wagner method, which linearizes the free surface condition and body boundary condition. The difference is that the MLM applies a nonlinear Bernoulli equation to obtain pressure distribution, which has been proven to be helpful to enhance the accuracy of hydrodynamic loads. The Wagner condition in this paper was generalized to solve the problem of the water entry of a wedge body with rotational velocity. The comparison of wet width between the MLM and a fully nonlinear numerical approach was given, and they agree well with each other. The effect of angular velocity on the hydrodynamic loads of a wedge body was investigated.     

Numerical simulation of the anti-shock performance of a gear case

Both sea battles and testing of ship in underwater explosions reveal unacceptably poor anti-shock performance of important shipboard equipment. Anti-shock performance of shipboard equipment is a significant factor determining fighting strength and survivability. The anti-shock performance of a shipboard gear case based on BV043/85 was investigated using numerical simulation. A geometric model of the gear case was built using MDT software and meshed in HyperMesh software, and then the finite element model of the gear case was formed. Using ABAQUS software, the anti-shock performance of the gear case was simulated. First, shock response of typical regions of gear case was determined. Next, some generalizations were made about the anti-shock performance of the gear case by analyzing the Mises stress of typical regions varied with shock inputs. Third, weak regions were determined from simulation results. The threshold values of shock resistance of the gear case at different impulse widths were obtained through interpolating the numerical simulation results selected from the most dangerous spot. This research provides a basis for further optimization of the design of gear cases.     

Numerical simulation of the water entry of a structure in free fall motion

To solve the problems concerning water entry of a structure, the RANS equations and volume of fluid （VOF） method are used. Combining the user-defined function （UDF） procedure with dynamic grids, the water impact on a structure in free fall is simulated, and the velocity, displacement and the pressure distribution on the structure are investigated. The results of the numerical simulation were compared with the experimental data, and solidly consistent results have been achieved, which validates the numerical model. Therefore, this method can be used to study the water impact problems of a structure.     

能量耗散效应的多域边界元法(英文)

The wave diffraction and radiation around a floating body is considered within the framework of the linear potential theory in a fairly perfect fluid.The fluid domain extended infinitely in the horizontal directions but is limited by the sea bed,the body hull,and the part of the free surface excluding the body waterplane,and is subdivided into two subdomains according to the body geometry.The two subdomains are connected by a control surface in fluid.In each subdomain,the velocity potential is described by using the usual boundary integral representation involving Green functions.The boundary integral equations are then established by satisfying the boundary conditions and the continuous condition of the potential and the normal derivation across the control surface.This multi-domain boundary element method(MDBEM) is particularly interesting for bodies with a hull form including moonpools to which the usual BEM presents singularities and slow convergence of numerical results.The application of the MDBEM to study the resonant motion of a water column in moonpools shows that the MDBEM provides an efficient and reliable prediction method.     

Collision analysis of the spar upper module docking

In order to assess the possible collision effect, a numerical simulation for the upper module and spar platlbrm docking at the speed of 0.2 m/s was conducted by using the software ANSYS/LS-DYNA, and the time history of the collision force, energy absorption and structural defonamtion during the collision was described. The purpose was to ensure that the platlbrm was safely put into operation. Furthermore, this paper analyzes different initial velocities and angles on the Von Mises stress and collision resultant force during the docking collision. The results of this paper showed that the docking could be conducted with higher security. The data in this paper can provide useful references for the determination of the upper module＇s offshore hoisting scheme and practical construction by contrasting the numerical simulation results of the parameters on the docking collision.     

逆压边界层流动的湍流模式研究（英文）

In this paper, a numerical study of flow in the turbulence boundary layer with adverse and pressure gradients(APGs) is conducted by using Reynolds-averaged Navier-Stokes(RANS) equations. This research chooses six typical turbulence models, which are critical to the computing precision, and to evaluating the issue of APGs. Local frictional resistance coefficient is compared between numerical and experimental results. The same comparisons of dimensionless averaged velocity profiles are also performed. It is found that results generated by Wilcox(2006) k-ω are most close to the experimental data. Meanwhile, turbulent quantities such as turbulent kinetic energy and Reynolds-stress are also studied.     

Turbulence Model Investigations on the Boundary Layer Flow with Adverse Pressure Gradients

In this paper, a numerical study of flow in the turbulence boundary layer with adverse and pressure gradients(APGs) is conducted by using Reynolds-averaged Navier-Stokes(RANS) equations. This research chooses six typical turbulence models, which are critical to the computing precision, and to evaluating the issue of APGs. Local frictional resistance coefficient is compared between numerical and experimental results. The same comparisons of dimensionless averaged velocity profiles are also performed. It is found that results generated by Wilcox(2006) k-ω are most close to the experimental data. Meanwhile, turbulent quantities such as turbulent kinetic energy and Reynolds-stress are also studied.     

基于B样条应用的船舶波浪增阻计算（英文）

Making an exact computation of added resistance in sea waves is of high interest due to the economic effects relating to ship design and operation. In this paper, a B-spline based method is developed for computation of added resistance. Based on the potential flow assumption, the velocity potential is computed using Green's formula. The Kochin function is applied to compute added resistance using Maruo's far-field method, the body surface is described by a B-spline curve and potentials and normal derivation of potentials are also described by B-spline basis functions and B-spline derivations. A collocation approach is applied for numerical computation, and integral equations are then evaluated by applying Gauss–Legendre quadrature. Computations are performed for a spheroid and different hull forms; results are validated by a comparison with experimental results. All results obtained with the present method show good agreement with experimental results.     

波浪中船体梁后极限强度特性的实验评估(英文)

Experimental investigations into the collapse behavior of a box-shape hull girder subjected to extreme wave-induced loads are presented.The experiment was performed using a scaled model in a tank.In the middle of the scaled model,sacrificial specimens with circular pillar and trough shapes which respectively show different bending moment-displacement characteristics were mounted to compare the dynamic collapse characteristics of the hull girder in waves.The specimens were designed by using finite element(FE)-analysis.Prior to the tank tests,static four-point-bending tests were conducted to detect the load-carrying capacity of the hull girder.It was shown that the load-carrying capacity of a ship including reduction of the capacity after the ultimate strength can be reproduced experimentally by employing the trough type specimens.Tank tests using these specimens were performed under a focused wave in which the hull girder collapses under once and repetitive focused waves.It was shown from the multiple collapse tests that the increase rate of collapse becomes higher once the load-carrying capacity enters the reduction path while the increase rate is lower before reaching the ultimate strength.     

The asymptotic field of a dynamically growing crack in a viscoelastic materia

A mechanical model of a fracturing viscoelastic material was developed to investigate viscous effects in a dynamically growing crack-tip field. It was shown that in the stable creep-growing phase,elastic deformation and viscous deformation are equally dominant in the near-tip field,and stress and strain have the same singularity,namely,(σ,ε ) ∝ r?1 /( n?1) . The asymptotic solution of separating variables of stress,stain and displacement in the crack-tip field was obtained by asymptotic analysis,and the resulting numerical value of stress and strain in the crack-tip field was obtained by the shooting method and the boundary condition of a mode I crack. Through numerical calculation,it was shown that the near-tip fields are mainly governed by the creep exponent n and Mach number M . When n →∞,the asymptotic solution of a viscoelastic material can be degenerated into that of Freund's elastic-ideally plastic material by analyzing basic equations.