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

АНАЛИТИЧЕСКАЯ ОЦЕНКА ЧАСТОТЫ СОБСТВЕННЫХ КОЛЕБАНИЙ ФЕРМЫ С ПРОИЗВОЛЬНЫМ ЧИСЛОМ ПАНЕЛЕЙ Текст научной статьи по специальности «Физика»

CC BY
176
18
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
БАЛОЧНАЯ ФЕРМА / СОБСТВЕННЫЕ КОЛЕБАНИЯ / НИЖНЯЯ ОЦЕНКА ЧАСТОТЫ / ВЕРХНЯЯ ОЦЕНКА ЧАСТОТЫ / МЕТОД ДОНКЕРЛЕЯ / ИНДУКЦИЯ / MAPLE / BEAM TRUSS / NATURAL OSCILLATIONS / LOWER FREQUENCY ESTIMATE / UPPER FREQUENCY ESTIMATE / DUNKERLEY METHOD / INDUCTION

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

Цель исследования - вывод формулы зависимости первой частоты собственных колебаний плоской статически определимой балочной фермы с параллельными поясами от числа панелей, размеров и одинаковых масс, сосредоточенных в узлах нижнего пояса фермы. Решетка фермы треугольная с вертикальными стойками. В решении использованы операторы системы компьютерной математики Maple. Методы. Основой для верхней оценки искомой частоты колебаний регулярной фермы является энергетический метод. В качестве формы прогиба фермы взят прогиб от действия равномерно распределенной нагрузки. Предполагаются только вертикальные перемещения грузов. Амплитудные значения прогиба фермы вычисляются по формуле Максвелла - Мора. Усилия в стержнях определяются в символьной форме методом вырезания узлов. Зависимость решения от числа панелей получается индуктивным обобщением серии решений для ферм с последовательно увеличивающимся числом панелей. Для последовательностей коэффициентов искомой формулы составляются и решаются однородные линейные рекуррентные уравнения четвертого порядка. Результаты. Решение сравнивается с численным решением, полученным из анализа всего спектра собственных частот колебаний системы масс, расположенных в узлах фермы. Частотное уравнение составляется и решается с помощью операторов поиска собственных значений в системе Maple. Показано, что полученная аналитическая оценка отличается от численного решения на доли процента. При этом с увеличением числа панелей погрешность энергетического метода монотонно уменьшается. Приведена более простая нижняя оценка частоты колебаний по методу Донкерлея. Точность оценки снизу значительно меньше оценки сверху, зависит от размеров и числа панелей.

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

ANALYTICAL ASSESSMENT OF THE FREQUENCY OF NATURAL VIBRATIONS OF A TRUSS WITH AN ARBITRARY NUMBER OF PANELS

The aim of the work is to derive a formula for the dependence of the first frequency of the natural oscillations of a planar statically determinate beam truss with parallel belts on the number of panels, sizes and masses concentrated in the nodes of the lower truss belt. Truss has a triangular lattice with vertical racks. The solution uses Maple computer math system operators. Methods. The basis for the upper estimate of the desired oscillation frequency of a regular truss is the energy method. As a form of deflection of the truss taken deflection from the action of a uniformly distributed load. Only vertical mass movements are assumed. The amplitude values of the deflection of the truss is calculated by the Maxwell - Mohr’s formula. The forces in the rods are determined in symbolic form by the method of cutting nodes. The dependence of the solution on the number of panels is obtained by an inductive generalization of a series of solutions for trusses with a successively increasing number of panels. For sequences of coefficients of the desired formula, fourth-order homogeneous linear recurrence equations are compiled and solved. Results. The solution is compared with the numerical one, obtained from the analysis of the entire spectrum of natural frequencies of oscillations of the mass system located at the nodes of the truss. The frequency equation is compiled and solved using Eigenvalue search operators in the Maple system. It is shown that the obtained analytical estimate differs from the numerical solution by a fraction of a percent. Moreover, with an increase in the number of panels, the error of the energy method decreases monotonically. A simpler lower bound for the oscillation frequency according to the Dunkerley method is presented. The accuracy of the lower estimate is much lower than the upper estimate, depending on the size and number of panels.

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

2020. 16(5). 351-360 Строительная механика инженерных конструкций и сооружений Structural Mechanics of Engineering Constructions and Buildings

HTTP://JOURNALS.RUDN.RU/STRUCTURAL-MECHANICS

DOI 10.22363/1815-5235-2020-16-5-351-360 УДК 624.04:531.391.3

НАУЧНАЯ СТАТЬЯ / RESEARCH PAPER

