Stochastic Response Characteristic and Equivalent Damping of Weak Nonlinear Energy Dissipation System under Biaxial Earthquake Action

The random response characteristic of weak nonlinear structure under biaxial earthquake excitation is investigated. The structure has a SDOF (single degree of freedom)with supporting braces and viscoelastic dampers. First, it adopts integral constitutive relation and establishes a differential and integral equations of motion. Then, according to the principle of energy balance, the equation is linearized. Finally, based on the stochastic averagingmethod, the general analytical solution of the variance of the displacement and velocity response and the equivalent damping is deduced and derived. At the same time, the joint probability density function of the amplitude and phase and displacement and velocity of the energy dissipation structure are also given. The dynamic characteristics of a structure with viscoelastic dampers are determined as a solution to the variance of displacement response, so the equivalent damping is taken into consideration as a solution to replace the original nonlinear damping. It means it has established a unified analytical solution of stochastic response analysis and equivalent damping of a SDOF nonlinear dissipation structure with the brace under biaxial earthquake action in this paper.


Introduction
In addition to the two seismic force components in the horizontal direction, there is still vertical seismic force components.Actual earthquake structure is always subject to vertical and horizontal earthquake actions.Under a larger action of earthquake, the response of structure is further increased.At this time, the vertical earthquake action can not be ignored.Therefore, it is important to research structure response in the horizontal and vertical earthquake.As land in big cities is limited, buildings are located close to each other.To reduce the seismic responses of buildings, adjacent buildings are linked together by connecting dampers, such as the Triple Towers in Downtown Tokyo [1].Researchers have proposed different types of connecting devices to connect adjacent buildings.These devices include passive dampers [2][3][4][5], semiactive dampers [6][7][8], and active dampers [9,10].It is now well recognized that seismic responses of adjacent buildings can be mitigated by connecting dampers.Biaxial earthquake action will aggravate the vibration of the structure; the devices of passive control will reduce the structural vibration.The passive control techniques, such as viscous and viscoelastic dampers, have been widely used [11].Linear viscoelastic damper is a kind of excellent performance of energy dissipation device and is widely used in seismic engineering.The integral model is the most general model of viscoelastic dampers [12].Other models, such as the complex modulus model [13], the fractional derivative model [14][15][16], and the general differential model, are all approximate model.Analytical modeling of a novel type of passive friction damper for seismic hazard mitigation of structural systems is present [17]; numerical results show that the proposed damper is more efficient in dissipating input seismic energy than a passive linear viscous damper with same force capacity.The equivalent linearization of the motion equations with Maxwell dampers will effectively solve the problem of nonlinear equations.A system with nonlinear dampers is usually replaced by an equivalent linear system, with its properties determined by using different methods, like equating the energy dissipated [18], equating power consumption [19], replacing the nonlinear viscous damping by an array of frequency and amplitude-dependent linear viscous models [20], and other random vibration theories [21].Malone and Connor [22] have reported a method setting a new degree of freedom at mass-less point between a dashpot and a stiffness spring of the Maxwell model and then apply a common numerical integration scheme.In this method, it is necessary to consider twice the degrees of freedom of the original system.It has been mainly applied to the analyses of material stress-strain relationships.Kitagawa et al. [23] have reported the analysis of reinforced concrete elements by considering the effect of strain speed.They treated the Maxwell model as a supplementary restoring force on the equation of motion of the system discretized by a central difference method, which is categorized into an explicit integration scheme.It can play a better role of shock absorption by adding brace to the viscoelastic damper.The brace is widely used in damper.The integral model is a typical type; this kind of damper [12] can be used to describe the instantaneous elasticity, creep, relaxation, and strain memory of viscoelastic dampers.Park et al. [24] and Singh et al. [25] describe the use of gradientbased optimization algorithms to obtain the optimal parameters of dampers and their supporting braces in structures subjected to seismic motions.More recently, Chen and Chai [26] also proposed a gradient-based numerical procedure for determining the minimum brace stiffness together with a set of optimal damper coefficients to meet a target response reduction.They used Maxwell model-based brace-damper systems and concluded that brace stiffness equal to the first storey stiffness would be adequate for the desirable levels of response reduction in typical applications.Since the structure is installed with the damper and then turned into the energy dissipation structure, its design can not be directly applied to the response spectrum method.At the same time, this makes the design of actual engineering very inconvenient.The damping ratio of the dissipation system is the sum of the damping ratio of the structure itself and the equivalent damping ratio of the damper.The linear response spectrum method can be used to calculate the equivalent damping ratio of the damper [27,28].Therefore, it is greatly significant to establish a equivalent structure.Then, the response spectrum method can be used directly used to structural analysis and engineering design.The relationships between equivalent damping and ductility for the direct displacement-based seismic design (DBSD) method are proposed [29].As the concept of the DBSD is addressed to highlight the importance of the proper determination of equivalent damping, in the DBSD, the equivalent stiffness is taken as the secant stiffness at maximum deformation, so the appropriate equivalent damping should be determined based on such a prescription.And twenty-one SDOF systems are designed according to the DBSD procedure and analyzed to indicate that the proposed equivalent damping relationships are suitable for the DBSD.In addition, stochastic averaging method is an effective approximation method for predicting the stochastic response of a structure.The basic assumption is small damping and weak broadband excitation.Compared with the modal strain energy method, it is easy to understand and obtain the general analytical solution under a close theoretical basis; the same result of decoupling method of the forced vibration mode under the case of linear small damping can be concluded.In fact, in recent years, the important theoretical results of linear and nonlinear random vibration are obtained by using the stochastic averaging method.It investigates the stochastic response of vibroimpact system with fractional derivative under Gaussian white noise excitation; the nonsmooth transformation and stochastic averaging method are used to obtain the analytical solutions of the equivalent stochastic system [30].The first-passage statistics of Duffing-Rayleigh-Mathieu system under wide-band colored noise excitation are studied by using stochastic averaging method.The motion equation of the original system is transformed into two time homogeneous diffusion Markovian processes of amplitude and phase after stochastic averaging [31].The equivalent linearization can solve the problem of nonlinear structure; a nonlinear stochastic optimal control strategy for single degree of freedom viscoelastic system with actuator saturation is proposed based on the stochastic averaging method and stochastic dynamical programing principle.As the viscoelastic system is converted into an equivalent nonlinear nonviscoelastic system by replacing the viscoelastic force with amplitude-dependent stiffness and damping [32], in this paper, it has used the equivalent linearization method and the stochastic averaging method, and it has also used the general integral model of viscous and viscoelastic dampers.Considering the comprehensive effect of brace, strain history of damper, dynamic characteristics of structure, and excitation, it establishes a complete analytical solution of stochastic response analysis and equivalent damping of a SDOF nonlinear dissipation structure with the brace under biaxial earthquake action.The new approach can be directly applied to damping engineering design with the response spectrum method.
The motion equation can be expressed as follows: where  is the mass,  is the damping,  is the stiffness, and {  } and { V } are horizontal and vertical inertial force vector.
The original calculation diagram of damper  1 () is the viscoelastic dampers force.Relevant parameters are listed as follows: . (2)

