Abstract

This paper addresses the importance of engineering asset management decisions and control. For this purpose, a Life-Cycle Cost (LCC) analysis is conducted for typical reinforced concrete (R/C) buildings located in Mexico City. The objective of this study is to develop an artificial neural network (ANN) model that can estimate the total expected cost of R/C buildings by using LCC functions. The total cost includes the initial cost and the cost of the damage caused by future possible ground motions at the site of interest. The present value of the cost includes: initial cost, repair or reconstruction cost, cost of damage to the contents, costs associated with the loss of life or injuries and economic losses. The structural performance is evaluated using probabilistic models, artificial neural networks models are used to obtain the seismic response of the buildings. The methodology is applied to a set of reinforced concrete buildings with 4, 8, and 12 stories which are located at the soft soil of Mexico City. Finally, it is concluded that the life-cycle cost is efficiently obtained using artificial neural network models for estimating the structural reliability of reinforced concrete buildings, in such a way that it can be used as an excellent planning tool that covers long spans of time.

1. Introduction

The life-cycle cost estimation has received abundant attention over the past decades. In that time, the challenge of estimating the total life-cycle costs of structures considering all variables involved in the problem has represented a difficult task. There is a vast amount of literature available regarding the estimation of the expected cost for different structural systems [119]; however, those studies generally are applied to a limited number of particular cases. Furthermore, one of the main limitations of those studies is the time consuming of the methodology to obtain the total expected cost. This study addresses this issue, to perform this task artificial neural networks are used aiming to minimize the time consuming of the methodology within low errors in the estimation. The great potential of ANNs is the high-speed processing provided in a massive parallel implementation [20]. ANNs can be developed and used for image recognition, natural language processing and so on. Nowadays, ANNs are vastly used for universal function approximation in numerical paradigms because of their excellent properties of self-learning, adaptivity, fault tolerance, nonlinearity, and advancement in input and output mapping [21]. In the engineering field ANN has been used for optimization and seismic code calibration [22], they have been applied for estimating concrete and reinforcement consumption in the construction of integral bridges [23]; moreover, an artificial neural network model has been applied to estimate the construction costs [24]. For bridge design an ANN was used to decrease the computational demand of box-girder element analysis [25]. The prediction and identification of seismic-induced damage in structures has been possible using neural networks [2628]. For earthquake engineering, the use of the power of ANN pattern recognition was applied to evaluate seismic risk problems [29].

The estimation of the life-cycle cost of buildings is a primary part of a construction project, it is considered one of the major criteria in building design for its importance in helping to choose economic structural configurations and to estimate future costs of ownership. In this paper, it is considered that the total expected life-cycle cost of buildings includes: the initial costs plus the expected costs of the damage caused to the structure by future earthquakes, including repair cost, cost of damage to the contents, the cost associated with the loss of life and injuries, and direct economic losses [30]. The loads considered for the structural design are: dead load, life load, and earthquake loads. The structural performance is evaluated using probabilistic analysis. The occurrence of earthquakes is described by a Poisson process [31]. The methodology is applied to a set of reinforced concrete buildings with 4, 8, and 12 stories (low, middle and high rise) located in soft soil with ground period of Ts = 2.0s of Mexico City, which are of main concern in Mexico City Building Code (MCBC). Finally, it is concluded that artificial intelligence can provide accurate results, minimize prediction errors, and it can be used as an excellent planning tool that covers long spans of time.

2. General Methodology

The general approach taken to obtain the life-cycle costs of buildings is summarized in the steps below and illustrated in Figure 1.

2.1. Structural Design of Buildings (Step 1)

The buildings are designed according to the Mexico City Building Code considering the seismic design and concrete regulations for buildings located in the soft soil of Mexico City.

2.2. Maximum Structural Capacity (Step 2)

The maximum structural capacity is obtained using incremental dynamic analysis (IDA) [32]. For this aim, the maximum inter-story drift (MID) is selected as the engineering demand parameter to conduct the analysis. The associated yielding and near collapse limit states of the buildings can be obtained from this analysis.

2.3. Structural Reliability (Step 3)

