Next Article in Journal
Sliding Mode Motion Control for AUV with Dual-Observer Considering Thruster Uncertainty
Previous Article in Journal
Mediterranean Coastal Lagoons Review: Sites to Visit before Disappearance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Probabilistic Damage Stability for Passenger Ships—The p-Factor Illusion and Reality

1
Maritime Safety Research Centre (MSRC), Department of Naval Architecture, Ocean and Marine Engineering, University of Strathclyde, Glasgow G4 0LZ, UK
2
Bureau Veritas, Marine and Offshore, 92937 Paris, France
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2022, 10(3), 348; https://doi.org/10.3390/jmse10030348
Submission received: 21 January 2022 / Revised: 24 February 2022 / Accepted: 25 February 2022 / Published: 1 March 2022
(This article belongs to the Section Ocean Engineering)

Abstract

:
The paper complements an earlier publication by the authors addressing the probability of survival in the IMO framework for damage stability assessment, the s-factor. The focus here is on the probability of occurrence of a certain damage scenario (breach), conditional on its dimensions and location (centre and port or starboard side), the p-factor. Pertinent assumptions and limitations are explained, following its evolution for specific application to passenger ships. Attempts to provide analytical descriptions of the damage breach distributions as tetrahedra shapes positioned along the ship length whilst accounting for changes in ship geometry, structural arrangements, and subdivision for consumption by the wider profession has led to misconceptions and misunderstandings of what exactly the p-factor is in the context of probabilistic damage stability calculations. This is evidenced by the fact that the same original damage breach distributions, derived in Project HARDER, based on largely cargo ships with the age spread over the last three decades of the previous century, are still being used today for all ship types, including modern passenger ships. Filling this gap, a new database for passenger ships developed in the EC-funded Project FLARE, is briefly presented, leading to new damage breach distributions specifically for passenger ships. It is believed that this paper will throw considerable light in enhancing understanding on the p-factor, which has been cluttered with unnecessary complexity from the outset.

1. Introduction

The probabilistic assessment of ship survivability after an accident should be a comprehensive process estimating the conditional probability of losing ship stability in the wake of a casualty. Even though the definition of a probabilistic framework has developed for the last 30 years, the actual regulations imposed by SOLAS 2009/2020, ref. [1], incorporate just a few elements of the provided research output. The only cause of accident included in the SOLAS framework is collisions, totally neglecting other sources of hazards for ships, such as groundings, that could be more frequent and dangerous for certain kinds of vessels as passenger ships. Furthermore, SOLAS provides a classification of the safety level of a ship based on the evaluation of indices instead of promoting a direct approach for the estimation of flooding risk. In such a case, the resulting probabilistic assessment neglects relevant aspects for ship survivability, such as the operational area and operating environment, the structural arrangements, the breaches definition and distributions and the vessel type. More importantly, focussing on indices, as a substitute for the direct assessment of flooding risk, deprives such assessments of the time element, hence crucial information on measures affects the improvements on the evolution of flooding leading to capsize as well as evacuation arrangements and associated Risk Control Options (RCOs), affecting evacuation in such scenarios.
On the other hand, SOLAS regulation provides a clear logic to evaluate ship survivability through an Attained Subdivision Index (A-Index):
A = j = 1 J i = 1 I w j p i s i
where,
  • j —represents the loading condition under consideration.
  • J —represents the total number of loading conditions considered in the calculation of A , usually three draughts covering the operational draught range of the vessel.
  • w j —represents a weighting factor applied to each initial draught.
  • i —represents each compartment or group of compartments under consideration for the loading condition, j .
  • I —is the total number of all feasible damage scenarios involving flooding of individual compartments or groups of adjacent compartments.
  • p i —is the probability that, for the loading condition, j , only the compartment or group of compartments under consideration are flooded, disregarding any horizontal subdivision.
  • s i —accounts for the conditional probability of survival following flooding of the compartment or group of compartments under consideration for loading condition j , weighted by the probability that the space above a horizontal subdivision may not be flooded.
The use of the A-Index as a safety measure gives a fully decoupled approach for the determination of the flooding probability (p-factor) and ship survivability (s-factor), as it was clear since the first studies of Wendel on the probabilistic damage stability assessment [2,3]. This simple but efficient distinction between casualty occurrence (p) and its consequence (s) can be used to incorporate research outcomes of the last decades in the field of ship safety. The present work gives a detailed overview of the enhancements provided within the FLARE project concerning the definition of p-factors, including relevant aspects of ship safety neglected or ignored by the current SOLAS regulation.

2. p-Factor Definition

Whilst the s-factor relates to the probability of a ship surviving a given damage (breach) in each loading condition and environment [4], the p-factor is used to define the probability of occurrence of a certain breach in each one of the pertinent hazards (collision, side and bottom grounding) conditional on its dimensions and location (centre and port or starboard side). This entails the need for probabilistic information pertaining to each of these elements, which is provided by the marginal distributions of the breach dimensions and location. Breaches are defined as 3-dimensional objects (location, side, and vertical position along the ship length). Deriving from this, damage breaches are often thought of and described as cuboids; however, this is not always the case. In areas where there is a curvature in the vessel waterline, i.e., outside of the parallel mid body, the damage breach ceases to be described as a cuboid. Instead, the penetration element of the damage breach follows the profile of the waterline corresponding to the draught being examined, offset by the penetration Ly. The p-factor is unaffected by this assumption as the dimensional properties of the damage remain the same. Instead, the geometrical properties of the breach are changed, see Figure 1. However, the spaces affected by the damage breach can vary.
All the clutter in the literature relates to how the p-factor is addressed in current SOLAS and how damage breaches are defined, concerning several pertinent characteristics:
  • SOLAS approach to defining and using the p-factor
“p-factor is the probability that, for a given loading condition, only the compartment or group of compartments under consideration are flooded, disregarding any horizontal subdivision”, Equation (1). SOLAS is still referring to “compartments” for collision damages only, using statistical data of breaches that relate to the last 3 decades of the previous century, the majority of which relate to cargo vessels, consequently disregarding all related information of modern passenger ships, ship size, speed and structural arrangements. Hence, material and speed for the vessel is under consideration. Moreover, the statistical database available in SOLAS includes collisions and contacts as part of the hazards. If there is no breach, then there is no p-factor as its definition is conditional on having a breach of given dimensions, location, and position.
  • Whether the distributions are marginal or conditional probabilities
Even though marginal distributions are supposed to be independent, attention should be paid to the damage penetration. The SOLAS framework implicitly assumes that for a collision damage breach, the ratio between dimensionless penetration and dimensionless length cannot exceed 15. Therefore, an upper limit should be introduced, having as the main consequence that damage length should be generated before damage penetration. Specific reference to this is made in the explanations provided for Figure 4 in Section 2.1.
  • Derivations of the breach distributions based on statistical or direct approaches
Crash analysis using verified numerical Finite Element codes, e.g., LS-Dyna or faster super-element codes, e.g., SHARP, as expanded in Section 4, are widely available, offering potential to address collision and grounding hazards for a specific ship in specified operational scenarios and environmental conditions. Yet, the profession continues to rely on statistical methods, using incomplete or, in the case of passenger ships, irrelevant statistical data, pertaining to cargo ships, for the definition of damage breaches.
  • Zonal or non-zonal approaches and definition of breaches in each approach
The reference of SOLAS to compartments, i.e., physical boundaries to be used in the integration of the probability distributions of breaches to derive the p-factors, is still creating problems between the traditionalists and modern naval architecture. The former believe that the p-factor should be calculated with the help of the law of total probability, resulting from Kolmogorov axioms, as it is in the SOLAS Convention. Using Monte Carlo (MC) sampling of the damage breach distributions is unable to calculate the true value of the A-index. As such, it is of no value for Naval Architects. This is the alienated view being referred to in the paper title. Notwithstanding the above, there are varying views on how to use the non-zonal approach, with confusion being the standard situation.
  • Sampling methods for numerical simulations/calculations of ship survivability
