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

Обложка

Цитировать

Полный текст

Аннотация

Рассмотрена проблема расчета первой собственной частоты колебаний и первой критической силы для балки с упругими опорами. Аналитический обзор литературы по решению таких задач показал, что в теории колебаний и теории устойчивости стержней учет условий закрепления основан на использовании коэффициентов опор, значения которых были получены после решения соответствующего дифференциального уравнения. В рассмотренной литературе содержится только ограниченный набор значений этих коэффициентов, в основном для идеальных опор простых типов: шарниры, заделка и др. Учет жесткости опор можно найти только в отдельных изданиях и только для ограниченного числа вариантов значений. В данной работе выполнен расчет коэффициентов опор в зависимости от жесткости закрепления балки для первой собственной частоты колебаний и первой критической силы. Полученные значения были разделены на три зоны жесткостей и аппроксимированы внутри каждой зоны квадратичными функциями. Использование квадратичной аппроксимации позволило получить простые аналитические зависимости, пригодные для инженерных прикладных расчетов, а разбиение жесткости на зоны обеспечило приемлемую погрешность получаемых значений. Также квадратичные зависимости позволили решать обратные задачи по определению жесткостей опор для заданного значения первой собственной частоты колебаний или первой критической силы. Проведено подробное исследование погрешности полученных аппроксимирующих функций по всему рассмотренному диапазону жесткостей, которое показало, что погрешность определения коэффициента опор при колебаниях составляет не более 2 %, а при потере устойчивости – 6 %. Погрешность зависит от сочетания жесткостей опор и может увеличиться, если жесткости различаются более чем на порядок. Также была установлена высокая чувствительность решения обратной задачи к входным данным, что является следствием высокой нелинейности зависимости коэффициентов опор от жесткости. Полученные результаты можно использовать при инженерных расчетах первой собственной частоты колебаний и первой критической силы балки с упругими опорами.

Полный текст

Введение*

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

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

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

С целью упрощения таких расчетов существует большое число книг и справочников по динамике балок, в которых предложены готовые решения по свободным колебаниям [3–27] и потери устойчивости [28–35] для простых случаев закрепления, которые полностью ограничивают или разрешают движение по какой-либо координате: шарнир, заделка, скользящая опора и др.,  

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

Однако реализовать на практике упомянутые выше идеальные опоры затруднительно, поскольку действительные опоры имеют некоторое промежуточное значение жесткости, которое влияет на динамическое поведение балки и вызывает расхождение между теоретическим расчетом и практическим результатом. В имеющейся литературе по динамике балок решения для промежуточных значений жесткостей представлены в очень ограниченном виде [36–41] или полученные решения не пригодны для прикладных инженерных расчетов [42–45]. Также справочные таблицы иногда содержат ошибочные значения, которые не всегда очевидны для неопытных расчетчиков. Например, в справочнике по динамике [41] в таблице коэффициентов опор на с. 168 указаны некорректные значения, из-за чего таблица даже потеряла симметричность относительно главной диагонали. Подобные неточности есть и в других изданиях [40; 46].

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

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

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

1. Математическая постановка задачи

Будем рассматривать однородную балку, закрепленную по краям в шарнирных опорах, жесткость которых k1 и k2 может изменяться от нуля (шарнир) до бесконечности (заделка).  

В дальнейших расчетах в качестве жесткости опор также будем использовать не ее абсолютные значения k1 и k2 [Н*м], а относительные безразмерные величины

C1=k1lEJminC2=k2lEJmin,                                                                              (1)

где l – длина; E – модуль Юнга материала; Jmin – минимальный момент инерции поперечного сечения.

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

1.1. Свободные колебания балки в упругих опорах

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

с жесткостью С1 и С2.

Уравнение свободных колебаний имеет вид

EJmin4yx4+m2yt2=0,                                                                                        (2)

где y = y(x) – функция изгибной формы балки; m – удельная масса.

 

Рис. 1. Балка с упругими опорами

Fig. 1. Spring-hinged beam.

 

Граничные условия закрепления балки имеют вид

y0,t=yl,t=0;EJmin2y0,tx2=k1y0,tx;EJmin2yl,tx2=k2yl,tx.                                                                                    (3)

Общее решение представим в виде

yx,t=i=1Yixsinωit,                                                                                       (4)

где  – i-я мода колебаний, общее решение для которой имеет вид

Yix=Aisinpix+Bicospix+Cisinhpix+Dicoshpix.                                                    (5)

