Parameters Optimization and Performance Evaluation of the Tuned Inerter Damper for the Seismic Protection of Adjacent Building Structures

Xiaofang Kang,Jian Wu,Xinqi Wang and Shancheng Lei

1Key Laboratory of Intelligent Underground Detection Technology,School of Civil Engineering,Anhui Jianzhu University,Hefei,230601,China

2Chengnan College,Changsha University of Science and Technology,Changsha,410015,China

ABSTRACT

In order to improve the seismic performance of adjacent buildings,two types of tuned inerter damper (TID)damping systems for adjacent buildings are proposed,which are composed of springs,inerter devices and dampers in serial or in parallel.The dynamic equations of TID adjacent building damping systems were derived,and the H2 norm criterion was used to optimize and adjust them,so that the system had the optimum damping performance under white noise random excitation.Taking TID frequency ratio and damping ratio as optimization parameters,the optimum analytical solutions of the displacement frequency response of the undamped structure under white noise excitation were obtained.The results showed that compared with the classic TMD,TID could obtain a better damping effect in the adjacent buildings.Comparing the TIDs composed of serial or parallel,it was found that the parallel TIDs had more significant advantages in controlling the peak displacement frequency response,while the H2 norm of the displacement frequency response of the damping system under the coupling of serial TID was smaller.Taking the adjacent building composed of two ten-story frame structures as an example,the displacement and energy collection time history analysis of the adjacent building coupled with the optimum design parameter TIDs were carried out.It was found that TID had a better damping effect in the full-time range compared with the classic TMD.This paper also studied the potential power of TID in adjacent buildings,which can be converted into available power resources during earthquakes.

KEYWORDS

Adjacent buildings;tuned inerter damper(TID);H2 norm optimization;vibration control;energy harvesting

1 Introduction

Since Lee et al.first proposed the method of vibration control for civil structures in 1972[1],the theory and method of vibration control have been developed rapidly.In order to deal with structural vibration problems caused by wind load or earthquake load,a series of control strategies have been formed,including active control,passive control,hybrid control,and intelligent control [2].Among them,passive vibration control is widely used in engineering practice because of its high stability,low cost,good control effect,and easy realization[3,4].The design of classic tuned mass damper(TMD)consists of a mass block,spring and viscous damper,which is a reliable structural vibration control device [5].Kareem et al.took high-rise buildings subjected to wind load as an example,introduced the application of TMDs in Australia,Canada,China,Japan,and the United States in recent years,and showed that TMD had superior seismic performance [6].However,for TMD,to obtain better and more robust control performance,it is necessary to use a larger auxiliary mass,that is,to increase the mass ratio of TMD to the controlled structure,which is often restricted by site conditions or costs[7].In order to solve the problems encountered in the practical application of TMD,researchers have developed many new dampers based on the existing TMD research,including electromagnetic inertial mass damper(EIMD)[8,9],rotational inertial damper(RID)[10,11],rotational inertial double-tuned mass damper (RIDTMD) [12–14],tuned mass damper inerter (TMDI) [15–17],tuned liquid inerter system (TLIS) [18],tuned inerter damper (TID) [19–23],tuned mass-inerter damper (TMID) [24],tuned fluid inerter(TFI)[25,26],clutching inerter damper(CID)[27],etc.These dampers introduced a promising passive vibration control device,namely inerter device.

Inerter device is a mechanical device with two ends which was first proposed by Smith in 2002[28].Its characteristic is that the applied inerter force is proportional to the relative acceleration between the two ends,and it usually consists of rack and pinion,ball screw and hydraulic mechanism[12,29,30].Inerter devices have significant mass amplification effect,and the effective mass produced can exceed 100 times of its physical mass[13,31,32].The principle of the inerter device is that it can transform the linear motion of both ends of the device into the rotational motion of the flywheel when an earthquake occurs,resulting in considerable inerter[33–35].For example,Ikago et al.realized an apparent mass of 300 kg with a ball screw device with a physical mass of only 2 kg[36],Javidialesaadi et al.designed and produced a ball screw mechanism with a physical mass of only 2 kg to produce an effective mass of 350 kg[13].It can be seen that in order to produce an equivalent effective mass,the actual physical mass of the inerter device is far less than that of the traditional TMD.In addition,the research also showed that when the TID is placed at the bottom of the main structure,the performance of the TID is the best,which is opposite to the TMD that must be placed at the position where the maximum displacement occurs,that is,the top of the structure [22].TID can obtain a large inerter mass ratio without increasing the physical mass,so it is considered that TID can replace the traditional TMD and be applied in engineering practice.Based on inerter device,Dai et al.explored the optimal design of Maxwell tuned-mass-damper-inerter (M-TMDI) to mitigate vortex-induced vibration (VIV) of bridges.Based on a two-DOF system,the optimal parameters of a specific M-TMDI,in which the end of the inerter is connected to the fixed ground,are analytically given using the inerter location as a design variable[37].In addition,Dai et al.aimed at the impact of the location of inerter on the control performance of tuned mass-damping-inerter(TMDI)in wind-induced vibration mitigation of flexible structures,discussed the control effect,optimum design parameters,and high-mode damping effects of the TMDI[38].Different types of control devices based on inerter devices are also used in base isolated structures.TIDs are installed at the bottom of the structure to further reduce the displacement of isolated floors and the inter story displacement of superstructure[35,39,40].Among them,the inerter device based on the generator device can collect energy while reducing the displacement of the upper structure,realizing the trade-off between the two contradictory goals of energy collection and vibration control[35].Gonzalez-Buelga et al.proposed a circuit containing both damping and energy collection effects,and connected it with the electromagnetic transducer to synthesize TID,which also achieved the dual goals of damping and energy collection[41].Unlike TMD driven by the absolute motion of the damper,TID is driven by the relative motion between the two ends of the device,which actually provides conditions for the corresponding energy collection mode of TID[42].