Even though there is some general guidance based on sampling error, there is no rigorous approach to define sample size for use in simulations/calculations, such numbers varying from 1000 samples to 100,000, based on how closely the breach distributions are represented, but without any reference to the reliability of data or the impact on the damage stability calculation in using different sample sizes.
  • Crashworthiness considerations
The question of using crashworthy ship structures to positively affect (reduce) damage breach distributions is another element where confusion prevails, in terms of what exactly this is, how it can be used to improve damage stability, how it is calculated and how it is applied optimally. In particular, the fact that the probability term implicit in the p-factor does not change, only the condition pertaining to the damage size in a given location in the ship. As a result, even though the concept has been around for decades it has not found any real application in ship design in so far as damage stability is concerned.
Each of these aspects will be further addressed in the following sections.

2.1. SOLAS Damage Breach Distributions

The derivation of p-factors, currently in use in SOLAS, originates from the HARDER project [5,6,7,8], during which collision damage statistics were processed to obtain probabilistic damage breach distributions, in terms of damage longitudinal position, longitudinal extent, transversal extent, the upper limit of vertical extent and side of damage (port/starboard). The mathematical integration of these distributions over box-shaped domains allows expressing the p-factors in the known analytical format of SOLAS on ship subdivision. The SOLAS underlying damage distributions have been obtained by pooling collision accidents of all types of ships available at the time, spanning the last 3 decades of the previous century. Moreover, the damage distributions do not explicitly consider the structural design or crashworthiness of the ship. Practically, this implies that even if a ship is designed with a high crashworthiness level, no gain is to be expected in terms of safety in the framework of the current regulations. A second consequence is that SOLAS damage distributions embody an ‘average’ crashworthiness level of the historically damaged ships, which is not necessarily representative of a specific type of ship, or applicable to any type of ship, in particular passenger ships, especially the modern giants populating the current fleet.
More specifically, it is acknowledged that the collision statistics include the main accidents involving cargo ships and tankers, Figure 2.
Geometrically, collision-type damage is idealised in SOLAS as a box with two faces parallel to the waterplane, two faces parallel to the ship transversal plane, and two faces following the hull longitudinal shape at the waterline. Furthermore, the damage box crosses the waterline as well as one side of the ship. In the general case, the damage is modelled using the 6 geometrical parameters ( X c , L x , L y , z U L , z L L , damage side), illustrated in Figure 3.
From a probabilistic point of view, the SOLAS underlying damage breach distributions associated with each potential damage parameter are exemplified in Figure 4. This figure provides an overview of the geometrical model of a collision breach, together with the independent marginal cumulative distributions of the breach characteristics in dimensional form. Starboard and portside damages are equiprobable. The damage is defined as potential, meaning that it could also extend outside the vessel limits. This aspect requires particular attention concerning the positioning of the damage at the ship extremities, keeping consistency with the analytical formulation of zonal p-factors. In case the potential damage is fully contained within the ship length L s , X c corresponds to the damage centre. If the damage partially extends outside the vessel, then the location of X c should be changed, as described in [7]. Even though marginal distributions are supposed to be independent, attention should be paid to the damage penetration, Ly. The SOLAS framework implicitly assumes that for a collision damage breach, the ratio between dimensionless penetration and dimensionless length cannot exceed 15. Therefore, an upper limit L y , max = 15 B L x / L c (where B =ship beam) should be introduced, having as a main consequence that damage length should be generated before damage penetration. As a last remark, the internal limit of the damage follows the waterline at z T (where T =ship draft) shifted by L y , then the collision damage is not always box-shaped.

2.2. Deriving p-Factors Using Zonal and Non-Zonal Damage Breach Distributions

In the zonal approach to probabilistic damage stability, currently adopted by IMO, collision damage cases are defined as three-dimensional cuboids, as outlined in Section 2.1. These are determined following discretisation of the vessel subdivision into zones, which can be conducted either in line with physical subdivision boundaries or “virtual” boundaries. Damage probabilities (p-factors) are then derived for each of these zonal damages and combinations thereof, using damage statistics in the form of marginal distributions, as provided in Figure 5. Damage breach p-factors are then generated by integrating the joint probability function of the non-dimensional damage location and non-dimensional damage length f ( x ¯ , y ¯ ) with respect to each damage zone and combination of zones. The resultant probability then accounts for the occurrence of all damage cases that would fall within the range of either a single zone or a combination of zones.
As the non-dimensional damage location and non-dimensional damage length are considered independent parameters, their joint probability density function can be expressed, as shown in Equation (2):
f ( x ¯ , y ¯ ) = a ( x ¯ ) b ( y ¯ )
The respective p-factor for a given damaged zone or combination of zones can then be calculated through the integration of the underlying probability functions for length and location as follows:
p ( x ¯ 1 , x ¯ 2 ) =   b ( y ¯ )   a ( x ¯ )   d x ¯ d y ¯
Once the probability of damaging a given zone or combination of zones is known, the zonal approach then requires two additional reduction factors r and v to account for the probability of differing degrees of damage penetration and height, respectively. The purpose of these factors is to weigh the damage probability in a manner reflective of the underlying damage distributions. Therefore, the final p-factor for a given zonal damage described by location, length, penetration, and height is given by Equation (4).
p i = p ( x ¯ 1 , x ¯ 2 ) · r ( x ¯ 1 , x ¯ 2 , b i ) · v ( H d )
In contrast, the non-zonal approach works by sampling pertinent damage probability distributions to produce a multitude of damage breaches, characterised by size and location. For this purpose, Monte Carlo sampling is generally adopted to create a damage sample from the marginal damage distributions previously described. The process utilises inverse transform sampling, which involves inverting the cumulative distribution function (CDF) of a given random variable, say F X , to produce F X 1 . Random numbers, u , are then generated from a uniform distribution in the interval [0, 1] and are inputted into the inverse CDF to solve for variable x , see Figure 6. This creates a sample population of the random variable being addressed that is representative of the underlying distribution.
By applying the above process to each of the damage distributions, shown in Figure 4, damage breaches can be constructed by combining the output from each sampled distribution. For a given vessel, each of these breaches will lead to a certain combination of rooms having been compromised. Inevitably, a number of damage breaches will lead to the same rooms being affected, thus creating a smaller subset of distinct damage cases for use in the Attained Index calculation. The frequency of repeated cases is then used to determine the damage probability (p-factor), corresponding to n/N, where n is the number of breaches damaging the same compartment (thus referring to a damage case) and N is the total number of breaches generated (sample size). Figure 7 below provides an illustrative example of the difference between zonal and non-zonal approaches. Here, on the left, the traditional zonal approach can be observed, where the probability of damaging the single zone highlighted is determined by integrating the joint probability function of damage location and damage length with respect to the extremities of the zone and its location. In contrast, on the right-hand side the non-zonal approach is demonstrated, where individual sampled damage scenarios appear as unique points all affecting the same single-zone domain. Here, the damage probability for this one zone damage would be determined as n/N which, if a damage sample size of 1000 scenarios is assumed, would lead to a p-factor of 6/1000 = 0.006. If a sufficient damage sample is taken, the zonal and non-zonal p-factors will converge to the same value.
The determination of the p-factor is detailed in the HARDER project, refs. [5,8,9,10], and amended in the SOLAS2009 probabilistic framework [11], catering for collision hazards only. However, collisions are not the only possible hazard constituting the flooding risk for a ship, especially for passenger ships. For the latter, a lack of due consideration for grounding (side and bottom) hazards at the IMO level over the past few decades, catering for these through deterministic requirements, has shifted the flooding risk focus with side and bottom groundings constituting now the majority of the flooding risk for passenger ships. Figure 8 is indicative of the current situation with flooding hazards for passenger ships.
Whilst SOLAS is still very relevant and the requirements for collision are still valid, there is a growing need for revision by adopting a more holistic regulatory framework accounting suitably for all pertinent hazards. Figure 8 from Project FLARE is indicative of the current situation with flooding hazards for passenger ships. The alarm for this state of affairs has been raised in the past [13,14].
Notwithstanding this, research on the topic of grounding hazards has been persistent and systematic, with significant contributions at the IMO level from Projects SAFEDOR [15], GOALDS [6], EMSA III [16] and eSAFE [17], but IMO rejected to include this in pertinent regulations. Specific developments include an accident database addressing all hazards [18] and leading to new damage breach distributions [19]. Directly related to p-factors determination in probabilistic damage stability calculations, a non-zonal approach for breach generation has been developed, e.g., [20,21] as well as calculations of all pertinent indices and their combination, based on the current IMO framework and accounting consistently for all hazards [22,23].