Подставляя общее решение в граничные условия, получим частотное уравнение в виде [27]

2αi4tanαi2tanhαi2+αi2lk1+k2EJmintanαi2tanhαi2k1k2l2E2Jmin211cosαi2coshαi2=0,     (6)

где αi=pil – коэффициент опор на i-й собственной моде колебаний.

Ограничиваясь рассмотрением первой собственной частоты колебаний (i=1) получаем искомую зависимость с учетом (1) в виде

α1=α1C1,C2.                                                                                                (7)

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

f1C1,C2=α12C1,C212πl2EJminm.                                                             (8)

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

1.2. Потеря устойчивости балки в упругих опорах

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

 

Рис. 2. Изгибная потеря устойчивости балки с упругими опорами

Fig. 2. Flexural buckling of the beam with spring-hinged supports

 

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

M1=k1α1;M2=k2α2.                                                                                (9)

Тогда уравнение устойчивости примет вид

EJmind2ydx2=PyM1lxlM2xl.                                                            (10)

Граничные условия определяются способом закрепления в виде ограничений на прогибы и углы поворота:

yx=0=0;  y'x=0=α1=M1k1;  yx=l=0;  y'x=l=α2=M2k2.            (11)

Общее решение примем в виде функции

y=Asinθx+BcosθxM1PllxM2Plx,                                               (12)

где

θ=PEJmin.                                                                                             (13)

Тогда уравнение упругой линии будет иметь вид

y=M2PsinθlM1Ptgθlsinθx+M1PcosθxM1PllxM2Plx,               (14)

в котором моменты в упругих шарнирах можно найти из системы двух уравнений:

M11+k1Plk1θPtgθl+M2k1θPsinθlk1Pl=0;M1k2Plk2θsinθlPk2θcosθlPtgθl+M2k2θPtgθlk2Pl1=0.           (15)

Ненулевое решение системы (15) найдем, приравняв к нулю ее определитель и решая трансцендентное уравнение [29],

Pik1+1lθitgθilθitgθil1lPik21lθisinθilθisinθil1l=0,         (16)

где i – номер критической силы.

Ограничиваясь рассмотрением первой критической силы (i = 1), получим искомую зависимость с учетом (1) в виде

μ1=μ1C1,C2=πθ1C1,C2l,                                                                           (17)

где µ1 – коэффициент приведения длины или коэффициент опор при устойчивости.

Первая критическая сила тогда определится по уравнению

PcrC1,C2=π2EJminμ1C1,C2l2.                                                                           (18)

Решение разрешающего уравнения задачи устойчивости (16) аналитическими способами затруднительно, поэтому воспользуемся численными методами.

2. Аппроксимация коэффициентов опор

Решение разрешающих уравнений (6) и (16) выполним численным итеративным методом последовательных приближений в специально созданных программах. Решение было выполнено для 500 вариантов сочетаний значений жесткостей опор балки в диапазоне от С = 0. 

до С = 1000. Ограничение значений диапазона жесткостей было вызвано тем, что с одной стороны расширение диапазона требует использование более сложных аналитических зависимостей для аппроксимации, а с другой стороны жесткость С = 1000 фактически соответствует жесткой заделке и дальнейшее увеличение значений жесткости практически не влияет на изменение величины коэффициентов опор.

2.1. Аппроксимация коэффициентов опор при свободных колебаниях

В табл. 1 приведен ограниченный набор результатов для коэффициентов опор a при свободных колебаниях балки для выбранного диапазона жесткостей, включая значение С = ∞.

 

Таблица 1. Значения коэффициентов опор a1(С1, С2)

С2

С1

0

0,1

1

5

10

20

30

40

50

100

200

500

1000

0

3,142

3,157

3,273

3,534

3,665

3,769

3,814

3,839

3,855

3,889

3,907

3,919

3,923

3,926

0,1

 

3,173

3,288

3,548

3,678

3,783

3,825

3,852

3,868

3,902

3,921

3,932

3,936

3,940

1

 

 

3,399

3,652

3,781

3,884

3,929

3,954

3,970

4,004

4,023

4,034

4,038

4,041

5

 

 

 

3,897

4,026

4,131

4,177

4,202

4,219

4,254

4,273

4,285

4,289

4,292

10

 

 

 

 

4,156

4,263

4,310

4,337

4,354

4,390

4,410

4,422

4,426

4,430

20

 

 

 

 

 

4,374

4,422

4,449

4,467

4,504

4,525

4,538

4,542

4,546

30

 

 

 

 

 

 

