Research Article  Open Access
Mohammad Parsaeimaram, Congqi Fang, Xianqi Luo, Chandan Shakya, "Seismic Performance Evaluation of DoubleSkin SemiBaseIsolated Building Using Incremental Dynamic Analysis", Advances in Civil Engineering, vol. 2018, Article ID 2747836, 19 pages, 2018. https://doi.org/10.1155/2018/2747836
Seismic Performance Evaluation of DoubleSkin SemiBaseIsolated Building Using Incremental Dynamic Analysis
Abstract
Base isolation is a widely accepted earthquake damage prevention technique. This method decouples the superstructure from the base by putting a flexible layer under each column, thereby elongating the time period of structures. Semibase isolation (SBI) in doubleskin structures is an effective technique to reduce the dynamic responses of structures due to earthquake motions, by utilizing the isolation devices in part of it instead of the entire base. This study presents a doubleskin structure consisting of outer fixed base frames with shear walls that have been detached from the inner core with a minor gap. The inner core of the structure has been dissociated from the base using the elastomeric bearings. Seismic response of 10story doubleskin structure with inner isolated core was compared to that of inner fixed base core to consider the yielding and collapse probability of the structure using the incremental dynamic analyses (IDA). The results showed that the time period in SBI buildings can be adjusted with the use of coupling beams between the inner and outer frames. Also, the time period and interstory drift ratio are both reduced as more floors are given coupling beams. However, these coupling beams are the most effective at the topmost floors.
1. Introduction
Seismic isolation is a design concept to dissociate the superstructure from its foundations completely or partially. Base isolation is an option to conventional structural design with fixed base and can be costeffective for buildings located in areas where a highly ground motion is likely to occur [1]. Conceptually, isolation systems minimize the earthquake responses of the superstructure by separating the building from the base. By lengthening the time period of the building using isolation systems, forces transmitted to the structure above the base are significantly reduced. In fact, buildings with prolonged time periods are subjected to lower seismic force and divert energy of the ground motion. In particular, highenergy seismic ground motions are diverted at higher mode frequencies [2, 3]. Energy dissipation capacity and lateral flexibility are the basic characteristics of an elastomeric isolation. Flexibility increases the time period of the structure, thereby reducing the earthquake force that has to be absorbed by the isolation system, while energy absorbing capacity itself reduces the seismic energy and increases damping of the building [4, 5]. During this procedure, a large amount of transmitted earthquake motion is dissipated while an adequate stiffness of the isolation device has to be provided to protect the structural integrity.
Uplift and compression force result from the overturning action due to the lateral forces on the lateral force resisting system (LFRS). In some cases, the isolators may be located under the LFRS (e.g., shear walls and bracing). In particular, the connection of the isolator to the LFRS should be carefully considered with respect to the transfer of tensile and compressive forces to the isolator under the columns. Although the horizontal and vertical stiffness of the rubber bearing can significantly withstand the transmitted compressive forces, the tensile forces (uplift) can cause instability in the isolation system.
The effects of the vertical component of earthquake acceleration on a structure had received relatively less attention in the past. Zhou et al. [6] studied several vertical and 3D isolation systems on the modern nuclear facilities. Their study consisted of a series of case study analysis of the modern nuclear power plant model to examine the benefits and challenges associated with 3D isolation, compared to that of the horizontal isolation. This study compared the performance of horizontal isolators alone with that of isolators capable of damping vertical frequencies and found that isolators having vertical frequencies under 3 Hz can effectively reduce the vertical responses for the studied nuclear power plant model.
Takahashi et al. [7] proposed a 3D seismic base isolation system in order to apply vertical and horizontal earthquake motions to an actual threestory building. In their proposed 3D seismic isolation, they used a device consisting of laminated rubber bearings as a horizontal isolation device and air springs as a vertical isolation device. The study showed that the vertical and horizontal maximum absolute acceleration drastically decreases with the use of 3D seismic base isolators under each supporting point.
Furthermore, 3D seismic isolation models have been developed in the past, by adding limited vertical isolation systems to particular portions of a horizontally isolated structure [8–12], and in these studies, it was found that vertical isolation systems could better protect the structure against the vertical component of the response. However, the studies stated above were restricted to the vertical component of excitation, and the uplift force behavior was limited to the effectiveness of the vertical component of the earthquake acceleration. Accounting for the uplift force due to overturning moment could have resulted in a different response of the isolated structures when compared to the vertical component of excitation alone. As previously stated, isolated structures have a longer fundamental time period due to a flexible layer between the base and superstructure. In addition, the base isolation system is appropriate for structures that possess short to medium heights and whose dominant modes are within a specific frequency range. However, in highly seismic zones and in structures with longer period components, it is difficult to build components flexible enough for the energy due to ground motion to be reflected [13].
The presented study proposes a doubleskin structure called SBI building to address the aforementioned issues and to reduce the effectiveness of uplift in isolation devices due to LFRS. The proposed model consists of isolators in only the inner part of the structure, in which the LFRS is located only at the perimeter of the building with fixed base as outer frames, and in which the inner gravity columns (inner core) have been supported by isolators beneath the columns with sufficient separation from the outer frame. To avoid uplift in the isolation devices, LFRS is located in the outer frames with fixed bases. Furthermore, in this system, the inner core is connected to the outer frames in every floor level in order to increase the stiffness of the structure and to decrease fundamental time period.
Several research projects have been conducted to account for various cases of doubleskin facade (DSF) buildings in which the lateral loads have been assigned to and are resisted by the primary building structure, while the exterior skin is employed as a curtain wall. For such cases, a significant amount of studies were carried out by Bedon and Amadio [14–16] for the mitigation and enhancement of glass curtains under extreme loads based on the potential of special mechanical connectors inserted between the primary structure and exterior glazing facade. And also they investigated the behavior of a cablesupported glazing facade subjected to high and mediumlevel air blast loads [17, 18].
Recently, a variety of studies have focused on the properties and applications of doubleskin system regarding environmental conditions and energy saving [19–21], while little research exists on the seismic behavior of this system. Moon [22, 23] studied the structural capacities of the doubleskin facade (DSF) systems in tall buildings with fixed base subjected to dynamic wind loads. Furthermore, Pipitone et al. [24] employed several layouts of doubleskin facade (DSF) with mass dampers to reduce the seismic vibration in structures subjected to variety of earthquake excitations. They proposed an optimized design configuration of doubleskin facade (DSF) based on numerical analysis. Fu and Zhang [25] also proposed an integrated control system to combine doubleskin facade and mass dampers in buildings. They show that by motorizing the DFCs, energy efficiency can be improved, and their results indicate that mass damper systems can significantly reduce structural motions under earthquake excitation. In [26], the effect of distributed mass dampers in multistory doubleskin facade building under a set of 20 earthquake excitations was investigated. They employed both genetic algorithms (GA) and particle swarm optimization (PSO) techniques to find the optimal design parameters for doubleskin facade panels. It was concluded that doubleskin facade panels can reduce the seismic response up to 35% of the uncontrolled case.
The aforementioned literature reviews have been carried out with a primary focus on the environmental applications and energy aspects and as a vibration absorber under earthquake excitation of the DSF system. In fact, no substantial studies on DSF structures focusing on the application of base isolation in the DSF systems could be found in the approachable references.
Here, the seismic response of 3D nonlinear finite element analysis of a doubleskin semibaseisolated (SBI) multistory building, composed of fixed base reinforced concrete (RC) frames with shear walls as outer frames and an inner isolated core with a sufficient gap from the outer frames, is investigated and the results are compared to that of inner fixed base core, under bidirectional earthquake acceleration. Incremental dynamic analyses (IDA) is used to estimate the performance of structures under seismic loading using OpenSees [27], for the modeling and analyses. Fragility functions of SBI and fixed base on the doubleskin structures are generated and compared and the results obtained are discussed.
2. Structure and Model Description
The prototype reinforced concrete (RC) building presented herein is a 10story doubleskin structure with moment resisting connections. The typical floor plan and elevations of the structural models used in this study are illustrated in Figure 1. The design of all buildings satisfy UBC97 [28] and ACI 31811 [29] requirements. Six prototype buildings are assumed to be located in seismic zone 3 on the stiff soil () based on UBC97 [28] classification. For the purposes of this study, four doubleskin structures and two conventional structures with identical configurations were considered as described below and results were compared for better understanding.
The doubleskin structures in consideration can be split into two parts, namely, inner core and outer frames. Among these, the outer frames were considered to be fixed. Furthermore, the two aforementioned conventional structures consist of the inner core of the doubleskin structure which has been embedded with the inside of the exterior skin.
Each model is given a name, starting with the letters “IFBC” or “IIC” to represent inner fixed base core and inner isolated core, respectively, followed by numbers or with the word “All” to denote the floor number in the inner core that are connected to the outer frames with coupling beams as illustrated in Detail2 of Figure 1. Detail2 shows the connection of the beam element between the inner core and the outer frames.
Based on the above statements, six prototype models are defined as (IIC10th), (IFBC10th), (IICAll), (IFBCAll), (IIC), and IFBC. Configurations of all the models are explained in the following:(i)Model IIC10th stands for a doubleskin structure with inner isolated core in which the inner core has been connected to the outer frames at the tenthfloor level, as illustrated in Figure 1(b).(ii)Model IFBC10th stands for a doubleskin structure with inner fixed base core in which the inner core has been connected to the outer frames at the tenthfloor level, as illustrated in Figure 1(c).(iii)Model IICAll stands for a doubleskin structure with inner isolated core in which the inner core has been connected to the outer frames at all floor levels.(iv)Model IFBCAll stands for a doubleskin structure with inner fixed base core in which the inner core has been connected to the outer frames at all floor levels.(v)Model IIC stands for a structure with an inner isolated core only, which has no outer frames.(vi)Model IFBC stands for a structure with an inner fixed base core only, which has no outer frames.
Four RC shear walls are located in the outer frames in both longitudinal and transverse directions to withstand lateral earthquake forces. Moment resisting connections are used at the intersection of beams and columns. The cross sections of the beam elements in the inner core and outer frames are assumed to be equal in all stories. The plan view of the 10story square buildings with three 4.0 m bays is depicted in Figure 1(a). All stories have typical heights equal to 3.0 meters. The dead and live loads are assumed to be 6.50 kN/m^{2} and 2.50 kN/m^{2}, respectively. The compressive strength of concrete is assumed to be 25 MPa, and the tensile strength of rebar is 400 MPa. Concrete slabs are at story levels with 150 mm thickness. In the analysis and design processes, all lateral displacement limitations and strength requirements, as mandated, have been checked in reference with UBC97 [28] and ACI 31811 [29]. The details along with the crosssectional dimensions of the elements are presented in Table 1.

