Аннотация -Данная работа посвящена исследованию нелинейной динамики контактного взаимодей-ствия двух гибких балок Тимошенко, находящихся под действием поперечной знакопеременной нагруз-ки. Учет контактного взаимодействия балок осуществлен по модели Кантора. Геометрическая нелиней-ность учтена по модели Т. фон Кармана. Система дифференциальных уравнений в частных производных двенадцатого порядка методом конечных разностей второго порядка сводится к системе обыкновенных дифференциальных уравнений. Полученная система решается методами типа Рунге -Кутты второго, четвертого и восьмого порядков. Методами нелинейной динамики исследовались хаотические колебания двух гибких балок Тимошенко и найдены оптимальные значения шага по пространственной координате, шага по времени для проведения численного эксперимента. Показано, что получаемые хаотические сиг-налы являются истинными, так как обеспечивается сходимость по всем применяемым численным мето-дам. В работе доказана и обоснована достоверность численных результатов решения задачи контактного взаимодействия двух балок, описываемых кинематической гипотезой Тимошенко, при малом зазоре между ними.Ключевые слова: Контактное взаимодействие, балка Тимошенко, хаос, метод конечных разностей, ме-тоды типа Рунге -Кутты, геометрическая нелинейность I. ВВЕДЕНИЕ Исследование нелинейной динамики и контактного взаимодействия балок является одной из наиболее акту-альных задач современной науки, ввиду широкого спектра применения балочных структур в различных отрас-лях современной промышленности. Известно несколько гипотез, описывающих динамику балочного элемента, это теории первого, второго, третьего приближений и другие. К гипотезе первого приближения относится гипо-теза Эйлера -Бернулли [1]. Гипотеза второго приближения -гипотеза С.П. Тимошенко [2], позволяет учиты-вать поворот нормали к срединной линии после деформации. Данная теория позволяет более точно описывать динамику балочного элемента. Исследованиям нелинейной динамики балок, описываемых гипотезами различ-ных приближений, посвящено большое количество российской и иностранной литературы [3][4][5]. Настоящая работа ставит своей целью доказать истинность хаотических колебаний двух балок Тимошенко с малым зазо-ром, находящихся под действием поперечной знакопеременной нагрузки, с учетом контактного взаимодей-ствия. В отечественной и иностранной литературе решения таких задач нет. Это является принципиально важ-ным вопросом, при решении таких сложных нелинейных систем уравнений численными методами, т.к. может получиться неверное решение из-за погрешности численных методов.В работе будем опираться на определение хаоса, данного Гуликом [6]. Гулик считает, что хаос существует тогда, когда либо имеется существенная зависимость от начальных условий, либо функция имеет положительный показатель Ляпунова в каждой точке области ее определения и поэтому не является в конечном итоге периодиче-ской. В качестве начальных условий, помимо условий, накладываемых на функции, входящие в систему диффе-ренциальных уравнений, будем понимать количество раз...