A simple formulation for predicting the ultimate strength of ships

The aim of this study is to derive a simple analytical formula for predicting the ultimate collapse strength of a single- and double-hull ship under a vertical bending moment, and also to characterize the accuracy and applicability for earlier approximate formulations. It is known that a ship hull will reach the overall collapse state if both collapse of the compression flange and yielding of the tension flange occur. Side shells in the vicinity of the compression and the tension flanges will often fail also, but the material around the final neutral axis will remain in the elastic state. Based on this observation, a credible distribution of longitudinal stresses around the hull section at the overall collapse state is assumed, and an explicit analytical equation for calculating the hull ultimate strength is obtained. A comparison between the derived formula and existing expressions is made for largescale box girder models, a one-third-scale frigate hull model, and full-scale ship hulls.List of symbols A _B total sectional area of outer bottom - A _B total sectional area of inner bottom - A _D total sectional area of deck - A _S half-sectional area of all sides (including longitudinal bulkheads and inner sides) - a _s sectional area of a longitudinal stiffener with effective plating - b breadth of plate between longitudinal stiffeners - D hull depth - D _B height of double bottom - E Young's modulus - g neutral axis position above the base line in the sagging condition or below the deck in the hogging condition - H depth of hull section in linear elastic state - I _s moment of inertia of a longitudinal stiffener with effective plating - l length of a longitudinal stiffener between transverse beams - M _E elastic bending moment - M _p fully plastic bending moment of hull section - M _u ultimate bending moment capacity of hull section - M _uh ,M _us ultimate bending moment in hogging or sagging conditions - r radius of gyration of a longitudinal stiffener with effective plating [=(I _s /a _s )^1/2] - t plate thickness - Z elastic section modulus at the compression flange - Z _B ,Z _D elastic section modulus at bottom or deck - slenderness ratio of plate between stiffeners [= (b/t)(_y/E)^1/2] - slenderness ratio of a longitudinal stiffener with effective plating [=(l/r)(_y/E)^1/2] - _y yield strength of the material - _yB , _yB , _yD yield strength of outer bottom, inner bottom - _yS deck, or side - _u ultimate buckling strength of the compression flange - _uB , _uB , _uD ultimate buckling strength of outer bottom - _uS inner bottom, deck, or side     

Surf-riding and oscillations of a ship in quartering waves

The behavior of a ship encountering large regular waves from astern at low frequency is the object of investigation, with a parallel study of surf-riding and periodic motion paterns. First, the theoretical analysis of surf-riding is extended from purely following to quartering seas. Steady-state continuation is used to identify all possible surf-riding states for one wavelength. Examination of stability indicates the existence of stable and unstable states and predicts a new type of oscillatory surf-riding. Global analysis is also applied to determine the areas of state space which lead to surf-riding for a given ship and wave conditions. In the case of overtaking waves, the large rudder-yaw-surge oscillations of the vessel are examined, showing the mechanism and conditions responsible for loss of controllability at certain vessel headings.List of symbols c wave celerity (m/s) - C(p) roll damping moment (Ntm) - g acceleration of gravity (m/s²) - GM metacentric height (m) - H wave height (m) - I _x ,I _z roll and yaw ship moments of inertia (kg m²) - k wave number (m^–1) - K _H ,K _W ,K _R hull reaction, wave, rudder, and propeller - K _p forces in the roll direction (Ntm) - m ship mass (kg) - n propeller rate of rotation (rpm) - N _H ,N _W ,N _R hull reaction, wave, rudder, and propeller - N _P moments in the yaw direction (Ntm) - p roll angular velocity (rad/s) - r rate-of-turn (rad/s) - R(,x) restoring moment (Ntm) - Res(u) ship resistance (Nt) - t time (s) - u surge velocity (m/s) - U vessel speed (m/s) - v sway velocity (m/s) - W ship weight (Nt) - x longitudinal position of the ship measured from the wave system (m) - x _G ,z _G longitudinal and vertical center of gravity (m) - x _S longitudinal position of a ship section (S), in the ship-fixed system (m) - X _H ,X _W ,X _R hull reaction, wave, rudder, and propeller - X _P forces in the surge direction (Nt) - y transverse position of the ship, measured from the wave system (m) - Y _H ,Y _W ,Y _R hull reaction, wave, rudder, and propeller - Y _p forces in the sway direction (Nt) - z _Y vertical position of the point of action of the lateral reaction force during turn (m) - z _W vertical position of the point of action of the lateral wave force (m) Greek symbols angle of drift (rad) - rudder angle (rad) - wavelength (m) - position of the ship in the earth-fixed system (m) - water density (kg/m³) - angle of heel (rad) - heading angle (rad) - _e frequency of encounter (rad/s) Hydrodynamic coefficients K roll added mass - N _v ,N _r yaw acceleration coefficients - N _v N _r N _rr N _rrv ,N _vvr yaw velocity coefficients K. Spyrou: Ship behavior in quartering waves - X _u surge acceleration coefficient - X _u X _vr surge velocity coefficients - Y _v ,Y _r sway acceleration coefficients - Y _v ,Y _r ,Y _vv ,Y _rr ,Y _vr sway velocity coefficients European Union-nominated Fellow of the Science and Technology Agency of Japan, Visiting Researcher, National Research Institute of Fisheries Engineering of Japan     

