Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 125856, 15 pages

http://dx.doi.org/10.1155/2013/125856

Research Article

Time Intervals for Maintenance of Offshore Structures Based on Multiobjective Optimization

Dante Tolentino and Sonia E. Ruiz

Instituto de Ingeniería, Universidad Nacional Autónoma de México, 04510 Coyoacán, DF, Mexico

Received 22 March 2013; Revised 28 June 2013; Accepted 5 July 2013

Academic Editor: Dan Simon

Copyright © 2013 Dante Tolentino and Sonia E. Ruiz. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

With the aim of establishing adequate time intervals for maintenance of offshore structures, an approach based on multiobjective optimization for making decisions is proposed. The formulation takes into account the degradation of the mechanical properties of the structures and its influence over time on both the structural capacity and the structural demand, given a maximum wave height. The set of time intervals for maintenance corresponds to a balance between three objectives: (a) structural reliability, (b) damage index, and (c) expected cumulative total cost. Structural reliability is expressed in terms of confidence factors as functions of time by means of closed-form mathematical expressions which consider structural deterioration. The multiobjective optimization is solved using an evolutionary genetic algorithm. The approach is applied to an offshore platform located at Campeche Bay in the Gulf of Mexico. The optimization criterion includes the reconstruction of the platform. Results indicate that if the first maintenance action is made in 5 years after installing the structure, the second repair action should be made in the following 7 to 10 years; however, if the first maintenance action is made in 6 years after installing the structure, then the second should be made in the following 5 to 8 years.

1. Introduction

Structures are continually exposed to different environmental loads that modify both their mechanical properties and their structural performance over time. Due to these modifications, it is useful to formulate maintenance plans for the structural systems. Several approaches of maintenance plans have been proposed in the literature considering different numbers of optimization objectives: that is, (a) minimizing the expected cumulative maintenance cost over a time interval (as a single objective) [1], (b) minimizing maintenance cost and maximizing load-carrying capacity and durability (as two independent objectives) [2], and (c) more powerful approaches that optimize several objectives simultaneously [3–6].

In this paper a multiobjective optimization approach is proposed. Such approach considers a balance among the following objectives: (1) reliability, expressed in terms of the confidence factor that takes into account the structural deterioration over time by means of closed-form mathematical expressions, (2) damage index, expressed in terms of the ratio between the structural capacity and the structural demand, and (3) expected cumulative total cost, which considers costs of inspection, repair, and failure. The optimization problem is solved by means of the elitist Nondominated Sorting Genetic Algorithm (NSGA II) [7] which gives place to a set of nondominated solutions that are used to recommend time intervals for making decisions for maintenance of an offshore structure.

For the particular case of offshore structures, the structural deterioration phenomenon is mainly due to fatigue caused by waves acting continuously on the steel elements. The decrease of the resistance of the structural capacity over a time interval is caused by the cracking of the tubular joints. The inspection of an offshore structural system has two objectives: (a) to detect both the presence and the size increase of cracks and (b) to perform the necessary maintenance actions to the structure accordingly. Many authors have proposed inspection and maintenance plans for offshore structures using different approaches such as the following: (a) the analysis of risk and reliability of welded connections subject to fatigue [8–10]; (b) the use of methodologies that take into account a fatigue sensitivity analysis in steel joints [11]; (c) the application of probabilistic detection methods in order to study the influence of repetitive inspection for fatigue in joints [12]; (d) the implementation of simplified approaches and the use of practical design parameters such as fatigue design factors [13] and/or reserve strength ratios [14]; (e) the consideration of damage caused by fatigue, buckling, and dents on structural elements [15], and (f) the use of Bayesian techniques [16].

However, neither a multiobjective criterion for making decisions of maintenance actions on “jacket” marine platforms, nor simplified closed-form mathematical expressions for taking into account the structural deterioration of both structural capacity and structural demand over a time interval are considered in the studies mentioned above, as it is done in the present study. The advantage of using closed-form mathematical expressions is the possibility of evaluating the structural reliability without having to perform numerical integrations.

2. Basic Definitions of Structural Reliability, Damage Index, and Expected Cumulative Total Cost

2.1. Structural Reliability

The structural reliability can be expressed by means of different indicators, such as probability of failure in a time interval, expected number of failures per unit time, confidence factors, and confidence levelss. In this paper, the confidence factor, , is used as a measure of the reliability implicit in the structure. The confidence factor was originally introduced in order to evaluate the reliability of steel building structures [17–20]. Later, the original methodology was extended to the case in which the mechanical properties of the structural elements present degradation (i.e., due to corrosion, fatigue, etc.), and as a consequence, either the structural capacity or the structural demand, or even both change over time [21–24]. In the present study the confidence factor is evaluated under the assumption that the structural capacity deteriorates linearly over time. For this case, the confidence factor, , is calculated as mentioned in Appendix A.

