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

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

CC BY
36
9
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
БАЛКА / УПРУГИЕ ОПОРЫ / КОЛЕБАНИЯ / УСТОЙЧИВОСТЬ / ПЕРВАЯ СОБСТВЕННАЯ ЧАСТОТА / ПЕРВАЯ КРИТИЧЕСКАЯ СИЛА / АНАЛИТИЧЕСКОЕ РЕШЕНИЕ

Аннотация научной статьи по физике, автор научной работы — Рабецкая О.И., Кудрявцев И.В., Митяев А.Е.

В работе рассматривается проблема обеспечения требуемой первой собственной частоты изгибных колебаний балки при действии продольной силы за счет введения необходимой жесткости опор. Рассматривая и объединяя уравнения свободных колебаний балки и уравнения, описывающие потерю ее устойчивости, было получено условие работоспособности на основе обеспечения минимально заданного значения первой собственной частоты колебаний с учетом действия продольной силы. При этом достижение нулевой частоты собственных колебаний соответствует потере устойчивости, что позволяет решать обе задачи. Данная задача математически сложна и в известной научной литературе ее решение обычно приводится только в графическом или табличном виде. Проблема заключается в нелинейной зависимости коэффициентов опор от жесткости при колебаниях и потере устойчивости. Для решения этой проблемы использовалась аппроксимация нелинейных зависимостей коэффициентов опор методом наименьших квадратов и получения квадратичных аппроксимирующих функций. В результате задача определения требуемой жесткости опор свелась к разрешающему алгебраическому уравнению четвертой степени, для которого существует аналитическое решение. Полученное решение позволяет определить жесткость опор балки, которая обеспечивает требуемое значение первой собственной частоты колебаний балки и ее первой критической нагрузки в виде внешней сжимающей силы или температурных воздействий. Замена нелинейных зависимостей коэффициентов опор от жесткости опор более простыми квадратичными функциями привела к относительно простым аналитическим зависимостям, которые позволяют преобразовывать разрешающее уравнение в соответствии с конкретной решаемой задачей. Вместе с тем, квадратичные функции повлияли на погрешность расчета, для снижения которой было произведено ограничение рассматриваемого диапазона жесткостей опор и разбиение его на три зоны. Проведено сравнение результатов расчетов по предложенному аналитическому решению с численными расчетами методом конечных элементов. Сравнение результатов расчета показало погрешность не более 5 % для рассматриваемого диапазона жесткостей опор, что вполне достаточно для инженерных расчетов балочных конструкций. Для ограничения погрешности результата рекомендуется, чтобы жесткости обоих опор были равны или же одного порядка.

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

Похожие темы научных работ по физике , автор научной работы — Рабецкая О.И., Кудрявцев И.В., Митяев А.Е.

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

BEAM SUPPORT STIFFNESS ANALYTIC SOLUTION FOR THE FIRST EIGENFREQUENCY AND CRITICAL FORCE

The work discusses the problem of providing the required first natural frequency of bending vibrations of the beam under the action of a longitudinal force by introducing the necessary stiffness of the supports. Considering and combining the equations of free vibrations of the beam and the equations describing the loss of its stability, the operability condition was obtained because of providing a minimum given value of the first natural frequency of vibrations considering the action of the axial force. In this case, the achievement of the zero frequency of natural vibration corresponds to the loss of stability, which allows solving both problems. This problem is mathematically complex, and in the known scientific literature its solution is usually given only in graphical or tabular form. The problem lies in the nonlinear dependence of the coefficients of supports on the stiffness during vibrations and loss of stability. To solve this problem, the approximation of the nonlinear coefficients of the supports by the least squares method and the obtaining of quadratic approximating functions was used. As a result, the problem of determining the required stiffness of the supports brought to a fourth-degree resolving algebraic equation, for which an analytic solution exists. The obtained solution allows the stiffness of the beam supports, which provides the required value of the first natural frequency of vibrations of the beam and its first critical load in the form of external compressive force or temperature effects. Replacing the nonlinear dependencies of the support coefficients with the stiffness of the supports with simpler quadratic functions led to relatively simple analytic dependencies that allow the resolution equation to be transformed according to the particular problem being solved. At the same time, quadratic functions influenced the calculation error, to reduce which, the range of the support stiffness under consideration was limited and divided into three zones. The results of calculations using the proposed analytical solution were compared with numerical calculations using finite element method. The comparison of the calculation results showed an error of not more than 5 % for the considered range of stiffness of the supports, which is quite enough for engineering calculations of beam structures. To limit the error of the result, it is recommended that the stiffnesses of both supports be equal or of the same order.

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

УДК 534.11 + 539.3

Doi: 10.31772/2712-8970-2022-23-4-708-720

Для цитирования: Рабецкая О. И., Кудрявцев И. В., Митяев А. Е. Аналитический расчет жесткости опор балки для обеспечения первой собственной частоты колебаний и критической силы // Сибирский аэрокосмический журнал. 2022. Т. 23, № 4. С. 708-720. Doi: 10.31772/2712-8970-2022-23-4-708-720.

For citation: Rabetskaya O. I., Kudryavtsev I. V., Mityaev A. E. [Beam support stiffness analytic solution for the first eigenfrequency and critical force]. Siberian Aerospace Journal. 2022, Vol. 23, No. 4, P. 708-720. Doi: 10.31772/2712-8970-2022-23-4-708-720.