In order to achieve the best control effect of the vibration control device,different types of optimization methods have been proposed,including the most commonly used method,namely H2norm optimization.The H2optimization criterion was first proposed by Crandall and Mark in 1963 and applied to the design of dynamic vibration absorbers (DVAs),with the goal of reducing the total vibration energy of the system at all frequencies[43].It can be understood that if the vibration control system is subjected to random excitation,that is,the excitation contains an infinite number of frequencies,then it is unnecessary to consider only the resonance frequency of the system.Under this optimization criterion,the area under the frequency response curve of the system (called H2norm) is the smallest [44].In the study of vibration control of adjacent buildings based on inerter,Palacios-Quiñonero et al.proposed a computational strategy to design inerter-based multi-actuation systems for the seismic protection of adjacent structures based on an H∞cost-function and used a constrained global-optimization solver to compute parameter configurations with high-performance characteristics.The research results fully demonstrate the flexibility and effectiveness of the proposed design methodology,and clearly show the superior performance and robustness of the TID actuation systems [45].Djerouni et al.proposed a new configuration of TMDIs with the goal of studying the seismic response of two adjacent buildings and the pounding distance between them.The results show that the proposed method is effective and better than the existing TID configuration,filling a gap in the research of investigating the effectiveness of TMDI in controlling seismic pounding and the impact of ground motion variability[46].Taking adjacent buildings as controlled structures,this paper investigated the vibration control performance of serial TID and parallel TID in adjacent buildings.Also,the analytical solutions of the optimum frequency ratio and damping ratio of structural vibration control under random excitation were obtained by using the H2norm optimization method.In addition,for the convenience of analysis,the random excitation used in this paper is white noise excitation when analyzing the frequency response of the system.

The main contents of this paper are summarized as follows:First,the analytical solutions of the optimum frequency ratio and the optimum damping ratio of adjacent building structures under the coupling of serial TID and parallel TID were obtained by using the H2optimization method.Secondly,the influence of different structural parameters on the optimum solutions was studied.Then,the classic TMD was introduced to compare with the two types of TID,and the frequency response of adjacent building structures under different damper coupling was analyzed.At the same time,the influence of damper mass was also considered.Then,by changing the optimum frequency ratio and the optimum damping ratio,the robust performance of vibration of adjacent buildings under TID coupling was analyzed.Finally,four types of seismic waves were used to analyze the time history of vibration control of adjacent building structures under two types of TIDs and classic TMD coupling,and the time history of energy collection of adjacent building structures under two types of TIDs coupling was analyzed.

2 Mathematical Model

2.1 Simplified Model of Adjacent Buildings

As shown in Fig.1a is the simplified model of adjacent buildings,wheremLi(i=1,2,...,n),cLi(i=1,2,...,n)andkLi(i=1,2,...,n)are the corresponding mass,damping and stiffness of theifloor of the left building respectively;mRi(i=1,2,...,n),cRi(i=1,2,...,n)andkRi(i=1,2,...,n)are the corresponding mass,damping and stiffness of theifloor of the right building respectively;csj(j=1,2,...,n)andksj(j=1,2,...,n)are the corresponding damping and stiffness of thejfloor of adjacent buildings.Under the excitation of seismic acceleration,in order to better analyze the vibration control effect of dampers on adjacent building structures,the left and right buildings are simplified into single degree of freedom structures,with their mass,stiffness,and damping values corresponding to the main modes of the structure[47,48].The simplified model is shown in Fig.1b.After simplification,the mass,damping and stiffness of the left building arem1,c1andk1,respectively,and the mass,damping and stiffness of the right building arem2,c2andk2,respectively.The damping and stiffness between the two buildings arecsandks,respectively.

Figure 1 :Simplified schematic graph of adjacent building structures

2.2 Inerter Device Concept

Adding ball screw mechanism(including shell and linear guide rail),flywheel and rotary generator inside the damper can transform the relative linear motion between the structure and the ground into rotation when an earthquake occurs,and drive the flywheel and generator to collect energy[35].The rotating generator and flywheel drive the inerter mass to rotate,and generate resistance to the driving structure.The damper is elastically installed with the structure,which is equivalent to a spring.Fig.2a shows a simplified model of the inerter device.It is characterized in that the generated force is proportional to the difference of acceleration at both ends [22].The mechanical equation can be expressed as:

wherexl,xrare the displacement of the left and right ends of the inerter device,the symbol “..” denotes the calculation of the second derivative,andbdenotes the effective mass.

Figure 2 : Inerter device and two types of inerter-based damper models: (a) Inerter device; (b) Serial TID;(c)Parallel TID

2.3 Damping System

According to Section 2.2,considering different connection modes,two damper models with different configurations based on inerter are formed.According to the different positions ofmdandkd,serial connection ofmdandcdis called serial TID,and parallel connection ofmdandcdis called parallel TID.The two TID models are shown in Figs.2b and 2c,wheremdis the effective mass of TID,cdandkdis the damping and stiffness of TID,respectively.Couple the different types of TIDs mentioned above with the simplified adjacent buildings,and the formed models are shown in Fig.3,wherex1,x2andxdare the displacement of left building,right building and TID under earthquake excitation,respectively.

Figure 3 :The damping system of adjacent buildings coupled with TID:(a)Serial TID damping system;(b)Parallel TID damping system

3 Parameter Optimization Based on H2 Norm Criterion

3.1 Equation Establishment

According to the coupled damping system in Fig.3,the dynamic balance equation of the system can be listed.We stipulate that the structure on the left which connected to TID is the left structure,and the structure on the right of adjacent buildings is the right structure,and they will be distinguished by left structure and right structure in subsequent expressions.

3.1.1SerialTID

According to the balance principle of force,the dynamic balance equation can be obtained as follows:

After dimensionless transformation by Laplace transform,the following results are obtained:

Define some new parameters for the next transformation.The defined parameters are shown in Table 1.After substituting them into Eq.(3),there are:

Table 1 : Parameters and their definitions during derivation

3.1.2ParallelTID

According to the balance principle of force,the dynamic balance equation can be obtained as follows:

Similar to the serial TID,it can be obtained from Table 1 and Eq.(5)as follows:

3.2 Solution of the Transfer Function

Before H2optimization,it is necessary to obtain the normalized displacement transfer functions of the left structure and the right structure under different connection modes of TID.Solve Eqs.(4)and(6)to obtain the following results.

3.2.1SerialTID

(1) Left structure

where the superscript “s” denotes serial connection,and the subscript “1” denotes the left structure,and its numerator and denominator coefficients are shown in Appendix A.

(2) Right structure