The structural reliability is obtained from the seismic demand hazard curve that estimates the mean of the annual rate of exceeding a certain MID. The curve is calculated by the following expression [33,34]:Where : represents the average number of times per year that MID exceed a given value of d; d: a given value of MID; D: structural demand, presented by the maximum inter-story drift; Sa: pseudo-acceleration associated with the fundamental period of vibration of the building; : the seismic hazard curve for the site of interest; andwhere : is the fragility curve; : the standard normal cumulative distribution function; : the median of the logarithmic value of the seismic demand; : logarithmic standard deviation of the demand.

2.4. Artificial Neural Network Training (Step 4)

In this study, the artificial neural network technique has been used for the structural design of the buildings and to estimate the capacity and seismic demand hazard curves. From the previous steps, and in order to obtain the structural capacity and seismic demand hazard curve, numerous nonlinear analyses that are time consuming are required. Therefore, it is possible to reduce the time for the analysis of the building using ANN models. The Matlab [35] software was used in this study to develop the ANNs.

2.5. Seismic Demands Simulation (Step 5)

Different values of seismic demands are simulated from the seismic demand hazard curve (equation (1)). For the numerical simulation, the average number of times per year that a seismic event occurs with a magnitude equal or greater than 6 is considered ν = 3. The events are modeled by a Poisson process; therefore, the earthquakes follow an exponential distribution. By using the inverse transformation method [36] different values of MID are simulated for the life-cycle of the buildings.

2.6. Damage Index (Step 6)

There are different approaches to estimate the structural damage [37,38]. In the present study, the level of structural degradation is estimated from a measure of physical damage, which is represented as a damage index, DI, expressed as the ratio between the global structural capacity of the building and structural demand. The values taken by the DI are between 0 and 1, when DI equals 0 represents no damage in the structure, and DI equals 1 represents the total damage. The DI is defined by equation (3) [39]:Where : is the maximum inter-story drift demand, which is obtained from the simulation of the seismic demands of the seismic demand hazard curve (step 5); : is the maximum inter-story drift associated with serviceability limit state (structure without damage), which is obtained from the statistic value of the limit state considered in the incremental dynamic analyses; : is the maximum inter-story drift associated with incipient collapse.

2.7. Life Cycle Cost Functions (Step 7)

The total cost of a structure with an expected life cycle of 50 years is the combination of the initial cost plus the expected damage cost:Where represents the total cost of the structure; : is the initial cost; is the cost associated with the damage.

2.7.1. Initial Cost

The initial cost is calculated using an approximated procedure proposed by De Leon [14]. The cost includes direct cost, indirect costs and the utility paid to the contractor. The direct cost is estimated from the material cost plus the cost of labor, approximately 40% of , . The indirect cost due to the nonstructural work the contractor has (i.e.: insurance, training, office expenses, etc.) is estimated to be 20% of , . The constructor fee is 15% of the summation of the direct cost and indirect cost, . The total initial cost can be estimated as the summation of and , equation (5);Where; : is the initial cost; : is the direct cost; : is the indirect cost; : is the utility.

The unit cost of the construction materials is shown in Table 1. For the cost analyses the average cost of concrete and steel in Mexico City is used.

2.7.2. Damage Cost

The damage costs of the structure during its useful life can be express as the following: repair or reconstruction, ; loss of content, ; economic losses, ; injury, ; and loss of life, . The damage cost can be expressed as Eq. (6);

The damage costs are obtained by means of the simulated level of structural damage present in the building. The damage index (Step 6) is used to measure the structural damage simulated.

(1) Repair Cost or Reconstruction. There are several methods to repair buildings damaged by previous earthquakes, depending on factors such as the level of damage, type of structure, configuration, logistics. In this study, it is assumed that the R/C buildings will be repaired using the jacketing technique of structural members, this technique allows the structure to be repaired effectively and return to practically its initial condition before the damage. Based on the repair costs of buildings that were damaged by previous earthquakes a relationship between the repair cost and the damage index was established. The structural damage in some cases is extremely severe such structures can no longer be repaired, and it is necessary to be demolished. De León and Ang [40] considered that for DI > 0.70 the R/C structure needs to be demolished. The costs of reconstruction are assumed to be equal to 1.2 , which includes costs of demolition, cleaning and redesigning the structure. For DI < 0.70 the repair cost is a function of the initial cost times the DI to the second power. The repair or reconstruction costs are defined by equation (7) and (8):