2.2. Damage Index

The damage index (DI) is defined as the ratio of the structural capacity, , and the structural demand, . The damage index varies between 0 and 1, where 0 indicates that the structure does not have any damage, and 1 refers to total damage, which means that the structural system has reached its failure condition.

2.3. Expected Cumulative Total Cost

The expected cumulative total cost function is defined as the sum of the expected cost of inspection, repair, and failure, at the end of a time interval , as derived in [25] and presented in Appendix B.

3. Basic Definitions of Genetic Algorithms

Genetic algorithms (GAs) are adaptive methods used to solve, search, and optimize problems. Based on the principle of the survival of the fittest [26], GAs have been applied successfully in many problems and different areas (i.e., economics, medicine, engineering, etc.). GAs are commonly composed of three basic operators: (1) selection, (2) crossover, and (3) mutation. Recently, GAs have been used by civil engineers due to the ease of handling multiple goals directly and simultaneously [27, 28]. In this study, GAs are selected as a tool of optimization in order to solve a problem in which three objectives must be satisfied. Commonly, in multiobjective optimization a set of nondominated solutions is obtained, all of which form a Pareto front. The concept of dominance is as follows [29].

Assuming that all objective functions are to be maximized (see Figure 1), it is said that the solution is dominated by solution if conditions 1 and 2 are true.

(1)Solution is not worse than considering all the objectives.(2)Solution is strictly better than in at least one objective.

Figure 1: Dominated and nondominated solutions.

In order to find a Pareto front that helps the decisions makers to formulate plans of maintenance of offshore structures which present structural degradation over a time interval, the elitist NSGA II method [7] was used in the present study. The approach consists of the following steps:

(1)Randomly create an initial population of size , where generation index (). (2)Selection, crossover, and mutation are made to create the offspring population of size .(3)Combine parent, , and offspring, , to create the population of size .(4)A fast nondominated sort of size is made to create the population of size .(5)Repeat the process from the second step using the population selected above to create the offspring population, .

The offspring population is obtained using binary tournament selection, intermediate crossover, and Gaussian mutation. These concepts are explained below.

Binary Tournament Selection. individuals (commonly 2) are selected randomly from the size population. The selected individuals compete against each other. The individual with the highest fitness win, and it will be included for crossover and mutation [30].

Intermediate Crossover. It is the process of creating offsprings by a weighted average of the parents, and it is controlled by the following ratio [31, 32]:

The value of “ratio” is commonly assumed equal to 1.2 [32].

Gaussian Mutation. It consists in adding a random number taken from a Gaussian distribution with mean 0 to each entry of the parent vector. The standard deviation is determined by the parameters “scale” and “shrink,” as follows:

where the “scale” parameter determines the standard deviation at the first generation; the “shrink” parameter controls how the standard deviation shrinks as generations go by; is the current generation, and represents the number of total generations under consideration. The values of “scale” and “shrink” are commonly assumed equal to 1.0 [32].

Fast Nondominated Sort. A “naive and slow” procedure [29] is applied in order to sort a population into different nondomination levels. In this procedure each solution is compared with every other solution in order to find out if it is dominated by any other solution in the population. When the comparisons of each solution are completed, all individuals in the first nondominated front are found. All individuals obtained in the first nondominated front are set apart temporarily in order to find the second nondominated front, and the procedure is repeated. The third nondominated front is obtained by setting apart temporarily all the individuals in the first and second front, and the procedure mentioned above is repeated. In summary, different fronts are obtained in this approach, where is the best nondominated set, is the next best nondominated set, and so on. In case of , all members of the set are chosen to create the new population ; the remaining members are chosen from the subsequent nondominated sets in the order of their ranking (). If , the members are chosen using a crowded-comparison operator [29], in descending order and the best solutions needed to fill the population slots are selected.

4. General Approach of the Multiobjective Optimization for Structural Maintenance

The multiobjective optimization criterion for maintenance of offshore structures formulated in the present study considers three objectives: confidence factor , damage index (DI), and expected cumulative total cost (E[CTC]). The first two objectives are subjected to the following constraints:

(i)the confidence factor is subject to .(ii)the damage index is subject to .

That is, the structural reliability expressed in terms of the confidence factor, , must be greater than unity [18] (values less than 1 consider that the structure is in the unsafe side). The damage index (DI) should be smaller or equal to a given threshold which is related to a permissible level of structural damage.

The procedure to obtain the set of nondominated solutions is divided here in two general steps: (a) the simulation of values of the three objectives (, DI, and E[CTC]) over a time interval, as indicated in the flowchart of Figure 2, and (b) the codification of the simulations, to perform the multiobjective optimization by means of NSGA II.

Figure 2: Flowchart corresponding to the three objectives simulation.

The simulation profiles are codified in order to perform the multiobjective optimization by using the GA mentioned in the previous section. Then, the optimization problem is solved as follows:

(1)maximize the minimum confidence factor;(2)minimize the maximum damage index;(3)minimize the expected cumulative total cost.

As a result, a set of nondominated solutions is found for making decisions and for establishing different time interval for maintenance. In this study, the criterion for selecting the maintenance time interval is based on a threshold value of the confidence factor. When the confidence factor takes values between a certain range of values of the nondominated solutions, maintenance actions should be performed. The range of values to select the repairing time interval is as illustrated in Figure 3, where the nondominated solutions in the bracket are those selected for repairing actions.

Figure 3: Criterion to select repairing time intervals.

Here it is assumed that after several nondominated solutions are selected for maintenance actions, the critical structural elements are repaired to recover their initial capacities, and a new set of nondominated solutions can be obtained after reconstruction. Then, a second set of nondominated solutions is obtained following the flowchart shown in Figure 2 (as done previously); however, an initial cumulative damage (ICD) is now taken into account for the second set of critical structural elements to be analyzed. It is noticed that in the example presented in Section 5 the second set of critical structural elements is different from the first set analyzed in the first step. For illustrative purposes, Figures 4(a), 4(b), and 4(c) show projections of one selected solution (indicated with circles) corresponding to the first maintenance time interval, as well as one solution (indicated with triangles) corresponding to the second time interval for structural maintenance.

(a)

(b)

(c)

Figure 4: Projections of selected solutions corresponding to the first and second time intervals for repairing actions.

5. Illustrative Example

The multiobjective approach proposed above is applied to an offshore “jacket” platform with the aim of finding its first and its second time interval for maintenance. The platform is supposed to be located in Campeche Bay. The structure is 48 m high and the water depth at the site is 45.11 m (see Figure 5). The expected weight of the deck was assumed equal to 500 ton. The “jacket” platform is represented by a simplified 2D model with mean mechanical properties. The joints were selected according to their global structural capacity contribution. Joints 1, 2, 3, 4, and 6 were selected for finding the first time interval for maintenance, whereas joints 5, 7, 8, and 9 were selected for finding the second time interval (once the joints 1, 2, 3, 4, and 6 were totally reconstructed).

Figure 5: 2D model of the “jacket” platform.

5.1. Evaluation of Structural Capacity over Different Time Intervals

The lateral resistance of the structure was obtained by means of nonlinear static analyses using twenty different simulated profiles of lateral loads. These profiles are related to the acting forces corresponding to the simulated waves that produced the maximum base shear on the structure. The appearance and growth of cracks, which give place to a reduction of the structural capacity, are considered as damage condition. Two crack points were considered for each structural element connected to the joints.

The average crack size simulation process of the selected points subject to random load was performed using the Monte Carlo technique. Operational and storm waves were considered. The arrival times between storms were assumed with exponential distribution. The storm wave heights follow a Gumbel distribution. The statistical parameters used for simulation of the crack size were taken from Silva and Heredia [33]. The crack size was obtained using the modified differential equation by Paris and Ergodan [34] and Sobczyk and Spencer [35]. In the differential equation, the random load was replaced by an equivalent cyclic load whose amplitude and frequency were expressed in terms of the mean properties of the random process. After obtaining the time history of the crack size, the capacity of the joint was reduced by a linear correction factor as follows [36, 37]:

where is the capacity of the cracked joint; is the capacity of the intact joint which can be calculated by means of the API recommendations [38]; and are the areas of the transverse section and of the crack, respectively, corresponding to all the elements connected to the joint. As an example, percentages of simulated cracked areas at the end of different time intervals, corresponding to joint 2, are shown in Figure 6.

Figure 6: Percentages of simulated cracked areas. Joint 2.

The median of structural capacity, , represented in terms of the global displacement (in meters, m) of the platform, and the standard deviation of its natural logarithm, are shown in Table 1 for different time intervals ( years). It is assumed that both the structural capacity and the structural demand follow a lognormal distribution in any time [17–21].

Table 1: Median of structural capacity, , and standard deviation of the natural logarithm of the structural capacity at different time intervals.

(years)
(m)

0
0.2520
0.0140

1
0.2516
0.0574

2
0.2492
0.0901

3
0.2462
0.1054

4
0.2429
0.1130

5
0.2378
0.1500

6
0.2322
0.1813

7
0.2297
0.1902

8
0.2272
0.1918

9
0.2206
0.1971

10
0.2190
0.1978

11
0.2179
0.2079

12
0.2151
0.2122

13
0.2116
0.2137

14
0.2106
0.2305

15
0.2092
0.2334

5.2. Evaluation of the Structural Demand over Different Time Intervals