where the superscript “s” denotes serial connection,the subscript “2” denotes the right structure,the denominator coefficient is the same as the left structure,and the numerator coefficient is shown in Appendix B.

3.2.2ParallelTID

(1) Left structure

where the superscript “p” denotes parallel connection,the subscript “1” denotes the left structure,and its numerator and denominator coefficients are shown in Appendix C.

(2) Right structure

where the superscript “p” denotes parallel connection,the subscript “2” denotes the right structure,the denominator coefficient is the same as the left structure,and the numerator coefficient is shown in Appendix D.

3.3 Parameters Optimization

In order to verify the damping effect of different types of TID coupled with adjacent buildings,H2optimization method is used to optimize the parameters of vibration control system,and the external excitation used in this paper is white noise excitation.If the evaluation index is defined asPI,the H2optimization problem can be expressed as the optimization of the corresponding variance.The form of evaluation indexPIis:

To simplify the analysis,take the damping between the left structure and the right structure as 0,that is,from Table 1:ξ1=0,μ2=0,μs=0.After solving the normalized transfer function,the parameters of the left structure and the right structure can be optimized for adjacent buildings with different types of TID coupling.

3.3.1OptimizationProcess

(1) Serial TID

After substituting Eqs.(7)and(8)into the evaluation indexPIof Eq.(11),using the formula in Appendix E,the evaluation indexPIare calculated as follows:

where the superscript “s” denotes serial connection,the subscripts “1,2” denote the left structure and the right structure,respectively.Psn(n=1,2,3,4),Qsn(n=1,2,3,4),Msn(n=1,2,3,4)andNsn(n=1,2,3,4)are polynomials composed ofβd,β2,λ2andλs,respectively,see Appendix F and Appendix G for the specific form.

(2) Parallel TID

After substituting Eqs.(9)and(10)into the evaluation indexPIof Eq.(11),using the formula in Appendix E,the evaluation indexPIare calculated as follows:

where the superscript “p” denotes parallel connection,the subscripts “1,2” denote the left structure and the right structure,respectively.Pp1,Qpn(n=1,2,3,4),MPn(n=1,2,3,4,5)andNpn(n=1,2,3,4,5,6)are polynomials composed ofβd,β2,λ2andλs,respectively,see Appendix H and Appendix I for the specific form.

3.3.2OptimizationResults

In order to obtain the optimum solutions,the evaluation indexPIis used to calculate the partial derivatives of the optimization parameters respectively,and make them equal to 0.The following results are obtained:

(1) Serial TID

By solving Eqs.(16)and(17),respectively,we can get the optimization results(Eq.(18)represents the left structure,and Eq.(19)represents the right structure):

Here,it is necessary to point out:

(2) Parallel TID

By solving Eq.(21)and by Eq.(20),we can get the optimization results of the left structure under the parallel TID coupling are:

When solving Eq.(22),this equation group is converted into the solution of Eq.(24),which includes:

where

see Appendix J for the specific form ofRn(n=1,2,...,10).

By solving Eq.(24),and with Eq.(20),the optimization results of the right structure under parallel TID coupling are obtained as follows:

whereUn(n=1,2,3,4) is a polynomial composed ofRn(n=1,2,...,10).See Appendix K for the specific form.

4 Parameters Analysis

4.1 Parameters Determination

Before parameters analysis,it is necessary to give specific values of parameters that have not been optimized.In order to facilitate the optimization analysis,the mass,damping and stiffness of the left structure and the right structure are determined to be consistent in this section,that is,Table 1 shows:β2=1.0,λ2=1.0.In the following analysis,we first determine the mass ratioβdand the stiffness ratioλs,and then conduct relevant parameters analysis.

4.1.1Determinationofβd

In Section 3.3.2,the process and results of parameter optimization are given.However,many of the optimization results have root signs.According to the mathematical knowledge,all values under the root sign should be greater than or equal to 0.After considering all the optimization results,it is found that the conditions affecting the value ofβdare shown in Eq.(27):

To determine the value ofβd,let the left side of inequality(27)be 0,that is,there is Eq.(28):

Combining Eqs.(H1)and(H2),the image of Eq.(28)is drawn as shown in Fig.4.It can be seen from the figure that ifγ≥0,the value ofβdshould be between 0 and 0.244.Considering that the value ofβdis often small in real life,we can temporarily take the mass ratioβd=0.1 for the next analysis.

4.1.2Determinationofλs

It was found during parameter optimization that ifλsis the parameter to be optimized,the equation set could not be solved.So,the determination ofλsin advance is a key problem,and it is also a different point and difficulty in dealing with the vibration control of adjacent buildings.In the actual project,ksis generally less than the stiffness of the left structurek1.Therefore,it can be determined in advance that the value range ofλsis 0.1～0.7.The determination ofβdhas been given in Section 4.1.1,and the value is 0.1.According to the above conditions,we can get the function of normalized displacement,which only contains normalized frequencyαand stiffness ratioλs.The two parameters,then,can draw a three-dimensional graph for analysis,and then determine the value ofλs.Fig.5 shows a three-dimensional graph of normalized displacement under different conditions,considering different values ofλs.In Fig.5,αis between 0 and 2,andλsis between 0.1 and 0.7.

As can be seen from Fig.5,whenλstakes different values,the normalized displacement graph of the left structure coupled with two connection modes TID is similar (Figs.5a and 5c),and the normalized displacement graph of the right structure is also similar(Figs.5b and 5d).The peak value of the normalized displacement(for the convenience of description,the peak value of the displacement transfer function graph is called the wave peak)shows that for the left structure,the value ofλsdoes not affect the number of wave peaks.But for the right structure,whenλs≤0.3 or so,the number of wave peaks will degenerate into two.It can also be seen from Fig.5 that whenλs≤0.2,the normalized displacement will increase sharply,while whenλsis between 0.2 and 0.7,the displacement image tends to be stable and at a low level.In order to facilitate the next research,and considering the better damping performance of the system whenλsis large and the constraints of the actual engineering structure,it is decided to setλsas 0.5.

4.1.3VerificationofOptimizationResults

Figure 5 : Three-dimensional graphs of normalized displacement when the stiffness ratio λs take different values:(a)The left structure under serial TID coupling;(b)The right structure under serial TID coupling;(c)The left structure under parallel TID coupling;(d)The right structure under parallel TID coupling

