Методы исследования нелинейных колебаний систем, связанные с преобразованием нелинейных уравнений движения систем с конечным числом степеней свободы к нелинейным интегральным уравнениям второго рода, достаточно широко распространены в научно-технической литературе [1–3]. В основном речь идет о построении амплитудно-частотных характеристик (АЧХ) и их анализе.

Между тем, воспользовавшись импульсными переходными функциями (ИПФ) исходных линейных «порождающих» уравнений, нелинейные дифференциальные уравнения могут быть записаны в виде систем нелинейных интегральных уравнений второго рода, которые, по существу, и определяют полное решение исходной системы в виде суммы двух интегралов типа свертки. Один из возможных алгоритмов вычисления подобных интегралов дан, в частности, в статье [4].

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

Общая схема и алгоритм расчета нелинейных систем с конечным числом степеней свободы проиллюстрированы на примере расчета системы с двумя степенями свободы (рис. 1).

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

Стоит отметить, что принятая расчетная схема достаточно широко используется при расчетах различных типов виброзащитных систем: виброизоляции, динамических гасителей колебаний и т. п. – с линейными и нелинейными характеристиками. Подобный алгоритм позволил записать решения линейных уравнений практически в замкнутом виде [2] и в качестве второго примера – решение уравнений плоских горизонтально-вращательных колебаний массивных тел (фундаментов), в том числе заглубленных в грунт [1] (рис. 2).

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

На начальном этапе дифференциальные уравнения колебаний систем при нелинейной зависимости «жесткость–перемещение» следует представить в виде (1), записав нелинейные составляющие в правой части:

; (1)

где для дальнейшего сокращения записей обозначено:

(2)

k1, k2 – жесткости упругих связей;

v1, v2 – коэффициент затухания в упругих связях;

Δy = y1 − y2.

Передаточные (ПФ) и импульсные переходные (ИПФ) функции этой системы определяют вид интегральных уравнений, которые, по существу, и являются решением системы нелинейных дифференциальных уравнений. Значение этих функций вычисляют из расчета линейной («порождающей») системы уравнений при:

qi(t) = Qieiωt – силовых;

qi(t) = – 0 Meiωt – кинематических воздействиях.(3)

где Qi – амплитудные значения усилий;

0 – закон смещения основания;

ω – частота вынужденных колебаний.

В матричном виде эту систему следует записать так:

(4)

где

M =  (5)

– матрицы масс и жесткости.

Используя известные зависимости [2][5], решения подобных систем могут быть представлены в виде разложения по формам собственных колебаний при том, что каждая составляющая полного решения определяется из решения систем уравнений – по структуре систем с одной степенью свободы, диссипативные силы в которых можно учесть, добавив в каждое из них дополнительные слагаемые пропорциональные скорости колебаний.

В работе, в частности, используется модифицированная модель Фойгта с основным параметром γ – коэффициентом неупругого сопротивления. Такая модель принята во многих нормативных документах [6][7], где приводятся также значения этого параметра для различных материалов и сред.

Подобная схема построения алгоритмов расчета более детально дана на примере расчета системы (1).

Следующий этап расчета – вычисление ПФ и ИПФ системы линейных уравнений (4).

Следуя общей схеме вычисления передаточных функций линейной системы [2], подставим в (4):

y→ = Y→eiωt = Q→H→eiωt – при силовом и y→к = Y→кeiωt = – 0 MH→eiωt – при кинематическом воздействии, где H→– передаточная функция, и, после сокращения на eiωt, запишем систему уравнений относительно элементов вектора Y→: [(1 + iω2v1)k1 – m1ω2]Y1 – (1 + iω2v1)k1Y2 = Q1; − (1 + iω2v1)k1Y1 + [(1 + iω2v1)k1 + (1 + iω2v2)k2 – m2ω2]Y2 = Q2. Из решения которой следует: (7) где Φij – приближенная (без учета диссипативных сил) амплитуда перемещений i-й массы при действии единичной силы 1 × eiωt на j-ю массу; D – определитель системы (4), который для сокращения выкладок удобно после некоторых преобразований записать в виде:

(8)

где

(9)

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

(10)

где φ1 = h1 + h1s1 + 1.

Следуя общей схеме построения ПФ [3][5], определитель системы (4) и его производную по ω2 запишем в виде:

(11)

(12)

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

(a)

где L(pi) – некоторая функция частоты.

Следует преобразовать и представить как сумму отдельных составляющих:

(б)

где pr – корни уравнения (8) (частоты собственных колебаний), т. е., по существу, в виде разложения по формам собственных колебаний «порождающей» системы, в которой коэффициенты Kr, следуя известной схеме [2], следует определять по формуле:

(в)

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

, (13)

где

. (14)

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

Передаточные и импульсные переходные функции для общих случаев линейных систем определены в [2].

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

,

при r = 1...n, (15)

где Dr – некоторая функция частоты pr.

ИПФ уравнения (15) удобно записать, воспользовавшись одним из решений для свободных колебаний и теоремой сохранения количества движения S = my·, в окончательном виде так

. (16)

Из (13) и (16) следуют такие общие зависимости для ИПФ системы (4)

, (17)

где

(18)

Определив ИПФ из решения системы линейных дифференциальных уравнений (4), решения нелинейной системы (1) можно записать в виде нелинейных интегральных уравнений второго рода. В частности, для несколько упрощенного варианта – при c1(Δy) = K1; c2(y2) ≠ K2, общее решение системы (1) будет определяться интегралами:

y1(t) = y01 (t) + W01(t)

где

(19)

(20)

(21)

(22)

где y0i и W0i – соответственно решение линейной системы (4) и составляющие решений от «фиктивной» нагрузки. Φ2(τ) – по формуле (2).

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

и соответственно ИПФ в виде сумм слагаемых, каждое из которых – произведения функций от t и τ, алгоритм расчета подобных систем можно записать в более компактном виде.

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

(23)

могут вычисляться по формулам [3]:

(24)

где

(25)

Алгоритмы построения решений (перемещений) системы нелинейных дифференциальных уравнений (1), записанных в виде сумм двух интегралов типа свертки (19)–(22), можно проиллюстрировать на примере уравнения (19). Воспользовавшись зависимостями (24), (25), составляющие решения по первой форме собственных колебаний системы линейных «порождающих» уравнений следует записать в виде

(г)

где

(д)

(е)

При вычислении перемещений по второй форме в (г) и (е) p1* следует заменить на p2* и поменять знак перед общим решением (г).

Точность построенных по формулам (24), (25) решений подобных систем оценивалась на примере уравнения (23) [8] при кубической зависимости «жесткость–перемещение» на нагрузки, возбуждаемые при работе виброактивного оборудования с вращающимися частями (насосы, вентиляторы и т. п.) при всех режимах – рабочем,

пуска

и остановки (26)

где постоянные a и b определяют время пуска и остановки.

В принятом алгоритме расчета при вычислении интегралов и, следовательно, полного решения системы (1) в качестве основного параметра, определяющего, главным образом, точность решений, был принят интервал по времени , где Tr – период колебаний, Nr – число разбиений. При погрешность решений по отношению к решениям при Δt, существенно меньше принятого, не превышала одного процента.

При достаточно малых значениях Δt возможно не уточнять величину нелинейной составляющей (не выполнять итерации), а при вычислении интегралов F1, F2 на каждом этапе принимать их средние значения.

Подобный алгоритм был использован и при расчете системы с одной и двумя степенями свободы (рис. 2), дополнительная связь в которых включает демпфер вязкого трения [2]. Зависимости, построенные по результатам расчетов обеих систем, приведены на рис. 3–5.

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