4,471

4,499

4,517

4,555

4,576

4,589

4,593

4,597

40

 

 

 

 

 

 

 

4,527

4,545

4,583

4,604

4,617

4,622

4,626

50

 

Symmetric

 

 

 

 

 

4,563

4,602

4,623

4,636

4,640

4,644

100

 

 

 

 

 

 

 

 

 

4,641

4,663

4,676

4,681

4,685

200

 

 

 

 

 

 

 

 

 

 

4,684

4,698

4,702

4,707

500

 

 

 

 

 

 

 

 

 

 

 

4,711

4,716

4,720

1000

 

 

 

 

 

 

 

 

 

 

 

 

4,721

4,725

 

 

 

 

 

 

 

 

 

 

 

 

 

4,730

 

Полученные численные данные графически показаны на рис. 3 в виде плоскости a1(С1, С2) и графика a1(С) для случая опор с одинаковой жесткостью: С = С1 = С2.

       

Рис. 3. Графическое представление коэффициентов опор при свободных колебаниях: a – график функции α1(С1, С2);  б – график функции α1(С)

Fig. 3. Graphs of support coefficients at free vibration: agraph of the function α1(С1, С2);  bgraph of the function α1(С)

 

Результаты показали, что зависимость α1(С1, С2) является сильно нелинейной, что не позволяет аппроксимировать ее одной простой аналитической функцией даже для принятого диапазона С = 0…1000. Исследования установили, что высокая точность аппроксимации может быть получена только при использовании полиномов высокой степени, однако это создает большие вычислительные трудности и не позволяет решать обратную задачу по определению требуемой жесткости опор для заданного значения коэффициента опор α1.

Учитывая симметрию данных в табл. 1, необходимость сохранения простых аналитических зависимостей и возможность проведения обратного расчета, примем в качестве аппроксимирующих зависимостей квадратичные функции. Неизбежное при этом снижение точности аппроксимации компенсируем путем разбиения выбранного диапазона жесткостей на три зоны: I – С = 0…10, II – С = 10…100, III – С = 100…1000. Внутри каждой зоны аппроксимируем полученные числовые значения коэффициентов опор a отдельной квадратичной функцией, коэффициенты которых определим методом наименьших квадратов [47; 48]. В результате получаем решение для каждой зоны в виде:

α1iC1,C2=a0i+a1iC1+C2+a2iC12+C22,                                                 (19)

где i – номер зоны, i= I, II, III.

Коэффициенты квадратичных функций для каждой зоны приведены в табл. 2.

 

Таблица 2. Коэффициенты функций α1i(С1, С2) 

Коэффициент

Зона i

I

II

III

C = 0…10

C = 10…100

C = 100…1000

a0i

3,1415

4,03

4,62

a1i

0,107

0,0075

1,25E-04

a2i

–0,00563

–4,35E-05

–7,400E-08

 

Полученные коэффициенты квадратичных аппроксимирующих функций (19) позволяют рассчитывать значения коэффициента опор a для уравнения (8) практически с любым сочетанием жесткостей опор. Однако погрешность определения коэффициента существенно увеличивается при большой разнице в жесткостях опор. Если же значения жесткостей опор одного порядка, то погрешность расчета коэффициента опор составляет не более 2 % для всего рассматриваемого диапазона жесткостей.

2.2. Аппроксимация коэффициентов опор при потере устойчивости

Используя аналогичный подход, выполним аппроксимацию коэффициентов опор µ при  

потере устойчивости балки. В табл. 3 приведен ограниченный набор результатов для коэффициентов опор µ1(С1, С2) для выбранного диапазона жесткостей, включая значение С = ∞.

Полученные численные данные графически показаны на рис. 4 в виде плоскости µ1(С1, С2) и графика µ1(С) для случая опор с одинаковой жесткостью: С1 = С2.

Зависимость µ(С1, С2) также является сильно нелинейной, что требует разбиения рассматриваемого диапазона жесткости на зоны. Внутри каждой зоны аппроксимируем значения квадратичными функциями

μ1iC1,C2=b0i+b1iC1+C2+b2iC12+C22,                                                    (20)

где i – номер зоны, i = I, II, III.

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

 

Таблица 3. Значения коэффициентов опор µ1(С1, С2) 

С2

С1

0

0,1

1

5

10

20

30

40

50

100

200

500

1000

0

1,000

0,987

0,920

0,803

0,760

0,732

0,722

0,716

0,713

0,706

0,703

0,701

0,700