After the values of two key parametersβdandλsare determined,they can be substituted into the optimization result equations under different conditions,and the values of the optimum solutionsλdoptandξdoptcan be obtained.At the same time,according to the evaluation indexPI(Eqs.(12)–(15))given earlier,a three-dimensional graph ofPI(Fig.6)can be drawn,and the accuracy of the optimization results can be verified.

It can be seen from Figs.6a and 6b that the evaluation indexesPIof the left structure and the right structure under the serial TID coupling are similar,the lowest point appears nearλd=0.09～0.11 andξd=0.10～0.14.For the three-dimensional graph of evaluation indexPIof left structure and right structure under the coupling of parallel TID,it can be seen from Figs.6c and 6d that whenλd≥0.28 andξd≥0.3,the value ofin three-dimensional graph is relatively stable and the minimum value ofPIcannot be obtained within this range,but whenλd=0.10～0.11 andξd=0.11～0.15,can get the minimum value.Similarly,whenλd=0.11～0.13 andξd=0.10～0.13,PIp2can get the minimum value.In order to further determine the values ofλdandξd,the contour graphs are drawn in detail as shown in Fig.7.

Figure 6 :Three-dimensional graphs of evaluation index PI by H2 optimization method(λd is the ratio of stiffness of TID to stiffness of left structure,ξd is the nominal damping ratio of TID):(a)The left structure under serial TID coupling; (b) The right structure under serial TID coupling; (c) The left structure under parallel TID coupling;(d)The right structure under parallel TID coupling

From the closure of the contour graphs in Figs.6 and 7,we can roughly judge the coordinates whenPI=0,and the results are listed in Table 2.As mentioned above,the calculated value of the optimum solutions can be obtained by substituting the previously determined parameters:β2=1.0,λ2=1.0,λs=0.5 andβd=0.1 into Eqs.(18),(19),(23) and (26),and the results are shown in Table 2.As can be seen from Table 2,compared with the calculated value of the optimum solutions,the error of the predicted optimum solutions from Fig.7 is small,and the maximum error is only 3.31%.Therefore,the optimization results obtained by the process in Section 3 are accurate.Due to the error in the estimated value obtained through figure judgment,the calculated values are used for the optimum solution during subsequent analysis.

Table 2 :Comparison between the estimated values and the calculated values of the optimum solutions

Figure 7 : The contour graphs of evaluation index PI of H2 optimization method (a,b,c and d correspond to Fig.6).The blue pentagram marks the point where PI=0,and its corresponding λd and ξd are the optimum solutions λdopt and ξdopt

4.2 Influence of Parameters on Optimum Solutions

In 1998,Rana et al.put forward that the frequency ratio is an important index that affects the damping performance when designing TMD through the fixed-points theory [49].Considering the impact of actual engineering on the optimum solution results is conducive to further understanding the laws of TID for vibration control of adjacent building structures.The influence of different parameters on the optimum solution is analyzed as follows:

4.2.1InfluenceofMassRatiosβdandβ2

Figs.8 and 9 are three-dimensional graphs of optimum solutionsfdoptandξdoptof left structure and right structure under the coupling of serial TID and parallel TID,respectively.Figs.10 and 11 are two-dimensional graphs corresponding to Figs.8 and 9,respectively.

Figure 8 :Three-dimensional graphs of optimum solutions fdopt and ξdopt under serial TID coupling(λ2=1.0,λs=0.5).In the figure,(a1)and(a2)correspond to the left structure,and(b1)and(b2)correspond to the right structure

Figure 9 : Three-dimensional graphs of optimum solutions fdopt and ξdopt under parallel TID coupling(λ2=1.0,λs=0.5).In the figure,(c1) and (c2) correspond to the left structure,and (d1) and (d2)correspond to the right structure

In Figs.8(a1) and 8(a2),the optimum frequency ratiofdoptof the left structure under serial TID coupling decreases with the increase of mass ratioβd,while the optimum damping ratioξdoptincreases with the increase of mass ratioβd.It can be seen from Fig.10(a1)that whenβ2is small,fdoptdecreases rapidly,and with the increase ofβ2,fdoptdecreases more and more slowly.In Fig.10(a2),the value ofβ2has little effect on the rising speed ofξdopt.In Figs.8(b1),and 8(b2),for the right structure coupled by serial TID,β2has obvious influence onfdoptandξdopt.Whenβ2is small,bothfdoptandξdoptdecrease with the increase of mass ratioβd;While whenβ2is larger,bothfdoptandξdoptincrease with the increase ofβd.The above results show that the mass ratioβ2has a great influence on the right structure in the adjacent building structures coupled by serial TID.In order to achieve the best control effect,with the increase ofβd,the required spring stiffnesskddecreases and the dampingcdincreases for the left structure.However,for the right structure,when the mass ratioβ2is small,with the increase ofβ2,the required spring stiffnesskdand dampingcdbecome smaller and smaller.When the mass ratioβ2is large,with the increase ofβ2,the required spring stiffnesskdand dampingcdbecome larger and larger.

Figure 10 :The curves of the optimum solutions(fdopt,ξdopt)of the left structure and the right structure under the coupling of serial TID when βd takes different values(a and b correspond to Figs.8a and 8b)

Under parallel TID coupling,it can be seen from Fig.9 that thefdoptandξdoptof the left structure and the right structure both increase with the increase ofβd.According to the observation in Fig.11,whenβ2is small,fdoptandξdoptincrease rapidly,while whenβ2is large,the increasing trend offdoptandξdoptslows down.Therefore,the mass ratioβ2has obvious influence on both the left structure and the right structure,and the influence trends are the same.For the left structure and the right structure,in order to achieve the best control effect,with the increase ofβd,the required spring stiffnesskdand dampingcdbecome larger and larger,and the largerβ2is,the slower the increase speed is.

4.2.2InfluenceofStiffnessRatiosλsandλ2

Figs.12 and 13 are three-dimensional graphs of optimum solutionsfdoptandξdoptof left structure and right structure under the coupling of serial TID and parallel TID,respectively.Figs.14 and 15 are two-dimensional graphs corresponding to Figs.12 and 13,respectively.

Figure 11 :The curves of the optimum solutions(fdopt,ξdopt)of the left structure and the right structure under the coupling of parallel TID when βd takes different values (c and d correspond to Figs.9c and 9d)