3. Description of Analytical Models
3.1. Material Properties and Elements Type
The analytical model behavior of the building is approximated by a continuous 3D finite element model using the OpenSees [27] nonlinear analysis program, as shown in Figure 2. Incremental dynamic time history analysis was conducted to determine the performance of the prototype buildings. To develop the analytical model inelastic forcebased formulations were used for the simulation of nonlinear beam and column elements. Beamcolumn elements were modeled as fiber sections with distributed material plasticity along the element length. Geometric transformation to account for the PDelta effect was used in the OpenSees model. Concrete01 was used for simulating uniaxial Kent–Scott–Park concrete material behavior taking into account unloading and reloading stiffness based on the work of Karsan and Jirsa [30], with zero tensile strength. The reinforcing bars were modeled using the uniaxial Steel01 to create a uniaxial bilinear steel material based on kinematic hardening defined by a nonlinear evolution equation.
3.2. Modeling of Base Isolation System
The analytical model of the lead rubber bearing (LRB) was generated using the elastomeric bearing element, which has a bilinear forcedeformation response envelope as implemented in the OpenSees program. This isolation system consists of layers of thin steel shim and rubber, to provide vertical stiffness and horizontal flexibility. For dissipating the seismic energy and reducing the bearing displacement, a lead core has been inserted into the holes in the elastomeric bearing [31, 32]. Sixteen equally sized LRBs were considered to be located under the inner core columns between the foundation and beamslab system, consisting of an orthogonal mesh of RC beams (650/650 mm) [33, 34]. The slabs at the base were identical to those of the other floors, as illustrated in Detail1 of Figure 1. For dynamic analysis, the mass of the structure was lumped at the beamcolumn joint at each floor level. A rigid diaphragm was considered at each floor level in order to obtain the same lateral displacement across the floor. Applying rigid diaphragm constraints in the OpenSees models also reduced modeling complexity and analysis time [35, 36]. The design of superstructure is based on the design basis earthquake (DBE), while isolation devices should be designed for displacement obtained from maximum considered earthquake (MCE) [37]. The mechanical properties of the isolators were set to comply with the recommendation of the UBC building code subjected to critical reaction force in inner core columns. Moreover, the first mechanical characteristics of the isolators can be identified by the manufacturer [38], and the preyield stiffness () can be either evaluated from available hysteresis loops obtained from elastomeric bearing tests or assigned as a factor of postyield stiffness () [39]. Three parameters were used to simulate the nonlinear forcedeformation behavior of the isolators through the bilinear hysteresis loop as shown in Figure 3(b), namely, (i) preyield stiffness () and (ii) postyield stiffness (), and characteristic strength (). is taken to be 10 times [39–41]. and were initially calculated using estimated deformation of the bearings and the effective period of the building. The characteristic strength is estimated from the hysteresis loop for the elastomeric bearings, and for LRB systems, is given bywhere and are the area and yield stress of the lead, respectively, and is assumed to be 10 MPa. The total effective stiffness () of the isolation system is generally estimated by the fundamental isolation period of the building (). In the case of isolated structures, base isolation is most appropriate when the isolation period was selected to be in the range between 1.5 and 3.0 s [42–44], which in this study was considered to be 2.0 s. The effective stiffness is expressed aswhere is the total mass of the building, while the postyield stiffness is given by the following formula:where is the design displacement. The compression stiffness of the bearing is determined by the axial stiffness and is given by the following formula:where is the crosssectional area of the bearing; , the total thickness of the rubber, is given by and taken as , where is maximum shear strain taken as 1.5, for all LRBs; is the compression modulus of the rubbersteel and is given by , where the shape factor of the bearing is assumed to take the value of 10 [39]; and is the shear modulus of the elastomeric bearing amounts to 0.80 MPa [39, 45]. The properties used in the design LRB bearings are presented in Figure 3. For this study, the effective damping of an isolator is assumed to be equal to 15%. The compression stiffness of the bearing was determined as axial stiffness and was set to 1301463.65 kN/m.
(a)
(b)
Base isolation requires a gap between the inner core and its surroundings to afford enough distance for the displacement of isolators. This gap should be appropriately preserved and must remain free through the entire life of the structure [46, 47]. The separation gap should be enough to allow the movement of isolation devices in the horizontal plane without colliding with its surroundings. In addition, special architectural details must be included to limit pounding during strong earthquakes. In order to evaluate the gap size, the total maximum displacement of the isolation system (), should be considered with regard to both the maximum translational displacement () and the torsional displacement. This relation is given bywhere is the short dimension of the structure measured perpendicular to the long dimension of the building () and is measured in the direction normal to the seismic loading from the center of rigidity of the isolation system and a corner. And is the eccentricity between the center of mass of the superstructure and the center of rigidity of the isolation system considering 5% accidental eccentricity, is the gravitational acceleration, is the seismic coefficient, and is a numerical coefficient related to the effective damping of the isolation system which is taken as 1.35, corresponding to damping coefficient of [28]. The analytical gap distance is shown in Detail1 of Figure 1.
3.3. Modeling of Multilayer Shell Element for Shear Wall
The RC shear walls were simulated using the multilayered shell element ShellMITC4, proposed by Dvorkin et al. [48] based on the theory of mixed interpolation of tensorial components (MITC) with four nodes, which was already provided in OpenSees. The multilayered shell element is based on the principles of composite material mechanics and has several layers, of different thicknesses, with different multidimensional concrete [49, 50] as shown in Figure 4. In this study, the section of solid shear walls was modeled with eightlayer shell elements, including the multidimensional concrete, reinforcement material, with corresponding thicknesses.
(a)
(b)
4. Ground Motion Ensemble
Nonlinear response history analysis was carried out on 3D OpenSees models of sample buildings. In the analytical models, 22 component pairs of horizontal ground motions were uniformly applied at the base of the buildings. Ground motion records selected were obtained from the next generation attenuation (NGA) project database developed by Pacific Earthquake Engineering Research (PEER) [51] The records had a moment magnitude between 6.5 and 6.9, and the closest distance to fault rupture ranged from 13.9 to 37.7 km, all of which were recorded on stiff soil. For brevity, only the acceleration response spectrum (5% damped) for the selected set of ground motion records are presented in Figure 5. Characteristics of the selected records set are tabulated in Table 2 including the directional component (), which shows the direction of seismic input in E and N directions, where the letters E and N indicate the global direction of structures in X and Y directions, respectively, as illustrated in Figure 1(a), moment magnitude (M), closest distance to fault rupture (R), and peak ground acceleration (PGA) values. Ground motions are scaled to represent a range of earthquake intensities up to collapse level ground motions.