0,699

0,1

 

0,979

0,913

0,797

0,754

0,727

0,717

0,711

0,708

0,702

0,698

0,696

0,695

0,695

1

 

 

0,854

0,750

0,711

0,686

0,676

0,672

0,669

0,663

0,660

0,658

0,657

0,657

5

 

 

 

0,661

0,627

0,605

0,596

0,592

0,590

0,585

0,582

0,581

0,580

0,580

10

 

 

 

 

0,594

0,572

0,564

0,560

0,558

0,553

0,550

0,549

0,548

0,548

20

 

 

 

 

 

0,551

0,543

0,539

0,537

0,532

0,529

0,528

0,527

0,527

30

 

 

 

 

 

 

0,535

0,531

0,529

0,524

0,521

0,520

0,520

0,519

40

 

 

 

 

 

 

 

0,527

0,525

0,520

0,518

0,516

0,516

0,515

50

 

Symmetric

 

 

 

 

 

0,522

0,517

0,515

0,514

0,513

0,513

100

 

 

 

 

 

 

 

 

 

0,513

0,510

0,509

0,508

0,508

200

 

 

 

 

 

 

 

 

 

 

0,508

0,506

0,506

0,505

500

 

 

 

 

 

 

 

 

 

 

 

0,505

0,504

0,504

1000

 

 

 

 

 

 

 

 

 

 

 

 

0,503

0,503

 

 

 

 

 

 

 

 

 

 

 

 

 

0,500

 

Рис. 4. Графическое представление коэффициентов опор при потере устойчивости: аграфик функции µ1(С1, С2);  бграфик функции µ1(С)

Fig. 4. Graphs of support coefficients at buckling: agraph of the function µ1(С1,С2);  bgraph of the function µ1(С)

 

Таблица 4. Коэффициенты функций µ1i (С1, С2) 

Коэффициент

Зона i

I

II

III

C = 0…10

C = 10…100

C = 100…1000

b0i

1

0,623

0,516

b1i

–0,0504

–0,0016

–1,81E-05

b2i

0,003

1,05E-05

1,2E-08

 

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

2.3. Решение обратной задачи

Использование аппроксимирующих квадратичных функций позволяет решать и обратную задачу по определению жесткости опор, обеспечивающих требуемые коэффициенты опор a или µ. Рассмотрим определение жесткости опор для известного значения коэффициента опор a1 при свободных колебаниях, которое можно найти согласно (8) при заданном значении первой собственной частоты колебаний f1 как

α1C1,C2=lf124π2mEJmin4.                                                                 (21)

Равенство (21) возможно при различном сочетании значений жесткостей опор С1 и С2, что делает задачу нахождения жесткости опор неоднозначной. Для определенности свяжем эти жесткости через заданный коэффициент n, определяющий отношение жесткостей опор:

C=C1=C2n.                                                                               (22)

Подставив эти данные в уравнение (19), получим разрешающее уравнение относительно общей жесткости опор С:

ai21+n2С2+ai11+nC+ai0α1=0.                                         (23)

Значения коэффициентов аij принимаются из табл. 2 для требуемой зоны i и решение с учетом знаков имеет вид

С=a1i1+n+a1i21+n24a2i1+n2a0iα12a2i1+n2.                       (24)

В случае решения обратной задачи для задачи устойчивости, требуемое значение коэффициента опор µ при заданном значении критической силы определится по (18) как

μ1C1,C2=π2EJminl2Pcr.                                                                      (25)

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

С=b1i1+nb1i21+n24b2i1+n2b0iμ12b2i1+n2,                         (26)

где коэффициенты bij определяются для i-й зоны, согласно табл. 4.

В обоих случаях решение будет корректным при определении жесткости опор внутри одной i-й зоны. После определения неизвестной С, требуемые значения жесткости опор определятся как:

C1=C;C2=Cn                                                                       (27)

и

k1=CEJminlk2=CEJminln .                                                          (28)

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

Обсуждение

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

Значения рассчитанных коэффициентов опор в табл. 1 и 3 сравнивались со значениями в известной литературе [36–41]. Сравнение показало, что найденные значения полностью совпадают по первым двум значащим цифрам и иногда различаются в третьей, что, по-видимому, вызвано способом округления. Прямое сравнения всех найденных в литературе значений показало, что отклонения не превышают 0,36 % для коэффициента опор a и 0,27 % для коэффициента µ.

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

коэффициентов опор по всему диапазону жесткостей при равной жесткости опор С = С1 = С2.

 