Figure 12 : (Continued)

Figure 12 :Three-dimensional graphs of the optimum solutions fdopt and ξdopt under serial TID coupling(βd=0.1,β2=1.0).In the figure,(a1) and (a2) correspond to the left structure,and (b1) and (b2)correspond to the right structure

Figure 13 : Three-dimensional graphs of the optimum solutions fdopt and ξdopt under parallel TID coupling(βd=0.1,β2=1.0).In the figure,(c1)and(c2)correspond to the left structure,and(d1)and(d2)correspond to the right structure

Figure 14 :The curves of the optimum solutions(fdopt,ξdopt)of the left structure and the right structure under the coupling of serial TID when λs takes different values (a and b correspond to Figs.12a and 12b)

Figure 15 :The curves of the optimum solutions(fdopt,ξdopt)of the left structure and the right structure under the coupling of parallel TID when λs takes different values (c and d correspond to Figs.13c and 13d)

For the adjacent buildings coupled by serial TID,it can be seen from Figs.12,14(a1)and 14(a2)that the optimum frequency ratiofdoptof the left structure increases with the increase of stiffness ratioλ2,and it is not obviously affected by stiffness ratioλs.The optimum damping ratioξdoptis obviously affected by the stiffness ratioλs.Whenλsis small,ξdoptfirst decreases and then increases with the increase ofλ2,while whenλsis large,ξdoptalways increases with the increase ofλ2.The above results show that,in order to achieve the best control effect,with the increase of stiffness ratioλsandλ2,the required stiffness provided by TID will be larger for the left structure coupled with serial TID.As for damping,whenλsis small,with the increase ofλ2,the damping provided by TID decreases first and then increases;whenλsis large,with the increase ofλs,the damping required by TID is larger.It can be seen from Figs.12,14(b1)and 14(b2)that the optimum frequency ratiofdoptand the best damping ratioξdoptof the right structure are significantly affected by the stiffness ratioλs.Whenλsis small(λs≤0.3),fdoptincreases first and then decreases and then increases with the increase of stiffness ratioλ2,andξdoptincreases first and then decreases with the increase of stiffness ratioλ2.Whenλsis larger(0.3＜λs≤0.7),bothfdoptandξdoptincrease with the increase of stiffness ratioλ2,and the increasing speed is relatively slow.The above results show that,in order to achieve the best control effect,whenλsis small,with the increase ofλ2,the stiffness provided by TID increases first and then decreases and then increases,and the required damping increases first and then decreases.Whenλsis larger,the largerλ2is,the greater the stiffness and damping that TID needs to provide.

For adjacent buildings with parallel TID coupling,it can be seen from Figs.13,15(c1)and 15(c2)that the optimum frequency ratiofdoptof the left structure is obviously affected by the stiffness ratioλs.Whenλsis small,fdoptfirst decreases and then increases with the increase ofλ2,while whenλsis larger,fdoptalways increases with the increase ofλ2.ξdoptis not obviously affected by the stiffness ratioλs.Under the same stiffness ratioλs,ξdoptdecreases with the increase ofλ2.The above results show that,in order to achieve the best control effect,whenλsis small,with the increase ofλ2,the stiffness provided by TID decreases first and then increases,and the required damping also decreases first and then increases; whenλsis larger,with the increase ofλ2,the stiffness required by TID is larger,and the required damping also decreases first and then increases.For the right structure,it can be seen from Figs.13,15(d1) and 15(d2) that the optimum frequency ratiofdoptand the optimum damping ratioξdoptare significantly affected by the stiffness ratioλs,similar to the serial TID.Whenλsis small(λs≤0.3),fdoptincreases first and then decreases and then increases with the increase of stiffness ratioλ2,andξdoptincreases first and then decreases with the increase of stiffness ratioλ2.Whenλsis larger(0.3＜λs≤0.7),bothfdoptandξdoptincrease with the increase of stiffness ratioλ2,and the increasing speed is relatively slow.This shows that,in order to achieve the best control effect,with the change of stiffness ratioλsandλ2,the change trend of required stiffness and damping of the right structure coupled with parallel TID is similar to that of the right structure coupled with serial TID.

5 Frequency Response Analysis

After parameter analysis,we introduce classic tuned mass damper (TMD) into the earthquake resistance of adjacent building structures,and analyze its vibration control performance.Under the condition that all the conditions of the specified structure are unchanged,the vibration control performance of two types of TID and TMD are compared.

5.1 Classic TMD System

5.1.1ClassicTMDEquations

Fig.16 shows the damping model of adjacent buildings coupled with TMD,which is the same as the definition in Section 3.1.It is defined that the structure on the left side is the left structure and the structure on the right side is the right structure.The mass,damping and stiffness of TMD aremt,ctandkt,respectively.m1,m2are the mass of the left structure and the right structure,respectively,c1,c2andcsare the damping of the left structure,the right structure and the damping between the left and right structures,respectively,k1,k2andksare the stiffness of the left structure,the right structure and the stiffness between the left and right structures,respectively.is the acceleration of seismic waves,andx1,x2andxtare the displacement of the left structure,right structure and TMD,respectively.

Figure 16 :The damping system of adjacent buildings coupled with classic TMD

According to the principle of force balance,the dynamic equilibrium equations are:

Using Laplace transform to make dimensionless,and after finishing,get:

whereβt=mt/m1is the ratio of the classic TMD mass to the left structure mass,ξt=is the nominal damping ratio of classic TMD,λt=kt/k1is the ratio of classic TMD stiffness to left structural stiffness,ft=ωt/ω1is the ratio of the nominal frequency of the classic TMD to the frequency of the left structure,Xtis the form ofxtafter Laplace transformation.The definitions of other parameters are the same as those in Section 3.1.1,see Table 1.

5.1.2OptimizationResults

Same as the steps in Section 3,the optimization results obtained by H2optimization method are as follows:

(1) Left structure

wherePn(n=1,2,3,4)andQn(n=1,2,3,4)are polynomials composed ofβt,β2,λ2andλs,respectively,see Appendix L for the specific form.

Similarly,for the classic TMD,there are:

(2) Right structure

whereMn(n=1,2,3,4)andNn(n=1,2,3,4)are polynomials composed ofβt,β2,λ2andλs,respectively,see Appendix M for the specific form.

