Next Article in Journal
In Search of Sustainability and Financial Returns: The Case of ESG Energy Funds
Next Article in Special Issue
Taxonomy and Indicators for ESG Investments
Previous Article in Journal
A Newly Acidophilic Bacterium Acidomyces acidothermus Was Isolated to Efficiently Bioleach Copper from Waste Printed Circuit Boards (WPCBs)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simplified Life Cycle Cost Estimation of Low-Rise Steel Buildings Using Fundamental Period

Department of Civil & Architectural Engineering, Sungkyunkwan University, Suwon 16419, Republic of Korea
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(3), 2706; https://doi.org/10.3390/su15032706
Submission received: 28 December 2022 / Revised: 30 January 2023 / Accepted: 30 January 2023 / Published: 2 February 2023

Abstract

:
In the current study, a simplified seismic life cycle cost (LCC) estimation procedure is proposed utilizing the mean values of the structure’s main input variables. The main input variables of the building are used for constructing a relationship between the structural fundamental period (T) and an average estimation of the LCC (LCCavg). Using the actual building properties related to damage probability, the T–LCCavg relationship is used to obtain the final LCC (LCCfin). The equivalent single degree of freedom (ESDOF) model and SAC-FEMA framework are utilized for damage probability calculation. The dispersion measure in demand is approximately calculated based on the mean plus one standard deviation of the seismic hazard response spectrum, and, then, verified through nonlinear time history (NLTH) analyses of the original structure. Five and three-story steel buildings are used as case studies for verification of the proposed method. The analysis results indicate that the proposed procedure provides reasonable LCC estimations for low-rise buildings dominated by the fundamental mode of vibration.

1. Introduction

The advancement in PBEE (performance-based earthquake engineering) in the last two decades has made great leaps in terms of quantifying and assessing the performance of buildings in a manner suitable for decision-makers. This can be attributed to the robust, yet complicated, probabilistic models, and the advancements in structural analysis and computational capabilities. On the other hand, the role of the structural engineer is challenged because of the complexity inherited in the process itself [1]. To this end, the need for simple and reliable methods suitable for loss assessment and LCC estimation is justified. Limited studies have been conducted on simplifying the LCC of building structures, which makes it an attractive field to investigate for researchers and engineers. Many of the recent studies in LCC are dedicated to the assessment and optimization of LCC. For example, Mahmoud and Cheng [2] provided a probabilistic framework for evaluating different design approaches when considering LCC. Two steel buildings with different heights subjected to different seismic and wind intensities were used as case studies. Cutfield et al. [3] presented a life cycle cost analysis of a conventional and base-isolated steel braced frame office building. The study uses the FEMA P-58 approach to assess the impact of moat wall pounding and business interruption. Matta [4] proposed a multi-hazard LCC approach as an efficient alternative for Cost-effective assessment of a TMD. Noureldin et al. [5] evaluated the cost-effectiveness of retrofitting steel structures using steel slit dampers and shape memory alloy-based hybrid dampers. Recently, Torti et al. [6] proposed an approach based on life cycle cost analysis to assess the effectiveness of equipping bridges with seismic structural health monitoring systems. Gidaris and Taflanidis [7] described a probabilistic, simulation-based framework for estimating life cycle cost and developed a stochastic search approach to support an efficient optimization under different design scenarios. Esteva et al. [8] provided an overview of the general framework supporting the optimum life cycle-based engineering decisions to examine the relationship between the target safety level and the required maintenance of the structure. Vitiello et al. [9] developed a semi-probabilistic approach to assess structures economically. The approach addresses the cost-effective retrofitting strategies of structures during their lifetime. Noureldin and Kim [10] proposed a parameterized LCC evaluation method of structures considering stiffness, strength, and ductility.
Risk assessment is considered an important step in any LCC methodology and simplifying this step will reduce the computational cost required for LCC calculations. Recently, some simplified methods have been proposed for risk assessment. For example, Fragiadakis and Vamvatsikos [11] used pushover analysis and approximate incremental dynamic analysis to propose a fast performance uncertainty-estimation method. Another important study was conducted by Dolsek [12], who proposed an approximate method for seismic risk assessment taking into account various uncertainties. The method was demonstrated using an example of a four-story reinforced concrete frame. Fajfar and Dolsek [13] proposed a practice-oriented probabilistic approach for the seismic performance assessment of building structures using the SAC-FEMA approach and the pushover-based N2 method. Sullivan et al. [14] discussed some simplified tools for risk assessment, such as the simple definition of collapse fragility, based on FEMA P-58 guidelines [15], and applied it to a benchmark 4-story building. Welch et al. [1] proposed a simplified probabilistic loss assessment methodology that builds on a direct displacement-based framework. The methodology was tested via an examination of two RC frame buildings. Kosič et al. [16] utilized a deterministic structural model with a probabilistic equivalent SDOF model to introduce an approximate seismic risk assessment framework. Kosič et al. [17] proposed dispersion values for a practice-oriented pushover-based method for the estimation of failure probability for eight selected examples. Noureldin and Kim [18] developed a simple methodology for seismic life cycle cost (LCC) estimation for a steel jacket offshore platform structure. In their study, a localized incremental dynamic analysis and an approximate fragility curve are used. In their study, Liu and Mi [19] proposed a new approach for analyzing the costs of energy-efficient buildings over their lifespan, called LCC analysis. The method probabilistically predicts the levels of seismic damage in energy-efficiency building features and includes the costs of repairing or replacing those features in the analysis. Recent studies in earthquake engineering [e.g., [20,21,22,23,24,25]] places an emphasis on the importance on the seismic assessment of different types of structures. This necessitates proposing new approaches of seismic assessment.
The objective of the current study is to estimate the life cycle cost of a building structure using its fundamental period as the main input parameter. The probability assessment proposed by the SAC-FEMA method in a closed-form is utilized. The proposed LCC procedure is used for estimating the seismic LCC of three- and five-story steel buildings. The accuracy of the proposed method is assessed by comparison with the results obtained from the conventional method based on NLTHA. Establishing a relationship between the fundamental structural period and the expected LCC will be useful both in new and existing structures. For new structures, during the preliminary structural design stage, where different lateral resisting systems are examined, it is required to estimate LCC quickly to finalize the structure system. The simplest parameter that can be used in this case is the building fundamental period (T). In this case, using an LCC–T relationship will be useful for judging the cost-effectiveness of each structural system. For existing structures, the natural period of the structure can be changed during its service life; for example, in the case of the installation of new non-structural elements, experiencing inelastic deformation, retrofitting of the existing structural system, or considering the effect of soil-structure interaction. Mwafy and Elnashai [26] showed that the natural period can be doubled when regular frames experience inelastic deformation. Soil-structure interaction may reduce the global system stiffness [27] of the building, which may cause period elongation of the superstructure. The change in the fundamental period of the structure will, consequently, lead to a change in the total seismic LCC. The proposed procedure will provide a simple and reliable estimation tool for the seismic life cycle cost of the structures that will be readily available for designers and decision-makers in the preliminary design of structures.
The current study addresses a crucial need in the field of LCC estimations of structures by focusing on simplification rather than assessment. This study stands out from previous studies by proposing a novel approach of using only one factor, the fundamental period, to estimate LCC approximately, instead of the multiple factors used in previous studies. Furthermore, this study establishes a graphical relationship between LCC and the fundamental period of the structure, providing a valuable tool for structural designers, particularly at the early stages of design where information about the structure is limited. This makes the current study a significant contribution to the field of LCC and sets it apart from previous studies on LCC.

