Система автоматизированного проектирования и расчета
+7 (499) 922  00  02

Реализация модели контактного слоя при расчете адгезионного соединения с использованием метода конечных элементов

6 Июня 2019 Реализация модели контактного слоя при расчете адгезионного соединения с использованием метода конечных элементов

Н.Ю. ЦЫБИН1,2
1 НИУ МГСУ, 129337, Москва, Ярославское шоссе, 26
2 ООО «ЦПР Форма», 115114, Москва, 1-й Кожевнический пер., 10
science@nikitatsybin.ru

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


Постановка задачи

Рассмотрим многослойный стержень, состоящих из трех одинаковых слоев субстрата, скрепленных двумя слоями адгезива. Эскизный вид модели приведен на рисунке ниже.

01.png

Рис. 1. Эскизный вид модели

В данной модели центральный слой субстрата закреплен левой гранью от горизонтальных смещений, правая его грань свободна. К двум внешним слоям справа приложена растягивающая нагрузка f01.png .

Возникшая необходимость детального изучения напряженно-деформированного состояния потребовала от исследователей использования различных современных вычислительных методов. Одним из таких является метод конечных элементов.

Решение с использованием метода конечных элементов

Для решения данной задачи с использованием программного комплекса ЛИРА 10.8, реализующим метод конечных элементов, была создана плоская модель. На рисунке ниже приведен общий вид данной модели.

02.png

Рис. 2. Общий вид конечно-элементной модели без контактного слоя

На данном рисунке красные элементы ­– слои адгезива, зеленые – слои субстрата. Решалась упругая задача с признаком схемы №2 (плоская рама). Моделирование слоев адгезива и субстрата осуществлялось с использование пластинчатых конечных элементов №24, 27 (КЭ плоского напряженного состояния). В угловых точках было выполнено сгущение сети.

Ниже записаны физико-механические и геометрические параметры модели

f02.png

При моделировании использовался масштабный коэффициент 100.

Ниже приведены результаты вычисления нормальных и касательных напряжений в модели.

03.png

Рис. 3. Результат вычисления нормальных напряжений f03.png

04.png

Рис. 4. Результат вычисления касательных напряжений f04.png

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

05.png

Рис. 5. Результат вычисления плотности энергии

Красным обведены зоны концентрации напряжений.

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

06-1.png06-2.png

Рис. 6. Зависимость величины максимальных касательных напряжений и плотности энергии в зависимости от размера конечных элементов

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

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

  1. Использование физически-нелинейных моделей деформирования материалов.
  2. Добавление в угловые точки мениска [3], [4].

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

Другие подходы, как правило, являются математическими и их реализация в конечно-элементной модели либо невозможна, либо требует от авторов программного комплекса разработки отдельных маловостребованных расчетных модулей. Одним из дополнительных примеров является добавление по торцам модели «нерастяжимых» мембран.

Для решения данной задачи предлагается использовать модель контактного слоя [5], предложенную [6] и реализованную во многих работах Р.А. Турусовым. Контактный слой это трансверсально ортотропная среда, в которой происходит межмолекулярное взаимодействие слоев адгезива и субстрата.

07.png

Рис. 7. Взаимодействие адгезива и субстрата через контактный слой

Параметры этой среды таковы, что ее можно абстрактно представить, как щетку коротких упругих стерженьков, не контактирующих между собой и ориентированных нормально к поверхности контакта. Так как стерженьки не соприкасаются, модули упругости данной среды в направлениях, перпендикулярных их боковой поверхности и коэффициент Пуассона равны нулю. В результате контактный слой имеет две основных характеристики, толщину f05.png и модуль сдвига f06.png, причем, так как коэффициент Пуассона равен нулю, модуль сдвига связан с модулем Юнга простым соотношением f07.png.

Аналитическое решение рассматриваемой задачи с использованием модели контактного слоя и гипотез Эйлера-Бернулли для слоев адгезива и субстрата можно найти в [7]. Здесь же основной целью является демонстрация способа реализации модели контактного слоя в программном комплексе ЛИРА 10.8.

08.png

Рис. 8. Общий вид конечно-элементной модели с контактными слоями

На представленном рисунке контактные слои синего цвета. В остальном модель полностью соответствует рассмотренной ранее.

09.png

Рис. 9. Результат вычисления нормальных напряжений f03.png

