Lin Hua, Fan Wu*, Jinlei Mu and Xiaolong Ma

(1. School of Naval Architecture & Ocean, Naval University of Engineering, Wuhan 430033, China;2. Naval Vessel Technical Evaluation Center for Hull Strucutre Strength, Wuhan 430033, China)

Abstract: In order to solve the problem of calculating the ultimate strength of marine plate under pitting corrosion, the plate taking from marine with pitting corrosion was studied. Based on the data of pitting corrosion recorded, the concept of pitting corrosion characteristic element matrix was proposed to describe the distribution and the forms of pitting. Moreover, the model of sensitivity calculation was established to analyse the sensitivity of pitting corrosion characteristic elements for the ultimate bearing capacity of ship structural plate. A new approach was proposed to calculate the ultimate strength of pitting damaged ship structural plate under combined loads based on the factor of rigidity reduction. Simultaneously this new approach was proven to be validated by finite element simulation. Finally, formula of ultimate strength of pitting damaged ship hull plate was established by series of numerical calculations based on the approach，and the reliability of the formula was validated as well.

Keywords: ship structure; pitting damage; ultimate strength; combined loads; sensitivity analysis

1 Introduction

Under the influences of electrochemical corrosion, ocean organism adhered and non-metal impurities in metal materials, pitting corrosion damage usually exists in the hull plates. As the main structures of ship, the pitting corrosion of hull plates will lead to the ultimate bearing capacity of hull reduce dramatically, thus impacting the security of ship structures[1-2]. In recent years, as the cognition of pitting corrosion′s failure mechanism has been improved, more attentions have been paid to the impact of pitting corrosion to plate strength[3-5]. The statistical theory was used to analyse the impact of partial corrosion to plate buckling load, and the methods to calculate the critical pressure of plate′s instability were obtained[6]. The ultimate strength of ship structures under the condition of pitting corrosion was analysed through finite element analysis[7-8]. The relationships between distribution of pitting corrosion, positions, plate thickness and ultimate strength were obtained. Some researches on ultimate strength of steel plates with pitting corrosion carried out under experiment and finite element method[9-10]. The results showed that, in terms of plates owning just a single pitting corrosion, the position of pitting corrosion is unimportant for the ultimate compression strength, but extremely significant for post-buckling strength. Low carbon steel plates with pitting corrosion were analysed through nonlinear finite element analysis. One-side pitting corrosion resulted more in the decrease of compressive ultimate power when the thickness of pitting corrosion plate was the same as discovered, and many methods were taken to represent the level of pitting corrosion is imperfect[11-13].

Currently, in terms of the research on impacts of pitting corrosion to ultimate strength of plates, many researchers analyse the situations under axial load and many achievements have been obtained, while the research considering the influence of lateral load will be more useful[14]. As the axial ability bearing compression of bottom plate under horizontal hydrostatic pressure represents the safety level of the whole hull structure under hogging state. Therefore, the plate of real ship with pitting corrosion was taken as the object to analyse pitting characteristic parameters′ impact to the sensitivity of ultimate strength of plate by establishing mechanical and finite element model in this paper. Moreover, the concept of factor of rigidity reduction was proposed to appraise the level of plate′s pitting corrosion. Numerical finite element calculations were used to obtain the equation calculating the axial ultimate strength of plate under the condition of combined load.

2 Pitting Characteristic Matrix

There is still no effective method to precisely describe the characteristics of pitting corrosion due to the complicated distribution. To some extent, it is detrimental for research on pitting corrosion[15-16]. In order to solve this the plates of actual ship were taken as an example, the forms and positions of pitting corrosion were analysed in this paper. The concept of pitting characteristic matrix is proposed to study on the calculating methods of the ultimate strength of plates.

2.1 Pitting Corrosion Forms

The statistical data from actual ship of pitting corrosion shows that most of the shapes of pitting corrosion in bulk freighter are cone, while most are hemispheric in tankers, and most are echelon cone in military vessels[17]. The shapes of pitting corrosion of components such as ribs are usually spherical. Although the shapes are different in different structures and ships, they can be mainly divided into cone and hemisphere(as shown in Fig.1).