2.3. Sampling Breach Distributions for Damage Stability Assessment

When generating scenarios by sampling probability distributions, it is important to ensure that the sample is a fair and accurate representation of the underlying distributions. The magnitude of the error incurred here is predominantly a function of the sample size and, as such, it is of great importance to ensure that a statistically valid sample is considered. However, as the sample size increases, so does the calculation time and computational cost, so one must seek to strike a balance between these two competing objectives. In order to make this determination, two approaches could be pursued; one using a commonly adopted engineering approach, and a more rigorous mathematical approach in the sampling process. In this respect, there are two points to consider. The first relates to how accurately the sampling process represents a given distribution, in which case a rigorous approach will produce better results. The second relates to the uncertainty in the determination of the statistical distributions being derived from limited accident data (the marginal distributions expanded upon in Section 2.1). In the latter case, given the fact that accident data is limited and unlikable, unwarranted accuracy in the sampling process might not bear improved accuracy in the calculation of damage stability and survivability (Attained Index of Subdivision from static calculations and Survivability Index from time-domain simulations), using suitable numerical models [24].
Adopting a pragmatic approach, the Standard Error (SE) of the mean may be used to ascertain sample quality and is a measure of the accuracy in which the sample mean x ¯ reflects the actual distribution mean μ, calculated in accordance with Equation (5). Two examples of distributions are considered here. One addressing a single parameter (SOLAS distribution of sea states, characterised by Hs) and the second, flooding risk aggregation, in this case represented by the time it takes a specific ship to capsize in pertinent critical flooding scenarios, and the CDF of time to capsize (TTC). By assessing the magnitude of the Standard Error as a function of sample size (N), the relationship between these two parameters can be derived, as shown in Figure 9, for the first case. Here, it can be observed that there are diminishing returns in error reduction for sample sizes greater than 750 samples. Similar tendencies were identified when assessing other parameters in this way, with a variation ±50 samples found across all cases. This would indicate an optimal sample size of 700–800 samples, in this particular case. However, the sampling process itself, provides a subset of all probable cases with proportional representation of various extents, but fails to capture all possible scenarios. This is particularly true in the case of low probability events (the rail-end of such distributions), which are often poorly represented within small samples. To provide an example, if one were to compare a random damage sample to zonal damages, the ratio of 2-compartment to 4-compartment damages would most likely be the same in each case, however, the sample would only consider a fraction of all probable 2 and 4-compartment cases. As such, by increasing sample size, a greater number of these “black swan” events would be captured, even though the error may remain for the most part unchanged.
S E = σ n
where,
  • σ = sample standard deviation,
  • n = number of samples.
Figure 9. Standard Error (SE) relative to sample size (N) for Hs.
Figure 9. Standard Error (SE) relative to sample size (N) for Hs.
Jmse 10 00348 g009
In addition to considering the Standard Error, confidence intervals are normally derived for each sample in order to illustrate the range of confidence across the sample CDF. For this purpose, the Dvoretzky–Kiefer–Wolfowitz inequality [25], is being utilised, which allows different rates in violation to be identified across the range of the distribution, see Equations (6) and (7). An example of how this error varies relative to sample size is also provided in Figure 10.
F   n ( x ) ε F ( x )   F n ( x ) + ε
ε = l n 2 α 2 n
where,
  • F(x) = the true sample CDF,
  • F n ( x ) = lower and upper bounds,
  • 1 − α = level of confidence, i.e., α = 0.05 for 95% confidence.
Figure 10. Confidence intervals: (a) 100 samples; (b) 1500 samples.
Figure 10. Confidence intervals: (a) 100 samples; (b) 1500 samples.
Jmse 10 00348 g010
Considering the sampling process from a more mathematical perspective, studies and applications in computer science suggest that Latin Hypercube (LH), Quasi-Monte Carlo (QMC) or Randomised Quasi-Monte Carlo (RQMC) methods ensure a faster convergence rate than the traditional Monte Carlo approach when addressing complex functions [26]. Considering this in the particular case of application to damaged ship stability/survivability, a preliminary study, limited to Cruise RoPax bottom groundings, has been carried out for the non-zonal approach, implementing an RQMC sampling method on a reference barge [27]. Traditionally, the application of Monte Carlo sampling of pertinent distributions in assessing ship survivability is well documented [24]. However, such a method introduces randomness in the process, leading to a dispersion of the attained survivability index within multiple sets of generated damages. To this end, recent work in [27] investigates sampling methods alternative to Monte Carlo, based on Latin Hypercube and Randomised Quasi-Monte Carlo processes. The sampling methods application for collisions, side and bottom groundings on a reference barge, available in the literature for benchmark purposes, shows that the Randomised Quasi-Monte Carlo method based on multidimensional Sobol sequences grants a lower dispersion of the final survivability index data within samples of an equivalent size. The application on a sample Cruise ship of Monte Carlo and Randomised Quasi-Monte Carlo methods highlights the possibility of reducing the number of damages necessary to evaluate the survivability index within an engineering consistent confidence interval. The sampling process of damages within the SOLAS probabilistic framework has been analysed, proposing three alternative sampling processes useful to reduce uncertainties and A-index variability, whilst adopting a non-zonal approach. More specifically, the performance of LH and RQMC sampling with a standard MC approach is addressed. The test case for collision, bottom, and side grounding damages on a simple reference barge highlights how the RQMC method, based on multi-dimensional Sobol sequences (SMPL-3), gives more benefits than other procedures in the reduction of variability for partial and total A-indices calculations. A detailed analysis of the evaluated p-factors highlights that the reduction of variability in the A-Index is strongly related to the reduction of the p-values evaluated per each unique damage case among multiple repetitions. Moreover, the SMPL-3 method is capable of detecting a higher number of unique damage cases compared to other methods. Therefore, it could significantly reduce the number of samples to be generated to achieve a target confidence level on the results. The benefits provided by SMPL-3 have been further highlighted by testing the sampling process on a complex internal layout, more granular than traditional geometries used for static calculations. Comparing results with traditional MC sampling, it has been found that the SMPL-3 method grants the same Confidence Interval (CI) on the final A-index, using approximately 1/3 of the total breach samples. However, to clearly identify a suitable lower limit for the sample size needed for damage stability assessment, a more extensive study on a wider number of ships with different sizes is needed. Nevertheless, the results on the reference barge and on the sample cruise ship indicate that the adoption of the SMPL-3 method could be very effective with different internal layouts and sizes. The same procedure can be extended also for dynamic analysis, where the benefits in terms of calculation reduction could be even higher than for static calculations.

2.4. Structural Crashworthiness

2.4.1. General Considerations

