Научная статья на тему 'ДЕФОРМАЦИИ И СОБСТВЕННАЯ ЧАСТОТА КОЛЕБАНИЙ ФЕРМЫ ПРОСТРАНСТВЕННОЙ МОДЕЛИ МНОГОЭТАЖНОГО ЗДАНИЯ'

ДЕФОРМАЦИИ И СОБСТВЕННАЯ ЧАСТОТА КОЛЕБАНИЙ ФЕРМЫ ПРОСТРАНСТВЕННОЙ МОДЕЛИ МНОГОЭТАЖНОГО ЗДАНИЯ Текст научной статьи по специальности «Строительство и архитектура»

CC BY
105
10
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Вестник МГСУ
ВАК
RSCI
Ключевые слова
МНОГОЭТАЖНОЕ ЗДАНИЕ / КАРКАС / ВЕТРОВАЯ НАГРУЗКА / СОБСТВЕННЫЕ КОЛЕБАНИЯ / АНАЛИТИЧЕСКОЕ РЕШЕНИЕ / ФОРМУЛА МАКСВЕЛЛА - МОРА

Аннотация научной статьи по строительству и архитектуре, автор научной работы — Кирсанов Михаил Николаевич

Введение. Расчет деформаций под действием различных нагрузок и собственных частот колебаний составляет необходимую часть общего расчета многоэтажного здания. Как правило, такие расчеты выполняются численно в различных системах, основанных на методе конечных элементов в специализированных программах. Аналитические решения редки, но они требуются для оценки численных решений и расчета упрощенных моделей сооружения. Ставится задача найти аналитическую зависимость прогиба каркаса многоэтажного здания и его первой собственной частоты от числа панелей. Масса модели фермы предполагается равномерно распределенной по узлам. Учитываются только горизонтальные смещения масс. Материалы и методы. С целью расчета усилий в стержнях статически определимой фермы методом вырезания узлов в символьной форме применяется программа, составленная в системе компьютерной математики Maple. Расчет смещения узлов производится с помощью интеграла Мора. Использование метода Донкерлея позволяет найти аналитическую форму нижней оценки первой частоты собственных колебаний для произвольной этажности здания. Анализ серии решений в символьной форме для ферм с последовательно увеличивающимся числом панелей дает возможность вывести итоговую расчетную формулу для частоты колебаний. При определении собственных чисел характеристической матрицы использовался оператор Eigenvalues из пакета LinearAlgebra. Результаты. Полученные формулы для прогиба сооружения (горизонтального смещения) под действием горизонтальной боковой нагрузки и формула для оценки первой частоты колебаний имеют простой вид, удобный для быстрой проверки численного решения и предварительной оценки работоспособности проектируемой конструкции. Выводы. Предложенная модель пространственной статически определимой модели многоэтажного здания дает возможность вывести расчетные формулы для деформаций и собственных частот сооружения при произвольной его этажности. Полученные формулы могут быть использованы в задачах оптимизации и быть основой при расчете более сложной статически неопределимой модели конструкции.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

DEFORMATIONS AND NATURAL VIBRATION FREQUENCY OF A THREE-DIMENSIONAL TRUSS IN A MULTI-STOREY BUILDING

Introduction. The analysis of deformations arising under the action of various loads and natural vibration frequencies is an integral part of the structural analysis of a multi-storey building. As a rule, this analysis is performed numerically in various systems using the finite element method and specialized software programmes. Analytical solutions are rare, but they are needed to evaluate the numerical ones and calculate simplified models of structures. The task is to find the analytical dependence of the frame deflection in a multi-storey building and its first natural frequency on the number of panels. The mass of the truss model is assumed to be distributed uniformly over the nodes. Only horizontal displacements of masses are taken into account. Materials and methods. To calculate the forces arising in the rods of a statically determinate truss using the method of joint isolation in the alphameric form, the software compiled in the Maple math system is applied. The calculation of displacements of nodes is performed using the Mohr integral. The Dunkerley method allows finding the analytical form of the lower estimate of the first frequency of natural vibrations for an arbitrary number of storeys in a building. The analysis of a series of solutions, made in the alphameric form for trusses with a successively growing number of panels allows deriving the final vibration frequency formula. When determining the natural vibrations of the characteristic matrix, the Eigenvalues operator from the LinearAlgebra package was used. Results. The resulting formulas used to analyze the deflection (horizontal displacement) of a structure under the action of a horizontal lateral load and the formula for estimating the first vibration frequency have simple forms, convenient for the quick verification of numerical solutions and preliminary evaluation of the behaviour of a designed structure. Conclusions. The proposed spatial statically determinate model of a multi-storey building allows deriving formulas for the analysis of deformations and natural frequencies of a structure with an arbitrary number of storeys. Resulting formulas can be applied to solve optimization tasks and serve as the basis for calculating a more complex statically indeterminate model of a structure.

Текст научной работы на тему «ДЕФОРМАЦИИ И СОБСТВЕННАЯ ЧАСТОТА КОЛЕБАНИЙ ФЕРМЫ ПРОСТРАНСТВЕННОЙ МОДЕЛИ МНОГОЭТАЖНОГО ЗДАНИЯ»

НАУЧНАЯ СТАТЬЯ / RESEARCH PAPER УДК 624.074

DOI: 10.22227/1997-0935.2022.5.580-588

Деформации и собственная частота колебаний фермы пространственной модели многоэтажного здания

