Seismic design of RC structures: A critical assessment in the framework of multi‐objective optimization

The assessment of seismic design codes has been the subject of intensive research work in an effort to reveal weak points that originated from the limitations in predicting with acceptable precision the response of the structures under moderate or severe earthquakes. The objective of this work is to evaluate the European seismic design code, i.e. the Eurocode 8 (EC8), when used for the design of 3D reinforced concrete buildings, versus a performance‐based design (PBD) procedure, in the framework of a multi‐objective optimization concept. The initial construction cost and the maximum interstorey drift for the 10/50 hazard level are the two objectives considered for the formulation of the multi‐objective optimization problem. The solution of such optimization problems is represented by the Pareto front curve which is the geometric locus of all Pareto optimum solutions. Limit‐state fragility curves for selected designs, taken from the Pareto front curves of the EC8 and PBD formulations, are developed for assessing the two seismic design procedures. Through this comparison it was found that a linear analysis in conjunction with the behaviour factor q of EC8 cannot capture the nonlinear behaviour of an RC structure. Consequently the corrected EC8 Pareto front curve, using the nonlinear static procedure, differs significantly with regard to the corresponding Pareto front obtained according to EC8. Furthermore, similar designs, with respect to the initial construction cost, obtained through the EC8 and PBD formulations were found to exhibit different maximum interstorey drift and limit‐state fragility curves. Copyright © 2007 John Wiley & Sons, Ltd.     

Seismic performance and new design procedure for chevron‐braced frames

The paper is concerned with the seismic design of steel‐braced frames in which the braces are configured in a chevron pattern. According to EuroCode 8 (EC8), the behaviour factor q, which allows for the trade‐off between the strength and ductility, is set at 2.5 for chevron‐braced frames, while 6.5 is assigned for most ductile steel moment‐resisting frames. Strength deterioration in post‐buckling regime varies with the brace's slenderness, but EC8 adopts a unique q value irrespective of the brace slenderness. The study focuses on reevaluation of the q value adequate for the seismic design of chevron‐braced frames. The present EC8 method for the calculation of brace strength supplies significantly different elastic stiffnesses and actual strengths for different values of brace slenderness. A new method to estimate the strength of a chevron brace pair is proposed, in which the yield strength (for the brace in tension) and the post‐buckling strength (for the brace in compression) are considered. The new method ensures an identical elastic stiffness and a similar strength regardless of the brace slenderness. The advantage of the proposed method over the conventional EC8 method is demonstrated for the capacity of the proposed method to control the maximum inter‐storey drift. The q values adequate for the chevron‐braced frames are examined in reference to the maximum inter‐storey drifts sustained by most ductile moment‐resisting frames. When the proposed method is employed for strength calculation, the q value of 3.5 is found to be reasonable. It is notable that the proposed method does not require larger cross‐sections for the braces compared to the cross‐sections required for the present EC8 method. Copyright © 2005 John Wiley & Sons, Ltd.     

Displacement capacity of masonry buildings as a basis for the assessment of behavior factor: an experimental study

The results of shaking table tests of a series of 1:5 scale masonry building models have been used for the assessment of values of structural behavior factor q for masonry structures, seismic force reduction factors proposed for the calculation of design seismic loads by Eurocode 8, European standard for the design of structures for earthquake resistance. Six models have been tested, representing prototype buildings of two different structural configurations and built with two different types of masonry materials. The study indicated that the reduction of seismic forces for the design depends not only on the type of masonry construction system, but also on structural configuration and mechanical characteristics of masonry materials. It has been also shown that besides displacement and energy dissipation capacity, damage limitation requirement should be taken into account when evaluating the values of behavior factor. On the basis of analysis of experimental results a conclusion can be made, that the values at the upper limit of the proposed range of values of structural behavior factor q for unreinforced and confined masonry construction systems are adequate, if pushover methods are used and the calculated global ductility of the structure is compared with the displacement demand. In the case where elastic analysis methods are used and significant overstrength is expected, the proposed values are conservative. However, additional research and parametric studies are needed to propose the modifications.     

Analytical prediction of the collapse earthquake for R/C buildings: Case studies

