#### Abstract

This paper conducts a parametric study on the seismic response of multistorey bidirectional eccentric structures from elastic stage to inelastic stage. Based on a simplified multistorey bidirectional eccentric model composed of bidirectional lateral load-resisting members, a general law is proposed for three-stage natural frequency variation behaviour from elastic stage to inelastic stage of eccentric frame structures with different layers. Different simplification treatments are conducted on each stage and the three stable parameter analysis stages are defined. The corresponding dynamic stiffness matrices and motion equations in different loading stages are derived. On this basis, a parametric analysis of seismic response of a three-storey bidirectional regular eccentric structure from elastic stage to inelastic stage is conducted. Effects of the uncoupled torsion to lateral frequency ratios () and bidirectional eccentricities on the seismic responses are investigated. The results reveal that as increases, translational displacement in the load direction first decreases and then increases; meanwhile, the displacement perpendicular to load direction and torsion displacement first rise and then decrease sharply. When , the coupling effect between the translation in the load direction and the torsion is at its strongest condition. Increasing the eccentricities leads to a decrease in the displacement in the load direction as well as an increase in the displacement perpendicular to load direction and torsion displacement. Variation regularity of inelastic seismic response is remarkably different from that in elastic stage. The lateral-torsional coupling effect of the bidirectional eccentric structure is closely related to both the period ratio and the bidirectional eccentricities.

#### 1. Introduction

The trends toward developing modern buildings with complex and irregular shapes have led to serious concerns with problems associated with structural irregularities. In these irregular structures having an uneven distribution of mass or stiffness, the centers of resistance do not coincide with the centers of mass. When these structures, so-called eccentric buildings, are subjected to earthquake ground motions, they undergo both translational and rotational motions, even though the earthquake excitation is purely unidirection translational. Observations following seismic events have revealed that these eccentric structures tend to suffer from structural damage more seriously than noneccentric ones because of the influence of lateral-torsional coupling [1, 2]. The lateral-torsional coupling behaviours of eccentric structures have since attracted continuous attention from scholars and engineers alike, and many research results have been reported [3–6].

The initial research focused on the single-storey eccentric structure. The system parameters of a single-storey, one-way eccentric structure in an elastic earthquake response can be analytically determined. These parameters include the uncoupled torsion to lateral frequency ratio, the stiffness eccentricity of the floor, and the natural period of the structure [7, 8]. Kan and Chopra [9] conducted the first parametric studies on the inelastic response of single-storey one-way eccentric structure. They declared that the lateral-torsional coupling effect in inelastic stage depends significantly on uncoupled torsion to lateral frequency ratio. The system in inelastic stage is more and more close to the translation motion of a system with a single degree of freedom. In contrast to these findings, scholars such as Tso, Sadek, and Bzorgnia [10, 11] reported that eccentric system does not respond primarily translation in inelastic stage. The different or even contradictory conclusions provided by these studies are mainly attributed to the fact that the researchers adopted different analysis models or parameter calculation methods. Kan and Chopra [9] adopted a very simple model that consisted of a single lateral load-resisting member. Meanwhile, studies by Tso and Sdek [10] and Bozorgnia and Tso [11] utilized models consisting of three lateral load-resisting members. Kan and Chopra [9] and Tso and Sadek [10] defined and calculated the uncoupled torsion to lateral frequency ratio Ω based on a vertical axis around the mass center, while Bozorgnia and Tso [11] selected a vertical axis around the stiffness center. From the above literature, it is known that the study of single-storey eccentricity structure is more mature. As a result of the comprehensive literature, experts have arrived at a universal conclusion about lateral-torsional coupling in the elastic stage of single-storey eccentric structure. However, the studies of lateral-torsional coupling in the inelastic stage have not always got generally applicable conclusions. And there are more research and discussion on this aspect [12–14].

The seismic responses of multistorey eccentric structures in elastic and inelastic stages are more complex than the responses of single-storey, one-way asymmetrical buildings subjected to earthquake motion. Currently, the available investigations primarily concentrate on the dynamic analysis of these regularly asymmetrical multistorey buildings. Maruić et al. [15] studied the elastic and inelastic seismic response of a multistorey regular eccentric steel frame structure under bidirectional earthquake motions. They found that the reduction in lateral stiffness due to plastic deformations in a certain direction may considerably influence the lateral load-resisting member in the orthogonal direction. Magliulo and Ramasco [16] compared the seismic response of a multistorey one-way eccentric frame structure under the influence of unidirectional and bidirectional earthquakes. Their results indicated that the maximum base shear and the top displacement are not very sensitive to the presence of the orthogonal component, which, conversely, leads to a large increase in the column plastic excursions. Jiang and Kuang [17] investigated the influence of uncoupled torsion to lateral frequency ratio and bidirectional eccentricities on the normalized frequencies of multistorey bidirectional stiffness eccentric structures from elastic stage to inelastic stage. Bosco et al. [18] proposed new corrective eccentricity aiming at generalizing and simplifying the previous corrective eccentricity method, and these equations of corrective eccentricities are defined based on the response of a set of single-storey systems representative of buildings. This research lays the foundation for the application of nonlinear static method in seismic assessment of eccentric structures. Wang [8] took strength eccentricity as the inelastic impact factor and studied the inelastic seismic response of multistorey eccentric structures. They declared that under the influence of earthquakes with different intensities, the maximum displacement ductility requirement mostly occurred in stiff side members, where it was about two to four times the value found in the flexible side members. Scholars including Wu et al. [19] analyzed the seismic responses of asymmetric regular multistorey buildings in both the elastic and inelastic stages. They found that the ductility responses of members at the stiff side increase with the increase of eccentricities. By adopting a 3D plastic hinge model, which closely approximates actual frame structures, Stathopoulos et al. [20] conducted inelastic dynamic time history analysis of multistorey bidirectional eccentric structures. In contrast to related noneccentric systems, the so-called flexible side frames had more ductility requirements than that of stiff side frames; this finding was also different from previous results using the simplified one-storey model [21]. Rizwan et al. [22] utilized the SAP2000 to construct a plastic hinge model of a 3D frame structure and conducted nonlinear dynamic history analysis on two kinds of asymmetric structures. The first kind was a mass symmetric system, with strength, stiffness, and both strength and stiffness eccentricities. The second kind was a mass eccentric system. Rizwan’s results indicated that, for the first kind of system, the response can be better corelated to the strength eccentricity. And the second kind of system exhibited almost uniform ductility demand in beams of various frames for smaller eccentricities. Halabian et al. [23] studied mass eccentricity and stiffness eccentricity in terms of the inelastic seismic response of a multistorey bidirectional eccentric frame-shear structure and evaluated the torsion effect using the unified ductility requirement distributed along the structure’s height direction. The results demonstrated that, as for a system consisting of dual lateral load-resisting members, column at the corner of the flexible side has greater ductility requirements than that of a similar frame system.