Михаил Николаевич Кирсанов

Национальный исследовательский университет «Московский энергетический институт»

(НИУ «МЭИ»); г. Москва, Россия

АННОТАЦИЯ

Введение. Расчет деформаций под действием различных нагрузок и собственных частот колебаний составляет необходимую часть общего расчета многоэтажного здания. Как правило, такие расчеты выполняются численно в различных системах, основанных на методе конечных элементов в специализированных программах. Аналитические решения редки, но они требуются для оценки численных решений и расчета упрощенных моделей сооружения. Ставится задача найти аналитическую зависимость прогиба каркаса многоэтажного здания и его первой собственной частоты от числа панелей. Масса модели фермы предполагается равномерно распределенной по узлам. Учитываются только горизонтальные смещения масс.

Материалы и методы. С целью расчета усилий в стержнях статически определимой фермы методом вырезания узлов в символьной форме применяется программа, составленная в системе компьютерной математики Maple. Расчет смещения узлов производится с помощью интеграла Мора. Использование метода Донкерлея позволяет найти аналитическую форму нижней оценки первой частоты собственных колебаний для произвольной этажности здания. Анализ серии решений в символьной форме для ферм с последовательно увеличивающимся числом панелей дает возможность вывести итоговую расчетную формулу для частоты колебаний. При определении собственных чисел g g характеристической матрицы использовался оператор Eigenvalues из пакета LinearAlgebra.

еч сч Результаты. Полученные формулы для прогиба сооружения (горизонтального смещения) под действием горизон-

njj' „у тальной боковой нагрузки и формула для оценки первой частоты колебаний имеют простой вид, удобный для быс-

^ Ф трой проверки численного решения и предварительной оценки работоспособности проектируемой конструкции.

Выводы. Предложенная модель пространственной статически определимой модели многоэтажного здания дает возможность вывести расчетные формулы для деформаций и собственных частот сооружения при произвольной его этажности. Полученные формулы могут быть использованы в задачах оптимизации и быть основой при расчете

ш г*

. т- более сложной статически неопределимой модели конструкции.

£ •

2 IE КЛЮЧЕВЫЕ СЛОВА: многоэтажное здание, каркас, ветровая нагрузка, собственные колебания, аналитическое

О — решение, формула Максвелла - Мора

N N

ф ф

Благодарности. Работа выполнена при финансовой поддержке фонда РНФ 22-21-00473.

= 5 ДЛЯ ЦИТИРОВАНИЯ: Кирсанов М.Н. Деформации и собственная частота колебаний фермы пространственной мо-

О ф дели многоэтажного здания // Вестник МГСУ. 2022. Т. 17. Вып. 5. С. 580-588. DOI: 10.22227/1997-0935.2022.5.580-588

—■

о Автор, ответственный за переписку: Михаил Николаевич Кирсанов, C216@ya.ru.

со <Т

Deformations and natural vibration frequency of a three-dimensional truss in a multi-storey building

.E o _

£ " Mikhail N. Kirsanov

m National Research University "Moscow Power Engineering Institute" (MPEI);

o E Moscow, Russian Federation

-

^ ABSTRACT

z £

OT g Introduction. The analysis of deformations arising under the action of various loads and natural vibration frequencies is an in— 2 tegral part of the structural analysis of a multi-storey building. As a rule, this analysis is performed numerically in various sys-

tems using the finite element method and specialized software programmes. Analytical solutions are rare, but they are needed O jj to evaluate the numerical ones and calculate simplified models of structures. The task is to find the analytical dependence of

O the frame deflection in a multi-storey building and its first natural frequency on the number of panels. The mass of the truss

x S model is assumed to be distributed uniformly over the nodes. Only horizontal displacements of masses are taken into account.

S Materials and methods. To calculate the forces arising in the rods of a statically determinate truss using the method of joint

¡E £ isolation in the alphameric form, the software compiled in the Maple math system is applied. The calculation of displace-

jj jj ments of nodes is performed using the Mohr integral. The Dunkerley method allows finding the analytical form of the lower

U > estimate of the first frequency of natural vibrations for an arbitrary number of storeys in a building. The analysis of a series

of solutions, made in the alphameric form for trusses with a successively growing number of panels allows deriving the final

© М.Н. Кирсанов, 2022

Распространяется на основании Creative Commons Attribution Non-Commercial (CC BY-NC)

vibration frequency formula. When determining the natural vibrations of the characteristic matrix, the Eigenvalues operator from the LinearAlgebra package was used.

Results. The resulting formulas used to analyze the deflection (horizontal displacement) of a structure under the action of a horizontal lateral load and the formula for estimating the first vibration frequency have simple forms, convenient for the quick verification of numerical solutions and preliminary evaluation of the behaviour of a designed structure. Conclusions. The proposed spatial statically determinate model of a multi-storey building allows deriving formulas for the analysis of deformations and natural frequencies of a structure with an arbitrary number of storeys. Resulting formulas can be applied to solve optimization tasks and serve as the basis for calculating a more complex statically indeterminate model of a structure.

KEYWORDS: multi-storey building, frame, wind load, natural vibrations, analytical solution, Maxwell - Mohr formula Acknowledgements.

This work was financially supported by the Russian Science Foundation 22-21-00473.

