Abstract

This study aims to stabilize the unwanted fluctuation of buildings as mechanical structures subjected to earth excitation as the noise. In this study, the ground motion is considered as a Wiener process, in which the governing stochastic differential equations have been presented in the form of Ito equation. To stabilize the vibration of the system, the ATMD system is considered and located on the upmost story of the building. A sliding mode controller has been utilized to control the ATMD system, which is a robust controller in the presence of uncertainty. For this purpose, the design of a sliding mode controller for the general dynamic system with Lipschitz nonlinearity and considering the Ito relations has been accomplished. The mentioned design has been implemented considering the presence of the Weiner process and existence of uncertainty in the structure and actuator. Then, the obtained general control law has been generalized to control the ATMD system. The results show that the designed controller is effective to reduce the effect of the unwanted impused vibrations on the building.

1. Introduction

Stochastic factors are one of the inherent aspects of most dynamical systems that occur in different ways such as external force and changes in the inherent parameters of the system. Stochastic Differential Equations (SDE) in financial mathematics and economics have been extensively investigated [16]. Also, the effect of stochastic factors in science and mechanical and electronic engineering has been considered [720].

Earthquake is an example of a natural phenomenon that influences the dynamics of a structures and building. Dynamic behavior of structures in the presence of earthquakes has been extensively studies in [2127]. In the mentioned studies, the dynamic behavior of the structures has been investigated under the influence of a particular earthquake. But in this paper, we aim to examine the dynamic behavior of the structure subjected to White Gaussian Noise (WGN) because the density function of its power spectrum is constant at all frequencies. So, unlike the previous studies, a new formulation should be considered.

For this purpose, Itô formulation is considered to solve governing SDE of the structure [2830]. The studied structure is an 11-story building equipped with an ATMD system at the upmost story. ATMD has been used to reduce unwanted vibrations. This system has been controlled by means of the sliding mode controller. For this reason, the sliding mode controller has been designed for the general and nonlinear Lipschitz dynamic system in the presence of actuator and system uncertainties and based on the Itô theory. Finally, the designed controller has been generalized to control the ATMD system.

The dynamic behavior of the structure in the active and passive mode and considering the uncertainties has been studied. The effect of various controller parameters on the dynamic behavior of the structure has been also investigated. Moreover, the effect of the controlled system and different parameters of the controller on the basin of attraction was also studied. At the end, the controller’s robustness to structural and actuator uncertainty has been studied and the obtained results have been presented.

2. Model Description

The physical and geometrical models of the studied system are presented in this section. The studied model is an 11-story building, the schematic view of which is shown in Figure 1. The mass of each story is denoted by , where i represents the story number. Moreover, the stiffness and damping coefficient for each story are represented by and , respectively. The degree of freedom (DoF) of the system is considered along the horizontal direction, since the ground displacement effects are horizontally applied to the structure. As shown in Figure 1, the displacement of each story is denoted by . Upon ground excitation, each story experiences a vibration and, given that the first mode of vibration is more likely to occur, the largest displacement occurs in the topmost story. To alleviate this deformation caused by earthquakes, the 11th story is equipped with an ATMD system. The mass of this system, its stiffness, and damping coefficient are denoted by , , and . After installation, the DoF of system increases to 12.

The dynamic equation governing the system behavior is presented in the following equation [31]:where , , and denote mass, stiffness, and damping matrices, respectively. Moreover, , , and represent the structural displacement, velocity, and acceleration, respectively. The term is the control force vector of the ATMD system responsible for controlling the vibrations of the 11th story. The term denotes the effects of base excitation on the system, where is the acceleration of base excitations and . The nominal matrix values for , , and are also presented in the following. However, note that uncertainty will be also considered for these values later in the control design process.

3. Mathematical Modeling

Without loss of generality, the controller is initially designed for a general system, after which it is applied to the studied structural control.

Consider the mathematical model presented in the following equation:where , , and represent continuous functions satisfying the Lipschitz condition. In this case, , , and indicate the control force, standard Wiener process, and white Gaussian noise, respectively. The objective of the controller design is for the to track . To this end, the dynamic error of should be defined. However, since is considered zero in structural control, the error is defined as and the corresponding dynamic error is defined as follows:

The Lyapunov function is considered as . From a mathematical viewpoint and based on Ito’s theory, equation (4) can be reformulated in the form of a differential equation as follows:

Given that , the abovementioned equation is presented as follows for simplification purposes:

The abovementioned equation represents an Itô SDE used instead of equation (4) and considering the Wiener process. The terms and represent the drift function and diffusion, respectively [28, 29].

The sliding surface was considered as in the design of the sliding mode controller, from which and can be derived. In the final form and according to equation (4), the equation can be rewritten to obtain equation (5):

Assuming and employing Itô ’s differentiation formula for , we have [28, 29]

Based on Ito’s formula, differentiating produces . Therefore, the value emerges as follows:

By including the following relations proposed by [32, 33] in the calculations,

And also considering the properties of the Wiener process [28],

According to equations (10) and (11), the expected value of is determined as follows:

Dividing the abovementioned relation by , is obtained as follows:

The stability condition for the sliding mode controller is defined as based on the Lyapunov second method for stability [34]. Assuming the systems involve no uncertainty, the controller stability is only guaranteed by considering , limiting the region of attraction associated with the sliding surface. However, given the presence of uncertainty in most of dynamical systems, the system uncertainties are included in the controller design equations in the following part.

The structural and actuator uncertainties have been considered in the controller model in this study. To this end, some of the inequalities associated with the and functions along with their nominal values should be taken into account.

Assuming as the nominal values for and is a positive function expressed as follows:

As a result, the following equations hold true:

If is defined aswhere and are positive values representing the upper and lower bounds of the function, the following equations are necessary for the controller design:

The following inequality holds true considering equations (15), (17), and (19):

Equation (21) also holds true considering equations (18) and (20):

By defining as follows, the relation holds true:where represents a positive value and . In fact, , where and , which are obtained as follows, cause the nominal and uncertainty terms to be negative definite, respectively:

Theory. assume as the following set:where is a positive value. In fact, represents the attraction set for the trajectory . Assuming at , the trajectory lies outside the attraction set, i.e., , we obtainwhich indicates a declining rate for as it tends to enter the set and ultimately remains in this region.
If it is located in the domain of for the moment , it will continue to be in this area.
The designed sliding mode controller is then applied to the system, as shown in Figure 1, based on the abovementioned equations. The equation governing the dynamic behavior of the 11th story is expressed as follows:The for this equation is defined as follows:where , , , , and denote nominal values and , , , , and represent the maximum value for the respective uncertainty term.
In this case, is defined as follows:Moreover, is assumed as in this case, meaning that the system’s actuator includes the uncertainty . Under these conditions, and are expressed as follows:Finally, the term associated with the control force is conveniently determined using equation (22).

4. Results and Discussion

In this section, the obtained results will be presented and the control of horizontal displacement of the 11-story building will be discussed. Given that the first vibrational mode of the building is easily excited, the largest displacement in this mode is experienced by the topmost story. Therefore, vibration analysis is conducted on the 11th story. The physical and geometrical specifications of the studied building are given in Table 1. In the first scenario, the uncertainty was assumed as , , , and . The resulting horizontal displacement for story 11 is demonstrated in Figure 2 for both active and passive cases. As shown in this figure, the vibrational amplitude of the controlled system is considerably smaller than that of the uncontrolled system.

The variations of the control force with respect to time are demonstrated in Figure 3. As shown, the chattering phenomenon is apparent within some time intervals. Chattering is mainly caused by the term in the control force relation. In fact, the discontinuity and undifferentiability of this function at point is responsible for the chattering phenomenon. The chattering around the zero point is highly harmful, since, in addition to the force magnitude, its sign also changes. However, at nonzero points, the chattering only causes a decrease or increase in the force magnitude. In our case, as shown, chattering occurs at nonzero points.

To resolve the chattering problem and satisfy the continuity and Lipschitz condition, the term was used as an approximation of the function for the function , where the continuity and differentiability conditions are satisfied at point . The diagram of horizontal displacements for different values is demonstrated in Figure 4. As shown, the displacement amplitude is increased by increasing . However, this increase is negligible compared to the displacement of the uncontrolled system.