The structural demand (given a maximum wave height) at the end of a time interval was obtained by means of dynamic “step by step” analysis, using simulated time histories of waves that correspond to different maximum wave heights. The crack growth simulations and the capacity reductions at selected joints considered to evaluate the capacity over time (Section 5.1) were also used in this section as damage condition. Figures 7(a) and 7(b) show the median structural demand as a function of the platform global displacement, and the corresponding standard deviation of the natural logarithm of the demand, for different intervals , respectively. The expression that fits the median values of the structural demand, , given a maximum wave height, , at the end of different time intervals , is given by (a similar expression was proposed by Tolentino et al. [24]). According to the above, the median structural demand, given a maximum wave height, is equal to . Similarly, the corresponding standard deviation of the natural logarithm of the demand is given by .

(a)

(b)

Figure 7: (a) Median of structural demand and (b) standard deviation of the natural logarithm of the demand, for different time intervals.

5.3. Damage Index

In this study the damage index is defined as follows:

where is the maximum global displacement of the structure without damage, represents the maximum global displacements that consider structural deterioration at instant , and is the ultimate global displacement demanded by the structure at instant .

5.4. Confidence Factor and Expected Cumulative Total Cost

The confidence factor values, , were calculated using . The parameters and were fitted according to , based on the median values of the structural capacity, , shown in Table 1. The epistemic uncertainties associated with the structural demand, , and with the structural capacity, , were assumed equal to 0.15. The parameters and fitted to the wave hazard curve were equal to and , which correspond to a maximum wave height equal to  m and to a return interval of 1485 years for the failure condition [39].

The expected cumulative total cost was evaluated using . An annual discount rate was assumed. The inspection cost by joint was 3,518 USD [40] and the repair cost was 20,000 USD for each critical point [41]. The cost of failure includes costs due to equipment damage, pollution, deferred production, and indirect losses [25].

5.5. Simulation of the Three Objectives

1024 profiles of damage index, confidence factor, and expected cumulative total cost, considering cumulative damage at the end of a time interval , were simulated, following the procedure presented in Figure 2. The simulations were based on the median values and the standard deviations of structural capacity and structural demand. Lognormal distribution functions were assumed.

5.6. Multiobjective Optimization before Reconstruction (First Time Interval for Maintenance)

The simulations were codified in order to perform the optimization by means of NSGA II. A population size of 50 was considered. A set of 12 optimal solutions in 24 generations were found. Different views of the Pareto front that include damage index, confidence factor, and the expected total cost over a time interval are shown in Figures 8(a), 8(b), and 8(c). Figure 8(d) shows the same results in 3D view. It can be observed that the damage index varies between 0.23 and 0.29, the confidence factor presents values between 1.10 and 1.24, and the expected cumulative total cost takes values between and USD. In order to compare different time intervals to perform maintenance actions, six optimal solutions (A, B, C, D, E, and F) were selected (see Figures 8(a)–8(d)). A comparison among the solutions selected is shown in Table 2.

Table 2: Comparison between solutions A, B, C, D, E, and F.

Solution
Time interval (years)
Damage index (DI)
Confidence factor ()
Expected cumulative total cost E[CTC]

(USD)

A
6
0.2345
1.0965
1105415

B
5
0.2387
1.1375
930566

C
5
0.2421
1.1452
850778

D
4
0.2762
1.2246
786761

E
4
0.2879
1.2381
818191

F
3
0.2939
1.2378
665032

(a)

(b)

(c)

(d)

Figure 8: Pareto front corresponding to the three objective functions and optimal solutions A, B, C, D, E, and F.

With the aim of establishing time intervals of structural maintenance, solutions associated with values of were selected for repairing actions. Then, in accordance to Table 2, maintenance actions on joints 1, 2, 3, 4, and 6 should be performed between 5 and 6 years (corresponding to solutions A, B, and C).

5.7. Multiobjective Optimization after Reconstruction (Second Time Interval for Maintenance)

In this section it is assumed that the maintenance (repairing) actions on joints 1, 2, 3, 4, and 6 are performed in accordance with solution A, and alternatively, with solution . It means that the first repair is done after 6 years of installing the platform (solution A), or alternatively, after 5 years if solution is selected (see Table 2). Also, it is considered that after maintenance, the joints recover their full capacity. Taking this into account, the joints that should be repaired in the next maintenance plan are joints 5, 7, 8, and 9 (which were selected in accordance with their contribution to global structural capacity).

The evaluation of the structural capacity and the structural demand (given a maximum wave height) at the end of a time interval , after the first maintenance action, was performed in a similar way as in Sections 5.1 and 5.2; however, it was now considered that joints 5, 7, 8, and 9 presented an initial cumulative damage (ICD). The ICD takes into account the deterioration due to fatigue that the joints suffered during the time interval between the installation of the structure and its first repair. Figure 9 shows an example of the percentages of cracked areas of joint 5 at the end of different time intervals, considering an ICD corresponding to 5 years (which is associated with solution C). It is noticed that in Figure 9 there is an initial cumulative damage (represented by /  × 100) at years (on the contrary, Figure 6 does not indicate any initial cumulative damage in joint 2).