The relationship between the normalize repair cost and the damage index is shown in Figure 2.

(2) Cost Due to Damage of Contents. The estimation of costs due to damage of contents depends on the type and use of the structure; for example, the damage of contents in a school does not have the same impact as that of an office building or shopping mall. In this study, it is assumed that the buildings are offices, it is considered as suggested by Suranhman and Rojiani [41], that the maximum content loss if DI > 1 is equal to 50% of the initial cost. For values of DI < 1, the expected cost is assumed to be a linear function of DI. The cost due to loss of content is defined by equations (9) and (10):

The relationship between the normalized cost due to damage of contents and the damage index is shown in Figure 3.

(3) Cost Due to Economic Loss. The cost due to economic loss depends on the economic activity for which the building is used. This study assumed the usage of the building as offices; therefore, the economic loss is associated with the loss of income due to rental during the time of repair or reconstruction.

The maximum cost for economic loss (DI > 1) is a function of the period of reconstruction (PR), the average rent in dollars per square meter (R) and the build area in meters of the structure (A). For values of DI < 1 the variation of cost is assumed to be a function of DI to the second power [14]. The average monthly rental for offices in Mexico City is equal to $19 US dollars/m2 and it is assumed that the maximum period of reconstruction of a building is equal to 24 months. The function due to economic loss is defined by equations (11) and (12).

The relationship between the normalized economic cost per square meter of building area and the damage index is shown in Figure 4.

(4) Cost Due to Loss of Life. The estimation of the cost associated with life loss is a difficult task because it is a sensitive issue [4245]. Two postures regarding the value of a person life can be found. The first one, human life has a statistical value based on its income; and the second one, the cost of human life is invaluable. For this work, the value of a life is considered by the average income of a person.

For the present study, the cost due to the loss of life is estimated relying on previous catastrophic seismic events. For this reason, it is necessary to estimate the average number of dead people inside a building during a seismic event. A nonlinear regression analysis [40] was obtained based on the total area of buildings that collapse during the 1985 Mexico City earthquake [46] and the number of deaths [47]. The result of this analysis is summarized in equation (13) and shown in Figure 5:Where : is the number of deaths; A: is the of collapsed buildings area in 1000 m2.

The cost due to the loss of life function (DI > 1) is expressed as the total number of deaths, , in the case of incipient collapse, multiplied by the value of the statistical life-time income of a person, . For smaller vales of DI, it is assumed that the cost function can be represented by the total number of deaths, multiplied by the statistical life-time income of a person multiplied by DI to the fourth power, Equation (15). It is considered that the average annual income of a person in Mexico is considered equal to $10000 US dollars [48], and that the useful working life per person is equal to 25 years, so the cost per death of a person is equal to $250,000 US dollars. The cost functions due to loss of life are summarized in equations (14) and (15) [14]:

(5) Cost Due to Injuries. The evaluation of cost due to injuries, , refers to expenses required for hospitalization of people who became injured during an earthquake [14]. The average number of people injured per unit area of collapse is equal to 0.0168/m2 estimated as the relation of the number of injuries reported in the 1985 Mexico City earthquake [47] with the total area of buildings collapsed [46]. The cost of injuries without disability, CLS, is considered equal to $2000 US dollars. Notice that costs for minor injuries, medical expenses and medicine, including a small stay in the hospital are considered. Injured without disability are considered to represent 90% of all injured people [40]. The cost of injuries resulting in disability, , is considered equal to $250,000 US dollars (equal to the cost of death). People with disabilities represent 10% of total injuries.

For DI < 1 the cost due to injuries is assumed to be a function of the total people injured in the building times the DI raised to the second power, Equation (16) [14]. For DI > 1 the cost function is assumed to be a function of the total people injured in the building with disability and without disability in the building, Equation (17).

The relationship between the normalized cost due to injuries per person injured and the damage index is shown in Figure 6.

3. Illustrative Example

The illustration of the methodology is presented by an example in this section. The methodology is applied to a 12-story, 3-bays reinforced concrete building.

3.1. Selected Ground Motion Records