Аналитический расчет жесткости опор балки для обеспечения первой собственной частоты колебаний и критической силы

1* 2 2 О. И. Рабецкая , И. В. Кудрявцев , А. Е. Митяев

1Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева Российская Федерация, 660037, г. Красноярск, просп. им. газ. «Красноярский Рабочий», 31

2Сибирский федеральный университет Российская Федерация, 660041, г. Красноярск, просп. Свободный, 79 E-mail: olga_rabez@mail.ru

В работе рассматривается проблема обеспечения требуемой первой собственной частоты из-гибных колебаний балки при действии продольной силы за счет введения необходимой жесткости опор. Рассматривая и объединяя уравнения свободных колебаний балки и уравнения, описывающие потерю ее устойчивости, было получено условие работоспособности на основе обеспечения минимально заданного значения первой собственной частоты колебаний с учетом действия продольной силы. При этом достижение нулевой частоты собственных колебаний соответствует потере устойчивости, что позволяет решать обе задачи. Данная задача математически сложна и в известной научной литературе ее решение обычно приводится только в графическом или табличном виде. Проблема заключается в нелинейной зависимости коэффициентов опор от жесткости при колебаниях и потере устойчивости. Для решения этой проблемы использовалась аппроксимация нелинейных зависимостей коэффициентов опор методом наименьших квадратов и получения квадратичных аппроксимирующих функций. В результате задача определения требуемой жесткости опор свелась к разрешающему алгебраическому уравнению четвертой степени, для которого существует аналитическое решение. Полученное решение позволяет определить жесткость опор балки, которая обеспечивает требуемое значение первой собственной частоты колебаний балки и ее первой критической нагрузки в виде внешней сжимающей силы или температурных воздействий. Замена нелинейных зависимостей коэффициентов опор от жесткости опор более простыми квадратичными функциями привела к относительно простым аналитическим зависимостям, которые позволяют преобразовывать разрешающее уравнение в соответствии с конкретной решаемой задачей. Вместе с тем, квадратичные функции повлияли на погрешность расчета, для снижения которой было произведено ограничение рассматриваемого диапазона жесткостей опор и разбиение его на три зоны. Проведено сравнение результатов расчетов по предложенному аналитическому решению с численными расчетами методом конечных элементов. Сравнение результатов расчета показало погрешность не более 5 % для рассматриваемого диапазона жесткостей опор, что вполне достаточно для инженерных расчетов балочных конструкций. Для ограничения погрешности результата рекомендуется, чтобы жесткости обоих опор были равны или же одного порядка.

Ключевые слова: балка, упругие опоры, колебания, устойчивость, первая собственная частота, первая критическая сила, аналитическое решение.

Beam support stiffness analytic solution for the first eigenfrequency

and critical force

O. I. Rabetskaya1*, I. V. Kudryavtsev2, A. E. Mityaev2

1Reshetnev Siberian State University of Science and Technology 31, Krasnoyarskii Rabochii prospekt, Krasnoyarsk, 660037, Russian Federation 2Siberian Federal University 79, Svobodniy prospekt, Krasnoyarsk, Russian Federation, 660041 *E-mail: olga_rabez@mail.ru

The work discusses the problem of providing the required first natural frequency of bending vibrations of the beam under the action of a longitudinal force by introducing the necessary stiffness of the supports. Considering and combining the equations of free vibrations of the beam and the equations describing the loss of its stability, the operability condition was obtained because ofproviding a minimum given value of the first natural frequency of vibrations considering the action of the axial force. In this case, the achievement of the zero frequency of natural vibration corresponds to the loss of stability, which allows solving both problems. This problem is mathematically complex, and in the known scientific literature its solution is usually given only in graphical or tabular form. The problem lies in the nonlinear dependence of the coefficients of supports on the stiffness during vibrations and loss of stability. To solve this problem, the approximation of the nonlinear coefficients of the supports by the least squares method and the obtaining of quadratic approximating functions was used. As a result, the problem of determining the required stiffness of the supports brought to a fourth-degree resolving algebraic equation, for which an analytic solution exists. The obtained solution allows the stiffness of the beam supports, which provides the required value of the first natural frequency of vibrations of the beam and its first critical load in the form of external compressive force or temperature effects. Replacing the nonlinear dependencies of the support coefficients with the stiffness of the supports with simpler quadratic functions led to relatively simple analytic dependencies that allow the resolution equation to be transformed according to the particular problem being solved. At the same time, quadratic functions influenced the calculation error, to reduce which, the range of the support stiffness under consideration was limited and divided into three zones. The results of calculations using the proposed analytical solution were compared with numerical calculations using finite element method. The comparison of the calculation results showed an error of not more than 5 % for the considered range of stiffness of the supports, which is quite enough for engineering calculations of beam structures. To limit the error of the result, it is recommended that the stiffnesses of both supports be equal or of the same order.

Keywords: beam, elastic supports, vibrations, stability, first eigenfrequency, first critical force, analytical solution.

Введение

Протяженные балочные конструкции, работающие в условиях вынужденных колебаний и подвергающиеся действию сжимающих сил, должны соответствовать определенным условиям работоспособности [1; 2]. Одними из таких условий работоспособности являются обеспечение значений первой собственной частоты колебаний и первой критической силы:

f >[ f ], Pcrl >[P], (1)