Figure 9: Percentages of the simulated cracked area on joint 5, with an ICD corresponding to solution C.

Tables 3 and 4 show the median of the structural capacity, , represented in terms of global displacement of the platform, corresponding to solutions and C, respectively, for different time intervals after the first maintenance actions. The standard deviations of the natural logarithm of the structural capacity are also shown in Tables 3 and 4.

Table 3: Median of the structural capacity, , and standard deviation of the natural logarithm, at the end of different time intervals, corresponding to solution A.

(years)
(m)

0
0.2530
0.0129

1
0.2515
0.0336

2
0.2499
0.0885

3
0.2474
0.1141

4
0.2448
0.1412

5
0.2420
0.1592

6
0.2373
0.1803

7
0.2365
0.1870

8
0.2325
0.1966

9
0.2304
0.2026

10
0.2296
0.2033

11
0.2286
0.2193

12
0.2249
0.2176

13
0.2246
0.2247

14
0.2227
0.2362

15
0.2213
0.2578

Table 4: Median of the structural capacity, , and standard deviation of the natural logarithm, at the end of different time intervals, corresponding to solution C.

(years)
(m)

0
0.2534
0.0146

1
0.2518
0.0383

2
0.2503
0.0863

3
0.2495
0.1196

4
0.2476
0.1250

5
0.2468
0.1398

6
0.2448
0.1433

7
0.2417
0.1509

8
0.2392
0.1602

9
0.2375
0.1668

10
0.2364
0.1702

11
0.2342
0.1841

12
0.2315
0.1915

13
0.2300
0.1920

14
0.2285
0.2026

15
0.2273
0.2332

The median of the structural demand (given a maximum wave height, ) was calculated as in Section 5.2. The expressions obtained are and , for solutions and C, respectively. The standard deviations of the natural logarithm of the demand, for a maximum wave height, at the end of different time intervals were fitted as and for solutions A and C, respectively.

A total of 1024 simulations corresponding to solution A, and alternatively, to solution were performed following the procedure presented in Figure 2. The simulations were codified for each case in order to realize the optimization procedure using NSGA II. A population size of 50 was considered. A set of 13 optimal solutions with 23 generations for solution was found, and 12 optimal solutions with 25 generations for solution were obtained. Figures 10 and 11 show the set of solutions (Pareto fronts) for cases and C, respectively. Next, some optimal solutions are selected in order to handle different options for structural maintenance. A comparison of two sets of solutions selected is shown in Tables 5 and 6 (see Figures 10 and 11).

Table 5: Comparison between solutions a₁, b₁, c₁, d₁, e₁, and f₁. ICD corresponding to 6 years.

Solution
Time interval (years)
Damage index

(DI)
Confidence factor ()
Expected cumulative total cost [CTC]

(USD)

8
0.2277
1.0859
2319250

7
0.2388
1.0935
2190760

6
0.2435
1.1119
2095360

5
0.2448
1.1136
2014650

e₁
3
0.2489
1.1630
1794310

f₁
3
0.2535
1.1751
1756520

Table 6: Comparison between solutions a₂, b₂, c₂, d₂, e₂, f₂, and g₂. ICD corresponding to 5 years.

Solution
Time interval (years)
Damage index

(DI)
Confidence factor ()
Expected cumulative total cost [CTC]

(USD)

10
0.2284
1.0827
2133420

9
0.2309
1.0886
2040560

8
0.2371
1.1200
1870170

7
0.2375
1.1255
1813980

e₂
5
0.2680
1.2029
1665650

f₂
4
0.2734
1.2284
1523060

g₂
4
0.2980
1.2458
1554500

Figure 10: Pareto front with 13 optimal solutions considering solution for the first time interval for maintenance and a set of solutions selected (, , , , , and ) for the second interval for maintenance.

Figure 11: Pareto front with 12 optimal solutions considering solution for the first time interval for maintenance and a set of solutions  selected (, , , , , , and ) for the second interval for maintenance.

In a similar way to the previous section, the time intervals for the second maintenance actions are performed when the structural confidence values are lower than 1.15. Tables 5 and 6 show that maintenance actions could be done at time intervals between 5 and 8 years (solutions , , , and ) if the first maintenance actions were done after 6 years of installing the platform, or alternatively, at intervals between 7 and 10 years (solutions , , , and ) if the first maintenance actions were performed after 5 years of installing the platform.

Figures 12(a), 12(b), and 12(c) show projections of solution A (see Table 2) together with projections of solutions , , , and corresponding to the three objectives (damage index, confidence factor, and expected cumulative total cost). Similarly, Figures 13(a), 13(b), and 13(c) show projections of solution together with projections of solutions  , , , and .