2. Framework for Constructing LCC Curve

The proposed framework is based primarily on establishing a relationship between the fundamental structural period (T) and the expected LCC for a certain hazard area. This relationship will be configured as a curve relating the T and LCC (LCC curve). This LCC curve is similar to the response spectrum curve because it gives an average estimate of the LCC based on the average values of the main parameters affecting the LCC.

2.1. The Steps of the Proposed Procedure

Figure 1 shows the proposed procedure, which is summarized as follows:
  • Modal analysis of the original MDOF structure is conducted and the basic dynamic properties are extracted. Three different hazard levels corresponding to 75, 500, and 2500 years return periods are selected based on the site-specific hazard curve along with their response spectra. The inelastic responses of different SDOF systems with different periods (T = 0.5, 1.0, 1.5, 2.5, etc.) are estimated using a nonlinear static technique such as the capacity spectrum method [28] or the N2 method [29,30].
  • Damage state probability is computed to estimate the LCC. Cornell et al. [31] provide the following equation to compute the damage state probability P L s :
P L s = H S a c ^ e x p 1 2 k 2 b 2   β D s a 2 + β C 2
where S a c ^ is median drift capacity spectral acceleration (which can be obtained from the approximate fragility curve); H S a c ^ is the probability that the intensity in a site will exceed or equal to S a c ^ on an annual basis; k is the slope of the hazard curve on a log–log plot; b is the slope of the demand–intensity relationship on a log–log plot; βD|s is the standard deviation of the natural logarithm of the drift demands for a given Sa; and βc is the standard deviation of natural logarithm for the median capacity of drift C (assumed 0.3 considering previous studies, e.g., [31]).
After obtaining the probabilities of damage states, the estimated LCC can be given as [32]
E C L c = C o + 0 L E C S D 1 1 + λ t d t = C o + L   E C S D
where Co is the initial cost of the construction, L is the expected life span (service life) of the building, λ is the rate of return (or annual discount rate), and E[CSD] is the damage cost of the building on an annual basis. A and E[CSD] can be formulated as:
α = 1 e x p q L q L
q = l n 1 + λ
C S D = i = 1 N C i P i
where N represents the considered limit states, Pi is the ith damage state probability of the building during its service life, and Ci is the ith damage state cost including damage and its repair. Three different limit states are used, which are IO, LS, and CP [33]. Ci is taken as a fraction of the initial cost (0.3, 0.7, and 1.0, respectively, for IO, LS, and CP limit states; [34,35]). Pi is given by:
P i = P Δ D > Δ C , i   P Δ D > Δ C , i + 1  
where ΔD is the maximum inter-story drift that a building experiences under an earthquake excitation (demand), and ΔC,i is the maximum inter-story drift that corresponds to the capacity of the structure at a the ith damage state (capacity). The probability of demand is greater than capacity, Δ D > Δ C , i , is evaluated as discussed in the previous step.
After that, an average LCC graph is constructed, which relates T (the fundamental period of the structure) with the corresponding LCC. T represents the stiffness of the structure and is usually used for the design or retrofit of structures [36]. The LCC graph is based on calculating the mean inter-story drift ratio (MIDR) of the multi-degree-of-freedom (MDOF) system using the building characteristics (structural variables) such as height, modal participation factor, and modal amplitude at the roof. All the parameters are adjusted to the most common mean values of building structures. Each point in the graph represents the average LCC for a particular fundamental period. This graph can be adjusted for any other building characteristics using a simple algorithm or electronic spreadsheet tool. The following equation is used for adjusting the MIDR of the equivalent single degree of freedom (ESDOF) system:
D M D O F = s S D O F i n e l Γ . φ r H
where Γ is the modal participation factor; φ r is the modal amplitude at the roof of the fundamental mode shape; and H is the height of the structure.
Some parameters such as drift, seismic hazard regression coefficients, and damping ratio, have a lognormal distribution considering previous studies [37,38]. However, other parameters such as initial cost and limit state repair costs have ranges and coefficients of variations mentioned in previous studies [37,39,40], but no specific statistical distribution is found in literature for them.
A probabilistic relationship between the response of the MDOF system and its corresponding ESDOF system needs to be established to quantify the uncertainty in obtaining the system response using the SDOF system. This probabilistic relationship is investigated by Jeong and Elnashai [41] for different structural systems and can be expressed as
D ¯ m a x M D O F = 1.13   D ¯ m a x E S D O F
where D ¯ m a x M D O F and D ¯ m a x E S D O F are the mean maximum seismic drift demand of the MDOF and ESDOF systems, respectively.
Next, similar graphs (T–LCC) for the main cost and damage probability variables, e.g., dispersion in demand, dispersion in capacity, hazard-curve slope factor, service time, etc., are constructed. In each graph, the mean, the upper bound (mean plus one standard deviation) and lower bound (mean minus one standard deviation) of a variable are considered, while the other input variables remain at their mean values. Some variables can be assumed based on previous studies such as b and βc. For moment frames, b and βc can be assumed as 1.0 and 0.3, respectively [31], and k and βD|s can be obtained from the hazard curve and the hazard response spectra, respectively. Based on these parameter values, the average LCC (LCCavg) is estimated using the fundamental period of the structure and the average LCC graph. Similarly, from the (T–LCC) curves of the main variables, the difference in LCC between the LCCavg and the corresponding LCC of the value of the variable can be obtained at the same T. If the LCC is larger than LCCavg, the difference is added to the LCCavg.. If the difference is smaller than LCCavg, then it is subtracted from LCCavg. This difference will be used for estimating the final LCC of the building.
After obtaining the difference in LCC (LCCdif) from the (T–LCC) graphs of the main variables, the final estimate of the LCC is determined by
LCCfinal = LCCavg ± LCCdif.

2.2. Median Drift Demand and the Drift Capacity Intensity

