## 1. Introduction

Past destructive earthquakes (e.g. the 1994 Northridge earthquake and the 1995 Kobe earthquake) have left a clear signature on the engineering community worldwide, changing thinking of structural engineers [1-2]. As such, after holding several workshops and conferences, an innovative approach namely Performance-Based Design (PBD) was presented by modern guidelines [3-5]. In principle, a structure designed using PBD approach should meet performance objectives in accordance with a set of specified reliabilities over its service life. This is aimed to reach structural design candidates associated with more predictable seismic behavior, quantifying and controlling the risk at an engineered acceptable level.

Both seismic demands and capacity parameters, that are inherently uncertain, are highly influential on the acceptable performance level of a structure. Furthermore, due to the fact that a structure on underlying soil is not rigid, soil-structure interaction (SSI) affects the responses of structures during an earthquake. Obviously, ignoring the SSI effects could lead to unrealistic structural responses and seismic demands. Hence, the effects of SSI should be considered in the seismic responses of structures [6]. Therefore, soil type, material properties of the structure, and ground motion characteristics randomly affect the seismic structural responses.

Deterministic structural optimization without considering the uncertainties in design manufacturing and operating processes may lead to unreliable design resulting in inappropriate balance between cost and safety. A proper design procedure must reasonably account for the inherent uncertain nature of a structural system and associated external load [7]. In structural optimization, non-deterministic performance of structures can be taken into account using robust design optimization (RDO) [8] and reliability-based design optimization (RBDO) [9]. RDO aims to minimize variation of the objective function, but RBDO optimizes the structural cost under reliability of the constraints.

A few studies have been implemented in a structural optimization problem, where RBDO is incorporated into PBD concept. Foley *et al*. [10] proposed a state-of-the-art model code and a PBD methodology. The methodology was applied to multiple-objective optimization problems for single storey and multi-storey structural frameworks with fully and partially restrained connections. Lagaros *et al.* [11] introduced a tool for the solution of realistic structural optimization problems that incorporate PBD under seismic loading into RBDO. The tool consisted of two distinctive methodologies based on artificial neural networks (ANNs). Fragiadakis *et al*. [12] presented optimum seismic design of reinforced (RC) structures considering the reliability constraints and PBD, where RBDO was implemented by the evolution strategies (ES). In the study by Moller *et al.* [13], the concept of PBD with RBDO was implemented by consideration of the uncertainties in the structural demands and capacities in order to evaluate reliability associated with each of the required performance levels. Khatibinia *et al*. [14] introduced RBDO of RC structures including SSI effects. In their study, the uncertainty of the structural demand was in terms of the random properties of the structure, underlying soil and the uncertain characteristics of artificially generated earthquakes. Also, the structural capacity associated with each of the required performance level in the concept of PBD was treated as an uncertain quantity.

Nonlinear dynamic analysis of structures using finite element method requires much computational effort. This drawback may accentuate when the nonlinear dynamic structural responses are required in RBDO of the structure using the Monte-Carlo Simulation (MCS) method, the importance sampling technique and the response surface method. In order to obtain an acceptable confidence within probabilities of the order close to 10^{-4} - 10^{-6}, the MCS method requires a large number of structural analyses. Based on Lagaros *et al*. [11], for such results, the large number of analyses ranges from *et al.* [19] stipulated that there is no difference between the performance of the SVM-based response surface method (SVM-RSM) and the radial basis function neural network-based response surface method (RBFN-RSM). As a comparative study, Moura *et al*. [20] assessed the SVM effectiveness in forecasting time-to-failure and reliability of engineered components based on time series data. The efficacy of SVM with respect to other learning methods was shown. In the context of structural reliability assessment, a sampling method based on the adaptive Markov chain simulation and support vector density estimation was also developed by Dai *et al*. [21]. In their study, the application of SVM was proposed as a density estimator for structural reliability analysis.

In this chapter, RBDO of RC structures with considering SSI effects under time-history earthquake loading is presented in accordance with the PBD concept of SEAOC guidelines [3]. In this work, a new discrete gravitational search algorithm (DGSA) and an efficient proposed meta-model were introduced for performing RBDO of RC structures [22]. The objective function is the total cost of the structure while the constraints are treated as deterministic and probabilistic. The annual probability of non-performance for each performance level is considered as the probabilistic constraint in RBDO procedure. The new DGSA based on the fundamental concept of the standard GSA [23] is introduced for finding the optimal designs in the RBDO procedure. In DGSA, the position of each agent is presented in positive integer numbers. Also, the velocity of each agent is modified based on the particle swarm optimizer with passive congregation (PSOPC) which was proposed by He *et al*. [24]. The modifications can improve the global exploration ability of DGSA and overcome the shortcomings of the Binary GSA (BGSA) model introduced by Rashedi *et al*. [25]. A meta-model-based MCS method is also presented herein to take into account the probabilistic constraint in conjunction with the nonlinear finite element analysis (FEM) of SSI system. Due to the fact that the computational cost of MCS for structural reliability analysis is high, the meta-model is proposed to predict the structural seismic responses of SSI system, and significantly reduce the computational effort of the RBDO process. The meta-model is a combination of weighted least squares support vector machine (WLS-SVM) [26] and Morlet wavelet kernel function, which is called WWLS-SVM [14, 22, 27]. The selection of WWLS-SVM parameters efficiently affects the prediction accuracy of WWLS-SVM. Hence, the parameters of WWLS-SVM and wavelet kernel are assigned by using the standard gravitational search algorithm (GSA).