In order to compare the damping performance of two types of TID and classic TMD,this section adopts the parameters determined in Section 4.1,namelyβ2=1.0,λ2=1.0,λs=0.5.When analyzing the parameters of TID,we take the mass ratioβd=0.1,so when analyzing the damping effect of adjacent buildings under the classic TMD coupling,we also take the mass ratioβt=0.1,and then make the next comparative analysis.

5.2 Comparison of Vibration Control Effect

In the previous analysis,the values of relevant parameters have been determined as follows:β2=1.0,λ2=1.0,λs=0.5,theβdof adjacent buildings under TID coupling is 0.1,and theβtof adjacent buildings under classic TMD coupling is also 0.1.The normalized displacement figure under different conditions is shown in Fig.17.

As can be seen from Fig.17,the similarity of normalized displacement curves under different conditions is high,and there are three peaks.According to the observation of Fig.17a,compared with the highest point of the whole curve,the parallel TID value is the smallest; according to the observation of Fig.17b,compared with the highest point of the whole curve,the classic TMD value is the smallest.The above results show that when damping adjacent buildings,the best control effect can be achieved by selecting parallel TID for the left structure and classic TMD for the right structure.For further comparison,the relevant data of Figs.17a and 17b are summarized in Table 3.It can be seen from the table that compared with the highest point of displacement curve,the left structure is smaller than the right structure for the adjacent buildings coupled by serial and parallel TIDs,which indicates that TID has better vibration control effect for the left structure.However,for the adjacent buildings coupled with classic TMD,the highest point of the displacement curve of the left structure is larger than that of the right structure,which indicates that the classic TMD has a better damping effect on the right structure.From the results given in Table 3,it can be seen that if the size of the areaS(i.e.,H2norm) enclosed by the displacement curve and the coordinate axis in the normalized frequency range of 0～2.0 is taken as the basis(hereinafter referred to as areaSfor convenience of analysis),the serial TID has a better vibration control effect for the left structure,and the classic TMD has a better vibration control effect for the right structure.

Table 3 : Normalized displacement correlation coefficients of adjacent buildings under different damper coupling

Figure 17 : Normalized displacement graphs of the left structure and right structure of adjacent buildings with different types of dampers(β2=1.0,λ2=1.0,λs=0.5,βd=0.1,βt=0.1)

As far as the damping effect of TID is concerned,only the specific case of mass ratioβd=0.1 has been analyzed above,which is obviously not well considered.Fig.18 is a three-dimensional graph of the normalized displacement of the left structure and the right structure under the coupling of serial TID and parallel TID.From this graph,it can be found that the normalized displacement will gradually decrease with the increase of the mass ratioβd.For better comparison,the peak value and areaSof the displacement curve whenβdtakes different values are calculated and summarized in Fig.19.

Figure 18 : Three-dimensional graphs of normalized displacement when the mass ratio βd takes different values:(a)The left structure under serial TID coupling;(b)The right structure under serial TID coupling;(c)The left structure under parallel TID coupling;(d)The right structure under parallel TID coupling(β2=1.0,λ2=1.0,λs=0.5)

It can be seen from Fig.19a that increasing the mass ratioβdcan obviously enhance the damping effect.From the peak value of displacement,whenβd≤0.1,increasingβdwill obviously reduce the peak value of displacement.However,whenβd＞0.1,the peak value of displacement will decrease more and more slowly.By observing the change of the areaSin the bar graph,the height difference between the two bar graphs is also slowly decreasing.According to the line graph and bar graph,it can be found that parallel TID has obvious advantages in controlling the peak displacement,and this advantage becomes greater with the increase ofβd.Comparing the size of areaS,it can be seen that serial TID has a significant advantage in controlling the overall vibration control effect,and this advantage increases with the increase ofβd.It can be seen from the line graph and curve graph in Fig.19b that the value ofβdhas the same effect on the peak displacement and areaSof the right structure and the left structure.At the same time,for the right structure,the serial TID has a significant advantage in controlling the overall vibration control effect,and this advantage increases with the increase ofβd.However,compared with the peak value of the displacement of the right structure,the serial TID and the parallel TID are similar in size,and the influence ofβdis not obvious.

Figure 19 : Considering the different values of mass ratio βd,the related calculation results of normalized displacement graphs of left structure and right structure under the coupling of serial TID and parallel TID(In the figure,the line graph represents the maximum value of the normalized displacement,and the column graph represents the area S enclosed by the displacement curve and the coordinate axis,that is,the H2 norm)

The above analysis results show that when the mass of the damper is small compared with that of the left structure(βd≤0.1),increasing the mass of the damper can really enhance the damping effect.However,when the mass of the damper increases gradually compared with that of the left structure(βd＞0.1),the constraints of actual cost and site conditions should be considered,and the damping effect of the damper with large mass compared with the left structure is not obviously improved.In addition,when controlling the peak displacement of adjacent buildings under earthquake,parallel TID should be given priority.If H2norm is taken into account,it is better to choose serial TID.

5.3 Robustness Performance Analysis

Optimization parameters in practical engineering often change.In order to make adjacent buildings achieve more robust vibration control effect,it is necessary to analyze their robustness.Fig.20 is the robustness analysis curve of the left structure and the right structure under the coupling of serial TID and parallel TID.The optimum frequency ratiofdoptand the optimum damping ratioξdoptare the influencing factors.

As shown in Fig.20,slightly change the damping ratio the value ofξdopthas little influence on the vibration control effect of adjacent buildings under serial and parallel TID coupling,and only causes slight changes at the peak position,which indicates that the normalized displacement is not sensitive to the change of damping ratio.However,if the frequency ratiofdoptis slightly changed,it can be seen from Fig.20 that the normalized displacement change is significant.In addition,it can be seen from the Fig.20d that there is a special case,that is,when the frequency ratiofdoptof the right structure under parallel TID coupling decreases,the amplitude of change is so large that the change of wave peak height cannot be accurately seen.

Figure 20 : Robustness performance analysis curves of left structure and right structure of adjacent buildings under the coupling of serial TID and parallel TID: (a) Serial TID is coupled with the left structure;(b)Serial TID is coupled with the right structure;(c)Parallel TID is coupled with the left structure;(d)Parallel TID is coupled with the right structure(β2=1.0,βd=0.1,λ2=1.0,λs=0.5)