An analytical methodology for estimating the values of behaviour factors (q-factors), suggested in a companion paper, is applied to two typical reinforced concrete structures, a frame and a dual system, designed according to current European trends in seismic design. The two case studies clearly indicated the need to consider both local and overall collapse criteria, since either one can be critical for the determination of the behaviour factor, depending on the type of the structure. The calculated values of q-factors were close to, but higher than those specified by modern seismic codes, the larger discrepancy occurring in the case of the dual system.     

Life-cycle cost analysis of design practices for RC framed structures

The objective of this study is to perform life-cycle cost analysis on three design practices namely weak ground storey, short and floating columns and their combinations. Life-cycle cost analysis is recognized as the only suitable tool for assessing the structural performance when the structure is expected to be functional for a long period of time. Life-cycle cost analysis is considered in this study assessing the behaviour of the three design practices against earthquake hazard. Although, a number of checks are performed in order to reduce the influence of these design practices on the seismic behaviour of reinforced concrete (RC) framed structures, it was found that the total life-cycle cost of partially infilled RC designs is significantly increased compared to that of the fully infilled one. Through the test example examined in the framework of this study general conclusions are obtained regarding the behaviour of the three design practices.     

The effects of repeated earthquake ground motions on the non‐linear response of SDOF systems

In many parts of the world, the repetition of medium–strong intensity earthquake ground motions at brief intervals of time has been observed. The new design philosophies for buildings in seismic areas are based on multi‐level design approaches, which take into account more than a single damageability limit state. According to these approaches, a sequence of seismic actions may produce important consequences on the structural safety. In this paper, the effects of repeated earthquake ground motions on the response of single‐degree‐of‐freedom systems (SDOF) with non‐linear behaviour are analysed. A comparison is performed with the effect of a single seismic event on the originally non‐damaged system for different hysteretic models in terms of pseudo‐acceleration response spectra, behaviour factor q and damage parameters. The elastic–perfect plastic system is the most vulnerable one under repeated earthquake ground motions and is characterized by a strong reduction of the q‐factor. A moment resisting steel frame is analysed as well, showing a reduction of the q‐factor under repeated earthquake ground motions even larger than that of an equivalent SDOF system. Copyright © 2002 John Wiley & Sons, Ltd.     

A unified approach for the design of high ductility steel frames with concentric braces in the framework of Eurocode 8

Eurocode 8 (EC8) stipulates design methods for frames with diagonal braces and for chevron braced frames, which differ as regards the numerical model adopted, the value of the behavior factor q and the estimation of the lateral strength provided by braces. Instead, in this paper, the use of the same design method is suggested for both types of concentrically braced frames. The design method is a generalization of the one proposed for chevron braced frames in a previous study. A numerical investigation is conducted to assess the reliability of this design method. A set of concentrically braced frames is designed according to the EC8 and proposed design methods. The seismic response of these frames is determined by nonlinear dynamic analysis. Finally, it is demonstrated that the proposed design method is equivalent to those provided by EC8, because it can ensure the same level of structural safety which would be expected when using EC8. Copyright © 2013 John Wiley & Sons, Ltd.     

Design method and behavior factor for steel frames with buckling restrained braces

Buckling restrained braces (BRBs) are very effective in dissipating energy through stable tension–compression hysteretic cycles and have been successfully experimented in the seismic protection of buildings. Their behavior has been studied extensively in the last decades and today the level of performance guaranteed by these devices and the technological constrains that have to be fulfilled to optimize their behavior are well known. Furthermore, several companies in the world have developed their own BRBs and are now producing them. In spite of this, many seismic codes (for instance, the EuroCode 8) do not stipulate provisions for the design and construction of earthquake‐resistant structures equipped with BRBs. This discourages the structural engineering community from using these devices and seriously limits their use in structural applications. In this paper a procedure for the seismic design of steel frames equipped with BRBs is proposed. Furthermore, the paper presents a numerical investigation aimed at validating this design procedure and proposing the value of the behavior factor q that should be used for this structural type. To this end, a set of frames with BRBs is first designed by means of several values of q. Then, the obtained frames are subjected to a set of accelerograms compatible with the elastic response spectrum considered in design. The seismic response of the frames is determined by nonlinear dynamic analysis and represented in terms of the ductility demand of BRBs and the internal force demand of nondissipative members (beams and columns). Finally, the largest value of q that leads to acceptable seismic performance of the analyzed frames is assumed as adequate. The value of q is given in the paper as a continuous function of the assumed ductility capacity of the BRBs. Copyright © 2012 John Wiley & Sons, Ltd.     