Numerical examples show that the wavelet as a kernel function is much better than those of the common kinds as kernel function in WLS-SVM. The accuracy and generalization of WWLS-SVM is improved using GSA. Furthermore, numerical results demonstrate the efficiency and computational advantages of the proposed DGSA for RBDO of structures.

## 2. RBDO of RC structures

### 2.1. Formulation of optimization

Seismic design optimization of RC structures under time-history earthquake loads is an ongoing research topic and has received great attention among researchers [14, 27-33]. As such, RBDO of RC structures with the consideration of SSI effects was investigated in accordance with PBD concept of SEAOC guidelines [3] under seismic loading. This work incorporates the acceptable performance levels and the RBDO theory to compare the achieved annual probability of non-performance with target values for each performance level. The objective of the RBDO problem is to minimize the total cost whereas the deterministic and probabilistic constraints should not exceed a specified target.

The RBDO problem of RC structures can be formulated in the following form:

where *C*_{TOT} is the total cost; *j*th performance level (performance function), *N*_{g} is the number of deterministic constraints, *X*, is expressed by *R*^{d} and design variables, *x*_{k,} can take values only from this set. The vector

### 2.2. Life- cycle cost assessment of RC structure

The total cost, *C*_{TOT}, of a structure is the initial structural cost for a new structure construction and the repair cost from an earthquake and different levels of damage that may occur during the life of structure. This cost can be expressed as a function of the design vector *X* and the time *t* as follows:

where *C*_{IC} is the initial cost of structure; *C*_{RC} is the present value of the repair cost. In this study, the initial cost is considered as the sum of the total cost of concrete, *C*_{C}, and the total cost of reinforcing steel, *C*_{S}, is given [14]:

where *b*_{i}, *h*_{i,} and *L*_{i} are the section dimensions of *i*th element and its length; and *A*_{Si} is the total reinforcement in the section of element. *N*_{e} is the number of the elements of the structure.

The repair cost refers to the cost of damage level from earthquake that may occur during the life of a structure. In this study, the overall damage index, *DI*_{overall}, is considered as an indicator of structural damage. The total expected cost of repair based on the overall damage index is expressed as follows [13]:

where *DI*_{overall}.

where *r* are the mean occurrence of earthquakes and the discount rate, respectively. *C*_{0} is the complete replacement cost, with *r*.

### 2.3. Constraint handling approach

A comprehensive overview of the most popular constraint handling approaches used in conjunction with meta-heuristic optimization methods was presented in the literature review by Coello Coello [35]. In the present study, the external penalty function method as one of the most common forms of the penalty function in the structural optimization [15, 27-29, 36-39] is employed to transform constrained RBDO problem into unconstrained one as follows:

where *PF* and *r*_{p} are the modified function, the penalty function, and an adjusting coefficient, respectively. The penalty function based on the violation of normalized constraints [36] is defined as the sum of all active constraints violations as indicated:

This formulation allows solutions with violated constraints, and the objective function is always greater than the non-violated one.

## 3. Reliability assessment of RC structure

In order to evaluate the system reliability index corresponding to each of the performance levels, RC structures should be assessed in the RBDO procedure [14, 22]. The system reliability index corresponding to each of the performance levels are estimated by MCS method. In the following subsections, the procedure of assessment of RC structures is explained.

### 3.1. Required database

In PBD approach, many uncertain variables influence the structural seismic responses. In the studies by Khatibinia *et al*. [14, 22], material properties of concrete, steel and soil, as well as earthquakes are considered as intervening uncertain variables. Because of a few historical records of earthquake for a selected site, selection of a proper ground motion record for a site is often difficult, even impossible in some cases. To overcome this problem, artificial earthquakes, statistically influenced by desired properties of a selected site, are utilized in seismic design of structures. The spectral representation method based on time domain procedure can be used for generation of artificial earthquakes [40]. The proper parameters for generation of artificial earthquakes are selected according by values proposed by Möller et al. [13]. In order to perform RBDO of RC structures using the proposed meta-model-based MCS, samples are generated randomly, and are used to train and test the meta-model. Inputs of the meta-model include the random combinations of intervening variables. In order to generate the database, seven random Peak Ground Acceleration (PGA) values are chosen. The PGA values are equal to 260, 350, 400, 550, 650, 700 and 800 (cm/sec^{2}). Accordingly, by using the Latin Hypercube Design (LHD) sampling method [41], and considering the 120 combinations based on the intervening variables for each PGA value, the total number 840 combinations are generated. Then, corresponding to the each PGA of each 840 combinations, five artificial earthquakes, as sub-combinations, with random phase angles are generated corresponding to each PGA value. For each of the 840 combinations, in the first phase, the structure is analyzed subjected to the combination of gravity loads according to ACI code [42]. Steel reinforcement ratios of longitudinal bars of structural elements’ cross-sections, *u*_{max}, maximum inter-storey drift, *DR*_{max}, maximum local damage index, *DIL*_{max}, and overall damage index, *DI*_{overall}, are considered as the structural seismic responses. After that, the mean, *i*th seismic response, *R*_{i}, corresponding to each of the 840 combinations subjected to five artificial earthquakes, are achieved. In this work, the modified Park-Ang damage index [43] as one of the most acceptable indices in seismic damage analysis of structures is used for calculating *DIL*_{max} and *DI*_{overall}.