The relationship between the drift demand (D) and the hazard intensity (Sa) can be obtained accurately using NLTHA. Based on that, the responses of the building are obtained for different levels of intensity measured and a regression relationship is established to consider record-to-record variability. Using power-law distribution [42], the relationship can be expressed as,
D ^ = a ( S a ) b
where D   ^ is the median drift demand; S a is the elastic spectral acceleration (a measure of ground motion intensity); and a and b are the regression coefficients for linear regression of drift demand D at the intensity S a in logarithmic space. The most accurate method to determine the regression coefficients a and b for an ESDOF system is by conducting a nonlinear dynamic time history response analysis (NLTHA) using all available ground motion records. This analysis should be used to plot the relationship between the drift demand and the intensity, and, then, a regression analysis should be performed on the natural logarithm of the drift demand versus the natural logarithm of the intensity. A simpler alternative may be used such that the mean response spectrum for each set of records is used to calculate the drift demand by using a non-linear static approach (e.g., the N2 method or the capacity spectrum method) based on the equivalent SDOF system. This means that three drift demands will be available to draw the relationship between the drift demand (D) and hazard intensity (Sa) in logarithmic space. From this relationship, a and b are obtained. It needs to be mentioned that two conditions are required to guarantee the accuracy of the procedure at this stage. The first is having an elastic-perfectly plastic force-deformation relationship for the equivalent SDOF system (i.e., bilinear idealization of the capacity curve of the original MDOF system). The second is the applicability of the equal displacement rule (i.e., the seismic demand of the inelastic SDOF system is equal to the seismic demand of the corresponding elastic SDOF system with the same period). It is worth mentioning that the precision of determining regression coefficients by utilizing the mean response spectrum based on a nonlinear static method may be low due to the following reasons:
  • Limited applicability: Nonlinear static analysis is typically only applied to regular structures and may not be suitable for irregular or complex structures;
  • Simplistic input motion: Nonlinear static analysis uses a simplified input motion, typically a lateral load pattern that simulates the expected ground motion. This simplification may not accurately represent the actual ground motion that the structure will experience during an earthquake;
  • The actual drift of the structure during an earthquake may be different from the performance-based drift obtained from non-linear static analysis, as the actual ground motion may not be as severe as the load pattern used in the analysis, and the structure may not reach its maximum capacity. Additionally, the actual ground motion may be more complex than the load pattern used in the analysis, and the structure may experience multiple load cycles.
The median drift capacity intensity S a c ^ corresponding to each limit state can be calculated from the relationship between the drift demand (D) and the hazard intensity (Sa), explained above. Using the assumed limit state values (for IO, LS, and CP) and the slope of the (D–Sa) relationship, the corresponding median drift capacity intensity S a c ^ can be obtained for each limit state. Assuming that the hazard curve is linear on a log–log space, the relationship between s a and H s a can be approximated by the following form [31]:
H s a = P S a s a = k o s a k
where s a is the elastic spectral acceleration (the intensity measure); H s a is the corresponding annual probability of exceeding s a , k o and k are the linear regression coefficients on a log–log graph.

2.3. Dispersion in the Seismic Demand

In the proposed method, two alternatives are proposed for calculating the dispersion in demand β D | s a . One is based on generating the elastic mean and the mean plus one standard deviation response spectra, which can be generated from the smooth elastic response spectrum for the case of a firm soil profile using amplification factors recommended by Newmark and Hall [43] and illustrated by Chopra [44]. Using the N2 method [18,45] and the elastic spectral accelerations from both the elastic mean and the mean plus one standard deviation response spectra, the inelastic responses of the SDOF system can be calculated. The difference between the two responses indicates the standard deviation in demand, which is βD|Sa. This alternative is used in the proposed method because of its simplicity, where only the natural period of the structure is required. The second alternative is based on conducting a response spectrum analysis on the MDOF structure using the elastic mean and mean plus one standard deviation response spectra. The difference between the two responses defines βD|Sa.
In the current study, a deterministic SDOF model is utilized and the other epistemic and aleatoric uncertainties are treated simply without using sampling techniques such as Monte Carlo simulations and Latin hypercube sampling. This adds an advantage to the proposed framework in the current study compared, for example, with other studies such as work from Kosič et al. [16,17], which requires many Monte Carlo simulations with Latin hypercube sampling at the level of the SDOF model. Another advantage of the proposed framework is through using the main input variables of the building and implementing them to adjust the mean LCC to obtain the final LCC. This is different than the study conducted by Fajfar and Dolsek [13], where all the variables required for the probabilistic loss estimation calculation are assumed.

3. Case Study Structures and Seismic Hazards

3.1. Design and Modelling of the Case Study Structures

To validate the proposed simplified LCC procedure, three- and five-story steel buildings were used as the case study structures. The steel structures were special moment resisting frames, having a typical story height of 4000 mm and a first story height of 5500 mm. A span of 6000 mm is considered for each bay in both directions. The configuration of the case-study models is shown in Figure 2 (one of the exterior frames was selected for the analysis). The gravity loads used for design were 4.1 kN/m2 and 2.5 kN/m2, respectively, for dead and live loads. The beams and columns of both structures were W-shaped sections. The steel yield stresses used for beams and columns were 250 and 345 MPa, respectively.
Plastic hinges used for conducting NLTHAs were introduced at the end of the columns and beams to account for the inelastic behavior of the members according to FEMA-356 [46]. Figure 3a shows the bending moment vs. rotation angle relationship of the flexural members. The level of deformations corresponding to each limit state according to FEMA-356 [46] is shown in Figure 3b. The figure shows that between the first yield (point B) and collapse (point C) there were three different levels, which were immediate occupancy (IO), life safety (LS), and collapse prevention (CP). SAP2000 [47] was used for conducting NLTHAs of the model structures including the automatic definition of the limit states related to each cross section of the structure based on ‘Tables 5 and 6’ in FEMA-356 [46]. A modal damping ratio of 2% of the critical damping was used for nonlinear static and dynamic analyses.

3.2. Seismic Hazards and Earthquake Ground Motions

A site-specific hazard curve with various peak ground accelerations was used. A location in Los Angeles with SDS = 1.46 and SD1 = 0.737 with soil type D was used in the study. The response spectra and the hazard curve of the site were obtained from the USGS [48] and shown in Figure 4.
Thirty earthquake records, three sets of ground motions corresponding to 75, 500, and 2500 years return periods as shown in Table 1, were selected from the Pacific Earthquake Engineering Research NGA database [49] and used in previous studies [50,51,52,53,54]. The earthquake records were scaled such that the peak ground accelerations of response spectra geometric mean and the target spectra are the same (Figure 5). Scale factors less than four were used to maintain the original record characteristics [55]. Some earthquake events were repeated for different return periods in Table 1 because the same earthquake event in the NGA data base can be used again by multiplying it by a different scale factor to obtain similar peak ground accelerations of both the earthquake response spectrum and the target spectrum for a specific return period.

4. Seismic LCC Curves of the Case Study Structures

