Wind-Induced Dynamic Response of Inter-Story Isolated Tall Buildings with Friction Pendulum Bearing Based on an Enhanced Simplified Model

: Isolation technology, especially for base isolation, is increasingly being applied in earthquake-prone areas. To satisfy some special demands (such as prevention from seawater erosion of an isolation layer, story-adding retrofit of existing buildings, avoidance of collision between base-isolated tall buildings, and so on), the isolation layer sometimes has to be set in the middle of a building to constitute inter-story isolated buildings. This special structural form inevitably encounters strong wind loads during service life, and its wind-resistant performance deserves to be investigated. This study conducts the wind-induced vibration analysis of inter-story isolated tall buildings with friction pendulum bearing (FPB). The nonlinear time domain analysis model and statistical linearization method to compute the wind-induced response of FPB inter-story isolated tall buildings are addressed based on an enhanced simplified model. Considering the independence of the upper and lower structures, two structural design schemes for inter-story isolated tall buildings are provided. Their dynamic characteristics are analyzed, and wind-induced responses are compared. Finally, the accuracy of the statistical linearization method is verified. This study provides an important theoretical basis for the structural design and wind resistance of inter-story isolated tall buildings.


Introduction
As one of the most effective earthquake-resistant technologies, base isolation has demonstrated its superior seismic performance after more than half a century of theoretical development and practical application [1][2][3].However, in engineering practice, the isolation layer sometimes has to be installed on the higher floor to prevent the seawater erosion of buildings located in coastal areas, or reduce excessive isolation layer displacement of those in high-density building areas [4].Additionally, the isolation system is adopted to serve as a "soft transfer" between the rigid chassis and flexible towers in order to avoid the stiffness discontinuity.Furthermore, the isolation technology is applied in retrofit strategy on existing buildings [5][6][7].These buildings are divided into two parts by the isolation layer to realize multi-functionality, such as commercial and residential functions [8].Thus, it has gained more attention and rapid development.Japan is at the forefront of the world in the field of building isolation and has built a few inter-story isolated tall buildings [8].The highest one, Roppongi Grand Tower, is up to 230 m [9].
Over the past two decades, extensive theoretical analyses have been conducted on the seismic resistance of inter-story isolated buildings [10][11][12][13][14][15].Recently, some research based on engineering cases has also been carried out, including shaking table tests on over-track multi-tower building [16,17] and story-adding retrofit of existing buildings by isolation bearings [18].The structural design of an inter-story isolated building is complex due to more dynamic parameters being involved.By specifying optimal control objectives, singleor multi-objective optimization methods based on genetic algorithms are used to determine these dynamic parameters [19][20][21][22].
The above study on the seismic performance of inter-story isolated buildings has obtained enough attention, and these efforts have laid the foundation for practical application of this type of structure.However, taller inter-story isolated buildings inevitably encounter strong wind loads during service life.The mainstream research of classical wind engineering is the wind-induced vibrations of tall buildings and long-span bridges [23][24][25].The research on the wind resistance of base-isolated tall buildings has become a hot topic in recent years.Siringoringo and Fujino [26] performed the real-time monitoring of a base-isolated building subject to typhoon.Feng and Chen [27,28] systematically studied the nonlinear alongwind and crosswind responses of base-isolated tall buildings with lead rubber bearings.Li et al. [29] provided the equivalent static wind load of the base-isolated tall building from the perspective of the load specification.Li et al. [30,31] investigated the uncoupled and three-dimensional coupled wind-induced response of a base-isolated tall building with friction pendulum bearing (FPB).Nevertheless, there is almost no report on the wind resistance of inter-story isolated tall buildings.
As the more generalized version of the base-isolated building, inter-story isolated tall buildings consist of three substructures: the lower structure (LS), isolation layer, and upper structure (US).The dynamic parameters of substructures have significant influences on the dynamic responses of those three parts.An efficient and accurate evaluation of the wind-induced response for this unique structural form is significant.In fact, the computation of wind-induced vibration of an isolated tall building is time-consuming, as the nonlinear mechanical characteristics of the isolation layer, multiple degrees of freedom (DOFs) involved in the dynamic model, and a large scale of sample of wind load time history are necessary to obtain the response statistic.The nonlinear time domain analysis of the structural system with a large number of DOFs could be inefficient for optimization problems in which the massive iterative computation is involved.Seeking an accurate and reasonable simplified model, and developing an effective statistical linearization method are important to improve the efficiency of nonlinear wind vibration analysis.
This paper extends the previous works [30,31] on the wind resistance of the FPB baseisolated tall building to the FPB inter-story isolated tall building.An analytical framework of the wind-induced vibration of an FPB inter-story isolated tall building by the nonlinear time domain analysis and statistical linearization method is first proposed.To improve computational efficiency, a reasonable simplified model with a few DOFs of the LS and US is used in the time domain analysis and statistical linearization method.Secondly, considering the independence of the LS and US, two structural design schemes for the inter-story isolated building are suggested, their modal characteristics are analyzed, and crosswind responses are compared.Then, the accuracy of the statistical linearization method is verified.Some important conclusions are highlighted at the end of the paper.

Theoretical Framework 2.1. Equations of Motion
As shown in Figure 1, an inter-story isolated building with friction pendulum bearing (FPB) consists of an M-story LS, isolation layer, and N-story US.The mass, damping, stiffness, and external load on the i-th story of the LS are m l,i , c l,i , k l,i , and P l,i (t) (i = 1, 2, . . . M), respectively.The counterparts on the j-th story of the US are denoted as m u,j , c u,j , k u,j , and P u,j (t) (j = 1, 2, . . . N), respectively.For the isolation layer, its mass, damping, and external load are m b , c b , and P b (t), respectively.The soil-structure interaction and dynamic interaction between adjacent structures are neglected in this study.

