^{1}

^{*}

^{1}

Continuum plate model in the form of a cantilever anisotropic plate developed in the framework of the bimoment theory of plates describing seismic oscillations of buildings is proposed in this paper as a dynamic model of a building. Formulas for the reduced moduli of elasticity, shear and density of the plate model of a building are given. Longitudinal oscillations of a building are studied using the continuum plate and box-like models of the building with Finite Element Model. Numerical results are obtained in the form of graphs, followed by their analysis.

Among numerous objects of study in mechanics of a deformed rigid body, multi-storey buildings and structures occupy a special place. Development of dynamic models of buildings and structures with consideration of strain of spatial nature is one of the most difficult and urgent problems in mechanics.

A universal model of the building has not been developed yet. This is due to the complex structure, diversity and numerous elements of the building, especially tall large-sized buildings. Therefore, one of the important tasks of modern mechanics and seismic stability of structures is the development of universal models of a building that adequately describes its spatial behavior.

There are numerous articles and monographs devoted to the development of the theory of seismic resistance of a building. Various methods for calculating buildings and structures for seismic actions have been developed, taking into account important factors, such as seismic loading, soil conditions of terrain and structural features of buildings [

This need arises due to the fact that the existing methodology for calculating and designing buildings does not solve the problems of the optimal ratio of the sizes of a building box, the diaphragm, floor elements rigidity and so on.

In this paper, a technique for calculating the continuum plate model of buildings, developed in [

In [

V b o x = d 1 d 2 d 3 , (1)

where d 1 , d 2 ―are the dimensions of the building box in plan, d 3 ―the height of the box.

Then, to determine the mass of the boxes, the following formula is used

m bd = ρ 0 V 0 = ρ bd V box , (2)

here V 0 ―is the sum of the volumes of bearing and room partitions and floors

V 0 = 2 ( d 1 d 2 h 3 2 + d 1 d 3 h 2 2 + d 2 d 3 h 1 2 ) , (3)

where h 1 , h 2 ―is the thickness of bearing and room partitions. h 3 ―is the thickness of the floor.

From relations (1)-(3) a formula for determining the reduced density of a plate model of the building is obtained

ρ bd = ρ 0 ( h 3 d 3 + h 2 d 2 + h 1 d 1 ) . (4)

Remaining reduced elastic characteristics and densities of the building are determined by the formulas:

E 1 ( bd ) = ξ 11 E 0 , E 2 ( bd ) = ξ 22 E 0 , E 3 ( bd ) = ξ 33 E 0 , G 12 ( bd ) = ξ 12 E 0 , G 13 ( bd ) = ξ 13 E 0 , G 23 ( bd ) = ξ 23 E 0 , ρ bd = ρ 0 ξ 0 . (5)

where E 0 ―is the modulus of elasticity of building material

The values of the coefficients ξ 11 , ξ 22 , ξ 33 , ξ 12 , ξ 13 , ξ 23 , ζ 0 for each building cell (room) are defined as functions of two spatial variables, E 0 , G 0 ―are the moduli of elasticity and shear of the strongest bearing panel of the building. Reduced modulus of elasticity of a building is determined by the formula

ξ 11 = S 11 S 01 , ξ 22 = S 22 S 02 , ξ 33 = S 33 S 03 E 0 , ξ 12 = S 11 S 01 , ξ 13 = h fpl b 1 λ * , ξ 23 = h 2 b , (6)

ζ 0 = V 1 V 0 . (7)

where S 01 , S 02 , S 03 ―are the cross-sectional areas of the building in three coordinate planes of one floor of the building; S 11 , S 22 , S 33 ―are the total areas of cross sections of the plates in the coordinate planes, forming one floor of the building; λ * ―is a coefficient characterizing the voids in the cross-section of the floor plate, equal to

λ * = a 1 h fpl − s * a 1 h fpl . (8)

Here a 1 ―is a distance between two transverse room partitions plates; s * ―is a sum of the voids areas in the cross-section of the floor.

Note that the above volumes and areas are determined, depending on the dimensions of plates, rooms and the building itself, as follows:

S 01 = b 1 H , S 02 = a H , S 03 = a b 1 . (9)

S 11 = 2 b 1 h 1 + b 1 h 2 + H h fpl , S 22 = 2 a h 1 + a h 2 + 2 H h 1 + ( k − 2 ) H h 2 , S 33 = 2 b 1 h 1 + a h 2 + ( k − 2 ) b 1 h 2 . (10)

To represent the values of dynamic characteristics of the plate model of a building, the following formula obtained in [

E bd = ( ω bd ω pl ) 2 ρ bd ρ pl E pl . (11)

The values of coefficients ξ 11 , ξ 22 , ξ 33 , ξ 12 , ξ 13 , ξ 23 are determined for each cell (room) of the building. In general, these coefficients are variables and are the functions of two spatial coordinates, which should be determined for the building in question from multiple numerical theoretical experiments and existing experimental data. The moduli of elasticity of continuum plate model of a building are determined by formulas (11).

The equations of longitudinal oscillations of a thick plate [

∂ N 11 ∂ x 1 + ∂ N 12 ∂ x 2 = ρ H ψ ¯ ¨ 1 , ∂ N 21 ∂ x 1 + ∂ N 22 ∂ x 2 = ρ H ψ ¯ ¨ 2 , (12)

∂ T 11 ∂ x 1 + ∂ T 12 ∂ x 2 − 4 p ¯ 13 = ρ H β ¯ ¨ 1 , ∂ T 12 ∂ x 1 + ∂ T 22 ∂ x 2 − 4 p ¯ 23 = ρ H β ¯ ¨ 2 , (13)

∂ p ¯ 13 ∂ x 1 + ∂ p ¯ 23 ∂ x 2 − 2 p ¯ 33 H = ρ r ¯ ¨ , ∂ τ ¯ 13 ∂ x 1 + ∂ τ ¯ 23 ∂ x 2 − 6 τ ¯ 33 H = ρ γ ¯ ¨ . (14)

The equations of motion for determining the displacements of external longitudinal walls, obtained by meeting the boundary conditions on the faces of the plate z = − h and z = + h using the method of displacements expansion into the Maclaurin infinite power series, are constructed in [

W ¯ = 1 2 ( 21 γ ¯ − 7 r ¯ ) − 1 30 H ( E 31 E 33 ∂ u ¯ 1 ∂ x 1 + E 32 E 33 ∂ u ¯ 2 ∂ x 2 ) , (15)

u ¯ k = 1 4 ( 21 β ¯ k − 3 ψ ¯ k ) − 1 20 H ∂ W ¯ ∂ x k + 1 20 H q ¯ k G k 3 ( k = 1 , 2 ) . (16)

Longitudinal, tangential forces and bimoments are determined with respect to the following nine unknown kinematic functions:

u ¯ k = u k ( + ) + u k ( − ) 2 , ψ ¯ k = 1 2 h ∫ − h h u k d z , β ¯ k = 1 2 h 3 ∫ − h h u k z 2 d z , ( k = 1 , 2 )

W ¯ = u 3 ( + ) − u 3 ( − ) 2 , r ¯ = 1 2 h 2 ∫ − h h u 3 z d z , γ ¯ = 1 2 h 4 ∫ − h h u 3 z 3 d z .

where u i , ( i = 1 , 3 ) ―are the displacements of internal points of the plate; u i ( − ) , u i ( + ) , ( i = 1 , 3 ) ―are the displacements of points on the faces of the plate z = − h and z = + h . Function r expresses an average value of normal displacements of the plate model of the building; W ¯ ―is a half-difference of normal displacements of two outer layers of the plate model; functions u ¯ 1 and u ¯ 2 ―are the half-sums of longitudinal displacements, representing the displacements of the outer layers of the plate model of the building along horizontal and vertical directions.

Expressions of longitudinal and tangential forces are written in the form:

N 11 = E 11 H ∂ ψ ¯ 1 ∂ x 1 + E 12 H ∂ ψ ¯ 2 ∂ x 2 + 2 E 13 W ¯ , N 22 = E 12 H ∂ ψ ¯ 1 ∂ x 1 + E 22 H ∂ ψ ¯ 2 ∂ x 2 + 2 E 23 W ¯ , N 12 = N 21 = G 12 ( H ∂ ψ ¯ 1 ∂ x 2 + H ∂ ψ ¯ 2 ∂ x 1 )

Expressions of longitudinal and tangential bimoments have the form:

T 11 = H ( E 11 ∂ β ¯ 1 ∂ x 1 + E 12 ∂ β ¯ 2 ∂ x 2 + E 13 2 W ¯ − 4 r ¯ H ) , T 22 = H ( E 12 ∂ β ¯ 1 ∂ x 1 + E 22 ∂ β ¯ 2 ∂ x 2 + E 23 2 W ¯ − 4 r ¯ H ) , T 12 = T 21 = H G 12 ( ∂ β ¯ 1 ∂ x 2 + ∂ β ¯ 2 ∂ x 1 ) .

Bimoments are expressed in the form

σ ¯ 11 = E 11 * ∂ u ¯ 1 ∂ x 1 + E 12 * ∂ u ¯ 2 ∂ x 2 , σ ¯ 12 = G 12 ( ∂ u ¯ 1 ∂ x 2 + ∂ u ¯ 2 ∂ x 1 ) , σ ¯ 22 = E 12 * ∂ u ¯ 1 ∂ x 1 + E 22 * ∂ u ¯ 2 ∂ x 2 .

where

E 11 * = E 11 − E 31 E 13 E 33 , E 12 * = E 12 − E 32 E 13 E 33 , E 22 * = E 22 − E 32 E 23 E 33

Intensities of transverse and normal bimoments p ¯ 13 , p ¯ 23 , τ ¯ 13 , τ ¯ 23 и p ¯ 33 , τ ¯ 33 have the form:

p ¯ k 3 = G k 3 ( ∂ r ¯ ∂ x k + 2 ( u ¯ k − ψ ¯ k ) H ) , τ ¯ k 3 = G k 3 ( ∂ γ ¯ ∂ x k + 2 ( u ¯ k − 3 β ¯ k ) H ) , ( k = 1 , 2 )

p ¯ 33 = E 31 ∂ ψ ¯ 1 ∂ x 1 + E 32 ∂ ψ ¯ 2 ∂ x 2 + E 33 2 W ¯ H , τ ¯ 33 = E 31 ∂ β ¯ 1 ∂ x 1 + E 32 ∂ β ¯ 2 ∂ x 2 + E 33 2 W ¯ − 4 r ¯ H .

Intensities of bimoments σ ¯ 11 * , σ ¯ 22 * are expressed as:

σ ¯ 11 * = − E 11 H ∂ 2 W ¯ ∂ x 1 2 − E 12 H ∂ 2 W ¯ ∂ x 2 2 + E 13 R ¯ H , σ ¯ 22 * = − E 12 H ∂ 2 W ¯ ∂ x 1 2 − E 22 H ∂ 2 W ¯ ∂ x 2 2 + E 23 R ¯ H ,

where

R ¯ = 420 ( W ¯ + 6 r ¯ − 15 γ ¯ ) .

The system of differential equations of motion (12)-(16) constitutes a joint system of seven equations with respect to nine unknown functions ψ ¯ 1 , ψ ¯ 2 , β ¯ 1 , β ¯ 2 , u ¯ 1 , u ¯ 2 , r ¯ , γ ¯ , W ¯ .

Let us write down the boundary conditions for the problem of bending and shear vibrations of buildings. Suppose that the foundation points move according to a given law u 0 ( t ) , and the lower part of the building moves horizontally with the foundation. From kinematic considerations it follows that the displacements will be written in the form:

u 1 ( x 1 , 0 , z , t ) = u 0 ( t ) , u 2 ( x 1 , 0 , z , t ) = u 3 ( x 1 , 0 , z , t ) = 0.

From kinematic conditions it follows that in the foundation of the building the following boundary conditions must be fulfilled:

ψ ¯ 1 = u 0 ( t ) , ψ ¯ 2 = 0 , β ¯ 1 = 1 3 u 0 ( t ) , β ¯ 2 = 0 , r ¯ = 0 , γ ¯ = 0 , u ¯ 1 = u 0 ( t ) , u ¯ 2 = 0 , W ¯ = 0.

On the free side faces of the building the conditions of zero force factors are fulfilled:

N 11 = 0 , N 12 = 0 , T 11 = 0 , T 12 = 0 , σ ¯ 11 = 0 , σ ¯ 12 = 0 , p ¯ 13 = 0 , τ ¯ 13 = 0 , σ ¯ 11 * = 0.

On the free top faces of the building the conditions of zero force factors are fulfilled:

N 12 = 0 , N 22 = 0 , T 12 = 0 , T 22 = 0 , σ ¯ 12 = 0 , σ ¯ 22 = 0 , p ¯ 23 = 0 , τ ¯ 23 = 0 , σ ¯ 22 * = 0.

The problem is solved by the Finite Differences Method. To approximate the first derivatives, the central difference schemes are taken.

To justify the expressions obtained with the continuous model of calculation, let us consider a specific problem with the use of necessary physical and mechanical characteristics of a box-like model. Calculation of the box-like models has been carried out by the FEM, the main point of which is the replacement of a real system (a multi-storey box-like building) with a discrete model of rectangular plane elements (walls and floors) connected at the node points and satisfying the equilibrium conditions of the converging system of forces at each node.

At each nodal point i, a part of box mass m_{i} is concentrated, to which convergent forces are applied: the weight m_{i}g, the inertia force Р i = m i у ¨ 0 = m i sin ( 2 π ω 0 t ) and the elastic force [k_{ij}]{q_{j}}, which is the reaction of this point connection to other points of finite element sampling. Here [k_{ij}] is the stiffness matrix and {q_{j}} is the displacement vector found as the result of solving the system of differential equations for the entire box

[ M ] { q ¨ } + [ K ] { q } = [ M ] { g } + [ M ] { q ¨ 0 }

The system of differential equations is solved by the Newmark method.

Mechanical characteristics and geometry of the rooms panels are accepted as follows: bending bearing panels have a modulus of elasticity E = 20000 МPа ; density ρ = 2700 kg / m 3 ; Poisson’s ratio ν = 0.3 . For the panel working on shear, the following mechanical characteristics are adopted: modulus of elasticity E = 7500 МPа , density ρ = 1200 kg / m 3 ; Poisson’s ratio ν = 0.3 .

Results of calculations of forced oscillations of a building within the framework of a thick plate model with the following dimensions of the floor slab of the building are: h 1 = 0.25 m , h 2 = 0.25 m , h fpl = 0.2 m , a 1 = 7.5 m , b 1 = 3 m , a = 30 m , H = 11 m .

In all cases the effect is given by the acceleration of the foundation at amplitude 0.1 g, i.e. at unit amplitude, which corresponds to the intensity of the 7-point earthquake u 0 = sin ( 2 π ω 0 t ) .

The effect frequency is chosen equal to ω 0 = 9.5 Hz . Such a high-frequency effect has a period of T = 0.1 s. The Gazli earthquake (1976) accelerogram had such predominant period. So, it can be said that the considered effect is in a certain sense analogous in predominant period to the Gazli earthquake (1976). The scheme of the building in question is shown in

Results obtained with FEM and FDM for the box-like and the continuum plate model of a building are shown in the form of graphs in

points of the slab on the basis of the box-like model with FEM of a two-storey building and their values are R m i d = 0.0000430 m and R e d g = 0.0000500 m , respectively.

Comparison of results based on the plotted graphs shows that the values of displacements based on the box-like model with FEM are 10% - 13% less than the results based on the plate continuum model.

When comparing results based on plotted graphs for the box-like model with FEM and the plate continuum model, it can be seen that the values of displacements do not differ significantly.

The values of the generalized displacements obtained from the calculations of two-storey and five-story buildings for seismic effects at different points of floor differ within the range of 10%.

In conclusion, it should be noted that on the basis of the bimoment theory of thick plates a dynamic plate-like model of the building has been developed, reflecting its spatial strains. Using the geometry of the building and its rooms, the reduced density, moduli of elasticity and shear of the plate model have been determined. The plate model of the building due to the choice of coefficients ξ 22 , ξ 33 , ξ 12 , ξ 13 , ξ 23 , adequately reflects the laws of displacements of the points of the building. So, the plate model is acceptable for determining the displacement fields of a tall building during an earthquake and the dynamic characteristics of a building.

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

Usarov, M.K. and Ayubov, G.T. (2018) Study of Longitudinal Oscillations of a Building on the Basis of Continuum Plate Model. Journal of Applied Mathematics and Physics, 6, 1928-1936. https://doi.org/10.4236/jamp.2018.69164