The Transfer of the Weak Nonlinear System Equation.
Considering the weak nonlinear SDOF system, the general energy dissipation structural equation can be expressed as follows (see [33,34]): where  is the mass,  is the damping,  is the stiffness, (, u ) is the weak nonlinear force including the nonlinear damping and the spring forces, ( 1  + ∫  0 ℎ 1 ( − ) u ()) is the modified damper with supporting forces, and −( ü  () + ü V ()) is a biaxial excitation.The main aim is to replace (3) with an equivalent linear one (see [35]).

𝑚 ü + 𝑐
According to article (see [35]),  0 can be expressed as follows: where   and   are the equivalent damping and stiffness, respectively; then the error between solutions of these two systems is minimized with the mean-square method.The difference between (3) and ( 4) is shown in the following: To get a relative precise result, the error  0 should be approximating to minimum.It is better to solve the following instead of ( 6): In order to choose the best equivalent damping   and the equivalent stiffness   , it is necessary to minimize the error with statistical procedure, which requires (7) to be approximating to minimum.
According to the method of multivariate function, the necessary and sufficient condition (see [36]) for the minimum of [ 0 2 ] is obtained; it requires that Equations ( 10) lead to two linear equations and determine the optimal values of   and   .
The required parameters can be obtained simultaneously as follows: It is known from the paper (see [37,38]) that   and   determined by the above formula lead to the minimum value of [ 0 2 ].It is important to note that it has to solve the linear random vibration system (4) to obtain the optimal values of   and   .