The selected ground motions in this study were scaled using the firstmode vibration period with 5% critical damping , with each ground motion scaled such that the scale factor is used to match the spectral acceleration at the firstmode period of each structure to the design spectral acceleration [37]. Once the scaled records were generated, each record was incrementally scaled in multiple levels by an increment of 0.05 g in spectral acceleration to cover the entire range of structural response, from yielding zone to dynamic instability.
A value of 5% Rayleigh damping, mass, and stiffness proportional damping, was used in the dynamic analysis [52, 53]. The gravity load combinations for nonlinear analysis were different from design gravity load combinations, and effective seismic weight for each floor was estimated using the load combination 1.05 DL + 0.25 LL [54, 55] and was converted to lumped masses at beamcolumn joint in the story levels for dynamic analysis, where DL is the nominal and superimposed dead load of the structure, and LL is the nominal live load.
5. Incremental Dynamic Analysis (IDA)
IDA is a parametric analysis technique that has been recently developed to carefully estimate the seismic performance of structures using various nonlinear dynamic analyses under several scaled ground motion records. The method has been described in detail by Vamvatsikos and Cornell [56]. IDA requires conducting a series of nonlinear dynamic time history analyses for each record for producing one or more curves of engineering demand parameter (EDP) as a damage measure (DM) versus intensity measure (IM) curve under the effect of several scaled ground motions. The ground motions are characterized by the intensity measure (IM), which can be expressed in terms of peak ground acceleration (PGA), peak ground velocity (PGV), or spectral acceleration (). Damage measure (DM) can be defined and related to performance limit states of the structure corresponding to various damage levels. For the present study, spectral acceleration at the firstmode period of vibration with 5% damping was used as the seismic intensity measure (IM), and the maximum interstory drift ratio was considered as damage measure (DM). Maximum interstory drift ratio is defined as the difference between the maximum displacements of floors in a story, divided by the story height (e.g., maximum, overall stories, and peak interstory drift ratio ).
In order to consider the limiting values of damage levels for the derivation of seismic fragility curves, limit states should be described on the IDA curves. Some studies considered different damage levels. For example, three limit states have been suggested in FEMA 350 [57], such as immediate occupancy (IO), collapse prevention (CP), and global dynamic instability collapse. The first limit state (IO) is defined at maximum interstory drift ratio , whereas the collapse prevention (CP) is named on the median IDA curve using the 20% tangent slope scale. Finally, global dynamic instability occurs by increasing levels of IM at constant steps, when the flatline on the IDA curve is reached.
DumovaJavanoska [58] and Altay et al. [59] defined different limit states, such as slight, moderate, extensive, and collapse as stated in HAZUSMH [60], whereas Kirçil and Polat [61] specified two different damage levels, that is, yield point at the end of the linear path and collapse which occurs in the nonconverging run. In this study, yielding and collapse were considered as essential damage levels for which yield and collapse points are displayed in Figure 6(f).
(a)
(b)
(c)
(d)
(e)
(f)
(g)
The yield structural capacity is defined by the point where the IDA curve starts to decline from the linear path, (i.e., the relationship between the maximum interstory drift ratio and spectral acceleration at the yield point of IDA curve is linear). Moreover, elastic IDA region can be defined as the region where the initial stiffness remains constant with increase in spectral acceleration with related slope remaining constant until a specific point (yield point) is reached. This is important since both slope and stiffness start changing with the increasing spectral acceleration. This attests to the elastic behavior of the structure. To identify the collapse capacity of structure, the selected earthquake records first need to be scaled with low values of IM and gradually increased to several higher IM levels until structural collapse occurs as a result of a numerical nonconverging run. In this study, the state of collapse capacity of the structure at the nonconverging run was considered as the limit state for failure.
The seismic performance of the six prototype buildings under IDA analysis was compared in terms of various performance criteria. Moreover, results from the IDA analysis were plotted as IDA curves in the IMDM plane. In each plot, locating the maximum drift observed in an analysis gives one point in the IM versus DM ( versus ) domain, and with increasing the values of IM, other points on each IDA curve are generated. Thus, IDA curves for six analytical models in the XX direction are shown in Figure 6. Therefore, IDA curves exhibit an extensive range of structural behavior and large recordtorecord variability. In such case, it is essential to summarize such dataset and quantify the randomness.
The results from IDA curves in the limitstate capacities (yielding and collapse) are tabulated in Tables 3–8 and also can be easily summarized into 16%, 50%, and 84% fractile values of DM and IM capacity for each limit state [62]. Summarized IDA curves were developed to estimate the variation of structural capacities as shown in Table 9. Figures 6(c), 6(d), and Table 9 show approximately similar performance for structures IFBCAll and IICAll for all the 16%, 50%, and 84% fractile IDA curves. The structural capacities of all prototype buildings at the collapse damage level were estimated through IDA on the median (50%) percentile, which defines the capacity of structures as the last IDA point in the nonconverging run. The median collapse and yielding capacities for models IFBCAll and IICAll are 2.788, 2.684 and 0.254, 0.264 with corresponding maximum interstory drift ratio () of 0.01275, 0.0115 and 0.001, 0.00115, respectively, with a minor difference of 0.104 in terms of Sa, which means that the SBI with inner isolated core, as connected to the outer frames at all floor levels, has somewhat similar performance compared to that of IFBCAll. In addition, based on an eigenvalue analysis as stated in Section 7, the time period of both IFBCAll () and IICAll () also has a minor difference in the firstmode vibration period. It can be explained by the fact that the behavior of the structures cannot be significantly changed by using the base isolation devices in the inner columns in the cases where the inner structure is connected to the outer fixed base frames at all floor levels (e.g., IICAll).