It can be seen that [8, 19, 20, 22, 23] drew different conclusions, which can be attributed to the differences among the analytical models, parameter calculation methods, eccentricity types, and input seismic waves adopted in these studies. Reference [19] calculated the approximate floor stiffness eccentricity based on the first order of the stiffness of lateral load-resisting members in each layer. Meanwhile, [20, 22] conducted push-over analysis of each frame and obtained related force-displacement curves. By fitting the results to a second-line model, the elastic lateral load-resisting stiffness of each frame was obtained, and the stiffness eccentricity was calculated. Reference [8] primarily analyzed the influence of mass eccentricity on the inelastic response, while [20, 22, 24] demonstrated that the inelastic response of the mass eccentricity structure differs from that of the stiffness eccentricity structure. Therefore, the findings of [8] require further validation in order to determine if they are applicable to multistorey structures with stiffness eccentricity.

In addition to the dynamic time history analysis, some experts have expanded the parameter analysis, but most of them are confined to the elastic stage. Chopra et al. [7, 25] conducted elastic parametric analysis of the seismic responses of multistorey eccentric structures using the mode-superposition response spectrum method. Their results indicated that the effects of lateral-torsional coupling on the response of a multistorey eccentric building and its associated one-storey system are similar. Duan and Chandler [26, 27] further reported that a single-storey structural model is not on its own sufficient to investigate completely the inelastic torsional effects in multistorey regularly asymmetric buildings which are unlike elastic studies. Anagnostopoulos et al. [28] pointed out that conclusions based on the widely used simplified, one story, eccentric systems of the shear-beam type, to actual, nonsymmetric buildings and consequent assessments of the pertinent code provisions, can be quite erroneous, unless special care is taken to match the basic properties of the simplified models to those of the real buildings. It can been seen that the lateral-torsional coupling behaviour is well understood within the limits of elastic deformations but is still far from being fully understood in inelastic stage even for a special multistorey regular eccentric structure. Studies show that the lateral-torsional coupling behaviour in inelastic stage is obviously different from that in elastic stage [17]. After entering the inelastic stage, the lateral load-resisting members experience buckling, unloading, and reverse loading under the influence of earthquake motion. Under these circumstances, the locations of stiffness centers continue to move, and the stiffness center of each layer no longer remains in the same vertical axis direction, creating difficulties for investigation and analysis. Therefore, the multistorey eccentric model should be adopted in the inelastic parameter analysis so as to obtain some regular conclusions.

A thorough review of the literature summarized above leads to the conclusion that no consensus has been reached concerning the inelastic behaviours of multistorey eccentric structures. Only by conducting systematic parametric analysis can we obtain the variation regularity of seismic response in the inelastic stage. The primary aim of this paper is to present a systematic parametric study on the variation of lateral-torsional coupling behaviour of multistorey bidirectional eccentric structure, thereby providing a basis for intensive studies leading to their effective improvement on other aspects. Therefore, based on the proposed multistorey eccentric simplified model composed of bidirectional lateral load-resisting members, this paper outlines a general law for the three-stage natural frequency variation behaviour from elastic to inelastic stage of eccentric structures with different layers. Different simplification treatments are conducted in these three stages; accordingly, three stable parameter analysis stages are defined, and the motion equations of different stages are deduced. The structural seismic response of each stage is gained by solving the motion equation in each stage. On this basis, parametric analyses on seismic responses of multistorey bidirectional regular eccentric structures are conducted from elastic to inelastic stage. Influences of uncoupled torsion to lateral frequency ratios and bidirectional eccentricities on the structural seismic response are investigated, and a general law of seismic response is obtained, changing with different parameters in all three stages. Finally, the variation regularity of torsion-translation displacement ratio with the period ratio under different eccentricities is discussed. The results can provide references for the design and earthquake-proofing of relevant structures.

#### 2. Analytical Model and Motion Equation

##### 2.1. Establishment of the Eccentric Simplified Model

A reinforced concrete bidirectional eccentric frame structure with spans and* N*-storey is illustrated in Figure 1. Only the stiffness eccentricities caused by the irregular planar distribution of columns are considered in this example, and no changes are presented along the vertical direction of the columns; in other words, the structure investigated herein is a regular eccentric structure.* CM* and* GC* stand for the mass center and geometric center of each layer, respectively;* CS* stands for the stiffness center, with a coordinate of (, ). A spatial shear-torsional type series rigid plate layer model with bidirectional lateral load-resisting members is introduced in order to conveniently conduct parametric analysis of structural models with different uncoupled torsion to lateral frequency ratios or stiffness eccentricities. The assumptions of this model, shown in Figure 2, are as follows:(1)The mass center (*CM*) of each floor distributes in one line and is located at the geometric center (*GC*) of the slab. The radius of gyration of the mass center of each floor is the same.(2)The masses of slab and columns of each floor are concentrated at the mass center of that floor.(3)The degrees of freedom of each floor are located at the mass center of the slab. Each floor has three degrees of freedom: translation degrees of freedom in the and directions as well as a torsion degree of freedom around the vertical axis.(4)The slab is rigid within its plane, and the out-of-plane stiffness is small and can be ignored.(5)Out-of-plane stiffness and torsional inertia of each lateral load-resisting member are assumed to be negligible and thus ignored in the analysis.