In this study three different buildings are located in soft soil with a dominant period TS = 2.0s. A set of 15 ground motion records were selected to perform the analyses of the buildings. Table 2 shows the main characteristics of the earthquake ground motions. The ground motion records were obtained from the IINGEN database [49]. The pseudo-acceleration elastic seismic response spectra of the records selected are shown in Figure 7 considering 5% of critical damping.

Figure 8 shows the seismic hazard curves associated with different fundamental structural periods. Those curves represent the probability that a structure will exceed a specified seismic intensity (e.g., pseudo-acceleration) in one year. The curves are used to solve Eq. 1, and corresponds to a site located at soft soil of Mexico City.

3.2. Reinforced Concrete Building Characteristics

Three types of reinforced concrete buildings, low, mid, and high rise were selected for the seismic analyses and design according with common practice in Mexico City. The buildings were represented by three-dimensional structural models with elevation between 4, 8, and 12 stories, an inter-story height equal to 4 meters. The buildings are constituted by orthogonal rigid reinforced concrete frames constituted by beams and columns, with three bays in both directions. While for the 4-story building the bay length is equal to 6m, for the 8 and 12-story buildings the bay length is equal to 8m. The frames are connected by floor slabs assumed to be rigid diaphragms. Figure 9 shows the geometric characteristics of the buildings.

3.3. Structural Design

The 12-story building is designed in accordance with the Mexico City Building Code [50]. The loads considered for the analysis are; dead load, live load and seismic load.

The 12- story building was divided into 3 groups of beams and 3 groups of columns to uniform the buildings sections. Table 3 shows the results of the structural design for the 12-story building. The sections characteristics shown in Table 3 are the parameters for beams; b: width, h: depth, Ase: longitudinal reinforcement inferior at the end region of beam, Ase’: longitudinal reinforcement superior at the end region of beam, Asm: longitudinal reinforcement inferior at the mid span, Asm’: longitudinal reinforcement superior at the mid span, se: stirrups spacing at the end of the beam, sm; stirrups spacing in the middle of the beam, for columns; b: width, h: depth, As: longitudinal reinforcement, s: stirrups spacing.

3.4. Building Structural Capacity

The structural capacity of the building is calculated using incremental dynamic analysis. For this aim, the software Ruaumoko3D [51] was used to apply a series of nonlinear structural analyses to the building in order to determine the structural behavior. Nonlinear structural analysis is performed by subjecting the structure to a set of records (Table 2) scaled at different increasing levels of pseudo-acceleration. The result extracted from each analysis is a discrete point used to form the IDA curve (see Figure 10). The modified Takeda hysteretic model was considered in the analyses to assume the degradation of stiffness and strength of the reinforced concrete elements.

The limit state near collapse is defined within the Ruaumoko computer program when one of the following conditions are reached: (a) Numerical instability, (b) decrement of the tangent lateral stiffness to 20% of the initial value or, c) ultimate rotation condition of a member. The yielding limit state was defined as the difference of 10% of the elastic slope.

3.5. Seismic Demand Hazard Curve

In order to obtain the seismic demand hazard curves, the fragility curves for the 12-story building are needed. The fragility curves represent the conditional probability of damage being exceeded a particular damage state, these curves are obtained from Eq. 2.

Figure 11 presents a set of 12 fragility curves. Those curves show the vulnerability of the building underground motions, it represents the conditional probability of exceeding a maximum inter-story drift (d) as stated in the Mexico City Building Code [50]. Two limit states are considered, the service limit state and near collapse state. This figure illustrates that when a ground motion intensity Sa = 900 cm/s2, the probability of reaching or exceeding d = 0.015 is approximately 30% for the 12-story building. From the figure it can be seen that for the same intensity Sa, the exceedance probability (P) decreases as the peak inter-story drift (d) value increases. The seismic demand hazard curve show in Figure 12 was obtained from Eq. 1, and it represents the mean annual rate of exceeding a maximum inter-story drift (d).

3.6. Reliability Analysis Using Artificial Neural Networks

In this study, the technique of artificial neural networks is used, first; to design buildings based on similar characteristics of the buildings used in the database, second; to estimate the structural capacity and the seismic demand hazard curve. Two different ANN models were trained using data from a set of buildings designed in previous studies [5254].