Statistical Characteristics of Displacement
and Velocity Response of Weak Nonlinear Energy Dissipation System under Biaxial Earthquake Action

The Transform of the Time Domain Dynamic Equation.
The motion equation of equivalent linear structure with viscoelastic dampers (4) could be written in the following form: where where the symbols  1 ,  1 , and  0 are structure self-vibration frequency, damping ratio, and the reciprocal of structure mass, respectively.Moreover,   and   are the equivalent damping and stiffness, respectively.According to the seismic code [39],   should be ascertained by the maximum between the following: So   can be determined by the following: where ü  and ü V are the horizontal and vertical acceleration, respectively.
Assume that So the time domain dynamic equation of the energy dissipation structure of a single degree of freedom with linear viscoelastic damper could be expressed in the following form:

Stochastic Averaging Equation.
According to the stochastic averaging theory, the standard Van-der-Pol transform is introduced: The stochastic averaging equations that fit the amplitude  1 () are shown in the following: where V 1 () and V 2 () are Wiener process of independent units and   1 ( 1 ) is the power spectrum function of  1 in the value of  1 ; the expression of  is shown in (22). where

The Transient Joint Probability Density Function of Each
Mode Shape of the Nonlinear Structure with Braces.Assume that the state variables of  1 () and Φ 1 () are  1 and  1 , respectively.Probability density function of  1 () is  1 ( 1 , ).The transient joint probability density function of  1 () and Φ 1 () is  1 ( 1 ,  1 , ) and the transient joint probability density function of () and u () is  1 (, u , ), where () is structure displacement and u () is the velocity.According to Itô equation (21), the transient joint probability density function  1 ( 1 ,  1 ,  |  0 ,  0 ,  0 ) that fits the FPK equation is shown in the following: Because (20) does not depend on Φ 1 (), the probability density function ( 1 ,  |  0 ,  0 ) determined by FPK equation is as follows: The initial conditions of ( 24) and ( 25) are, respectively, as follows: Comparing with ( 24) and ( 25), we obtain the relationship of solution under the static initial conditions the following: Meanwhile, we obtain the transient joint probability density function of the original weak nonlinear structure from transient displacement () and transient velocity u () under the static initial condition.
Assume that the form of  1 ( 1 , ) is described as follows: where  1 () is the undetermined function.Equation ( 31) is substituted into (28); we transform the system of (31) into the following form.Then ( 32) is substituted into (31); we can obtain the analytical solution of  1 ( 1 , ).
According to ( 29) and ( 32), we can obtain the response variance of the structural displacement and velocity, respectively.