The slab for each floor is supported by four massless, axially inextensible load-resisting members: frame A, frame B, frame C, and frame D. By combining the frame in the direction of the* j*th floor of the original structure to double sides, frames A and B are formed, which provide lateral stiffness in the direction. Similarly, by combining the frame in the direction of the* j*th floor of the original structure to double sides, frames C and D are formed, which provide lateral stiffness in the direction. Frame A and frame B are located at equal distance from, but on opposite sides of,* x*-axis; frame C and frame D are located at equal distance from, but on opposite sides of,* y*-axis. Compared to the* j*th floor of the original model, the* j*th floor of the simplified model has the same total lateral stiffness on the horizontal axis, the same torsion stiffness about the mass center, and equal bidirectional eccentricities. The calculation formulas of the simplified model are listed as (1)–(3). The simultaneous solution of six formulas in (1)–(3) yields the parameters of the simplified model, which are listed in (4)–(6).where* k*_{xji} and* k*_{yji} denote the lateral stiffness of the* i*th column of floor in the and directions, respectively; and are, respectively, the and coordinates of the* i*th column of floor* j*; = and = represent the total translational stiffness of floor in the and directions, respectively; = + denotes the torsional stiffness about the mass center of floor* j*; and designate the and coordinates of the stiffness center (*CS*) of floor . As in multistorey buildings, the center of rigidity cannot be exactly defined, an approximate* CS* of floor is computed herein for reference purposes, = = [23]. The normalized stiffness eccentricities of floor are defined as , ; is the radius of gyration of the mass center of floor.* K*_{xjb},* K*_{xja},* K*_{yjd}, and* K*_{yjc} are the lateral stiffness of each frame in the* j*th floor of the simplified model.

It can be seen from (1)-(2) that and of floor of the original model are kept constant in the elastic stage; thus, and also remain constant. Consequently, and of floor of the simplified model do not vary. When the structure enters the inelastic stage, and vary with the load, and correspondently, and will also vary. It is suggested that when conducting inelastic analysis by adopting a simplified model in replace of the original structure, the values of and can be adjusted so that the simplified model can be used to accurately replace the original model for inelastic analysis.

##### 2.2. Solving the Motion Equation

The stiffness matrix of an* N*-storey bidirectional eccentric structure can be acquired based on the direction stiffness method [29], which also allows its motion equation to be deduced. Take a three-storey bidirectional eccentric structure serving as an example, the motion equation under an earthquake ground motion in the direction can be expressed aswhere , , , and are the mass matrix, stiffness matrix, displacement matrix, and transformation matrix, respectively, as described in [22]. The damping matrix and share the same modality. That is to say, is assumed to be proportional to the stiffness matrix . It is defined byin which , is the first natural angular frequency, and is the corresponding damping ratio. According to the provisions by Code for Seismic Design of Building [30], is advisable to 0.05.

Equation (7) is a differential equation in the time domain, and its calculation process is tedious. Therefore, it is necessary to transform it into a linear equation in complex frequency domain. is defined as the Fourier transformation of the displacement . Thus, the Fourier transform is applied to (7), and the following results can be obtained:

Equation (9) can be solved using the modal analysis method. However, because there are nine degrees of freedom, the solution cannot be explicitly expressed. Therefore, software programming is required to calculate the natural frequency and related vibration mode via numerical method. By expanding the structural seismic response in complex frequency domain with its inherent vibration mode as the base, there iswhere is the column vector of modal displacements of structure and is the modal displacement of* k*th modal.

Assuming that a structure on a rigid foundation possesses classical normal modes, its natural frequency and vibration mode satisfy the following orthogonality conditions [31]:where is N x N identity matrix, and are, respectively, the natural frequency of* k*th mode and the damping ratio of* k*th mode, and denotes a diagonal matrix.

Introducing the transformation of (10) and applying the orthogonality conditions (9), when premultiplied by , become

Defining the modal structure transfer function

is also called the complex frequency response function. The modal displacement can be obtained from (12).

Substituting (15) into (10), the structural seismic response in complex frequency domain can be obtained.

The ratio of the displacements in complex frequency domain at the mass center of each floor to Fourier transform of the ground motion displacement, denoted as , is referred to as the displacement transfer function. And is a physical parameter representing its own characteristics and only related to the frequency of earthquake motion , regardless of the type of earthquake motion. As the ratio reflects the dynamic response of the system under unit harmonic excitation, the transfer function can represent the structural dynamic response when given the same input. The interstory displacement is defined as follows:

and are, respectively, the interstory translation displacement transfer functions in the and directions of floor* j*; is the interstory torsion displacement transfer function of floor* j*.

#### 3. Analysis Methods

##### 3.1. Natural Frequency Analysis of Multistorey Eccentric Frame Structures

This analysis takes a typical 5×3 spans three-storey bidirectional regular eccentric frame structure with spans of 6.9 m and 6.3 m as an example. Its plane layout is shown in Figure 1. And its layer height is 3.9 m. The columns in Figure 1 marked as KZ1 have dimensions of 0.5 m × 0.55 m (×, similarly hereinafter), with a total column number of 6. The remaining frame columns marked as KZ2 have dimensions of 0.40 m × 0.45 m, with a total column number of 18. The cross-sections of the longitudinal beams are 0.25 m × 0.6 m, the cross-sections of the transverse beams are 0.25 m × 0.55 m, and the dimensions of the secondary beams are 0.2 m × 0.4 m. The thickness of the floor slab is 0.1 m. The concrete strength of the floor, beam, and column is C30. Moreover, there is no variation in each column along the vertical direction. That is, the structure investigated herein is a regular eccentric structure. Applying the method described in Section 2.1 yields the simplified model of the eccentric frame structure. In order to carry out the full-process analysis from elastic to inelastic stage, the original structure is treated as a shear frame structure, and the corresponding column is considered to be a member with shear deformation as the main part. That is to say, the restoring force model used in this article is based on the member. The force-displacement relationship of each lateral load-resisting member-both in the and directions of the original and simplified models is simulated by the bilinear hysteretic model as shown in Figure 3. denotes the stiffness of the member in the elastic stage; designates the stiffness after yielding, taken as 10% of the elastic stiffness, and it is tangent stiffness [32]; and is the yield displacement, which can be calculated according to empirical equations [33].

The nonlinear static analyses are, respectively, carried out on the original three-storey bidirectional eccentric model and corresponding simplified model. Through gradual loading, inverted triangle static loads in the direction are applied on the mass centers of each floor. Nonlinear problems are settled based on the load increment method. The eccentricities and the displacement of each floor are calculated at each step. Then, the displacement of each member is calculated according to its coordinate position. Next, the stiffness of each member is also modified according to the its displacement and related restoring force model at each step in order to form a new stiffness matrix for the next step analysis. Therefore, the nonsynchronization of the displacements of the members, which is induced by the torsion of the floor, can be considered in detail during analysis. When the story drift angle of a certain layer reaches 1/50 of the inelastic story drift angle limit of frame structures dictated in the Code for Seismic Design of Building [30], loading is stopped. The entire analysis process is achieved by M language programming in MATLAB.