где f ] и [Р] - допускаемые значения частоты и силы.

Постановка задачи колебаний и устойчивости балки сводится к дифференциальным уравнениям с заданными граничными условиями, которые определяются условиями закрепления балки. Для простых условий закреплений решение такой задачи тривиально и его можно найти в научной литературе и справочниках по динамике [3-34]. Однако в существующей литературе опоры имеют идеализированные модели: шарниры, жесткая заделка и др. В динамике балок принято

учитывать опорные закрепления в виде так называемых коэффициентов опор, которые получаются из решения дифференциальных уравнений динамического состояния и для простых случаев закрепления, когда жесткость равна нулю или же бесконечности, их значения приведены в научной литературе. Промежуточные значения жесткостей опор и соответствующие им значения коэффициентов опор упоминаются редко и обычно в виде графиков или таблиц [35-40]. В действительности же, опоры практически всегда будут обладать некоторой конечной жесткостью, которая будет влиять как на значение первой частоты колебаний, так и на первую критическую силу. Учет жесткости опор усложняет решение задачи динамики балки, поскольку приводит к трансцендентному уравнению, для которого уже нет явного аналитического решения.

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

1. Уравнения динамического поведения балки

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

1.1. Получение разрешающего уравнения

Пусть прямолинейная балка, закрепленная в шарнирных опорах с некоторой жесткостью k2, подвергается внешнему воздействию в виде продольной силы Р (рис. 1). Динамическое состояние балки будем оценивать первой частотой собственных изгибных колебаний и первой критической силой.

Рис. 1. Расчетная схема балки с упругими опорами Fig. 1. Spring-hinged beam

Уравнение свободных изгибных колебаний балки с учетом действия сжимающей продольной силы Р имеет вид [3-26]:

д4 y д 2 y д 2 y

EJmin ^Г + P -TT + m -тг = °> (2)

dx дх dt

где y = y(x, t) - функция прогиба; E - модуль Юнга материала; Jmin - минимальный момент инерции поперечного сечения; m - удельная масса, m = pS, где р - плотность материала; S -площадь поперечного сечения; Р - продольная сила.

Продольная сила Р, при достижении ей некоторого критического значения, может привести к потере устойчивости [27-34]. Уравнение устойчивости балки имеет вид:

EJmn dXT + PtT = °. (3)

дх дх

Решение уравнения свободных колебаний (2) для первой частоты собственных колебаний при Р = 0 имеет вид

/ - ^ ^ (4)

где а - коэффициент опор для первой собственной частоты колебаний, который учитывает влияние способа закрепления балки.

Решение задачи устойчивости (3) для первой критической силы запишется как

п2 Е/ ■

р - 71 Е/Ш1П (5)

сг1 2 ?2 ' ^ '

Ц2 • I

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

Задача устойчивости связана со свободными колебаниями продольной силой Р, поскольку действие этой силы будет изменять весь спектр частот балки. Данное влияние описывается уравнением Галефа [41; 42]:

Л(р<>0) - /1(р-0) 1" — , (6)

где /1(Р <> 0) - первая частота свободных колебаний при действии продольной силы Р; /1(Р = 0) - первая частота свободных колебаний при отсутствии продольной силы Р; Рсг1 - первая критическая сила.

Учтем возможное воздействие температуры на балку (рис. 1), для этого разложим продольную силу на две компоненты, температурную и силовую:

Р - РАТ + РР-а,•АТ • ЕБ + РР, (7)

где РАТ - продольная сила от изменения температуры; Рр - внешнее силовое воздействие; а, -коэффициент температурного расширения материала балки.

Подставляя (4), (5), (7) в уравнение (6), получим уравнение для первой частоты собственных колебаний балки, учитывающее температурное и силовое воздействия и способ закрепления в виде

/ 1(р<>0) -I ~т) 'л/—|п2Е/шт-Ц2 • 12 •(а, -АТ• ЕБ + Рр)|. (8)

Если выражение (8) будет равно нулю, то это свидетельствует о динамической потере устойчивости и позволяет использовать данное выражения сразу для двух целей: обеспечение частоты собственных колебаний и критической силы. Значения коэффициентов опор а и ц для распространенных случаев закреплений и требуемой частоты собственных колебаний или критической силы можно найти в различных справочниках по динамике [35-40] или рассчитать из уравнений динамического состояния. В данной работе будем рассматривать только первую собственную частоту колебаний и первую критическую силу.

1.2. Учет влияния жесткости опор

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

а-а(¿1,к2) , ц-ц(кх,^) . (9)

Выразим коэффициенты опор через относительные жесткости С

С=к Е—, с2 = ^ Е—. (10)

Е/- Е/-

Ш1П Ш1П

Тогда для расчета необходимо определить функции:

« = «(1 С2 ) Ц^^ С2 ) . (11)

В справочниках функции (11) обычно приводятся в табличном или графическом виде, поскольку данные зависимости сильно нелинейные и не могут быть выражены простой аналитической зависимостью. Для решения этой проблемы в работе [43] были получены аналитические функции (11) для ограниченных диапазонов изменения жесткостей Сь С2. Используем эти результаты для решения задачи определения жесткости опор балки и обеспечения заданных значений ее первой собственной частоты колебаний и первой критической силы.

2. Определение требуемой жесткости опор