10.png

Рис. 10. Результат вычисления касательных напряжений f04.png

11.png

Рис. 11. Результат вычисления плотности энергии

Как видим, после добавления контактного слоя максимум касательных напряжений сместился от угловой точки. Для наглядности разместим рядом результаты вычисления касательных напряжений в окрестности данной точки для двух рассмотренных моделей.

 12.png  13.png

Рис. 12. Касательные напряжения в угловой точке для модели без контактного слоя

Рис. 13. Касательные напряжения в угловой точке для модели с контактным слоем

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

14-1.png14-2.png

Рис. 14. Зависимость величины максимальных касательных напряжений и плотности энергии в зависимости от размера конечных элементов

Здесь мы видим существенно различный ход кривых представленных на Рис. 6и на данном рисунке. После добавления контактного слоя удалось избежать сингулярности касательных напряжений.

Приложение 1

15.png

Рис. 15. Характеристики материала контактного слоя

16.png

Рис. 16. Направление осей выравнивания напряжений и осей ортотропии

Заключение

В работе показано, что расчет адгезионных соединений в предположении абсолютно жесткого контакта слоев приводит к возникновению бесконечных напряжений. Показано, что для избегания данной проблемы эффективным способом является использование модели контактного слоя. Продемонстрирован простой метод реализации данной модели при расчете с использованием метода конечных элементов в программном комплексе ЛИРА 10.

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

  1. Lertora E., Campanella D., Mandolfino C., Gambaro C., Fratini L., Buffa G. Comparison Between FSW and Bonded Lap Joints – a Preliminary Investigation // AIP Conference Proceedings, Vol. 1896(1):110004, 2017.
  2. Туразян А.В., Рабинович А.Л. // ДАН СССР. 1970. Т. 194. № 6. С. 1305-1307.
  3. Calic A., Yildirim S. Effect of adherend recessing on bi-adhesively bonded single-lap joints // Sadhana, Vol. 42, No. 3, 207. pp. 317-325.
  4. Calic A. Effect of adherend shape on stress concentration reduction of adhesively bonded single lap joint // Engineering Review, Vol. 16, 2016. pp. 29-34.
  5. Турусов Р.А. Адгезионная механика: монография. 2-е изд. Москва: НИУ МГСУ, 2015. 232 с.
  6. Турусов Р.А., Вуба К.Т., Фрейдин А.С. Исследование влияния температурно-влажностных факторов на прочностные и деформационные свойства клеевых соединений древесины со стальной арматурой // Труды ЦНИИСК им. В.А.Кучеренко. 1972. № 24. С. 86-124.
  7. Andreev V.I., Turusov R.A., Tsybin N.Y. The contact layer method in calculating of the shear compounds // MATEC Web of Conferences, No. 117, 2017. pp. 1-7.
  8. Турусов Р.А., Вуба К.Т. Напряженное состояние и особенности оценки прочности адгезионных соединений при отрыве // Физика и химия обработки материалов, Т. 2, 1980. С. 108-115.
  9. Турусов Р.А., Маневич Л.И. Метод контактного слоя в адгезионной механике. Одномерные задачи. Сдвиг соединения внахлестку // Клеи. Герметики. Технологии, № 8, 2009. С. 2-12.
  10. Турусов Р.А., Маневич Л.И. Метод контактного слоя в адгезионной механике // Клеи. Герметики, Технологии. 2009. № 6. С. 2-11.
  11. Турусов Р.А., Маневич Л.И. Введение в адгезионную механику // Клеи. Герметики. Технологии, № 5, 2009. С. 2-8.

Следите за нашими новостями в социальных сетях:

Возврат к списку


Материалы по теме:


Подписка

Вы хотите первыми узнавать о выходе новых версий, проводимых мероприятиях и акциях компании? Подписывайтесь!