Structural design has traditionally been exploited as a means of managing safety related to accidental loads and breaches of hulls. In the 20th century, nuclear-powered ships faced a clear danger if the reactors were to be physically damaged, e.g., by a ship-to-ship collision. This led to Woisin [28] describing some reconfiguration of the hull that would result in a higher tolerance in the collision energy of the side structures prior to undergoing breaching. These first investigations served the purpose of not only creating more crashworthy side structure designs, but also in capturing the mechanics of ship-to-ship collisions. From that period, the work of Minorsky [29] should be noted, which established the proportional relationship between the capacity to absorb collision energy and the volume of the structure involved in deformation. McDermott et al. [30] showed that the key element for ship structures to have an extended capacity to absorb energy is to allow the structure to undergo large membrane tension. Based on his conclusion, substantial work followed with Pedersen and Zhang [31], attempting to estimate collision energy and loads based on the Minorsky empirical formula, while Amdahl [32], Lützen [8], Wierzbicki and Abramowicz [33], and Kitamura [34,35], developed analytical methods using an upper-bound theorem, referred to super-element solutions, the latter addressing both collisions and groundings. Deriving from these findings, a series of novel designs of both side and bottom structures have been and are still being investigated by Lehmann and Peschmann [36], Ludolphy and Boon [37], Graaf et al. [38], Naar et al. [39], Klanac [40], and Klanac et al. [41]. What all these studies have in common is that their conceptual developments are focused on the definition of the topology of a novel crashworthy structure, such as shown here in Figure 11.
Based on these estimation methodologies, many studies have been conducted focusing on protecting certain regions of interest against external forces, such as offshore structures in Storheim and Amdahl [43], Mujeeb-Ahmed et al. [44] and LNG tanker in Wang and H. C. Yu [45]. More recently, Paik [46,47] and Wilson [48], proposed advanced techniques for finite element modelling to simulate structural crashworthiness with increased accuracy in collisions and groundings. Most of these studies conclude that the crashworthiness of ships can be controlled effectively with conventional double-bottom and double-sided structures. Concerning the latter, a detailed methodological approach has been presented in [49], with application on a Cruise ship operating in the Finland archipelago, which is further elaborated in Section 4.

2.4.2. Impact of Crashworthiness on p-Factors

As mentioned earlier, the damage probability distributions utilised within SOLAS are based on accident statistics without taking explicitly into account the structural design or crashworthiness of the ship. This implies that even if a ship is designed with a high crashworthiness level, no gain is to be expected in terms of safety in the framework of the current regulations. In principle, SOLAS damage distributions embody an ‘average’ crashworthiness level of the historically damaged ships, which is not representative of a specific type of ship, for example, modern passenger ships. However, in the same way, these distributions can also be formed on the basis of the crash analysis conducted on an area within the vessel having structural protection. This would yield local damage distributions (p-factors) to be used instead of the standard SOLAS assumptions in case of damages involving such protected spaces. The impact of this consideration is demonstrated heuristically in Figure 12 and expanded further in Section 4.

3. Statistical Approach to p-Factors Determination

Probabilistic distributions of ship collision and grounding breaches are an essential part in the design of crashworthy ships. As indicated earlier, current SOLAS damage distributions for collision are developed based on all ship types. However, large differences in ship design, operation, and safety regulations, may render such assumptions invalid. Moreover, the number of accidents leading to the flooding of large passenger ships are rare, which poses a statistical challenge to obtain desired and accurate distributions. Over the years, there has been a continuous research effort towards the development of damage breach distributions, through various EU-funded projects such as HARDER [5], SAFEDOR [15], GOALDS [6], EMSA III [16] and eSAFE [17]. However, what is currently adopted by SOLAS regulations still pertains only to the earliest of these projects, namely Project HARDER. To address this gap, a concerted effort in the EU Project FLARE [12], focused on devising new damage breach distributions, specifically for large passenger ships, whilst addressing both collision and grounding accidents. To this end, use is made of a newly developed accident database undertaken within this project, leading to the development of pertinent damage distributions for damage length, height, penetration, and damage location.

3.1. Overview of the FLARE Accident Database

This section provides a brief discussion of the flooding database developed in FLARE [18,19]. Figure 13 illustrates the distribution of flooding cases for different types of accidents, spanning the period 1999–2020 for Cruise and RoPax ships, extracted from the IHS Sea-web. The record shows that the number of hull/machinery damages and groundings dominate, followed by collisions. This study focuses mainly on ship flooding due to the initiating events developing external to the ship, namely collision and grounding, disregarding contact where only a few flooding events (8) are registered.

3.1.1. Data Filters

Focusing on the scope of the database development, the following filters are employed to extract the casualty and fleet at risk data:
  • Accident period: 1 January 1999 to 31 October 2020
  • Accident type: collision and grounding
  • Ship size: GT ≥ 3500
  • Ship length (overall): ≥80 m
  • Ship type: Cruise, RoPax, Pure passenger, and RoPax (Rail)
  • Location: worldwide
  • Class type: IACS and non-IACS (for the fleet at risk)
Keeping in mind the 1995 SOLAS Conference and scope of the FLARE project, worldwide accidents during the last 20 years have been investigated. The identification of different accidents into collision and grounding is in line with the definition of accident types mentioned in IMO MSC/Circ. 953, i.e., Collision: striking or being struck by another ship (regardless of whether underway, anchored, or moored); stranding (or grounding): being aground, or hitting/touching shore or sea bottom or underwater objects (wrecks, etc.). To filter large passenger ships from the database, a lower threshold value of 3500 GT is selected, representing an average value based on a simple comparison of Cruise and RoPax ships having an overall length of 100 m. It is, essentially, a compromise between having enough data in the database for meaningful statistical analysis while focusing on large passenger ships. For the same reason, the filter for the ship-built year in the accident period has not been applied in this study.

3.1.2. Data Sources

Figure 14 summarises various sources from which the data is collected. The FLARE database is built mainly on five sources, supplemented by data from ship operators and other public sources, namely:
  • Sea-web (by IHS Markit), IHS Sea-web [50].
  • IMO GISIS (Global Integrated Shipping Information System) [51].
  • EMSA EMCIP (European Marine Casualty Information Platform) [52].
Initially, all the collision and grounding accidents were thoroughly examined, based on the different accident categories defined in the Sea-web, whilst cross-checking accident data with available accident reports and other online sources.

3.1.3. Data Taxonomy

A well-structured taxonomy has been defined to ensure the data is captured and organised in a meaningful manner. The newly updated taxonomy results are evolved from the Sea-web, EMSA EMCIP, and IMO GISIS databases, with the addition of fields related to the natural light at the time of the accident, more explicit details on the weather conditions, damage component, and location.

3.2. Probabilistic Modelling of Breach Distributions

Based on the developed accident database, a detailed statistical analysis was undertaken to derive breach distributions for pertinent ship types and hazards. Most of the breach information mentioned in the database contains qualitative descriptions, for example, relating to the breach as the hole, gash, tear, crack, above/below the waterline, etc., with no real quantitative measures of the damage opening. Table 1 shows the number of samples provided for the different damage locations in collision and grounding accidents, whilst Table 2 indicates the total number of breach data (quantitative measures) available in the database for collision, side, and bottom grounding. The figures clearly indicate that the recorded number of cases is scarce, especially for damage penetration, where such information has been registered in only one case (bottom grounding).
Figure 15, Figure 16 and Figure 17 show the percentage of accidents in the longitudinal, transverse, and vertical damage positions of the ship hull related to collision accidents. For both RoPax and Cruise ships, a similar trend is observed for all the damage positions. Along the length (longitudinal) position of the ship, the bow of the ship dominates, which includes 42% RoPax and 52.6% Cruise. The majority of the collisions occurred above the waterline (84.6% RoPax and 77.8% Cruise). The collisions at the starboard side (52.9% RoPax and 54.5% Cruise) of the ship marginally dominate the port side.
The statistical characteristics of damage parameters (length and width) are analysed based on the best-fit probability distribution function (PDF). The goodness-of-fit (GOF) method, using the Kolmogorov–Smirnov (K–S) tests combined with probability plots for a 95% confidence interval, is used to verify the selected PDF compatibility. The most well-known PDFs were chosen based on their popularity and relevance. The selected PDF is further confirmed using the lowest test statistics, which are the difference between the data sample and the fitted empirical CDF. Based on the results of the statistical analysis, the PDF and CDF of the damage characteristics were established for collision, bottom grounding, and side grounding. Figure 18 shows the breach probability distributions for damage length and breadth for the 3 hazards (collision, side and bottom groundings). These distributions need to be normalised by accounting for the fleet at risk for each one of the hazards and, of course, as indicated earlier, 6 distributions are needed for each hazard to completely describe the breach distributions. The results presented here are early work in Project FLARE to be completed in due course. The purpose of presenting it here is to demonstrate the methodology that needs to be followed in the statistical approach for damage breach definitions.
Table 3 summarises the details of the types of distributions selected and their parameters, along with the corresponding p-factors.