A comparative application of different EC8-3 procedures for the seismic safety assessment of existing structures

An application of the EC8-3 procedure for safety assessment is presented herein. Besides testing the applicability of the code procedure, this application aims to assess the consistency of the safety assessment results that are obtained. Based on the application of the EC8-3 procedure, the study assesses if the different methods of analysis that are permitted lead to similar safety results and identifies the factors that may affect these results. Furthermore, the results obtained by the EC8-3 procedure are complemented by a probabilistic approach yielding their corresponding fragility values. By comparing the D/C ratios and the probabilistic results, the application aims to determine if similar D/C ratios lead to similar probabilistic results. Furthermore, the application tries to determine if a correlation can be established between the D/C ratios and the expected fragility values.     

Emulation of simulated earthquake catalogues

In earthquake occurrence studies, the so-called q value can be considered both as one of the parameters describing the distribution of interevent times and as an index of non-extensivity. Using simulated datasets, we compare four kinds of estimators, based on principle of maximum entropy (POME), method of moments (MOM), maximum likelihood (MLE), and probability weighted moments (PWM) of the parameters (q and τ ₀) of the distribution of inter-events times, assumed to be a generalized Pareto distribution (GPD), as defined by Tsallis (1988) in the frame of non-extensive statistical physics. We then propose to use the unbiased version of PWM estimators to compute the q value for the distribution of inter-event times in a realistic earthquake catalogue simulated according to the epidemic type aftershock sequence (ETAS) model. Finally, we use these findings to build a statistical emulator of the q values of ETAS model. We employ treed Gaussian processes to obtain partitions of the parameter space so that the resulting model respects sharp changes in physical behaviour. The emulator is used to understand the joint effects of input parameters on the q value, exploring the relationship between ETAS model formulation and distribution of inter-event times.     

A methodology for the probabilistic assessment of behaviour factors

Given the importance that traditional force-based seismic design still currently exhibits, studies addressing issues related to the definition of the behaviour factor values are considered to be of most interest. A probabilistic methodology is proposed for the calibration of the q-factor relating its value with two fundamental parameters, the displacement ductility capacity measured at a relevant location of the structure and the failure probability P _f . The general foundation of this procedure is based on the probabilistic quantification of the seismic action and, by applying a transformation procedure, of the structural seismic demand in terms of displacement ductility. By recalling well established structural reliability procedures and by making use of nonlinear analysis methods, both static and dynamic, a general probabilistic framework, which is able to relate the ductility capacity, the failure probability P _f and the behaviour factor, is defined. In order to illustrate some of the potentialities of the methodology, an application example is presented, addressing the q-factor assessment for a set of regular and irregular reinforced concrete frame structures, enforcing a given P _f and two different ductility levels.     

Analytical prediction of the collapse earthquake for R/C buildings: Suggested methodology

Criteria to be used for the analytical prediction of collapse of R/C buildings subjected to earthquake excitation are suggested. It is attempted to include in the criteria as many of the parameters involved as possible, while keeping the complexity, as well as the cost of analysis reasonably low. Two combined collapse criteria are finally suggested, corresponding to member failure and storey failure. Using these criteria, quantitative estimations can be obtained of the actual values of behaviour factors (q-factors) used in modern seismic codes. The ground motion dependence of these factors, as well as the sensitivity of their calculated values to the assumption made about the stiffness of R/C members, is illustrated with reference to two typical structures.     

Estimates of Horizontal Groundwater Flow Velocities in Boreholes