Figure 6 shows the relationship between the fundamental period and the average LCC for structures having a different number of stories with the same number of bays, story height, and mass per floor as the case study buildings. This graph was obtained after conducting a small-scale parametric study on these structures. To construct these graphs, the fundamental mode shapes were normalized to have a modal mass equal to 1.0 for all buildings. The modal participation factor and the fundamental modal amplitude on the roof were calculated from the modal analysis. Each building curve on the graph was adjusted based on the aforementioned characteristics. The figure shows that if the natural period of the structure increases, the LCC increases exponentially. In other words, if the frame system experiences more ductility demand, this causes the LCC to increase. This is an interesting observation because, commonly, allowing the system to be more ductile is considered an advantage for the system since it allows for more energy dissipation. However, on the other hand, it causes larger damage to the structure and increases overall. This necessitates the structural designer to make a preliminary LCC estimation in the early design stage using the basic structural information available (such as the fundamental period), as proposed in the current study. This requires designers to make an early estimate of the LCC using basic structural information, such as the fundamental period, as proposed in the current study.
Figure 7 shows the average LCC graph for the three-story and five-story buildings with an indication of the building’s fundamental period and the corresponding average LCC. Figure 8 and Figure 9 show the mean and the upper and lower bounds of the LCC curves corresponding to the variation in the selected variables for the three- and five-story buildings. As shown in these figures, the upper and lower bounds of the LCC curves are, sometimes, inversely proportional to the upper and lower bounds of the parameter value itself. This is true because some parameters have an inverse relationship with LCCs. For example, increasing the damping ratio will result in more energy dissipation in the system, which causes smaller lateral displacement and, consequently, leads to less LCC estimation. Table 2 shows the mean value and coefficient of variation (COV) used for all variables required to estimate the LCC. All references for these values can be found in [56].
Table 3 and Table 4 show the parameters used in the estimation of the LCC. One case is for estimation of the LCC while all parameters are maintained at their mean values. The other case is for estimation of the LCC with some parameters adjusted to the building under investigation (three or five-story building). For each adjusted parameter, the LCC is re-calculated, and the difference is shown in the table. If the difference has a positive sign, as in the case of the limit state, the calculated LCC resulting from this change gives a larger value than the LCCavg, and the opposite is true. The parameters different from the mean are indicated in the table in bold letters. These parameters have been adjusted to the original building case.
The dispersion in demand parameter βD|Sa has been calculated by conducting full NLTH analyses on the original structure using the 30 time history records. The two alternative calculations for βD|Sa, which were explained earlier, were compared with the case of the NLTH analysis of the MDOF system as shown in Table 5 for the three different return period spectra. As can be observed in the table, the alternative method provides acceptable results, especially for the 75 and 500-year return period earthquakes. For the three-story building, the percentage differences between the proposed method and the NLTH case using the MDOF system turned out to be 7.4% and 3.7% for the case of the 75 and 500-year return period earthquakes, respectively. The differences were 15.0% and 54.0%, respectively, for the five-story building. For the case of the 2500-year return period earthquake, the difference was found to be considerably large for the two buildings. It is worth mentioning that the probability of exceedance of the CP limit state is generally small, as shown in Table 3 and Table 4; the values being 0.49% and 0.55 for the five and three-story buildings, respectively. These small values will overshadow the inaccuracy in calculating the dispersion demand associated with the 2500-year return period.
Figure 10 shows the mean and the mean plus one standard deviation response spectra used for calculating βD|Sa. In this figure, the mean plus one standard deviation response spectra were obtained by multiplying the mean response spectra by amplification factors obtained from previous studies [43,44]. These factors were 2.7, 2.3, and 2.1, respectively, for acceleration, velocity, and displacement sensitive regions of the elastic response spectra with a 5% damping ratio assuming lognormal probability distribution for the spectral ordinates. Figure 11 shows the difference in the probability of exceeding the limit state P(LS|Sa) for both the proposed method and the NLTH reference method. It can be observed that the proposed method provides generally reasonable results for both buildings for the limit states one and two (i.e., IO and LS, respectively). The difference in the five-story building was found to be 0.96%, 0.25%, and 0.38%, respectively, for the IO, LS, and CP limit states. For the three-story building, the differences were 3.84%, 1.28%, and 1.19%, respectively, for the IO, LS, and CP limit states.
It can be observed in Table 3 and Table 4 that the effect of some parameters on the LCC was more significant than others. For example, in the five-story building, the LS and k had more effect on the LCC compared to βD|Sa and t, while in the three-story building the LS and βD|Sa had more effect on the LCC compared to the k and t. In the five-story building, the percentage difference in LCC between the proposed method and the reference method was 17.2%, while in the three-story building this difference reduced to 5.6% (refer to Figure 12). This indicates that the proposed method provides an upper bound for the LCC for both cases with less margin in the case of the three-story building and that the method is applicable for low-rise buildings. This may be attributed to transforming the MDOF to the ESDOF system, wherein accuracy depends on the dominance of the fundamental mode. In addition, it was observed that the change in a lower limit state, such as IO, has more effect on the LCC compared to the change in a higher limit state, such as CP. On the other hand, the structure service life (t) was found to have a limited effect on the LCC; 7.2% and 6.8% of the LCC for the three and five-story buildings, respectively.
It is worth mentioning that the main purpose of the proposed framework is to easily predict the LCC of any building with acceptable accuracy while considering some limitations. For example, it is suitable for low-to-medium-rise regular buildings in which dynamic responses are dominated by the first vibration mode. Additionally, it is suitable for buildings where the equal displacement rule can be applied. Moreover, the framework assumes common values for some variables required for calculating damage probability and cost items. If more specific values of these variables are known, then they can be easily implemented in the framework.
The outcome of this study has the potential to greatly benefit both researchers and end users in the field of structural design and LCC analysis. For researchers, the proposed simplification of LCC estimations using only one factor, the fundamental period, opens up new opportunities for further investigation and advancements in the field. Additionally, the graphical representation of the relationship between LCC and the fundamental period of the structure provides a clear and intuitive tool for researchers to use in future studies. For end users such as structural designers, the outcome of this study can be applied in the early design stages of a structure, where information about the structure is limited, to make informed decisions about the design and ultimately reduce the overall LCC of the building. This can result in cost savings and improved energy efficiency for the building. Additionally, the graphical representation of the relationship between LCC and the fundamental period of the structure can be useful for designers to make trade-off decisions between different design options.

5. Limitations of this Study and Recommendations for Future Research

As this study is focused on the simplification of the seismic LCC estimation process of building structures, it has some limitations. One limitation of this study is that it is based on a set of assumptions, such as the seismic hazard level and the failure probability of structural systems, which may not be accurate for all regions and buildings. Additionally, the study only considers a limited number of structural systems for simplification, and it may not be applicable to other types of systems such as high-rise buildings. Another limitation is that the proposed method is based on SAC-FEMA framework for damage probability calculation, which may not be accurate for all structures. This framework needs to be validated for different types of structures and for different seismic hazard levels. It is also worth mentioning that the proposed method is based on a deterministic approach, which might not be suitable for some situations where a probabilistic approach is needed. These limitations should be considered when interpreting the results and when applying the proposed method to real-world structures.
It also should be mentioned that more research still needs to be carried out for the validation of the proposed method in different types of structures. It would also be valuable to investigate the effect of different seismic hazard levels on the proposed method. The study results can be used to develop design guidelines for cost-effective buildings, and the authors recommend that future research should investigate the possibility of incorporating the proposed method in building design codes and software. This will enable designers to have an automated LCC analysis tool, which will help them to make better decisions in terms of cost and sustainability.

6. Conclusions

