CC BY
7
#
Discrete & Continuous Models
2019, 27 (2) 124-132 http://journals.rudn.ru/miph
& Applied Computational Science
Mathematical modeling
Research article
UDC 519.6;699.841
DOI: 10.22363/2658-4670-2019-27-2-124-132
Seismic stability of oscillating building on kinematic
supports
Sergei P. Karnilovich, Konstantin P. Lovetskiy, Leonid A. Sevastianov, Eugene L. Shchesnyak
Peoples' Friendship University of Russia (RUDN University) 6 Miklukho-Maklaya St., Moscow, 117198, Russian Federation
(received: February 4, 2019; accepted: November 14, 2019)
The design of kinematic supports is considered, which allows to damp the oscillation energy of seismic waves during earthquakes. The building rests on supports that have the geometry of straight cylinders. When horizontal ground oscillations occur, the supports are deflected at a small angle "0. At the same time, their centre of gravity rises and tends to return to its original position under the action of two forces on each support: the weight of the building evenly distributed over all the supports, and the weight of the support itself. The first force is applied to the highest point of the support, the second one is applied to the centre of gravity of the support, so that the rotational moments of two forces act on the support.
It should be noted that under very strong vibrations of the ground, the projection of the centre of gravity could move beyond the base of the support. In this case, the supports will begin to tip over. We confine ourselves to considering such deviations that the rotational moments of the forces of gravity still tend to return the supports to their initial state of equilibrium.
Key words and phrases: ensuring seismic stability of buildings during earthquakes, the equation of motion of a physical pendulum, vibration damping.
The amount of energy transferred to the building depends on the relation between the spectra of seismic effects and natural oscillations of the building. The closer the peaks of the spectra, the greater the energy transferred to the building under similar conditions. This energy is mainly absorbed by the inelastic deformations of the structure. Based on the above facts, two main tasks can be formulated aimed at ensuring the seismic resistance of a building:
— to separate the spectra and thereby reduce the amount of energy transmitted to the building and
© Karnilovich S.P., Lovetskiy K.P., Sevastianov L.A., Shchesnyak E.L., 2019
1. Introduction
© ®
This work is licensed under a Creative Commons Attribution 4.0 International License http://creativecommons.org/licenses/by/4-0/
— to provide the absorption of the remaining energy using special measures.
In the publications on deterministic analysis of the dynamics of constructions, describing the time-dependent motion of the system under the action of seismic load (strongly oscillating and irregular) [1-8], various types of the building supports are studied, aimed at efficient damping the energy of the spatial movement of the building caused by an earthquake. Among them the most successful solutions have been proposed by A. M. Kurzanov and Yu. P. Cherepinsky [6]. In this paper, we consider in detail the design and operation of Kurzanov's kinematic supports, which have been well-proven in experimental studies [1,7,8]. Our theoretical study of the functioning of kinematic supports is aimed at their mathematical modelling and the subsequent selection of parameters of the model of supports in order to solve both of the above problems of seismic stability of buildings. Note that in recent years the interest in the study of kinematic supports that provide the seismic resistance of buildings has noticeably increased [9-11]. See also [12-19].
2. Mathematical model of the functioning of kinematic supports
The design of kinematic supports is considered, which allows damping the oscillation energy of seismic waves during earthquakes. The building having the weight M g rests on n supports, each having the weight mg. The supports have the geometry of straight cylinders of height h and base diameter a.
Under the action of an incident seismic wave, the whole construction "building + n supports" comes into a complex movement. Of all the seismic waves, let us consider horizontal waves. When horizontal ground oscillations occur, the supports are deflected by an angle -0. At the same time, their centre of gravity rises and due to this fact the building, moving horizontally, rises and acquires additional potential energy. After that, the building tends to return to its original position and acquires additional kinetic energy, with the result that each of the supports acquires additional kinetic energy under the action of two forces on each support: the weight of the building evenly distributed on each support, and the weight of the support itself. The first force is applied to the highest point of the support (see Figure 1), the second one is applied to the centre of gravity of the support (see Figure 2), so that the rotational moments of two forces act on the support. Below we consider only a part of the system, namely, a support isolated from the seismic impact.
Each support has the shape of a cylinder with the height h and diameter a = 2r. It should be noted that under very strong vibrations of the ground, the projection of the centre of gravity can move horizontally beyond the base of the support, i.e., p > a, where a is the angle between the diagonal of the support and its height, so that a = arctan(^). In this case, the supports will begin to tip over. We confine ourselves to considering such deviations, when the rotational moments of the forces of gravity still tend to return the supports to the initial state of equilibrium, that is, we will consider the case p < a. Then a = p + 0, where 0 is the angle of deviation of the flat base of the support from the horizontal, while the supports perform non-linear oscillations. To compose the equation of motion of a support, we take into account that it is a physical pendulum.
Figure 1. Rotational force generated by the pressure of the building on the supports
Figure 2. Rotational force generated by the weight of the support
Before the appearance of seismic effects, the "building + n supports" system is in a state of stable equilibrium, i.e., the supports stand on the ground, the building rests on the supports, with the centre of gravity in the lowest position and the lowest potential energy. After the deviation of the "building + n supports" system from the equilibrium position of the supports (along with the building), they rise and fall, rotating around one of their edges during the first half-period of the oscillatory motion. Then the supports touch their bases to the ground with a blow, followed by the loss of a part e of the energy due to the inelastic impact of the ground. After that, due to the remaining (1 — e) kinetic energy of the support (along with the building), they rise and fall, rotating around its other edge during the second half-period of the oscillatory motion. Then the next blow occurs and the oscillatory motion continues with damping.
During each of the half-periods, the movement occurs under the action of gravity forces and their rotational moments.
The returning force of a uniformly distributed building weight acting on a support is FM = — ^ sin p. The distance from the point of application of the force FM to the axis of rotation is equal to LM = 21 = \J a2 + h2. The torque of the force is expressed as FM • LM = — a2 + h2 sin
The returning force generated by the weight of the support itself is calculated as Fm = mg sin (p. The distance from the point of application of the force FM to the axis of rotation is I. The torque of the force is calculated using the formula Fm * I.
The moment of inertia of the support is equal to J = ^r2 + ^h2. The equation of motion of a support, assuming that it is a physical pendulum, has the form
Fm + FM = —J(p. (1)
Taking into account the explicit form of the restoring forces, we obtain the relation
/■ 2Mg 12Z2 + h2 ..
mql sin --/sin p =--—-p, (2)
n 24 v 7
which is reduced to the Lagrange differential equation describing the dynamics of the motion of the supports after a seismic shock, bringing the entire system out of equilibrium:
=
( 12Z2 + h2 ' Si
n
Sin ip.
In the case of limited oscillations, i.e., when ^ < a, we get
•• (nm + 2M\ 24gl
— vT2]2+~h2J^" sin(«-V0=°.
(3)
(4)
The obtained Eq. (4) is the equation of free oscillations (not disturbed by the continuing seismic effect).
Figure 3. Position of the support at the moment of maximum lifting of the centre of gravity (t = 0)
Figure 4. The support of the building. Maximum deviation after changing the axis of rotation
For a complete description of the evolutionary process of support oscillations, we supplement Eq. (4) with the initial conditions, taking the maximum deflection of the support after the seismic shock shown in Figure 3 for the initial position at zero time. Then at the moment of time corresponding to the highest ascent of the centre of gravity of the support (the centre of rotation of the support is located at its lower right point, see Figure 3) the movement of an individual support is described by Eq. (4) with the initial
conditions (we assume the deflection angle positive)
A (0) = ^ const, A (°) : ^const > 0. (5)
Under the influence of the gravity of the support and the entire building, the centres of gravity of the supports will tend to return to their original position. In this case, the angle 01 >0, 01 ^ 0, in accordance with the solution of the Cauchy problem (4), will tend to zero. At that moment, when the magnitude of the angle equals zero, the inertial forces will force the building and the supports to move further in the same horizontal direction, which will lead (after the impact of the support on the base surface) to a change of the rotation axis of the supports, see Figure 4. The motion of the support system and the building itself is still described by exactly the same equations, but with the variable changed for 02 :
•• (nm + 2M\ 24gl .
^ - [up^)— sin(*-^ ) = 0 (6)
and with other initial conditions (again we consider the deflection angle 02 to be positive)
A (0) = 0
À (0)=À(0). ()
It is important to note that at the time of the collision of the support with the ground surface, the system consisting of the building and the supports loses some fraction of the kinetic energy. Setting the restitution coefficient to be equal to Cr = (1 — e) < 1 and considering the kinetic energy losses, we arrive at the relation
~2^2after = Cr^befor ^ Rafter = V^rfibefor, (8)
which, in turn, makes it possible to determine the new velocities of the supports after their collisions with the surface and recalculate the initial conditions of the problems (5) and (7) when going through the zero value of the rotation angle using the formulas
i>k+1 (t + 0) = ^k (t — 0), À+1 (t + 0) = —Vc~4k(t — 0), ( )
for the generalized equation of motion of the supports and the building:
. . .. / nm + 2M\ 24gl . , . , ,, ,, ip = sign(y) (12l2 + nm sin(a — sign(V№). (10)
The formulated system can be solved by means of any stable numerical method, for example, the 4ife-order Runge-Kutta method with automatic step selection. We get the following plots for the solutions (Figures 5-7).
Angle, in rad
Figure 5. The angle of deviation from the vertical. The solid line corresponds to the recovery coefficient Cr = 0.9, the dashed line — to Cr = 0.8
rad/sec
corresponds to the recovery coefficient Cr = 0.9, the dashed line — to Cr = 0.8
3. Conclusion
The model (4) is a nonlinear conservative dynamical system. The motion of the supports and the building body during the non-linear oscillations can be considered as oscillations of the coupled physical pendulums. At the same time, the coupling of pendulums is not conservative, but contains a factor proportional to the rolling friction of the supports on the base of the building body. Moreover, the rolling radius depends on the angle: the larger the angle the smaller the radius, which means greater friction force between the n supports and the base of the building.
Figure 7. Phase portrait of the "oscillatory system" (Cr = 0.8).
In addition, of course, this entire unified system, including the reaction of the soil, moves under the action of a seismic strongly oscillating and irregular disturbance generated by an earthquake. Only as a result of the joint consideration of all these factors and we can count on an adequate description of the deterministic dynamics of the building under the influence of seismic perturbations from an earthquake.
However, even the analysis of equation (4) allows detuning of the natural frequencies of free vibrations of a building system on supports by varying the weight of the building and the number of supports.
The equation obtained by us is valid for angles ^ < a, but it can be easily generalized to the case ^ > a, when the movement of the supports will cease to be oscillatory, and the whole structure will lose stability. However, this case, interesting while considering the driving forces of horizontal seismic vibrations, will be analysed elsewhere.
Acknowledgments
The publication has been prepared under support from the "RUDN University Program 5-100" and the Russian Foundation for Basic Research, grants No. 18-07-00567 and No. 18-51-18005.
References
1. A. M. Kurzanov, S. Yu. Semenov, Dynamic tests of a multi-storey monolithic house in Sochi [Naturnyye dinamicheskiye ispytaniya stroy-ashchegosya mnogoetazhnogo seysmo izolirovannogo monolitnogo doma v Sochi], Industrial and civil construction [Promyshlennoye i grazhdan-skoye stroitel'stvo] (3) (2005) 42-43, in Russian.
2. A. K. Chopra, Dynamic of structures. Theory and applications to earthquake engineering, Prentice-Hall, New Jersey, 2006.
3. A. Martelli, M. Forny, Seismic isolation: present application and perspectives, Yerevan, Armenia, 2006, pp. 1-26.
4. H. Masahiko, O. Shin, Response Control and Seismic Isolation of Buildings, Taylor & Francis, New York, 2006.
5. V. L. Kharlanov, Deterministic analysis of metal frames for high-intensity dynamic loads: a monograph [Determinirovannyy analiz metallicheskikh karkasov na dinamicheskiye nagruzki vysokoy intensivnosti], VolgGASU, 2006, in Russian.
6. G. E. Avidon, E. A. Karlina, Vibration propertie of buildings with base isolating foundations by A. M. Kurzanov and Yu.D.Cherepinsky [Osobennosti kolebaniy zdaniy s seysmoizoliruyushchimi fundamentami A. M. Kurzanova i Yu. D. Cherepinskogo], Earthquake engineering. Building safety [Seysmostoykoye stroitel'stvo. Bezopasnost' sooruzheniy] (1) (2008) 42-45, in Russian.
7. A. M. Kurzanov, Modern state of standardization of design of structures for seismic load [Sovremennoye sostoyaniye normirovaniya rascheta sooruzheniy na seysmicheskuyu nagruzku], Industrial and civil construction [Promyshlennoye i grazhdanskoye stroitel'stvo] (11) (2009) 52-53, in Russian.
8. E. L. Schesnyak, et al., STO RUDN 02066463-001-2011 "Standards for Designing Structures in Seismic Areas" [STO RUDN 02066463-0012011. Stroitel'stvo v seysmicheskikh rayonakh. Normy proyektirovaniya sooruzheniy], RUDN, Moscow, 2011, in Russian.
9. S. P. Karnilovich, K. P. Lovetsky, L. A. Sevastyanov, E. L. Schesnyak, Seismic systems based on A.M. Kurzanov's kinematic supports [Seysmoizolyatsiya zdaniy na osnove kinematicheskikh opor Kurzanova A.M.], in: XIX international Conference on Distributed Computer and Communication Networks: Control, Computation, Communications (DCCN-2016). Proceedings of the Nineteenth International Scientific Conference "Mathematical Modelling, Simulation and Control Problems", Vol. 2, 2016, pp. 159-164, in Russian.
10. O. V. Druzhinina, L. A. Sevastianov, S. A. Vasilyev, D. G. Vasilyeva, Lyapunov stability analysis for the generalized Kapitza pendulum, Journal of Physics: Conference Series 937, article number 012011 (2017). doi:10.1088/1742-6596/937/1/012011.
11. O. V. Druzhinina, L. A. Sevastianov, S. A. Vasilyev, D. G. Vasilyeva, Numerical analysis of kurzanov bearing oscillation, in: M. B. Kochanov (Ed.), 7th International conference "Problems of Mathematical Physics and Mathematical Modelling": Books of abstracts (Moscow, NRNU MEPhI, 25-27 June), National Research Nuclear University MEPhI, Moscow, 2018, pp. 164-166.
12. Y. M. Parulekar, G. R. Reddy, Passive response control systems for seismic response reduction: a state-of-the-art review, International Journal of Structural Stability and Dynamics 09 (01) (2009) 151-177. doi:10.1142/S0219455409002965.
13. M. Palermo, S. Silvestri, L. Landi, G. Gasparini, T. Trombetti, A "direct five-step procedure" for the preliminary seismic design of buildings with added viscous dampers, Engineering Structures 173 (2018) 933-950. doi:10.1016/j.engstruct.2018.06.103.
14. S. Wang, S. A. Mahin, High-performance computer-aided optimization of viscous dampers for improving the seismic performance of a tall building, Soil Dynamics and Earthquake Engineering 113 (2018) 454-461. doi:10.1016/j.soildyn.2018.06.008.
15. C.-M. Chang, S. Shia, C.-Y. Yang, Design of buildings with seismic isolation using linear quadratic algorithm, Procedia Engineering 199 (2017) 1610-1615. doi:10.1016/j.proeng.2017.09.069.
16. J. M. Kelly, The role of damping in seismic isolation, Earthquake Engineering and Structural Dynamics 28 (1) (1999) 3-20. doi:10.1002/(SICI)1096-9845(199901)28:1<3::AID-EQE801>3.0.C0;2-D.
17. I. Politopoulos, A review of adverse effects of damping in seismic isolation, Earthquake Engineering and Structural Dynamics 37 (3) (2008) 447-465. doi:10.1002/eqe.763.
18. G. P. Warn, K. L. Ryan, A review of seismic isolation for buildings : historical development and research needs, Buildings (2) (2012) 300-325. doi:10.3390/buildings2030300.
19. S. Moretti, A. Trozzo, V. Terzic, G. P. Cimellaro, S. Mahin, Utilizing base-isolation systems to increase earthquake resiliency of healthcare and school buildings, Procedia Economics and Finance 18 (2014) 969-976. doi:10.1016/S2212-5671(14)01024-7.
For citation:
S. P. Karnilovich, K. P. Lovetskiy, L. A. Sevastianov, E. L. Shchesnyak, Seismic stability of oscillating building on kinematic supports, Discrete and Continuous Models and Applied Computational Science 27 (2) (2019) 124-132. DOI: 10.22363/2658-4670-2019-27-2-124-132.
Information about the authors:
Sergei P. Karnilovich (Russian Federation) — Associate Professor, Ph.D., Assistant Professor Institute of Physical Research and Technology of Peoples' Friendship University of Russia (RUDN University) (e-mail: karnilovich-sp@rudn.ru, phone: +7(495)434-42-12, Scopus Author ID: 6505797810)
Konstantin P. Lovetskiy (Russian Federation) — Associate Professor, Ph.D., Assistant Professor of Department of Applied Probability and Informatics of Peoples' Friendship University of Russia (RUDN University) (e-mail: lovetskiy-kp@rudn.ru, phone: +7(495)955-07-83, ORCID: https://orcid.org/0000-0002-3645-1060, ResearcherID: A-5725-2017, Scopus Author ID: 18634692900)
Leonid A. Sevastianov (Russian Federation) — professor, Doctor of Physical and Mathematical Sciences, professor of Department of Applied Probability and Informatics of Peoples' Friendship University of Russia (RUDN University) (e-mail: sevastianov-la@rudn.ru, phone: +7(495)955-07-83, ORCID: https://orcid.org/0000-0002-1856-4643, ResearcherID: B-8497-2016, Scopus Author ID: 8783969400)
Eugene L. Shchesnyak (Russian Federation) — professor, Doctor of Economics, First Vice-Rector, Vice-Rector for Economic Activity of Peoples' Friendship University of Russia (RUDN University) (e-mail: shcheasnyak-el@rudn.ru, phone: +7(495)433-73-68)