(a) Cone form pitting

(b) Hemisphere form pitting

In addition, research shows that the shapes of pitting corrosion just contribute little to the ultimate strength of components on condition of the same total volume[18]. Therefore, the cone and hemisphere shape of pitting corrosion could be transferred to cylinder shape in the same volume. The transformation equations as follows:

Cone shap→Cylinder shap：

Hemisphere shap→Cylinder shap：

wherehcis the depth of cone pitting corrosion,hsis the depth of hemisphere pitting corrosion,hcyis the depth of cylinder pitting corrosion, anddis the diameter of pitting corrosion.

2.2 Distribution of Corrosion Pitting

The data from collecting the situation of plates pitting corrosion of actual ships show corrosion distributed randomly in the plates under the influences of plate material situation (such as non-metal impurities and weld crack) and outer environment (such as temperature, humidity and PH value). According to the data, the pitting corrosion of above waterline is always one-side distributed, mainly on the outer side and with rarely perforating phenomenon, is shown in Fig.2. While at the areas close to waterline, it is always two-side distributed. The pitting characteristic parameters could be divided into geometrical parameters and distribution parameters in order to describe random distribution of pitting corrosion accurately. Geometrical parameters include diameter of pitting corrosiondand depth of pitting corrosionh. Distribution parameters mean the position of pitting corrosion on the coordinate, as shown in Fig.3.

Fig.2 Pitting distribution on ship plate

Fig.3 Pitting distribution on coordinate system

2.3 Pitting Characteristic Matrix Establishing

Pitting characteristic matrix was proposed in order to describe the shapes, sizes and positions of pitting corrosion of actual ship plates more precisely. The established steps as following:

Step1Take the middle surface of hull plate as coordinate planes, the length of plate asxaxis, the width of plate asyaxis, thus establishing in Dartesian coordinate system.

Step2Take the largest diametercof pitting corrosion as a standard to zone the hull plates and show in Fig.2. There is the number offcells along the length and the number ofgcells along the width of plate.

where “int” means down round-off number.

Step4The distribution characteristic matrix of pitting corrosion could be gained after integrating the dates and equations.

where, ifaij=0, means no-pitting and ifaij=1, means pitting.

3 Sensitivity of Pitting Characteristic Factors on Ultimate Strength

When the hull plate loses its stability under combined loads, the corresponding strength is defined ultimate strength. As the mechanism of pitting characteristic parameters influencing the ultimate strength of plates is really complicated, it is necessary to analyze the sensitivity of pitting characteristic parameters to test the level of influence on the plate ultimate strength.

3.1 Sensitivity Analysis Model of Pitting Parameters

The accuracy of analysis of sensitivity based on statistics is effective and reliable. An approach of selecting sample within the scope of parameter value was introduced to achieve analysis and the specific principle, and the method can be found in Ref.[19]. In this approach, the sensitivity of ultimate strength under the condition of combined loads was calculated based on the mechanical model and finite element model. The main parameters contain pitting position, diameter and depth.

3.1.1Mechanicalmodel

There are always beams, ribs and other support structures around hull plates, sometimes also longitudinal girder, cross beam, etc. Generally, the hull plates are more vulnerable than these structures, so the boundary of hull plates with pitting corrosion can be seen as simply support. Although the load of hull plates is complicated, it can be divided into two types: 1)The loads perpendicular to plate surface (such as gravity of cargo and hydrostatic pressure); 2)The axial loads on the plate[20](such as interaction force between plates and plane stress). For the bottom plate, the main load is horizontal hydrostaticqand axial loadTunder total longitudinal bend. The mechanical model is shown in Fig.4. The hydrostatic load is determined by ship draught, and in this caseq=0.03 MPa, axial pre-compressing loadT′=1 N/mm.

3.1.2Finiteelementmodel

Fig.4 Mechanical model of hull plate