Using the base isolation system in the structures with long time period cannot be significantly effective. Whereas, as illustrated in Table 9, results from analysis of the isolated structure (e.g., IIC) embedded into the outer stiff frames with fixed bases while connected to the outer frames only at the topmost floor (e.g., IIC10th) show that time period and interstory drift ratio decrease from 0.87 s to 0.569 s and 0.20 to 0.0149, respectively, even in comparison to the IFBC10th structure. Comparing the elastic behavior of the SBI structure, Figures 6(b) and 6(e) show that IIC10th and IIC have a similar elastic behavior with a short elastic region and both reach yielding capacity quite early. From Table 9, we can see that in terms of collapse capacity, the IFBC model when compared to those with embedded outer stiff frames (e.g., IFBCAll) displayed a marked reduction of 16.42% and 94.40% in spectral acceleration and interstory drift ratio, respectively.
6. Fragility Analysis Procedure
Once the IDA is completed, fragility curves can be obtained. Based on the results of the distribution of nonlinear dynamic analysis, fragility curves can be employed to quantitatively estimate the collapse capacity of the structures. A fragility curve relies on a graphical relationship of the probability of exceeding specific damage levels in terms of yielding and collapse due to earthquake as a function of given intensity. It is generally expressed in the form of twoparameter lognormal distribution function. Based on this statement, the cumulative probability of exceeding damage levels modified on the IM is expressed aswhere is the standard normal distribution, IM is the intensity measure of the earthquake in term of , is the damage level, and and are the mean and standard deviation of . The mean and standard deviation of values corresponding to the prescribed limit states are determined based on a regression analysis of the IDA results for each structural model [63]. Table 10 summarizes the mean and standard deviation of lognormal distributed intensity measure in terms of for each damage level obtained using the median collapse intensity and lognormal distributed ground motion ().