For the training phase of the first model a feedforward backpropagation network was used to design a 12-story building. The architectural model of the ANN used was an input layer, hidden layers and the output layer, with 7 neurons for the input layer, 5 neurons for the hidden layer and 36 neurons for the output. The activation functions were hyperbolic tangent sigmoid (tansig) for the hidden and output layer. A set of 200 buildings were used as data for the training phase of the first ANN model. The 7 neurons in the input layer represents the general characteristics of the buildings; number of bays in X direction, number of bays in Y direction, story levels, story height, bays spacing X, bays spacing Y, soil period. The 36 output neurons represent the parameters of the beams and columns designed.

The results of the design obtained from the first ANN model are shown in Table 4, which are compared with the results obtained using traditional approaches for seismic design (see Table 3). The results showed excellent similarity with the actual values, attaining a good accuracy prediction with a relative error less than 9%. Moreover, the selected 12-story building was not used in the dataset of the training phase.

The second ANN model estimates the seismic demand hazard curve and capacity of the 12-story building. The input layer with 44 neurons, two hidden layers with 5 neurons and the output layer with 14 neurons, constitutes the architectural model used in the second ANN tool. The same set of 200 buildings was used in the training phase. The activation functions were hyperbolic tangent sigmoid (tansig) for the two hidden layers, and output layer. The input layer parameters are the buildings general and section characteristics and the output layer represents the capacity and seismic demand hazard curve. Table 5 summarizes the input vector parameters.

Figure 13 shows the results of the ANN model and the comparison to those obtained using equation (1). From the figure, it can be seen that the estimation of the seismic demand hazard curve simulated with the ANN; SIM curve, attained good accuracy on predicting equation (1), with a mean square error (MSE) of 2.42e-08 about 9% relative error.

3.7. Seismic Demand Simulation

The simulation of the seismic demands considering 50 years of service life is performed using the seismic demand hazard curve estimated by the ANN. For the simulation, the average number of seismic events per year is needed, in this study an average number of 3 seismic events with magnitude equal or greater than 6 was used [22]. Figure 14 shows the seismic demands in terms of MID, considering the service life for the 12-story building.

3.8. Damage Index in Service Life of the Building

The damages are evaluated using equation (3). The damage index is then obtained for every seismic demand simulated for the service life of the building. The limit state associated with serviceability, , and the limit state associated with incipient collapse,  = 0.0356, were obtained from the results of the ANN training. Figure 15 shows the damage index for the simulated demands scenario of Figure 14.

3.9. Cost of the 12 Story Building

The initial cost of the 12-story building is obtained from the amount of material (concrete and steel) necessary for its construction. The material is obtained from the quantification of the material necessary to build the beam and column elements obtained from the structural design, 557 tons of reinforced steel and 3660 cubic meters of concrete were quantified. Equation (5) is applied to obtain the initial costs.

The damage costs are obtained for the simulated seismic demand of the year 14.07, where the demand is 0.014 for the 12-story building. Applying equation (3) for the year 14.07, the damage index is 0.359 (Figure 15). The costs due to damage are then:

Repair cost or reconstruction (Equation (7)):

Cost due to damage of contents (Equation (9)):

Cost due to economic loss (Equation (12)):

Cost due to loss of life (Equations (13) and (14)):

Cost due to injuries (Equation (16))

The total costs due to damage for the year 14.07 was obtained by the summation of all the costs due to damage (Equation (6)).

4. Results and Discussion

This section shows the results of applying the methodology for the 4, 8 and 12 story buildings.

The ANN prediction results of the seismic demand hazard curve are compared to those of Equation (1). Figure 16 shows the predicted values of the seismic demand hazard curve for the 4-story building (SIM curve), based on the results of the predictions done by the ANN, the curve shows a deviation from the curve of Eq.1. The mean square error 2.97e-07 is estimated from the deviation, about a 12% of relative error; therefore, the usefulness of the seismic demand hazard curve to simulate the seismic demand during the life-cycle of the building is sufficient to predict accurate demand values. The seismic demand hazard curve for the 8-story building is shown in Figure 17, where the mean square error estimated from the deviation of the curve predicted by the ANN is 3.81Ee-07, about a 6% of relative error.