История статьи

Поступила в редакцию: 3 июля 2020 г. Доработана: 11 сентября 2020 г. Принята к публикации: 19 сентября 2020 г.

Аналитическая оценка частоты собственных колебаний фермы с произвольным числом панелей

М.Н. Кирсанов

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

Красноказарменная, 14

c216@ya.ru

Аннотация

Цель исследования - вывод формулы зависимости первой частоты собственных колебаний плоской статически определимой балочной фермы с параллельными поясами от числа панелей, размеров и одинаковых масс, сосредоточенных в узлах нижнего пояса фермы. Решетка фермы треугольная с вертикальными стойками. В решении использованы операторы системы компьютерной математики Maple. Методы. Основой для верхней оценки искомой частоты колебаний регулярной фермы является энергетический метод. В качестве формы прогиба фермы взят прогиб от действия равномерно распределенной нагрузки. Предполагаются только вертикальные перемещения грузов. Амплитудные значения прогиба фермы вычисляются по формуле Максвелла - Мора. Усилия в стержнях определяются в символьной форме методом вырезания узлов. Зависимость решения от числа панелей получается индуктивным обобщением серии решений для ферм с последовательно увеличивающимся числом панелей. Для последовательностей коэффициентов искомой формулы составляются и решаются однородные линейные рекуррентные уравнения четвертого порядка. Результаты. Решение сравнивается с численным решением, полученным из анализа всего спектра собственных частот колебаний системы масс, расположенных в узлах фермы. Частотное уравнение составляется и решается с помощью операторов поиска собственных значений в системе Maple. Показано, что полученная аналитическая оценка отличается от численного решения на доли процента. При этом с увеличением числа панелей погрешность энергетического метода монотонно уменьшается. Приведена более простая нижняя оценка частоты колебаний по методу Донкерлея. Точность оценки снизу значительно меньше оценки сверху, зависит от размеров и числа панелей.

Ключевые слова: балочная ферма, собственные колебания, нижняя оценка частоты, верхняя оценка частоты, метод Донкерлея, Maple, индукция

Для цитирования

Кирсанов М.Н. Аналитическая оценка частоты собственных колебаний фермы с произвольным числом панелей // Строительная механика инженерных конструкций и сооружений. 2020. Т. 16. № 5. С. 351-360. http://dx.doi.org/10.22363/1815-5235-2020-16-5-351-360

Кирсанов Михаил Николаевич, доктор физико-математических наук, профессор кафедры робототехники, мехатроники, динамики и прочности машин Института энергомашиностроения и механики; eLIBRARY SPIN-код: 8679-6853, ORCID iD: https://orcid.org/0000-0002-8588-3871, ResearcherlD: H-9967-2013, IstinaResearcherlD: 2939132.

© Кирсанов M.H.. 2020

li^v 0 I This work is licensed under a Creative Commons Attribution 4.0 International License https://creativec0mm0ns.0rg/licenses/by/4.Q/

Article history

Received: July 3, 2020 Revised: September 11, 2020 Accepted: September 19, 2020

Analytical assessment of the frequency of natural vibrations of a truss with an arbitrary number of panels

Mikhail N. Kirsanov

National Research University "Moscow Power Engineering Institute ", 14 Krasnokazarmennaya St, Moscow, 111250, Russian Federation c216@ya.ru

Abstract

The aim of the work is to derive a formula for the dependence of the first frequency of the natural oscillations of a planar statically determinate beam truss with parallel belts on the number of panels, sizes and masses concentrated in the nodes of the lower truss belt. Truss has a triangular lattice with vertical racks. The solution uses Maple computer math system operators. Methods. The basis for the upper estimate of the desired oscillation frequency of a regular truss is the energy method. As a form of deflection of the truss taken deflection from the action of a uniformly distributed load. Only vertical mass movements are assumed. The amplitude values of the deflection of the truss is calculated by the Maxwell - Mohr's formula. The forces in the rods are determined in symbolic form by the method of cutting nodes. The dependence of the solution on the number of panels is obtained by an inductive generalization of a series of solutions for trusses with a successively increasing number of panels. For sequences of coefficients of the desired formula, fourth-order homogeneous linear recurrence equations are compiled and solved. Results. The solution is compared with the numerical one, obtained from the analysis of the entire spectrum of natural frequencies of oscillations of the mass system located at the nodes of the truss. The frequency equation is compiled and solved using Eigenvalue search operators in the Maple system. It is shown that the obtained analytical estimate differs from the numerical solution by a fraction of a percent. Moreover, with an increase in the number of panels, the error of the energy method decreases monotonically. A simpler lower bound for the oscillation frequency according to the Dunkerley method is presented. The accuracy of the lower estimate is much lower than the upper estimate, depending on the size and number of panels.