Currently, monitoring tools can be deployed in observation boreholes to better assess groundwater flow, flux of dissolved contaminants and their mass discharge in an aquifer. The relationship between horizontal water velocity in observation boreholes and Darcy fluxes in the surrounding aquifer has been studied for natural flow conditions (i.e., no pumping). Interpretation of measurements taken with dilution tests, the colloidal borescope, the Heat Pulse Flowmeter, and other techniques require the conversion of observed borehole velocity u to aquifer Darcy flux q_∞ . This conversion is typically done through a proportionality factor α = u/q_∞ . In experimental studies as well as in theoretical developments, reported values of α vary almost three orders of magnitude (from 0.5 to 10). This large variability in reported values of α could be explained by: (1) unclear distinction between Darcy flux and water seepage velocity, (2) unclear definition of water velocity in the borehole, (3) effects of well screen and the presence of the measurement device itself on the observable velocities, and (4) hydraulic conditions in the borehole annulus. We address (1), (2) from a conceptual/theoretical perspective, and (3) by means of numerical simulations. We show that issue (1) in low porosity aquifers can yield to order-of-magnitude discrepancies in estimates of q_∞ ; (2) may result in discrepancies of up to 50%, and (3) can cause differences up to 20% of water velocity in the borehole void space compared to the theoretical case of an open borehole.     

Optimum seismic design of reinforced concrete frames according to Eurocode 8 and fib Model Code 2010

Traditional seismic design, like the one adopted in Eurocode 8 (EC8), is force‐based and examining a single level of seismic action. In order to provide improved control of structural damage for different levels of seismic action, the new fib Model Code 2010 (MC2010) includes a fully fledged displacement‐based and performance‐based seismic design methodology. However, the level of complexity and computational effort of the MC2010 methodology is significantly increased. Hence, the use of automated optimization techniques for obtaining cost‐effective design solutions becomes appealing if not necessary. This study employs genetic algorithms to derive and compare optimum seismic design solutions of reinforced concrete frames according to EC8 and MC2010. This is important because MC2010 is meant to serve as a basis for future seismic design codes. It is found that MC2010 drives to more cost‐effective solutions than EC8 for regions of low seismicity and better or similar costs for regions of moderate seismicity. For high‐seismicity regions, MC2010 may yield similar or increased structural costs. This depends strongly on the provisions adopted for selecting the set of ground motions. In all cases, MC2010 provides enhanced control of structural damage. Copyright © 2016 John Wiley & Sons, Ltd.     

Consistent inelastic design spectra: Hysteretic and input energy

This is the second of two companion papers on inelastic design spectra (for strength, displacement, hysteretic and input energy) for systems with a prescribed ductility factor. All the spectra are consistent (interrelated and based on the same assumptions). This paper deals with two quantities related to cumulative damage: hysteretic and input energy. The input data for the procedure are the characteristics of the expected ground motion in terms of a smooth elastic pseudo-acceleration spectrum and the time integral of the square of the ground acceleration ∫a² dt. Simple, approximate expressions for two dimensionless parameters (the parameter γ and the hysteretic to input energy ratio E_HE_I) have been proposed. The parameter 7, which controls the reduction of the deformation capacity of structures due to low-cycle fatigue, depends on the natural period of the system, the prescribed ductility factor, the hysteretic behaviour and the ground motion characteristics. The ratio E_H/E_I is influenced by damping, the ductility factor and the hysteretic behaviour. Very good approximations to the inelastic spectra for hysteretic and input energy can be derived from the elastic spectrum using the spectra for the reduction factor R, proposed in the companion paper, and the proposed values for γ and E_H/E_I     

Evaluation of ASCE-41, ATC-40 and N2 static pushover methods based on optimally designed buildings

