Full Text

Прямоугольные в плане панели в качестве элементов конструкций находят широчайшее применение в различных областях техники: авиастроении, космической промышленности, судостроении, строительстве, приборостроении и др. Элементы подобных конструкций находятся под действием внешнего динамического давления различного вида и работают в средах с переменными во времени свойствами (случайные ударные воздействия, звуковое давление). Поэтому возникает необходимость комплексного исследования поведения таких систем. Для современного строительства важнейшим аспектом является установление множества параметров воздействия, которые будут характеризовать безопасный и опасный режимы работы элементов ответственных конструкций. Для моделирования случайных во времени свойств внешней среды хорошо подходит белый шум. Как известно, поведение нелинейных динамических систем серьезно реагирует на малейшие изменения начальных условий, что может быть связано в том числе с влиянием шумовых полей. Несмотря на более чем полувековую историю, вопросы моделирования сложной динамики механических систем остаются актуальными [1-6]. Однако влияние шумовых полей на сложные колебания механических систем мало изучено. Показано, что случайные воздействия играют весьма существенную роль в поведении динамических систем в таких областях, как физика, химия, биология [7]. Для динамических систем с одной степенью свободы установлено, то внешние шумы способны вызывать качественную перестройку колебательных режимов системы [8-10]. В работах F. Jinqing, Zhao De-min, Gan Chaunbiao, Manohar C.S., и других серьезное внимание уделяется теоретическим исследования уравнений типа Дуффинга с учетом внешних шумов [11]. В статье В.Д. Потапова показано, что неустойчивая по Ляпунову детерминированная параметрическая система может быть стабилизирована путем наложения на параметрическую нагрузку случайного шума [12]. Статья [13] посвящена влиянию шума на контактное взаимодействие пластины, подкрепленной ребрами жесткости. В работе [14] излагаются результаты опытного изучения балочно-ударной системы при возмущениях типа гауссова шума. Показано, что при учете нескольких форм колебаний реакция нелинейной модели сильно отличается от реакции расчетной модели с одной степенью свободы. В статье [15] также показано, что учет количества степеней свободы существенно влияет на достоверность получаемых результатов. Данная работа посвящена исследованию влияния поля белого шума на нелинейную динамику гибких прямоугольных в плане панелей под действием внешней продольной нагрузки. Вопрос влияния шумовой составляющей на характер колебаний рассматриваемых панелей особенно интересен и важен, так как его касаются впервые и нигде в литературе ранее он освещен не был. Математическая модель колебаний прямоугольной в плане панели основывается на гипотезах Кирхгофа. Геометрическая нелинейность учитывается в форме Т. Кармана [16]. В прямоугольной системе координат трехмерная область имеет вид: (рис. 1). Рис. 1. Схема приложения нагрузки и шумового поля к пластине Математическая модель колебаний панели описывается системой нелинейных дифференциальных уравнений в частных производных, записанных в безразмерном виде [17]: (1) где - известные нелинейные операторы; , и - функция прогиба и усилия соответственно; - коэффициент Пуассона; - удельный вес материала пластины; - модуль Юнга материала пластины; - ускорение свободного падения; , где - размеры пластины в плане по и по соответственно. Система (1) приведена к безразмерному виду с использованием следующих безразмерных параметров: - безразмерные параметры пластины по и по соответственно; ; ; - время, ; - коэффициент демпфирования среды; - внешний белый шум. Черточка над безразмерными параметрами в уравнениях (1) для простоты опущена. К уравнениям (1) присоединяются граничные и начальные условия в зависимости от типа нагружения и закрепления пластины. Методы решения Система (1) с учетом граничных и начальных условий сводится к нелинейной системе обыкновенных дифференциальных уравнений (СОДУ) методом конечных разностей с аппроксимацией по пространственным переменным. Первое уравнение нелинейной СОДУ по времени решается методом Рунге-Кутта 4-го порядка относительно функции прогиба. Выбор метода обоснован в статье [18]. Далее значения функции прогиба подставляются в правую часть второго уравнения системы. Второе уравнение системы (1) становится линейным и решается методом обратной матрицы относительно функции усилия на каждом шаге по времени. Число разбиений в методе конечных разностей . Выбор количества разбиений МКР и обоснование достоверности получаемых результатов приведены в работе [19]. Применение метода конечных разностей позволяет рассматривать пластину как систему с большим числом степеней свободы, в данном численном эксперименте их число 196 (14×14). Для глубокого понимания поведения динамических систем под действием той или иной нагрузки необходима информация, скрытая в частотной области сигнала. Получаемые в результате исследований динамики нелинейных пластин сигналы зачастую имеют сложную структуру по временной оси. Поэтому в работе применялся аппарат вейвлет-преобразований. Основным преимуществом вейвлет-анализа перед Фурье является возможность отследить локализованные особенности нестационарного сигнала. В то время как в коэффициентах Фурье отражается поведение сигнала на всем временном интервале. Если с помощью Фурье-анализа рассматривать сигнал, частотные характеристики которого существенно меняются с течением времени, то на выходе будем иметь «сумму всех его поведений», наложенных друг на друга. Вейвлет спектр сигнала, в свою очередь, серьезно зависит от материнского вейвлета, на основании которого он строится. В работе [19] приведено обоснование выбора в качестве материнского вейвлета Gauss32, который дает лучшую локализацию по частоте по сравнению с ранее используемым вейвлетом Morlet. Следует отметить, что как бы ни была хороша локализация вейвлета по оси частот, на вейвлет-спектре хорошо различимы только первые три-четыре бифуркации системы. Если в спектре колебаний присутствует частота, амплитуда которой значительно превосходит амплитуду другой частоты, то на вейвлет-спектре последняя будет плохо различима. Фурье-спектр дает возможность видеть гораздо большее количество бифуркаций системы, регистрирует все, даже незначительные, движения частот и визуализирует частоты вне зависимости от их амплитуды. Поэтому при исследовании нелинейной динамики распределенных механических систем следует применять Фурье- и вейвлет-анализ совместно, опираясь на сильные стороны каждого из них. Численные эксперименты, проведенные в работе, направлены на изучение влияния интенсивности внешнего шумового поля, нормального к поверхности пластины, на характер колебаний рассматриваемой динамической системы. Продольная нагрузка задавалась в виде где , - частота и амплитуда внешнего воздействия соответственно. Колебания рассматривались на временном интервале . Коэффициент диссипации среды Частота продольной нагрузки что превосходит частоту собственных колебаний системы. Отношение размеров панели в плане . К системе (1) присоединялись неоднородные граничные условия опирания на гибкие несжимаемые (нерастяжимые) ребра [20] при при (2) и нулевые начальные условия (3) Аддитивный шум добавлен в систему в форме случайного слагаемого с постоянной интенсивностью где - интенсивность шума; rand() - стандартная функция языка C++, принимающая случайное целое число от 0 до RAND_MAX, RAND_MAX - константа, равная 65535. Выражение принимает произвольные дробные значения в диапазоне (-1;1). Данная математическая модель белого шума предложена Perry R. Cook и Gary P. Scavone (Центр компьютерных исследовании в области музыки и акустики (CCRMA) Стэндфордского университета). В рамках численного эксперимента были получены две последовательные бифуркации удвоения периода колебаний рассматриваемой динамической системы. После чего система резко переходила к хаосу. В работе [21] показано, что для динамических систем в виде прямоугольных в плане пластин первая бифуркация зачастую является бифуркацией Андронова-Хопфа. Следует отметить, что пластина в рассматриваемом случае (в диапазоне амплитуд нагрузки очень резко реагирует на малейшее изменение нагрузки. Так, двухчастотные колебания сменяются хаотическими при изменении значения амплитуды нагрузки на 0,001 %. Дальнейший рост управляющего параметра на 0,001 приводит колебания системы к двухчастотным Как видно из табл. 1, хаотические колебания системы при фиксированной амплитуде продольной нагрузки меняют свой характер с течением времени. Есть окна, где доминирует частота возбуждающей силы и области, где преобладает частота ( Из последней строки табл. 1 видно, что окна хаоса могут меняться областями с одной или двумя бифуркациями Хопфа. В табл. 1 вейвлет-спектры построены на основании материнского вейвлета Gauss 32, ось ординат - частота, ось абсцисс - время. Таблица 1 Вейвлет- и Фурье-спектры Вейвлет-спектры Фурье-спектры Вейвлет-спектры Фурье-спектры Окончание табл. 1 Вейвлет-спектры Фурье-спектры Анализ влияния внешнего, нормального белого шума на колебания системы проводился в том диапазоне амплитуд внешних продольных нагрузок, где система сильно неустойчива (табл. 2). Интенсивность белого шума в эксперименте бралась следующим образом: Были построены сигналы с шагом по амплитуде продольной нагрузки 0,001. Выявлено, что влияние внешнего шума интенсивности ухудшало состояние системы, при все сигналы были хаотическими. В диапазоне система становилась менее чувствительной к шумовому воздействию. Отдельно стоит отметить в данном случае наблюдалось улучшение состояния системы по сравнению с в диапазоне амплитуды продольной силы Интенсивность внешнего шума улучшала состояние системы в диапазоне Таблица 2 Влияние интенсивности поля белого шума на характер колебаний пластины Окончание табл. 2 Численный эксперимент показал, что амплитуда и характер сигналов совпадают там, где колебания системы в шумовом поле и без него отличны от хаотических. Шум интенсивности снижал амплитуду сигнала в области хаотических колебаний (рис. 2). Проводился анализ влияния интенсивности шумового поля, нормального к плоскости панели, на форму ее колебаний. В случае амплитуды внешней продольной нагрузки на всем временном интервале были получены несимметричные формы колебаний. Наложение шумового поля интенсивности уменьшало степень несимметрии форм панели. В эксперименте с полем белого шума интенсивности несимметричных форм колебаний не обнаружено (табл. 3). Рис. 2. Сигналы при амплитуде продольной нагрузки и Таблица 3 Поверхности панели при Проведенные численные эксперименты показывают, что поле белого шума способно уменьшать амплитуду колебаний панели, сокращать количество частот в спектре колебаний системы и переводить несимметричные формы колебаний к симметрии. Таким образом, можно утверждать, что воздействие на динамическую систему шумовым полем может приводить к более безопасным колебательным режимам. То есть с помощью белого шума можно управлять характером колебаний механической системы.

About the authors

E Y Krylova

NIIM Nizhny Novgorod State University N.I. Lobachevsky

Email: kat.krylova@bk.ru

T V Yakovleva

NIIM Nizhny Novgorod State University N.I. Lobachevsky

Email: yan-tan1987@mail.ru

V G Bazhenov

NIIM Nizhny Novgorod State University N.I. Lobachevsky

Email: bazhenov@mech.unn.ru

References

Statistics

Views

Abstract - 197

PDF (Russian) - 74