边坡滑动面三维空间有限元分析

提出了一种用三维潜在滑移面法分析三维边坡滑动面及稳定性的新方法,该方法根据弹塑性有限元的应力分析结果,用数值积分方法确定边坡的潜在滑移面、最危险潜在滑移面和边坡稳定性安全系数,具有充分的理论依据。通过对大量三维边坡实例的稳定性安全系数及最危险滑动面的分析得出,三维边坡按平面边坡处理会引起较大的分析误差,三维边坡应该按三维问题进行分析。     

基于生长算法构建Delaunay三角网的研究

提出了一种基于凸包构建Delaunay三角网的算法,它包括离散点数据的三角剖分及将约束边嵌入三角网。在对离散点数据进行三角剖分时,首先建立一个外轮廓为凸包的初始三角网,然后将剩余的点依次与既有三角网合并来生成包含所有数据点的三角网。为了将约束边嵌入三角网,找出约束边穿过的三角形,删除这些三角形并用约束边将这些三角形所围成的多边形分成两个简单多边形,再对这两个简单多边形进行Delaunay三角网剖分即可达到嵌入约束边的目的。     

基于自适应遗传算法的路堤边坡稳定性分析方法

基于圆弧滑动面假定,提出了一种用自适应遗传算法搜索最危险滑动面及其对应的最小安全系数的新方法。该方法是一种改进的遗传算法,采用自适应求取适值、动态调整交叉率和变异率、自适应区间收缩。自适应遗传算法有效克服了传统方法易陷入局部极小的缺陷,提高了算法的搜索效率、精度和稳定性。     

山区公路层间滑移破坏机理分析及对策研究

通过现场调查,初步分析了山区公路沥青路面基面层滑移的主要原因。应用有限元法,对基层和面层之间的剪应力进行计算分析,得出了基面层间剪应力变化的一般规律,认为层间最大剪应力主要取决于面层厚度和荷载条件。在前述计算分析的基础上,利用自行研制的DLG A路面材料剪切仪进行了大量的层间滑移室内模拟试验,考察了不同的层间处理措施、温度、基层类型等因素对层间抗剪能力的影响,提出了防止基面层滑移的技术措施。     

Simulation on progressive collapse behaviour of whole ship model under extreme waves using idealized structural unit method

To ensure the safety of navigating ship, working loads and structural load-carrying capacity are two important aspects. In the present paper, a total simulation system combing load calculation and structural collapse analysis is applied to simulate progressive collapse behaviour of a single-hull Kamsarmax type bulk carrier. A three dimensional singularity distribution method is adopted to calculate pressure distribution with time history. A mixed structural model, collapse part simulated by ISUM elements and remaining part by elastic FEM elements with relative coarse mesh, is proposed for collapse analysis. Progressive collapse behaviour obtained by ISUM is good agreement with that by nonlinear software package, MARC. However, the calculation time of ISUM analysis is about 1/70 of MARC analysis. The applicability to structure system, high accuracy and sufficient efficiency of ISUM had been demonstrated.     

Tuning of transfer functions for analysis of wave–ship interactions

Transfer functions are often used together with a wave spectrum for analysis of wave–ship interactions, where one application addresses the prediction of wave-induced motions or other types of global responses. This paper presents a simple and practical method which can be used to tune the transfer function of such responses to facilitate improved prediction capability. The input to the method consists of a measured response, i.e. time series sequences from a given sensor, the 2D wave spectrum characterising the seaway in which the measurements are taken, and an initial estimate of the transfer function for the response in study. The paper presents results obtained using data from an in-service container ship. The 2D wave spectra are taken from the ERA5 database, while the transfer function is computed by a simple closed-form expression. The paper shows that the application of the tuned transfer function leads to predictions which are significantly improved compared to using the transfer function without tuning.     

基于ANN的船舶撞击高桩码头群桩损伤位置预测

本文采用有限元软件ABAQUS建立了船舶撞击高桩码头群桩的空间有限元模型。通过计算评估了撞击力、桩体刚度、撞击位置和撞击角度下对群桩结构损伤位置的影响。基于人工神经网络(ANN)方法,对不同参数组合下的群桩结构损伤位置进行了预测,并对ANN方法的可行性进行了评估。     

考虑横向漂移的船舶路径跟踪自抗扰控制

针对船舶运动系统中内部动态不确定和外部干扰等问题,进行了欠驱动船舶路径跟踪的自抗扰方法研究。利用Backstepping设计参考航向角,并通过线性扩张状态观测器对流干扰和横向运动引起的横向漂移进行估计。其次,根据自抗扰算法对航向进行控制,采用线性扩张状态观测器对外界干扰及内部不确定项进行估计。最后仿真结果表明,在风流干扰下所设计的控制器仍能使船准确地跟踪上参考路径,验证了所提控制方案的有效性。