FOR CITATION: Kirsanov M.N. Deformations and natural vibration frequency of a three-dimensional truss in a multi-storey building. Vestnik MGSU [Monthly Journal on Construction and Architecture]. 2022; 17(5):580-588. DOI: 10.22227/19970935.2022.5.580-588 (rus.).

Corresponding author: Mikhail N. Kirsanov, C216@ya.ru.

ВВЕДЕНИЕ

Одним из направлений анализа напряженно-деформированного состояния (НДС) и колебаний многоэтажных зданий является расчет упрощенных модельных представлений каркасов зданий [1, 2]. Расчеты, как правило, производятся в численной форме в различных специализированных программах, основанных на методе конечных элементов [3-6]. В более сложных моделях при расчете динамики сооружения используется волновая теория [7]. Уровни вибрации зданий от движения поездов метрополитена изучены в исследовании [8]. Дискретно-континуальный подход к численному моделированию НДС высотных зданий применялся в работе [9]. Влияние свойств грунтов на состояние высотных зданий рассмотрено в публикации [10]. Альтернативой численным методам в расчете строительных конструкций служат аналитические методы. В трудах [11, 12] в системе символьной математики использовался метод представления решения в виде рядов. Для регулярных систем аналитические решения в виде конечных формул дает метод индукции [13]. Одними из первых, кто поднял проблему существования и расчета регулярных статически определимых стержневых конструкций, были R.G. Hutchinson и N.A. Fleck [14, 15]. Получены формулы для прогиба ряда плоских регулярных ферм с произвольным числом панелей [16-20]. Аналогичные решения для пространственных ферм, как правило, более сложны [21]. В исследовании [22] приведены различные схемы регулярных плоских ферм и формулы для вычисления их прогибов для различных нагрузок в зависимости от числа панелей. Метод индукции применялся также в решениях задач о колебании фермы [21, 23-25].

В настоящей работе рассматривается регулярная статически определимая пространственная ферма многоэтажного здания (рис. 1). Высота этажа h. Два крыла здания по n панелей по горизонтали и k этажей соединены центральной частью, содержащей боковые раскосы длиной Va2 + h2 и Vb2 + h2 . В общем случае крылья здания несимметричные, ширина одного из них — a, другого — b. Ставится

задача найти аналитическую зависимость прогиба каркаса многоэтажного здания от числа панелей при действии боковой равномерно распределенной по узлам ветровой стационарной нагрузки. Цель исследования — индуктивным методом, пользуясь нижней оценкой Донкерлея, вывести формулу для первой частоты собственных колебаний конструкции с произвольным числом панелей.

Число стержней фермы пг = 6^п + т + 2), включая 2(п + т) + 4 опорные стойки и по три опорных раскоса в средней части корпуса на каждом этаже. Число внутренних узлов фермы К = 2^п + т + 2).

МАТЕРИАЛЫ И МЕТОДЫ

Расчет усилий в стержнях. Для расчета усилий в статически определимых фермах не требуется

kh

Рис. 1. Схема здания. Боковая равномерная узловая нагрузка интенсивностью P

Fig. 1. The structural diagram of the building. The lateral uniform nodal load having intensity P

< П

IH

kK

G Г

S 2

0 CO § CO

1 S

У 1

J to

u-

^ I

n °

S 3 o

zs ( о §

E w § 2

n 0

S 6

r 6 t (

S )

ii

® 7 i

. DO

■ T

s □

s У с о <D * Ultt

2 2 О О 10 10 10 10

h

численный метод. Усилия рассчитываются методом вырезания узлов. Для этого составляется система уравнений равновесия всех узлов в проекции на три оси координат. В систему уравнений в качестве неизвестных входят и реакции опор. Матрица системы уравнений состоит из направляющих косинусов усилий, значения которых вычисляются по координатам узлов и порядку соединения стержней в узлах. Стержни и узлы нумеруются (рис. 2).

д = 2и~1)(т + п + 2),

=аЦ-\), У1+д =Ьт,

2!+ч г' = 1> •••> й + 1>

xi+n+Uq = an,

Номера концов стержней записываются в специальные списки Ф. i = 1, ..., п , задающие структуру конструкции. Стойки в здании, например, кодируются так:

Ф+*(4яг + л)+5)+в(]-1) = [ + 0(] " !)' »' + ]

] = 1, ..., к -1, I = 1, ..., 0,

где 0 = 2(т + п) + 4.

Направляющие косинусы имеют вид: gv. =

N N

N N

О О

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

N N

1П 10

К Ш