Alternative static pushover methods for the seismic design of new structures are assessed with the aid of advanced computational tools. The current state-of-practice static pushover methods as suggested in the provisions of European and American regulations are implemented in this comparative study. In particular the static pushover methods are: the displacement coefficient method of ASCE-41, the ATC-40 capacity spectrum method and the N2 method of Eurocode 8. Such analysis methods are typically recommended for the performance assessment of existing structures, and therefore most of the existing comparative studies are focused on the performance of one or more structures. Therefore, contrary to previous research studies, we use static pushover methods to perform design and we then compare the capacity of the outcome designs with reference to the results of nonlinear response history analysis. This alternative approach pinpoints the pros and cons of each method since the discrepancies between static and dynamic analysis are propagated to the properties of the final structure. All methods are implemented in an optimum performance-based design framework to obtain the lower-bound designs for two regular and two irregular reinforced concrete building configurations. The outcome designs are compared with respect to the maximum interstorey drift and maximum roof drift demand obtained with the Incremental Dynamic Analysis method. To allow the comparison, also the life-cycle cost of each design is calculated; i.e. a parameter that is used to measure the damage cost due to future earthquakes that will occur during the design life of the structure. The problem of finding the lower bound designs is handled with an Evolutionary type optimization algorithm.     

Confidence interval of real‐time forecast stages provided by the STAFOM‐RCM model: the case study of the Tiber River (Italy)

This study proposes a statistically based procedure to quantify the confidence interval (CI) to be associated to the stages forecast by a simple model called STAge FOrecasting Model‐Rating Curve Model (STAFOM‐RCM). This model can be used for single river reaches characterized by different intermediate drainage areas and mean wave travel times when real‐time stage records, cross section surveys and rating curves are available at both ends. The model requires, at each time of forecast, an estimate of the lateral contribution q_for between the two sections delimiting the reach. The CI of the stage is provided by analyzing the statistical properties of model output in terms of lateral flow, and it is derived from the CI of the lateral contribution q_for which, in turn, is set up by associating to each q_for the q_opt which allows STAFOM‐RCM to reproduce the exact observed stage. From an operative point of view, the q_for values are ranked in order of magnitude and subdivided in classes where the q_opt values can be represented through normal distributions of proper mean and variance from which an interval of selected confidence level for q_for is computed and transferred to the stage. Three river reaches of the Tiber river, in central Italy, are used as case study. A sensitivity analysis is also performed in order to identify the minimum calibration set of flood events. The CIs obtained are consistent with the level of confidence selected and have practical utility. An interesting aspect is that different CI widths can be produced for the same forecast stage since they depend on the estimate of q_for made at the time of forecast. Overall, the proposed procedure for CI estimate is simple and can be conveniently adapted for other forecasting models provided that they have physically based parameters which need to be updated during the forecast. Copyright © 2012 John Wiley & Sons, Ltd.     

Consistent inelastic design spectra: Strength and displacement

A procedure for the determination of inelastic design spectra (for strength, displacement, hysteretic and input energy) for systems with a prescribed ductility factor has been developed. All the spectra are consistent (interrelated and based on the same assumptions). This is the first of two companion papers which deals with the ‘classical’ structural parameters: strength and displacement. The input data are the characteristics of the expected ground motion in terms of a smooth elastic pseudo-acceleration spectrum. Simple, approximate expressions for the strength reduction factor R are proposed. The value of R depends on the natural period of the system, the prescribed ductility factor, the hysteretic behaviour, damping and ground motion. Fairly accurate approximations to the inelastic spectra for strength and displacement can be derived from the elastic spectrum using the proposed values for R.     

Maximum likelihood (ML) estimation of depth from the spectrum of aeromagnetic fields

Summary Starting from the assumption that the aeromagnetic field is a Gaussian random function the probability density function of log-radial spectrum is shown to be a slightly asymmetric non-Gaussian function, [2q/2 (q/2)]^–1 q( ^q/2)–1 exp(q/2(r-exp(r))). The depth to magnetic layer is determined by maximum likelihood (ML) technique and is compared with the least square (LS) estimate. The difference between the two is only marginal, about 15%. The least square estimate is lower than the maximum likelihood estimate.     

考虑模糊可靠度的抗震结构优化设计

在结构优化的目标函数中，不仅应考虑结构造价，而且应考虑结构服役期的损失期望。为此，我们定义了抗震结构模糊可靠度的概念。利用抗震设计原则和地震危险性分析给出了一个计算模糊可靠度的方法。然后提出了三步优化法；第一步，寻求一系列对应于不同设防烈度的最小造价设计；第二步，同时考虑造价和地震引起的损失期望，寻求最优设防烈度；第三步，寻找相应于此最优烈度的最小造价设计。