In this study, a simplified procedure for estimating the seismic LCC of building structures was presented. The procedure used the fundamental period of the structure as the main input to obtain an average seismic LCC estimation. This average LCC was adjusted using the actual value of the variables related to damage probability and LCC to obtain the final LCC of the building. In the present study, the dispersion measure in demand βD|Sa was calculated using the mean plus one standard deviation of the response spectrum, which eliminates the need for performing nonlinear dynamic time history analyses. The SAC-FEMA framework was utilized along with the transformation of the MDOF system to the ESDOF system. The procedure was applied to three- and five-story buildings located in LA. The results were verified through a comparison with the results obtained from NLTH analysis-based procedure. The main findings can be summarized as follows:
  • The analysis results showed that the proposed procedure provides an upper bound for the seismic LCC of the low-rise buildings, which makes it suitable for preliminary LCC estimations. The percentage difference in LCC between the proposed method and the reference method is 17.2% and 5.6%, respectively, for the five- and three-story buildings.
  • The structure service life (t) is found to have a limited effect on the LCC; 7.2% and 6.8% of the LCC for the three and five-story buildings, respectively. The contributions of the limit state and βD|Sa were found to have a significant impact on LCC estimation.
  • The percentage difference in seismic demand dispersion measure (βD|Sa) compared with the reference method for the 75 and 500 years return period spectra was found to be 7.4% and 3.7% for the three-story building, respectively. For the five-story building, these values were 15.0% and 54.0%, respectively.
  • The differences between the proposed and the reference method in the probability of exceeding a limit state P(LS|Sa) for the IO, LS, and CP limit states were found to be 0.96%, 0.25%, and 0.38%, respectively, for the five-story building. For the three-story building, the differences turned out to be 3.84%, 1.28%, and 1.19%, respectively, for the IO, LS, and CP limit states.
The proposed method eliminates the use of nonlinear dynamic analysis; instead, a simple nonlinear static method such as N2 or capacity spectrum method can be used. The method does not need the construction of fragility curves as the capacity spectral acceleration is calculated from the relationship between the spectral acceleration and the MIDR. The proposed method utilizes the response spectrum approach to calculate the dispersion measure in demand without the need for constructing the IDA curves. The results of the method provide a good approximation for the expected LCC of low-rise buildings and the framework can be generalized for other buildings, which can be used for preliminary LCC estimation needed for decision-makers and early design stage. Moreover, the framework can be easily customized for any particular structure dominated by the first mode shape. Additionally, the framework provides an insight into the effect of important variables on LCCs through a set of graphs showing the relationship between the LCC and the period of the structure for a different variation in the the main input variables of the LCC.

Author Contributions

Conceptualization and methodology, M.N. and J.K.; analysis and validation, M.N.; writing, M.N. and J.K.; supervision and funding acquisition, J.K. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Research Foundation of Korea (NRF), grant funded by the Korea government (MSIT) (No. 2021R1A2C2006631).

Data Availability Statement

The data and resources presented in this research are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