Physical Coordinates
The displacement of the isolation layer, with respect to the top story of the LS, is denoted as  .The displacement vectors of the LS and US are represented as   , ,  , , ⋯ ,  , and   , ,  , , ⋯ ,  , , respectively, where  is relative to the isolation slab.In physical coordinates, the equations of motion corresponding to three substructures of the inter-story isolated building are given by where  ,  , and  are the mass, damping, and stiffness matrices of the LS, respectively;  ,  , and  are the same metrices of the US;   ,  ,  ,  , …,  ,  and   ,  ,  ,  , …,  ,  are the wind load vectors along the stories of the LS and US, respectively; the force vectors  , and  , represent the inertial and external forces which act on the top story of the LS;  is the N-dimension unit column vector; the total weight of the US and isolation slab is   where       and  9.8 m/s ;  is the radius of the sliding plate. represents the friction force of the bearings, which is expressed as [32]     (3a) where M l , C l , and K l are the mass, damping, and stiffness matrices of the LS, respectively; M u , C u , and K u are the same metrices of the US; P l = [P l, 1 (t), P l, 2 (t), . . . ,P l, M (t)] T and P u = [P u, 1 (t), P u, 2 (t), . . . ,P u, N (t)] T are the wind load vectors along the stories of the LS and US, respectively; the force vectors F l,inertial and F l,load represent the inertial and external forces which act on the top story of the LS; r N is the N-dimension unit column vector; the total weight of the US and isolation slab is W = mg where m = r T N M u r N + m b and g = 9.8 m/s 2 ; R is the radius of the sliding plate.F f represents the friction force of the bearings, which is expressed as [32] .
where µ min and µ max are the friction coefficients when the velocity of isolation slab .u b = 0 and .u b = ∞, respectively.Factor a controls the variation of the friction coefficient from µ min to µ max ; and z b is the normalized hysteresis displacement.The parameters controlling the hysteresis loop shape need to meet the conditions γ + β = 1 (γ, β > 0) and η ≥ 1.A set of reasonable parameters γ = β = 0.5, η = 2 is used in this study.The low yield displacement is taken as Y = 0.13 mm.
In this study, the shear building model is used.The mass and stiffness matrices of the LS and US are represented as Because the LS can be regarded as a fixed-base building, its damping matrix C l adopts Rayleigh damping [33], i.e., C l = a l0 M l + a l1 K l , where coefficients a l0 and a l1 are determined by giving any two modal damping ratios.The US and isolation layer constitute a base-isolated building.A stiffness-proportional damping matrix is appropriate for the US, i.e., C u = a u1 K u [34].
The cross-power spectral density of fluctuating wind loads on the i-th and j-th stories can be derived from the base bending moment according to the Japanese code [28,29] where h i is the distance between the i-th story and ground; H and B are the height and width of the building; H 0 is the height of each story; α s is the wind profile coefficient; ρ is the air density; k y is the decay factor; U H is the mean wind speed at the building top; and S c M ( f ) is the power spectral density (PSD) of the coefficient of the base bending moment [35].Based on the PSD matrix of the story forces, the time history of fluctuating wind load is generated stochastically by the spectral representation method [36].When the load on stories is non-stationary wind, the corresponding non-stationary wind-induced responses are also obtained [37].

Modal Coordinates
The displacements of the LS and US can be expressed, respectively, by their first rand s-order modes (r ≤ M; s ≤ N): where the modal matrix composed of the first r-order mode shape of the LS and the first s-order mode shape of the US, respectively.They can be determined by the eigenvalue analysis using mass and stiffness matrices.For convenience, the last elements ϕ li,M of ϕ li (i = 1, 2, . .., r) and ϕ uj,N of ϕ uj (j = 1, 2, . .., s) are normalized as the unit.
Substituting Y l and Y u in Equation ( 6) into Equations (1a)-(1c), and then premultiplying Φ T l and Φ T u on both sides of Equation (1a) and Equation (1c), respectively, the equations of motion with dimensionless coefficients are derived as (r li + 1) .. .. .. where are the mass ratios of the LS and US, respec- are the modal participation mass factors; are the modal frequencies of the LS and US, respectively; are the modal damping ratios; are the frequency and linear damping ratio of the isolation layer; and the corresponding inputs for three parts are , and Furthermore, Equation (7) can also be written as the matrix form in which Λ l = diag(r l1 , . . ., r lr ), Λ u = diag(r u1 , . . ., r us ); where r r is the r-dimension unit column vector.The state-space form of the equation of motion is further given by This nonlinear equation of motion can be solved by the Runge-Kutta method to obtain the structural response history.
The structural wind-induced response is commonly controlled by low modes.A few low-order dynamic parameters of the LS and US and the mechanical parameters of the isolation layer determine the structural system.Theoretically, through machine learning, it is possible to establish a mapping relationship between structural response and those parameters.Subsequently, the wind-induced response of buildings can also be predicted accurately [38].

Enhanced Simplified Model with Reduced DOFs
The equations of motion of inter-story isolated buildings in physical coordinates involve the structural dynamics model with full DOFs.The nonlinear time domain analysis of a structural system with a large number of DOFs could be inefficient.By modal truncation of the LS and US, the equations of motion in the modal coordinates can greatly reduce the DOFs of a system.The simplified model LrUs is defined as the combination of LS with first r modes, US with first s modes, and an isolation system (the simplest model is L1U1 with 3 DOFs).There are M × N simplified models of the original structure when the different r and s are combined.A reasonable simplified model is capable of accurately estimating the response of the original structure with the number of DOFs as few as possible.The author [29] proposed a method to select DOFs of the LS and US based on the relationship of generalized stiffness between the original structure and substructure: where K eff, k is the k-th generalized stiffness of the original structure, k b,eff is the equivalent stiffness of the isolation layer, and K li = ϕ T li K l ϕ li and K uj = ϕ T uj K u ϕ uj are the i-th generalized stiffness of LS and the j-th one of US, respectively.q k = q l1, k . . .q li, k . . .q lM, k q bk q u1, k . . .q uj, k . . .q uN, k T is the mode vector corresponding to the undamped vibration of the original structure with full DOFs.
Using Equation (11), a reasonable reduced model is selected through the following procedure: (1) The number of modes considered for the original structure is labelled as Np.
(2) For those Np modes, their important generalized stiffness contributions from the LS, US, and isolation layer need to be up to 98%.The highest mode of the LS and US that participates the cumulative operation is denoted as r i and s i (i = 1, 2,. .., Np), respectively.
For the nonlinear wind-induced vibration of a FPB isolated tall building, the isolation layer is always in a non-sliding phase or sliding phase.The equivalent stiffness of the isolation layer k b,eff is taken as W[(1/R) + (µ min /Y)] or W/R to the corresponding two phases.The final simplified model is the result with the consideration of the above two equivalent stiffnesses.