Equivalent Damping of Weak Nonlinear Structure with the Viscoelastic Damping and the Braces
The actual ground motion is highly random characteristics.
Because of the rationality and practicality of the earthquake, the ground motion model still needs to be further improved.So the the response spectrum method is adopted in most countries.Once the structure is installed with the damper and it turns into an energy dissipation structure, the response spectrum method can not be directly applied to these structures.Therefore, it is greatly significant to establish the equivalent structure which can be used directly with the response spectrum method.The calculation diagram is shown in Figure 3.
Where  0 1 () = ∫  0 ℎ 1 ( − ) u () is the equivalent to a damping force of   u from (4), the motion equation of the structure can be described as follows: In this case, (35) may be written as the following form: where According to the stochastic averaging method, it is known that the probability density function of the amplitude response ( 1 ()) of the equivalent structure is  1 ( 1 , ).The probability density function fitting the FPK equation is as follows: The amplitude probability density function of the original structure can be applied to (30); the amplitude probability density function of the equivalent structure is appropriate for (37).We will know the difference by comparing with (30) and (37).After the following processing, the expression can be expressed as follows: where   is the equivalent damping ratio of damper; it is consistent with the equivalent damping ratio of the Maxwell damper with the general integral model.For arbitrary random biaxial earthquake excitations ü  () and ü V (), all stochastic response characteristics calculated with the proposed method in equivalent structure are the same as these of the original structure.The equivalent damping ratio of the whole weak nonlinear dissipation structure is established as follows: That is, the equivalent structure can be used as a total equivalent ratio of   instead of the original structure damping ratio  1 ; then we can use response spectrum method for structural analysis and engineering design.

Numerical Example
It shows a SDOF nonlinear generalized Maxwell damper energy dissipation structure and the equivalent structure in Figure 4; the earthquake intensity is 8 degrees (0.2 g); its mass, stiffness, damping, and damping ratio are, respectively,   = 2 kg,   = 100 N/m,   = 2 N⋅s/m, and  1 = 0.05.The nonlinear structure is subjected to transient forces under biaxial earthquake.  1 =  1 (0) = 500 × 10 −6 (m 2 /s 3 ),  = 0.2 s.The performance parameters of Maxwell damper in parallel are listed as follows: the brace  1 = 200 N/m, equilibrium modulus  1 = 200 N/m, ℎ 1 = 200 s −2 , element damping coefficient  0 = 30 N⋅s/m, and the stiffness  0 = 50 kN/m.The excellent frequency and damping ratio of the site are  1 = 9.67 s −1 and   1 = 0.9, respectively.Spectral intensity factor  0 = 0.01387 m 2 /s 3 .According to the equivalent damping ratio formula, when  = 0.2 s, the attached equivalent damping ratio   of damper and the response variance of equivalent structural displacement are calculated; the response variance of original structure is also obtained by the frequency domain method.
According to (2), (14), and (35), we can obtain the value of the following parameters: Mathematical Problems in Engineering Hence, The total coefficient of the parallel spring group is equal to the sum of the coefficients of each spring: According to (36),   can be calculated as follows: From ( 32) and ( 34), we can conduct the following calculations: (45) Hence, we can obtain the following parameters values: (46) According to the frequency domain method, frequency response function and the variance of displacement are obtained, respectively. where It is known that we have calculated the maximum displacement standard deviation by frequency domain method and equivalent structure.The results of maximum displacement standard deviation are given in Table 1.Results of the two methods are gradually approaching with the increase of the damping coefficient.The maximum displacement relative error is gradually reduced with the increase of the damping coefficient.When  0 increases to a certain value, the results have a higher precision accuracy.

Conclusions
In this paper, a weak nonlinear structural system with one degree of freedom is researched and a systematically research on the random response characteristic of structure was conducted, which is under biaxial earthquake action.First, integral constitutive relation is adopted; it then establishes a differential and integral equations of motion of SDOF weak nonlinear structure containing the general integral model viscoelastic dampers and the braces.And, then, the motion equation is linearized according to the principle of energy balance.Finally, based on the stochastic averaging method, the general analytical solution of the variance of the displacement, velocity response, and equivalent damping is deduced and derived.The joint probability density function of the amplitude and phase and displacement and velocity of the energy dissipation structure are also given at the same time.Numerical example shows the availability and accuracy of the proposed method.It means it has established a complete analytical solution of stochastic response analysis and equivalent damping of a SDOF nonlinear dissipation structure with the brace under biaxial earthquake action in this paper.The proposed method provides a beneficial reference for the engineering design of this kind of structure.

Table 1 :
Results comparison of the frequency domain method and the proposed method in this paper.