U 3 > (Л

С И

со

J

- £

ф Ф

о £ —■

о

о <£ со <f CD

8 «

Z ■ ^ от

от IE

Е О

CL ° ^ с

ю о

S ц

о Е

СП ^ т- ^

£

ОТ О

О (П

= 1,

где первый индекс v

= ^Ф, д - Ч2 ' 1 указывает на ось проекции V = х, у, 2; второй индекс — номер стержня; 11 = ^ ^ - уф^ ) —

yi+n+1+д =b(m-i),

Zi

j'+n+l+g

: hj, г = 1,..., т,

i+n+m+1+д

■ a(i-1),

У!+п+т+1+д =b(m + l),

i+n+m+l+q

= kj, 1=1, ..., И+ 2,

i+2n+m+l+q

= а(и +1),

yi+2n+m+3+q =b(m-i +1),

i+2n+m+3+g

= hj, i = 1, ..., 7И + 1.

длина стержня. Элементы матрицы системы уравнений равновесия узлов имеют вид:

^3Ф,Д-2+j,i = gj,i, ^3Ф,_2-2+j,i = ~Sj,i,

j = 1, 2, 3, i = 1, ..., иг,

где параметр j — это номер осей x, y, z соответственно.

Система уравнений равновесия узлов записывается в матричном виде GS = T, где S — вектор, содержащий неизвестные усилия и реакции опор. Вектор нагрузок T имеет длину n Нагрузки, приложенные к узлу i в проекции на ось x, записываются в элементы T3i _ на ось y — в элементы T3i _ Вертикальные внешние силы содержатся в элементах T3i , i = 1, ..., K. Для решения матричного уравнения в символьной форме в системе компьютерной математики Maple используется метод обратной

C

зо x

Рис. 2. Нумерация узлов и стержней фермы, k = 2, n = 3, m = 1 Fig. 2. The numbering of truss nodes and rods, k = 2, n = 3, m = 1

матрицы: S = G-1T. По сравнению со специальными операторами пакета LinearAlgebra этот метод более быстрый.

РЕЗУЛЬТАТЫ

Деформации от ветровой нагрузки. Многоэтажные здания, обладая большой парусностью, могут иметь заметные деформации от боковой ветровой нагрузки. Рассмотрим, например, равномерно распределенную узловую нагрузку интенсивностью Р, приложенную перпендикулярно к одной из стен здания (рис. 1). Вычислим перемещение по направлению действия нагрузки точки С. По формуле Максвелла - Мора получим:

Д = 2>ГS® la/(EF).

а2, = P

9(a2 + с2) + 39h2

Д2,2 = P

А2,З = P

Д 2, 4 = P

Д 2, 5 = P

EFa2 ?

,12(a2 + с2) ■+ -62Й2

EFa2 ?

,15(a2 + с2) + 90h2

EFa2 ?

,18(a2 + с2) + ■123Ä2

EFa2

> 21(a2 + с2) ■+ -161Й2

EFa

\п = Р

EFa2

К=Р

(n + 2)(a3 +с3 +h3)

EFa2

(6я +12) (а3 + с3)+ (13и2+ 71и + 90)/г3

~ЁЫ2 ' (3)

(10и + 20 )(а3+ с3) + (40и2+ 210и + 260)/г3 ЕЕа2 '

(15 п + 30) (а3+ с3) + (95 и2+ 490и + 600)й3 ЕЕа2

Находим общие члены последовательности коэффициентов при подобных слагаемых в выражениях (2), (3). Получим итоговую зависимость смещения точки С в зависимости от размеров фермы, числа этажей и панелей:

А3 , = Р

\п=Р

(1)

A t ,и = P

Сумма составлена по всем деформируемым стержням конструкции, включая опорные стойки и раскосы. Здесь введены обозначения: S¿■P> — усилие от действия внешней нагрузки в стержне с номером а; ^а(1) — усилие в этом же стержне от действия единичной горизонтальной силы, приложенной к узлу С, перемещение которого разыскивается; /а — длина стержня; Е — модуль упругости материала; F — площадь поперечного сечения стержня.

Проанализируем случай симметричной конструкции: т = п. Задача имеет два натуральных параметра: число панелей п и число этажей k. Обозначим Д горизонтальное смещение узла С конструкции с п панелями и k этажами. Фиксируем сначала k и находим прогиб для нескольких ферм с последовательно увеличивающимся числом панелей. При k = 2 в результате расчетов по формуле (1) имеем:

C1 (a3 + с3) + C2 h EFa2

(4)

где С = (n + 2)k(k + 1)/2,

С2 = k(k + 1)(n + 2)(3(n + 3)k2 + (n + 3)k - 4n)/24. (5)

Заметим, что решение не зависит от размера b. Введем величину безразмерного прогиба А' = = EFAk,nj(P0H), отнесенного к высоте фермы H = = kh и суммарной нагрузке P0 = (n + 2)kP. Построим графики зависимости (4) с коэффициентами (5) от числа этажей (рис. 3).

Кривые, построенные для различных размеров а, пересекаются. Начиная с высоты 2 этажа порядок кривых меняется. Для малой этажности перемещение тем больше, чем больше величина а, для зданий выше трех этажей зависимость от этой величины более естественная: чем больше ширина здания а, тем меньше смещение узла C. Кривые имеют квадратичную асимптотику: lim А'/ k2 = h2( n + 3)/(8a2 V k/ / \ /

Оценка низшей частоты собственных колебаний. Инерционные свойства здания моделируем равными сосредоточенными в узлах массами. Число степеней свободы в предположении, что массы смещаются только по горизонтали, равно числу узлов K = 2k(n + m +2). Решить задачу о собственных частотах системы с многими степенями свободы не представляется возможным. Метод Донкер-лея [21, 23-26] позволяет получить нижнюю аналитическую оценку первой собственной частоты.

Приближенная оценка частоты по Донкер-лею юв выражается через парциальные частоты ю :

Здесь c = Va2 + h2 . Пользуясь операторами Maple, обобщаем эти выражения по числу панелей n:

(Зи + 6)(а3+с3) + (5и2/2 + 31и/2 + 21)й3

-.(2)

Аналогично, при других значениях k имеем выражения:

< п

8 8 i Н

k к

G Г

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

0 со

n со

1 2 У 1

J со

^ I

n ° o

=! (

о =?

о n

CO co

0)

(6)

-2 V"1 -2

=Ь ®- •

Р=1

При вычислении парциальных частот необходимо составить отдельные дифференциальные уравнения для масс д:

пу + d у = 0,

Р Р Р

Р = 1, ..., к,

ю со о 2 § Г §6

c я

h о

С n

ф )

ii

® 7 л ' . DO

■ T

s □

s У с о <D *

JJ JJ

2 2 О О 2 2 2 2

a-1

30

25

20

15

10

1

2

3

4

5

6

7

8

сч N

N N

О О

N N

U) 10

