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

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.
"ЛИРА софт" открывает новое направление
Выгодные условия на приобретение российского ПО: nanoCAD, Renga, Pilot-BIM и др.
05 апреля 2023
Разбор примеров из практики по расчету на сейсмические воздействия
Рассмотрим реальные примеры уже построенных или проектируемых объектов
22 марта 2023
Приглашаем принять участие в конференции «Градостроительное развитие Алтая». Доступен онлайн формат
Темы конференции: Геологические условия. Высотное строительство. Современные технологии и методы обеспечения качества и безопасности объектов.
03 февраля 2023
Работа клиентской поддержки в новогодние праздники
с 31 декабря по 8 января наш офис будет на каникулах
30 декабря 2022
Все новости
Реализация модели контактного слоя при расчете адгезионного соединения с использованием метода конечных элементов
В большинстве опытов по испытанию адгезионных соединений измеряется средняя адгезионная прочность. Данная величина вычисляется как отношение разрушающей нагрузки к площади склейки. Подобный подход подразумевает равномерное распределение касательных напряжений. Исследователи давно обнаружили, что средняя адгезионная прочность соединения является сильной функцией геометрических [1] и физико-механических параметров модели и, следовательно, делает малоинформативными и несопоставимыми экспериментальные данные, выполненные на отличающихся образцах. Малочисленные результаты по измерению касательных напряжений по площади склейки с использованием преимущественно поляризационно-оптических методов [2] показывают, что распределение напряжений является нелинейной функцией. При этом наблюдается концентрация напряжений у торцов модели. В связи с этими фактами возникает необходимость детального изучения напряженно-деформированного состояния адгезионных соединений.
06 июня 2019
Оценка точности нелинейного статического метода анализа сейсмостойкости сооружений
В статье рассмотрено практическое применение методики нелинейного статического анализа сейсмостойкости зданий и сооружений. Произведен расчет одноэтажной стальной рамы нелинейным статическим и нелинейным динамическим методами. В результате анализа полученных результатов расчета показана значимость высших форм колебаний и необходимость анализа их влияния на реакцию системы.
06 февраля 2018
Напряженно-деформированное состояние коррозионно - поврежденных железобетонных элементов при динамическом нагружении
С помощью современного программно-вычислительного комплекса  ЛИРА 10.6 выполнена сравнительная оценка напряженно–деформированного состояния  не поврежденного и коррозионно-поврежденного железобетонного элемента при динамическом и статическом нагружении. Проанализировано   влияния ослабленного коррозией бетонного участка сжатой зоны на перераспределение напряжений в сечении.
25 января 2018
Применение технологий BIM при расчете зданий в условиях сложной геотехнической обстановки в связке программ Revit, ЛИРА 10.6 и PLAXIS 3D
В статье рассматривается методика совместной работы ПК ЛИРА 10.6 и PLAXIS 3D посредством API модуля. С позиции инженера-расчетчика рассматриваются возможности передачи моделей между различными программами с применением технологий информационного моделирования.
21 июня 2017
Все публикации
Разбор примеров из практики по расчету на сейсмические воздействия
Рассмотрим реальные примеры уже построенных или проектируемых объектов
22 марта 2023
Вебинар: ЛИРА 10 - API. Разработка пользовательских скриптов и плагинов
Вы увидите процесс пошагового создания скриптов и плагинов и поймете, насколько это просто!
02 декабря 2022
Вебинар: От каркаса до расчета в BIM
Приглашаем присоединиться к трансляции, особенно если вам интересна тема взаимодействия Revit и ЛИРА 10
03 ноября 2022
Способы соединения конечных элементов и узлов в ЛИРА 10
Приглашаем на вебинар по применению инструментов соединения элементов и узлов в ЛИРА 10. Будут рассмотрены примеры из инженерной практики: моделирование шарниров, учет несоосности, передача усилий между различными типами КЭ, автоматическое соединение элементов
30 сентября 2022
Все записи вебинаров