Statistical Linearization Method
In nonlinear random vibration analysis, the statistics of structural responses are computed by response history analysis.A large number of sample computations are required.Another frequency-domain-based statistical linearization method to obtain the statistics of structural responses directly is introduced as follows.
The nonlinear term µ .
where coefficients b 1 and b 2 are derived according to the Gaussian distribution assumption [32]: where )dx; and σ .u b , σ z b , and ρ .u b z b are the standard deviations (STDs) of .u b and z b , as well as their correlation coefficient.
Baber and Wen [39] provided an analytic expression for b 3 , b 4 , and b 5 : in which where Accordingly, the corresponding linearized equation for an FPB inter-story isolated building is given by where The statistics of the steady-state response can be determined by where S PP (ω) is the PSD matrix of the story wind load vector P = [P T l P b P T u ] T .
In the statistical linearization method, Equation ( 16

Building Description
For comparison, the fixed-base, base-isolated and inter-story isolated tall buildings, shown in Figure 2, are involved in this study.They have identical building dimensions, with a total height of H = 202 m and square section of B = D = 0.2 H.The story heights of those three building are 4 m.The 2 m rigid foundation and isolation layer are set for fixed-base or isolated buildings.The first story of the fixed-base and base-isolated buildings is supported by rigid foundation and isolation bearings, respectively.Thus, all the buildings have 51 stories.

Building Description
For comparison, the fixed-base, base-isolated and inter-story isolated tall buildings, shown in Figure 2, are involved in this study.They have identical building dimensions, with a total height of  202 m and square section of   0.2 .The story heights of those three building are 4 m.The 2 m rigid foundation and isolation layer are set for fixed-base or isolated buildings.The first story of the fixed-base and base-isolated buildings is supported by rigid foundation and isolation bearings, respectively.Thus, all the buildings have 51 stories.
For the fixed-base building, the mass of each story is  1.2288 10 kg according to the building density  192 kg/m .The first frequency is  0.20 Hz .The story stiffness is determined by use of the linear first mode shape, i.e.,   /.According to eigenvalue analysis, the first three order frequencies of the structure are 0.20, 0.49, and 0.77 Hz, respectively.The building adopts Rayleigh damping with the first and third modal damping ratios   1% to determine the damping matrix.The eigenvalue analysis indicates that the first three modal damping ratios of a fixed-base building are 1.00%, 0.83%, and 1.00%, respectively.A base-isolated building is regarded as a flexible fixed-base building which replaces the rigid foundation by the isolation layer.Ryan and Polanco [34] indicated that a stiffnessproportional damping matrix is appropriate for the US of base-isolated buildings.The first modal damping ratio  1% is calibrated to define this damping matrix.For the isolation system, the mass of the isolation slab is equal to that of one story of the US, i.e.,  For the fixed-base building, the mass of each story is m i = 1.2288 × 10 6 kg according to the building density ρ s = 192 kg/m 3 .The first frequency is f 01 = 0.20 Hz.The story stiffness is determined by use of the linear first mode shape, i.e., ϕ 01 (z) = z/H.According to eigenvalue analysis, the first three order frequencies of the structure are 0.20, 0.49, and 0.77 Hz, respectively.The building adopts Rayleigh damping with the first and third modal damping ratios ζ 01 = ζ 03 = 1% to determine the damping matrix.The eigenvalue analysis indicates that the first three modal damping ratios of a fixed-base building are 1.00%, 0.83%, and 1.00%, respectively.
A base-isolated building is regarded as a flexible fixed-base building which replaces the rigid foundation by the isolation layer.Ryan and Polanco [34] indicated that a stiffnessproportional damping matrix is appropriate for the US of base-isolated buildings.The first modal damping ratio ζ 01 = 1% is calibrated to define this damping matrix.For the isolation system, the mass of the isolation slab is equal to that of one story of the US, i.e., m b = m i .The radius of the circular sliding surface of the FPB is R = 2 m.The parameters related to the friction coefficient are taken as µ min = 0.01, µ min /µ min = 2 and a = 23.6 s/m.Additionally, it is assumed that the dissipative energy of the isolation system is only from the friction mechanism, and the linear damping ratio is ζ b = 0.
An inter-story isolated building consists of the isolation layer, LS, and US.The parameters of the isolation system are the same as those in a base-isolated building.In fact, there is no uniform regulation for the preliminary structural design for the LS and US.Two design schemes to define the LS and US are introduced as follows: Scheme A: Similar to the base-isolated building, the isolation layer is placed in any middle story of the fixed-base building, instead of the base.The story mass and stiffness of the LS and US are determined by the fixed-base building.
Scheme B: The LS and US are independent of each other and are designed separately.Specifically, the masses of the LS and US are calculated based on the building density, ρ l = ρ u = 192 kg/m 3 .The first modal frequencies of the LS and US are f l1 =10/M and f u1 = 10/N Hz, respectively.By use of the linear first mode shape, the story stiffness of each part is obtained.
The damping matrices of the LS and US are determined by their respective modal damping ratios in two design schemes.For the LS, the damping matrix C l is calculated according to In scheme A, a fixed-base tall building with a specified height is assumed, the to-bebuilt inter-story isolated tall building is designed by installing an isolation layer in the middle of this fixed-base building.Similar to base-isolated buildings, this structural design of an inter-story isolated tall building only needs to consider the dynamic characteristics of the fixed-base building and the mechanical parameters of the isolation layer.The design method is simple, and the integrity of the building is great.Scheme B is more flexible, especially when the characteristics of the LS and US are significantly different.It is suitable for isolated buildings on large bottom podiums or over tracks and adding-story-adding retrofit on existing buildings by isolation bearings.

Scheme A
The dynamic model of an inter-story isolated tall building has M + 1 + N = 51 degrees of freedom, denoted as 51DOFs.According to Section 2.2, a method to select the simplified model of the original building with a 20-story LS (M = 20) is illustrated here.Figure 3 presents the contributions of generalized stiffness from the LS, isolation layer, and US to the first three generalized stiffnesses of the original building when the isolation layer is in the non-sliding and sliding phases.It can be seen that the generalized stiffness of the original building is contributed by only a few low modes of the LS and US.For example, the first modes of the LS and US contribute 91.24% to the first generalized stiffness of the original building, as the isolation layer in the non-sliding phase.This data is up to 92.86% in the sliding phase.If the first four DOFs of the LS and US are taken, the first three modes of the original building have an estimation with an accuracy up to 98% in both the non-sliding and sliding phases.Therefore, the reasonable simplified building model is L4U4.
The location of the isolation layer is labeled as M, i.e., the number of stories of the LS.The base-isolated building is represented by M = 0. Table 1 gives selective simplified models of inter-story isolated buildings with different isolation locations.
Figures 4 and 5 show the variation in the first three frequencies and damping ratios for the simplified model LrUs with the different isolation location when the isolation layer is in the non-sliding or sliding phase.The results from 51DOFs and L1U1 models are also provided for comparison.In addition, Figures 6 and 7    Figures 4 and 5 show the variation in the first three frequencies and damping rat for the simplified model LrUs with the different isolation location when the isolation la is in the non-sliding or sliding phase.The results from 51DOFs and L1U1 models are a provided for comparison.In addition, Figures 6 and 7 depict the first three mode sha of the inter-story isolated building with M = 20.It is observed that scheme A does not alter the story stiffness of the LS and US with the elevation of the isolation layer.Therefore, the first modal frequency decreases slightly, while the second and third ones remain almost unchanged when the isolation layer is in the non-sliding phase.The modal frequencies in the sliding phase are always lower than those in the non-sliding phase.Their variation becomes more complex with the rise in the isolation layer.It is observed that scheme A does not alter the story stiffness of the LS and US with the elevation of the isolation layer.Therefore, the first modal frequency decreases slightly, while the second and third ones remain almost unchanged when the isolation layer is in the non-sliding phase.The modal frequencies in the sliding phase are always lower than those in the non-sliding phase.Their variation becomes more complex with the rise in the isolation layer.It is observed that scheme A does not alter the story stiffness of the LS and US with the elevation of the isolation layer.Therefore, the first modal frequency decreases slightly, while the second and third ones remain almost unchanged when the isolation layer is in the non-sliding phase.The modal frequencies in the sliding phase are always lower than those in the non-sliding phase.Their variation becomes more complex with the rise in the isolation layer.With the elevation in the isolation layer, the first modal damping ratio initially decreases and then increases, while the second and third modal damping ratios reduce on the whole.In fact, the modal damping ratios of the original building (51DOFs) are the combination of those of the LS, isolation layer, and US [15].
In the non-sliding phase, the structural mode shapes are similar to those of the fixedbase building, while in the sliding phase, the linear stiffness (/) provided by the bearings is always smaller than the adjacent story stiffness of the LS and US, resulting in a clear drift at the isolation layer.
It is evident that the simplified model, LrUs, can accurately estimate the first three modal characteristics of the original building.Strictly speaking, even higher modes (such as the fourth or fifth mode) also have a great prediction.Because the wind-induced response of an isolated tall building is primarily controlled by the first three modes, the simplified model LrUs can entirely represent the model 51DOFs to compute wind-induced responses.The model L1U1 overestimates the first frequency and damping ratio and provides a different first mode shape.The estimations of the second and third modal characteristics are usually inaccurate.

Scheme B
When structural design adopts scheme B, the selection for the DOF of an appropriate simplified model is presented in Table 2. Figures 8 and 9 illustrate the first three modal frequencies and damping ratios.Because the first frequencies of the LS and US are inversely proportional to their respective numbers of stories, it is the softest for the whole structure when the isolation layer is positioned near the middle of the building.Therefore, the frequency decreases first and then increases with the isolation location M. The minimum first frequency is observed at around M = 25.Compared with the first frequency, the change in the second frequency with M has the opposite trend, and the third frequency has a more complex variation.The damping matrices for the LS and US are set independently.Therefore, the variations in modal damping ratios are similar to those in scheme A.  It is observed that scheme A does not alter the story stiffness of the LS and US with the elevation of the isolation layer.Therefore, the first modal frequency decreases slightly, while the second and third ones remain almost unchanged when the isolation layer is in the non-sliding phase.The modal frequencies in the sliding phase are always lower than those in the non-sliding phase.Their variation becomes more complex with the rise in the isolation layer.
With the elevation in the isolation layer, the first modal damping ratio initially decreases and then increases, while the second and third modal damping ratios reduce on the whole.In fact, the modal damping ratios of the original building (51DOFs) are the combination of those of the LS, isolation layer, and US [15].
In the non-sliding phase, the structural mode shapes are similar to those of the fixedbase building, while in the sliding phase, the linear stiffness (W/R) provided by the bearings is always smaller than the adjacent story stiffness of the LS and US, resulting in a clear drift at the isolation layer.
It is evident that the simplified model, LrUs, can accurately estimate the first three modal characteristics of the original building.Strictly speaking, even higher modes (such as the fourth or fifth mode) also have a great prediction.Because the wind-induced response of an isolated tall building is primarily controlled by the first three modes, the simplified model LrUs can entirely represent the model 51DOFs to compute wind-induced responses.The model L1U1 overestimates the first frequency and damping ratio and provides a different first mode shape.The estimations of the second and third modal characteristics are usually inaccurate.

Scheme B
When structural design adopts scheme B, the selection for the DOF of an appropriate simplified model is presented in Table 2. Figures 8 and 9 illustrate the first three modal frequencies and damping ratios.Because the first frequencies of the LS and US are inversely proportional to their respective numbers of stories, it is the softest for the whole structure when the isolation layer is positioned near the middle of the building.Therefore, the frequency decreases first and then increases with the isolation location M. The minimum first frequency is observed at around M = 25.Compared with the first frequency, the change in the second frequency with M has the opposite trend, and the third frequency has a more complex variation.The damping matrices for the LS and US are set independently.Therefore, the variations in modal damping ratios are similar to those in scheme A.   In scheme B, the story stiffness of the LS and US varies with the elevation of the isolation layer, and the LS story stiffness near the top is typically smaller than the US story stiffness adjacent to the bottom.Therefore, there is a significant mode displacement corresponding to stories close to the top of the LS (shows in Figure 10a).Nevertheless, the isolation layer stiffness in the sliding phase, /, is always smaller than its adjacent story stiffness of the LS and US, which leads to an obvious mode displacement at the isolation layer (shows in Figure 11).
Compared to scheme A, the model L1U1 causes a larger error for the prediction on the structural modal characteristics.Fortunately, the LrUs model still provides a satisfactory estimation.In scheme B, the story stiffness of the LS and US varies with the elevation of the isolation layer, and the LS story stiffness near the top is typically smaller than the US story stiffness adjacent to the bottom.Therefore, there is a significant mode displacement corresponding to stories close to the top of the LS (shows in Figure 10a).Nevertheless, the isolation layer stiffness in the sliding phase, /, is always smaller than its adjacent story stiffness of the LS and US, which leads to an obvious mode displacement at the isolation layer (shows in Figure 11).
Compared to scheme A, the model L1U1 causes a larger error for the prediction on the structural modal characteristics.Fortunately, the LrUs model still provides a satisfactory estimation.In scheme B, the story stiffness of the LS and US varies with the elevation of the isolation layer, and the LS story stiffness near the top is typically smaller than the US story stiffness adjacent to the bottom.Therefore, there is a significant mode displacement corresponding to stories close to the top of the LS (shows in Figure 10a).Nevertheless, the isolation layer stiffness in the sliding phase, W/R, is always smaller than its adjacent story stiffness of the LS and US, which leads to an obvious mode displacement at the isolation layer (shows in Figure 11).

Wind-Induced Response
This study focuses on the buffeting response of an inter-story isolated building to a fluctuating wind load.The wind load on the stories in Equation ( 5) is based on a quasistatic assumption.The dynamic wind-structure interaction is ignored, and the vortex-induced vibration is not considered here.
Because of the non-zero mean component, the along-wind load of the building is obviously different from the crosswind load.The mean wind load affects the displacement at the isolation layer before arriving to the steady state.However, this influence is significant only in the case of small sliding of the bearing under the low wind speed.When a distinct sliding occurs at the isolation layer, its displacement quickly reaches a steady state, and the along-wind response can be separately computed for the static response under the mean load and the dynamic response under the fluctuating load [30].There are consistent solution procedures and similar structural response characteristics for alongwind and crosswind fluctuating responses.In the following discussion, the crosswind vibration of an inter-story isolated building will be focused on to analyze the wind-induced response characteristic and validate the effectiveness of the statistical linearization method.
The building locates in II terrain category with wind profile coefficient  0.

Wind-Induced Response
This study focuses on the buffeting response of an inter-story isolated building to a fluctuating wind load.The wind load on the stories in Equation ( 5) is based on a quasistatic assumption.The dynamic wind-structure interaction is ignored, and the vortex-induced vibration is not considered here.
Because of the non-zero mean component, the along-wind load of the building is obviously different from the crosswind load.The mean wind load affects the displacement at the isolation layer before arriving to the steady state.However, this influence is significant only in the case of small sliding of the bearing under the low wind speed.When a distinct sliding occurs at the isolation layer, its displacement quickly reaches a steady state, and the along-wind response can be separately computed for the static response under the mean load and the dynamic response under the fluctuating load [30].There are consistent solution procedures and similar structural response characteristics for alongwind and crosswind fluctuating responses.In the following discussion, the crosswind vibration of an inter-story isolated building will be focused on to analyze the wind-induced response characteristic and validate the effectiveness of the statistical linearization method.
The building locates in II terrain category with wind profile coefficient  0. Compared to scheme A, the model L1U1 causes a larger error for the prediction on the structural modal characteristics.Fortunately, the LrUs model still provides a satisfactory estimation.

Wind-Induced Response
This study focuses on the buffeting response of an inter-story isolated building to a fluctuating wind load.The wind load on the stories in Equation ( 5) is based on a quasi-static assumption.The dynamic wind-structure interaction is ignored, and the vortex-induced vibration is not considered here.
Because of the non-zero mean component, the along-wind load of the building is obviously different from the crosswind load.The mean wind load affects the displacement at the isolation layer before arriving to the steady state.However, this influence is significant only in the case of small sliding of the bearing under the low wind speed.When a distinct sliding occurs at the isolation layer, its displacement quickly reaches a steady state, and the along-wind response can be separately computed for the static response under the mean load and the dynamic response under the fluctuating load [30].There are consistent solution procedures and similar structural response characteristics for along-wind and crosswind fluctuating responses.In the following discussion, the crosswind vibration of an inter-story isolated building will be focused on to analyze the wind-induced response characteristic and validate the effectiveness of the statistical linearization method.
The building locates in II terrain category with wind profile coefficient α s = 0.15.A decay factor k y = 5 is used in crosswind PSDs of the story wind loads.The PSD of the crosswind base bending moment coefficient is determined from reference [35] with its STD σ CM = 0.16, as shown in Figure 12.Using the spectral representation method, 100 history samples with a duration of 15 min for the story wind loads are simulated.The PSD of the base moment coefficient obtained through stochastic simulation is consistent with that from the Japanese code.This study focuses on the statistics of structural responses in the steady state, and the first 5 min of the response history are discarded to eliminate the transient effect in the following analysis.
ings 2024, 14, x FOR PEER REVIEW the Japanese code.This study focuses on the statistics of structura state, and the first 5 min of the response history are discarded to eli in the following analysis.It is observed that the inter-story drift angle of the fixed-base buil < 1/550 (the limitation of the inter-story drift angle of the reinforc ture [40]).The drift at the isolation layer increases with its elevati of the LS and US remain almost at the same level.The story drifts isolation layer have a noticeable change, which leads to an obvio forces of corresponding stories.Apart from the isolation layer b top, story accelerations change a little with the isolation location.ergy dissipation of bearings increases the damping of the structu fore, the responses of isolated buildings reduce significantly com fixed-base building.For instance, the acceleration at the buildin  For the isolation location M = 10, 20, 30, and 40, the STDs of responses, including displacement, drift, acceleration, and shear force, of inter-story isolated buildings designed through scheme A to crosswind load with U H = 50 m/s are portrayed in Figure 13.The corresponding results for a fixed-base building and the base-isolated building (M = 0) are also plotted, which are comparable with those in reference [28].The STDs of the top displacements of the fixed-base building in this study and in reference [28] are 0.26 and 0.23 m, respectively; and the STDs of the base shear forces of the base-isolated buildings shown in Figure 13d and in that reference are 7.84 × 10 6 and 8.12 × 10 6 kN, respectively.It is observed that the inter-story drift angle of the fixed-base building is 0.263/202 = 0.0013 < 1/550 (the limitation of the inter-story drift angle of the reinforced concrete frame structure [40]).The drift at the isolation layer increases with its elevation, while the story drifts of the LS and US remain almost at the same level.The story drifts of the LS adjacent to the isolation layer have a noticeable change, which leads to an obvious variation in the shear forces of corresponding stories.Apart from the isolation layer being positioned near the top, story accelerations change a little with the isolation location.The friction-induced energy dissipation of bearings increases the damping of the structural system [30,31].Therefore, the responses of isolated buildings reduce significantly compared with those of the fixed-base building.For instance, the acceleration at the building top decreases by 44% and the shear force at the building base declines by 39% when the isolation location is M = 30.By use of the LrUs and L1U1 models, Table 3 summarizes the STDs of some important responses of isolated buildings with different isolation locations.These responses include top displacement of the LS, drift at the isolation layer, top displacement of the US relative to isolated layer, top accelerations of the LS and US, as well as base shear forces of the LS and US, which are denoted as Dl top, Di, Du top, Al top, Au top, Vl base, and Vu base, respectively.The ratios of structural response STDs obtained from the L1U1 and LrUs models are shown in Figure 14.By use of the LrUs and L1U1 models, Table 3 summarizes the STDs of some important responses of isolated buildings with different isolation locations.These responses include top displacement of the LS, drift at the isolation layer, top displacement of the US relative to isolated layer, top accelerations of the LS and US, as well as base shear forces of the LS and US, which are denoted as D l top , D i , D u top , A l top , A u top , V l base , and V u base , respectively.The ratios of structural response STDs obtained from the L1U1 and LrUs models are shown in Figure 14.It can be found that the structural displacement is primarily dominated by the first mode, while the acceleration is controlled by the first three modes.As the isolation layer is raised, Al top is gradually dominated by the first mode, while the first mode is always the main contribution to Au top.Due to the inaccurate estimation on modal characteristics, the structural responses provided by L1U1 show errors more or less.Specifically, Vl base is overestimated and other responses are underestimated.Dl top is more prone to underestimation than Du top.It is observed that the stories near the top of the LS exhibit a significant drift with the relatively low story stiffness, while the drift of the US can be effectively mitigated.With the rise in isolation location, the drift at the building top decreases greatly, and the largest drift at the isolation layer continues to grow.Because the structural response is dominated by the first mode, and the first modal frequency and damping ratio of the inter-story isolated building are always lower than those of the base-isolated building (as shown in Figures 8a and 9a), the displacement, acceleration, and shear forces of the former are consistently larger than those of the latter.The maximal wind-induced response appears at M = 30, at which the acceleration and shear force approach the results of the fixed-base building.Similarly, due to the additional damping from the isolation layer, the story acceleration and shear force of the isolated building are, on the whole, less than those of the fixedbase building.It can be found that the structural displacement is primarily dominated by the first mode, while the acceleration is controlled by the first three modes.As the isolation layer is raised, A l top is gradually dominated by the first mode, while the first mode is always the main contribution to A u top .Due to the inaccurate estimation on modal characteristics, the structural responses provided by L1U1 show errors more or less.Specifically, V l base is overestimated and other responses are underestimated.D l top is more prone to underestimation than D u top .

Scheme B
With scheme B, the STDs of displacement, drift, acceleration, and shear force along the stories of the inter-story isolated building are depicted in Figure 16.Table 4 summarizes the STDs of important responses computed by simplified models L1U1 and LrUs.The ratios of response STDs by those two models are illustrated in Figure 17.The normalized PSDs of D l top , D u top , A l top , and A u top are displayed in Figure 18.It should be pointed that for a known site and specified height of the inter-story isolated tall building, the location of the isolation layer directly affects the dynamic characteristics of the LS and US, thereby affecting the dynamic response of the building.The optimal location of the isolation layer can be further determined through genetic algorithm-based single or multi-objective optimization methods [21,22].The optimization objectives (related to the building's intended use and safety) such as minimizing the displacement of the isolation layer and the acceleration of the LS and US will directly affect the location of the isolation layer.It should be pointed that for a known site and specified height of the inter-story isolated tall building, the location of the isolation layer directly affects the dynamic characteristics of the LS and US, thereby affecting the dynamic response of the building.The optimal location of the isolation layer can be further determined through genetic algorithm-based single or multi-objective optimization methods [21,22].The optimization objectives (related to the building's intended use and safety) such as minimizing the displacement of the isolation layer and the acceleration of the LS and US will directly affect the location of the isolation layer.

Accuracy of the Statistical Linearization Method
Taking the inter-story isolated building (M = 20) designed by scheme A as an example, Figure 19  It is observed that the stories near the top of the LS exhibit a significant drift with the relatively low story stiffness, while the drift of the US can be effectively mitigated.With the rise in isolation location, the drift at the building top decreases greatly, and the largest drift at the isolation layer continues to grow.Because the structural response is dominated by the first mode, and the first modal frequency and damping ratio of the inter-story isolated building are always lower than those of the base-isolated building (as shown in Figures 8a and 9a), the displacement, acceleration, and shear forces of the former are consistently larger than those of the latter.The maximal wind-induced response appears at M = 30, at which the acceleration and shear force approach the results of the fixed-base building.Similarly, due to the additional damping from the isolation layer, the story acceleration and shear force of the isolated building are, on the whole, less than those of the fixed-base building.
In scheme B, the estimation of structural dynamic properties by L1U1 is worse than that in scheme A, resulting in less accurate predictions of structural responses.For example, D l top is underestimated by 50% when M = 30.With the comparison between Figures 15  and 18, scheme B provides a more flexible building structure, and its structural response is predominantly contributed by the first and second modes.
It should be pointed that for a known site and specified height of the inter-story isolated tall building, the location of the isolation layer directly affects the dynamic characteristics of the LS and US, thereby affecting the dynamic response of the building.The optimal location of the isolation layer can be further determined through genetic algorithm-based single or multi-objective optimization methods [21,22].The optimization objectives (related to the building's intended use and safety) such as minimizing the displacement of the isolation layer and the acceleration of the LS and US will directly affect the location of the isolation layer.

Accuracy of the Statistical Linearization Method
Taking the inter-story isolated building (M = 20) designed by scheme A as an example, Figure 19 compares the response STDs computed by response history analysis (RHA) and statistical linearization method (SLM) in the wind speed range U H = 20 ∼ 60 m/s.It can be observed that the SLM slightly overestimates displacement at the isolation layer at a high wind speed, while it underestimates the accelerations at the top of the LS, somewhat.The other structural responses have a great prediction.Through the eigenvalue analysis of state-space equation, i.e., Equation ( 16), the equivalent dynamic properties of the inter-story isolated building with nonlinearity can be obtained.Figure 20 provides the variation in the first three frequencies and damping ratios of the equivalent linear system with the wind speed.It is clear that the first and third frequencies of the building decrease slightly with the growth of   , while the second frequency is almost unchanged.This variation is consistent with the result in Figure 4a where frequencies change with the bearing state from the non-sliding phase to the sliding phase.The friction-induced energy dissipation brings about a significant increase in the first and third modal damping ratios.It is the additional damping at the isolation layer, the windinduced responses of inter-story isolated tall buildings are reduced, on the whole.Through the eigenvalue analysis of state-space equation, i.e., Equation ( 16), the equivalent dynamic properties of the inter-story isolated building with nonlinearity can be obtained.Figure 20 provides the variation in the first three frequencies and damping ratios of the equivalent linear system with the wind speed.It is clear that the first and third frequencies of the building decrease slightly with the growth of U H , while the second frequency is almost unchanged.This variation is consistent with the result in Figure 4a where frequencies change with the bearing state from the non-sliding phase to the sliding phase.The friction-induced energy dissipation brings about a significant increase in the first and third modal damping ratios.It is the additional damping at the isolation layer, the wind-induced responses of inter-story isolated tall buildings are reduced, on the whole.Through the eigenvalue analysis of state-space equation, i.e., Equation ( 16), the equivalent dynamic properties of the inter-story isolated building with nonlinearity can be obtained.Figure 20 provides the variation in the first three frequencies and damping ratios of the equivalent linear system with the wind speed.It is clear that the first and third frequencies of the building decrease slightly with the growth of  , while the second frequency is almost unchanged.This variation is consistent with the result in Figure 4a where frequencies change with the bearing state from the non-sliding phase to the sliding phase.The friction-induced energy dissipation brings about a significant increase in the first and third modal damping ratios.It is the additional damping at the isolation layer, the windinduced responses of inter-story isolated tall buildings are reduced, on the whole.It should be noted that the SLM provides an equivalent linear system to represent a nonlinear system by minimizing the mean square error between these two systems.The accuracy of the SLM depends on the degree of nonlinearity of the responses, i.e., the effectiveness of the SLM is related to the structural nonlinearity and the input load.It should be noted that the SLM provides an equivalent linear system to represent a nonlinear system by minimizing the mean square error between these two systems.The accuracy of the SLM depends on the degree of nonlinearity of the responses, i.e., the effectiveness of the SLM is related to the structural nonlinearity and the input load.

Conclusions
This study conducted the wind-induced vibration analysis of an inter-story isolated tall building with FPBs.The nonlinear time domain analysis framework and statistical linearization method to obtain the wind-induced response of the building are built by a simplified model with reduced DOFs.Considering the independence of the LS and US, two structural design schemes for FPB inter-story isolated tall buildings are provided.Their dynamic properties are analyzed, and the crosswind responses are compared.Finally, the accuracy of the statistical linearization method is confirmed.Some important conclusions are listed as follows: 1.The main factors affecting structural responses include low-order modal parameters of the LS and US (mass ratios, frequency ratios, damping ratios), and the friction coefficient and sliding plate radius of the FPB, as well as wind loads.
2. Although the wind-induced vibration of inter-story isolated tall buildings is controlled by low modes, the simplified model L1U1 consisting of the first mode of the LS and US cannot accurately model the first three dynamic characteristics of the original building, resulting in a significant over-and under-estimation of the structural responses.For instance, the top displacement of the LS is underestimated by nearly 50% in the case of structural design scheme B. The proposed simplified model LrUs can predict the first three modes precisely, and the structural response also has an accurate estimation.
3. The story stiffness of the LS and US in design scheme A is determined by the corresponding fixed-base building.The existence of the isolation layer affects the response of stories near the isolation layer.The displacement at isolation layer increases with its elevation, other responses are not significantly affected by the isolation location.The top acceleration and base shear force of the inter-story isolated building with M = 30 reduce by 44% and 30%, respectively, which attributes to the additional damping from friction-induced energy dissipation.
4. The story stiffness of the LS and US in design scheme B is defined independently.Its dynamic characteristic and wind-induced response vary greatly with the isolation location.When the isolation layer is near to middle position (M = 30), the smallest first modal frequency leads to the highest structural response and the worst isolation performance.Overall, scheme B makes the building structure more flexible than scheme A, and gives the greater wind-induced vibration.
5. Except for a slight underestimation of the top acceleration of the LS, and an overestimation of the displacement of the isolation layer at the high wind speed, the statistical linearization is accurate to predict other structural responses.
This study investigates the wind-induced response of FPB inter-story isolated buildings from a theoretical perspective.In future research, it is necessary to conduct a corresponding experimental study to confirm the accuracy and reliability of the proposed methods and design schemes.In addition, the response of a hybrid isolation system consisting of rubber bearings and friction pendulum bearings to three-dimensional synchronous wind loads also deserves investigation.

Figure 1 .
Figure 1.Shear model of inter-story isolated tall buildings.

Figure 1 .
Figure 1.Shear model of inter-story isolated tall buildings.2.1.1.Physical Coordinates The displacement of the isolation layer, with respect to the top story of the LS, is denoted as u b .The displacement vectors of the LS and US are represented as Y l = [y l,1 , y l,2 , . . ., y l,M ] T and Y u = [y u,1 , y u,2 , . . ., y u,N ] T , respectively, where Y u is relative to the isolation slab.In physical coordinates, the equations of motion corresponding to three substructures of the inter-story isolated building are given by

. u b z b in Equation ( 8 )
and the dimensionless hysteresis velocity term .z b in Equation (3c) are linearized as µ .u b z b = b 1 .u b + b 2 z b ) always contains the unknown parameters b 1 , b 2 , b 4 , and b 5 .These parameters depend on the STD of responses .u b and z b , as well as their correlation coefficient.Therefore, an iterative computation is needed.These initial values of parameters can be taken as b 1 = 0, b 2 = µ min , b 4 = 1/Y, and b 5 = 0.
depict the first three mode shapes of the inter-story isolated building with M = 20.dings 2024, 14, x FOR PEER REVIEW 10 o (a) Isolation layer in the non-sliding phase.(b) Isolation layer in the sliding phase.

Figure 3 .
Figure 3. Contribution to the generalized stiffness of an inter-story isolated tall building with M 20.

Figure 3 .
Figure 3. Contribution to the generalized stiffness of an inter-story isolated tall building with M = 20.

Figure 4 .Figure 5 .
Figure 4. First three modal frequencies of building designed by scheme A.

Figure 6 .Figure 4 .
Figure 6.First three mode shapes of a building in the non-sliding phase M = 20.

Figure 4 .Figure 5 .
Figure 4. First three modal frequencies of building designed by scheme A.

Figure 6 .Figure 5 .
Figure 6.First three mode shapes of a building in the non-sliding phase M = 20.

Figure 4 .Figure 5 .
Figure 4. First three modal frequencies of building designed by scheme A.

Figure 6 .Figure 6 .Figure 7 .
Figure 6.First three mode shapes of a building in the non-sliding phase M = 20.

Figure 7 .
Figure 7. First three mode shapes of a building in the sliding phase M = 20.

Figure 8 .Figure 9 .
Figure 8.First three modal frequencies of building designed by scheme B.

Figure 8 .
Figure 8.First three modal frequencies of building designed by scheme B.

Figure 8 .Figure 9 .
Figure 8.First three modal frequencies of building designed by scheme B.

Figure 9 .
Figure 9. First three modal damping ratios of building designed by scheme B.

Figure 10 .Figure 11 .
Figure 10.First three mode shapes of a building in the non-sliding phase (M = 20).

Figure 11 .
Figure 11.First three mode shapes of a building in the sliding phase (M = 20).

Figure 12 .
Figure 12.PSD of the crosswind base bending moment coefficient.

4. 1 .
Dynamic Response of Buildings 4.1.1.Scheme A For the isolation location M = 10, 20, 30, and 40, the STDs displacement, drift, acceleration, and shear force, of inter-stor signed through scheme A to crosswind load with  50 m/s 13.The corresponding results for a fixed-base building and the b = 0) are also plotted, which are comparable with those in referen top displacements of the fixed-base building in this study and and 0.23 m, respectively; and the STDs of the base shear forces o ings shown in Figure 13d and in that reference are 7.84 × 10 6 and 8