The normalized stiffness eccentricities* b*_{x1} and* b*_{y1} of the original model and the corresponding simplified model are evaluated in Figure 4, where the loading coefficient *α* of the horizontal axis is defined as the proportion of the applied load to the entire final load on the structure, with a maximum value of 1. It can be noted from Figure 4 that the variation curve of* b*_{x1} with *α* has no mutation points, indicating that the members perpendicular to the load direction on the first floor do not yield. Conversely, mutation points appear on the variation curve of* b*_{y1} with *α*, suggesting that the members in the load direction on the first floor gradually yield. Variations of* b*_{y2} and* b*_{x2} are respectively similar to those of the first floor. Moreover,* b*_{y3} and* b*_{x3} remain straight lines with *α*, which indicate that all the members on the third floor do not yield. It also can be seen from Figure 4 that, during the entire loading process, variation of* b*_{y1} of the simplified model matches favorably with that of the original model. Only during the yielding process of some members, there are some differences between the original model and the corresponding simplified model. For the overall analysis of lateral-torsional coupling behaviour of the eccentric structure in this paper, the impact of this difference is within the acceptance of engineering applications.

The first-order natural frequencies of the two models are illustrated in Figure 5. As shown in Figure 5, the first-order natural frequency of the eccentric frame decreases from 2.569 Hz to 0.873 Hz. At the same time, the second- and third-order natural frequencies also experience reductions in amplitude smaller than that of the first-order natural frequency. Throughout the entire loading process, the first three natural frequencies of the two models agree well with each other; the elastic stages of the first-order frequency of the two models have a difference of 0.20%, while the maximum deviation in the inelastic stage is 8.25%. These results indicate that, in the elastic stage, the simplified model can effectively fill the role of the original model during analysis. As for the inelastic stage, there are only minor differences during the yielding process of some lateral load-resisting members, but the deviation is still acceptable for engineering. The stable variation of the natural frequencies reflects stable variation of the structural integral stiffness, which indicates that the structural stiffness changes of the two models are consistent in the process from elastic to inelastic stage. Therefore, the proposed simplified model is also suitable for analysis in the inelastic stage, and it can reliably replace the original model for full-process parameter analysis.

It can be noted from Figure 5 that the natural frequency reduces gradually through the whole process from elastic to inelastic stage, with the variation occurring in three obvious stages. The first stage takes the form of a horizontal straight line, the second stage is a huge reduction segment like a diagonal line, and the third stage is an approximately horizontal segment with relatively small variation amplitude. When the loading coefficient is larger than 0.92, the structure is close to sustaining severe damage; therefore, the analysis does not consider this stage. To further validate the variation tendency of the natural frequency of multistorey eccentric frame structures, typical five- and seven-storey structures are selected to analyze vibration regularity. The planar layouts of five- and seven-storey regular eccentric frame structures are similar to that displayed in Figure 1, with a floor height of 3.9 m. The cross-section dimensions of the large columns (KZ1 frame columns) of five-storey eccentric frame structure are 0.60 m × 0.65 m, while those of the small columns (KZ2 frame columns) are 0.50 m × 0.55 m. The cross-section dimensions of the large and small columns of seven-storey eccentric frame structure are 0.65 m × 0.70 m and 0.55 m × 0.60 m, respectively. The simplified models of eccentric frame structures with five and seven layers can be obtained according to the method described in Section 2.1. The force-displacement relationship of each member of the original and simplified models of the five- or seven-storey structures is also described by the bilinear hysteretic model shown in Figure 3. The corresponding parameters (, , and ) are calculated by the same method as the three-storey model. Full-process nonlinear static analyses can then be conducted with these two models. The variation curves of the first-order natural frequencies are shown in Figures 6 and 7.

A comparison of Figures 5, 6, and 7 indicates that the frequency variation regularities of five- and seven-storey eccentric frame structures is similar to that of the three-storey eccentric frame structure; that is, all of the curves are characterized by three different stages. Although the three-, five-, and seven-storey eccentric structures have the same structure plane, their cross-sections of the members are different. The corresponding elastic lateral stiffness and yielding displacement values also differ from each other, as do the initial bidirectional eccentricities and uncoupled torsion to lateral frequency ratios. That is to say, the three examples analyzed above represent common multistorey frame structures, and the three-stage variation regularity represents the general law of multistorey eccentric frame structures. This finding agrees with the frequency variation regularity obtained in other types of analysis [34].

The frequency variation amplitude is largest in the second stage, during which time it also has its maximum influence on the structure. This indicates that the second stage is the most important component for analysis. According to the second stage variation tendency of the natural frequency in different layers, the frequency during this stage is fitted into a slant line with a linear relationship to the loading coefficient via the least-squares method, making the analysis procedure more convenient. In the third stage, the frequency variation amplitude is very small. During analysis, the average frequency value can be selected for the third stage to simplify the calculation process. The general form of the three-stage frequency can be expressed by where is the first-order natural frequency in the elastic stage of the structure;* k* and are, respectively, the monomial coefficient and constant term of the linear variation of frequency in the second stage; and* f*_{e3} is the third stage frequency, with its value is equal to the average frequency of this stage. The specific values of the first-order frequency in three stages for the three-storey eccentric frame structure are given in in which , , and are, respectively, the values of the simplified three-stage of the first-order natural frequency. Figure 8 illustrates the variation curves of the simplified three-stage frequency and the original frequency.

Similar methods can be adopted to obtain corresponding simplified three-stage frequencies for frame structures with other numbers of layers.

##### 3.2. Description of the Full-Process Three-Stage Analysis

Starting from the three-stage variation features of the natural frequency, three stable parametric analysis stages are defined as the first, second, and third stage. Using the first-order frequency curve in Figure 8 as an example, any one point in the first and third stages can represent the state of the entire segment. Meanwhile, because the second stage is a slant line, selecting any one point as a representative allows the rest of the points to be obtained based on the linear relationship. It can be noted that the structure is in elastic during the first stage parametric analysis, without any lateral load-resisting member yielding. During the second stage, the lateral load-resisting members on the first floor yield gradually along the load direction. For different values of *α*, the related frequency and the total lateral stiffness of the first floor are also different, implying different degrees of inelasticity development in this stage. In the second stage, a point near the center is selected for representativeness. In the third stage, members on the first floor all yield along the load direction, while those on the second floor yield gradually. At the same time, the lateral load-resisting members on the third floor remain elastic. Through the above analysis, we can know that the variation curves of* b*_{x1},* b*_{x2}, and* b*_{x3} with *α* all stay straights. It shows that, under the load in the direction, all the members perpendicular to the load direction of the system studied in this article are still in the elastic stage when the loads are limited by = 1/50. Hence, in order to carry out the full-process parametric analysis, the following analyses only consider the elastic stage of the members perpendicular to the load direction. The story stiffness of each floor in each stage together with other parameters of the entire system in that stage is substituted into the stiffness matrix so as to determine the stiffness matrix in each stage. The motion equation of each stage can also be acquired by substituting the stiffness matrix into (7). The motion equation of each stage can be solved by (9)–(16) and parametric analysis of each stage can be conducted.