To further determine the size of the change,the detailed drawing of Fig.20d is shown in Fig.21.It can be seen from the figure that the right structure in adjacent buildings under the coupling of parallel TIDs is extremely sensitive to the reduction of frequency ratiofdopt,and the displacement after the change is very large,and the vibration control effect decreases obviously.The above results show that,in the actual engineering structure design,in order to achieve more robust vibration control effect of adjacent building structures,the change of stiffness is more worthy of attention than the change of damping of TID.

Figure 21 : Robust performance analysis curve of the right structure of adjacent buildings under the coupling of parallel TID

6 Time Domain Simulation

6.1 Vibration Control Performance

In the foregoing analysis,we assume that the external excitation is white noise excitation with evenly distributed energy at all frequencies,but the ground motion caused by earthquake is actually not ideal white noise [35].Therefore,it is necessary to simulate and verify the seismic performance of adjacent buildings under the coupling of serial TID and parallel TID.Taking two ten-story frame structures as an example,the relevant parameters of the adjacent building structures are:m1=m2=5.897×104kg,c1=c2=67.4×103N·s/m,cs=33.7×103N·s/m,k1=k2=4.81×107N/m,ks=2.405×107N/m[50],βd=0.1,βt=0.1 and relevant parameters of dampers are determined in Table 4.The external excitation is four commonly used seismic waves,namely:El Centro(May 18,1940),Taft(July 21,1952),Loma Prieta(October 17,1989)and Chi-Chi(September 21,1999)[2,35].Fig.22 shows the acceleration time history of four types of seismic waves.Under the excitation of four types of seismic waves,the displacement time history curves of adjacent buildings coupled with different dampers are shown in Figs.23 and 24.

Table 4 : Relevant parameters of damper

Figure 22 :Acceleration time history curves of four types of seismic waves

Figure 23 : Displacement time history graphs of the left structure under four types of seismic waves excitation

It can be seen from Fig.23 that under the excitation of different types of seismic waves,compared with the displacement of the left structure in the uncontrolled state,the three dampers have good damping effects in the whole range.Comparing the displacement of three dampers,it is found that the displacement of classical TMD is slightly larger than that of serial or parallel TID.In order to make a more accurate comparison,the peak and root mean square values of the displacement of the left structure coupled with three different types of dampers are summarized,as shown in Table 5.It can be seen from Table 5 that,compared with the peak displacement in uncontrolled state,the peak control effect of the serial TID is at least 10.1% and at most 47.8%;the peak control effect of parallel TID is at least 9.5% and at most 45.9%;the peak control effect of classic TMD is at least 2.8% and at most 42.3%.As for the control effect of RMS,the control effect of serial TID is at least 21.6% and at most 53.2%;the control effect of parallel TID is at least 18.6% and at most 51.3%,the control effect of classical TMD is at least 14.7% and at most 49.4%.It can be found that under the vibration control of the above model,the root mean square control effect is often greater than the peak damping ratio.The above results show that for the left structure coupled with three types of dampers,all of them can achieve good damping effect,among which the serial TID has the best damping effect,the parallel TID is the second,and the classic TMD is the last.

Table 5 :Peak values and root mean square values of displacement time history under different seismic waves excitation

It can be seen from Fig.24 that,similar to the left structure,under the excitation of different types of seismic waves,compared with the displacement of the right structure in the uncontrolled state,the three dampers have good damping effects in the whole range.Comparing the displacement of three dampers,it is found that the displacement of classic TMD is slightly larger than that of serial or parallel TID.It can be seen from Table 5 that,for the right structure,compared with the peak displacement in uncontrolled state,the peak control effect of the serial TID is at least 12.6% and at most 51.1%;the peak control effect of parallel TID is at least 11.4% and at most 50.7%;the peak control effect of classic TMD is at least 7.9% and at most 48.0%.As for the control effect of RMS,the control effect of serial TID is at least 22.6% and at most 51.9%;the control effect of parallel TID is at least 21.6% and at most 50.6%,the control effect of classical TMD is at least 20.6% and at most 51.3%.The above results show that for the right structure coupled with three types of dampers,all of them can achieve good damping effect,among which the serial TID has the best damping effect,the parallel TID is the second,and the classic TMD is the last.In addition,it can be found that the right structure can achieve better damping effect than the left structure when the adjacent buildings are coupled with the same damper.

6.2 Energy Collection

As mentioned in Section 2.2,installing the roller generator inside the inerter device can collect part of the vibration energy during the earthquake and make use of it.Therefore,in this paper,we connect the left structure with the TID to collect the vibration energy of the building structure.Based on the fact that the power consumption of the external resistor is basically equivalent to that of the TID buffer element,the power calculation formula can be expressed asP=c1˙x21.Fig.25 shows the energy collection graph of the left structure under the excitation of four types of seismic waves,coupled with serial TID and parallel TID.It can be seen from Fig.25 that the energy power value of the left structure under the excitation of El Centro seismic waves is larger than that of the other three types of seismic waves,showing the best energy collection effect,and the energy collection effects of the other three types of seismic waves are similar.For more accurate comparison,the peak and average values of power in Fig.25 are summarized in Table 6.As can be seen from Table 6,under the excitation of the same type of seismic wave,in most cases,the peak value and average value of the collected power of the left structure with parallel TID coupling are higher than those of the left structure with serial TID,this shows that,on the premise of similar vibration control effect,the parallel TID should be preferred for the energy collection of the left structure of adjacent buildings.

Table 6 : Energy collection under different seismic waves excitation

Figure 25 :Energy collection graphs of the left structure under four types of seismic waves excitation

It can be seen from Table 6 that the left structure under the excitation of El Centro seismic wave has a remarkable energy collection effect,with its power peak value reaching 2525.56Wand power root mean square value reaching 44.66W.According to the calculation,the energy collected in the whole period of El Centro seismic wave is 1.44×105J.Compared with the classical TMD,it is known from the analysis results in Section 5.2 that the vibration reduction effect of TID is better than that of the classic TMD.At the same time,TID has the advantage of collecting energy,when an earthquake occurs,it can store part of the energy for the structure itself,and use it according to the actual situation.Therefore,in real life,when the economic cost is similar,compared with the classic TMD,TID should be preferred to achieve the dual effects of structural vibration control and energy collection.