COVcoefficient of variation
CPcollapse prevention limit state
Ddrift demand
ESDOFequivalent single degree of freedom
Hthe height of the structure.
IOimmediate occupancy limit state
Lthe expected life span (service life) of the building
LCCife cycle cost
LCCavgaverage estimation of the LCC
LCCdifthe difference in LCC
LCCfinalthe final estimate of the LCC
LSlife safety limit state
MDOFmulti-degree-of-freedom
MIDRthe mean inter-story drift ratio
Nthe number of the considered limit states
NLTHnonlinear time history
Tbuilding fundamental period
UHSuniform hazard spectrum
Γ the modal participation factor
S a c ^ median drift capacity spectral acceleration
H S a c ^ the probability that the intensity in a site will exceed or equal to S a c ^ on an annual basis
kthe slope of the hazard curve on a log–log plot
bthe slope of the demand–intensity relationship on a log–log plot
βD|sthe standard deviation of the natural logarithm of the drift demands for a given Sa
βcthe standard deviation of natural logarithm for the median capacity of drift C
P L s damage state probability
P(LS|Sa) the probability of exceeding the limit state,
Cothe initial cost of the construction,
λthe rate of return (or annual discount rate)
E[CSD] the damage cost of the building on an annual basis
Pithe ith damage state probability of the building during its service life
Cithe ith damage state cost including damage and its repair.
ΔDthe maximum inter-story drift that a building experiences under an earthquake excitation (demand),
ΔC,Ithe maximum inter-story drift that corresponds to the capacity of the structure at a the ith damage state (capacity).
φ r the modal amplitude at the roof of the fundamental mode shape;
D ¯ m a x M D O F the mean maximum seismic drift demand of the MDOF system
D ¯ m a x E S D O F the mean maximum seismic drift demand of the ESDOF system
Sathe elastic spectral acceleration (a measure of ground motion intensity
D ^ median drift demand
a and bthe regression coefficients for linear regression of drift demand D at the intensity S a in logarithmic space
H s a the corresponding annual probability of exceeding s a
k o and k the linear regression coefficients on a log–log graph
SDSthe design spectral acceleration at short periods
SD1the design spectral acceleration at one second

References

  1. Welch, D.P.; Sullivan, T.J.; Calvi, G.M. Developing direct displacement-based procedures for simplified loss assessment in performance-based earthquake engineering. J. Earthq. Eng. 2014, 18, 290–322. [Google Scholar] [CrossRef]
  2. Mahmoud, H.; Cheng, G. Framework for lifecycle cost assessment of steel buildings under seismic and wind hazards. J. Struct. Eng. 2017, 143, 04016186. [Google Scholar] [CrossRef]
  3. Cutfield, M.; Ryan, K.; Ma, Q. Comparative life cycle analysis of conventional and base-isolated buildings. Earthq. Spectra 2016, 32, 323–343. [Google Scholar] [CrossRef]
  4. Matta, E. Life-cycle cost assessment of inelastic buildings with tuned mass dampers in seismic areas. Int. J. Lifecycle Perform. Eng. 2016, 2, 1–21. [Google Scholar] [CrossRef]
  5. NourEldin, M.; Naeem, A.; Kim, J. Life-cycle cost evaluation of steel structures retrofitted with steel slit damper and shape memory alloy–based hybrid damper. Adv. Struct. Eng. 2019, 22, 3–16. [Google Scholar] [CrossRef]
  6. Torti, M.; Venanzi, I.; Laflamme, S.; Ubertini, F. Life-cycle management cost analysis of transportation bridges equipped with seismic structural health monitoring systems. Struct. Health Monit. 2022, 21, 100–117. [Google Scholar] [CrossRef]
  7. Gidaris, I.; Taflanidis, A.A. Performance assessment and optimization of fluid viscous dampers through life-cycle cost criteria and comparison to alternative design approaches. Bull. Earthq. Eng. 2015, 13, 1003–1028. [Google Scholar] [CrossRef]
  8. Esteva, L.; Díaz-López, O.J.; Vásquez, A.; León, J.A. Structural damage accumulation and control for life cycle optimum seismic performance of buildings. Struct. Infrastruct. Eng. 2016, 12, 848–860. [Google Scholar] [CrossRef]
  9. Vitiello, U.; Asprone, D.; Di Ludovico, M.; Prota, A. Life-cycle cost optimization of the seismic retrofit of existing RC structures. Bull. Earthq. Eng. 2017, 15, 2245–2271. [Google Scholar] [CrossRef]
  10. Noureldin, M.; Kim, J. Parameterized seismic life-cycle cost evaluation method for building structures. Struct. Infrastruct. Eng. 2021, 17, 425–439. [Google Scholar] [CrossRef]
  11. Fragiadakis, M.; Vamvatsikos, D. Fast performance uncertainty estimation via pushover and approximate IDA. Earthq. Eng. Struct. Dyn. 2010, 39, 683–703. [Google Scholar] [CrossRef]
  12. Dolšek, M. Simplified method for seismic risk assessment of buildings with consideration of aleatory and epistemic uncertainty. Struct. Infrastruct. Eng. 2012, 8, 939–953. [Google Scholar] [CrossRef]
  13. Fajfar, P.; Dolšek, M. A practice-oriented estimation of the failure probability of building structures. Earthq. Eng. Struct. Dyn. 2012, 41, 531–547. [Google Scholar] [CrossRef]
  14. Sullivan, T.J.; Welch, D.P.; Calvi, G.M. Simplified seismic performance assessment and implications for seismic design. Earthq. Eng. Eng. Vib. 2014, 13, 95–122. [Google Scholar] [CrossRef]
  15. Fema, P. 58 Seismic Performance Assessment of Buildings; Federal Emergency Management Agency: Washington, DC, USA, 2012.
  16. Kosič, M.; Fajfar, P.; Dolšek, M. Approximate seismic risk assessment of building structures with explicit consideration of uncertainties. Earthq. Eng. Struct. Dyn. 2014, 43, 1483–1502. [Google Scholar] [CrossRef]
  17. Kosič, M.; Dolšek, M.; Fajfar, P. Dispersions for the pushover-based risk assessment of reinforced concrete frames and cantilever walls. Earthq. Eng. Struct. Dyn. 2016, 45, 2163–2183. [Google Scholar] [CrossRef]
  18. El-Din, M.N.; Kim, J. Simplified seismic life cycle cost estimation of a steel jacket offshore platform structure. Struct. Infrastruct. Eng. 2017, 13, 1027–1044. [Google Scholar] [CrossRef]
  19. Liu, M.M.; Mi, B. Life cycle cost analysis of energy-efficient buildings subjected to earthquakes. Energy Build. 2017, 154, 581–589. [Google Scholar] [CrossRef]
  20. Mahdavi, G.; Nasrollahzadeh, K.; Hariri-Ardebili, M.A. Optimal FRP jacket placement in RC frame structures towards a resilient seismic design. Sustainability 2019, 11, 6985. [Google Scholar] [CrossRef]
  21. Fikri, R.; Dizhur, D.; Walsh, K.; Ingham, J. Seismic performance of reinforced concrete frame with masonry infill buildings in the 2010/2011 Canterbury, New Zealand earthquakes. Bull. Earthq. Eng. 2019, 17, 737–757. [Google Scholar] [CrossRef]
  22. Fikri, R.; Dizhur, D.; Ingham, J. Typological study and statistical assessment of parameters influencing earthquake vulnerability of commercial RCFMI buildings in New Zealand. Bull. Earthq. Eng. 2019, 17, 2011–2036. [Google Scholar] [CrossRef]
  23. Hasan, M.; Saidi, T.; Afifuddin, M.; Setiawan, B. The assessment and strengthening proposal of building structure after the Pidie Jaya earthquake in December 2016. J. King Saud Univ.-Eng. Sci. 2021, in press. [Google Scholar] [CrossRef]
  24. Hasan, M.; Afifuddin, M.; Putra, R. Large deformation of building frame structure due to December 7, 2016 Pidie Jaya earthquake. InIOP Conf. Ser. Mater. Sci. Eng. 2020, 933, 012034. [Google Scholar] [CrossRef]
  25. Riswandi, M.; Hasan, M.; Afifuddin, M.; Putra, R.; Maulana, M.; Bilqaini, G.R. Performance assessment of building structures due to Pidie Jaya earthquake using pushover analysis. InIOP Conf. Ser. Mater. Sci. Eng. 2021, 1087, 012022. [Google Scholar] [CrossRef]
  26. Mwafy, A.M.; Elnashai, A.S. Static pushover versus dynamic collapse analysis of RC buildings. Eng. Struct. 2001, 23, 407–424. [Google Scholar] [CrossRef]
  27. Mylonakis, G.; Gazetas, G. Seismic soil-structure interaction: Beneficial or detrimental? J. Earthq. Eng. 2000, 4, 277–301. [Google Scholar] [CrossRef]
  28. ATC, Applied Technology Council. Seismic Evaluation and Retrofit of Concrete Buildings; ATC 40 Report; Applied Technology Council: Redwood City, CA, USA, 1996. [Google Scholar]
  29. Fajfar, P. Capacity spectrum method based on inelastic demand spectra. Earthq. Eng. Struct. Dyn. 1999, 28, 979–993. [Google Scholar] [CrossRef]
  30. Fajfar, P. A nonlinear analysis method for performance-based seismic design. Earthq. Spectra 2000, 16, 573–592. [Google Scholar] [CrossRef]
  31. Cornell, C.A.; Jalayer, F.; Hamburger, R.O.; Foutch, D.A. Probabilistic basis for 2000 SAC federal emergency management agency steel moment frame guidelines. J. Struct. Eng. 2002, 128, 526–533. [Google Scholar] [CrossRef]
  32. Wen, Y.K.; Kang, Y.J. Minimum building life-cycle cost design criteria. I: Methodology. J. Struct. Eng. 2001, 127, 330–337. [Google Scholar] [CrossRef]
  33. Gencturk, B. Life-cycle cost assessment of RC and ECC frames using structural optimization. Earthq. Eng. Struct. Dyn. 2013, 42, 61–79. [Google Scholar] [CrossRef]
  34. Gencturk, B.; Elnashai, A.S. Life cycle cost considerations in seismic design optimization of structures. In Structural Seismic Design Optimization and Earthquake Engineering: Formulations and Applications; IGI Global: Hershey, PA, USA, 2012; pp. 1–22. [Google Scholar]
  35. Fragiadakis, M.; Lagaros, N.D.; Papadrakakis, M. Performance-based multiobjective optimum design of steel structures considering life-cycle cost. Struct. Multidiscip. Optim. 2006, 32, 1–11. [Google Scholar] [CrossRef]
  36. Nour Eldin, M.; Naeem, A.; Kim, J. Seismic retrofit of a structure using self-centring precast concrete frames with enlarged beam ends. Mag. Concr. Res. 2020, 72, 1155–1170. [Google Scholar] [CrossRef]
  37. Choun, Y.S.; ASElnashai, A.S. A simplified framework for probabilistic earthquake loss estimation. Probab. Eng. Mech. 2010, 25, 355–364. [Google Scholar] [CrossRef]
  38. Eldin, M.N. A Simplified Method for Seismic Life-Cycle Cost Estimation of Structures with Application on Sensitivity Analysis. Ph.D. Thesis, Sungkyunkwan University, Suwon, Republic of Korea, 2014. [Google Scholar] [CrossRef]
  39. RS Means Corp. Means Assemblies Cost Data; R.S. Means Company: Kingston, MA, USA, 1997. [Google Scholar]
  40. Lagaros, N.D.; Fotis, A.D.; Krikos, S.A. Assessment of seismic design procedures based on the total cost. Earthq. Eng. Struct. Dyn. 2006, 35, 1381–1401. [Google Scholar] [CrossRef]
  41. Jeong, S.H.; Elnashai, A.S. Probabilistic fragility analysis parameterized by fundamental response quantities. Eng. Struct. 2007, 29, 1238–1251. [Google Scholar] [CrossRef]
  42. Luco, N.; Cornell, C.A. Structure-specific scalar intensity measures for near-source and ordinary earthquake ground motions. Earthq. Spectra 2007, 23, 357–392. [Google Scholar] [CrossRef]
  43. Newmark, N.M.; Hall, W.J. Earthquake Spectra and Design. Engineering Monographs on Earthquake Criteria; Earthquake Engineering Research Institute: Oakland, CA, USA, 1982. [Google Scholar]
  44. Chopra, A.K. Dynamics of structures; Pearson Education: Noida, India, 2007. [Google Scholar]
  45. Abd-alla, M.N. Application of Recent Techniques of Pushover for Evaluating Seismic Performance of Multistory Building; Faculty of Engineering, Cairo University: Oula, Egypt, 2007. [Google Scholar]
  46. FEMA 356 FE. Pre Standard and Commentary for the Seismic Rehabilitation of Buildings; FEMA Publication No. 356. Federal Emergency Management Agency: Washington, DC, USA, 2000.
  47. SAP2000. Analysis Reference Manual; Computers and Structures, Inc.: Berkeley, CA, USA, 2015. [Google Scholar]
  48. USGS. U.S. Geological Survey. 2021. Available online: https://www.usgs.gov/ (accessed on 8 September 2021).
  49. Pacific Earthquake Engineering Research (PEER) Center. NGA and Strong Motion Database. 2021. Available online: https://ngawest2.berkeley.edu (accessed on 15 October 2021).
  50. Noureldin, M.; Ali, A.; Nasab, M.S.; Kim, J. Optimum distribution of seismic energy dissipation devices using neural network and fuzzy inference system. Comput.-Aided Civ. Infrastruct. Eng. 2021, 36, 1306–1321. [Google Scholar] [CrossRef]
  51. Eldin, M.N.; Dereje, A.J.; Kim, J. Seismic retrofit of framed buildings using self-centering PC frames. J. Struct. Eng. 2020, 146, 04020208. [Google Scholar] [CrossRef]
  52. Noureldin, M.; Memon, S.A.; Gharagoz, M.; Kim, J. Performance-based seismic retrofit of RC structures using concentric braced frames equipped with friction dampers and disc springs. Eng. Struct. 2021, 243, 112555. [Google Scholar] [CrossRef]
  53. Noureldin, M.; Ahmed, S.; Kim, J. Self-centering steel slotted friction device for seismic retrofit of beam-column joints. Steel Compos. Struct. 2021, 41, 13–30. [Google Scholar]
  54. Noureldin, M.; Ali, A.; Sim, S.; Kim, J. A machine learning procedure for seismic qualitative assessment and design of structures considering safety and serviceability. J. Build. Eng. 2022, 50, 104190. [Google Scholar] [CrossRef]
  55. Hancock, J.; Bommer, J.J.; Stafford, P.J. Numbers of scaled and matched accelerograms required for inelastic dynamic analyses. Earthq. Eng. Struct. Dyn. 2008, 37, 1585–1607. [Google Scholar] [CrossRef]
  56. Eldin, M.N.; Kim, J. Sensitivity analysis on seismic life-cycle cost of a fixed-steel offshore platform structure. Ocean Eng. 2016, 121, 323–340. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the proposed procedure.
Figure 1. Flowchart of the proposed procedure.
Sustainability 15 02706 g001
Figure 2. Prototype structures: (a) Structural plan; (b) Selected frames.
Figure 2. Prototype structures: (a) Structural plan; (b) Selected frames.
Sustainability 15 02706 g002
Figure 3. Nonlinear modelling of flexural members.
Figure 3. Nonlinear modelling of flexural members.
Sustainability 15 02706 g003
Figure 4. Seismic hazard curve and uniform hazard spectra used in the analysis: (a) Seismic hazard curve of Los Angeles; (b) Uniform hazard spectrum (UHS) with 75 years, 500 years, and 2500 years return periods.
Figure 4. Seismic hazard curve and uniform hazard spectra used in the analysis: (a) Seismic hazard curve of Los Angeles; (b) Uniform hazard spectrum (UHS) with 75 years, 500 years, and 2500 years return periods.
Sustainability 15 02706 g004aSustainability 15 02706 g004b
Figure 5. Response spectra of ground motions and the geometric mean scaled to PGA of UHS with 75, 500, and 2500 years return periods.
Figure 5. Response spectra of ground motions and the geometric mean scaled to PGA of UHS with 75, 500, and 2500 years return periods.
Sustainability 15 02706 g005
Figure 6. A relationship between the fundamental period and the average LCC for structures having a different number of stories.
Figure 6. A relationship between the fundamental period and the average LCC for structures having a different number of stories.
Sustainability 15 02706 g006
Figure 7. The average LCC graphs for the three-story and five-story buildings.
Figure 7. The average LCC graphs for the three-story and five-story buildings.
Sustainability 15 02706 g007
Figure 8. Variation in the LCC–T curves for the variation in the selected parameters (three-story building).
Figure 8. Variation in the LCC–T curves for the variation in the selected parameters (three-story building).
Sustainability 15 02706 g008aSustainability 15 02706 g008b
Figure 9. Variation in the LCC–T curves for variation in the selected parameters (five-story building). (a) Dispersion in structure capacity, (b) Dispersion in structure demand, (c) Cost of limit states, (d) Annual discount rate.
Figure 9. Variation in the LCC–T curves for variation in the selected parameters (five-story building). (a) Dispersion in structure capacity, (b) Dispersion in structure demand, (c) Cost of limit states, (d) Annual discount rate.
Sustainability 15 02706 g009aSustainability 15 02706 g009b
Figure 10. Mean and mean plus one standard deviation of the hazard response spectra.
Figure 10. Mean and mean plus one standard deviation of the hazard response spectra.
Sustainability 15 02706 g010
Figure 11. Difference in the probability of exceeding the limit state P(LS|Sa).
Figure 11. Difference in the probability of exceeding the limit state P(LS|Sa).
Sustainability 15 02706 g011
Figure 12. LCC obtained from the proposed and the NLTH analysis methods.
Figure 12. LCC obtained from the proposed and the NLTH analysis methods.
Sustainability 15 02706 g012
Table 1. Selected ground motions for different hazard levels.
Table 1. Selected ground motions for different hazard levels.
Return PeriodNo.NGA #Earthquake EventYearStationMag.Vs (m/Sec)
75 YR131Parkfield1952Cholame—Shandon6.19385.43
220Northern California1954Ferndale City Hall6.5219.31
368San Fernando1971LA—Hollywood Stor FF6.61316.46
4138Tabas, Iran1978Boshrooyeh7.35324.57
5161Imperial Valley1979Brawley Airport6.53208.71
6163Imperial Valley1979Calipatria Fire Station6.53205.78
7173Imperial Valley1979El Centro Array #106.53471.53
8172Imperial Valley1979El Centro Array #16.53237.33
9175Imperial Valley1979El Centro Array #126.53196.88
10176Imperial Valley1979El Centro Array #136.53249.92
500 YR116Imperial Valley1940El Centro Array #96.95213.44
12158Imperial Valley1979Aeropuerto Mexicali6.53259.86
13161Imperial Valley1979Brawley Airport6.53208.71
14165Imperial Valley1979Chihuahua6.53242.05
15169Imperial Valley1979Delta6.53242.05
16170Imperial Valley1979EC County Center FF6.53192.05
17171Imperial Valley1979El Centro—Meloland Geot6.53264.57
18164Imperial Valley1979Cerro Prieto”6.53202.85
19174Imperial Valley1979El Centro Array #116.53196.25
20178Imperial Valley1979El Centro Array #36.53162.94
2500 YR216Imperial Valley1940El Centro Array #96.95213.44
2230Parkfield1966Cholame—Shandon6.19289.56
2315Kern County1966Taft Lincoln School7.36256.82
2433Parkfield1966Temblor pre-19696.19527.92
2557San Fernando1971Castaic—Old Ridge Route6.61450.28
2695Managua Nicar1972Managua_ESSO6.24288.77
27125Friuli, Italy1976Tolmezzo6.5505.23
28126Gazli, USSR1976Karakyr6.8259.59
29139Tabas, Iran1978Dayhook7.35471.53
30158Imperial Valley1979Aeropuerto Mexicali6.53259.86
Table 2. Mean and covariance of the parameters used for the LCC estimation [56].
Table 2. Mean and covariance of the parameters used for the LCC estimation [56].
kbBD|s
(75, 500, 2500)yr
βcCotλCiηLimit State
(IO, LS, CP)
Mean 2.51.0(0.4,0.5,0.6)0.3-400.030.3,0.7,1.00.03(0.0,0.02,0.03)
COV0.10.050.30.30.30.30.30.30.50.3
Table 3. Summary of the parameters used in LCC estimation of the five-story structure.
Table 3. Summary of the parameters used in LCC estimation of the five-story structure.
ParameterLCCavg
Using an Average Value for All Parameters (1)
LCCfinal
Proposed Method Adjusted to 5 Story Building (2)
LCCdif
The Difference in LCC (3)–(1)
LCC Adjusted to a Particular Parameter (3)
Limit stateLS1LS2LS3LS1LS2LS3--
MIDR (%)12411.52+3532158
t (years)4030−1221683
S A c ^ (g)0.1070.2150.430---
βD|Sa0.40.50.60.530.370.27+2552060
H (Sa )3.82c 10−26.75 × 10−31.19 × 10−3---
ko0.0001444--
k2.502.26−3811424
b1.01.0--
βc0.30.3
P(LS |Sa ) %8.35%1.95%0.49%-
Pi (%)6.39%1.47%0.49%
Co ($)1000
(Assumed reference value)
Ci ($)3007001000Assuming 0.3, 0.7 and 1.0 of Co, respectively.
λ0.030.03
q0.030.03=ln(1 + λ)
α0.660.66=1 − exp(−ql)/ql
(LCC − Co) ($)805810additional cost
LCC ($)18051910
[=1805 + 353 – 122 + 255 − 381]
-
LCC from NLTH and MDOF (LCCMDOF)1630
Table 4. Summary of the parameters used in LCC estimation of the three-story structure.
Table 4. Summary of the parameters used in LCC estimation of the three-story structure.
ParameterLCCavg
Using an Average Value for all Parameters (1)
LCCfinal
Proposed Method Adjusted to 5 Story Building (2)
LCCdif
The Difference in LCC (3)–(1)
LCC Adjusted to a Particular Parameter (3)
Limit stateLS1LS2LS3LS1LS2LS3Difference in LCCCalculate LCC
MIDR (%)12411.52+3962299
t (years)4030−1371766
S A c ^ (g)0.1280.2560.512---
βD|Sa0.40.50.60.540.540.23+3402243
H (Sa)4.29 × 10−27.58 × 10−31.34 × 10−3---
ko0.000252--
k2.502.38−2351668
b1.01.0--
βc0.30.3
P(LS|Sa) %9.36%2.19%0.55%-
Pi (%)7.17%1.65%0.55%
Co ($)1000
(Assumed reference value)
Ci ($)3007001000Assuming 0.3, 0.7 and 1.0 of Co, respectively.
λ0.030.03
q0.030.03= ln(1 + λ)
α0.660.66= 1 − exp(− ql)/ql
(LCC − Co) ($)9031267additional cost
LCC ($)19032267
[=1903 + 396 – 137 + 340 − 235]
-
LCC from NLTH and MDOF2147
Table 5. Comparison of βD|Sa obtained from NLTHA, RSA, and the proposed method.
Table 5. Comparison of βD|Sa obtained from NLTHA, RSA, and the proposed method.
ModelLimit State
(Return Period)
βD|Sa Based on NLTHA Using the MDOFβD|Sa Based on RSA Using the MDOFβD|Sa Based on the Proposed Method Using the ESDOF
IO (75 yr)0.580.540.54
3 storyLS (500 yr)0.560.530.54
CP (2500 yr)0.560.210.23
IO (75 yr)0.610.550.53
5 storyLS (500 yr)0.570.370.37
CP (2500 yr)0.570.270.27
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Noureldin, M.; Kim, J. Simplified Life Cycle Cost Estimation of Low-Rise Steel Buildings Using Fundamental Period. Sustainability 2023, 15, 2706. https://doi.org/10.3390/su15032706

AMA Style

Noureldin M, Kim J. Simplified Life Cycle Cost Estimation of Low-Rise Steel Buildings Using Fundamental Period. Sustainability. 2023; 15(3):2706. https://doi.org/10.3390/su15032706

Chicago/Turabian Style

Noureldin, Mohamed, and Jinkoo Kim. 2023. "Simplified Life Cycle Cost Estimation of Low-Rise Steel Buildings Using Fundamental Period" Sustainability 15, no. 3: 2706. https://doi.org/10.3390/su15032706

APA Style

Noureldin, M., & Kim, J. (2023). Simplified Life Cycle Cost Estimation of Low-Rise Steel Buildings Using Fundamental Period. Sustainability, 15(3), 2706. https://doi.org/10.3390/su15032706

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

Article Metrics

Back to TopTop