### 3.2. Limit state functions

The operational, life safety and collapse prevention levels have been defined as the performance levels. A performance level depends on some limit state functions. A limit state function,

where

The limit state functions and their probability distribution function (PDF) for the performance levels, according to SEAOC guidelines (2000), are shown in Table 1.

Performance level | Limit sate function | G (X¯) | PDF |

Operational | |||

Elastic roof displacement | Lognormal | ||

Inter-storey drift | Lognormal | ||

Life safety | |||

Inter-storey drift | Lognormal | ||

Max. local damage index | Beta | ||

Global damage index | Beta | ||

Collapse | |||

Inter-storey drift | Lognormal | ||

Max. local damage index | Beta | ||

Global damage index | Beta |

In Table 1,

**3.3. Annual probability of non-performance**

The non-performance probability,

where

For each performance level, the total exceeding probability,

where

## 4. Seismic responses and SSI system

### 4.1. Seismic responses of SSI system

There are two main approaches for modeling and analyzing SSI systems, namely the direct method and the substructure method either in time domain or in frequency domain [6]. Considering the discretized dynamic equations of structure and soil simultaneously, the direct method models the soil and structure together, and the responses of soil and structure are determined simultaneously by analyzing SSI system in each time step [22].

In the direct method, the discretization of nonlinear dynamic equations can be expressed in FEM framework as:

where *t*. The term *x* direction. The column matrix,

Over the past two decades, the use of damage and energy concepts for the seismic performance evaluation and design of structures has attracted considerable attention among the researchers [30, 44-47]. These concepts can be simultaneously used through a combined damage index namely Park-Ang damage index. The index is taken into account as a combined index, defined as the linear combination of the maximum displacement and the hysteretic energy dissipation for a structural element. For this reason, the damage index [44] is one of the indices that have widely been used for damage assessment and damage-based design of RC structures [14, 22, 28-30, 45-46, 48]. As shown in Table 1, some limit states of the performance levels depend on the damage indices.

An improved version of the index namely modified Park-Ang damage index [43] is defined based on the cross-section deformation of structural elements as:

where

where *i*. *Ne* and *Ns* are the number of structural elements and stories, respectively.

### 4.2. Finite element model of SSI system

OpenSEES [49], as an open-source computational software framework, is used for by simulation of SSI system, and performing nonlinear dynamic analyses of SSI system depicted in Fig. 1. Assuming materials of constant properties over its depth, soil encompasses different layers, and the foundation is considered as rigid strip footing. Beams and columns of structure are modeled using force-based nonlinear beam-column element with considering the spread plasticity along the element’s length. The integration along each element is based on Gauss-Lobatto quadrature rule. Also, the infinite boundaries of soil are modeled using the artificial boundaries (Fig. 1). The model of soil-structure system shown in Fig. 1 was successfully used by [14, 22, 31].

The Kent-Scott-Park model [50] is utilized for modeling the confined and unconfined concrete of cross-sections of structural elements. The constitutive parameters of this model are: *f*_{c}=concrete peak strength in compression, *f*_{u}=residual strength,

where

In this study, the one-dimensional *J*_{2} plasticity model with linear hardening is utilized for modeling the constitutive behavior of the steel reinforcement. The material parameters defining *J*_{2} plasticity model are: *f*_{y}=yield strength, *H*=hardening modulus and *E*=Young’s modulus (Fig. 2(b)).

Soil layers are modeled using isoperimetric four-node quadrilateral finite elements and assuming bilinear displacement interpolation. The plane strain condition is assumed for the soil domain with considering a constant soil thickness corresponding to the inter-frame distance. The material of the soil is modeled using a modified pressure-independent multi-yield-surface *J*_{2} plasticity model [53]. As shown in Fig. 3, this nonlinear model of soil material is described by a shear stress-strain backbone curve. The detailed description of the parameters of the shear stress-strain backbone curve can be found in [53].

One of the major problems in SSI system for infinite media has been the modeling of the domain boundaries. Infinite boundaries have to absorb all outgoing waves and reflect no waves back into the computational domain. In this study, the standard viscous boundary proposed by Lysmer and Kuhlemeyer [54] is used for this purpose. This boundary can be described by two series of dashpots oriented normal and tangential to the boundary of a finite element mesh (Fig. 1) as follows:

where *a* and *b* are dimensionless parameters to be determined, and *i=*1,2*,...,n* representing the soil layers’ number. Other parameters of soil material depend on

The material damping matrix, *C*, of the SSI system is assembled by its corresponding damping matrices of structure and soil and considering the Rayleigh damping model. The factors of proportionality for damping matrices of structure and soil are computed based on 5% and 10% viscous damping respectively for structure and soil. The P-Δ effects are regarded in nonlinear time-history analyses. The accelerated Newton algorithm based on Krylov subspaces [55] is utilized for solving nonlinear equations of SSI system equilibrium.

## 5. Artificial earthquakes

For RBDO of structures, it is then necessary to utilize accelerograms of compatible characteristics with a desired site. It is often difficult or impossible in some cases to choose a proper record for a site, since historically recorded accelerograms for a given site could be limited or scare. Hence, artificial earthquakes that are statistically influenced by desired properties of the given site are very useful for seismic design of structures. In this work, spectral representation method based on time domain procedure is used for the generation of synthetic ground motion records. The non-stationary ground motion is simulated using this method as [13]:

where

where

where

Parameter | |||

T (sec) | 5.12 | 10.24 | 20.48 |

T1 (sec) | 0.50 | 1.50 | 2.00 |

T2 (sec) | 4.00 | 8.00 | 16.00 |

c | 2.0 | 1.0 | 0.7 |

d | 2.0 | 2.0 | 2.0 |

NFR | 100 | 200 | 300 |

fmax (Hz) | 12 | 15 | 15 |

δS | 0.40 | 0.40 | 0.40 |

The PGA values are obtained corresponding to hazard curves and produced for a specific region. As shown in Table 3, in this work the hazard curves presented by Möller et al. [13] are used. An artificial earthquake generated based on Eq. (19) is shown in Fig. 4.

Probability of exceedance | Recurrence interval | PGA (g) |

50% in 50 years | 73 | 0.27 |

10% in 50 years | 475 | 0.6 |

5% in 50 years | 975 | 0.8 |

## 6. The new discrete gravitational search algorithm

Based on the work presented by Khatibinia *et al*. [14], a new discrete gravitational search algorithm (DGSA) based on the standard gravitational search algorithm (GSA) is utilized to find the optimum designs in the RBDO procedure.

### 6.1. Gravitational search algorithm

GSA was introduced by Rashedi *et al*. [23] as a new stochastic population based search algorithm based on the law of gravity and mass interactions. In GSA, each agent of the population represents a potential solution of the optimization problem. The *i*th agent in *t*th iteration is associated with a position vector, *D* is dimension of the solution space. Based on [23], the mass of each agent is calculated after computing the current population fitness for a minimization problem as follows:

where *N*, *M*_{i}(*t*) and *i* at *t*th iteration, respectively; and

To compute the acceleration of an agent, total forces from a set of heavier masses applied to it should be considered based on the law of gravity (Eq. (23)). Afterwards, the next velocity of an agent is calculated as a fraction of its current velocity added to its acceleration (Eq. (24)). Then, its next position could be calculated using Equation (25):

where *i*th agent in dimension *d*, respectively. *i* and *j*. *kbest* is the set of first *k* agents with the best fitness value and biggest mass, which is a function of time, initialized to *N* and is decreased linearly to 1.0. *G* is a decreasing function of time.

### 6.2. The proposed discrete GSA

The binary GSA (BGSA) for solving discrete problem was developed by Rashedi *et al*. [25]. In the BGSA model, the positions of agents are indicated by one or zero. Hence, Eq. (24) can be used without any changes for updating the velocity of agents. Because of terms of probability, the velocity must be converted into interval [0, 1] by a transfer function,

Based on Eq. (26), a large computer memory is needed for the position of agents in BGSA. Also, coding and encoding of the position of agents is a time consuming process. In order to overcome the shortcomings of BGSA, a new DGSA based on the fundamental concept of the standard GSA with passive congregation is presented herein. The passive congregation strategy as perturbations operator can transfer information among agents in the optimization procedure [24]. Therefore, the search performance of DGSA can be improved using the passive congregation. To achieve this purpose, Khatibinia *et al*. [14] modified the velocity of agents based on the particle swarm optimizer with passive congregation (PSOPC) which is proposed by He *at al*. [24]. The modified velocity is expressed as follows:

where *i*th agent, the best previous position among all the agents and a randomly selected agent from the population, respectively. *r*_{i,1} and *r*_{i,2} are the random number in the interval [0, 1], respectively; *r*_{i,3} is the uniform random number in the interval (0, 1).

In DGSA, the scalar *A*={*A*_{1}, *A*_{2},...,*A*_{n}}. Hence, the position of agents is updated by the following equation instead of Eq. (26):

Therefore, the coding and encoding of the position of agents are omitted; and the position of agents is calculated as the integer value. The current position of agents may violate from the values of the set *A*. To avoid this problem, the current position of particles is limited as:

where

## 7. Approximating the structural seismic responses

MCS requires excessive computational cost for RBDO of structures in order to obtain an acceptable accuracy [11]. Because of the drawback, Khatibinia *et al*. [22] proposed a meta-model that predicted the mean and the standard deviation of the structural seismic responses in the RBDO process based on MCS. The meta-model consists of weighted least squares support vector machine (WLS-SVM) and wavelet kernel function, which is called WWLS-SVM [14, 22, 27].

### 7.1. Weighted least squares support vector machines

WLS-SVM was introduced as excellent machine learning algorithms in large-scale problems by Suykens *et al*. [26]. In fact, assigning weights to SVM as well as to the least squares version of SVM (LS-SVM) is resulted in more robust and precise prediction of functions [26].

WLS-SVM is described as the following optimization problem in primal weight space [26]:

Subject to the following equality constraints:

where

It is impossible to indirectly compute

Based on the Karush-Kuhn-Tucker (KKT) conditions, by eliminating

where

According to Mercer’s condition, a kernel

Consequently, the final WLS-SVM model for the prediction of functions becomes:

Weight

where *i*th primal diagonal element of inverse of matrix *D*, which is the matrix on the left-hand of the system of the linear Eq. (34) [59]. After that, weights

In WLS-SVM, Gaussian radial basis function (RBF) is frequently used as the kernel function, and it is expressed as:

where

Based upon the Suykens *et al*. [26], the WLS-SVM model based on RBF kernel function for predicting the output data is implemented using the following procedure:

**Step 1.** Assign training data *N=Ntot*.

**Step 2.** Find an optimum *Ntot* training data by 10-fold cross-validation, then solve linear system (Eq. (34)), and give the model (Eq. (37)).

**Step 3.** Sort the values

**Step 4.** Remove a small number of *M* points (typically 5% of the *N* points) which has the smallest values in the sorted

**Step 5.** Retain *N-M* points and set *N=N-M*.

**Step 6.** Go to 2 and retrain on the reduced training set.

### 7.2. The new meta-model-based wavelet kernel

Wavelets as kernel function have been introduced and developed in ANNs and SVMs [60-63]. It has been shown that wavelet kernel functions are superior to other kernel functions in the training ANN and SVM. Accordingly, the kernel function of WLS-SVM is substituted with a specific kind of wavelet functions proposed by Khatibinia et al. [22]. The meta-model based on wavelet kernel function is called WWLS-SVM. The cosine-Gaussian Morlet wavelet is used as the kernel function of WLS-SVM. The wavelet function is mathematically written as follows:

where *a* and *b* are the scale factor and the translation factor, respectively.

According to Zhang *et al*. [64], the translation-invariant wavelet kernels can be explained as follows:

where *n* is the number of samples;

Therefore, according to Eqs. (40) and (41), the wavelet kernel function of the cosine-Gaussian Morlet wavelet is given as follows:

The accuracy of WWLS-SVM prediction depends on the good selection of its parameters. Selecting appropriate values of these parameters is important for obtaining the excellent predicting performance. Hence, in this study, GSA is used to find the WWLS-SVM optimal parameter, *a* and *MAPE*) is used to evaluate the performance of the WWLS-SVM model as follows:

where *y* and

The WWLS-SVM training stage during GSA is performed according to the *k*-fold cross-validation (CV) [26]. Consequently, an optimization problem based on *MAPE* of the *k*-fold CV (*MAPE*_{k-fold CV}) is expressed as:

The converged solution is affected by the setting value of parameters in GSA. In this study, the values are selected based on the general recommendations by Rashedi *et al*. (2009). The flowchart of the meta-model based on WWLS-SVM [22] and GSA [14, 23] is shown in Fig. 5.

## 8. Predicting failure probability of structures

In the RBDO procedure, nonlinear time-history analysis of SSI system is used and it may be failed regarding a number of random structures [65]. In fact, a number of structures collapse and then lose their stability. Hence, these structures should be identified and eliminated from optimization process. For this purpose, a failure probability is considered as stability criterion. An efficient method is presented to train the failure probability with high performance [65]. This efficient method is consisted of a modified adaptive neuro fuzzy inference system (ANFIS) with a hybrid of fuzzy c-means (FCM) [66] and fuzzy particle swarm optimization (FPSO) [67]. To train the modified ANFIS, the input–output data are classified by a hybrid algorithm consisting of FCM-FPSO clustering. The optimum number of ANFIS fuzzy rules is determined by subtractive algorithm (SA).

### 8.1. Hybrid of FCM and FPSO for clustering

The FCM algorithm has been extensively studied and is known to converge to a local optimum in nonlinear problems. Moreover, the FPSO algorithm is robust method to increase the probability of achieving the global optimum in comparison with the FCM algorithm. The FCM algorithm is faster than the FPSO algorithm because it requires fewer function evaluations. This shortcoming of FPSO can be dealt with selecting an adequate initial swarm [65].

In this study, a hybrid clustering algorithm called FCM-FPSO is presented to use the merits of both FCM and FPSO algorithms and increase the procedure of convergence. In this way, the FCM algorithm finds an adequate initial swarm FPSO algorithm for commencing the FPSO. For this purpose, first, the FCM algorithm is utilized to find a preliminary optimization that shown by *N*_{FCM} times to create the some part of the initial swarm FPSO. Other particles of the initial swarm, i.e.

### 8.2. Modified ANFIS

An ANFIS model depends on the number of ANFIS fuzzy rules and membership functions. In other words, creating an ANFIS model with a minimum number of fuzzy rules can eliminate a well-known drawback. Therefore, for overcoming of this drawback, Khatibinia *et al*. [65] proposed a modified ANFIS to predict the probability of failure. In this model, the number of clusters, the cluster centers and membership grades are considered as parameters which optimized by subtractive algorithm (SA) and the hybrid of FCM-FPSO and used in FIS for tuning ANFIS. The algorithm flow of the proposed model is shown in Fig. 7. The proposed method is executed in the following steps [65]:

**Step 1.** SA finds the optimum number of the clusters (*nc)*.

**Step 2.** The hybrid FCM-FPSO algorithm partitions training data to *nc* clusters and determines membership grades each of clusters. This parameters is used for optimizing the center of rules and membership functions for the input and output data.

**Step 3.** The FIS structure with a minimum number of fuzzy rules and membership functions is generated by using SA and the hybrid FCM-FPSO algorithm. The FIS uses Gaussian function and linear function for membership function of input and output, respectively. These parameters are tuned for the ANFIS.

**Step 4.** The ANFIS is employed for training data. The ANFIS uses a hybrid learning algorithm to identify parameters of Sugeno-type fuzzy inference systems. a combination of the least-squares method and the back-propagation gradient descent method are applied for training FIS membership functions.

