ОБЩЕЕ РЕШЕНИЕ ЗАДАЧИ ОБ ИЗГИБЕ МНОГОСЛОЙНОЙ БАЛКИ В РЯДАХ ФУРЬЕ

Полный текст

Исходная система уравнений Исходная система дифференциальных уравнений, полученная в [1] для задачи изгиба многослойной балки (), записывается в виде В уравнения : ?k, Nk - перемещение и соответственно продольная сила в слое k; Dk, Bk - приведенные жесткость на изгиб и сжатие слоя определяемые из выражений: - модуль Юнга слоя ; - ширина балки; - расстояние от верхней и нижней грани слоя k до нейтральной оси соответственно, - переменная, отсчитываемая от верхней грани слоя k; - толщина слоя k; - приведенная жесткость контактного слоя k на сдвиг и сжатие: Все величины, относящиеся к контактному слою, отмечены символом *. Nf,k - вынужденные усилия в слое k, обусловленные наличием предварительного натяжения, действием температуры и т.д. Пример с учетом данных усилий можно найти, например, в работах [2 - 5]; qk - нормальная равномерно-распределенная нагрузка в слое k. Рис. 1 Модель многослойной балки Из системы уравнений путем подстановки , где - число слоев, можно получить частные случаи разрешающих уравнений. Граничные условия Будем рассматривать шарнирно-опёртую балку (), граничные условия для которой записываются в виде: Рис. 2 Расчетная схема Для данных граничных условий решение системы уравнений можно получить в виде разложения в ряды Фурье по синусам, так как в этом случае все граничные условия будут удовлетворены. Разложение исходных функций в ряды Фурье В общем случае разложение произвольной функции в ряд Фурье по синусам с периодом на интервале записывается в виде где Таким образом получим: Поперечную нагрузку и вынужденные усилия также необходимо разложить в ряды Фурье , где Преобразование исходной системы уравнений Подставим выражения и в систему уравнений . После преобразований получим систему из двух разрешающих алгебраических уравнений для слоя k: в которой Система разрешающих уравнений для балки в целом может быть получена из уравнений и представлена в матричной форме для -го члена разложения следующим образом: где - матрица коэффициентов размером ; - вектор неизвестных длиной ; - вектор приложенных усилий длиной . Решением системы уравнений является выражение В развернутой форме систему можно записать в виде: Стоит заметить, что при вычислении коэффициентов в первых и последних двух строках матрицы и вектора необходимо учесть, что . Пример расчёта. Модель трехслойной балки. Разрешающие уравнения для трехслойной балки. В качестве примера расчета рассмотрим трехслойную балку (), состоящую из двух внешних слоев углепластика и внутреннего слоя синтактика на основе полых стеклянных микросфер. Подобные конструкции все чаще находят себе применение в авиастроении и машиностроении, где необходимо обеспечить высокую жесткость и прочность, заменяя собой металлы, имея более низкую массу. Стоит заметить, что слои эпоксидной смолы, скрепляющие субстрат, в данной схеме не учитываются. В данном случае разрешающая система уравнений примет вид: Коэффициенты, входящие в систему уравнений , имеют вид: Рис. 3 Модель трехслойной балки Выражения для n-ых членов разложения функций в ряды Фурье, полученные из системы уравнений весьма громоздки, поэтому здесь не приводятся. При расчете использованы следующие геометрические и физико-механические характеристики Модуль Юнга для синтактика на основе стеклосфер с наполнением 50% получен из результатов испытания на растяжение 5 образцов стержней на испытательной машине MTS insight 100. Из решения системы определяются перемещения и продольные силы для каждого из слоев. Остальные параметры балки выражаются через них следующим образом Исследование сходимости аналитического решения в рядах, по сравнению с численным решением На рисунках 4, 5 приведены кривые, демонстрирующие погрешность вычисления прогибов в середине балки и касательных напряжений на краю для слоя 1, полученные из решения в рядах, по сравнению с численным решением исходной системы для различного числа учтенных членов разложения [4]. Для остальных параметров кривые имеют аналогичный характер. Для численного решения исходной системы уравнений использовался метод Рунге-Кутта 5-го порядка с максимально допустимой погрешностью . Сходимость данных величин в слое 0 будет проанализирована отдельно, в связи с наличием в них краевых эффектов. Погрешность вычислялась следующим образом: где - число учитываемых членов разложения. Рис. 4. Погрешность вычисления прогибов Рис. 5. Погрешность вычисления касательных напряжений в контактном слое Для слоя 0 произведено отдельное сравнение сходимости моментов и поперечных сил, в связи с тем, что сходимость данных величин в результате возникновения краевого эффекта и, соответственно, существенной нелинейности, намного хуже. На рисунках 6, 7 показаны зависимости, демонстрирующие данный факт. Рис. 6. Эпюры моментов для различного числа членов ряда и численного решения Рис. 7. Эпюры поперечных сил для различного числа членов ряда и численного решения На рис. 8, 9 приведены нормированные графики, отражающие характер изменения напряженно-деформированного состояния трехслойной балки в слое 1 при 20 учитываемых членах разложения и слое 0 при 500 членах разложения. Нормирование осуществлялось по следующей формуле: Рис. 8. Перемещения, углы поворота и усилия в слое 1 Рис. 9. Перемещения, углы поворота и усилия в слое 0, а также напряжения в контактном слое 1 На представленных графиках не отображены продольные силы в связи с тем, что характер их распределения совпадает с характером распределения моментов. Заключение Разработана методика решения задачи об изгибе многослойной-шарнирно опертой балки на основе разложения в ряды Фурье по синусам. Данный метод представлен в матричном виде, в результате чего решение может быть легко получено с использованием современным программных комплексов линейной алгебры. Стоит заметить, что ряды для перемещений, углов поворота, моментов и продольных сил достаточно быстро сходятся, в результате решение с допустимой погрешностью в 5% можно получить уже при 10 учитываемых членах разложения. Для поперечных сил и касательных напряжений аналогичную погрешность имеет решение со 100 учитываемыми членами разложения. Однако, для анализа напряженно-деформированного состояния в зоне краевого эффекта необходимо рассматривать решения с учетом порядка 200-500 членов разложения в зависимости от условий нагружения модели. Разработанный метод может использоваться также в задачах ползучести [7, 8] для уточнения краевого эффекта