Keywords: beam truss, natural oscillations, lower frequency estimate, upper frequency estimate, Dunkerley method, Maple, induction

For citation

Kirsanov M.N. Analytical assessment of the frequency of natural vibrations of a truss with an arbitrary number of panels. Structural Mechanics of Engineering Constructions and Buildings. 2020;16(5):351-360. (In Russ.) http://dx.doi.org/10.22363/1815-5235-2020-16-5-351-360

1. Введение

Наряду с жесткостью и прочностью конструкции значение первой частоты колебаний фермы является одной из ее важнейших эксплуатационных характеристик. Эта величина входит в динамические расчеты вынужденных и затухающих колебаний фермы, в решения задач вибрации и сейсмостойкости. Существуют численные расчеты спектров колебаний ферм, реализуемые в различных специализированных пакетах [1-6]. Несмотря на широкое распространение численных методов расчета, их высокую точность и большой выбор доступных моделей ферм (линейных и нелинейных) для анализа частот колебаний, аналитические методы, сводящиеся к сравнительно простым формулам, редки.

С появлением в арсенале инженера современных систем компьютерной математики (Maple, Mathe-matica, Reduce и др.) аналитические методы становятся все более распространенными [7-12]. Формулы, полученные в результате решения уравнений в символьном виде, удобны как для оценки точности решений в специализированных пакетах, так и для простого предварительного расчета частотных характеристик сооружения. Наиболее значимы решения применимые для широкого класса систем. Расширение класса рассчитываемых аналитически систем в регулярных фермах достигается введением в набор параметров формулы числа панелей. Возможности аналитических решений для регулярных статически определимых ферм (плоских и пространственных) рассматривались в трудах R.G. Hutchinson, N.A. Fleck,

Mikhail N. Kirsanov, Doctor of Physical and Mathematical Sciences, Professor of the Department of Robotics, Mechatronics, Dynamics and Machine Strength of the Institute of Power Machinery and Mechanics; eLIBRARY SPIN-code: 8679-6853, ORCID iD: https://orcid.org/0000-0002-8588-3871, ResearcherID:H-9967-2013, IstinaResearcherID: 2939132.

F.W. Zok, R.M. Latture, M.R. Begley [13-15]. В работах [16; 17] предложены алгоритмы расчета динамики в общем случае статически неопределимых стержневых конструкций в аналитической форме с применением методов строительной механики. В настоящем исследовании ставится цель получить аналитическое решение задачи о низшей частоте собственных колебаний балочной фермы (рис. 1) в виде формулы, включающей в себя не только размеры и инертные свойства фермы, но и число панелей.

2. Методы

Энергетический метод. Масса фермы с треугольной решеткой и стойками моделируется одинаковыми грузами в узлах (шарнирах) нижнего пояса. Рассматриваются только вертикальные колебания грузов. Число степеней свободы системы грузов фермы с n панелями равно N = 2n - 1. Длина каждой панели 2a, высота h. Рассматриваемая ферма относится к числу регулярных статически определимых конструкций, для которых доступны аналитические методы нахождения зависимостей усилий и деформаций от числа ячеек регулярности (здесь - панелей).

аааааааааа

I-1-1-1-1-1-1-1-1-1-

Рис. 1. Схема фермы, n = 5 [Figure 1. Scheme of the truss, n = 5]

Система дифференциальных уравнений динамики грузов имеет матричный вид:

mnY+DnY = 0, (1)

где Y = [yi,y2, ...,yN~\T - вектор вертикальных перемещений масс 1,..., N; DN - матрица жесткости; MN -матрица инерции размером N х N; Y - вектор ускорений.

В случае одинаковых масс грузов матрица инерции имеет диагональный вид MN = mIN . Элементы матрицы податливости BN, являющейся обратной к матрице жесткости DN , определяются по формуле Максвелла - Мора:

=S )S" )/а /(EF), (2)

а=1

где EF - жесткость стержней; S^) - усилие в стержне а от действия единичной вертикальной силы в узле i; la - длина стержня а ; ц = 8n - число стержней фермы вместе с тремя стержнями, моделирующими подвижную и неподвижную опоры.