4. Direct Approach to Modelling Breach Distributions

The first step of the methodology is to run a very large number of scenarios for each hazard in question, namely collision, side grounding and bottom grounding, considering a reference ship. The aim is to simulate a large range of representative breaches, adopting, for example, a design of experiment strategies and using suitable crash analysis software. In the example presented here, the Super-Element software SHARP is utilised and the hazard considered is a collision. In the particular example considered, a collision scenario is defined by the following parameters: (a) striking ship type, (b) striking ship initial surge velocity, (c) struck ship initial surge velocity, (d) impact longitudinal position, (e) collision angle, (f) striking ship draft and (g) struck ship draft. For each of these parameters, a range of values has been defined in order to build a load case matrix capable of inducing a large range of pertinent breaches. One thousand nine hundred and eighty collision scenarios have been defined by considering the combination of parameters presented in Table 4. As indicated earlier, current SOLAS damage distributions for collisions are developed based on all ship types. In this respect, considering that the number of accidents leading to the flooding of large passenger ships are rare, as demonstrated by the data presented in Section 3, this poses a statistical challenge to obtain desired and accurate distributions. This, in turn, makes a direct approach much more attractive, especially considering that the right tools are available for this purpose.
With respect to the definition of the collision scenarios, it is to be noted that:
  • Since SHARP considers the structural description of one half of the ship (collisions are modelled at the port side), the structure of the ship has been considered symmetrical, and hence a unique model is used.
  • In all simulations, the struck ship is supposed to be at rest (no initial surge velocity). This is in accordance with Lützen [9], who observed from the collision accident statistics that the most likely surge velocity of the struck ship would be zero. Furthermore, the struck ship considered for the case study, having very limited draft variability, was assumed to be at design draft.
  • According to the probabilistic damage analysis model, the longitudinal position is independent of all other damage variables. On this basis, only impacts at the midship section are modelled. However, the actual longitudinal position varied so that transverse bulkheads can also be directly hit.
In simulating collision scenarios, a large range of striking ships is considered, as it drives the damage size obtained, but also the relationship between the damage longitudinal, transverse and vertical extents. For the analysis presented here, 11 striking ships of various types and dimensions were modelled. The general characteristics of the striking ships considered represent the world fleet and are shown in Table 5.
For this case study, all the calculations have been carried out considering a reference ship the FLOODSTAND SHIP B Cruise ship [53], the main particulars of which are given in Table 6. The super-element structural description has been modelled for a section that is 100 m long along the ship parallel body and centred on the mid-ship section. All materials have been modelled as rigid-perfectly plastic with S235 mild steel properties (see Table 7). The failure strain—which in SHARP is compared to the averaged tension stress within the super-elements—has been considered equal to 10%. Similar values have been observed by other authors to provide a good fit between super-element predictions and experimental results, such as Zhang [54], Lützen [9], and Buldgen et al. [55]. The SHARP super-element model of the struck ship is shown in Figure 19. Its hydrodynamic properties, as required by MCOL, have been obtained using the BV Hydrostar software [56].
As far as the striking ships are concerned, the bow shape has been modelled in SHARP and the ships have been assumed to be rigid. For the studied ship, this assumption is supported by the finite element analysis (FEA) computations carried out during the benchmark of SHARP using striking ship 8, which showed a good agreement between the FEA and SHARP results.
After a simulation of all collision scenarios and filtering damages not compatible with the SOLAS description (i.e., mainly damages with a lower vertical limit above the waterline), it was examined to which extent potential SOLAS damages can be practically simulated. This is demonstrated in Figure 20, Figure 21, Figure 22, Figure 23 and Figure 24, where the main damage parameters ( L x , L y , z UL , z LL ) are presented by pair plots. Overall, it is deemed that the SOLAS domains are well populated by the simulation results. Some unpopulated areas are discussed below:
  • Figure 20 shows that no damages of a length higher than 50 m are obtained. A potential explanation is that the calculation matrix lacks very severe scenarios. Another explanation would be that for the reference ship considered, the SOLAS damage limit of 60 m cannot be physically reached when considering realistic scenarios.
  • Figure 20 also shows that longitudinal damages higher than 20 m ( L x > 20   m ) with low penetrations ( L y < 2.5   m ) cannot be simulated. This may be due to the fact that no initial surge velocity was considered for the struck ship. It could also come from the underlying SOLAS model, which considers that for such types of damages, the longitudinal and transverse extents are independent.
  • Figure 21 shows that the domain is well-populated due to the large striking ships database. No damages have been simulated, with the damage upper limit slightly above the waterline and the damage lower limit slightly below. The simulation of such damages would typically require that the damage is due to the bulb of the striking ship only and that the combination of striking ship draft and bulb height is adequate.
  • From Figure 22, it can be noted that no longitudinal damage can be simulated with a vertical position just above the waterline. However, this was expected since long damages mainly correspond to the more massive striking ships with a high freeboard.
  • Figure 23 shows that simulated damages with large penetration have a lower vertical limit close to the ship bottom. This was expected given the bow shapes of the striking ships.
In Figure 24, the results from Figure 20 are shown after clustering the data into either the striking ship’s initial velocity or collision angle. It is observed, as expected, that the striking ship’s initial velocity has a significant influence on the damage extent and that the collision angle has a strong impact on the damage length.

5. Concluding Remarks

Despite a late start and slow early development in the subject of probabilistic damage stability, the past three decades have seen remarkable progress in the evolutionary development of this subject. Such progress covers specific elements in the probabilistic damage stability calculation/simulation process as well as the process itself. Focussing on the requisite data for such calculation/simulation, no input is more important than the damage breaches for each related hazard (collision, side grounding, bottom grounding) and associated probabilistic content, the so-called p-factors. Pursuing clarification in such determination, the following areas and concerns have been addressed, leading to specific conclusions and recommendations for further work to improve knowledge in this specific subject.:
  • Clarification on what exactly p-factors are and how they are defined in terms of marginal distributions of six parameters: length, breadth, height, location, side of the ship, upper and lower location.
  • How to sample such distributions in order to ensure sufficient accuracy in the damage sample.
  • Explanation of what constitutes zonal or non-zonal methods in damage breach generation.
  • Derivation of the marginal breach distributions based on statistical methods, describing, and using a new accidents database, specific for passenger ships and addressing all pertinent hazards (collision, side grounding, bottom grounding).
  • Explanation and demonstration of a direct approach to deriving pertinent p-factors, using a passenger ship operating in the Gulf of Finland.

Author Contributions

Conceptualization, D.V.; Data curation, D.V., M.P.M.-A., D.P., F.M. and F.C.; Formal analysis, D.V., M.P.M.-A., D.P., F.M. and F.C.; Investigation, D.V., M.P.M.-A., D.P., F.M. and F.C.; Methodology, D.V., M.P.M.-A., D.P., F.M. and F.C.; Resources, D.V., M.P.M.-A., D.P., F.M. and F.C.; Software, M.P.M.-A., D.P., F.M. and F.C.; Supervision, D.V.; Validation, D.V., M.P.M.-A., D.P., F.M. and F.C.; Visualization, D.V., M.P.M.-A., D.P., F.M. and F.C.; Writing—original draft, D.V., M.P.M.-A., D.P., F.M. and F.C.; Writing—review & editing, D.V., M.P.M.-A., D.P., F.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by European Union project FLooding Accident REsponse (FLARE) number, grant number 814753, under H2020 program.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The views presented in the paper are those of the authors alone.