## 9. Numerical examples

In this work, two RC frame structures shown in Fig. 8 are selected as illustrative examples. Three layers of sand associated with material properties varying over its depth are considered as the soil under the frames. The depth of each soil layer and the entire width of soil domain are considered to be 10 m and 100 m, respectively. The soil is also assumed to have plane strain condition with a constant thickness of 5.0 m in proportion to the inter-frame distance. Inertia properties of the soil mesh are considered using lumped mass matrices modeling with soil mass density of 17 kN/m^{3} for all soil layers. The values of the dead and live loads are considered to be *DL*=5.884 N/mm^{2} (600 kg/m^{2}) and *LL*=1.961 N/mm^{2} (200 kg/m^{2}), respectively.

For vertical continuity on the dimensions along the height of a column, the section database of columns is divided into three types in the height of RC frame. Hence, a database shown in Table 4 is generated. Similarly, the section database of beams is divided into three types in the height of RC frame. Distribution of beam dimensions along the height of the frame is shown in Table 4. The diameter of longitudinal bars for beams and columns is laid between 12 mm and 32 mm in the databases.

The initial cost is calculated for *DI*_{overall} adjusting a Beta distribution. Then, the density function is assigned by the derivative of the cumulative distribution. The target values of reliability indices corresponding to the three performance levels (Table 1) considered in RBDO of RC structures are equal to 1.276, 2.326 and 2.697, respectively [52]. For RBDO of RC structures, the probability density function (PDF), the mean and standard deviation (SD) values for each random constitutive parameter are listed in Table 5.

The concrete material parameters shown in Table 5 are considered for the cover of column cross-sections. The strain corresponding to the peak strength, *f*_{uo,} for unconfined concrete are selected as 0.002 and 0.0, respectively. Shear-wave velocity, *V*_{s,} and friction angle, *a*_{g}, and the central frequency, *f*_{g}, for soil filter are considered as random variables in the RBDO process. The PGA value is also taken into account as follows [13, 22]:

where

Number | Column (h × h)a | Beam (h × b)a | ||||

Type 1 | Type 2 | Type 3 | Type 1 | Type 2 | Type 3 | |

1 | 65 | 55 | 45 | 55 | 50 | 45 |

2 | 60 | 50 | 40 | 55 | 50 | 45 |

3 | 55 | 50 | 40 | 55 | 50 | 45 |

4 | 55 | 45 | 40 | 55 | 50 | 45 |

5 | 55 | 45 | 35 | 50 | 45 | 40 |

6 | 50 | 45 | 40 | 50 | 45 | 40 |

7 | 55 | 45 | 35 | 50 | 45 | 40 |

8 | - | - | - | 50 | 45 | 35 |

The presented DGSA requires the user to specify several internal parameters that can affect convergence behavior at the search space. It is found that a population of 50 agents can be adequate. Higher values are not recommended, as this will increase significantly computation time in RBDO. In addition, different optimization runs are carried out for RBDO model in this study, so optimum designs are found by DGSA about 150 iterations. Due to the effect of decreasing gravity, the actual value of the gravitational constant, *G*(*t*), depends on the actual age of the universe. In this study, *G*(*t*) is considered as a linear decreasing function [23] in DGSA. Since the initial value of the gravitational constant is found to affect the optimization results significantly, the fixed value of *G*_{0} = 50 is utilized in this study. In order to consider the stochastic nature of the optimization process, ten independent optimization runs are performed and the best solution is considered as the final results.

### 9.1. Example 1: Six-storey RC frame

Six-storey RC frame is shown in Fig. 8(a). In the frame, the length of each bay and the height of stories are 5 m and 3 m, respectively. The members of the structure are divided into four groups for the columns C_{1}, C_{2}, C_{3}, C_{4} and four groups B_{1}, B_{2}, B_{3} and B_{4} for the beams. The groups of structural elements are presented in Fig. 8(a).

Material | Variable | PDF | Mean | SD |

Concrete | ||||

fco (Mpa) | Lognormal | 28 | 2 | |

Normal | 0.0035 | 0.00035 | ||

Steel | ||||

fy (Mpa) | Lognormal | 340 | 25 | |

E (Mpa) | Lognormal | 210000 | 8000 | |

H | Normal | 0.015 | 0.0015 | |

Soil | ||||

Vs,1 (m/s) | Normal | 375 | 10 | |

φ1 | Normal | 37.5 | 1.0 | |

Vs,2 (m/s) | Normal | 300 | 20 | |

φ2 | Normal | 37.5 | 1.0 | |

Vs,3 (m/s) | Normal | 200 | 10 | |

φ3 | Normal | 32.5 | 1.0 |

Variable | PDF | Mean | SD |

^{2}) | Lognormal | 300 | 100 |

Normal | - | 0.15 | |

Normal | 2.5 | 0.375 |

#### 9.1.1. Training and testing the meta-model

To predict the mean, *i*th seismic response during RBDO of the frame, the proposed meta-models are trained based on the generated database. The WWLS-SVM training during GSA is performed according to five-fold cross-validation. The lower and upper bonds of the parameters required in the optimization process are selected as