Опорные стержни принимаются недеформируемыми, суммирование в (2) на эти стержни не распространяется. Усилия в стержнях фермы определяются одновременно с реакциями опор из решения системы уравнений равновесия всех узлов фермы. Матрица системы строится по значениям направляющих косинусов усилий, найденных по данным координат концов стержней. Решение ищется в символьной форме с использованием возможностей компьютерной математики Maple [18; 19\.

Умножая (1) слева на BN с учетом замены Y = —o2Y, соответствующей колебаниям вида

Ук = uk sin(o + ФоХ (3)

сводим задачу к проблеме собственных чисел матрицы BN : BNY = XY, где X = 1/ (дао2) - собственное число матрицы BN ; о - собственная частота колебаний.

Усилия S^) в стержнях фермы, входящие в элементы матрицы BN , определяются из решения системы уравнений узлов фермы, в которую входят также и реакции опор. Для этого в системе компьютерной математики Maple составляется матрица уравнений равновесия узлов, состоящая из направляющих косинусов усилий в стержнях.

Получить аналитические решения для фермы с произвольным числом панелей не удается. Возможностей системы компьютерной математики Maple хватает только для анализа ферм с n =1,2,...,6. Например, для n = 2 оператором Eigenvalues получены следующие характеристические числа:

X = (a3 + с3 + 2h3)/(2h2 EF ), X2,3 = (2с3 + h3 + 6a3 ± (2с6 + 16a3с3 + h6 + 32a6)(1/2))/(2h2EF).

Соответствующие частоты вычисляются по формуле ю = ф/(тХ). Эти частоты включены также

во все спектры частот ферм рассматриваемого вида с четным числом панелей, что подтверждает свойство вложения спектров регулярных систем [20; 21].

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

Tmax max * (4)

Кинетическая энергия масс в узлах нижнего пояса:

N

Tmax = maX S mVl /2

k=1

Согласно (3) вертикальная скорость массы k имеет вид

vk = У к =юик sin(rat + Фо). Отсюда, полагая max(sin(ra t + ф0)) = 1, получаем

N

Tmax =ю2 S mu2k/2, (5)

к=1

где амплитуда вертикального смещения вычисляется по формуле Максвелла - Мора:

Uk = S Sf)sS^X / (EF) =PS Sf) S^x / (EF) = Puk.

a=1 a=1

Введены обозначения S^ - усилие в стержне a = 1,...,ц от действия нагрузки P, распределенной по узлам нижнего пояса, и Sak) - усилие в этом же стержне от единичной (безразмерной) нагрузки, приложенной к массе с номером k, S^ ) = Sa(P) / P. Выбор такой нагрузки определяется близостью изогнутой оси нижнего пояса форме колебаний системы грузов с первой частотой. Таким образом, (5) приобретает вид

Tmax = PW S ^/2, (6)

k=1

где йк = ик / Р = ^ ))1а / (ЕР) - амплитуда смещений массы с номером к под действием распреде-

а=1

ленной нагрузки (рис. 2), отнесенная к величине P.

Рис. 2. Равномерно распределенная нагрузка, n = 6 [Figure 2. Uniformly distributed load, n = 6\

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

П max = £ Sa^) Ala / 2 =£ ^ la / (2EF).

a=1 a=1

N

В силу линейности задачи S^) = P£ S^) . Отсюда

k=1

ц N N ц N

Птах = P2£S^OP)£S^% /(2EF) = P2££SOp)S^% / (2EF) = P2£uk. (7)

a=1 k=1 k=1 a k=1

Из (4) - (6) следует верхняя оценка первой частоты колебаний фермы (формула Релея):

N N

°R =Е uk/ S muk• (8)

k=1 k=1

Для получения требуемой зависимости частоты от числа панелей перемещения uk должны быть получены так же, как функции n. План решения следующий: 1) расчет перемещения массы с номером 1 при различных значениях n с последующим определением общего члена последовательностей u1 (n) ; 2) расчет перемещений масс 2, 3, 4, ...; 3) обобщение формул u1(n), u2(n), u3(n),... по номеру массы и получение искомой зависимости uk (n) .

Выведем формулу для смещения первой массы от действия распределенной нагрузки.

Расчет перемещения для ферм с различным числом панелей показывает, что вид решения не зависит от n:

uk = £ Cr3/(h2 EF ),

i=1

