Abstract and keywords
Abstract (English):
The paper presents a mathematical model of the size-dependent behavior of a linear Bernoulli beam based on the gradient deformation theory. The study of nonlinear dynamics of contact interaction is carried out. The synchronization of oscillations on the basis of wavelet analysis is investigated.

Keywords:
Bernoulli-Euler beam, size-dependent behavior, gradient theory of elasticity, synchronization of oscillations
Text
Publication text (PDF): Read Download

Введение

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

Впервые задача о контакте упругих тел, первоначально соприкасающимися в точке была сформулирована и решена Г.Герцем [1] в 1882. Развитие техники поставило проблему контактных задач в ряд актуальных в механике деформируемого твердого тела. Сложность этих задач обусловило большое число подходов и математических методов при их решении. Принимая во внимание различные условия контакта, связанные с более сложными законами поведения, приводят к внедрению новых и нестандартных моделей, выражаемых с помощью эволюционных вариационных неравенств. Современное состояние исследований зарубежных ученых по данному вопросу можно найти в [2]. В работе [3] рассмотрена контактная задача без трения для конструкции, состоящей из двух балок Бернулли – Эйлера с различными длинами и различными переменными толщинами. Один конец короткой балки защемлен; совпадающий с ним конец длинной балки закреплен шарнирно. Второй конец каждой балки свободен. К длинной балке приложена заданная нагрузка. Под действием этой нагрузки балки испытывают слабый совместный изгиб с возможным отставанием, то есть имеется односторонний контакт. В статье [4] представлены результаты использования методов теории оболочек для решения задач определения напряженно-деформированного состояния и долговечности тонкостенных конструкций. Показано, что одной из разновидности задач для оболочек, которые часто встречаются при расчете различных технических объектов, являются задачи с подвижными границами – контактные задачи. В статье рассмотрен один из эффективных способов решения контактных задач для оболочек вращения – метод виртуальных элементов. В отличие от метода конечных элементов увеличение количества виртуальных элементов не увеличивает размер матрицы жесткости конструкции. Исследованию контактного взаимодействия нано балок и нано пластин по модифицированной моментной теории посвящены работы [5-7].

 

1. Постановка проблемы

 

Основываясь на теории напряжений более высокого порядка Mindlin[7] и Lametal. [8] построена математическая модель нано балок Бернулли-Эйлера , соединенных через краевые условия, под действием поперечной нагрузки. Для учета размерных эффектов вводятся три параметра шкалы длины материала . Для учета контакта между балками, используется винклерова связь между обжатием и контактным давлением между двумя балками [9]:

 

,                                                           (1)

где Y = 1, если W1 > w2 + hk то есть имеется контакт между пластиной и балкой, иначе Y = 1.

Математическая модель контакта двух нано балок, на основе кинематической гипотезы Бернулли-Эйлера, описывается системой разрешающих уравнений:

,                              (2)

где m – номер балки (m=1,2),  - зазор между балками. Система (2) с помощью соотношений:

 

        (3)

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

 

.                                                                                            (4)

 

 

Начальные условия имеют вид:

 

.                                                                           (5)

Система (2, 4,5) сводилась к задаче Коши с помощью метода конечных разностей O(h2). Задача Коши решалась методом Рунге-Кутта 4 порядка. Проведено исследование сходимости по методу конечных разностей, на основании которого выбрано оптимальное количество разбиений. Шаг разбиений для метода Рунге-Кутта определялся по принципу Рунге.

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

Применение вейвлет-преобразования хаотического сигнала  рассматриваемой механической системы, характеризуемого плохой определимостью фазы, дает возможность разработать новый метод выявления фазовой синхронизации. Запишем вейвлетную поверхность , которая характеризует поведение механической системы на каждом временном масштабе s в любой момент времени t. Введем так же интегральное распределение энергии вейвлетного спектра по временным масштабам . Использование вейвлетного преобразования дает возможность вычислить фазу  для каждого временного масштаба s, т.е. возможно характеризовать поведение каждого временного масштаба с помощью ассоциированной  с ним фазы. Следует отметить, что фазовая синхронизация означает, что происходит захват хаотических сигналов, в то время как амплитуды этих сигналов остаются несвязанными друг с другом и выглядят хаотическими. Захват фаз влечет за собой совпадение частот сигналов. Частота хаотического сигнала определяется как средняя скорость изменения фазы . Синхронизированными оказываются те временные масштабы, на которые приходятся наибольшие доли энергии  Остальные временные масштабы по-прежнему остаются несинхронизированными. При этом, фазовая синхронизация ведет к появлению и захвату фаз на синхронизированных временных масштабах s: , здесь ,  - непрерывные фазы 1 и 2 балки, соответствующие синронизированным временным масштабам s. В численных примерах графически эти зоны отмечаются черным цветом.

 

2. Численные результаты

 

Для рассмотренных задач использовались следующие значения параметров: .

В таблице 1 представлены результаты исследований для случая когда все безразмерные параметры  равны 0.3. Приведены следующие характеристики колебаний: ф) сигналы 1 и второй балок; б) фазовая хаотическая синхронизация колебаний; в) и е) спектры Фурье для 1 и 2 балок соответственно; г) и ж) 2Dвейвлет спектры на базе вейвлета Морле; д) и з) фазовые портреты.

График сигналов показывает, что балки контактируют на всем временном интервале и начиная со времени  сигналы имеют зеркальный вид. Таким образом энергии колебаний обеих балок равны. На это указывает и график фазовой хаотической синхронизации – начиная с этого времени синхронизация (темный цвет) происходит на частоте возбуждения. До момента времени синхронизация колебаний происходила на более широком диапазоне частот и как показывают вейвлет спектры (г и ж) энергия первой балки существенно больше чем энергия второй. Спектр Фурье обоих балок показывает, что колебания первой и второй балок существенно хаотичны, но энергия имеет «пики» на одних и тех же частотах: . Вейвлет спектр второй балки указывает на то, что наибольшее значение энергия колебаний (соответствующая квадрату вейвлет коэффициентов ) сосредоточено на низких частотах . На более высоких частотах имеются одномоментные всплески активности, соответствующие переходным процессам, хорошо заметным на графике сигналов а). Вид фазовых портретов для обоих балок несимметричен относительно нуля по оси w и имеет сдвиг для первой балки в сторону отрицательных прогибов, а для второй в сторону положительных прогибов, что указывает опять же на симметричность колебаний балок. Отметим интересную особенность колебательного процесса. В момент времени  на вейвлет спектре обеих балок на частотах  происходит всплеск энергии, который приводит к усложнению сигнала, т.е. его вид становится более хаотичным – это видно графике сигналов а). После  сигнал становится более стабильным. Хотя колебания так же остаются хаотическими.

 

Таблица 1. Характеристики нелинейной динамики контактного взаимодействия нано балок Бернулли-Эйлера

Выводы

 

В работе построена математическая модель контактного взаимодействия нанобалок Бернулли-Эйлера на основе градиентной теории упругости. Проведено исследования динамики при одном значении трех параметров шкалы длины материала. Выявлены интересные особенности контактного взаимодействия нано балок. Изучение фазовой хаотической синхронизации выявило связь между всеми характеристиками колебаний (Фурье и вейвлет спектры, фазовый портрет и сигнал) и фазовой синхронизацией. Это позволяет сделать вывод о необходимости дальнейших углубленных исследований в области фазовой хаотической синхронизации контактного взаимодействия нано балок.

Работа выполнена при финансовой поддержке Российского научного фонда, грант № 16-11-10138-П.

References

1. Hertz, H. GesammelteWerke. Bd. 1. Leipzig, 1895.

2. H. T. Banks, S. Hu, Z. R. Kenz, A brief review of elasticity and viscoelasticity for solids // Advances in Applied Mathematics and Mechanics. 2011. 3(1). P. 1-51

3. Osipenko M. A., Kasatkin A. A. Torque load at the contact unilateral contact beams // proceedings of the Saratov University. New series. Series: Mathematics. Mechanics. Informatics. 2019. T. 19. No. 1. Pp. 69-81.

4. Emelyanov I. G. Mechanics of shells and its applications / / Diagnostics, Resource and Mechanics of materials and structures. 2018. 2. P. 6-28.

5. Krysko V. A., Awrejcewicz J., Papkova I. V., Saltykova O. A., Krysko A.V. Chaotic contact dynamics of two microbeams under various kinematic hypotheses // International Journal of Nonlinear Sciences and Numerical Simulation. 2019. 20(3-4). P. 373-386

6. Krysko V. A., Awrejcewicz J., Saltykova O. A., Papkova I. V., Krysko A.V. On the contact interaction of a twolayer beam structure with clearance described by kinematic models of the first, second and third order approximation. Mechanical Systems and Signal Processing. 2019. 115. P.°696-719

7. Mindlin R.D. Second gradient of strain and surfacetension in linear elasticity. Int J Solids Struct. 1965;1:417-38.

8. Yang F, Chong ACM, Lam DCC, Tong P. Couple stress based strain gradient theory for elasticity. Int J Solids Struct. 2002;39:2731-43.

9. Kantor B. ya. Method for solving contact problems of nonlinear shell theory / B. Ya. Kantor, T. L. Bogatyrenko / / Dokl. UKRAINIAN ACADEMY OF SCIENCES. ser. A. 1986. No. 1. Pp. 18-21.

Login or Create
* Forgot password?