(a)

(b)

(c)

Figure 12: Different projections of solution and of solutions , , , and .

(a)

(b)

(c)

Figure 13: Different projections of solution and of solutions , , , and .

In summary, in accordance with the criterion selected for performing maintenance actions (values of ), Figures 12(a) and 12(b) show that the first optimal time interval for reconstructing joints 1, 2, 3, 4, and 6 can be after 6 years of installing the offshore “jacket” platform (see Table 2, solution A). If this solution is selected, then the second optimal time interval for repairing actions on joints 5, 7, 8, and 9 should be after 5, 6, 7, or 8 years (solutions , , , or ) of reconstructing the structure. On the other hand, if the maintenance actions for the first time are done after 5 years of installing the marine structure (solution in Table 2), then the repairing actions on joints 5, 7, 8, and 9 should be performed after 7, 8, 9, or 10 years after reconstructing the structure (solutions  , , , or ), as shown in Figures 13(a) and 13(b).

6. Conclusions

An approach based on multiobjective optimization in order to find adequate time intervals for maintenance, by means of genetic algorithms, is applied. A set of nondominated solutions is found based on three objectives: (1) confidence factor, which is represented by closed-form mathematical expressions that take into account the structural deterioration over time, (2) damage index, expressed in terms of the ratio between the structural capacity and structural demand, and (3) expected cumulative total cost, which considers inspection, repair, and failure costs.

The evolutionary Nondominated Sorting Genetic Algorithm NSGA II was selected as a tool to solve the optimization problem. Different sets of nondominated solutions were found in order to recommend different time intervals for maintenance actions of a fixed steel “jacket” platform subjected to structural deterioration over time. The implemented GA is very efficient since it only requires a few seconds to find a set of nondominated solutions based on different objectives; it also requires less than 30 generations in order to obtain the Pareto front.

The results show that, for the offshore “jacket” platform analyzed, it is recommended to make the maintenance actions on the first set of selected joints, between 5 and 6 years after installing the platform. In the case of performing the first maintenance action after 5 years, the second time interval for repairing the second set of joints is recommended to be done between 7 and 10 years. On the other hand, if the first maintenance actions are performed after 6 years of installing the platform, the second time interval for maintenance actions should be between 5 and 8 years.

The multiobjective optimization approach presented, together with the closed-form mathematical expression used, constitutes an efficient tool for making maintenance plans for offshore structures; moreover, it can be adapted to consider different kinds of structures (buildings, transmission towers, bridges, etc.). The criterion can be extended to take into account different types of maintenance actions in order to keep the structural system within prescribed performance levels.

A limitation of the approach presented here is that it is adequate only for making decisions about structural maintenance actions for systems subjected to ordinary excitations; it does not apply to the identification of repair actions that should be taken after the occurrence of extraordinary environmental events, for example, after the action of hurricanes with a return period much longer than that considered here (see Section 5.4).

Appendices

A. Confidence Factor That Takes into Account the Structural Deterioration

The confidence factor over the time interval can be expressed as follows [21, 24]:

where

where and are the capacity reduction factor and the demand amplification factor at instant , respectively; and are the variances of the natural logarithm of the demand, , given a maximum wave height, , and of the structural capacity, , respectively; and represent the epistemic uncertainties associated with the structural demand and the structural capacity respectively; is the median of the structural capacity at the limit state of interest, which vary linearly in the interval ; is the median of the structural demand caused by environmental loads related to the tolerable annual exceedance rate , at time ;   represents the correction function of the expected number of failures which considers the variation in the interval of the structural demand (given a maximum wave height) and the structural capacity simultaneously.

Hypergeometric functions implicit in (A.6) can be solved using a simpler method by means of the following hypergeometric series [42, 43]:

where

where , , and represent real numbers, and must be an integer. The value will depend on the problem; for example, for the case analyzed in this study, the hypergeometric series is developed up to , in order to obtain a good approximation.

B. Expected Cumulative Total Cost

The expected cumulative total cost function is defined as the summation of the total expected cost of inspection, repair, and failure at the end of a time interval as follows [25]:

where ,  , and are the expected cost of inspection, repair, and failure, respectively; is the discount rate; is the expected number of failures at the end of a time interval considering the variation of the structural capacity and structural demand over time [24]; is the number of elements to be repaired; is the conditional probability that element with certain level of cumulative damage , at time , is higher or equal to an established damage level ; is the number of time intervals of interest; is the expected number of failures for the time intervals of interest, and years. Equation (B.1) implies that the structure survives up to the end of the time interval of interest.

Acknowledgments

This research project had the support of DGAPA-UNAM (PAPIIT IN107011-3). The first author thanks CONACYT for the economic support to develop his PhD research.

References

J. S. Kong and D. M. Frangopol, “Life-cycle reliability-based maintenance cost optimization of deteriorating structures with emphasis on bridges,” Journal of Structural Engineering, vol. 129, no. 6, pp. 818–828, 2003.