The initial cost of the 4-story building is obtained by applying Equation (5). The amount of material necessary for the construction was obtained from the quantification of the material necessary to build the structural members given by the structural design of the ANN, a total of 39.5 tons of steel and 353.75 cubic meters of concrete were obtained. By applying Equation (5) and multiplying by unit cost of the material (Table 1) the initial cost US dollars.

Table 6 shows the cost due to damage expected over the life cycle of the 4-story building. From the table it can be seen the damage scenario simulated throughout the 50 year life-cycle of the building, and the damage cost associated with that damage. The total life-cycle cost for the 4-story building is the sum of the initial cost plus the cost due to damage US dollars.

To estimate the initial cost for the 8-story building, the ANN was used to design the building and quantify the material necessary for its construction. A total of 230 tons of steel and 1905 cubic meters of concrete were quantified. The initial cost is obtained applying equation (5), the cost is US dollars. The damage cost for the simulated structural demand scenario of the life-cycle of the 8-story building is shown in Table 7. The life-cycle cost is obtained as sum of the initial cost plus the damage cost, resulting in .

The damage cost for the simulated structural demand scenario of the life-cycle of the 12-story building is shown in Table 8. The life-cycle cost of the 12-story structure is obtained from the sum of the costs associated with damage during the useful life of the structure and the initial cost, for this scenario of simulated demands, the total cost is: .

5. Conclusions

In this paper, a general methodology to obtain the life-cycle costs for buildings located at the soft soil of Mexico City was presented. The costs considered in the analysis are the initial costs and the cost due to damage. The structural damage was obtained by simulating the damage using the seismic demand hazard curve; however, the number of nonlinear analyses to obtain the curve can be complex and time demanding. This paper proposes a methodology to reduce the analysis time by using ANNs. The proposed method uses ANNs to model a relationship between the parameters describing the structure and the seismic demand hazard curve. Two feedforward backpropagation network models were selected, the first model was used to obtain the structural design of the buildings, and the second model to obtain the seismic demand hazard curve. The analyzes were conducted considering 3D R/C buildings with varying story height, and bay length subjected to soft soil ground motion records. The performance of the ANN was evaluated by obtaining the seismic demand hazard curve for the different buildings and it was found to be capable of successfully predicting the curves. The accuracy of the developed model was tested by comparing the results of the traditional method and those obtained by the ANNs. The results show that the relative error of the seismic demand hazard curve for the different buildings were between 6% and 12%. It can be considered that the results were acceptable, taking into consideration the computational time reduced from about a few hours when obtaining an exact solution to a few seconds when the ANN model is used. A limitation in the applicability of the presented method is due to the nature of the ANNs that learn the relationships using examples data sets, the particular implementation proposed is then limited to the range of parameters, such as types of buildings, types of soil, properties of materials in construction, etc. Used for the training of the ANNs, any extension of its applicability in other types of buildings will require the generation of additional data to expand the parameters used in the training of the network. Based on the results obtained from the ANN, the seismic demand hazard curves were used to simulate demands in the life-cycle of the structures, from the simulated demands the damage scenario for the life-cycle of the structures was obtained. The total costs of the structures that are given by the initial costs and costs for damage during the useful life were obtained. The total cost obtained in the analysis for the 4-story building was $323,586 US dollars, for the 8-story building the total cost obtained was $4.64 million US dollars, and for the 12-story building the total cost during its life-cycle was $4.91 million US dollars. The proposed approach to obtain the life-cycle cost would be of great use to investors and builders in helping them to consider the future costs of constructing buildings on soft ground in Mexico City. Finally, it is concluded that life-cycle cost is efficiently obtain using artificial neural networks models for the reliability analysis of reinforced concrete building, so it can be used as an excellent planning tool that covers long spans of time.

Data Availability

The data used to support the findings of this study are available from the corresponding authors upon request.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

The first author would like to thank CONACYT (Consejo Nacional de Ciencia y Tecnología from México), for the scholarship provided to pursue his PhD studies. The second and sixth author appreciate the support of CONACYT under grant Ciencia Básica 287103. Thanks are given to the PASPA program of DGAPA-UNAM for the support provided to the fourth author. Financial support also was received from the Universidad Autónoma de Sinaloa under grant PROFAPI 2022.