#### 4. Parametric Analyses of the Structural Seismic Responses

Scholars pay special attention to the amplitude of interstory displacement of multistorey eccentric structures, and the interstory displacement can determine the deformation distribution of each floor. Therefore, the peak values of the interstory translational displacement transfer functions in the and directions, that is, and , are studied. And the peak value of the interstory torsion displacement transfer function is also investigated. This section discusses the variation regularity of these peak values with uncoupled torsion to lateral frequency ratios and bidirectional eccentricities. For the convenience of discussion, the peak value of the displacement transfer function is called “displacement.” Here, Ω is defined as the ratio of the first-order pure torsional frequency to the first-order pure translational frequency in the direction of the corresponding symmetric system. For a certain and , different ratios can be generated by adjusting the distribution of the lateral load-resisting members in the direction and the value of . At the same time, the distribution of the lateral load-resisting elements in the direction and the total lateral stiffness are kept as constant. That is to say, a set of Ω can be generated by changing* K*_{yjc},* K*_{yjd}*,* and . Meanwhile,* K*_{xja},* K*_{xjb},* b*, , and and the sum of* K*_{yjc} and* K*_{yjd} are kept invariant. There are nine degrees of freedom for the three-storey bidirectional eccentric simplified model analyzed in this paper, and the explicit solution cannot be given directly when Ω is adjusted. Therefore, it is necessary to use the numerical calculation method of developed MATLAB program. With MATLAB’s M program language, the corresponding* K*_{yjc},* K*_{yjd}, and can be easily obtained when some other parameters are kept constant. Then, the value of Ω of the symmetric system corresponding to the eccentric model can be obtained by using the “eig” function trial calculation in MATLAB. The corresponding parameters are constantly adjusted according to this method until the desired value of Ω is transferred. The simplified models with different values of and can be established by adopting the approach specified in study [17].The range of Ω is 1.0 to 1.8, and the range of and is 0.1 to 0.4 [35].

##### 4.1. The Effect of on Interstory Displacement

Variations of interstory translation displacements in the direction with Ω from the first to the second and then to the third stage are analyzed in Figure 9. It is evident from Figure 9(a) that the interstory translation displacement in the direction decreases firstly and then increases with the adding of Ω in the first stage. The turning point is at Ω=1.1. When Ω is greater than 1.4, the increase amplitude is only 1.71% to 5.33%, tending to stable. Massive calculations reveal that, within the range of 1.0 to 1.4, Ω plays the dominant role in the translation displacement along the direction when is less than or equal to 0.2. When is greater than 0.2, the key range of Ω is between 1 and 1.6. Compared to the first stage, the second and third stages (Figures 9(b) and 9(c)) see smaller magnitudes of changes in the translational displacements, and the variation of the interstory translation displacement in the direction with Ω is no longer consistent. The influence of Ω on the translation displacement is reduced by the in-depth development of structural inelasticity, and the variation tendencies of displacements of each floor with Ω are also affected to a certain degree.

**(a) The first stage**

**(b) The second stage**

**(c) The third stage**

The effects of Ω on the interstory translation displacement in the direction and the torsion displacement in different stages are, respectively, expressed in Figures 10 and 11. As Figures 10(a) and 11(a) imply, the translation displacement in the direction and the torsion displacement increase firstly and then decrease sharply with the increase of Ω in the first stage, with the turning point occurring at Ω=1.1. From the first to the second and then to the third stage, the reduction amplitude of the translation displacement in the direction becomes smaller and smaller, while that of torsion displacement in the second stage remains relatively large (Figure 11(b)). Moreover, the torsion displacement demonstrates a tendency toward slow increase with relatively large values of Ω in the third stage.

**(a) The first stage**

**(b) The second stage**

**(c) The third stage**

**(a) The first stage**

**(b) The second stage**

**(c) The third stage**

Comparing Figures 9, 10, and 11 reveals that, in the first stage, increasing Ω causes the torsion displacement to increase when Ω is less than or equal to 1.1. When Ω increases to 1.1, the torsion displacement reaches its maximum value, and the translation displacement in the direction is at its minimum. The coupling effect between the translation in the direction and the torsion becomes quite strong, which is extremely unfavorable to the torsion resistance design. With the increase of Ω (Ω>1.1), the structural seismic response gradually approaches the pure translation along the load direction, and the coupled torsional displacement and the translational displacement perpendicular to the load direction become increasingly smaller. The coupling effect between the translation in the load direction and the torsion is substantially reduced. The above finding is due to the fact that the eccentric structure with a larger Ω has a larger torsional stiffness and a larger ratio of torsional stiffness to lateral stiffness; thus, the lateral-torsional coupling effect and the torsional effect will be weak. That is to say, the structural seismic response will gradually approach the pure translation along the load direction with the increase of Ω. Therefore, it will be good practice in design to achieve a large Ω value, thereby avoiding torsion failure.

##### 4.2. The Effect of on Interstory Displacement

Figure 12 presents the effects of on the interstory translation displacements in the direction from the first to the second and then to the third stage. It can be noted from Figure 12(a) that increasing causes translation displacement in the direction to decrease uniformly. Specifically, the average reduction percentage in the translation displacement in the direction achieves 5.43% when increases 0.1. After entering the second stage (Figure 12(b)), translation displacement in the direction on the first floor decreases gradually, while those on the second and top floors still increase. In the third stage (Figure 12(c)), translation displacement in the direction on the first floor slightly increases, while those on the second and top floors decrease slowly, in total contrast to the second stage. These findings indicate that, in the first stage, variation tendencies of the interstory translation displacements in the direction with are similar, while different inelasticity development stages have comparably large influences on the variation regularity.