A. Miyamoto, K. Kawamura, and H. Nakamura, “Bridge management system and maintenance optimization for existing bridges,” Computer-Aided Civil and Infrastructure Engineering, vol. 15, no. 1, pp. 45–55, 2000.

M. Liu and D. M. Frangopol, “Optimal bridge maintenance planning based on probabilistic performance prediction,” Engineering Structures, vol. 26, no. 7, pp. 991–1002, 2004.

H. Furuta, T. Kameda, and M. Erami, “A practical bridge management system using new multi-objective genetic algorithm,” in Bridge Maintenance, Safety, Management, Life-Cycle Performance and Cost, P. J. S. Cruz, D. M. Frangopol, and L. C. Neves, Eds., pp. 131–132, Taylor & Francis Group, London, UK, 2006.

L. A. C. Neves, D. M. Frangopol, and P. J. S. Cruz, “Probabilistic lifetime-oriented multiobjective optimization of bridge maintenance: single maintenance type,” Journal of Structural Engineering, vol. 132, no. 6, pp. 991–1005, 2006.

N. M. Okasha and D. M. Frangopol, “Lifetime-oriented multi-objective optimization of structural maintenance considering system reliability, redundancy and life-cycle cost using GA,” Structural Safety, vol. 31, no. 6, pp. 460–474, 2009.

K. Deb, A. Pratap, S. Agarwal, and T. Meyarivan, “A fast and elitist multiobjective genetic algorithm: NSGA-II,” IEEE Transactions on Evolutionary Computation, vol. 6, no. 2, pp. 182–197, 2002.