References

  1. IMO. SOLAS-International Convention for the Safety of Life at Sea; IMO: London, UK, 2009. [Google Scholar]
  2. Wendel, K. Die Wahrscheinlichkeit des Uberstehens von Verletzungen. Schiffstechnik 1960, 7, 47–61. [Google Scholar]
  3. Wendel, K. Subdivision of ships. In Proceedings of the 1968 Diamond Jubilee International Meeting—75th Anniversary, SNAME, New York, NY, USA, 4–9 November 1968; Paper No 12. p. 27. [Google Scholar]
  4. Vassalos, D.; Mujeeb-Ahmed, M.P. Conception and Evolution of the Probabilistic Methods for Damage Stability and Flooding Risk Assessment. J. Mar. Sci. Eng. 2021, 9, 667. [Google Scholar] [CrossRef]
  5. HARDER. Harmonisation of Rules and Design Rationale; EC Contact No. GDRB-CT-1998-00028, Final Technical Report; IMO Publishing: London, UK, 1999–2003.
  6. GOALDS. Goal-Based Damage Stability; Project funded by the European 13th Commission, FP7-DG Research, Grant Agreement 233876; EC: Brussels, Belgium, 2009–2012.
  7. Bulian, G.; Francescutto, A. Probability of Flooding Due to Grounding Damage Using a p-factor Formulation; GOALDS Project; EC: Brussels, Belgium, 2010. [Google Scholar]
  8. Lützen, M. Ship Collision Damage. Ph.D. Thesis, Technical University of Denmark, Lyngby, Denmark, 2001. [Google Scholar]
  9. Lützen, M. Damage Distributions; HARDER Document 2-22-D-2001-01-4; Version:4; EC Contact No. GDRB-CT-1998-00028. 29 July 2002. [Google Scholar]
  10. Pawlowski, M. Subdivision and Damaged Stability of Ships; Euro-MTEC Book Series; Fundacja Promocji Przemyslu: Gdansk, Poland, 2004. [Google Scholar]
  11. IMO. Resolution MSC.216(82)—Amendments to the International Convention for the Safety of Life at Sea 1974; IMO: London, UK, 8 December 2006. [Google Scholar]
  12. FLARE. EU H2020—MG2.2 Flooding Accident Response, 2019–2022. Available online: https://flare-project.eu/ (accessed on 28 January 2022).
  13. Papanikolaou, A.; Spanos, D.; Boulougouris, E.; Eliopoulou, E.; Alissafaki, A. Investigation into the sinking of the Ro-Ro passenger ferry Express Samina. Int. Shipbuild. Prog. 2004, 51, 95–120. [Google Scholar]
  14. Papanikolaou, A.; Eliopoulou, E. On the development of the new harmonized damage stability regulations for dry cargo and passenger ships. Reliab. Eng. Syst. Saf. 2008, 93, 1305–1316. [Google Scholar] [CrossRef]
  15. SAFEDOR. Design, Operation and Regulation for Safety; Integrated Project, FP6_2 Contract TIP4-CT-2005-516278; SAFEDOR: Hamburg, Germany, 2005–2009; Available online: www.safedor.org (accessed on 28 January 2022).
  16. EMSA III. Study to Assess Acceptable and Practicable Risk Levels for Damage Stability of Passenger Ships; EMSA III: Lisbon, Portugal, 2013–2015.
  17. eSAFE. Damage Stability of Cruise Ships; Joint Industry Project; Cruise Ship Safety Forum (CSSF); 2019; Available online: https://www.google.com.hk/url?sa=t&rct=j&q=&esrc=s&source=web&cd=&cad=rja&uact=8&ved=2ahUKEwiXrpbD-9XrAhXHc94KHZ82CC4QFjAAegQIBRAB&url=http%3A%2F%2Fwww.shipstab.org%2Ffiles%2FProceedings%2FISSW%2FISSW_2019_Helsinki_Finland%2FPapers%2F7_0_Luhmann_ISSW2019.pdf&usg=AOvVaw2vaPE7xWjn4JuLJkMPIU9S (accessed on 6 September 2020).
  18. Mujeeb-Ahmed, M.P.; Vassalos, D.; Boulougouris, E. Development of collision and grounding accident database for large passenger ships. In Proceedings of the First International Conference on the Stability and Safety of Ships and Ocean Vehicles (STAB&S 2021), Glasgow, UK, 7–11 June 2021. [Google Scholar]
  19. Mujeeb-Ahmed, M.P.; Vassalos, D.; Boulougouris, E. Probabilistic damage distribution and risk modelling of collision and grounding accidents for large passenger ships. In Proceedings of the First International Conference on the Stability and Safety of Ships and Ocean Vehicles (STAB&S 2021), Glasgow, UK, 7–11 June 2021. [Google Scholar]
  20. Zaraphonitis, G.; Bulian, G.; Lindroth, D.; Hamann, R.; Luhmann, H.; Cardinale, M.; Routi, A.L.; Bertin, R.; Harper, G. Evaluation of Risk from Raking Damages Due to Grounding; Final report; DNVGL Report 2015-0168 Rev. 2, Project EMSA/OP/10/2013; Technical Report; European Maritime Safety Agency: Lisbon, Portugal, 2015. [Google Scholar]
  21. Bulian, G.; Lindroth, D.; Ruponen, P.; Zaraphonitis, G. Probabilistic assessment of damaged ship survivability in case of grounding: Development and testing of a direct non-zonal approach. Ocean Eng. 2016, 120, 331–338. [Google Scholar] [CrossRef]
  22. Zaraphonitis, G.; Bulian, G.; Hamann, R.; Eliopoulou, E.; Cardinale, M.; Luhman, H. eSAFE-D2.2.1-Description of Methodology; Joint Industry Project eSAFE. 2017. [Google Scholar]
  23. Bulian, G.; Cardinale, M.; Dafermos, G.; Lindroth, D.; Ruponen, P.; Zaraphonitis, G. Probabilistic assessment of damaged survivability of passenger ships in case of grounding or contact. Ocean Eng. 2020, 218, 107396. [Google Scholar] [CrossRef]
  24. Vassalos, D.; Paterson, D. Towards Unsinkable Ships. Ocean Eng. 2021, 232, 109096. [Google Scholar] [CrossRef]
  25. Dvoretzky, A.; Kiefer, J.; Wolfowitz, J. Asymptotic Minimax Character of the Sample Distribution Function and of the Classical Multinomial Estimator. Ann. Math. Statist. 1956, 27, 642–669. [Google Scholar] [CrossRef]
  26. Cools, R.; Nuyens, D. Monte Carlo and Quasi-Monte Carlo Methods; Springer: Cham, Switzerland, 2014. [Google Scholar] [CrossRef] [Green Version]
  27. Mauro, F.; Paterson, D.; Michalec, R.; Boulougouris, E.; Vassalos, D. A damage sampling method to reduce A-index standard deviation in the probabilistic assessment of ship survivability using a non-zonal approach. In Proceedings of the First International Conference on the Stability and Safety of Ships and Ocean Vehicles (STAB&S 2021), Glasgow, UK, 7–11 June 2021. [Google Scholar]
  28. Woisin, G. Design against collision. In Proceedings of the International Symposium on Advances in Marine Technology, Trondheim, Norway, 2 June 1979; pp. 309–336. [Google Scholar]
  29. Minorsky, V. An analysis of ship collision with reference to protection of nuclear power ships. J. Ship Res. 1959, 3, 208–214. [Google Scholar]
  30. McDermott, J.F.; Kline, R.G.; Jones, E.L.; Maniar, N.M.; Chiang, W.P. Tanker structural analysis for minor collisions. SNAME Trans. 1974, 82, 382–414. [Google Scholar]
  31. Pedersen, P.T.; Zhang, S. Absorbed energy in ship collisions and grounding: Revising Minorsky’s empirical method. J. Ship Res. 2000, 44, 140–154. [Google Scholar] [CrossRef]
  32. Amdahl, J. Energy Absorption in Ship-Platform Impacts. Ph.D. Thesis, Norweigian Institute of Technology, Trondheim, Norway, 1982. [Google Scholar]
  33. Wierzbicki, T.; Abramowicz, W. On the crushing mechanics of thin-walled structures. J. Appl. Mech. 1983, 50, 727–734. [Google Scholar] [CrossRef]
  34. Kitamura, O. Comparative study on collision resistance of side structure. Mar. Technol. SNAME News 1997, 34, 293–308. [Google Scholar] [CrossRef]
  35. Kitamura, O. FEM approach to simulate of collision and grounding damage. In Proceedings of the Second International Conference on Collision and Grounding of the Ships, Copenhagen, Denmark, 1–3 July 2001. [Google Scholar]
  36. Lehmann, E.; Peschmann, J. Energy absorption by the steel structure of ships in the event of collision. Mar. Struct. 2002, 15, 429–441. [Google Scholar] [CrossRef]
  37. Ludolphy, J.; Boon, B. Collision resistant side shell structures for ships. In Proceedings of the International Maritime Design Conference, Kyongju, Korea, 21–24 May 2000. [Google Scholar]
  38. Graaf, B.V.; Vredeveldt, J.; Broekhuijsen, J. Construction aspects for the Schelde Y-shape crashworthy hull structure. In Proceedings of the Third International Conference on Collision and Grounding of Ships, Izu, Japan, 25–27 October 2004; pp. 229–233. [Google Scholar]
  39. Naar, H.; Kujala, P.; Simonsen, B.; Ludolphy, H. Comparison of the crashworthiness of various bottom and side structures. Mar. Struct. 2002, 15, 443–460. [Google Scholar] [CrossRef]
  40. Klanac, A. Design Methods for Safe Ship Structures; Aalto University: Aalto, Finland, 2011. [Google Scholar]
  41. Klanac, A.; Ehlers, S.; Tabri, K.; Rudan, S.; Broekhuijsen, J. Qualitative design assessment of crashworthy structures. In Proceedings of the International Maritime Association of Mediterranean (IMAM), Lisbon, Portugal, 26–30 September 2005; pp. 461–469. [Google Scholar]
  42. Ehlers, S.; Broekhuijsen, J.; Alsos, H.; Biehl, F.; Tabri, K. Simulating the collision response of ship side structures: A failure criteria benchmark study. Int. Shipbuild. Prog. 2008, 55, 127–144. [Google Scholar]
  43. Storheim, M.; Amdahl, J. Design of offshore structures against accidental ship collisions. Mar. Struct. 2014, 37, 135–172. [Google Scholar] [CrossRef]
  44. Mujeeb-Ahmed, M.P.; Ince, S.T.; Paik, J.K. Computational models for the structural crashworthiness analysis of a fixed-type offshore platform in collisions with an offshore supply vessel. Thin-Walled Struct. 2020, 154, 106868. [Google Scholar] [CrossRef]
  45. Wang, B.; Yu, H.C.; Basu, R. Ship and ice collision modeling and strength evaluation of LNG ship structure. In Proceedings of the 27th International Conference on Offshore Mechanics and Arctic Engineering, ASME, Estoril, Portugal, 15–20 June 2008. [Google Scholar]
  46. Paik, J.K. Practical techniques for finite element modeling to simulate structural crashworthiness in ship collisions and grounding Ships and Offshore Structures. Ships Offshore Struct. 2007, 2, 69–80. [Google Scholar] [CrossRef]
  47. Paik, J.K. Advanced Structural Safety Studies with Extreme Conditions and Accidents; Springer: Singapore, 2020. [Google Scholar] [CrossRef]
  48. Wilson, P. Basic Naval Architecture; Springer: Berlin/Heidelberg, Germany, 2018. [Google Scholar] [CrossRef]
  49. Conti, F.; Le Sourne, H.; Vassalos, D.; Kujala, P.; Lindroth, D.; Kim, S.J.; Hirdaris, S. A comparative method for scaling SOLAS collision damage distributions based on ship crashworthiness—Application to probabilistic damage stability analysis of a passenger ship. Ships Offshore Struct. 2021. [Google Scholar] [CrossRef]
  50. Sea-Web. IHS Markit. Available online: https://ihsmarkit.com/products/sea-web-maritime-reference.html (accessed on 1 January 2021).
  51. IMO GISIS. Available online: https://webaccounts.imo.org/Common/weblogin.aspx?App=GISISPublic&ReturnUrl=https%3a%2f%2fgisis.imo.org%2fPublic%2fMCI%2fDefault.aspx&error_message=login_required (accessed on 1 January 2021).
  52. EMSA EMCIP. Available online: https://portal.emsa.europa.eu/emcip-public/#/dashboard (accessed on 1 January 2021).
  53. Luhmann, H. FLOODSTAND WP1: Concept Ship Design B; MEYER WERFT GmbH: Papenburg, Germany, 2009. [Google Scholar]
  54. Zhang, S. The Mechanics of Ship Collision. Ph.D. Thesis, Technical University of Denmark, Limby, Denmark, 1999. [Google Scholar]
  55. Buldgen, L.; Le Sourne, H.; Besnard, N.; Rigo, P. Extension of the super-element method to the analysis of the oblique collision between two ships. Mar. Struct. 2012, 29, 22–57. [Google Scholar] [CrossRef] [Green Version]
  56. Bureau Veritas. Hydrostar for Experts User Manual; Bureau Veritas: Paris, France, 2019. [Google Scholar]