к <v

U 3 > 1Л С И

U I»

Рис. 3. Зависимость смещения узла C от числа этажей: I — a = 8 м; II — a = 12 м; III — a = 16 м

Fig. 3. The dependence of the displacement of node C on the number of storeys: I — a = 8 m; II — a = 12 m;

III — a = 16 m

у где ур — ускорение; — коэффициент жесткости; ур — смещение массы в узле. Отсюда частота колебаний отдельного груза имеет вид

щ * ■ (7)

Коэффициент жесткости вычисляется по формуле Максвелла - Мора:

8р = 1ар =£ № ) У(Ер). (8)

J

- £

ф Ф

О ё —■

о

о <£

со <f

CD i;

8 «

Z ■ ^ ОТ

от IE

Е О

CL ° ^ с

ю о

S «

о Е ст> ^

т-

£

ОТ о

£ w

í!

о (ñ

i=i

Из формул (6) и (7) следует юй2 = 3

ю 2 =

p=i

= ц ^ —. Так как податливость есть величина, об-

p=i d

ратная жесткости 1/d = 5, то

г p p'

р=1 ct=l

Последовательно вычислим суммы Ък = h2

^=1 а=1

ZZ((^) ^а при n = m = 1, a = b, k = 1, 2, 3, ...

=(iW2 +16 )а3 + 30c3 + 30A3, S2 = (28л/2 + 62) а3 + 90c3 + 266A3,

^з =

+138) + 180c3 + 1060A3, S4 = (56^2 + 244) а3 + 300c3 + 2940A3, S5 + 380) а3 + 450c3 + 6610A3,...

В общем случае £k = C1 аъ + C2С + C3hi, где коэффициенты в этом выражении определяются как общие члены коэффициентов в выражении (10). Методами Maple с использованием операторов составления и решения рекуррентных уравнений получим:

C = к (V2 +1 + 15k), C2 = 15k (к +1), C3 = к (к +1) (22k2 + 22к +1)).

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

В результате имеем нижнюю оценку для первой ч астоты по Донкерлею:

(üd = h

EF

ц {С1аъ + C2c3 +

c3h у

=Ф, =^ТЕ(^Р])\/(ЕР)= ^/{h2EF). (9)

Для оценки степени приближения этого решения вычислим численно весь спектр частот конструкции и сравним частоту юв с наименьшей частотой спектра.

Динамика системы с К степенями свободы описывается системой дифференциальных уравнений движения грузов в матричном виде:

MKY + DKY = 0,

(11)

(10)

где MK — матрица инерции размером К*К; У — вектор ускорений; БК — матрица жесткости; У — вектор вертикальных перемещений масс 1, ..., К. В случае одинаковых масс матрица инерции выражается через единичную матрицу MK = д1К. Матрица податливости ВК, являющаяся обратной к матрице

5

k

жесткости DK, определяется по формуле Максвелла - Мора:

bt,j = XS° S(j lJ(EF),

(12)

где S (а — усилие в стержне а от действия единичной вертикальной силы в узле I. Задача сводится к проблеме собственных чисел матрицы Вк. Умножая слева матричное уравнение (11) на Вк с учетом замены У = -ю2У, соответствующей гармоническим колебаниям z. = и. sm(юt + ф0), получим В^ = = где 1 = 1/(дю2) — собственное число матрицы В; ю — собственная частота колебаний. Отсюда

с» = ^1/ (|А). Для определения собственного числа матрицы можно использовать оператор Eigenvalues из пакета LinearAlgebra системы Maple.

Рассмотрим пример здания с k этажами с размерами h = 3,25 м, a = 6 м, д = 600 кг. Жесткость стальных стержней фермы примем EF = 1,62 х х 106 кН. На рис. 4 построены зависимости оценки юс Донкерлея и первой частоты собственных колебаний системы с K степенями свободы, полученной численно. Заметно некоторое повышение точности оценки с увеличением этажности здания.

Найдем зависимость относительной погрешности оценки Донкерлея от числа панелей по высоте сооружения и высоты каждой панели (этажа) h. Введем обозначение для величины погрешности е = |юв - Mjl/fflj. На рис. 5 три кривые, рассчитанные для разных высот, показывают, что наибольшая точность расчетов достигается для высоких зданий с большим числом этажей.

70

60

50

40

30

20

10

\

\

ш

ш D

-

0,54 0,52 0,50 0,48 0,46 0,44 0,42 0,40

\ I

II

III

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 1718

Рис. 4. Сравнение аналитической нижней оценки первой собственной частоты по методу Донкерлея с первой частотой спектра частот, полученного численно Fig. 4. Using the Dunkerley method to compare the analytical lower estimate of the first natural frequency with the numerically derived first frequency of the spectrum

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18

Рис. 5. Относительная погрешность оценки Донкерлея: I — h = 3,25 м; II — h = 3,5 м; III — h = 3,75 м Fig. 5. The relative error of the Dunkerley estimation: I — h = 3.25 m; II — h = 3.5 m; III — h = 3.75 m

ЗАКЛЮЧЕНИЕ

Построена простая математическая модель стального каркаса многоэтажного здания в виде статически определимой фермы. Принятые упрощения реальной статически неопределимой конструкции, содержащей не только стержни, но и плиты, позволили получить простые расчетные формулы для величины горизонтального смещения угловой точки здания на верхнем этаже в зависимости от числа панелей при действии на здание ветровой горизонтальной нагрузки. Эта задача решена в два этапа для двух натуральных параметров регулярности, которые определяют число этажей k и число панелей n в каждом из двух отсеков. Сначала проводилась индукция по одному параметру, в результате которой получалась некоторая промежуточная общая формула, затем серия этих формул обобщалась по другому параметру. Для проверки решения очередность параметров менялась. Решена и другая задача для предложенной конструкции — расчет низшей собственной частоты в приближении Донкерлея. Точность полученного решения оказалась не очень высокой, но с увеличением высоты h и этажности k здания точность росла. Это особенно ценно потому, что каркасные модели высоких зданий содержат большое число стержней и в численном решении в этом случае может сказываться эффект накопления ошибок округления. Разработанный алгоритм может быть использован и для других каркасных сооружений как основа для более сложных расчетов, например, методом сил.

< п 18

i Н

k к

G Г

o

n со

1 2

У 1

J со I

n °

2 3 o

zs ( n

co co

Q)

l\J CO О

2§ > §

c я

h о

cn

CD )

ii

® 7 л '

. DO ■

s □

s у с о <D *

JJ JJ 22 О О 10 10 10 10

ш, c

СПИСОК ИСТОЧНИКОВ

1. Shatilov Y.Y., Lyapin A.A. Vibration-based damage detection techniques for health monitoring of construction of a multi-storey building // Materials Science Forum. 2018. Vol. 931. Pp. 178-183. DOI: 10.4028/www.scien-tific.net/msf.931.178

2. Rafezy B., Zare A., Howson W.P. Coupled lateral-torsional frequencies of asymmetric, three-dimensional frame structures // International Journal of Solids and Structures. 2007. Vol. 44. Issue 1. Pp. 128-144. DOI: 10.1016/j.ijsolstr.2006.04.019

3. Alam Z., Sun L., Zhang C., Su Z., Samali B. Experimental and numerical investigation on the complex behaviour of the localised seismic response in a multi-storey plan-asymmetric structure // Structure and Infrastructure Engineering. 2021. Vol. 17. Issue 1. Pp. 86-102. DOI: 10.1080/15732479.2020.1730914

4. Wang Y., ThambiratnamD.P., Chan T.H.T., Nguyen A. Method development of damage detection in asymmetric buildings // Journal of Sound and Vibration. 2018. Vol. 413. Pp. 41-56. DOI: 10.1016/j.jsv.2017.10.015

5. Юзвенко А.Ю. Анализ возможностей программного комплекса Autodesk Robot Structural

Jj су Analysis Professional для динамического расчета мно-° ° гоэтажного здания // Вестник Гродненского государ-uj uf ственного университета имени Янки Купалы. Серия 6. g ® Техника. 2019. Т. 9. № 1. С. 40-49. j? ¡я 6. Малахова А., Мухин М. Проектирование желе-

2 зобетонных конструкций с использованием программ. т- ного комплекса ЛИРА : учебное пособие. М. : НИУ £ g МГСУ, 2016. 119 с.

7. Назаров Ю.П., Позняк Е.В., Филимонов А.В. Н ¡2 Анализ вида волновой модели и получение расчетных ^ параметров сейсмического воздействия для высотно-Л го здания // Промышленное и гражданское строитель-о| ство.2014.№ 5.С. 40-45.

§ о 8. Дашевский М.А., Мондрус В.Л. Прогноз уров-

со ^ ней вибрации зданий от движения поездов метрополи-

§ § тена // Промышленное и гражданское строительство.

2013.№11.С. 52-54.

от 9. Акимов П.А. О развитии дискретно-контину-

от Е

— -{д ального подхода к численному моделированию состо-

.Е о яния несущих систем высотных зданий // Промышлен-

6Ъ с ное и гражданское строительство. 2015. № 3. С. 16-20. со 10. Тер-Мартиросян З.Г., Тер-Мартиросян А.З.

с? | Исследования грунтов оснований высотных зданий //

g ° Основания, фундаменты и механика грунтов. 2009.

£ № 5. с. 2-12. от g 11. GoloskokovD.P.,MatrosovA.V. A superposition

T ^ method in the analysis of an isotropic rectangle // Applied

>J Э Mathematical Sciences. 2016. Vol. 10. Pp. 2647-2660.

| g DOI: 10.12988/ams.2016.67211

^ S 12. Goloskokov D.P., Matrosov A.V. Comparison

I of two analytical approaches to the analysis of grillages //

¡3 In 2015 International Conference "Stability and Control

£ £ Processes" in Memory of V.I. Zubov (SCP). 2015. DOI: 10.1109/SCP.2015.7342169

13. Бука-ВайвадеК., КирсановМ.Н., Сердюк Д.О. Calculation of deformations of a cantilever-frame planar truss model with an arbitrary number of panels // Вестник МГСУ. 2020. Т. 15. № 4. С. 510-517. DOI: 10.22227/1997-0935.2020.4.510-517

14. HutchinsonR.G., FleckN.A. Microarchitecture cellular solids — the hunt for statically determinate periodic trusses // ZAMM. 2005. Vol. 85. Issue 9. Pp. 607-617. DOI: 10.1002/zamm.200410208

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

15. Hutchinson R.G., FleckN.A. The structural performance of the periodic truss // Journal of the Mechanics and Physics of Solids. 2006. Vol. 54. Issue 4. Pp. 756-782. DOI: 10.1016/j.jmps.2005.10.008

16. Dai Qiao. Analytical dependence of planar truss deformations on the number of panels // AlfaBuild. 2021. Vol. 17. Issue 2. P. 1701. DOI: 10.34910/ALF.17.1

17. Комерзан Е.В., Свириденко О.В. Аналитический расчет прогиба плоской внешне статически неопределимой фермы с произвольным числом панелей // Строительная механика и конструкции. 2021. № 2 (29). С. 29-37.

18. Ovsyannikova V.M. Dependence of deformations of a trapezous truss beam on the number of panels // Structural Mechanics and Structures. 2020. № 3 (26). Pp. 13-20.

19. Овсянникова В.М. Зависимость прогиба плоской внешне статически неопределимой фермы от числа панелей // Строительная механика и конструкции. 2020. № 4 (27). С. 16-25.

20. Суд И.Б. Вывод формул для прогиба шпрен-гельной балочной фермы с произвольным числом панелей в системе Maple // Строительная механика и конструкции. 2020. № 2 (25). С. 25-32.

21. Kirsanov M., Vorobyev O. Calculating of a spatial cantilever truss natural vibration frequency with an arbitrary number of panels: analytical solution // Construction of Unique Buildings and Structures. 2021. Vol. 94. Issue 1. P. 9402. DOI: 10.4123/CUBS.94.2

22. Kirsanov M.N. Planar Trusses: Schemes and Formulas. UK : Cambridge Scholars Publishing, 2019. 206 p.

23. Petrenko V. The natural frequency of a two-span truss // AlfaBuild. 2021. Vol. 20 Issue 5 (20). P. 2001. DOI: 10.34910/ALF.20.1

24. Петренко В.Ф. Оценка собственной частоты двухпролетной фермы с учетом жесткости опор // Строительная механика и конструкции. 2021. № 4 (31). С. 16-25. DOI: 10.36622/VSTU.2021.31.4.002

25. Петриченко Е.А. Нижняя граница частоты собственных колебаний фермы Финка // Строительная механика и конструкции. 2020. № 3 (26). С. 21-29.

26. Levy C. An iterative technique based on the Dun-kerley method for determining the natural frequencies of vibrating systems // Journal of Sound and Vibration. 1991. Vol. 150. Issue 1. Pp. 111-118. DOI: 10.1016/0022-460X(91)90405-9

Поступила в редакцию 30 марта 2022 г. Принята в доработанном виде 5 мая 2022 г. Одобрена для публикации 5 мая 2022 г.

Об авторе: Михаил Николаевич Кирсанов — доктор физико-математических наук, профессор кафедры робототехники, мехатроники, динамики и прочности машин; Национальный исследовательский университет «Московский энергетический институт» (НИУ «МЭИ»); 111250, г. Москва, ул. Красноказарменная, д. 14; SPIN-код: 8679-6853, Scopus: 16412815600, ResearcherID: H-9967-2013, ORCID: 0000-0002-8588-3871, Google Scholar: FfoNGFwAAAAJ, IstinaResearcherID: 2939132; C216@ya.ru.

REFERENCES

1. Shatilov Y.Y., Lyapin A.A. Vibration-based damage detection techniques for health monitoring of construction of a multi-storey building. Materials Science Forum. 2018; 931:178-183. DOI: 10.4028/www. scientific.net/msf.931.178

2. Rafezy B., Zare A., Howson W.P. Coupled la-teral-torsional frequencies of asymmetric, three-dimensional frame structures. International Journal of Solids and Structures. 2007; 44(1):128-144. DOI: 10.1016/j. ijsolstr.2006.04.019

3. Alam Z., Sun L., Zhang C., Su Z., Samali B. Experimental and numerical investigation on the complex behaviour of the localised seismic response in a multi-storey plan-asymmetric structure. Structure and Infrastructure Engineering. 2021; 17(1):86-102. DOI: 10.1080/15732479.2020.1730914

4. Wang Y., Thambiratnam D.P., Chan T.H.T., Nguyen A. Method development of damage detection in asymmetric buildings. Journal of Sound and Vibration. 2018; 413:41-56. DOI: 10.1016/j.jsv.2017.10.015

5. Yuzvenko A.Yu. Analysis of the capabilities of autodesk robot structural analysis professional software package for dynamic calculation of a multistorey building. Vesnik of Yanka Kupala State University of Grodno. Series 6. Engineering Science. 2019; 9(1):40-49. (rus.).

6. Malakhova A., Mukhin M. Design of reinforced concrete structures using the LIRA software package: textbook. Moscow, NRU MGSU, 2016; 119. (rus.).

7. Nazarov Y.P., Poznyak E.V., Filimonov A.V. Wave model analysis and calculation of seismic impact parameters for high-rise building. Industrial and Civil Engineering. 2014; 5:40-45. (rus.).

8. Dashevskij M.A., Mondrus V.L. Prognosis of vibration levels of buildings from underground train traffic. Industrial and Civil Engineering. 2013; 11:5254. (rus.).

9. Akimov P.A. About development of discrete-continual approach to numerical modeling of load-bearing structures of high-rise buildings. Industrial and Civil Engineering. 2015; 3:16-20. (rus.).

10. Ter-Martirosyan Z.G., Ter-Martirosyan A.Z. Soil beds of high-rise buildings. Soil Mechanics and Foundation Engineering. 2009; 5:2-12. (rus.).

11. Goloskokov D.P., Matrosov A.V. A superposition method in the analysis of an isotropic rectangle. Applied Mathematical Sciences. 2016; 10:2647-2660. DOI: 10.12988/ams.2016.67211

12. Goloskokov D.P., Matrosov A.V. Comparison of two analytical approaches to the analysis of grillages. 2015 International Conference "Stability and Control Processes " in Memory of V.I. Zubov (SCP). 2015. DOI: 10.1109/SCP.2015.7342169

13. Buka-Vaivade K., Kirsanov M.N., Ser-djuk D.O. Calculation of deformations of a cantilever-frame planar truss model with an arbitrary number of panels. Vestnik MGSU [Monthly Journal on Construction and Architecture]. 2020; 15(4):510-517. DOI: 10.22227/1997-0935.2020.4.510-517

14. Hutchinson R.G., Fleck N.A. Microarchitec-tured cellular solids — the hunt for statically determinate periodic trusses. ZAMM. 2005; 85(9):607-617. DOI: 10.1002/zamm.200410208

15. Hutchinson R.G., Fleck N.A. The structural performance of the periodic truss. Journal of the Mechanics and Physics of Solids. 2006; 54(4):756-782. DOI: 10.1016/jjmps.2005.10.008

16. Dai Qiao. Analytical dependence of planar truss deformations on the number of panels. AlfaBuild. 2021; 17(2):1701. DOI: 10.34910/ALF.17.1

17. Komerzan E.V., Sviridenko O.V. Analytical calculation of the deflection of a plane external statically undeterminated truss with an arbitrary number of panels. Structural Mechanics and Structures. 2021; 2(29):29-37. (rus.).

18. Ovsyannikova V.M. Dependence of deformations of a trapezous truss beam on the number of panels. Structural Mechanics and Structures. 2020; 3(26):13-20. (rus.).

19. Ovsyannikova V.M. Dependence of the deflection of a planar external statically undeterminable truss on the number of panels. Structural Mechanics and Structures. 2020; 4(27):16-25. (rus.).

20. Sud I.B. Derivation of formulas for deflection of the girder truss with an arbitrary number of panels in the Maple system. Structural Mechanics and Structures. 2020; 2(25):25-32. (rus.).

21. Kirsanov M., Vorobyev O. Calculating of a spatial cantilever truss natural vibration frequency

< n 88 i Н

k к

G Г

S 2

o n

1 2

У 1

J to

u-I

n

2 3 o

=! ( n

u

n 2

n 0

2 § > §

t (

ф )

ii

i7 i

. DO

■ T

s □

(Л у с о ii

JJJJ 22 о о 10 10 10 10

with an arbitrary number of panels: analytical solution.

Construction of Unique Buildings and Structures. 2021; 94(1):9402. DOI: 10.4123/CUBS.94.2

22. Kirsanov M.N. Planar Trusses: Schemes and Formulas. UK, Cambridge Scholars Publishing, 2019; 206.

23. Petrenko V. The natural frequency of a two-span truss. AlfaBuild. 2021; 20(5):(20):2001. DOI: 10.34910/ALF.20.1

24. Petrenko V.F. Estimation of the natural frequency of a two-span truss, taking into account the sup-

Received March 30, 2022.

Adopted in revised form on May 5, 2022.

Approved for publication on May 5, 2022.

port stiffeness. Structural Mechanics and Structures. 2021; 4(31):16-25. DOI: 10.36622/VSTU.2021.31.4.002 (rus.).

25. Petrichenko E.A. Lower bound of the natural oscillation frequency of the fink truss. Structural Mechanics and Structures. 2020; 3(26):21-29. (rus.).

26. Levy C. An iterative technique based on the Dunkerley method for determining the natural frequencies of vibrating systems. Journal of Sound and Vibration. 1991; 150(1):111-118. DOI: 10.1016/0022-460X(91)90405-9

Bionotes: Mikhail N. Kirsanov — Doctor of Physical and Mathematical Sciences, Professor of the Department of Robotics, Mechatronics, Dynamics and Strength of Machines; National Research University "Moscow Power Engineering Institute" (MPEI); 14 Krasnokazarmennaya st., Moscow, 111250, Russian Federation; SPIN-code: 8679-6853, Scopus: 16412815600, ResearcherlD: H-9967-2013, ORCID: 0000-0002-8588-3871, Google Scholar: FfoNGFwAAAAJ, IstinaResearcherlD: 2939132; C216@ya.ru.

N N

N N

О О

tv N

in in

К (V

U 3

> (Л

с и

to I»

i

<D <u

о ё

о

о о СО <

8 « §

(Л "

со E

E О

CL °

^ с

ю о

S «

о E

со ^

T- ^

CO CO

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

■s

iE 3s

О (П

i Надоели баннеры? Вы всегда можете отключить рекламу.