Рис. 5. Погрешности аппроксимации квадратичными функциями

Fig. 5. Approximation errors by quadratic functions

 

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

распределяется по синусоиде и отклонения коэффициента a не превышают 2 %, а коэффициента µ  – менее 6,2 % во всем рассмотренном диапазоне жесткостей. Этого вполне достаточно для инженерных расчетов балочных конструкций в первом приближении, поскольку используемые в теории балок допущения могут приводить к таким же отклонениям в решении.

Исследование погрешности аппроксимации при различных сочетаниях жесткостей С1 и С2 показало, что если жесткости опор одного порядка, то отклонения значений коэффициентов опор a и µ не превышают значения, представленные на рис. 5.

При разнице в жесткостях опор более чем в 10 раз погрешность увеличивается до 15 %, что может потребовать несколько итераций при определении требуемого значения коэффициента опор. Также следует отметить, что из-за высокой нелинейности функций a1(С1, С2) и µ1(С1, С2) результат решения обратной задачи (24), (26) напрямую зависит от точности входных данных.

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

Заключение

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

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

 

______________________________

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

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

×

Об авторах

Илья Владимирович Кудрявцев

Сибирский федеральный университет

Автор, ответственный за переписку.
Email: ikudryavtsev@sfu-kras.ru

кандидат технических наук, доцент кафедры прикладной механики

Россия, 660041, Красноярск, Свободный, 79

Ольга Ивановна Рабецкая

Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева

Email: olga_rabez@mail.ru

кандидат технических наук, доцент кафедры технической механики

Россия, 660037, Красноярск, проспект имени газеты «Красноярский Рабочий», 31

Александр Евгеньевич Митяев

Сибирский федеральный университет

Email: pi-prm@mail.ru

кандидат технических наук, заведующий кафедрой прикладной механики