Fig.5 Relation between strain and stress of material

3.2 Sensitivity Analysis

Sensitivity analysis represents the level of input variables′ influence as well as the contribution to output. The reason to study on sensitivity analysis is to confirm the influence level of pitting characteristics to the ultimate strength of plates under the setting condition of loading, thus gaining the fundamental knowledge for future research. Sensitivity analysis of pitting characteristic parameters is performed based on finite element analysis and Isight software through programming Umat program of Abaqus, by changing pitting positions, diameters and depths. The sensitivity calculation process is shown in Fig.6.

Fig.6 Flowchart of sensitivity analysis

The model of sensitivity analysis is that: length of platea=640 mm, widthb=480 mm, thicknesst=8 mm, elasticity modulusE=2.06×105MPa, Poisson ratiou=0.3. According to ships′ actual records, the diameter changes within [0 mm, 28 mm] and the depth changes within [0 mm, 4 mm], and the pitting corrosion position moves from the central to the edge. Owing to parameters have been normalized, the calculation results are not influenced by parameters′ dimensions. Therefore, the precision of sensitivity analysis has been improved a lot. The sensitivity analysis result of pitting corrosion parameters is shown in Fig.7.

As shown in Fig.7, the depth of pitting corrosion is the most influential factor, followed by position and diameter. If considering pitting characteristic factors′ interaction, the influence of diameter and depth of pitting corrosion is the primary factor, followed by the effect between diameter and position of pitting corrosion.

Fig.7 Results of sensitivity analysis

4 Calculation of Ultimate Load

Corrosion volume and degree of pitting intensity are usually taken as the criteria to evaluate the ultimate strength of plate with pitting corrosion. Although those methods are applicable to some extent, rigidity reduction is more reasonable if considering the thin-wall effects and position of pitting corrosion. Therefore, the concept of rigidity reduction factor is proposed to calculate the ultimate load of hull plates with pitting corrosion based on pitting characteristic matrix.

4.1 Calculation of Rigidity Reduction Factor

The essence of pitting corrosion influences on the strength of plates is reducing the stiffness of plates. Given the random distribution of pitting corrosion can decentralize the stiffness and the flexibilityβof plates is an important parameter to reflect the buckling strength of plates[21]. The hull plate with pitting corrosion is taken as thickness-changing plate along with width direction in this paper. The method to calculate rigidity reduction factorλis shown as follows:

Step1Based on pitting characteristic matrix, the volume of plate in a cell isV=tc2, the equivalent thicknesstiin each unit along the width of hull plate is

whereaiis the element in distribution characteristic matrix,hiis pitting depth, anddiis pitting diameter.

Step2The stiffness of hull plate after pitting corrosion is

whereZiis the height centroid of unit width.

Step3The initial stiffness of hull plate is

D=EI

whereIis the inertia moment of marine plate before pitting corrosion.

Step4The rigidity reduction factor of hull plate is

4.2 Calculation of Ultimate Strength

It is necessary to verify the validity of rigidity reduction factor method before it is used to calculate the ultimate strength of hull plates with pitting corrosion. In this case, the axial ultimate load of actual hull plate is calculated in the way of rigidity reduction and actual pitting corrosion situation respectively, based on established mechanical model and finite element model. The specific parameters of hull plates are shown in Table 1.

Distribution characteristic matrix of pitting corrosion of hull plates isA′=[aij],i=1,…,24,j=1,…,32. The value of elementsa5,5,a12,2,a15,2,a11,9,a20,8,a5,25,a19,24is 1, and the remainder is 0. These seven corrosion pitting′ geometrical characteristic elements are shown in Table 2.

Table 1 Parameters of actual ship structural plate

Table 2 Geometrical elements of pitting corrosion (mm)

Based on the data from Tables 1 and 2, the axial load-displacement curve for the node which displacement is the largest in pitting plate has been obtained by different calculations methods through finite element under the horizontal uniform loadq=0.03 MPa, and is shown in Fig.8.

Fig.8Axialcompressionload-displacementcurvesofpittingplate