2 Figure 12 .
Figure 12.PSD of the crosswind base bending moment coefficient.

Figure 13 .
Figure 13.Height-wise variation in the response STD of inter-story isolated tall buildings designed by scheme A.

Figure 13 .
Figure 13.Height-wise variation in the response STD of inter-story isolated tall buildings designed by scheme A.
Figure 15 provides the PSD of D l top , D u top , A l top , and A u top normalized by their variances as M = 0, 10, 20, 30, and 40.

Figure 14 .
Figure 14.Ratios of the response STDs according to model L1U1 to those according to model LrUs (scheme A).

Figure 14 .
Figure 14.Ratios of the response STDs according to model L1U1 to those according to model LrUs (scheme A).
Buildings 2024, 14, x FOR PEER REVIEW 18 of 25 (a) Top displacement of LS (b) Top displacement of US (c) Top acceleration of the LS.(d) Top acceleration of the US.

Figure 15 .
Figure 15.PSDs of the building top displacement and acceleration in scheme A. 4.1.2.Scheme B With scheme B, the STDs of displacement, drift, acceleration, and shear force along the stories of the inter-story isolated building are depicted in Figure 16.Table 4 summarizes the STDs of important responses computed by simplified models L1U1 and LrUs.The ratios of response STDs by those two models are illustrated in Figure 17.The normalized PSDs of Dl top, Du top, Al top, and Au top are displayed in Figure 18.It is observed that the stories near the top of the LS exhibit a significant drift with the relatively low story stiffness, while the drift of the US can be effectively mitigated.With the rise in isolation location, the drift at the building top decreases greatly, and the largest drift at the isolation layer continues to grow.Because the structural response is dominated by the first mode, and the first modal frequency and damping ratio of the inter-story isolated building are always lower than those of the base-isolated building (as shown in Figures 8a and 9a), the displacement, acceleration, and shear forces of the former are consistently larger than those of the latter.The maximal wind-induced response appears at M = 30, at which the acceleration and shear force approach the results of the fixed-base building.Similarly, due to the additional damping from the isolation layer, the story acceleration and shear force of the isolated building are, on the whole, less than those of the fixedbase building.