7 Conclusions

In this paper,two types of tuned inerter damper (TID) adjacent building damping systems composed of springs,inerters and dampers in serial or in parallel were proposed.H2norm criterion was adopted to optimize and adjust the damping system of adjacent buildings,so that the system had the best damping performance under the random excitation of white noise.The parameter analysis,frequency response analysis,robustness analysis,time history response analysis and energy collection analysis of the system were carried out successively.Through the above analysis,the following conclusions are drawn:

(1)As for the H2norm optimization results of adjacent buildings,when considering the value of mass ratiosβ2andβd,only the change trend offdoptandξdoptof the right structure under the coupling of serial TID was fluctuating,and the rest changed smoothly.When considering the values of stiffness ratioλ2andλs,the change trend offdoptandξdoptof the right structure under the coupling of serial TID and parallel TID were fluctuating,while the rest changed smoothly.For the left structure,theξdoptof serial TID and thefdoptof parallel TID were fluctuating,and the rest of the situation changed smoothly.

(2)When studying the displacement frequency response of adjacent buildings,it was found that whenβd=0.1,based on the peak displacement,compared with the classic TMD,the parallel TID had a better damping effect on the left structure,but the improvement effect of TID on the right structure was not obvious.Based on the H2norm in the normalized frequency range of 0～2.0,compared with the classic TMD,the serial TID had a better control effect on the left structure,but the improvement effect of TID on the right structure was not obvious.Changing the mass ratio,whenβd≤0.1,increasing the damper mass could obviously enhance the damping effect,but whenβd＞0.1,the enhancement effect was not obvious.In addition,parallel TID had obvious advantages in controlling the peak displacement of the left structure,and the largerβdis,the more obvious it is,but it had no obvious superiority for the right structure.Serial TID had obvious advantages in controlling H2norm of the left structure and the right structure,and the largerβdis,the more obvious it is.

(3) As for the seismic robustness of adjacent building structures under TID coupling,the normalized displacement was less affected by the optimum damping ratioξdopt,but significantly affected by the optimum frequency ratiofdopt,especially for the right structure under parallel TID coupling.In actual engineering,when facing the vibration problem of adjacent buildings under white noise excitation,if we want to achieve more robust control,the first thing we need to pay attention to was the stiffness of the damper,followed by the damping.

(4)The time history analysis results of adjacent buildings excited by four types of seismic waves showed that all three types of dampers could achieve a good damping effect,among which the serial TID was the best,the parallel TID was the second,and the classic TMD was the last.In addition,for adjacent buildings coupled with the same damper,the right structure could often achieve a better damping effect than the left structure.According to the time history analysis of energy collection of adjacent buildings excited by four types of seismic waves,it was found that the left structure excited by El Centro seismic waves had the best energy collection effect,with its instantaneous peak power reaching 2525.56Wand the root mean square power reaching 44.66Win the whole period.According to the calculation,the energy collected in the whole period of El Centro seismic wave was 1.44×105J.In addition,when collecting the vibration energy of adjacent building structures under earthquake excitation,it was found that parallel TID should be given priority.

Based on the above analysis results,in general,the proposed TIDs with different connection modes had significant effects on the vibration control of adjacent buildings,and their performance was superior to that of classic TMDs.In addition,this type of TID also had the function of energy collection,which further expands its application range.The research results of this article could provide some reference significance for the vibration control research of adjacent buildings,however,the application of TIDs proposed in this paper in practical engineering needs further research and demonstration to achieve better vibration control effects.

Acknowledgement:The authors are very grateful to the editors and all anonymous reviewers for their insightful comments.

Funding Statement: This research was funded by the Natural Science Research Project of Higher Education Institutions in Anhui Province(Grant No.2022AH040045),the Anhui Provincial Natural Science Foundation (Grant No.2008085QE245),the Project of Science and Technology Plan of Department of Housing and Urban-Rural Development of Anhui Province(Grant No.2021-YF22).

Author Contributions:The authors confirm contribution to the paper as follows: study conception and design: Xiaofang Kang and Jian Wu; data collection: Jian Wu and Xinqi Wang; analysis and interpretation of results:Jian Wu and Shancheng Lei;draft manuscript preparation:Xiaofang Kang and Jian Wu.All authors reviewed the results and approved the final version of the manuscript.

Availability of Data and Materials:No new data and materials were created or analyzed in this study.Data and material sharing are not applicable for this article.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

Appendix A.The numerator and denominator coefficients of Eq.(7)are as follows:

Appendix B.The numerator coefficient of Eq.(8)is as follows:

Appendix C.The numerator and denominator coefficients of Eq.(9)are as follows:

Appendix D.The numerator coefficient of Eq.(10)is as follows:

Appendix E.Using the Residue Theorem,the evaluation indexPIcan be calculated by the following formula:

where

Appendix F.In Eq.(12),the specific forms ofPsn(n=1,2,3,4)andQsn(n=1,2,3,4)are:

Appendix G.In Eq.(13),the specific forms ofMns(n=1,2,3,4)andNns(n=1,2,3,4)are:

Appendix H.In Eq.(14),the specific forms ofPp1,Qpn(n=1,2,3,4)are:

Appendix I.In Eq.(15),the specific forms ofMnP(n=1,2,3,4,5)andNnp(n=1,2,3,4,5,6)are:

Appendix J.The specific forms ofRn(n=1,2,...,10)are:

Appendix K.The optimization results Eq.(26)can be obtained by solving Eq.(24),where the specific form of polynomialRn(n=1,2,...,10)are as follows:

Here,it is necessary to point out the calculation process of the optimum solutionsfdoptandξdoptunder parallel TID coupling in Table 4.

Substituting Eqs.(I1),(I2),(20),(J1)–(J10)and(K1)–(K5)into Eq.(26),the calculated result is:

In the above Equation,iis the imaginary unit,omit the smallest imaginary part and keep three decimal places,and get the optimum solutions in Table 4.

Appendix L.In Eq.(31),the specific forms ofPn(n=1,2,3,4)andQn(n=1,2,3,4)are:

Appendix M.In Eq.(33),the specific forms ofMn(n=1,2,3,4)andNn(n=1,2,3,4)are: