Multi-Attribute Decision-Making Ship Structural Design

: This study develops a procedure for performing multi-attribute decision-making ship structural design of a multi-purpose ship. The already designed ship is further structurally designed to comply with the requirements of the Classiﬁcation Societies. The ship hull and structural components are veriﬁed against yielding, buckling, and ultimate strength. Based on the ultimate limit state (ULS), the ﬁrst order reliability method (FORM) is employed to analyse the structural risk in reducing the probability of failure. The costs associated with materials, manufacturing, and labour are estimated. The structural risk analysis is performed, accounting for different hazard issues related to loss of ship, loss of cargo, loss of human life, and accidental spill of fuel and oil. The risk-based analysis is used to identify an optimum level of ship structural safety, i


Introduction
The ship structural design consists of an iterative decision-making process involving various aspects, such as the type of service, cargo transported, velocity, etc., to determine the optimal ship and structural configuration.One of the objectives of ship design is to identify an efficient and environmentally friendly design solution.
The International Maritime Organization (IMO) developed rules and guidelines that regulated the ship's design to maintain safety in the maritime sector.In 1969, IMO gave responsibility for applying maritime safety rules and standards to the International Association of Classification Societies (IACS).Due to increased competition in the maritime sector, which already accounts for about 90% of world trade, it is necessary to design and optimise more efficient ships with a reasonable safety level for a lower construction, operation, and dismantling expenditure and not only to comply with the minimum requirements, as defined by IACS, but to create safety and environmentally friendly ships.
The ship design goes through a series of evolutionary stages converging to a single design solution, where the most traditional method is the spiral design [1].Different optimisation techniques have been employed to solve the ship design problem, traditionally being a sequential and iterative process, allowing for developing more competitive design solutions.Nowacki, et al. [2] developed computer-aided design and optimisation algorithms, and Hughes [3] and Hughes, et al. [4] established the essential steps in optimising the ship structures.
Optimisation tools are currently taking a more general approach and becoming more reliable in contrast to what was done in the past, where the optimisation was directed to a single objective, becoming a limiting factor.To solve these limitations, several studies [5][6][7][8][9] developed design and optimisation techniques, which were later [10,11] incorporated into multi-criteria optimisation models that include the structural weight and cost of production.IMO adopted risk-based assessment procedures for ship design by defining the formal safety assessment (FSA) [12][13][14][15] as a method to enhance maritime safety, protecting human life in the sea, maritime environment, and cargo and ship integrity by employing cost-benefit assessments.Concerns are raised because many already built ships cannot adequately meet the new requirements introduced by IMO, which impacted the risk-based design performance [16] and short sea shipping [17].
Many studies have been performed using the holistic ship design approach in [16,18,19].A design risk-based framework covering the ultimate limit state and hazards from accidental spills, cargo loss, ship, and crew members was introduced in [18].
A risk-based framework has been developed in a series of studies [19] for ship and structural design and maintenance planning, where the risk analysed covered the structural failure.For structural failure, the time profiles of performance incorporating the structural degradation for structural integrity assessment were also part of the framework [20].
Several aspects related to the ship's structural critical failure modes, type of cargo to which the structure will be subjected during its operational life, and weight and centre of gravity when developing the structural configuration need to be considered when creating a mathematical model that represents the structure, validates the design criteria for different types of failure modes that may occur, and modifies the structural configuration for proper performance and optimisation levels, thus avoiding unnecessary costs.The modern ship designer considers various solutions across performance metrics using the multiple criteria decision-making (MCDM) approach, which involves the selection of the best alternative from pre-specified alternatives [21][22][23].The multiple objective decision-making (MODM) alternatives can be set with multiple equality and inequality constraints in the decision space [24,25].
Recently, several studies [26] analysed the operational characteristics of several shiptype design solutions.The studies focused on the main dimensions of ships, intact stability, and seakeeping performance.Particular attention was paid to the vertical bow acceleration, deck wetness, hydrodynamic impact, propeller emersion, motion sickness, and waveinduced water resistance, which leads to additional brake horsepower.
The present study develops a procedure for performing a multi-attribute decisionmaking multi-purpose ship design.Ship hull design magnification factors (SHDMF)are introduced to generate alternative ship design solutions for an already designed ship.The introduced ship design magnification factors reflect the ship hull structural capacity but also associated risk related to the loss of ship, cargo, human life in the open sea, sea pollution due to accidental fuel and oil spill, construction and operating cost, energy efficiency, and cargo capacity of the ship.Employing MADM, several alternative design solutions are evaluated.

Ship Design
Ship design assesses the owner's specification requirements (ship type, dead weight, speed, data of a similar ship, etc.), life cycle cost, and shipyard capabilities.Specification requirements consider ship hull descriptors, shaft horsepower, lightship weight, dead weight, cargo capacity, free board, stability, seakeeping, midship section design, and ultimate strength assessment of ship hull subjected to still water and wave-induced loads, as well as the progressive structural failure.Due to the substantial number of items that need to be considered, the Pareto optimisation algorithm [27,28] is employed to conclude the best design choice in defining the design parameters.
The work presented in [29] analysed the prerequisites for the environmental pollution driven by maritime transportations in adopting a new ship design concept.The new developed risk-based ship design minimises the risk related to environmental pollution while considering the life cycle assessment and energy efficiency of the ship propulsion system.The optimal design solution for the main dimensions and operation characteristics of a fleet of multi-purpose ships was based on the energy efficiency design index accounting for shipbuilding, operation cost, and resale costs at the end of the service life, where the acceptable design solution involved the requirements for ship resistance, propulsion and stability, free-board, seakeeping, and manoeuvrability.
In the present study, one of the multi-purpose ships that was a part of the fleet designed and analysed in [29] is further developed, concerning the midship section scantling satisfying the additional service grab loading [30] accounting for the structural risk and operational characteristics using MADM [31].The main dimensions of the ship used for the analysis here are given in Table 1.The main dimensions of the initially designed ship will be kept constant during the study.The ship design solution will be evaluated by introducing ship hull design modification factors.SHDMFs are employed as a risk control option to allow the design to be modified realistically and identify the effect on the ultimate capacity of the ship hull structure for estimating the reliability beta index, life cycle cost assessment, and resulting risk.

Midship Structural Design
The midship section of the multi-purpose ship is designed according to the rules of the Classification Society [30], following the ship structural design principle, loads, motions and accelerations, internal and external pressures, and forces needed for the process of the ship hull structural scantling satisfying the yielding, buckling, and ultimate strength requirements.
The width of plates and the spacing of ordinary stiffeners and girder span is defined, considering that the longitudinal span is two to three times the frame spacing.Normal strength steel is used for building the bottom and close to the neutral axis structural components.The choice of high-tensile steel for the structural elements distant from the neutral axis, in this case, the deck area and the hatch coaming, is due to the significant stresses in those locations (see Figure 1).
The vertical bending moments imposed on the ship hull include the still water, waveinduced bending moments, and respective shear forces estimated as stipulated in [23].The load cases used for the structural design are load cases "a" and "b" and load cases "c" and "d", where load cases "a" and "b" refer to the ship in upright conditions, i.e., at rest or having surge, heave, and pitch motions and load cases "c" and "d" refer to the ship in inclined conditions, i.e., sway, roll, and yaw motions [30].The vertical bending moments imposed on the ship hull include the still water, waveinduced bending moments, and respective shear forces estimated as stipulated in [23].The load cases used for the structural design are load cases "a" and "b" and load cases "c" and "d", where load cases "a" and "b" refer to the ship in upright conditions, i.e., at rest or having surge, heave, and pitch motions and load cases "c" and "d" refer to the ship in inclined conditions, i.e., sway, roll, and yaw motions [30].
The scantling procedure involves decision variables, constraints, and objective functions according to [30].The decision variables involved in the optimisation process include the ship hull plate thickness, , , defined as: the material yield stress, , as a discrete variable, according to the existing shipbuilding steels: The constraints applied are related to the local strength, represented by the minimum plate thickness: the minimum section modulus (longitudinal strength) for the deck and bottom are defined as: and the minimum critical buckling stress is defined as: where is the usage factor, , is the acting stress, and , is the critical buckling stress.The scantling of the midship section defines the net thickness of the plate panel subjected to in-plane normal stresses acting on the shorter side, which needs to be not less than the one as defined in [30] in satisfying the yield criteria: = 14.9 $ $ % &'( ) ( ( *+ , * + ( .+ , ./ 0 1 ≥ , The scantling procedure involves decision variables, constraints, and objective functions according to [30].The decision variables involved in the optimisation process include the ship hull plate thickness, t c,i , defined as: the material yield stress, σ y , as a discrete variable, according to the existing shipbuilding steels: σ y, min ≤ σ y, lower shell ≤ σ y, max The constraints applied are related to the local strength, represented by the minimum plate thickness: the minimum section modulus (longitudinal strength) for the deck and bottom are defined as: and the minimum critical buckling stress is defined as: where η i is the usage factor, σ a,i is the acting stress, and σ c,i is the critical buckling stress.The scantling of the midship section defines the net thickness of the plate panel subjected to in-plane normal stresses acting on the shorter side, which needs to be not less than the one as defined in [30] in satisfying the yield criteria: where: p S is the still water pressure and p w is the wave-induced pressure, s is the shorter side of plating, and l is the longer side of plating, C a is the aspect ratio of the plate panel, C r is the coefficient of curvature R y is the minimum yield stress, and γ R , γ m , γ S2 , γ W2 are utilisation factors.
The minimum net shear sectional area A Sh and the net section modulus W, for ordinary longitudinal stiffener subjected to lateral pressure, are to be not less as defined in [30]: where β S , β b , m are factors defined in [32].The ship's structural elements are also subjected to a compressive load, which may lead to a buckling failure.The type of the load and their combinations employed in the current analysis is axial, bending, and shear loadings.The common causes of plate buckling of ship structural components are due to the high compressive and residual stresses, high shear stresses, combined stresses, lack of flexural rigidity, lack of stiffening, significant initial imperfections, extensive and improper use of high tensile steel, and excessive material wastage due to general and local pitting corrosion [32][33][34][35][36].
The general modes of failure of stiffened panels are identified as the lateral buckling of stiffeners, torsional buckling of stiffeners, flexural buckling of stiffeners, flexural buckling for plate stiffener combination, and buckling of plate panels between stiffeners.The combined critical stress σ comb , for plate panels subjected to compressive axial, bending, and shear loads are defined as [30]: where F, R eH are defined in [30].
The critical buckling stress, σ c for compressive axial and bending loads is defined as [30]: where σ E is the Euler buckling stress.
The critical buckling stress of the ordinary stiffeners is estimated as [30]: where σ b is defined in [30].Bureau Veritas software MARS 2000 [37] is employed to analyse the ship's structural hull girder's ultimate strength.The software adopted the progressive collapse method [31] to analyse the ultimate strength of the hull girder between two adjacent frames [32].In this analysis, the midship section is divided into structure elements: stiffener attaching plating element and hard corner element, acting independently in their failure modes.
The incremental-iterative approach obtains the moment-curvature relationship [38].Each iteration's bending moment acting on the hull girder transverse section increases due to the imposed curvature.Each structural component has an axial strain due to the angle of rotation of the hull girder transverse section about its horizontal neutral axis.In the sagging conditions, the structural elements above the neutral axis are shortened, while those below the neutral axis are lengthened.The location of the neutral axis and the cross-section of the ship are calculated based on the failure mode of each structural element as the external moment is applied.The tensile structural elements represent a single mode of elastic-plastic failure, while in compression, they present the mode of buckling or yielding.
The pink dashed line in Figure 2 shows the applied bending moment the hull girder needs to support [23].As seen in Figure 2, the higher ultimate strength occurs in the hogging condition and lower in sagging loading conditions, respectively.The worst possible scenario for structure failure appears in the sagging loading condition.angle of rotation of the hull girder transverse section about its horizontal neutral axis.In the sagging conditions, the structural elements above the neutral axis are shortened, while those below the neutral axis are lengthened.The location of the neutral axis and the crosssection of the ship are calculated based on the failure mode of each structural element as the external moment is applied.The tensile structural elements represent a single mode of elastic-plastic failure, while in compression, they present the mode of buckling or yielding.The pink dashed line in Figure 2 shows the applied bending moment the hull girder needs to support [23].As seen in Figure 2, the higher ultimate strength occurs in the hogging condition and lower in sagging loading conditions, respectively.The worst possible scenario for structure failure appears in the sagging loading condition.A change in the ultimate strength is most effectively achieved by modifying the thickness of structural components [39,40].By employing SHDMF in the structural analysis, the midship section structure is redesigned, keeping the thickness of the structure closest to the neutral axis constant (inner side and side shell) and changing the thickness of the structure farthest from the neutral axis (bottom, double bottom, and deck structure) where the global failure most probably is expected to occur.