In order to validate the performance and accuracy of the proposed meta-model, relative root-mean-squared error, i.e. *RRMSE*, and *R*^{2} as the absolute fraction of variance, during testing the meta-model and WLS-SVM, are defined using the following equations:

Parameter | Mean | Standard deviation | ||||

umax | 375.73 | 2.892 | 4.274 | 284.62 | 1.647 | 4.048 |

DRmax | 386.04 | 3.804 | 6.347 | 314.42 | 2.104 | 4.702 |

DILmax | 400.48 | 3.615 | 5.895 | 373.67 | 2.548 | 3.692 |

DIoverall | 365.37 | 4.052 | 6.329 | 308.38 | 1.947 | 4.082 |

The smaller *RRMSE* and *MAPE* and the larger *R*^{2} are indicative of better performance generality. The comparison of WWLS-SVM and WLS-SVM, with respect to *MAPE*, *RRMSE*, and *R*^{2} is shown in Table 8.

Method | Statistical parameters | Mean | Standard deviation | ||||||

umax | DRmax | DILmax | DIoverall | umax | DRmax | DILmax | DIoverall | ||

WWLS-SVM | |||||||||

MAPE | 2.008 | 2.417 | 2.337 | 2.172 | 2.163 | 2.075 | 2.532 | 3.406 | |

RRMSE | 0.018 | 0.022 | 0.020 | 0.021 | 0.023 | 0.024 | 0.045 | 0.028 | |

R2 | 0.9999 | 0.9998 | 0.9998 | 0.9999 | 0.9999 | 0.9998 | 0.9998 | 0.9997 | |

WLS-SVM | |||||||||

MAPE | 5.328 | 6.302 | 5.392 | 5.862 | 5.386 | 6.017 | 6.007 | 5.737 | |

RRMSE | 0.057 | 0.074 | 0.059 | 0.061 | 0.079 | 0.094 | 0.091 | 0.064 | |

R2 | 0.9988 | 0.9979 | 0.9987 | 0.9987 | 0.9985 | 0.9984 | 0.9988 | 0.9988 |

As given in Table 8, the proposed meta-model trained for the mean and the standard deviation of seismic responses has proper performance generality. Thus, the approximating performance of the meta-model based on WWLS-SVM and GSA is better than the WLS-SVM with RBF kernel in predictive ability and precision.

#### 9.1.2. Results of RBDO

In this example, RBDO of the RC frame is performed using DGSA associated with WWLS-SVM-based MCS. In the reliability process, the reliability indices, ^{6} samples generated with the LHD method. The cross-section of beams and columns are selected from Types 2 and 3, which are shown in Table 4. The optimum designs of the RC frame are listed in Table 9. Furthermore, the optimal solutions of DGSA are also compared with those of BGSA in Table 9.

As shown in Table 9, the optimal solutions of DGSA are better than those of BGSA in terms of the total cost and the number of iterations. The minimum reliability index,

The convergence histories of the optimum objective function are shown in Fig. 9 for DGSA and BGSA models. As can be seen in Fig. 9, DGSA method is more efficient than BGSA method. Optimum designs are found by DGSA and BGSA in 4450 and 5900 required approximate analyses by the meta-model, respectively.

Element groups no. | DGSA | BGSA | ||

Cross-section | ρ (%) | Cross-section | ρ (%) | |

C_{1} | 55 | 2.67 | 55 | 3.23 |

C_{2} | 55 | 2.64 | 55 | 3.20 |

C_{3} | 45 | 2.56 | 45 | 2.60 |

C_{4} | 45 | 2.33 | 45 | 2.48 |

B_{1} | 50 | 1.98 | 50 | 2.38 |

B_{2} | 50 | 2.13 | 50 | 2.30 |

B_{3} | 45 | 1.81 | 45 | 1.88 |

B_{4} | 45 | 1.69 | 45 | 1.78 |

CIC (Euro) | 3448 | 3542 | ||

CRC (Euro) | 1080 | 1045 | ||

CTOT (Euro) | 4528 | 4587 | ||

Iterations | 89 | 118 | ||

1.4435 | 1.6901 | |||

2.4885 | 2.7342 | |||

2.7986 | 3.1384 |

### 9.2. Example 2: Nine-storey RC frame

Nine-storey RC frame is shown in Fig. 8(b). In the frame, the length of each bay and the height of stories are 5 m and 3 m, respectively. The members of the structure are divided into six groups for the columns and six groups for the beams. The groups of structural elements are presented in Fig. 8(b).

#### 9.2.1. Training and testing the meta-model

After training database using the presented WWLS-SVM optimal parameters of the meta-model associated with the mean and the standard deviation of seismic responses are shown in Table 10. Furthermore, the performance generality of the proposed meta-model and WLS-SVM is given in Table 10 in terms of *MAPE,* *RRMSE* and *R*^{2}.

Method | Statistical parameters | Mean | Standard deviation | ||||||

umax | DRmax | DILmax | DIoverall | umax | DRmax | DILmax | DIoverall | ||

WWLS-SVM | |||||||||

MAPE | 2.837 | 3.028 | 3.127 | 2.689 | 3.024 | 2.682 | 3.008 | 3.105 | |

RRMSE | 0.026 | 0.038 | 0.035 | 0.023 | 0.032 | 0.024 | 0.031 | 0.036 | |

R2 | 0.9999 | 0.9998 | 0.9998 | 0.9999 | 0.9998 | 0.9999 | 0.9998 | 0.9998 | |