4.3 Calculating Formula of Ultimate Strength

The validity of rigidity reduction factor method has been demonstrated, and the evaluation formula of ultimate strength of hull plates with pitting corrosion under combined load is proposed based on finite element numerical calculation. According toClassificationRulesofSteelSea-goingVesselfrom CCS, hull plates would be fixed or replaced when the level of pitting corrosion is reached. So that the actual range of rigidity reduction factorλis [0, 0.48) and the thickness range of actual hull plates is always [6 mm, 12 mm]. The changing rule between axial ultimate loadTcand rigidity reduction factorλunder different flexibility ofβwhen the horizontal hydrostatic pressureq=0.03 MPa is obtained based on series of finite element numerical calculation, and shown in Fig.9.

Fig.9RelationbetweenTcandλunderdifferentflexibility

As shown in Fig.9, the dimension of plate has a remarkable impact on the axial ultimate load. The ultimate strength reduction ratioTc/Tu(whereTcis the axial ultimate load considering pitting andTuis the axial ultimate load without pitting) was taken to weigh the damage degree of plate pitting to remove the effect of dimension. The changing trend betweenTc/Tuandλis shown in Fig.10. As been shown, the geometrical sizes of plate are less influential to the ultimate strength reduction ratioTc/Tu, thus the purpose of removing the size effect of plates has been achieved.

Fig.10RelationbetweenTc/Tuandλunderdifferentflexibility

The relationship between ultimate strength reduction ratioTc/Tuand rigidity reduction factorλis obtained based a series of finite element numerical calculations, and the rule was shown in Fig.11. The evaluating formula of ultimate strength of hull plates with pitting corrosion under condition of horizontal hydrostatic pressureq=0.03 MPa was obtained by regression analysis, and as follows:

(1)

In terms of the formula of ultimate strength of hull plates with pitting corrosion, the correlation coefficientR2is 0.999 8, and the standard deviationSis 9.75×10-5. The maximum deviation between Eq.(1) and finite element calculations is 3.25%. Therefore, applying the Eq.(1) to evaluate the ultimate strength of hull plates in that case is reliable and applicable. The evaluating formula of ultimate strength of hull plates with pitting corrosion under different horizontal hydrostatic load would be obtained by the method presented in the paper. However, for those heavy pitting corrosion with holes, more researches need to be performed because of the hull plate′s failure mode has been transformed.

Fig.11RelationbetweenTc/Tuandλundernumericalcalculationandfittingformula

5 Conclusions

The concept of pitting characteristic matrix has been proposed based on the analysis of the characteristic factors of actual hull plates with pitting corrosion as well as sensitivity of pitting parameters analyzed. The method to estimate the ultimate strength of hull plates with pitting corrosion under certain horizontal hydrostatic load has been carried out by the way of ultimate strength reduction ratio based on the sensitivity analysis and pitting characteristic matrix. The formula to evaluate the ultimate strength of hull plates with pitting corrosion which is proven to be feasible has been presented. The conclusions can be drawn as follows:

1) The depth of pitting is the most influential factor to the ultimate strength of hull plates with pitting corrosion under the condition of single influence of factor, while the inter-influence between diameter and depth of pitting becomes the most influential factor, followed by inter-influence between diameter and distribution of pitting when considering interactive influence between factors.

2) The method to estimate the pitting corrosion′s influence on ultimate strength of hull plates in the way of rigidity reduction is feasible and rational. The formula of calculating ultimate strength of hull plate with pitting corrosion that has been obtained under combined loads in the way of rigidity reduction is reliable and applicable in engineering. Due to the standard deviationSis 9.75×10-5, and the correlation coefficientR2is 0.999 8, the formula could be the criteria and relevant standard to calculate the ultimate strength of hull plates with pitting corrosion.

3) The concepts of pitting characteristic matrix and rigidity reduction factor are proposed in this paper, which develop a new way to study the pitting damage. The formula is established to evaluate the degree pitting damage could make a great contribution to further research in this area.