Жесткость опор, обеспечивающая первую частоту колебаний и первую критическую силу, является неявной переменной в уравнении (8), для определения которой необходимо подставить в аналитическом виде функции (11) и разработать методику решения полученного разрешающего уравнения.

2.1. Получение разрешающего уравнения

Преобразуем уравнение (8) таким образом, чтобы из него можно было выразить искомую жесткость опор и разработать метод ее нахождения:

\2

^С) • ^ [У / - Ц2 (С, С2 ) • /2 • (а, -АГ • ЕБ + РР )] ] (12)

где/ - требуемое значение первой собственной частоты колебаний. Преобразуем уравнение (12) к виду

С

А-ц2СС2) = 4, ), (13)

а (C1, С2)

где

А = л2Е/ш1п; В = /2^ (а, ■ АТ • ЕБ + РР ); С = 4т[/!]2 л4/4. (14)

Решение уравнения (13) относительно жесткостей С^ требует наличия не самих функций коэффициентов опор (11), а их 2-й и 4-й степеней. Данные функции получим методом наименьших квадратов [44; 45]. Для этого на основе исходных функций (11) создадим набор данных для ограниченного диапазона жесткостей С = 0 - 1000, возведем их в соответствующую степень и вновь аппроксимируем полученные результаты квадратичными функциями вида

ц2 (С,, С2 ) = « •( С2 + С22 ) + «2 •(С! + С2 ) + «3, (15)

(С,,С2) = ь •(с2 + С22) + Ь2 •(( + С2) + Ь3. (16)

а4

С целью повышения точности аппроксимации, разделим рассматриваемый диапазон жесткостей на три зоны. Полученные значения коэффициентов в уравнениях (15), (16) приведены в табл. 1.

Уравнение (13) имеет неопределенность решения в виде бесконечного сочетания жесткостей С\ и С2, которые удовлетворяют этой задаче. Для устранения данной определенности примем С2 = п • Сь тогда вместо (15), (16) получим

ц2 (С1) = « • (( + п2 С2 ) + а2 •(С1 + п • С) + а3, (17)

а

(с1 ) = Ь1 •(с2 + п2 • с2) + Ь2 •(С1 + п• С1) + Ь3.

(18)

Таблица 1

Значения коэффициентов аппроксимирующих функций

Зона Коэф( жциенты функции ц2 Коэффициенты функции а4

а! а2 а3 Ь1 Ь2 Ь3

I: С = 0 - 10 0,0053 -0,0855 1 -0,725 17,29 97,4

II: С = 10 - 100 1,21Е-05 -1,82Е-03 0,386 -1,84Е-02 2,945 243

III: С = 100 - 1000 1,3Е-08 -1,93Е-05 0,2664 -4,5Е-05 0,0675 451,5

Для краткости записи примем обозначения:

°1п = а1 •(1+п2);

а2 п = а2

•(1 + п ),

К = Ь + п2); Ь2п = ъ2 •( + п).

(19)

(20)

Во избежание путаницы в обозначениях, далее обозначим искомую жесткость как х. Тогда с учетом (17)-(20) уравнение (13) примет вид

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

л-(

-(х + а2пх + а3 )-Б =

с

Ь1пХ + Ь2 пХ + Ь3

(21)

Введем новые обозначения коэффициентов:

с1 =-Ба1 -(1 + п2 ); с2 =-Ба2-(1 + п2); с3 =-Ба3 + А; * Ь1 ^ + п2) • * Ь2 ^ + п2) • * Ьз

а = с ; *2 = с ; = с •

В результате получим разрешающее уравнение задачи 4-й степени в виде

(22)

С х + С2 х + С3 —

1

х2 + *2 х + *3

или

с1*1 х4 + (с1*2 + с2) х3 + (схй3 + с2*2 + с3йх) х2 + (с2й3 + съй2) х + (с3*3 -1) = 0.

(23)

(24)

Решение алгебраического уравнения 4-й степени (24) позволит определить искомое значение жесткости опор, которым будет один из его 4-х корней.

2.2. Аналитическое решение разрешающего уравнения

Как известно, 4-я степень уравнения является наивысшей, для которой существуют аналитические методы решения. В данной работе для решения разрешающего уравнения (24) использовался метод, разработанный Ю. А. Несмеевым в работах [46; 47]. Согласно этому методу, преобразуем уравнение (24) так, чтобы коэффициент при старшей степени стал равным единице:

х4 + а3 х3 + а2 х2 + а1 х + а0 = 0, (25)

где

с*-} + с-2 сл<*3 + с-2 + с 3*1 с-2 + с3*2 7

а3 = -а2 = -—-а1 = -а0 = с3а3 -1. (26)

с1а1 схах сгаг

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

а = 1; Ь =-а2; с = а1а3-4а0; а = -а2-а0а^ + 4а0а2 (27)

Решаем это кубическое уравнение в канонической форме, для этого находим его коэффициенты:

9a 2

Рассчитываем коэффициент r как

Далее находим характерный угол ф/

3ac - b2 (ab) bc d 3 2

p=——; q——+~; k = p

27 6a2 2a

= sign (q)• v|p|.

ф = arctg

3 -1

V

Из трех корней кубического уравнения выбираем первый:

л-ф^ b

щ = 2r•cos

(28)

(29)

(30)

(31)

Далее решение сводится к двум квадратным уравнениям, вспомогательные коэффициенты которых равны