Россия, 660041, Красноярск, Свободный, 79

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

  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. Бабаков И. М. Теория колебаний. М. : Дрофа, 2004. 591 с.
  4. Ильин М. М., Колесников К. С., Саратов Ю. С. Теория колебаний. М. : МГТУ имени
  5. Н. Э. Баумана, 2001. 272 с.
  6. Журавлев В. Ф., Климов Д. М. Прикладные методы в теории колебаний. М. : Наука, 1988. 328 с.
  7. Тимошенко С. П., Янг Д. Х., Уивер У. Колебания в инженерном деле. М. : Машиностроение, 1985. 472 с.
  8. Пановко Я. Введение в теорию механических колебаний. М. : Наука, 1991. 256 с.
  9. Яблонский А. А., Норейко С. С. Курс теории колебаний. СПб. : Лань, 2003. 254 с.
  10. Доев В. С. Поперечные колебания балок. М. : КНОРУС, 2016. 412 с.
  11. Блехман И. И. Вибрационная механика. М. : Физматлит, 1994. 400 с.
  12. Клаф В. К. Динамика сооружений. М. : Стройиздат, 1979. 320 с.
  13. Balachandran B. Vibrations. Toronto : Cengage Learning, 2009. 737 p.
  14. Benaroya H., Nagurka M., Han S. Mechanical vibration. London : CRC Press, 2017. 602 p.
  15. Bottega W. J. Engineering vibrations. New York : CRC Press, 2006. 750 p.
  16. Clough R. E. Dynamics of Structures. New York : McGraw-Hill College, 1995. 752 p.
  17. Geradin M., Rixen D. J. Mechanical vibrations. London : John Wiley & Sons, 2015. 617 p.
  18. Hagedorn P. Vibrations and waves in continuous mechanical systems. New Jersey: John Wiley & Sons, 2007. 388 p.
  19. Inman D. J. Engineering vibration, Pearson Education : NJ, 2014. 720 p.
  20. Kelly S. G. Advanced vibration analysis. New York : CRC Press, 2007. 650 p.
  21. Jazar R. N. Advanced vibrations. A modern approach. Springer : NY, 2013. 695 p.
  22. Kelly S. G. Mechanical vibrations. Theory and applications. Cengage Learning: NY, 2012. 896 p.
  23. Leissa A. W. Vibration of continuous systems, McGraw-Hill : New York, 2011. 524 p.
  24. Meirovitch L. Fundamentals of vibrations. McGraw-Hill,Book Co : New York, 2001. 826 p.
  25. Rades M. Mechanical vibrations II. Structural dynamic modeling. Printech Publisher: Turin, 2010. 354 p.
  26. Rao S. Mechanical vibrations. Pearson Education Limited: London, 2018. 1295 p.
  27. Shabana A. S. Theory of vibration. Springer-Verlag : New York, 2019. 382 p.
  28. Hibbeler R. C. Free vibration of a beam supported by unsymmetrical spring-hinges. J. Appl. Mech, 1975. Vol. 42(2), P. 501–502. doi: 10.1115/1.3423612.
  29. Алфутов Н. А., Колесников К. С. Устойчивость движения и равновесия. М. : МГТУ
  30. имени Н. Э. Баумана, 2003. 256 с.
  31. Тимошенко С. П. Устойчивость стержней, пластин и оболочек. М. : Наука, 1971. 807 с.
  32. Farshad M., Stability of Structures. Elsevier Science B.V. Amsterdam, 1994. 434 p.
  33. Ziemian R. D., Guide to Stability Design Criteria for Metal Structures. John Wiley & Sons: NY, 2010. 1117 p.
  34. Timoshenko S. P., Gere J. M., Theory of Elastic Stability. Dover Publications: New York, 2009. 560 p.
  35. Jerath. S., Structural Stability Theory and Practice: Buckling of Columns, Beams, Plates, and Shells. John Wiley & Sons: Chichester, 2020. 672 p.
  36. Thomsen J. J. Vibrations and stability. New York, 2003. 420 p. 404 p.
  37. Yoo C. H. Stability of structures. Elsevier : London, 2011. 529p.
  38. Биргер И. А., Пановко Я. Прочность, устойчивость, колебания. Т. 3. М. : Машиностроение, 1988. 567 с.
  39. Бидерман В. Л. Теория механических колебаний. М. : Высшая школа, 1980. 408 с.
  40. Коренев Б. Справочник по динамике сооружений. М. : Стройиздат, 1972. 511 с.
  41. Уманский А. А. Справочник проектировщика. Т. 2. М. : Стройиздат, 1973. 415 с.
  42. Wang C. M. Exact solutions for buckling of structural members. CRC Press: New York, 2005. 212 p.
  43. Blevins R. D. Formulas for dynamics, acoustics and vibration. John Wiley & Sons, Ltd: Chichester, 2016. 458 p.
  44. Lin Y. K. Free vibration of a continuous beam on elastic supports. International Journal of Mechanical Sciences. 1962. Vol. 4. P. 409–423.
  45. Lin H., Chang S. C. Free vibration analysis of multi-span beams with intermediate flexible constraints. Journal of Sound and Vibration. 2005. Vol. 281(1–2). P. 155–169. doi: 10.1016/j.jsv.2004.01.010.
  46. Luo J., Zhu S., Zhai W. Exact closed-form solution for free vibration of Euler-Bernoulli and Timoshenko beams with intermediate elastic supports. International Journal of Mechanical Sciences. 2022. Vol. 213. Р. 106842. doi: 10.1016/j.ijmecsci.2021.106842.
  47. Free vibrations of Bernoulli-Euler beams with intermediate elastic support / M. J. Maurizi,
  48. D. V. Bambill, P. M. Bellés et al. // Journal of Sound and Vibration. 2005. Vol. 281(3–5). P. 1238–1239. doi: 10.1016/j.jsv.2004.06.014.
  49. Hibbeler R. C. Erratum: “Free Vibration of a Beam Supported by Unsymmetrical Spring Hinges” (Journal of Applied Mechanics, 1975, 42, pp. 501–502)". ASME // J. Appl. Mech. 1981. Vol. 48(2). P. 449. doi: 10.1115/1.3157647.
  50. Ивченко И., Медведев Ю. И. Математическая статистика. М. : URSS, 2014. 352 с.
  51. Дрейпер Н., Смит Г. Прикладной регрессионный анализ. М. : Вильямс, 2016. 912 с.

Дополнительные файлы

Доп. файлы
Действие
1. JATS XML
2. Рис. 1. Балка с упругими опорами

Скачать (19KB)
3. Рис. 2. Изгибная потеря устойчивости балки с упругими опорами

Скачать (26KB)
4. Рис. 3. Графическое представление коэффициентов опор при свободных колебаниях:

Скачать (216KB)
5. Рис. 4. Графическое представление коэффициентов опор при потере устойчивости

Скачать (203KB)
6. Рис. 5. Погрешности аппроксимации квадратичными функциями

Скачать (111KB)

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

Creative Commons License
Эта статья доступна по лицензии Creative Commons Attribution 4.0 International License.