WLS-SVM | |||||||||

MAPE | 5.538 | 6.346 | 6.483 | 6.006 | 6.305 | 5.396 | 5.843 | 5.579 | |

RRMSE | 0.094 | 0.138 | 0.162 | 0.105 | 0.1057 | 0.0987 | 0.0998 | 0.0924 | |

R2 | 0.9987 | 0.9985 | 0.9981 | 0.9988 | 0.9987 | 0.9988 | 0.9988 | 0.9988 |

The results of Table 10 demonstrate that the meta-model is better than the WLS-SVM method in terms of performance generality. Therefore, the meta-model is reliably employed to predict the necessary responses during the RBDO process.

#### 9.2.2. Results of RBDO

As the first example, in this example RBDO of the RC frame is performed using DGSA and BGSA associated with WWLS-SVM-based MCS. In this example, the cross-section of beams and columns are selected from Types 1, 2 and 3, which are shown in Table 4. The best optimum designs of the RC frame are listed in Table 11.

As revealed in Table 11, the optimal solutions of DGSA are better than those of BGSA in terms of the total cost and the number of iterations. The minimum reliability index,

The convergence histories of the optimum objective function are shown in Fig. 10 for DGSA and BGSA models. As can be seen in Fig. 10, DGSA method is more efficient than BGSA method. Optimum designs are found by DGSA and BGSA in 4150 and 6150 required approximate analyses by the meta-model, respectively.

Element groups No. | DGSA | BGSA | ||

Cross-section | ρ (%) | Cross-section | ρ (%) | |

C_{1} | 65 | 2.49 | 65 | 2.90 |

C_{2} | 60 | 2.50 | 60 | 3.08 |

C_{3} | 55 | 2.53 | 55 | 2.87 |

C_{4} | 50 | 2.49 | 50 | 2.79 |

C_{5} | 45 | 2.27 | 45 | 2.68 |

C_{6} | 40 | 2.29 | 40 | 2.89 |

B_{1} | 55 | 1.88 | 55 | 2.18 |

B_{2} | 55 | 1.80 | 55 | 1.98 |

B_{3} | 50 | 1.82 | 50 | 2.01 |

B_{4} | 50 | 1.87 | 50 | 1.98 |

B_{5} | 45 | 1.68 | 45 | 1.79 |

B_{6} | 45 | 1.66 | 45 | 1.85 |

CIC (Euro) | 5571 | 5736 | ||

CRC (Euro) | 1004 | 987 | ||

CTOT (Euro) | 6575 | 6723 | ||

Iterations | 83 | 123 | ||

1.4985 | 1.7101 | |||

2.4953 | 2.7996 | |||

2.8514 | 3.1837 |

## 10. Conclusions

In general, the optimum design of structures depends on a number of parameters that are inherently uncertain. Reliability-based design optimization (RBDO) has been employed as the only method that assesses the influence of uncertain parameters and balance both cost and safety of structures. To account for all necessary uncertain and random parameters in RBDO of RC structures and to achieve the realistic optimum design of RC structures, the uncertain material properties of soil and structure, as well as the characteristics of ground motions should be considered as random parameters. Furthermore, the realistic seismic responses of RC structures can be account by consideration of soil-structure interaction (SSI) effects. In this work, a new discrete gravitational search algorithm (DGSA) and a new meta-modeling framework were incorporated for RBDO of RC structures with Performance-Based Design (PBD) under seismic loading. The objective of the RBDO problem was to minimize the total cost whereas the deterministic constraints and the system reliability index corresponding to each of the performance levels should not exceed a specified target. Based on this study, the following conclusions can be drawn:

To reduce the computational effort and computational cost of the Monte-Carlo Simulation (MCS) method, a new meta-model based on a wavelet weighted least squares support vector machine (WWLS-SVM) and gravitational search algorithm (GSA) was utilized in the RBDO procedure. Therefore, the proposed meta-model, as a substitute for the nonlinear dynamic analysis of SSI system, can estimate the reliability index through MCS with a small computational cost.

The WWLS-SVM and kernel parameters were simultaneously optimized in the proposed meta-model in order to improve performance generality of WWLS-SVM. Numerical results of training and testing the meta-model indicated that performance generality of the meta-model was higher in comparison to WLS-SVM. Hence, the proposed meta-model can predict the nonlinear dynamic analysis of SSI system in terms of accuracy and flexibility.

The proposed DGSA was presented based on the standard GSA with passive congregation. The passive congregation strategy can be considered as perturbations operator in the optimization procedure. Therefore, the presented DGSA using the passive congregation can transfer information among agents avoiding local minima. Furthermore, the coding and encoding of the position of agents as a time consuming process is omitted in DGSA. To eliminate this drawback, the position of agents was calculated as the integer value. The optimum designs obtained by DGSA were compared with those produced by BGSA model. Numerical examples showed that the proposed DGSA can converge and reach the optimum design more quickly than the BGSA model.

Future extension of current research could include reducing the computations involved in the PBD by replacing MCS with the response surface method or the importance sampling technique. The constraints imposed on the objective function could be also treated as random quantities (see [68]).

## Acknowledgments

Our special thanks go to Dr. Eysa Salajegheh (Distinguished Professor of Structural Engineering) in Department of Civil Engineering at Shahid Bahonar University of Kerman, Iran, for his cooperation in this research work.