где для сокращения записи введены обозначения размеров фермы r1 = a, r2 = c = Va2 + h2, r3 = h. Коэффициенты Ci, i = 1,2,3 зависят только от n. Для выявления закономерности в последовательности решений потребовалось рассчитать прогибы ферм с числом панелей n = 1,...,8. Для значений коэффициента C1 получена следующая последовательность: 5, 35/2, 42, 165/2, 143, 455/2, 340, 969/2. Для ее элементов оператор rgf_findrecur системы Maple возвращает рекуррентное уравнение Cn 1 = 4Cn—11 — 6Cn—21 + 4Cn—31 — Cn—41. Оператор rsolve с использованием начальных значений 5, 35/2, 42, 165/2 дает решение этого уравнения: C1 = n(4n2 —1)/6. Число рассчитываемых ферм n = 1,...,8 определяется возможностью получить рекуррентное уравнение. Если длина последовательности недостаточно велика, оператор rgf_findrecur не дает решения. Аналогично при других значениях номеров масс k имеем следующие выражения:

k = 2, C1 = (n — 1)(4n2 + 4n — 3) / 3, k = 3, C1 = (2n — 3)(2n2 + 3n — 4)/2,

k = 4, Cj = 2(n-2)(4n2 + 8n-15)/3, k = 5, C1 = 5(2n - 5)(2n2 + 5n -12)/6,

Обобщение этих формул по k с применением операторов Maple дает окончательное выражение зависимости коэффициента C1 от номера массы и числа панелей:

C (k) = k(k - 2n)(k2 - 2nk -1 - 4n2) /12. (9)

Также в два этапа, обобщением сначала по n, затем по k, получаются и другие коэффициенты:

C2(k) = -k(k-2n)/2, C3(k) = (1 -(-1)k)/2. (10)

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

Отдельно вычислим числитель и знаменатель в (8). Числитель представим в виде

N 3

Z Uk = Z g fJ(h2EF), (11)

k=1 a=1

где g =ZC(k), i = 1,2,3.

k=1

Конечные суммы коэффициентов (9) - (0) дает оператор sum системы Maple:

g = (16n4 -1)/ 30, g2 = (4n2 -1)/6, g3 = 1. (12)

Знаменатель (8):

N 3 3

S тй1 = mSS fa pr > 3/(Ä 4 E 2 F2), (13)

к=1 a=1 ß=1

N

где обозначено Zap = Z Ca (k)CP (k), ^ P = 1, 2,3•

k=1

Посредством суммирования оператором sum получаем следующие коэффициенты: fu = n(4n2 - 1)(496n6 + 328n4 + 103n2 +18)/11340,

f22 = n(16^ - 1)/60, f33 = n,

f12 = n(4n2 - 1)(68n4 + 31n2 + 6) /1260,

f13 = n(8n4 + 5n2 + 2)/30, f23 = n(2n2 +1)/6. (14)

Таким образом, верхнюю оценку первой частоты фермы в зависимости от числа панелей можно получить по формуле

h2 EF S g,

3

r3

a a

= 3 3 a=1--(15)

R 3 3

, ,3r3

'aß' a' р

a=1 ß=1

mSS fap rar ;

с коэффициентами (12), (14).

Оценка Донкерлея. Нижнюю оценку первой частоты колебаний получим по формуле Донкерлея:

N

ш-2

(16)

D

к=1

гдешк - частота колебания одной массы m , расположенной в узле k + 1 нижнего пояса, принимая нумерацию узлов от левой опоры.

Уравнение (1) в случае колебаний одной массы будет скалярным:

шУк + dkУk = 0

где yk - вертикальное перемещение массы; ук - вектор ускорений; dк - коэффициент жесткости (к -номер массы).

Частота колебаний груза ш* = ■^й^Гш . Коэффициент жесткости, обратный коэффициенту податливости, определяется по формуле Максвелла - Мора:

5 * = 1/йк = £() ) /в/( ЕЕ),

а=1

где ¿О*) - усилия в стержне с номером а от действия единичной вертикальной силы, приложенной к узлу к + 1 (рис. 3).

I-1-1-1-1-!-1-1-1-1-1-h