Structural Reliability Analysis
The ship hull structural progressive collapse due to the loss of structural stiffness and strength is related to the loss of equilibrium, attainment of the maximum capacity of the resistance by yielding, rupture or fracture and the instability resulting from the buckling or elastoplastic collapse of plating, stiffened panels, and support members, and it is defined as the ultimate limit state.The limit state function in the reliability assessment of the ship hull structure is based on the ultimate strength of the ship hull, as defined in [41]: where O P is the ultimate bending moment and O P is the still water bending moment fitted to a normal distribution [42].The regression equations define the statistical descriptors of the still water bending moment as a function of the length of the ship, Q and deadweight ratio, 3 = JR3S/Full loadL [43,44] as: A change in the ultimate strength is most effectively achieved by modifying the thickness of structural components [39,40].By employing SHDMF in the structural analysis, the midship section structure is redesigned, keeping the thickness of the structure closest to the neutral axis constant (inner side and side shell) and changing the thickness of the structure farthest from the neutral axis (bottom, double bottom, and deck structure) where the global failure most probably is expected to occur.

Structural Reliability Analysis
The ship hull structural progressive collapse due to the loss of structural stiffness and strength is related to the loss of equilibrium, attainment of the maximum capacity of the resistance by yielding, rupture or fracture and the instability resulting from the buckling or elastoplastic collapse of plating, stiffened panels, and support members, and it is defined as the ultimate limit state.The limit state function in the reliability assessment of the ship hull structure is based on the ultimate strength of the ship hull, as defined in [41]: where M u is the ultimate bending moment and M sw is the still water bending moment fitted to a normal distribution [42].The regression equations define the statistical descriptors of the still water bending moment as a function of the length of the ship, L and dead-weight ratio, W = (DWT/Full load) [43,44] as: where M sw,CS is the still water bending moment as given in [30] and Mean(M sw,max ) and StDev(M sw,max ) are defined as a function of L and W as suggested in [43].
The wave-induced bending moment, M w is fitted to the Gumbel distribution [41], x u is the uncertainty on ultimate strength, x sw is the uncertainty prediction on the still water bending moment, x w is the uncertainties in the wave-induced bending moment due to linear seakeeping analysis, and x s is a factor that accounts for the nonlinearities in sagging load.The uncertainty coefficients x sw , x w , x s are fitted to a normal distribution function, with a mean value of 1.00 and a standard deviation of 0.1.The model uncertainty on ultimate strength x u is fitted to a normal distribution function with a mean value of 1.05 and a standard deviation of 0.1.It is assumed that the MARS 2000 software calculates the ultimate bending moment at the confidence level of 5%, where M 5% u = M C u and it is fitted to a log-normal probability density function, where the covariance, COV, is 0.08 [45].
The reliability index of the midship hull structure for the net and gross designs is estimated at the end of the service life of the ship, considered here as τ s = 25 years, assuming that the structure of the ship is entirely corroded, i.e., there is no corrosion margin as determined by the Classification Society rules, concerning the net ship hull structural design, and the gross structural design is considered when the ship structure is not corroded, and the corrosion protection is appropriately functioning.It is also assumed that the structure of the midship section is subjected to general corrosion degradation, where over the years, its degradation has occurred for all structural components.The mean value and standard deviation of the corrosion depth as a function of time is taken as defined in [46], and the statistical descriptors of the vertical wave-induced bending moment are in [41].
Using the first order reliability method [47,48], failure may occur when g(X) fails, leads to P f = P(g(X) < 0), which may be evaluated by using the standard normal distribution function as P f = P(g(X) < 0) = Φ(−β), where β is the beta reliability index [49].
The time-dependent non-linear beta reliability index, β(t), where t ∈ [0, τ S ] is defined as [19]: where τ c = 6.50 years is the coating life and τ t = 11 years is the transition life, β gross is the beta reliability index related to the gross thickness of structural components and β net is related to the net thickness.The time-dependent non-linear Beta reliability index as a function of SHDMFs can be seen in Figure 3.
J. Mar.Sci.Eng.2022, 10, x FOR PEER REVIEW 7 of 18 where O ,]* is the still water bending moment as given in [30] and OZ[\9O , ; and ^ RZ_9O , ; are defined as a function of Q and 3 as suggested in [43].The wave-induced bending moment, O P is fitted to the Gumbel distribution [41], M N is the uncertainty on ultimate strength, M N is the uncertainty prediction on the still water bending moment, M N is the uncertainties in the wave-induced bending moment due to linear seakeeping analysis, and M N is a factor that accounts for the nonlinearities in sagging load.The uncertainty coefficients M N , M N , M N are fitted to a normal distribution function, with a mean value of 1.00 and a standard deviation of 0.1.The model uncertainty on ultimate strength M N is fitted to a normal distribution function with a mean value of 1.05 and a standard deviation of 0.1.It is assumed that the MARS 2000 software calculates the ultimate bending moment at the confidence level of 5%, where O `% = O ] and it is fitted to a log-normal probability density function, where the covariance, COV, is 0.08 [45].
The reliability index of the midship hull structure for the net and gross designs is estimated at the end of the service life of the ship, considered here as b = 25 years, assuming that the structure of the ship is entirely corroded, i.e., there is no corrosion margin as determined by the Classification Society rules, concerning the net ship hull structural design, and the gross structural design is considered when the ship structure is not corroded, and the corrosion protection is appropriately functioning.It is also assumed that the structure of the midship section is subjected to general corrosion degradation, where over the years, its degradation has occurred for all structural components.The mean value and standard deviation of the corrosion depth as a function of time is taken as defined in [46], and the statistical descriptors of the vertical wave-induced bending moment are in [41].
Using the first order reliability method [47,48], failure may occur when IJKL fails, leads to g h = gJIJKL i 0L, which may be evaluated by using the standard normal distribution function as g h = gJIJKL i 0L = ΦJ−4L, where 4 is the beta reliability index [49].
The time-dependent non-linear beta reliability index, 4J L, where ∈ l0, b * m is defined as [19]: where b = 6.50 years is the coating life and b = 11 years is the transition life, 4 n is the beta reliability index related to the gross thickness of structural components and 4 is related to the net thickness.The time-dependent non-linear Beta reliability index as a function of SHDMFs can be seen in Figure 3.

Cost-Benefit Analysis
The ship is subjected to corrosion degradation and structural failures during the service life due to progressive structural collapse.To control the risk associated with the structural collapse of the ship's hull, accounting for the existing uncertainties, based on an identified failure scenario, the risk is measured as a product of the likelihood of structural failure and its consequences defined as: where P f ,j P g X 1,j t ≤ 0 is the probability of failure, C f ,j (X 2 |t) is the consequence cost of failure, X 1,j and X 2,j are the vectors of parameters involved in the probability of failure and consequence analyses that occur during the service life of the ship τ s .
The objective is to reduce the risk to an acceptable level by optimising the ship's hull structure and evaluating alternative decision-making options.The method used to define the acceptable risk is through the reliability index that minimises the cost of the structure's design, where various failure modes may result in economic, environmental, and human losses and other consequences.
The risk, Risk(t n |SHDMF, β) is defined as [19]: where C P f (t n |SHDMF, β) is the cost associated with the structural failure over the service life of the ship and C me (SHDMF, β) is the cost of implementing a structural safety measure accounting for SHDMF, including the cost of material and labour needed to redesign the ship hull structure.The life cycle cost (LCC) of the ship is defined as the total cost of all distinct phases of life of the ship and its equipment, including design, acquisition, operation, maintenance, upgrade, and dismantling.It is defined as a sum of cost estimates from the beginning to the end of the life cycle and is used in the design process of all engineering systems, including ships [20,29,31] and offshore structures.The total costs of the ship related to the distinct phases of the ship's life cycle are divided into three diverse groups, including capital expenditure (CAPEX), operating expenditure (OPEX), and dismantling expenditure (DECEX).
The MARAD [50] system is used for the group-specific ship systems and their associated costs.The systems groups are based on the different components in constructing the ship's life cycle, and their costs are included in the construction project.
The initial capital cost estimate for constructing multi-purpose ships is based on several design parameters, such as main dimensions, dead-weight, lightship weight, propulsive power, etc.A regression analysis makes it possible to estimate CAPEX using a mathematical relationship between the input parameters and the cost of construction [51].Construction costs are divided into four components: material, labour, overheads, and profits.It is assumed that the steel price is 580 EUR/tonne, and the equipment is 1500 EUR/tonne [52].
The annual operating expenditure is a sum of the salary of crew members, costs related to the stores and supplies, insurance, port expenses, and annual fuel cost.The dismantling cost is evaluated as a function of the lightship weight.
The cost associated with the structural failure over the service life of the ship is estimated as a function of SHDMF, reliability index, and time [19]: where P f t j SHDMF, β is the probability of failure, C S t j SHDMF, β is the total cost of the ship in the year t j ∈ [0,τ s ], C C is the cost associated with the loss of the cargo, C d is the cost associated with the accidental spill, C v is the cost associated with the loss of human life, and γ = 5% is the assumed value of the discount rate (see Figure 4).
where g h 9 v y^|RO=, 4; is the probability of failure, $ * 9 v y^|RO=, 4; is the total cost of the ship in the year v ∈ [0,b * ], $ ] is the cost associated with the loss of the cargo, $ is the cost associated with the accidental spill, $ • is the cost associated with the loss of human life, and ( = 5% is the assumed value of the discount rate (see Figure 4).The cost of the ship, $ * 9 v y^|RO=, 4;, is a function of the ship's age, that is, the initial cost of the ship, ‚ = 0 years, and the scrapping cost = 25 ƒears accounting for corrosion degradation [21] is estimated as [46]: where $ J ‚ |^|RO=, 4L is the initial cost of the ship, $ J |^|RO=, 4L is the scrapping cost, v is the year of the operation and $ = 270 EUR/tonne is the assumed value of the scrap cost.
The cost associated with the loss of cargo, $ ] , is estimated as [39]: where $ n = 1200 EUR/tonne is the assumed cost of a tonne of cargo, … n = 20% is the considered partial factor of the cargo lost, and g n is the total amount of the cargo of the ship, estimated as 7200 tonnes.
The cost of the accident spill, $ is estimated as [39]: where … = 10% is the considered a partial factor of the fuel oil spill, g = 10% is the probability that the fuel oil split reaches the shoreline [36], CATS = 60,000 EUR/tonne [53] is the cost of one tonne of accidentally spilt fuel oil and 3 h is the total weight of fuel oil in tonnes.
The cost associated with the loss of human life, $ • , related to the loss of crew members is estimated as: where \ is the number of crew members, … is the probability of loss of the life of a crew member (considered 25%), and ICAF is the implied cost of averting the fatality.The cost of the ship, C S t j SHDMF, β , is a function of the ship's age, that is, the initial cost of the ship, t 0 = 0 years, and the scrapping cost t n = 25 years accounting for corrosion degradation [21] is estimated as [46]: where C s (t 0 |SHDMF, β) is the initial cost of the ship, C s (t n |SHDMF, β) is the scrapping cost, t j is the year of the operation and C scrap = 270 EUR/tonne is the assumed value of the scrap cost.The cost associated with the loss of cargo, C C , is estimated as [39]: where C cargo = 1200 EUR/tonne is the assumed cost of a tonne of cargo, f cargo = 20% is the considered partial factor of the cargo lost, and P cargo is the total amount of the cargo of the ship, estimated as 7200 tonnes.
The cost of the accident spill, C d is estimated as [39]: where f spill = 10% is the considered a partial factor of the fuel oil spill, P sl = 10% is the probability that the fuel oil split reaches the shoreline [36], CATS = 60,000 EUR/tonne [53] is the cost of one tonne of accidentally spilt fuel oil and W f uel oil is the total weight of fuel oil in tonnes.
The cost associated with the loss of human life, C v , related to the loss of crew members is estimated as: where n crew is the number of crew members, f crew is the probability of loss of the life of a crew member (considered 25%), and ICAF is the implied cost of averting the fatality.The cost associated with the fatality is based on the implied cost of avoiding a fatality ICAF, which uses a risk model obtained from the average of the OCDE countries [54].
The cost of implementing the structural safety measure accounts for the SHDMF, which is also associated with the reliability level, including the cost of material and labour.
The cost of structural safety measures, C me (SHDMF, β) is positive or negative depending on if the value of SHDMF is larger or smaller than 1, respectively [39]: where ∆W steel (SHDMF, β) is the weight of steel, in tonnes, because of the design modifications factor, C steel is the cost of steel, and C labour is the cost of labour.The weight of steel is estimated as [39]: where W steel is the weight of steel related to the ship hull structural design (see Figure 5).The cost associated with the fatality is based on the implied cost of avoiding a fatality ICAF, which uses a risk model obtained from the average of the OCDE countries [54].
The cost of implementing the structural safety measure accounts for the SHDMF, which is also associated with the reliability level, including the cost of material and labour.The cost of structural safety measures, $ J^|RO=, 4L is positive or negative depending on if the value of SHDMF is larger or smaller than 1, respectively [39]: where ∆3 J^|RO=, 4L is the weight of steel, in tonnes, because of the design modifications factor, $ is the cost of steel, and $ is the cost of labour.The weight of steel is estimated as [39]: where 3 is the weight of steel related to the ship hull structural design (see Figure 5).The cost-benefit analysis is used to identify an optimum level of ship safety, i.e., the optimum reliability index, controlling the risk associated with changing the initial design.The cost-benefit analysis of the redesigned structure, according to the ship hull design magnification factors related to the scantlings of the midship section, is conducted based on the expected total risk.
The target reliability level β reflects the structural failure cost associated with risk control since each cost is a function of the reliability index, see Figure 6.According to the study performed in [55], the range of the target reliability index over the ship's service life can be between 2.79 and 5.38.The most economical target reliability index here is 4.29.The cost-benefit analysis is used to identify an optimum level of ship safety, i.e., the optimum reliability index, controlling the risk associated changing the initial design.The cost-benefit analysis of the redesigned structure, according to the ship hull design magnification factors related to the scantlings of the midship section, is conducted based on the expected total risk.
The target reliability level β reflects the structural failure cost associated with risk control since each cost is a function of the reliability index, see Figure 6.According to the study performed in [55], the range of the target reliability index over the ship's service life can be between 2.79 and 5.38.The most economical target reliability index here is 4.29.There are some difficulties in defining the cost-related parameters, which is explained by the fact that it is necessary to account for factors that impact the cost forecast.This becomes a critical issue due to the long in-service time of ships and the evolution of prices in such an extended period.Time series analysis may be used in economic forecasting.To analyse and forecast prices, an average, naive, random walk or decomposition, exponential smoothing, or autoregressive integrated moving average [56] process may be There are some difficulties in defining the cost-related parameters, which is explained by the fact that it is necessary to account for factors that impact the cost forecast.This becomes a critical issue due to the long in-service time of ships and the evolution of prices in such an extended period.Time series analysis may be used in economic forecasting.To analyse and forecast prices, an average, naive, random walk or decomposition, exponential smoothing, or autoregressive integrated moving average [56] process may be used.

Decision-Making Ship Design
The ship hull structural design is a multiple-objective decision problem, which requires identifying an optimal structural design solution to fulfil the established primary descriptors of the ship and satisfy the constraints.The goal here is to achieve fast and accurate design solutions finalised at minimising the total expected cost and risk, which is a function of numerous factors, including the ship hull strength and cargo capacity, lightship weight, and energy efficiency.Any acceptable design solution needs to follow the ship design criteria related to the free-board, stability, seakeeping, manoeuvrability, etc.
In this regard, additionally to the risk of structural collapse, estimated in monetary terms, it is essential to include in the multiple attribute's decision analysis the annual operating cost, annual construction cost, cargo capacity measured by the number of containers transported, and attained energy efficiency design index as a function of SHDMF.The lightship weight is defined as the sum of the weight of the hull structure, equipment, outfitting, and machinery [56,57].
Since the ship weight and buoyancy forces are in balance, any change in the structural design due to the application of the SHDMF will reflect the lightship weight LW.To keep the predefined draft, trim, and the longitudinal centre of gravity location, the dead weight DWT needs to be adjusted in the part of the cargo and ballast content, which will lead to new acting still water bending moment Msw.The estimated LW, DWT and Msw for the generated design solutions as a function of SHDMF are presented in Figure 7.The acting still water bending moments are shown as red circles in Figure 8, and the design still water bending moment, for hogging and sagging, by the green and blue envelopes.As can be seen, the attained still water bending moments as a function of SHDMF at the midship section are much below the design one as stipulated in [30].The difference between the design and acting steel water bending moment will enhance the already established structural reliability index estimated based on the design steel water and waveinduced bending moment and ultimate strength.The acting still water bending moments are shown as red circles in Figure 8, and the design still water bending moment, for hogging and sagging, by the green and blue envelopes.As can be seen, the attained still water bending moments as a function of SHDMF at the midship section are much below the design one as stipulated in [30].The difference between the design and acting steel water bending moment will enhance the already established structural reliability index estimated based on the design steel water and wave-induced bending moment and ultimate strength.
The acting still water bending moments are shown as red circles in Figure 8, and the design still water bending moment, for hogging and sagging, by the green and blue envelopes.As can be seen, the attained still water bending moments as a function of SHDMF at the midship section are much below the design one as stipulated in [30].The difference between the design and acting steel water bending moment will enhance the already established structural reliability index estimated based on the design steel water and waveinduced bending moment and ultimate strength.The operating costs are ongoing expenses connected with the vessel's day-to-day running (excluding fuel, which is included in voyage costs) and an allowance for everyday repairs and maintenance.Operating costs' principal components include manning costs, stores, routine repair and maintenance, insurance, and administration.The annual operating cost, transported containers, and construction costs, as a function of SHDMF, are given in Figure 9.
As seen in Figure 9, the ship construction and operating cost for transporting the same cargo in 20-foot containers are of the same trends, and by increasing SHDMF, they are also growing.However, going into a more detailed analysis of different structural design solutions as a function of SHDMF, it seems that the number of containers the ship may transport is growing with the reduction in the SHDMF.The operating costs are the ongoing expenses connected with the vessel's day-to-day running (excluding fuel, which is included in voyage costs) and an allowance for everyday repairs and maintenance.Operating costs' principal components include manning costs, stores, routine repair and maintenance, insurance, and administration.The annual operating cost, transported containers, and construction costs, as a function of SHDMF, are given in Figure 9.As shown in Figure 10, for all structural design solutions as a function of SHDMF, the relatively more energy-efficient ships indicate a lower level of EEDI for the lower SHDMF.This can be explained by reducing SHDMF, the ship's lightship weight is diminished due to the reduced plate thicknesses.This reduction is compensated by increasing the DWT and the respective cargo capacity by keeping the draft constant.The encountered design solution does not involve enhancing the propulsion system and any engineering solution in reducing EEDI.As seen in Figure 9, the ship construction and operating cost for transporting the same cargo in 20-foot containers are of the same trends, and by increasing SHDMF, they are also growing.However, going into a more detailed analysis of different structural design solutions as a function of SHDMF, it seems that the number of containers the ship may transport is growing with the reduction in the SHDMF.
As shown in Figure 10, for all structural design solutions as a function of SHDMF, the relatively more energy-efficient ships indicate a lower level of EEDI for the lower SHDMF.This can be explained by reducing SHDMF, the ship's lightship weight is diminished due to the reduced plate thicknesses.This reduction is compensated by increasing the DWT and the respective cargo capacity by keeping the draft constant.The encountered design solution does not involve enhancing the propulsion system and any engineering solution in reducing EEDI.
As shown in Figure 10, for all structural design solutions as a function of SHDMF, the relatively more energy-efficient ships indicate a lower level of EEDI for the lower SHDMF.This can be explained by reducing SHDMF, the ship's lightship weight is diminished due to the reduced plate thicknesses.This reduction is compensated by increasing the DWT and the respective cargo capacity by keeping the draft constant.The encountered design solution does not involve enhancing the propulsion system and any engineering solution in reducing EEDI.MCDM identifies the "best" alternative from pre-selected options, considering a set of several objectives [21][22][23], which can be a set with multiple equality and inequality type constraints in the decision space [24,25].
The impact of the design parameters on the ship's design solution is analysed using the multi-attribute decision-making Technique of Order Preference by Similarity to Ideal Solution (TOPSIS) [21].The multi-attribute decision-making aims to identify and choose the best alternatives that fit the design goals.The decision algorithms induce an order on a set of options based on the information [58] about the design descriptions of possible solutions, criteria to measure the performance of solutions, and the preference statements to indicate relative importance.MCDM identifies the "best" alternative from pre-selected options, considering a set of several objectives [21][22][23], which can be a set with multiple equality and inequality type constraints in the decision space [24,25].
The impact of the design parameters on the ship's design solution is analysed using the multi-attribute decision-making Technique of Order Preference by Similarity to Ideal Solution (TOPSIS) [21].The multi-attribute decision-making aims to identify and choose the best alternatives that fit the design goals.The decision algorithms induce an order on a set of options based on the information [58] about the design descriptions of possible solutions, criteria to measure the performance of solutions, and the preference statements to indicate relative importance.
Eleven ship design alternatives as a function of SHDMF are generated and compared.The ideal option, the best level for all attributes considered, and the perfect negative option, which has the worst attribute values, are identified.The TOPSIS method identifies the closest to the ideal solution and farthest from the excellent damaging alternative.
In any ship design solution, i = 1, . . .n is a function of the design criteria, which are scored, x ij concerning the criterion j = 1, . . .m, where a matrix X = x ij of n × m is developed.C 11 to C 2m are defined based on the ship design, scantling, and strength analysis [21].J + is the set of benefit criteria, where the more significant score represents a better condition.J − is the set of negative criteria, where the lower score represents the better condition.
The TOPSIS calculation procedure covers several steps.The first one is the construction of a normalised decision matrix [59]: In the next step, the weighted normalised decision matrix is constructed using a set of weights r ij for each criterion w j : The positive ideal solution is defined by: where v * j = max v ij if j ∈ J + , min v ij if j ∈ J + and the negative ideal solution is defined as: where v j = min v ij if j ∈ J − , max v ij if j ∈ J − and the negative ideal solution is defined as: The separation measures for each alternative are estimated as [59]: The relative closeness to the ideal solution C * i is calculated as: The option C * i , that is closest to 1, is the best-suited solution.In the present study, an inverse multi-attribute decision-making calculation is also performed for any design solution to identify the importance of different criteria, quantified by the maximum relative closeness to the ideal solution.Initially, the weights of all design criteria are assumed to be equal and uniformly distributed.At the end of the calculations, the weights are estimated according to the maximum relative closeness to the ideal design solution and identify the significance of any individual criterion in any specific design solution.
Assuming that the five design criteria have the same weight, 0.2 in the multi-attribute decision-making analysis, the relative closeness to the ideal solution makes the sixth design solution the most acceptable solution at 97.1% significance, followed by the seventh and fifth design solutions (see Figure 11).However, the weight of the design criteria may differ from the assumed uniform distribution here, and a deeper analysis needs to be performed.The inverse calculation is performed to identify the importance of different criteria for the maximum relative closeness to the ideal solution.As can be seen from Figure 12, the best-suited design criterion estimated for the design solution, DS : .Associated with SHDMF = 1 is EEDI = 37%, followed by the cargo capacity, represented by the number of containers transported, 22%, the construction cost of 21% and the operating cost of 20%.The risk is not seen as a significant factor for DS1.However, the DS : design solution has the minimum construction and operating cost and maximum cost associated with the risk.The inverse calculation is performed to identify the importance of different criteria for the maximum relative closeness to the ideal solution.As can be seen from Figure 12, the best-suited design criterion estimated for the design solution, DS 1 .Associated with SHDMF = 1 is EEDI = 37%, followed by the cargo capacity, represented by the number of containers transported, 22%, the construction cost of 21% and the operating cost of 20%.The risk is not seen as a significant factor for DS 1 .However, the DS 1 design solution has the minimum construction and operating cost and maximum cost associated with the risk.
The inverse calculation is performed to identify the importance of different criteria for the maximum relative closeness to the ideal solution.As can be seen from Figure 12, the best-suited design criterion estimated for the design solution, DS : .Associated with SHDMF = 1 is EEDI = 37%, followed by the cargo capacity, represented by the number of containers transported, 22%, the construction cost of 21% and the operating cost of 20%.The risk is not seen as a significant factor for DS1.However, the DS : design solution has the minimum construction and operating cost and maximum cost associated with the risk.The design solution, DS2, is dominated by the significance of the cargo capacity of 100%.Starting from DS3, where the significance of the risk dominates, where RISK = 96% and a significance of 4% is estimated for the resting parameters.The significance of the risk is estimated as 84% for DS4, with an operating cost of 18% and construction cost of 1%.For DS5 to DS11, the significance of the risk is 100%.
The multi-attribute decision-making approach employed here identifies the best design solution from the alternative design solutions, accounting for the annual operating cost, annual construction cost, risk, cargo capacity, and energy efficiency design index.Design solution 6 is the optimum design solution associated with a beta reliability index of 4.29, annual construction cost of EUR 527,732, operating cost of EUR 604,603, cargo capacity of 412 containers, and energy efficiency design index of 29.24 gCO2/tonne-mile.The design solution, DS 2, is dominated by the significance of the cargo capacity of 100%.Starting from DS 3 , where the significance of the risk dominates, where RISK = 96% and a significance of 4% is estimated for the resting parameters.The significance of the risk is estimated as 84% for DS 4 , with an operating cost of 18% and construction cost of 1%.For DS 5 to DS 11, the significance of the risk is 100%.
The multi-attribute decision-making approach employed here identifies the best design solution from the alternative design solutions, accounting for the annual operating cost, annual construction cost, risk, cargo capacity, and energy efficiency design index.Design solution 6 is the optimum design solution associated with a beta reliability index of 4.29, annual construction cost of EUR 527,732, operating cost of EUR 604,603, cargo capacity of 412 containers, and energy efficiency design index of 29.24 gCO 2 /tonne-mile.

Conclusions
The present study developed a design procedure to integrate the innovative riskbased ship hull structural design, capital and operational expenditure, cargo capacity, and energy efficiency for a multi-purpose ship in identifying an optimum level of ship safety, i.e., the optimum reliability index, controlling the risk associated with the ship hull structure.The identified design solution is associated with the minimum expected total cost, leading to lower construction and operational costs, energy efficiency, risk, and maximum cargo capacity.The study used MADM, considering several objectives simultaneously for different scenarios, employing the TOPSIS method, and identifying the best ship structural design solution for all scenarios as a function of SHDMF.It can be noticed that the procedure developed in this work is flexible enough to accommodate different design criteria and possible hazards during the service life of the ship.
The next step of the present study is to adapt the current procedure for the structural design of lightweight ships made of composite materials, including honeycomb sandwich panels replacing the conventional steel structures, to gain benefits of more cargo capacity and cost reduction in building modern energy-efficient and environment-friendly ships.

Figure 2 .
Figure 2. Bending moment-curvature relationship, SHDMF = 1, vertical (blue) is the bending moment axis, horizontal (green) is the curvature axis, resultant (red) is the hull girder strength and applied (pink) is the bending moment the hull girder needs to support.

Figure 2 .
Figure 2. Bending moment-curvature relationship, SHDMF = 1, vertical (blue) is the bending moment axis, horizontal (green) is the curvature axis, resultant (red) is the hull girder strength and applied (pink) is the bending moment the hull girder needs to support.

Figure 3 .
Figure 3. Beta reliability index β for net and gross scantling.

Figure 4 .
Figure 4. Cost of structural failure, safety measure, and total cost.

Figure 4 .
Figure 4. Cost of structural failure, safety measure, and total cost.

Figure 5 .
Figure 5. Safety measures, steel and labour cost, and weight of steel.

Figure 5 .
Figure 5. Safety measures, steel and labour cost, and weight of steel.

Figure 7 .
Figure 7. LW, DW and steel water vertical bending moment, Msw at midship section.

Figure 7 .
Figure 7. LW, DW and steel water vertical bending moment, Msw at midship section.

Figure 8 .
Figure 8. Steel water vertical bending moment at midship section, Msw, regular scale (left) and magnified scale (right).

Figure 8 .
Figure 8. Steel water vertical bending moment at midship section, Msw, regular scale (left) and magnified scale (right).

Figure 9 .
Figure 9. Transported containers, ship annual operating and construction costs.

Figure 9 .
Figure 9. Transported containers, ship annual operating and construction costs.

Table 1 .
Main Dimensions of Multi-Purpose Ship.