Figure 15 .
Figure 15.PSDs of the building top displacement and acceleration in scheme A.

Buildings 2024 ,Figure 16 .
Figure 16.Height-wise variation of the response STD of inter-story isolated tall buildings designed by scheme B.

Figure 16 .
Figure 16.Height-wise variation of the response STD of inter-story isolated tall buildings designed by scheme B.

Figure 17 .
Figure 17.Ratios of response STDs with the model L1U1 compared to those with the model LrUs (scheme B).

Figure 18 .
Figure 18.PSDs of the building top displacement and acceleration in scheme B.

Figure 17 . 25 Figure 17 .
Figure 17.Ratios of response STDs with the model L1U1 compared to those with the model LrUs (scheme B).

Figure 18 .
Figure 18.PSDs of the building top displacement and acceleration in scheme B.

Figure 18 .
Figure 18.PSDs of the building top displacement and acceleration in scheme B.

Figure 19 .
Figure 19.Variation in the response STD with wind speed.
(a) Equivalent frequency (b) Equivalent damping ratio

Figure 20 .
Figure 20.First three frequencies and damping ratios of a linearized system.

Figure 19 .
Figure 19.Variation in the response STD with wind speed.

Buildings 2024 ,
14, x FOR PEER REVIEW 22 of 25 (g) Base shear force of the US.

Figure 19 .
Figure 19.Variation in the response STD with wind speed.
(a) Equivalent frequency (b) Equivalent damping ratio

Figure 20 .
Figure 20.First three frequencies and damping ratios of a linearized system.

Figure 20 .
Figure 20.First three frequencies and damping ratios of a linearized system.

Table 1 .
Selected DOFs of the simplified models for scheme A.

Table 1 .
Selected DOFs of the simplified models for scheme A.

Table 2 .
Selected DOFs of the simplified models for scheme B.

Table 2 .
Selected DOFs of the simplified models for scheme B.

Table 3 .
Structural response STDs by the LrUs and L1U1 models (scheme A).

Table 4 .
Structural response STD according to the LrUs and L1U1 models (scheme B).

Table 4 .
Structural response STD according to the LrUs and L1U1 models (scheme B).