Об авторах

РОБЕРТ АЛЕКСЕЕВИЧ ТУРУСОВ

Национальный исследовательский Московский государственный строительный университет

Автор, ответственный за переписку.
Email: rob-turusov@yandex.ru

ТУРУСОВ РОБЕРТ АЛЕКСЕЕВИЧ родился в 1940 году в г. Москва, окончил Московский физико-технический институт в 1965 году, доктор физико-математических наук, профессор, Национальный исследовательский Московский государственный строительный университет, 129337, г. Москва, Ярославское шоссе, д. 26. Область научных интересов: физика и механика композитов и полимеров

129337, г. Москва, Ярославское шоссе, д. 26

ВЛАДИМИР ИГОРЕВИЧ АНДРЕЕВ

Национальный исследовательский Московский государственный строительный университет

Email: asv@mgsu.ru

АНДРЕЕВ ВЛАДИМИР ИГОРЕВИЧ родился в 1941 году в г. Москва, окончил Московский физико-технический институт в 1964 году, академик Российской академии архитектуры и строительных наук, доктор технических наук, профессор, заведующий кафедрой сопротивления материалов, главный ученый секретарь Ассоциации строительных вузов, Национальный исследовательский Московский государственный строительный университет, 129337, г. Москва, Ярославское шоссе, д. 26. Область научных интересов: механика неоднородных тел, теория толстых оболочек

129337, г. Москва, Ярославское шоссе, д. 26

НИКИТА ЮРЬЕВИЧ ЦЫБИН

Национальный исследовательский Московский государственный строительный университет

Email: science@nikitatsybin.ru

ЦЫБИН НИКИТА ЮРЬЕВИЧ родился в 1992 году в г. Новомосковск, окончил Национальный исследовательский Московский государственный строительный университет в 2015 году, аспирант кафедры сопротивления материалов, Национальный исследовательский Московский государственный строительный университет, 129337, г. Москва, Ярославское шоссе, д. 26. Тема кандидатской диссертации «Определение напряженно-деформированного состояния многослойных конструкций с применением метода контактного слоя»

129337, г. Москва, Ярославское шоссе, д. 26

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