2 2 a3 Щ

ix = + Щ -a2; d2 = a0; d3 = a3u1 -2a1.

Первое квадратное уравнение имеет вид

k"i X + ^2 x + кз — 0 ,

где

Его корни равны

к1 = 1; к2 = Ol. 1 2 2

• к = U1 -1; кз = 2

X1 = 2Г(-к2 Wк2 - 4к1к3

Аналогично решаем второе квадратное уравнение:

l1 x +l2 x +l3 = 0 ,

где

/ = Г / = ^3- d. ! = Ü1 h ~ ^ l2 ~ 2 V 1' 'з _ 2

Корни равны:

X4 = (-12 ^l22 - 4/1/з ) .

Искомая жесткость опор будет наименьшей из положительных корней (35) и (38):

С1 = x = min (xt: xt > 0, i = 1,2,3,4) . Абсолютные жесткости опор согласно (10) определятся, как

FI

к = C1--min-

FI

к2 = min

(32)

(33)

(34)

(35)

(36)

(37)

(38)

(39)

(40)

/ 2 1 /

Жесткости опор (40) будут обеспечивать требуемые значения коэффициентов опор (11) и, следовательно, значение первой собственной частоты балки и/или первой критической силы. Для автоматизации расчетов по приведенной выше методике была написана программа для ЭВМ.

3. Пример расчета

В известной научной литературе задача расчета колебаний или устойчивости балки с учетом жесткости опор имеет весьма ограниченное распространение и зачастую представлена в общем виде или в виде численных расчетов, что затрудняет сравнение. Поэтому методика была верифицирована путем сравнения результатов расчета с численным решением методом конечных элементов в программе Ansys. Использовалась балочная модель из конечных элементов типа ВЕАМ189, жесткость опор задавалась конечными элементами СОМВ1Ш4.

Расчетная схема балки показана на рис. 1 и имеет следующие характеристики: длина I = 0,5 м; круглое поперечное сечение диаметром Б = 25 мм; материал - сталь Е = 2*105 МПа, плотность р = 7800 кг/м3, КТР aí = 1,2*10-5 1/°С. Шарнирные опоры балки имеют равную жесткость: к = к = к2 или С = С = С2 (10). Необходимо определить жесткость опор, которая обеспечит заданное значение первой собственной частоты колебаний Л при вариации температуры АТ.

При заданных геометрических параметрах и материале балки, а также при отсутствии продольной силы такая конструкция обеспечивает первую собственную частоту Л = 200,7 Гц для жесткости опор С = 0 и Л = 455 Гц для С = да. Зададим значение Л = 300 Гц и будем поднимать температуру с шагом 50 °С, компенсируя падение значения первой собственной частоты с помощью увеличения жесткости опор.

Результаты расчетов сведены в табл. 2, а на рис. 2 показан один из вариантов расчета в программе Ansys. В табл. 2, помимо заданной температуры, приведено эквивалентное значение продольной силы в балке Рэкв. Сравнение результатов расчета выполнено по значению первой собственной частоты колебаний балки при заданной температуре и жесткости опор.

Таблица 2

Результаты расчетов

Исходные данные расчета

f1, Гц 300

ДТ, °С 0 50 100 150 200 250 282

-^зев кН 0 58,9 117,8 176,7 235,6 294,5 332,2

Результаты по предлагаемой методике

C 4,26 5,97 8,47 21,37 34,72 51,11 1000

к, Н*м 33296 46638 66195 166975 271287 399301 7812500

Результаты МКЭ

fi, Гц 295,9 289,4 284,4 311,6 304,2 286,4 296,5

д/1, % 1,37 3,53 5,2 -3,87 -1,4 4,53 1,17

fiOCAL SOLUTION 5TSP-1 3VB -1 'J5W* <АТО) RSY3-0 ШХ -1.16081 5МХ »1, If ОМ ANSYS ММ ASP. 20 2022 16(12:51

.120963 . ЗвбМВ +71Э89В .М4913 .9026?« i.UDM

Рис. 2. Пример расчета при ДТ = 282 °С и С = 1000 Fig. 2. Example of calculation at ДТ = 282 °С and C = 1000

Полученные результаты показывают хорошую сходимость по всему рассмотренному диапазону жесткостей опор.

4. Обсуждение

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

В работе получено аналитическое решение задачи определения жесткостей опор балки для обеспечения значений ее первой собственной частоты колебаний и первой критической силы. Постановка задачи свелась к алгебраическому уравнению четвертого порядка (24), решение которого можно графически представить как пересечение двух кривых, составленных из разных частей исходного уравнения (13), как показано на рис. 3. Первая положительная точка пересечения функций определяет требуемую жесткость опор С1 (39).

Рис. 3. Графическая интерпретация решения Fig. 3. Graphical interpretation of the solution

Точность расчета жесткости сильно зависит от точности аппроксимирующих зависимостей, что весьма сложно обеспечить применением квадратичных функций (15), (16), поскольку исходные кривые (11) сильно нелинейные. Необходимое повышение точности достигнуто разбиением диапазона жесткостей на три зоны, аналогично тому, как это было сделано в работе [43].

Особенностью решения по предложенной методике является необходимость задания ненулевого значения осевой силы (7), поскольку в противном случае в разрешающем уравнении (23) левая часть окажется равной нулю и вместо уравнения 4-й степени мы получим квадратное уравнение, решение которого гораздо проще, но невозможно по предлагаемой методике. То же касается и частоты /¡, нулевое значение приведет к делению на ноль в правой части уравнения (23). При необходимости такой расчет может быть проведен, если задать малое значение соответствующей величины.

Оценить погрешность предложенного аналитического решения в целом весьма затруднительно, поскольку каждая из аппроксимирующих функций (15), (16) имеет свои отклонения по диапазону, которые накладываются друг на друга при решении разрешающего уравнения (23). Проведенные сравнительные расчеты с вариативностью значений продольной силы и требуемого значения первой собственной частоты колебаний показали точность расчета 5 %, что вполне приемлемо для инженерных расчетов в первом приближении. Погрешность существенно увеличивается при п > 10, т. е. жесткости опор сильно различаются между собой.

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

Заключение

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

Благодарности

Исследование выполнено при финансовой поддержке Российского фонда фундаментальных исследований, Правительства Красноярского края и Красноярского краевого фонда науки в рамках научного проекта № 20-48-242922.

Acknowledgements

The research was funded by RFBR, Krasnoyarsk Territory and Krasnoyarsk Regional Fund of Science, project number 20-48-242922.

Библиографические ссылки

1. Некоторые аспекты моделирования динамики трансформируемых космических конструкций / Ц. Джан, В. Н. Зимин, А. В. Крылов, С. А. Чурилин // Сибирский журнал науки и технологий. 2019. Т. 20, № 1. С. 68-73. DOI: 10.31772/2587-6066-2019-20-1-68-73.

2. Кудрявцев И. В. Обеспечение динамического состояния прямолинейных волноводных трактов при нагреве с помощью расстановки опор // Вестник Московского авиационного института. 2021. Т. 28, № 4. С. 92-105. DOI: 10.34759/vst-2021-4-92-105.

3. Тимошенко С. П., Янг Д. Х., Уивер У. Колебания в инженерном деле. М. : Машиностроение, 1985. 472 с.

4. Бабаков И. М. Теория колебаний. М. : Дрофа, 2004. 591 с.

5. Журавлев В. Ф., Климов Д. М. Прикладные методы в теории колебаний. М. : Наука, 1988. 328 с.

6. Ильин М. М., Колесников К. С., Саратов Ю. С. Теория колебаний. М. : МГТУ, 2001. 272 с.

7. Яблонский А. А., Норейко С. С. Курс теории колебаний. СПб. : Лань, 2003. 254 с.

8. Пановко Я. Г. Введение в теорию механических колебаний. М. : Наука, 1991. 256 с.

9. Блехман И. И. Вибрационная механика. М. : Физматлит, 1994. 400 с.

10. Клаф В. К. Динамика сооружений. М. : Стройиздат, 1979. 320 с.

11. Доев В. С. Поперечные колебания балок. М. : КНОРУС, 2016. 412 с.

12. Balachandran B. Vibrations. Toronto: Cengage Learning, 2009. 737 p.

13. Benaroya H., Nagurka M., Han S. Mechanical vibration. CRC Press: London, 2017. 602 p.

14. Leissa A. W. Vibration of continuous systems, McGraw-Hill: New York, 2011. 524 p.

15. Bottega W. J. Engineering vibrations. CRC Press: New York, 2006. 750 p.

16. Meirovitch L. Fundamentals of vibrations. McGraw-Hill, Book Co: New York, 2001. 826 p.

17. Clough R. E. Dynamics of Structures. McGraw-Hill College: New York, 1995. 752 p.

18. Shabana A. S. Theory of vibration. Springer-Verlag: New York, 2019. 382 p.

19. Geradin M., Rixen D. J. Mechanical vibrations. John Wiley & Sons: London, 2015. 617 p.

20. Rao S. Mechanical vibrations. Pearson Education Limited: London, 2018. 1295 p.

21. Hagedorn P. Vibrations and waves in continuous mechanical systems. John Wiley & Sons: New Jersey, 2007. 388 p.

22. Kelly S. G. Mechanical vibrations. Theory and applications. Cengage Learning: NY, 2012. 896 p.

23. Rades M. Mechanical vibrations II. Printech Publisher: Turin, 2010. 354 p.

24. Inman D. J. Engineering vibration, Pearson Education: NJ, 2014. 720 p.

25. Jazar R. N. Advanced vibrations. A modern approach. Springer: New York, 2013. 695 p.

26. Kelly S. G. Advanced vibration analysis. CRC Press: New York, 2007. 650 p.

27. Тимошенко С. П. Устойчивость стержней, пластин и оболочек. М. : Наука, 1971. 807 с.

28. Алфутов Н. А., Колесников К. С. Устойчивость движения и равновесия. М. : МГТУ, 2003. 256 с.

29. Farshad M., Stability of Structures. Elsevier Science B. V.: Amsterdam, 1994. 434 p.

30. Jerath. S., Structural Stability Theory and Practice: Buckling of Columns, Beams, Plates, and Shells. John Wiley & Sons: Chichester, 2020. 672 p.

31. Timoshenko S. P., Gere J. M., Theory of Elastic Stability. Dover Publications: New York,

2009. 560 p.

32. Thomsen J. J. Vibrations and stability. New York, 2003. 420 p.

33. Yoo C. H. Stability of structures. Elsevier: London, 2011. 529 p.

34. Ziemian R. D., Guide to Stability Design Criteria for Metal Structures. JohnWiley&Sons: NY,

2010. 1117 p.

35. Бидерман В. Л. Теория механических колебаний. М. : Высшая школа, 1980. 408 с.

36. Биргер И. А., Пановко Я. Г. Прочность, устойчивость, колебания. Т. 3. М. : Машиностроение, 1988. 567 с.

37. Коренев Б. Г. Справочник по динамике сооружений. М. : Стройиздат, 1972. 511 с.

38. Уманский А. А. Справочник проектировщика. Т. 2. М. : Стройиздат, 1973. 415 с.

39. Blevins R. D. Formulas for dynamics, acoustics and vibration. John Wiley & Sons, Ltd: Chichester, 2016. 458 p.

40. Wang C. M. Exact solutions for buckling of structural members. CRC Press: New York, 2005, 212 p.

41. Galef A. E. Bending frequencies of compressed beams // Journal of the Acoustical Society of America. 1968. Vol. 44(2). P. 643. DOI: 10.1121/1.1911144.

42. Bokaian A. Natural frequencies of beams under compressive axial loads // Journal of Sound and Vibration. 1988. Vol. 126(1). P. 49-65. DOI: 10.1016/0022-460X(88)90397-5.

43. Кудрявцев И. В., Рабецкая О. И., Митяев А. Е. Аппроксимация значений коэффициентов опор балки при колебаниях и потери устойчивости // Сибирский аэрокосмический журнал. 2022. Т. 23, № 3. С. 461-474.

44. Дрейпер Н. Смит Г. Прикладной регрессионный анализ. М. : Вильямс, 2016. 912 с.

45. Ивченко Г. И., Медведев Ю. И. Математическая статистика. М. : URSS, 2014. 352 с.

46. Несмеев Ю. А. Об одном подходе к решению алгебраических уравнений 3-й и 4-й степеней // Вестник Томского гос. ун-та. Математика и механика. 2011. № 1(13). С. 26-30.

47. Несмеев Ю. А. Развитие одного подхода к решению алгебраического уравнения 4-й степени // Вестник Томского гос. ун-та. Математика и механика. 2013. № 4(24). С. 29-38.

References

1. Zhang Z., Zimin V. N., Krylov A. V., Churilin S. A. [The definite questions of simulation of transformable space structures dynamics]. Siberian Journal of Science and Technology. 2019, Vol. 20, No. 1, P. 68-73 (In Russ.). DOI: 10.31772/2587-6066-2019-20-1-68-73.

2. Kudryavtsev I. V. Ensuring dynamic state of straight waveguide paths at heating by supports arrangement. Aerospace MAI Journal. 2021, Vol. 28, №. 4, P. 92-105. DOI: 10.34759/vst-2021-4-92-105. (In Russ.).

3. Timoshenko S. P., Yang D. Kh., Uiver U. Kolebaniya v inzhenernom dele [Vibrations in Engineering]. Moscow, Mashinostroenie Publ., 1985, 472 p.

4. Babakov I. M. Teoriya kolebaniy [Theory of vibrations]. Moscow, Drofa Publ., 2004, 591 p.

5. Zhuravlev V. F., Klimov D. M. Prikladnye metody v teorii kolebaniy [Applied methods in vibration theory]. Мoscow, Nauka Publ., 1988, 328 p.

6. Il'in M. M., Kolesnikov K. S., Saratov Yu. S. Teoriya kolebaniy\Theory of vibrations]. Мoscow, MGTU Publ., 2001, 272 p.

7. Yablonskiy A. A., Noreyko S. S. Kurs teorii kolebaniy [Oscillation theory course]. Sankt-Peterburg, Lan' Publ., 2003, 254 p.

8. Panovko Ya. G. Vvedenie v teoriyu mekhanicheskikh kolebaniy. [Introduction to the theory of mechanical Vibrations]. Moscow, Nauka Publ., 1991, 256 p.

9. Blekhman I. I. Vibratsionnaya mekhanika [Vibration mechanics]. Moscow, Fizmatlit Publ., 1994, 400 p.

10. Klaf V. K. Dinamika sooruzheniy [Dynamics of structures]. Moscow, Stroyizdat Publ., 1979, 320 p.

11. Doev V. S. Poperechnye kolebaniya balok [Transverse vibrations of beams]. Moscow, KNORUS Publ., 2016, 412 p.

12. Balachandran B. Vibrations. Toronto: Cengage Learning, 2009. 737 p.

13. Benaroya H., Nagurka M., Han S. Mechanical vibration. CRC Press: London, 2017, 602 p.

14. Leissa A. W. Vibration of continuous systems, McGraw-Hill: New York, 2011, 524 p.

15. Bottega W. J. Engineering vibrations. CRC Press: New York, 2006, 750 p.

16. Meirovitch L. Fundamentals of vibrations. McGraw-Hill,Book Co: New York, 2001, 826 p.

17. Clough R. E. Dynamics of Structures. McGraw-Hill College: New York, 1995, 752 p.

18. Shabana A. S. Theory of vibration. Springer-Verlag: New York, 2019, 382 p.

19. Geradin M., Rixen D.J. Mechanical vibrations. John Wiley & Sons: London, 2015, 617 p.

20. Rao S. Mechanical vibrations. Pearson Education Limited: London, 2018, 1295 p.

21. Hagedorn P. Vibrations and waves in continuous mechanical systems. John Wiley & Sons: New Jersey, 2007, 388 p.

22. Kelly S. G. Mechanical vibrations. Theory and applications. Cengage Learning: NY, 2012, 896 p.

23. Rades M. Mechanical vibrations II. Printech Publisher: Turin, 2010, 354 p.

24. Inman D. J. Engineering vibration, Pearson Education: NJ, 2014, 720 p.

25. Jazar R. N. Advanced vibrations. A modern approach. Springer: New York, 2013, 695 p.

26. Kelly S. G. Advanced vibration analysis. CRC Press: New York, 2007, 650 p.

27. Timoshenko S. P. Ustoychivost' sterzhney, plastin i obolochek [Stability of rods, plates and shells]. Moscow, Nauka Publ., 1971, 807p.

28. Alfutov N. A., Kolesnikov K. S. Ustoychivost' dvizheniya i ravnovesiya [Stability of movement and equilibrium]. Moscow, MGTU Publ., 2003, 256 p.

29. Farshad M. Stability of Structures. Elsevier Science B. V.: Amsterdam, 1994, 434 p.

30. Jerath S., Structural Stability Theory and Practice: Buckling of Columns, Beams, Plates, and Shells. John Wiley & Sons: Chichester, 2020, 672 p.

31. Timoshenko S. P., Gere J. M. Theory of Elastic Stability. Dover Publications: New York,

2009, 560 p.

32. Thomsen J. J. Vibrations and stability. New York, 2003, 420 p.

33. Yoo C. H. Stability of structures. Elsevier: London, 2011, 529p.

34. Ziemian R. D., Guide to Stability Design Criteria for Metal Structures. JohnWiley&Sons: NY,

2010, 1117 p.

35. Biderman V.L. Teoriya mekhanicheskikh kolebaniy [Theory of mechanical vibrations]. Moscow, Vysshaya shkola, 1980. 408 p.

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

36. Birger I. A., Panovko Ya. G. Prochnost', ustoychivost', kolebaniya [Strength, stability, vibrations]. Vol. 3. Moscow, Mashinostroenie Publ., 1988, 567 p.

37. Korenev B. G. Spravochnik po dinamike sooruzheniy [Structure dynamics guide]. Moscow, Stroyizdat Publ., 1972, 511 p.

38. Umanskiy A. A. Spravochnik proektirovshchika [Designer's Handbook]. Vol. 2. Moscow, Stroyizdat Publ., 1973, 415 p.

39. Blevins R. D. Formulas for dynamics, acoustics and vibration. John Wiley & Sons, Ltd: Chichester, 2016, 458 p.

40. Wang C. M. Exact solutions for buckling of structural members. CRC Press: New York, 2005, 212 p.

41. Galef A. E. Bending frequencies of compressed beams. Journal of the Acoustical Society of America. 1968, Vol. 44(2), P. 643. DOI: 10.1121/1.1911144.

42. Bokaian A. Natural frequencies of beams under compressive axial loads. Journal of Sound and Vibration. 1988, Vol. 126(1), P. 49-65. DOI: 10.1016/0022-460X(88)90397-5.

43. Kudryavtsev I. V., Rabetskaya O. I., Mityaev A. E. [Approximation of beam support coefficient values at vibrations and buckling]. Siberian Aerospace Journal. 2022, Vol 23, No. 3, P. 461-474 (In Russ.).

44. Dreyper N. Smit G. Prikladnoy regressionnyy analiz [Applied regression analysis]. Moscow: Vil'yams Publ., 2016, 912 p.

45. Ivchenko G. I., Medvedev Yu. I. Matematicheskaya statistika [Mathematical statistics]. Мoscow, URSS, 2014, 352 p.

46. Nesmeev Yu. A. An approach to solution of algebraic equations of the third and fourth degrees. Tomsk State University Journal of Mathematics and Mechanics. 2011, No. 1(13), P. 26-30. (In Russ.).

47. Nesmeev Yu. A. The development of an approach for the solution of the fourth degree algebraic equation. Tomsk State University Journal of Mathematics and Mechanics. 2013, No. 4(24), P. 29-38 (In Russ.).

(© Рабецкая О. И., Кудрявцев И. В., Митяев А. Е., 2022

Рабецкая Ольга Ивановна - кандидат технических наук, доцент кафедры технической механики; Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева. E-mail: olga_rabez@mail.ru.

Кудрявцев Илья Владимирович - кандидат технических наук, доцент кафедры прикладной механики; Сибирский федеральный университет. E-mail: ikudryavtsev@sfu-kras.ru.

Митяев Александр Евгеньевич - кандидат технических наук, заведующий кафедрой прикладной механики; Сибирский федеральный университет. E-mail: pi-prm@mail.ru.

Rabetskaya Olga Ivanovna - Cand. Sc., assistant professor, department of technical mechanics; Reshetnev Siberian State University of Science and Technology. E-mail: olga_rabez@mail.ru.

Kudryavtsev Ilya Vladimirovich - Cand. Sc., assistant professor, department of applied mechanics; Siberian Federal University. E-mail: ikudryavtsev@sfu-kras.ru.

Mityaev Alexander Evgenievich - Cand. Sc., head of the department of applied mechanics; Siberian Federal University. E-mail: pi-prm@mail.ru.

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