**(a) The first stage**

**(b) The second stage**

**(c) The third stage**

The effects of on the interstory translation displacement in the direction and the torsion displacement in different stages are, respectively, evaluated in Figures 13 and 14. It can be seen from Figures 13(a) and 14(a) that the translation displacements in the direction and the torsion displacement increase with the increase of in the first stage, and the increase percentages reduce slowly. Compared to the first stage, increasing exerts greater effects on the translational displacement in the direction and torsional displacement in the second stage when is approximately greater than or equal to 0.3. For example, as increases from 0.3 to 0.4, the variation percentages of the translational displacement in direction and the torsional displacement are, respectively, 4.1% and 17.9% in the first stage. In the second stage, the variation magnitudes increase to 57.9% and 65.6%, respectively. In the third stage, the displacement variation amplitude is clearly reduced. Figures 13 and 14 also indicate that the variation tendencies for interstory translation displacement in the direction and torsion displacement with are similar on all floors. During the second stage, the development of structural inelasticity enhances the increase amplitude of the displacements. In the third stage, the influence of on displacements is reduced.

**(a) The first stage**

**(b) The second stage**

**(c) The third stage**

**(a) The first stage**

**(b) The second stage**

**(c) The third stage**

##### 4.3. The Effect of on Interstory Displacement

The effects of on the interstory translation displacements in the direction in the different stages are shown in Figure 15. It is evident from Figure 15(a) that the translation displacement in the direction gradually decreases as increases in the first stage, and the reduction amplitude rises. Figures 15(b) and 15(c) show that increasing exerts smaller effects on the variation of translation displacement in the direction in the second and third stages than in the first stage. It also can be noted from Figure 15(c) that the translation displacement in the direction of the second floor slowly increases and those of the other two floors slightly decrease with the increase of . These findings demonstrate that the variation tendencies of interstory translation displacements in the direction with are similar in the first stage. With the yielding of the lateral load-resisting members, the influences of on the displacements are weakened in the second and third stages. And the variation regularity of interstory translation displacements in the direction is no longer consistent.

**(a) The first stage**

**(b) The second stage**

**(c) The third stage**

The influences of on the interstory translation displacement in the direction and the torsion displacement in different stages are given in Figures 16 and 17. It is evident from Figure 16(a) that increasing causes translation displacement in the direction to increase in the first stage, and the corresponding increasing percentage gradually decreases. Meanwhile, the torsion displacement gradually increases with the increase of and tends to be stable implied in Figure 17(a). Compared to the first stage, the second stage sees smaller variation percentages of translation displacement in the direction and the torsion displacement. It can be seen from Figure 16(c) that the translation displacement in the direction increases slightly with the increase of in the third stage. Figure 17(c) shows that the torsion displacements of the first and second floors first reduce and then increase slowly, and the torsion displacement of the top floor decreases gradually, but the variation amplitude is relatively small in this stage. These findings imply that the development of structural inelasticity weakens the influences of on displacements and changes the variation tendency of interstory torsion displacement with .

**(a) The first stage**

**(b) The second stage**

**(c) The third stage**

**(a) The first stage**

**(b) The second stage**

**(c) The third stage**

The parametric analysis of lateral-torsional coupling effect of multistorey bidirectional eccentric structures has indicated that the first-order natural frequency of the eccentric structure becomes closer to that of the corresponding symmetric system from the first to the second and then to the third stage, and the lateral-torsional coupling effect in the first vibration mode is significantly reduced [17]. Accordingly, the first-order vibration mode tends toward the pure translation vibration mode in the load direction. This also determines that the translation displacement in the load direction becomes greater during this process; the coupled torsional displacement and the translational displacement in the direction become increasingly smaller (Figures 9–17). In particular, the translation displacement in the direction of the first floor presents the largest increase amplitude from the first stage to the second stage. For a system with Ω=1.4 and ==0.2, its translation displacement in the direction of the first floor increases 53.6% from the first to the second stage and only 14.9% from the second to the third stage. The above finding is due to the fact that the members in the load direction of floor one firstly yields in the second stage and further yield in the third stage, while members on other floors in the load direction and the members perpendicular to the load direction do not yield. The yielding of the members leads to a considerable reduction on the overall lateral stiffness in the load direction of floor one. As a result, the first floor becomes increasingly weak compared to other floors. The peak values of the interstory translational displacement transfer function in the load direction of floor one increase significantly. In the third stage, all the lateral load-resisting members in the load direction of floor one enter yield, and the members in the load direction of floor two begin to yield; therefore, the increase amplitude of the translation displacement of floor one in the load direction reduces. The overall structural seismic response gradually approaches the pure translation in the load direction, and the coupled torsional displacement and the translational displacement perpendicular to the load direction become relatively smaller. The coupling effect between the translation in the load direction and the torsion is greatly reduced.

Through the parametric analysis of the structural displacement, a general law for the seismic response of multistorey bidirectional eccentric structures varying with different parameters in the elastic and inelastic stages is obtained. These research results can provide references for structural and seismic design in engineering.

In the process of inelastic development, the actual damping of the structure is increased [36]. The damping problem in inelastic stage needs further research and discussion in basic theory and experimental verification, and this is another complex research topic. So the damping is assumed to be in direct proportion only to the linear portion of stiffness in this paper. When the stiffness-proportional damping is adopted, the actual damping in the second and third stages will be lowered, which will increase the response of the structure. This paper focuses on conducting a parametric analysis of the whole process from elastic to inelastic stage. And the influence regularities of uncoupled torsion to lateral frequency ratios and bidirectional eccentricities on the structural seismic response are mainly investigated. Due to the reason of damping change, the responses of the structure in the second and third stages are larger than the actual reactions. However, for the study of the variation law in the parametric analysis, this effect is not particularly sensitive.

#### 5. Discussion of Torsion-Translation Displacement Ratio under Different Eccentricities

The parametric analysis above indicates that the bidirectional eccentricities have a considerable influence on the translation and torsion displacements of the eccentric structures, and they also affect the ratio of torsion displacement to translation displacement. This section primarily focuses on the ratio of torsion displacement to translation displacement in the load direction, and torsion-translation displacement ratio of the top floor is defined as follows:

The torsion-translation displacement ratio of the top floor varying with coupled period ratio in the first stage is shown in Figure 18. is the first natural period of vibration mode which is chiefly torsion, and is the first natural period of vibration mode which is primarily translational movement. According to the findings in Figure 18, an increased with any and causes the displacement ratio to increase. And the lateral-torsional coupling effect of the eccentric structure becomes stronger. It also can be noted from Figure 18 that the displacement ratio changes significantly with the bidirectional eccentricities and when is the same. In other words, values of of eccentric structures with different bidirectional eccentricities are remarkably different from each other when the displacement ratio is the same. The results reveal that the structure with larger values of and has a smaller . As implied in Figures 18(a) and 18(b), the displacement ratio gradually increases as increases from 0.1 to 0.4 when the coupled period ratio is constant. Similarly, it can be seen from Figure 18(c) that increasing causes the coupled period ratio to decrease when the displacement ratio is fixed. These findings indicate that the value of corresponding to a specific displacement ratio gradually decreases with the increase of and . Similar phenomena not discussed here also occur under other values of and .

**(a)**= 0.1

**(b)**= 0.2

**(c)**= 0.4Reference [37] studied the elastic lateral-torsional coupling effect of high-rise buildings with one-way mass eccentricity and the approximated calculation formula for the top floor relative torsion effect (torsion-translation displacement ratio) associated with eccentricity and was obtained. Based on these results, the authors proposed recommended period ratio limit values to control the relative torsion effect of one-way mass eccentricity structures. Related results have been adopted in Technical Specification for Concrete Structures of Tall Building (JGJ 3-2010) [38]. The results of this paper indicate that the lateral-torsional coupling effect of bidirectional eccentric structures is not only related to the period ratio, but also affected by the bidirectional eccentricities. High-rise buildings, which are more sensitive to the torsion effect, require additional attention.

Figure 19 displays the variation of coupled period ratios and noncoupled period ratios corresponding to different values of Ω. Obviously, the noncoupled period is the reciprocal of Ω, while the coupled period is obtained via calculation. It is evident from Figure 19 that the coupled period ratios of the structures with various eccentricities differ from the noncoupled period ratios. The greater the bidirectional eccentricities, the greater the difference between the coupled period ratio and noncoupled period ratio. The difference is most significant in structures with Ω=1, where the coupled period ratio is as much as 31.8% greater than the noncoupled period ratio. When Ω is greater than or equal to 1.1, the coupled period ratio is smaller than the noncoupled period ratio. As Ω increases, however, the difference reduces, although it remains notable. Even for eccentric structure with Ω=1.6, the difference still reaches 8.47% when =0.2 and =0.4.

**(a)**= 0.2

**(b)**= 0.2The preceding analysis clearly demonstrates the relationship between the coupled period ratio, noncoupled period ratio, and bidirectional eccentricities in multistorey bidirectional eccentric structures. Reference [37] suggested that, under the premise of satisfying the required displacement ratio when the period ratio is between 0.8 and 0.9 [38], the difference between the noncoupled period ratio and coupled period ratio is within 5%. Obviously, this one-way mass eccentricity analysis result is not appropriate for multistorey structures with bidirectional stiffness eccentricities.

#### 6. Validation of the Dynamic Response of Eccentric Simplified Model

##### 6.1. Construction of the 3D Refined Model

Section 3.1 has already analyzed and discussed the applicability of the simplified model in the inelastic stage. To further validate the accuracy and reliability of the proposed simplified model in the parametric analysis of the seismic response, a corresponding 3D refined model for the three-storey eccentric structure with its plane layout as Figure 1 is established. And the dynamic responses of the proposed simplified model are validated by the nonlinear dynamic time history analysis. The beam element of BEAM188 which can define the cross-section dimension is applied to simulate the frame beams and columns, and the shell element of SHELL63 is used to simulate the floor slab. The related structural parameters of the three-storey eccentric structure are the same as those in 3.1 Section. When the nonlinear dynamic time history analysis of the 3D refined model is conducted by ANSYS, the simulation process becomes difficult if the structure directly utilizes the force-displacement resilience model. Hence, the elastic-plastic constitutive models of the beam and column elements adopt the bilinearity isotropic hardening model (BISO) which follows flow rules related to Von Mises yielding. Only the elastic modulus before yielding, the tangent modulus after yielding, and the yield strength should be defined in ANSYS, and the tangent modulus after yielding is also set as 10% of the elastic modulus [39]. The stiffness matrix coefficient calculation of damping defined by ANSYS is the same as the damping coefficient defined in Section 2.2.

Since this study focuses on the overall analysis of the structure and considers the coupling effect between translational displacement and torsional displacement of the eccentric structure, an integral analysis model of reinforced concrete is selected. It is assumed that rebars are dispersed in the whole element and the element is considered as a continuous homogeneous material. The contribution of rebars to the whole structure is realized by increasing the material elastic modulus via the stiffness* EI* equivalence principle [40].in which* E* and are, respectively, the elastic modulus and cross-section inertia moment; superscripts and refer to the concrete part and rebar part, respectively. The moment of inertia of the concrete is calculated considering the uncracked cross-section.

The dimension of mesh should also be taken into consideration during the construction of the finite element model. When the mesh size is small, the calculation accuracy is high, but due to the high degrees of freedom, the required calculation time, calculation amount, and memory space are great accordingly. Meanwhile, when the mesh size is smaller than 5 cm, the calculation may become difficult. Similarly, the mesh size cannot be too oversized; otherwise, the calculation accuracy will be affected. Based on above aspects, the 3D refined model is constructed, as shown in Figure 20.

##### 6.2. Selection of Seismic Waves

The selection of seismic waves should meet the requirements of three seismic oscillation elements, i.e., spectral characteristics, time of duration, and amplitude. The spectral characteristics can be described by the seismic influence coefficient curves and determined according to the site classification and designed seismic grouping. Reference [38] stipulates that the last duration of seismic waves should not be less than five times of the building basic period and 15 s, and the time interval of seismic waves can be set as 0.01 s or 0.02 s. The NS component of the 1976 Tianjin earthquake record, the NS component of the 1940 EL-Centro earthquake record, and an artificial Tianjin wave are selected as the external load of dynamic analysis. The selected seismic waves include a main energy section lasting duration of 15 s and an interval of 0.02 s. The seismic fortification intensity in Tianjin is 7 degrees. According to the Code for Seismic Design of Building [30], the basic seismic acceleration is designed as 0.15g. To conduct analyses under frequently occurred earthquake and rare earthquake, the seismic waves are proportionally increased to maximum values of 0.55 m/s^{2} and 3.1 m/s^{2}.