Подписаться
Приглашаем на конференцию Строительство, проектирование объектов капитального строительства Республики Крым и города Севастополь
15-16 августа 2019 г.
Республика Крым, г. Симферополь
Научно-практическая конференция "Строительство, проектирование объектов капитального строительства Республики Крым и города Севастополь"
05 августа 2019
Акция: Скидка 50% на ПК ЛИРА 10
Специальное предложение от "ЛИРА софт". Скидка 50% на расчетный комплекс ПК ЛИРА 10 для организаций, работающих в альтернативном программном комплексе
31 июля 2019
Обучись расчетному комплексу ЛИРА 10 бесплатно
До 31 июля вы можете пройти бесплатно базовый курс, состоящий из 5-ти видео-уроков, с пошаговым объяснением выполнения задач и самостоятельной работой для отработки полученных навыков.
09 июля 2019
Семинар в Душанбе «Современные методы моделирования и расчета зданий и сооружений на статические и сейсмические нагрузки и воздействия»
Компания ЛИРА софт приглашает инженеров-расчетчиков, конструкторов, ГИПов, сотрудников экспертизы, сотрудников ВУЗов, а также других специалистов, имеющих отношение к расчетам строительных конструкций, принять участие в семинаре, посвященном расчету зданий и сооружений на статические и сейсмические нагрузки и воздействия.
04 июля 2019
Все новости
Реализация модели контактного слоя при расчете адгезионного соединения с использованием метода конечных элементов
В большинстве опытов по испытанию адгезионных соединений измеряется средняя адгезионная прочность. Данная величина вычисляется как отношение разрушающей нагрузки к площади склейки. Подобный подход подразумевает равномерное распределение касательных напряжений. Исследователи давно обнаружили, что средняя адгезионная прочность соединения является сильной функцией геометрических [1] и физико-механических параметров модели и, следовательно, делает малоинформативными и несопоставимыми экспериментальные данные, выполненные на отличающихся образцах. Малочисленные результаты по измерению касательных напряжений по площади склейки с использованием преимущественно поляризационно-оптических методов [2] показывают, что распределение напряжений является нелинейной функцией. При этом наблюдается концентрация напряжений у торцов модели. В связи с этими фактами возникает необходимость детального изучения напряженно-деформированного состояния адгезионных соединений.
06 июня 2019
Оценка точности нелинейного статического метода анализа сейсмостойкости сооружений

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

06 февраля 2018
Напряженно-деформированное состояние коррозионно - поврежденных железобетонных элементов при динамическом нагружении
С помощью современного программно-вычислительного комплекса  ЛИРА 10.6 выполнена сравнительная оценка напряженно–деформированного состояния  не поврежденного и коррозионно-поврежденного железобетонного элемента при динамическом и статическом нагружении. Проанализировано   влияния ослабленного коррозией бетонного участка сжатой зоны на перераспределение напряжений в сечении.
25 января 2018
Применение технологий BIM при расчете зданий в условиях сложной геотехнической обстановки в связке программ Revit, ЛИРА 10.6 и PLAXIS 3D
В статье рассматривается методика совместной работы ПК ЛИРА 10.6 и PLAXIS 3D посредством API модуля. С позиции инженера-расчетчика рассматриваются возможности передачи моделей между различными программами с применением технологий информационного моделирования.
21 июня 2017
Все публикации
Вебинар "Расчет зданий и сооружений на устойчивость против прогрессирующего обрушения"
На вебинаре рассматривается практическая реализация требований СП 385.1325800.2018, демонстрируются и подробно разбираются расчеты зданий и сооружений различных конструктивных систем на устойчивость против прогрессирующего обрушения
09 августа 2019
Вебинар "Расчет зданий и сооружений в сейсмоопасных районах"
Расчет зданий в сейсмоопасных районах является одним из самых сложных типов расчета и обычно вызывает массу вопросов как по требованиям норм проектирования, так и по применяемым методам расчета. В рамках вебинара мы рассмотрим все расчетные положения, а также «подводные камни» норм проектирования в сейсмических районах.
02 августа 2019
"ПК ЛИРА 10 в задачах". Тема 25. Расчет фильтрации
Учёт фильтрации грунтовых вод является важной задачей во многих областях строительства (геотехника, гидротехника, экология).
В ПК ЛИРА 10.8 реализован тип задачи, который позволяет проводить расчет поля давления и скорости фильтрации грунтовых вод в расчетных схемах с произвольной геометрией для дальнейшего определения напряженно-деформированного состояния от действия вычисленного дополнительного давления.
05 июля 2019
ЛИРА 10. Базовый. Урок 1. Расчет железобетонного каркаса
Цели и задачи: Овладеть навыками создания моделей железобетонных конструкций с нуля. Произвести расчет на основе сочетания нагрузок, научиться анализировать результаты.
03 июля 2019
Все записи вебинаров