The fragility curves of the model structures constructed for yielding and collapse together are depicted in Figure 7. From the yield probability curves in Figure 7(a), it can be observed that the movement of yielding curves from right to left could be attributed to the reduction in yielding capacity of structures. Generally, flatter curves (e.g., IFBC10th, IFBCAll, and IICAll) exhibit a longer elastic range, whereas when the slope of the yielding curves continually increase (e.g., IIC10th, IIC, and IFBC), elastic range becomes short and maximum yielding capacity is reached early. As one would expect, the yielding phenomenon is generally higher when the inner core is completely connected to the outer frame in more floor levels. This means that the elastic range in model IIC10th can be changed from higher to lower values by connecting the inner core in some floor to the outer frames. As mentioned earlier, Figures 6(b) and 6(e) indicate that the elastic behavior of IIC10th and IIC are very close to each other.
(a)
(b)
Figure 7(b) depicts the collapse probability of sample buildings in terms of . From this figure, it is evident that the IIC structure has significantly lower collapse probability compared to that of IFBC10th, IIC10th, IFBCAll, IFBC, or IICAll. Moreover, fragility curves for IFBCAll and IICAll indicate that when the inner isolated core is connected to the outer frame at all floor levels, the behavior of the structure is very close to that of IFBCAll structure. Moreover, Figure 7(b) shows that collapse probability of IIC10th is higher, and it continues to decrease with the increase in the IM. This pattern can be numerically simulated at any specified higher IM levels. For example, given the spectral acceleration in (also indicated in Figure 7(b)), the corresponding collapse probability are 18.25%, 56.37%, 67.30%, 94.26%, 70.92%, and 41.58% for the IIC, IFBCAll, IICAll, IFBC10th, IIC10th, and IFBC, respectively. From this analysis, it can be concluded that the IIC10th structure has reduced collapse probability at higher IM ranges of approximately 25% when compared to that of IFBC10th. Furthermore, Figure 7(b) shows that all prototype models, in contrast to model IIC10th, have lower collapse probability. Moreover, models IFBCAll and IICAll with all connected floors have approximately identical performance pertaining to collapse phenomenon. Therefore, in doubleskin SBI structures, the probability of collapse can be reduced by connecting more floor levels to the outer frames.
7. Seismic Responses of Prototype Buildings
An eigenvalue analysis was conducted to better understand the response of the prototype buildings. This analysis elicits the calculation of time periods and vibration modes of the structures. To validate the consistency of the analytical models, structural periods and vibration modes of sample models simulated in ETABS [39] and OpenSees were compared. Therefore, in this study, time period of first three modes and vibration modes built in OpenSees and ETABS were compared, as shown in Tables 11–13. Results from Tables 11–13 are compared in this section to have a better understanding and also to corroborate the capabilities of the SBI system. Table 13 shows the time period of first three modes for all prototype models. For brevity, only the first three modes of vibration for the model inner fixed base core (IFBC) is tabulated in Tables 11 and 12. Results indicate that the time period and vibration modes obtained from OpenSees program matched very well with those of ETABS software.



Table 13 indicates that the time period in doubleskin SBI structure can be changed by increasing the number of connected floor levels. Except IIC10th, all other kinds of doubleskin structures (e.g., IFBCAll, IICAll, IIC10th, and IFBC10th) have a lower time period. This indicates that as the number of connected floors increases, the time period significantly decreases. Based on this statement, it can be also noted that when an isolated structure with long time period is encased and embedded inside a stiff outer frame to create a doubleskin SBI structure, the time period needs to be limited to a reasonable range of 1.5–3.0 s [42–44], by using more coupling beams in any floor level. It means that the dynamic characteristics of SBI structure such as time period can be drastically adjustable.
To evaluate the effectiveness of LRB isolation system in the SBI building, nonlinear dynamic analyses were performed with selected ground motions mentioned above. The seismic responses of prototype buildings were compared in terms of maximum absolute roof acceleration, maximum absolute roof displacement, and maximum absolute base shear force. The resulting curves extracted from a singlerecord IDA study subjected to Northrige01, 1994 (no. 14) ground motion data were compared for each model. The no. 14 ground motion possess 0.621 PGA with duration of 24 (s). Figures 8 and 9 show the predicted time history of the earthquake responses in the nonlinear time history analysis on six comparative models subjected to no. 14 ground motion. For better understanding, the maximum absolute responses such as accelerations, displacement, and base shear force are tabulated in Table 14.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)
(a)
(b)
(c)
(d)
(e)
(f)