The acceleration time history records of these three seismic waves (amplitude increased to 3.1 m/s^{2}) and the corresponding spectra are shown in Figures 21 and 22, respectively. The frequency components of these three seismic waves are quite different. The Tianjin wave is mainly composed of low frequency components, and the high frequency content is low. The frequency components of EL-Centro waves are moderate, while the high frequency components of Tianjin artificial wave are relatively high. The frequency components and preeminent frequencies of seismic waves have a great influence on the dynamic responses of structures. Therefore, the conclusions obtained by the selected waves with different frequencies components will be more representative.

**(a) Tianjin wave**

**(b) EL-Centro wave**

**(c) Tianjin artificial wave**

**(a) Tianjin wave**

**(b) EL-Centro wave**

**(c) Tianjin artificial wave**

##### 6.3. Comparison of the Dynamic Response of 3D Refined Model and Simplified Model

The dynamic time history analyses of both simplified model and 3D refined model are carried out under the action of the ground motion in the direction. And the dynamic responses of the two models are compared. Under the action of frequently occurred earthquake, the structure is in the elastic stage. And under the action of rare earthquake, the structure enters the inelastic stage. The time history curves of the top floor displacement at the center of mass under the action of EL-Centro wave are shown in Figures 23 and 24.* u*_{x3} and* u*_{y3} denote the top floor translational displacement in the and directions, respectively. And represents the top floor torsional angle. The corresponding peak values of the top floor displacement are given in Tables 1 and 2. The peak value of the top floor torsional displacement is equal to the peak value of the torsional angle multiplying the corresponding radius of rotation.

**(a)**The displacement in the direction

**(b)**The displacement in the direction**(c) The torsional angle**

**(a)**The displacement in the direction

**(b)**The displacement in the direction**(c) The torsional angle**

It is evident from Figure 23 that the displacement time history curve of the simplified model is in good agreement with that of the 3D refined model under the action of EL-Centro wave. Under the action of the other two kinds of seismic waves, the displacement time history curves of two models are also coincided with each other. Restricted to the length of the article, the corresponding displacement time history curves are not presented in the paper. It also can be seen from Table 1 that the peak values of displacement in the direction of the two models deviate within 5%. Meanwhile, the peak values of translational displacement in the direction and the torsional displacement of the two models have differences within 10%. Figure 24 displays that the displacement time history curves of the two models differ a little under the action of rare earthquake but still agree with each other generally. As evidenced in Table 2, the peak values of translation displacement in the direction differ within 10%. For the rest displacement peak values, beside two values deviate greater than 13%, the rest values all vary within 10%.

In this paper, the overall response of the structure under the action of ground motion in the direction is mainly analyzed. In the inelastic stage, the peak values of displacement in the direction and torsional displacement have relatively large errors, but their actual values are very small and have limited influence on the structure. The main structural response is the translation displacement in the direction. And the translation displacements in the direction of the two models all have minor errors in both the elastic and inelastic stages, which are within the acceptance of engineering applications. Hence, for the parametric analysis of the coupling effect between translation and torsional displacement of the eccentric structure under the external load in the direction conducted in this study, the simplified eccentric model proposed is accurate and reliable.

#### 7. Conclusions

Based on the simplified multistorey bidirectional eccentric model, this study conducted a full-process parametric analysis of structural seismic response from elastic stage to inelastic stage. It can be clearly seen that the lateral-torsional coupling behaviours in the inelastic stage are remarkably different from that in the elastic stage. The variations of structural seismic response with parameters are related to the degree of inelasticity development. The following specific conclusions can be drawn:(1)Taking a typical multistorey bidirectional regular eccentric frame system as an illustrative example, a simplified bidirectional asymmetric multistorey dual lateral load-resisting model is constructed. The simplified model is applicable to inelastic parametric analysis. The development of simplified model laid the foundations for the successive parametric analyses in the following work.(2)During the full process from elastic stage to inelastic stage, the natural frequency of eccentric frame structures with different layers exhibits three-stage variation regularity. Different simplification treatments were conducted on these three stages. The second stage, which plays a key role during analysis, is fitted into a slant line via the least-squares method with a linear relation to the loading coefficient. Based on this, three-stage analyses are defined to provide convenience for parametric analysis of different inelastic development stages. The corresponding dynamic stiffness matrices and motion equations in different loading stages are derived.(3)Based on the three defined analysis stages, the full-process parametric analysis of lateral-torsional coupling behaviour of the multistorey bidirectional eccentric structure is conducted for the first time from elastic to inelastic stage. The variation tendencies of the lateral-torsional coupling behaviours under the major parameters, including the uncoupled torsion to lateral frequency ratios and eccentricities with different values and stages, are analyzed. The general laws of each parameter affecting the lateral-torsional coupling behaviours in three stages are concluded. The analysis results reported that the coupling effect between the translation in the load direction and the torsion is the strongest when Ω is equal to 1.1, which is extremely unfavorable for the torsion resistance design.(4)From the first to the second and then to the third stage, the overall structural seismic response gradually approaches the pure translation in the load direction, and the coupled torsional displacement and the translational displacement perpendicular to the load direction become relatively smaller. The coupling effect between the translation in the load direction and the torsion is greatly reduced.(5)Compared to noncoupled period ratios, coupled period ratios are different under different eccentricities. The greater the bidirectional eccentricities, the greater the difference between the coupled period ratio and noncoupled period ratio. Analytical results indicate that the lateral-torsional coupling effect of bidirectional eccentric structures is closely related to both the period ratio and the bidirectional eccentricities. High-rise buildings, which are more sensitive to the torsion effect, require additional attention. It is suggested that the effect of eccentricities should be considered when the coupled period ratio is limited by the relevant specifications.(6)Through the nonlinear dynamic time history analysis of the proposed simplified model and the corresponding 3D refined model established on ANSYS software, the dynamic response of the simplified model is verified. Results show that the displacement time history curves of the two models agree well with each other, and the differences in the displacement peak values of the elastic and inelastic stages are within 5% and 10%, respectively. It is proved that the simplified model proposed in this paper is completely feasible for the parametric analysis of the overall response of the eccentric structure along the force direction, and the results are accurate and reliable.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

This work was financially supported by the National Science Foundation Project (Research Project no. 51278335), National Youth Science Foundation Project (Research Project no. 51208356), and the National Science Foundation Project (Research Project no. 51478312), China.