The variations of force corresponding to Figure 4 for different values are shown in Figure 5. As indicated, the amplitude of force at is smaller compared to other values. Moreover, the chattering phenomenon is also fully resolved in this case, while it is still observed at other cases, for example, . As shown in Figures 4 and 5, the vibrational amplitudes for both and functions are consistent, which is due to the accurate approximation of the function by .

Moreover, the phase diagram for the 11th story is also demonstrated in Figure 6 for different cases. As shown, by increasing , the region of attraction is also extended. This region was almost similar for both and functions.

The effect of on the dynamic behavior of the structure and the region of attraction is discussed in this part. The vibrational amplitude of the uncontrolled system for different values is shown in Figure 7. As demonstrated, the displacement amplitude is linearly increased by increasing .

Also, the attraction domain for the mentioned figure is presented in Figure 8. As is clear from this figure, the domain of attraction set is extended with increasing .

The horizontal displacement for the controlled case is demonstrated in Figure 9 with respect to different values. As shown, by increasing , the vibrational amplitude is also increased, but its value is substantially smaller compared to the uncontrolled case. The force variations corresponding to Figure 9 is demonstrated in the diagram of Figure 10. As shown, the amplitude of force variations also experiences an increase as increases.

Moreover, the region of attraction corresponding to Figure 9 is demonstrated in Figure 11. As shown, similar to the uncontrolled case, the region of attraction increases by increasing .

The effects of and on the behavior of a controlled system is discussed in this section. The variations in horizontal displacement of the system are demonstrated in Figure 12 for different values. As shown, for values lower than 1000, increasing the values is not significantly effective, but increasing this parameter to 100,000 causes the horizontal displacement of the system to increase substantially.

Additionally, the region of attraction for an value of 100,000 is shown in Figure 13. As shown, the region of attraction is considerably expanded at this value.

The displacement amplitude of the system for different values and a fixed value of 1 is demonstrated in Figure 14. As shown, the displacement amplitude is substantially decreased by increasing lambda. The region of attraction for a value of 5 is shown in Figure 15. Consider this figure suggests that the region of attraction becomes more limited for a lambda value of 5.

This section discusses the controller robustness to system uncertainty. In the first scenario, only the system uncertainties were considered and the actuator uncertainties were neglected. The variations of horizontal displacement for different values are demonstrated in Figure 16. As shown, the controller is highly robust to uncertainties and even within some time intervals, and the horizontal displacement of the system is decreased by increasing .

This is due to the fact that is also increased by increasing . However, as expressed in equation (22), presence of in this equation virtually increases the term , and as shown in Figure 12, increasing increases the controller robustness. The regions of attraction for and are plotted in Figure 17. As shown, despite its larger extent for compared to , the region of attraction has become more compact.

The results for the case with actuator uncertainty and were extracted and presented in Figure 18. In this case, . The results for different values are shown in this figure, which indicates the significant controller robustness in presence of actuator uncertainty.

The region of attraction in this case is depicted in Figure 19. As shown, the extent of region of attraction is increased by increasing .

The results in presence of actuator and system uncertainties are shown in Figure 20. In this case, is set to be 0.2. As shown, the controller exhibits a high robustness in presence of actuator and system uncertainties. The region of attraction for this case is demonstrated in Figure 21. As shown, the extent of region of attraction is increased by increasing .

5. Conclusion

This study has examined the dynamic behavior of an 11-story building equipped with an ATMD system. The stimulation force has been applied on the structure in the form of a Wiener process and an earthquake. The sliding mode scheme has been used to control the ATMD system. So, considering Lipschitz nonlinearity and based on the Ito formulation, a sliding mode controller in the presence of the uncertainty for the general dynamical system with the second-order governing stochastic differential equation has been designed.

The designed controller has been further developed to control the ATMD system. The dynamic behaviors of the structure in the active and passive modes have been simulated. The presented results demonstrate the high ability of the sliding mode controller to reduce unwanted vibrations of buildings under stimulation of the ground. Also, the results show that the designed controller has a good robustness in the presence of structural and actuator uncertainties. In addition, the effects of the controller parameters on system behavior have been studied. The results show that reduction of and increase of reduce the structural vibration amplitude.

Data Availability

No data were used to support this study.

Conflicts of Interest

The author declares that there are no conflicts of interest.

Acknowledgments

This work was supported by The Research Council of Oman with Grant no. ORG/CBS/14/008.