As stated in previous sections, seismic behavior of model IFBCAll and that of model IICAll are very close to each other. Furthermore, results from this section significantly corroborate the above statement. From Table 14, maximum absolute floor acceleration and displacement for the IFBCAll are 17.95 (g) and 0.128 (m), respectively, which is somewhat close to the IICAll with maximum absolute floor acceleration and displacement of 17.80 (g) and 0.126 (m), respectively.
By comparing the model IIC with other five models, it is deduced that the maximum absolute roof displacement experienced in the IIC (0.19 m) is larger than the maximum absolute roof displacement in other five. Thus, embedding model IIC in the outer frames (as in IIC10th and IICAll), roof displacement is decreased by approximately 32%, while roof acceleration in model IIC (8.66 (g)) is drastically reduced.
As can be seen from Table 14, the roof acceleration and roof displacement of model IIC10th were (resp.) 0.70 and 0.72 times that of model IFBC10th, which means that the base isolation devices can drastically reduce earthquake responses of building.
Figure 9 shows base sheardisplacement hysteresis curves. The total base shear forces are estimated by summing the horizontal reaction of the force developed in the cross section of column bases. The model IIC showed lower total base shear forces compared to those of the other four models. Figure 9(b) and Table 14 show that IIC10th has slightly better performance in reducing the seismic responses of building in comparison to models IFBC10th, IFBCAll, IICAll, and IFBC.
In general, the purpose of base isolation systems is to minimize the damage of the structural and nonstructural elements due to ground shaking. As seen in this section, isolated structures (e.g., model IIC) has the largest roof displacement when compared with other prototype models, with large displacement causing damage in the structural elements. Therefore, using SBI systems in this study and encasing the isolated structures with stiff outer frames resulted in model IIC10th showing substantially reduced displacement, while acceleration and base shear could not be reduced the increase in the stiffness of the structure caused by connecting the topmost floor to the outer frames, when compared to that of the IIC model. However, SBI system and specifically model IIC10th reduce both the acceleration and base shear in comparison to the IFBC10th, IFBCAll, IICAll, and IFBC with comparably similar displacements. It is clear that when using the coupling beams in more floor levels (e.g., IICAll), the seismic behavior of SBI structures becomes closer to that of the structure with entire fixed base (e.g., IFBCAll). It means that the effectiveness of using the base isolation devices in doubleskin SBI structure when all floor levels are connected to the outer frames is somewhat trivial.
It is therefore evident that in an SBI system, a model configuration between models IICAll and IIC10th, for example, with two connected floor levels, can be appropriate.
As mentioned earlier, due to the overturning action due to lateral forces, a tensile force (uplift) arises in the lateral force resisting system (LFRS). In some cases, the uplift due to the isolators located under the LFRS may cause instabilities in the isolation systems. To avoid this uplift force in the isolation devices, doubleskin SBI systems can have their entire LFRS removed from the inner isolated structure into the outer frames with fixed bases.
8. Conclusions
A numerical model has been presented for simulating the behavior of 3D 10story doubleskin SBI structure, using incremental dynamic analysis (IDA) subjected to 22 pairs of bidirectional scaled ground motion records. The results of the analysis of the sample structures were compared to those of the structures with added coupling beam at every floor level to investigate the seismic behavior of doubleskin SBI buildings. Timehistory analysis was also conducted using representative Northrige01, 1994 (no. 14) ground motion data in terms of maximum absolute roof acceleration, maximum absolute roof displacement, and maximum absolute base shear force. Based on the results obtained from the IDA and fragility and timehistory analysis, the following conclusions are drawn:(1)By comparing values from the elastic region obtained from the IDA, it was found that both the IIC10th and IIC were subject to a lower elastic region and reached yielding capacity much earlier than that of IFBC10th, IFBCAll, and IICAll.(2)Results showed that for models with inner isolated core connected to the outer frames at all floor levels (model IICAll), the earthquake responses of structures such as collapse and yielding capacity, interstory drift ratio , and also spectral acceleration are very similar to the model with the inner fixed base core (model IFBCAll). This implies that the benefits of base isolation in SBI system are reduced if the inner isolated core is connected to the outer frames at all floor levels.(3)In contrast to IIC (with time period and interstory drift ratio ), IIC10th (with and ) has both shorter time period and reduced interstory drift ratio, it can thus be proposed that isolated structures that have long time period and large displacement when encased and embedded inside a stiff outer frame result in reduced responses. It means that doubleskin SBI system can be significantly more effective than isolated structures with long time period and large lateral displacement.(4)Results from the collapse probability curves show that except model IIC10th, all prototype models have approximately lower collapse level, and models IFBCAll and IICAll with all floors connected with the outer frame have approximately identical performance in relation to collapse. It can be observed that for doubleskin SBI system with increased number of connected floor levels, the probability of collapse decreases. Therefore, with respect to model IIC10th, it is better to increase the number of connected floors, while with respect to model IICAll it is better to decrease the number of connected floors.(5)Eigenvalue analysis shows that the time period in doubleskin SBI structure can be changed by increasing the number of connected floors. Except IIC10th, all other kinds of doubleskin structures (e.g., IFBCAll, IICAll, IIC10th, and IFBC10th) have lower time period. This indicates that, as the number of connected floors increase, the time period decreases. It means that by coupling beams in any floor levels, dynamic behavior of SBI system can be adjusted. This expression indicates that the dynamic behavior of SBI structure such as time period can be drastically adjusted by adding or removing coupling beams at various floor levels.(6)Timehistory analysis with representative Northrige01, 1994 (no. 14) ground motion data was used to compare the seismic responses of prototype buildings. Results showed that the response of IFBCAll subjected to Northrige01 is very close to IICAll in terms of maximum absolute roof acceleration and maximum absolute roof displacement. Finally, it is concluded that in SBI system, when inner core is connected to the outer frame only at the top (e.g., model IIC10th), results show that model IIC10th has substantially reduced displacement, while acceleration and base shear could not be reduced with the increase of the stiffness of the structure caused by connecting the topmost floor to the outer frames, when compared to that of the IIC model. In contrast to models IFBC10th, IFBCAll, IICAll, and IFBC, model IIC10th has both acceleration and base shear reduced.(7)According to the results obtained from model IICAll, base isolation devices in doubleskin SBI structure with all connected floor levels have somewhat trivial efficiency. Moreover, results show that the model IIC10th exhibits a higher probability of collapse. Therefore, a model configuration of SBI system between models IICAll and IIC10th with for example two connected floor levels can be appropriate.(8)Doubleskin SBI system was found to be effective in reducing the time period and displacement in conventional isolated structure with long time period and large displacement. Results reveal that by increasing or reducing the number of connected floors in the SBI system, dynamic behaviors of the SBI system can be changed. Therefore, an adjustable structure can be created by using coupling beams at every floor level. As stated earlier, the SBI system can provide significant structural integrity by transferring the lateral force resisting system (LFRS) from the inner core (isolated structure) to the outer frames with fixed base by eliminating the tensile forces (uplift) from LFRS which can cause instability in isolation devices.
Conflicts of Interest
The authors declare no conflicts of interest.
Authors’ Contributions
All the authors made contributions to writing, revising, and proofreading the manuscript.
Acknowledgments
The authors wish to acknowledge the supports of the National Natural Science Foundation of China (NSFC) (Grant nos. 51048001 and 51178264) and National Railway Ministry of China (Grant no. J2011G003).
References
 A. K. Chopra, Dynamics of Structures: Theory and Applications to Earthquake Engineering, PrenticeHall, Upper Saddle River, NJ, USA, 2001.
 B. S. Taranath, Wind and Earthquake Resistant Buildings: Structural Analysis and Design, CRC Press, Boca Raton, FL, USA, 2004.
 T. K. Datta, Seismic Analysis of Structures, John Wiley & Sons, Hoboken, NJ, USA, 2010.
 R. S. Jangid and T. K. Datta, “Nonlinear response of torsionally coupled base isolated structure,” Journal of Structural Engineering, vol. 120, no. 1, pp. 1–22, 1994. View at: Publisher Site  Google Scholar
 S. Nagarajaiah, A. M. Reinhorn, and M. C. Constantinou, “Torsion in baseisolated structures with elastomeric isolation systems,” Journal of Structural Engineering, vol. 119, no. 10, pp. 2932–2951, 1993. View at: Publisher Site  Google Scholar
 Z. Zhou, J. Wong, and S. Mahin, “Potentiality of using vertical and threedimensional isolation systems in nuclear structures,” Nuclear Engineering and Technology, vol. 48, no. 5, pp. 1237–1251, 2016. View at: Publisher Site  Google Scholar
 O. Takahashi, H. Aida, J. Suhara, R. Matsumoto, Y. Tsuyuki, and T. Fujita, “Construction of civil building using three dimensional seismic isolation system: part 1, design of building using three dimensional seismic isolation system,” in Proceedings of the 14th World Conference on Earthquake Engineering, Beijing, China, October 2008. View at: Google Scholar
 M. Morishita, S. Kitamura, S. Moro, Y. Kamishima, and T. Somaki, “Study on 3dimensional seismic isolation system for next generation nuclear power plant–vertical component isolation system with coned disk spring,” in Proceedings of the 13th World Conference on Earthquake Engineering, Vancouver, BC, Canada, August 2004. View at: Google Scholar
 S. Ogiso, “Development of 3D seismic isolator using metallic bellows,” in Transactions of the 17th International Conference on Structural Mechanics in Reactor Technology, Prague, Czech Republic, August 2003. View at: Google Scholar
 K. Inoue, M. Fushimi, S. Moro, M. Morishita, S. Kitamura, and T. Fujita, “Development of threedimensional seismic isolation system for next generation nuclear power plant,” in Proceedings of the 13th World Conference on Earthquake Engineering, Vancouver, BC, Canada, August 2004. View at: Google Scholar
 T. Somaki, Development of 3Dimensional Base Isolation System for Nuclear Power Plants, SMIRT, Washington, DC, USA, 2001.
 J. Wong, Z. Zhou, and S. Mahin, Seismic Isolation of Nuclear Power Plants, Electric Power Research Institute Inc., Palo Alto, CA, USA, 2013, Tech. Rep. 3002000561.
 J. Marko, Influence of Damping Systems on Building Structures Subject to Seismic Effects, School of Civil Engineering, Physical Infrastructure Center, Queensland University of Technology, Brisbane, Australia, 2006, Ph.D. thesis.
 C. Bedon and C. Amadio, “Numerical assessment of vibration control systems for multihazard design and mitigation of glass curtain walls,” Journal of Building Engineering, vol. 15, pp. 1–13, 2018. View at: Publisher Site  Google Scholar
 C. Bedon and C. Amadio, “Enhancement of the seismic performance of multistorey buildings by means of dissipative glazing curtain walls,” Engineering Structures, vol. 152, pp. 320–334, 2017. View at: Publisher Site  Google Scholar
 C. Bedon and C. Amadio, “Passive control systems for the blast enhancement of glazing curtain walls under explosive loads,” Open Civil Engineering Journal, vol. 11, no. 1, pp. 396–419, 2017. View at: Publisher Site  Google Scholar
 C. Amadio and C. Bedon, “Viscoelastic spider connectors for the mitigation of cablesupported façades subjected to air blast loading,” Engineering Structures, vol. 42, pp. 190–200, 2012. View at: Publisher Site  Google Scholar
 C. Bedon and C. Amadio, “Exploratory numerical analysis of twoway straight cablenet façades subjected to air blast loads,” Engineering Structures, vol. 79, pp. 276–289, 2014. View at: Publisher Site  Google Scholar
 P. C. Wong, D. Prasad, and M. Behnia, “A new type of doubleskin façade configuration for the hot and humid climate,” Energy and Buildings, vol. 40, no. 10, pp. 1941–1945, 2008. View at: Publisher Site  Google Scholar
 S. Chou, K. Chua, and J. Ho, “A study on the effects of double skin façades on the energy management in buildings,” Energy Conversion and Management, vol. 50, no. 9, pp. 2275–2281, 2009. View at: Publisher Site  Google Scholar
 D.W. Kim and C.S. Park, “Difficulties and limitations in performance simulation of a double skin façade with EnergyPlus,” Energy and Buildings, vol. 43, no. 12, pp. 3635–3645, 2011. View at: Publisher Site  Google Scholar
 K. S. Moon, “Tall building motion control using double skin façades,” Journal of Architectural Engineering, vol. 15, no. 3, pp. 84–90, 2009. View at: Publisher Site  Google Scholar
 K. S. Moon, “Structural design of double skin facades as damping devices for tall buildings,” Procedia Engineering, vol. 14, pp. 1351–1358, 2011. View at: Publisher Site  Google Scholar
 G. Pipitone, G. Barone, and A. Palmeri, “Optimal design of doubleskin façades as vibration absorbers,” Structural Control and Health Monitoring, vol. 25, no. 2, p. e2086, 2018. View at: Publisher Site  Google Scholar
 T. S. Fu and R. Zhang, “Integrating doubleskin façades and mass dampers for structural safety and energy efficiency,” Journal of Architectural Engineering, vol. 22, no. 4, p. 04016014, 2016. View at: Publisher Site  Google Scholar
 G. Barone, A. Palmeri, and A. Khetawat, “Passive control of building structures using doubleskin façades as vibration absorbers,” in Proceedings of the Fifteenth International Conference on Civil, Structural and Environmental Engineering Computing, Prague, Czech Republic, September 2015. View at: Google Scholar
 S. Mazzoni, F. McKenna, M. H. Scott et al., OpenSees Command Language Manual, University of California, Berkley, MI, USA, 2007.
 Uniform Building Code, (UBC97), Part. 2, 1997.
 ACI Committee, AC Institute, and I.O.f. Standardization, Building Code Requirements for Structural Concrete (ACI 31811) and Commentary, American Concrete Institute, Farmington Hills, MI, USA, 2011.
 I. D. Karsan and J. O. Jirsa, “Behavior of concrete under compressive loadings,” Journal of the Structural Division, vol. 95, no. 12, pp. 2543–2563, 1969. View at: Google Scholar
 C. Providakis, “Effect of LRB isolators and supplemental viscous dampers on seismic isolated buildings under nearfault excitations,” Engineering Structures, vol. 30, no. 5, pp. 1187–1198, 2008. View at: Publisher Site  Google Scholar
 R. S. Jangid, “Optimum lead–rubber isolation bearings for nearfault motions,” Engineering Structures, vol. 29, no. 10, pp. 2503–2513, 2007. View at: Publisher Site  Google Scholar
 R. Han, Y. Li, and J. van de Lindt, “Seismic risk of base isolated nonductile reinforced concrete buildings considering uncertainties and mainshock–aftershock sequences,” Structural Safety, vol. 50, pp. 39–56, 2014. View at: Publisher Site  Google Scholar
 V. Kilar and D. Koren, “Seismic behaviour of asymmetric base isolated structures with various distributions of isolators,” Engineering Structures, vol. 31, no. 4, pp. 910–921, 2009. View at: Publisher Site  Google Scholar
 E. Kalkan and A. K. Chopra, “Evaluation of modal pushover–based scaling of one component of ground motion: tall buildings,” Earthquake Spectra, vol. 28, no. 4, pp. 1469–1493, 2012. View at: Publisher Site  Google Scholar
 M. Fairhurst, X. Zhang, and T. Tannert, “Nonlinear dynamic analysis of a novel timber–steel hybrid system,” in Proceedings of the World Conference on Timber Engineering, Quebec City, QC, Canada, August 2014. View at: Google Scholar
 A. Masroor and G. Mosqueda, Seismic Response of Base Isolated Buildings Considering Pounding to Moat Walls, MCEER, 2013.
 C. Bedon and A. Morassi, “Dynamic testing and parameter identification of a baseisolated bridge,” Engineering Structures, vol. 60, pp. 85–99, 2014. View at: Publisher Site  Google Scholar
 F. Naeim and J. M. Kelly, Design of Seismic Isolated Structures: from Theory to Practice, John Wiley & Sons, New York, NY, USA, 1999.
 T. K. Datta, Seismic Analysis of Structures, Indian Institute of Technology, John Wiley & Sons (Asia) Pte Ltd., Delhi, India, 1994.
 O. Gomase and S. Bakre, “Performance of nonlinear elastomeric baseisolated building structure,” International Journal of Civil and Structural Engineering, vol. 2, no. 1, p. 280, 2011. View at: Google Scholar
 A. TenaColunga and C. ZambranaRojas, “Dynamic torsional amplifications of baseisolated structures with an eccentric isolation system,” Engineering Structures, vol. 28, no. 1, pp. 72–83, 2006. View at: Publisher Site  Google Scholar
 A. TenaColunga and J. L. EscamillaCruz, “Torsional amplifications in asymmetric baseisolated structures,” Engineering Structures, vol. 29, no. 2, pp. 237–247, 2007. View at: Publisher Site  Google Scholar
 C. P. Providakis, “Pushover analysis of baseisolated steel–concrete composite structures under nearfault excitations,” Soil Dynamics and Earthquake Engineering, vol. 28, no. 4, pp. 293–304, 2008. View at: Publisher Site  Google Scholar
 R. Tyler and W. Robinson, “Highstrain tests on leadrubber bearings for earthquake loadings,” Earthquake Engineering, vol. 17, no. 2, 1984. View at: Google Scholar
 H. C. Tsai, “Dynamic analysis of baseisolated shear beams bumping against stops,” Earthquake Engineering & Structural Dynamics, vol. 26, no. 5, pp. 515–528, 1997. View at: Publisher Site  Google Scholar
 J. Chen, C. Zhao, Q. Xu, and C. Yuan, “Seismic analysis and evaluation of the base isolation system in AP1000 NI under SSE loading,” Nuclear Engineering and Design, vol. 278, pp. 117–133, 2014. View at: Publisher Site  Google Scholar
 E. N. Dvorkin, D. Pantuso, and E. A. Repetto, “A formulation of the MITC4 shell element for finite strain elastoplastic analysis,” Computer Methods in Applied Mechanics and Engineering, vol. 125, no. 1–4, pp. 17–40, 1995. View at: Publisher Site  Google Scholar
 X. Lu, L. Xie, H. Guan, Y. Huang, and X. Lu, “A shear wall element for nonlinear seismic analysis of supertall buildings using OpenSees,” Finite Elements in Analysis and Design, vol. 98, pp. 14–25, 2015. View at: Publisher Site  Google Scholar
 L. Xie, X. Lu, X. Lu, Y. Huang, and L. Ye, “Multilayer shell element for shear walls in OpenSees,” Computing in Civil and Building Engineering, pp. 1190–1197, 2014. View at: Publisher Site  Google Scholar
 The Pacific Earthquake Engineering Research Center, PEER Database, University of California, Berkeley, CA, USA, 2013.
 M. S. Zadeh and M. S. Saiidi, PreTest Analytical Studies of NEESRSG 4Span Bridge Model Using OpenSees, University of Nevada, Reno, NV, USA, 2007, Tech. Rep. CCEER07.
 P. Jehel, P. Léger, and A. Ibrahimbegovic, “Initial versus tangent stiffnessbased Rayleigh damping in inelastic time history seismic analyses,” Earthquake Engineering & Structural Dynamics, vol. 43, no. 3, pp. 467–484, 2014. View at: Publisher Site  Google Scholar
 FEMA P695, Quantification of Building Seismic Performance Factors, Federal Emergency Management Agency, Washington, DC, USA, 2009.
 Engineers, A.S.o.C., Seismic Rehabilitation of Existing Buildings, vol. 41, ASCE Publications, Reston, VI, USA, 2007.
 D. Vamvatsikos and C. A. Cornell, “Incremental dynamic analysis,” Earthquake Engineering & Structural Dynamics, vol. 31, no. 3, pp. 491–514, 2002. View at: Publisher Site  Google Scholar
 FEMA 350, Recommended Seismic Design Criteria for New Steel MomentFrame Buildings, FEMA, Washington, DC, USA, 2000.
 E. DumovaJovanoska, “Fragility curves for reinforced concrete structures in Skopje (Macedonia) region,” Soil Dynamics and Earthquake Engineering, vol. 19, no. 6, pp. 455–466, 2000. View at: Publisher Site  Google Scholar
 G. Altay, G. Deodatis, G. Franco et al., “Benefitcost analysis for earthquake mitigation: evaluating measures for apartment houses in Turkey,” in Proceedings of the 2nd Annual IIASADPRI Meeting, Integrated Disaster Management, Laxenburg, Austria, July 2002. View at: Google Scholar
 Federal Emergency Management Agency, HAZUSMH MR4 Earthquake Model Technical Manual, Federal Emergency Management Agency, Washington, DC, USA, 2009.
 M. S. Kirçil and Z. Polat, “Fragility analysis of midrise R/C frame buildings,” Engineering Structures, vol. 28, no. 9, pp. 1335–1345, 2006. View at: Publisher Site  Google Scholar
 D. Vamvatsikos, F. Jalayer, and C. A. Cornell, “Application of incremental dynamic analysis to an RCstructure,” in Proceedings of the FIB Symposium on Concrete Structures in Seismic Regions, Athens, Greece, May 2003. View at: Google Scholar
 A. Gunduz, Probability, Statistics, Risk and Reliability in Engineering, Kure publications, Istanbul, Turkey, 1996.
Copyright
Copyright © 2018 Mohammad Parsaeimaram et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.