R. Skjong, “Reliability based optimization of inspection strategies,” in Proceedings of the 4th International Conference on Structural Safety and Reliability (ICOSSAR '85), vol. 3, pp. 614–618, Kobe, Japan, 1985.

H. O. Madsen, J. D. Sørensen, and R. Olesen, “Optimal inspection planning for fatigue damage of offshore structures,” in Proceedings of the 5th International Conference on Structural Safety and Reliability (ICOSSAR '89), pp. 2099–2106, San Francisco, Calif, USA, August 1989.

J. D. S. Sørensen, M. H. Faber, R. Rackwitz, and P. Thoft-Christensen, “Modeling in optimal inspection and repair,” in Proceedings of the 10th International Conference on Offshore Mechanics and Arctic Engineering (OMAE '91), pp. 281–288, Stavenger, Noruega, 1991.

P. Thoft-Christensen and J. D. Sørensen, “Optimal strategies for inspection and repair of structural systems,” Civil Engineering Systems, vol. 4, no. 2, pp. 94–100, 1987.

T. Moan, O. T. Vardal, N.-C. Hellevig, and K. Skjoldli, “In-service observations of cracks in North Sea jackets. A study on initial crack depth and POD values,” in Proceedings of the 16th International Conference on Offshore Mechanics and Arctic Engineering (OMAE '97), pp. 189–197, Yokohama, Japan, April 1997.

D. Straub, J. D. Sørensen, J. Goyet, and M. H. Faber, “Benefits of risk based inspection planning for offshore structures,” in Proceedings of the 25th International Conference on Offshore Mechanics and Arctic Engineering (OMAE '06), pp. 1–10, Hamburg, Germany, June 2006.

M. H. Faber, “Reliability based inspection planning of fatigue damaged offshore platforms,” in Proceedings of the 1st International Symposium on Analysis of Structural Risk and Reliability for Offshore Facilities (IMP '01), Cancun, Mexico, 2001.

C. E. O. Estrada and D. de León Escobedo, “Development of a cost-benefit model for inspection of offshore jacket structures in Mexico,” in Proceedings of the 22nd International Conference on Offshore Mechanics and Arctic Engineering; Offshore Technology Ocean Space Utilization, pp. 87–98, Cancun, Mexico, June 2003.

J. D. Sørensen and G. Ersdal, “Risk based inspection planning of ageing structures,” in Proceedings of the 27th International Conference on Offshore Mechanics and Arctic Engineering (OMAE '08), pp. 399–408, Estoril, Portugal, June 2008.

C. A. Cornell, F. Jalayer, R. O. Hamburger, and D. A. Foutch, “Probabilistic basis for 2000 SAC federal emergency management agency steel moment frame guidelines,” Journal of Structural Engineering, vol. 128, no. 4, pp. 526–533, 2002.

C. A. Cornell, “Calculating building seismic performance reliability: a basis for multi-level design norms,” in Proceedings of the 11th World Conference on Earthquake Engineering, 1996, paper no. 2122.

Federal Emergency Management Agency (FEMA), “Recommended seismic evaluation and upgrade criteria for existing welded steel moment-frame buildings,” Tech. Rep. FEMA 351, SAC Joint Venture, Washington, DC, USA, 2000.

S.-Y. Yun, R. O. Hamburger, C. A. Cornell, and D. A. Foutch, “Seismic performance evaluation for steel moment frames,” Journal of Structural Engineering, vol. 128, no. 4, pp. 534–545, 2002.

M. A. Torres and S. E. Ruiz, “Structural reliability evaluation considering capacity degradation over time,” Engineering Structures, vol. 29, no. 9, pp. 2183–2192, 2007.

D. Celarec, D. Vamvatsikos, and M. Dolšek, “Simplified estimation of seismic risk for reinforced concrete buildings with consideration of corrosion over time,” Bulletin of Earthquake Engineering, vol. 9, no. 4, pp. 1137–1155, 2011.

D. Vamvatsikos and M. Dolšek, “Equivalent constant rates for performance-based seismic assessment of ageing structures,” Structural Safety, vol. 33, no. 1, pp. 8–18, 2011.

D. Tolentino, S. E. Ruiz, and M. A. Torres, “Simplified closed-form expressions for the mean failure rate of structures considering structural deterioration,” Structure and Infrastructure Engineering, vol. 8, no. 5, pp. 483–496, 2012.

D. Tolentino, Optimización multiobjetivo para el mantenimiento de estructuras considerando la influencia del daño acumulado [Ph.D. thesis], Graduate Program in Engineering, Universidad Nacional Autónoma de México, 2013 (Spanish).

D. E. Goldberg, Genetic Algorithms in Search, Optimization and Machine Learning, Addison Wesley Longman Inc., Boston, Mass, USA, 1989.

M. Fragiadakis, N. D. Lagaros, and M. Papadrakakis, “Performance-based multiobjective optimum design of steel structures considering life-cycle cost,” Structural and Multidisciplinary Optimization, vol. 32, no. 1, pp. 1–11, 2006.

O. Lavan and G. F. Dargush, “Multi-objective evolutionary seismic design with passive energy dissipation systems,” Journal of Earthquake Engineering, vol. 13, no. 6, pp. 758–790, 2009.

K. Deb, Multi-Objective Optimization Using Evolutionary Algorithms, John Wiley & Sons, Chichester, UK, 2001.

D. E. Goldberg and K. Deb, “A comparative analysis of selection schemes used in genetic algorithms,” in Foundations of Genetic Algorithms (Bloomington, IN, 1990), G. J. E. Rawlins, Ed., pp. 69–93, Morgan Kaufmann, San Mateo, Calif, USA, 1991.

H. M. Mühlenbein and D. Schlierkamp-Voosen, “Predictive models for the breeder genetic algorithm,” Evolutionary Computation, vol. 1, no. 1, pp. 25–49, 1993.

The MathWorks, Global Optimization Toolbox 3 User’s Guide, Version 3.1, 2010.

F. L. Silva and E. Heredia, “Effect of uncertainties on the reliability of fatigue damaged systems,” in Proceedings of the 23th International Conference on Offshore Mechanics and Arctic Engineering (OMAE '04), vol. 2, pp. 427–434, Vancouver, Canada, 2004.

P. Paris and F. A. Erdogan, “Critical analysis of crack propagation laws,” Journal of Basic Engineering, vol. 85, pp. 258–534, 1963.

K. Sobczyk and F. B. Spencer, Random Fatigue: From Data Theory, Academic Press, Boston, Mass, USA, 1992.

A. Stacey, J. V. Sharp, and N. W. Nichols, “Static strength assessment of cracked tubular joints,” in Proceedings of the 15th International Conference on Offshore Mechanics and Arctic Engineering (OMAE '96), vol. 3, pp. 211–224, June 1996.

F. M. Burdekin, “The static strength of cracked joints in tubular members,” Offshore Technology Report 2001/080, Health and Safety Executive, London, UK, 2002.

American Petroleum Institute (API), Recommended Practice for Planning Designing and Construction of Fixed Offshore Platforms Load and Resistance Factor Design, API RP 2a-LRFD, Washington, DC, USA, 1993.

PEMEX, “Design and evaluation of fixed marine platforms in Campeche Bay,” NRF-2003-Pemex-2000, Mexico, 2000

A. Raine, “The development of Alternating Current Field Measurement (ACFM) technology as a technique for the detection of surface breaking defects in conducting material and its use in commercial and industrial applications,” in Proceedings of the 15th World Conference on Non-Destructive Testing, Rome, Italy, 2000.

The Marine Technology Directorate Limited, Review of Repairs to Offshore Structures and Pipelines, vol. 94–102, 1994.

L. R. Ford, Differential Equations, McGraw-Hill Book Company, New York, NY, USA, 1955.

E. D. Rainville, Intermediate Course in Differential Equations, John Wiley & Sons, New York, NY, USA, 1961.