Рис. 3. К определению коэффициента жесткости, n = 6, k = 2 [Figure 3. On the determination of stiffness coefficient, n = 6, k = 2\

Методом индукции, последовательно рассчитывая смещение сначала массы 1 от действия нагрузки на этот узел при различных числах панелей n = 1,2,3,... и обобщая эти решения на произвольное число n, затем рассчитывая смещение других масс, прикладывая всякий раз единичную силу к месту положения массы, получаем решение для произвольных значений n и k:

8 k =Е Cr3/(h2 EF),

i=1

где

C1 = k (2k2 — 4nk — 1)(k — 2n) / (6n), C2 = — k (k — 2и)/(2и),

C3 = (1 — (—1)k )/2.

Из (16) суммированием этих выражений по k получаем нижнюю оценку для первой частоты по Донкерлею:

N

o—D = m£8k = m((32n4 + 20n2 — 7)a3 + 15(4n2 — 1)c3 + 90h3n)/(90h2 EF). (17)

k=1

3. Результаты

Оценить точность полученных решений можно графически. Рассмотрим для примера ферму с n панелями высотой h = 5 м, длиной панели a = 3 м и грузами m = 100 кг в узлах нижнего пояса. Жесткость

8 "ТТ

стальных стержней фермы примем EF = 0,2 -10 Н. Первая частота Or собственных колебаний фермы (15),

полученная по энергетическому методу Рэлея, и оценка Юр Донкерлея (17) с увеличением числа панелей сближаются (рис. 4). Численное же значение Ш1 собственной частоты системы с N = 2п — 1 степенями

свободы, найденная как минимальная частота полного спектра частот, почти полностью совпадает с оценкой Релея, кривые полностью сливаются. Для уточненной оценки степени полученного приближения введем относительные погрешности =| Юр —Ш1 | /Ш1, гк =| —Ш1 | /ш^

СО,

r-1

80- \ \

\ \ \ \

V \

\ Л \ Y

a

10- \

9 10 11 12

Рис. 4. Первая частота колебаний, полученная тремя способами в зависимости от числа панелей [Figure 4. The first oscillation frequency obtained by three methods depending on the number of panels]

Рис. 5. Погрешность оценки Донкерлея в зависимости от числа панелей: 1 - h = 3 м; 2 - h = 4 м; 3 - h = 5 м [Figure 5. Dunkerley's estimation error depending on the number of panels: 1 - h = 3 m; 2 - h = 4 m; 3 - h = 5 m]

0.002-

Рис. 6. Погрешность оценки Рэлея в зависимость от числа панелей:

1 - h = 3 м; 2 - h = 4 м; 3 - h = 5 м [Figure 6. Rayleigh's estimation error depending on the number of panels: 1 - h = 3 m; 2 - h = 4 m; 3 - h = 5 m]

В зависимости от числа панелей погрешность решения Донкерлея меняется от 4 % при большом числе панелей до 22 % при п = 2 (рис. 5). Графики построены для тех же параметров фермы, что и рис. 4. Для фермы с меньшей высотой степень приближения оценки Донкерлея больше.

Оценка Рэлея, так же как и оценка Донкерлея, лучше при большем числе панелей (рис. 6), за исключением начального участка кривой, имеющей резкий всплеск при п = 3. Главная же особенность оценки, полученной по Рэлею, - весьма высокая ее точность. Точность порядка 0,2 % сопоставима с точ-

ностью входных данных: размеров, упругих характеристик и масс. На рис. 3 кривая Ш1, найденная численно из анализа спектра системы со многими степенями свободы, практически сливается с кривой зависимости от числа панелей частоты по Рэлею Or .

4. Заключение

Для решения поставленной задачи в работе использована система компьютерной математики Maple. Выбор этой программы определялся только предпочтением автора. С не меньшим успехом здесь можно использовать и другие аналогичные программы, например Mathematica или Derive. В основном операторы Maple использовались для составления и решения в символьной форме уравнений равновесия узлов фермы. Определение же общих членов последовательностей можно выполнить в каком-либо онлайн-сервисе, например [22\. По сравнению с решениями задач о прогибе ферм с произвольным числом панелей, полученными методом индукции, формулы для оценок первой частоты оказались заметно сложнее. Однако они имеют замкнутую форму, не содержат сумм и рядов, которые надо вычислять, обрывая суммирование на некотором слагаемом, что вносит элемент субъективности в решение. Полученные формулы также не связаны с применением специальных функций и итеративными вычислениями. Выведенные формулы можно использовать при довольно большом числе панелей, то есть именно в тех случаях, где наиболее вероятно накопление ошибок вычислений. Кроме того, замкнутый аналитический вид решений позволяет применять все средства математического анализа для выявления их особенностей.

Список литературы

1. Ufimtcev E. Dynamic Calculation of Nonlinear Oscillations of Flat Trusses. Part 2. Examples of Calculations // Procedia Engineering. 2017. Vol. 206. Pp. 850-856. DOI: 10.1016/jproeng.2017.10.561.

2. Tejani G.G., Savsani V.J., Patel V.K., Mirjalili S. Truss optimization with natural frequency bounds using improved symbiotic organisms search // Knowledge-Based Systems. 2018. Vol. 143. Pp. 162-178. DOI: 10.1016/j.knosys.2017.12.012.

3. Ufimtsev E., Voronina M. Research of Total Mechanical Energy of Steel Roof Truss during Structurally Nonlinear Oscillations // Procedia Engineering. 2016. Vol. 150. Pp. 1891-1897. DOI: 10.1016/j.proeng.2016.07.188

4. Jalbi S., Bhattacharya S. Closed form solution for the first natural frequency of offshore wind turbine jackets supported on multiple foundations incorporating soil-structure interaction // Soil Dynamics and Earthquake Engineering. 2018. Vol. 113 (May). Pp. 593-613. https://doi.org/10.1016/j.soildyn.2018.06.011

5. Kilikevicius A., Fursenko A., Jurevicius M., Kilikeviciene K., Bureika G. Analysis of parameters of railway bridge vibration caused by moving rail vehicles // Measurement and Control (United Kingdom). 2019. Vol. 52. No. 9-10. Pp. 1210-1219. DOI: 10.1177/0020294019836123.

6. Алдушкин Р.В., Савин С.Ю. Исследование работы треугольных ферм при статических и динамических воздействиях // Строительство и реконструкция. 2010. № 3 (29). С. 3-6. https://elibrary.ru/download/elibrary_ 15503961_71987283.pdf (дата обращения: 01.05.2020).

7. Bolotina T.D. The deflection of the flat arch truss with a triangular lattice depending on the number of panels // Вестник научных конференций. 2016. № 4-3 (8). С. 7-8.

8. Тимофеева Т.А. Формулы для расчета прогиба плоской решетчатой рамы с произвольным числом панелей // Строительная механика и конструкции. 2019. № 4 (23). С. 26-33.

9. Бойко А.Ю., Ткачук Г.Н. Вывод формул зависимости прогиба плоской шарнирно-стержневой рамы от числа панелей в системе Maple // Строительная механика и конструкции. 2019. № 4 (23). С. 15-25.

10. Белянкин Н.А., Бойко А.Ю. Формулы для прогиба балочной фермы с произвольным числом панелей при равномерном загружении // Строительная механика и конструкции. 2019. № 1 (20). С. 21-29.

11. Ткачук Г.Н. Формула зависимости прогиба несимметрично нагруженной плоской фермы с усиленными раскосами от числа панелей // Строительная механика и конструкции. 2019. № 2 (21). С. 32-39.

12. Тиньков Д.В. Сравнительный анализ аналитических решений задачи о прогибе ферменных конструкций // Инженерно-строительный журнал. 2015. № 5 (57). С. 66-73.

13. Hutchinson R.G., FleckN.A. Microarchitectured cellular solids - the hunt for statically determinate periodic trusses // ZAMM Z. Angew. Math. Mech. 2005. Vol. 85. No. 9. Pp. 607-617.

14. Hutchinson R.G., Fleck N.A. The structural performance of the periodic truss // Journal of the Mechanics and Physics of Solids. 2006. Vol. 54. No. 4. Pp. 756-782.

15. Zok F.W., Latture R.M., Begley M.R. Periodic truss structures // Journal of the Mechanics and Physics of Solids. 2016. Vol. 96. Pp. 184-203.

16. Рыбаков Л.С., Мишустин И.В. Собственные колебания плоских регулярных упругих ферм ортогональной структуры // Механика композиционных материалов и конструкций. 1999. Т. 5. № 2. С. 3-16.

17. Мишустин И.В., Рыбаков Л.С. Колебания плоских упругих ферм ортогональной структуры // Известия Академии наук. Механика твердого тела. 2003. № 2. С. 168-184.

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

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

19. Kirsanov M. Analytical Solution of a Spacer Beam Truss Deflection with an Arbitrary Number of Panels // Construction of Unique Buildings and Structures. 2020. Vol. 88 Article No. 8802. DOI: 10.18720/CUBS.88.2.

20. Кирсанов М.Н., Тиньков Д.В. Аналитические выражения частот малых колебаний балочной фермы с произвольным числом панелей // Строительная механика и конструкции. 2019. Т. 1. № 20. С. 14-20.

21. Тиньков Д.В. Аналитические решения задач о собственных частотах колебаний регулярных стержневых систем: дис. ... к.т.н. М., 2019. 113 с.

22. WolframAlpha System. URL: https://www.wolframalpha.com/examples/mathematics/ (accessed: 03.07.2020).

References

1. Ufimtcev E. Dynamic Calculation of Nonlinear Oscillations of Flat Trusses. Part 2. Examples of Calculations. ProcediaEngineering. 2017;206:850-856. DOI: 10.1016/j.proeng.2017.10.561.

2. Tejani G.G., Savsani V.J., Patel V.K., Miijalili S. Truss optimization with natural frequency bounds using improved symbiotic organisms search. Knowledge-Based Systems. 2018;143:162-178. DOI: 10.1016/j.knosys.2017.12.012.

3. Ufimtsev E., Voronina M. Research of Total Mechanical Energy of Steel Roof Truss during Structurally Nonlinear Oscillations. Procedia Engineering. 2016;150:1891-1897. DOI: 10.1016/j.proeng.2016.07.188.

4. Jalbi S., Bhattacharya S. Closed form solution for the first natural frequency of offshore wind turbine jackets supported on multiple foundations incorporating soil-structure interaction. Soil Dynamics and Earthquake Engineering. 2018;113(May):593-613. https://doi.org/10.1016/j.soildyn.2018.06.011

5. Kilikevicius A., Fursenko A., Jurevicius M., Kilikeviciene K., Bureika G. Analysis of parameters of railway bridge vibration caused by moving rail vehicles. Measurement and Control (United Kingdom). 2019;52(9-10):1210-1219. DOI: 10.1177/0020294019836123.

6. Aldushkin R.V., Savin S.Yu. Investigation of the operation of triangular farms with static and dynamic effects. Construction and reconstruction. 2010;3(29):3-6. (In Russ.) Available from: https://elibrary.ru/download/elibrary_ 15503961_71987283.pdf (accessed: 01.05.2020).

7. Bolotina T.D. The deflection of the flat arch truss with a triangular lattice depending on the number of panels. Bulletin of Scientific Conferences. 2016;(8)4-3:7-8.

8. Timofeeva Т.А. Formulas for calculating the deflection of a flat lattice frame with an arbitrary number of panels. Structural mechanics and structures. 2019;4(23):26-33. (In Russ.)

9. Boyko A.Y., Tkachuk G.N. Derivation of the formulas for the deflection of a flat hinged-rod frame in the of symbol mathematics Maple system. Structural mechanics and structures. 2019;4(23):15-25. (In Russ.)

10. Belyankin N.A., Boyko A.Y. Formula for deflection of a girder with an arbitrary number of panels under the uniform load. Structural mechanics and structures. 2019;1(20):21-29. (In Russ.)

11. Tkachuk G.N. The formula for the dependence of the deflection of an asymmetrically loaded flat truss with reinforced braces on the number of panels. Structural mechanics and structures. 2019;2(21):32-39. (In Russ.)

12. Tinkov D.V. Comparative analysis of analytical solutions to the problem of deflection of truss structures. Magazine of Civil Engineering. 2015;5(57):66-73. (In Russ.)

13. Hutchinson R.G., Fleck N.A. Microarchitectured cellular solids - the hunt for statically determinate periodic trusses. ZAMMZ. Angew. Math. Mech. 2005;9(85):607-617.

14. 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.

15. Zok F.W., Latture R.M., Begley M.R. Periodic truss structures. Journal of the Mechanics and Physics of Solids. 2016;96:184-203.

16. Rybakov L.S., Mishustin I.V. Natural vibrations of flat regular elastic trusses of orthogonal structure. Mechanics of composite materials and structures. 1999;2(5):3-16. (In Russ.)

17. Mishustin I.V., Rybakov L.S. Oscillations of flat elastic trusses of orthogonal structure. News of the Academy of Sciences. Solid Mechanics. 2003;2:168-184. (In Russ.)

18. Buka-Vaivade K., Kirsanov M.N., Serdjuks D.O. Calculation of deformations of a cantilever frame planar truss model with an arbitrary number of panels. VestnikMGSU. 2020;15(4):510-517.

19. Kirsanov M. Analytical Solution of a Spacer Beam Truss Deflection with an Arbitrary Number of Panels. Construction of Unique Buildings and Structures. 2020;88:8802.

20. Kirsanov M.N., Tinkov D.V. Analytical expressions of the frequencies of small vibrations of a beam truss with an arbitrary number of panels. Structural Mechanics and Structures. 2019;20(1):14-20. (In Russ.)

21. Tinkov D.V. Analytical solutions to problems on natural frequencies of oscillations of regular rod systems (Thesis of Candidate of Technical Sciences). Moscow; 2019. (In Russ.)

22. WolframAlpha System. Available from: https://www.wolframalpha.com/examples/mathematics/ (accessed: 03.07.2020).

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