Figure 1. Breach definition along the ship side.
Figure 1. Breach definition along the ship side.
Jmse 10 00348 g001
Figure 2. Ship-type breakdown in collision statistics, Project GOALDS [6].
Figure 2. Ship-type breakdown in collision statistics, Project GOALDS [6].
Jmse 10 00348 g002
Figure 3. Geometric properties of a damage breach.
Figure 3. Geometric properties of a damage breach.
Jmse 10 00348 g003
Figure 4. (a) Damage centre longitudinal position cumulative distribution function; (b) Damage longitudinal extent cumulative distribution function; (c) Damage transversal extent conditional cumulative distribution function; (d) Damage vertical position upper limit cumulative distribution function; (e) Damage vertical position lower limit cumulative distribution function; (f) Damage side index probability mass function.
Figure 4. (a) Damage centre longitudinal position cumulative distribution function; (b) Damage longitudinal extent cumulative distribution function; (c) Damage transversal extent conditional cumulative distribution function; (d) Damage vertical position upper limit cumulative distribution function; (e) Damage vertical position lower limit cumulative distribution function; (f) Damage side index probability mass function.
Jmse 10 00348 g004
Figure 5. 1-Zone damage example using the zonal method.
Figure 5. 1-Zone damage example using the zonal method.
Jmse 10 00348 g005
Figure 6. Inverse transform sampling.
Figure 6. Inverse transform sampling.
Jmse 10 00348 g006
Figure 7. Zonal vs non-zonal 1-compartment damage.
Figure 7. Zonal vs non-zonal 1-compartment damage.
Jmse 10 00348 g007
Figure 8. Recent statistics on the flooding risk of passenger ships, Project FLARE [12].
Figure 8. Recent statistics on the flooding risk of passenger ships, Project FLARE [12].
Jmse 10 00348 g008
Figure 11. Concepts of crashworthy structures: (a) Longitudinal structure on-board an inland waterway gas carrier [37]; (b) Transverse structure on board a RoPax vessel, Ehlers et al. [42]; (c) Corrugated structure on board an inland waterway [42].
Figure 11. Concepts of crashworthy structures: (a) Longitudinal structure on-board an inland waterway gas carrier [37]; (b) Transverse structure on board a RoPax vessel, Ehlers et al. [42]; (c) Corrugated structure on board an inland waterway [42].
Jmse 10 00348 g011
Figure 12. Impact of a crashworthy ship’s structural section (blue line) on damage breach distribution for a typical ship structure (black line).
Figure 12. Impact of a crashworthy ship’s structural section (blue line) on damage breach distribution for a typical ship structure (black line).
Jmse 10 00348 g012
Figure 13. Number of flooding cases for different accidents registered over the last 20 years for Cruise and RoPax ships.
Figure 13. Number of flooding cases for different accidents registered over the last 20 years for Cruise and RoPax ships.
Jmse 10 00348 g013
Figure 14. Main information sources of the FLARE database.
Figure 14. Main information sources of the FLARE database.
Jmse 10 00348 g014
Figure 15. Longitudinal distribution of damage breaches.
Figure 15. Longitudinal distribution of damage breaches.
Jmse 10 00348 g015
Figure 16. Transverse distribution of damage breaches.
Figure 16. Transverse distribution of damage breaches.
Jmse 10 00348 g016
Figure 17. Vertical distribution of damage breaches.
Figure 17. Vertical distribution of damage breaches.
Jmse 10 00348 g017
Figure 18. (a) CDF of damage length (L) for collision; (b) CDF of damage width (W) for collision; (c) CDF of damage length (L) for bottom grounding; (d) CDF of damage width (W) for bottom grounding; (e) CDF of damage length (L) for side grounding; (f) CDF of damage width (W) for side grounding.
Figure 18. (a) CDF of damage length (L) for collision; (b) CDF of damage width (W) for collision; (c) CDF of damage length (L) for bottom grounding; (d) CDF of damage width (W) for bottom grounding; (e) CDF of damage length (L) for side grounding; (f) CDF of damage width (W) for side grounding.
Jmse 10 00348 g018aJmse 10 00348 g018b
Figure 19. Struck ship SHARP model.
Figure 19. Struck ship SHARP model.
Jmse 10 00348 g019
Figure 20. Penetration versus damage length for simulated damages.
Figure 20. Penetration versus damage length for simulated damages.
Jmse 10 00348 g020
Figure 21. Damage vertical position upper limit versus lower limit for simulated damages.
Figure 21. Damage vertical position upper limit versus lower limit for simulated damages.
Jmse 10 00348 g021
Figure 22. Damage vertical positions versus damage length for simulated damages: (a) upper limit; (b) lower limit.
Figure 22. Damage vertical positions versus damage length for simulated damages: (a) upper limit; (b) lower limit.
Jmse 10 00348 g022
Figure 23. Damage vertical positions versus penetration for simulated damages: (a) upper limit; (b) lower limit.
Figure 23. Damage vertical positions versus penetration for simulated damages: (a) upper limit; (b) lower limit.
Jmse 10 00348 g023
Figure 24. Penetration versus damage length for simulated damages: (a) data clustered by collision angle; (b) data clustered by striking ship initial velocity.
Figure 24. Penetration versus damage length for simulated damages: (a) data clustered by collision angle; (b) data clustered by striking ship initial velocity.
Jmse 10 00348 g024
Table 1. Number of accident cases providing qualitative measures of hull damage positions.
Table 1. Number of accident cases providing qualitative measures of hull damage positions.
Damage PositionCollisionGrounding
Bow567
Stern48
Port4812
Starboard8446
Above the waterline663
Below the waterline1985
Table 2. Number of accident cases providing quantitative measures of hull breaches.
Table 2. Number of accident cases providing quantitative measures of hull breaches.
Damage ExtentCollisionSide GroundingBottom Grounding
Length (L)321412
Width (W)10105
Penetration (D)001
Table 3. The probability distribution of breach extents for collision, bottom, and side grounding.
Table 3. The probability distribution of breach extents for collision, bottom, and side grounding.
Accident TypeDamage
Characteristics
PDFParameter
CollisionDamage length (L)3-P Log-logistic α = 1.2086
β = 3.64
γ = 0.0042
Damage width (W)3-P Log-logistic α = 1.5891   β = 2.6846   γ = 0.1695
Bottom groundingDamage length (L)2-P Weibull α = 0.5055   β = 13.22
Damage width (W)3-P Weibull α = 0.4146   β = 4.939   γ = 0.008
Side groundingDamage length (L)3-P Log-logistic α = 0.5635   β = 1.219
γ = 0.07
Damage width (W)General extreme value α = 0.9275
β = 0.4160   γ = 0.3089
Note: 3-P and 2-P denotes three-parameter and two-parameter, respectively.
Table 4. Parameters used in collision crash analysis.
Table 4. Parameters used in collision crash analysis.
ParameterUnit Values
Striking ship type11 ships (see Table 5)
Striking ship initial surge velocity [m/s]2, 4, 6, 8, 10
Impact longitudinal position [m]95.2, 103.6, 112
Collision angle [degrees]30, 45, 60, 90
Table 5. Striking ship’s general characteristics.
Table 5. Striking ship’s general characteristics.
IDVessel TypeLength Overall [m]Moulded Breadth [m]Draft
[m]
Depth
[m]
Displacement
[tonnes]
1Cargo92.214.04.910.03500
2CSV80.017.66.813.83500
3Chemical carrier110.019.57.610.611,064
4Gas carrrier155.022.706.9217.9516,006
5Cargo145.015.878.0011.1515,415
6RoRo180.030.506.8015.8022,062
7Passenger251.028.806.6019.3529,558
8RoPax221.030.006.9015.3230,114
9Bulk carrier180.030.0010.0015.0050,000
10Container300.048.2012.5024.60119,130
11Tanker274.042.0014.9021.00140,000
Table 6. Reference ship main particulars.
Table 6. Reference ship main particulars.
ParameterValue
LPP [m]216.8
Breath moulded B [m]32.2
Depth D [m]16
Draft T [m]7.2
Displacement [tonnes]33,923
Table 7. Material parameters considered.
Table 7. Material parameters considered.
ParameterValue
Yield strength [MPa]235
Tensile strength [MPa]400
Flow stress [MPa]317.5
Failure strain [-]10%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vassalos, D.; Mujeeb-Ahmed, M.P.; Paterson, D.; Mauro, F.; Conti, F. Probabilistic Damage Stability for Passenger Ships—The p-Factor Illusion and Reality. J. Mar. Sci. Eng. 2022, 10, 348. https://doi.org/10.3390/jmse10030348

AMA Style

Vassalos D, Mujeeb-Ahmed MP, Paterson D, Mauro F, Conti F. Probabilistic Damage Stability for Passenger Ships—The p-Factor Illusion and Reality. Journal of Marine Science and Engineering. 2022; 10(3):348. https://doi.org/10.3390/jmse10030348

Chicago/Turabian Style

Vassalos, Dracos, M. P. Mujeeb-Ahmed, Donald Paterson, Francesco Mauro, and Fabien Conti. 2022. "Probabilistic Damage Stability for Passenger Ships—The p-Factor Illusion and Reality" Journal of Marine Science and Engineering 10, no. 3: 348. https://doi.org/10.3390/jmse10030348

APA Style

Vassalos, D., Mujeeb-Ahmed, M. P., Paterson, D., Mauro, F., & Conti, F. (2022). Probabilistic Damage Stability for Passenger Ships—The p-Factor Illusion and Reality. Journal of Marine Science and Engineering, 10(3), 348. https://doi.org/10.3390/jmse10030348

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop