Напряженное состояние двухслойной прямоугольной пластинки при сдвиге. Упрощенная двумерная модель

Аннотация


Исследованию напряженного состояния клеевых соединений уделяется значительное внимание, ибо несущая способность конструкции, как правило, определяется прочностью соединений, где напряженное состояние нерегулярно. Большинство существующих математических моделей соединений являются одномерными, т.е. подразумевают равномерное распределение напряжений по ширине соединения. Однако существуют конструкции, для которых классические модели неприменимы. Для расчета напряженного состояния таких соединений необходимо учитывать неравномерность напряжений не только по длине, но и по ширине соединения. Для решения таких задач предложена упрощенная двумерная модель нахлесточного клеевого соединения прямоугольных пластинок. Упрощение состоит в том, что рассматриваются перемещения слоев только вдоль одной из осей координат. Модель является двумерным обобщением классической модели соединения Фолькерсена. Напряжения полагаются распределенными равномерно по толщине слоев, а клеевой слой работает только на сдвиг. Данные упрощения позволили получить аналитическое решение задачи. Решена задача о напряженном состоянии клеевого соединения двух прямоугольных пластин, одна из которых жестко закреплена вдоль одной из сторон, а вторая пластина нагружена неравномерной сдвиговой нагрузкой на противоположной стороне. Задача сведена к системе дифференциальных уравнений второго порядка в частных производных относительно продольных перемещений двух несущих (внешних) слоев. Решение строится при помощи метода разделения переменных и представляет собой функциональный ряд, состоящий из собственных функций. Краевые условия на ненагруженных торцах удовлетворяются точно. Удовлетворение краевых условий на боковых сторонах приводит к системе линейных уравнений относительно неизвестных коэффициентов функционального ряда. Доказана сходимость полученного решения. Решена модельная задача, и проведено сравнение численных результатов с результатами расчетов, выполненных при помощи метода конечных элементов. Показано, что предложенный подход обладает точностью, достаточной для задач проектирования.

Полный текст

Введение Большинство математических моделей клеевых соединений внахлест, позволяющих описать напряженно-деформированное состояние соединения в аналитической форме, являются одномерными [1-5]. Другими словами, предполагают равномерное распределение напряжений по ширине соединения и априорно заданное распределение напряжений по толщине клеевого слоя и соединяемых слоев (как правило, равномерное или линейное). Данный подход успешно применяется для моделирования напряженного состояния двухсрезных клеевых соединений [6, 7]. Наличие симметричных накладок с двух сторон соединения позволяет уменьшить влияние изгибающих моментов в соединяемых слоях на напряженное состояние клеевого слоя. Одним из направлений развития и уточнения математических моделей соединений является изучение распределения напряжений по толщине соединения, т.е. создание двумерных моделей соединения внахлест [8-13]. При этом распределение напряжений по ширине соединения полагается равномерным. Однако при расчете напряженного состояния некоторых конструкций, например соединений силовых элементов или ремонтных накладок с обшивкой, необходимо учитывать неравномерность напряженно-деформированного состояния склеиваемых пластин не только по длине, но и по ширине соединения. Как правило, для исследования двумерного напряженного состояния соединений используются различные численные методы, такие как метод конечных элементов или метод конечных разностей [14-20]. Аналитическое решение задачи о напряженном состоянии соединения в общей двумерной (в плоскости склейки) постановке пока не известно. В связи с этим для решения задач предложено две упрощенные модели, которые позволяют получить решение задачи в аналитической форме: 1) для изучения влияния на напряженное состояние соединения поперечных деформаций, обусловленных коэффициентами Пуассона соединяемых пластин, касательные напряжения в них полагаются равными нулю [21-24]. Приложенная к соединению нагрузка в этом случае полагается равномерно распределенной вдоль боковых сторон пластинок, а соединяемые пластины полагаются абсолютно податливыми на сдвиг; 2) для решения задачи о напряженном состоянии соединения при неравномерной нагрузке перемещения несущих слоев в плоскости соединения в поперечном направлении полагаются равными нулю, т.е. соединяемые пластины полагаются абсолютно жесткими в направлении, поперечном приложенной продольной нагрузке [25]. Гипотеза об отсутствии поперечных перемещений ранее была использована при решении ряда задач для прямоугольных пластин и клеевых соединений (двумерная по толщине модель соединения) [9, 26]. В рамках этой модели также получено аналитическое решение задачи о напряженном состоянии клеевого соединения пластинок разной ширины [27]. В перечисленных выше работах [25, 27] упрощенная двумерная математическая модель клеевого соединения была использована для описания напряженного состояния соединения, вызванного продольной нормальной нагрузкой (растяжение-сжатие). Однако клееные конструкции также могут быть нагружены касательными напряжениями в плоскости соединения [28-31]. Примерами таких конструкций являются соединения полок балок со стенкой балки, соединения панелей обшивки фюзеляжа или крыла самолета и т.д. Целью данной работы является построение аналитического решения и исследование напряженного состояния клеевого соединения внахлест прямоугольных пластин, нагруженных касательной нагрузкой. Данная задача в двумерной постановке решается впервые. 1. Постановка задачи Рассмотрим клеевое соединение двух прямоугольных пластин (), показанное на рис. 1. К противоположным боковым сторонам и приложена сдвиговая нагрузка. Торцы и свободны от нагрузки. Обозначим толщины первого и второго несущих слоев и соответственно. Толщину соединительного слоя обозначим . Считаем, что несущие слои деформируются только в плоскости соединения (плоскости ), клеевой слой работает только на сдвиг, и напряжения равномерно распределены по толщине слоев. Рис. 1. Схема клеевого соединения Fig. 1. A scheme of the glued joint Полагаем, что , и изгибом пластин в плоскости соединения можно пренебречь. Поскольку поперечные (в направлении оси y) перемещения несущих слоев считаем равными нулю, усилия в поперечном направлении и соответствующие касательные напряжения в клее также будут равны нулю. Уравнения равновесия элементов несущих слоев в данном случае имеют вид [16, 25] , , (1) где , - нормальные (в продольном направлении) и касательные усилия в несущем в слое k, k = 1, 2; τ - касательные напряжения в клеевом слое в продольном направлении. Соотношения Коши при условии равенства нулю поперечных перемещений имеют вид , , , (2) где - продольные перемещения слоя k. Напряжения в клеевой прослойке полагаем пропорциональными разности перемещений слоев: , (3) где - модуль сдвига клеевого слоя. Подставив приведенные выше соотношения в уравнения (1), получим систему (4) где ; ; . Краевые условия на торцах имеют вид , т.е. с учетом (2) . (5) На боковых сторонах заданы либо перемещения , (6) либо касательные усилия ; . (7) 2. Построение решения Из первого уравнения системы (3) следует, что . (8) Подставив (8) во второе уравнение системы (4), найдем (9) где , , , , . Обозначим (10) и определим, при каких значениях коэффициентов и уравнение (11) эквивалентно уравнению (9). Подставим (10) в (11) и приравняем коэффициенты при производных в полученном уравнении коэффициентам уравнения (10), получим зависимости , , а также . Последнее условие, очевидно, ограничивает область применения предлагаемого метода решения. Однако если соединяемые материалы одинаковы (что характерно для сэндвич-панелей и ремонтных накладок) или же обладают близкими коэффициентами Пуассона, то предложенный подход оправдан. Обозначим . В таком случае получим . Выражение (10) можно представить в виде . (12) Перемещения второго слоя (8) соответственно представим в виде . (13) Учитывая краевые условия (5), из соотношения (12) следует . (14) Уравнение (11) решим методом разделения переменных. Частные решения уравнения будем искать в виде . После подстановки в уравнение (11) и удовлетворения граничных условий (14) получим Поскольку обе части равенства зависят от разных переменных, то полученное равенство возможно только в случае, когда обе его части постоянны: (15) В результате получим спектральную задачу (16) Система собственных функций данной задачи будет иметь вид где Из второй части равенства (15), учитывая, что получим уравнение , решением которого является семейство функций . Следовательно, общее решение уравнения (11) можно представить в виде (17) где ; , - семейства произвольных коэффициентов. На следующем этапе построения решения необходимо найти общее решение уравнения (12). Решение будем искать в виде линейной суперпозиции решений V и W: (18) где V является общим решением однородного уравнения (12) с однородными граничными условиями (5), (19) Функция W представляет собой частное решение неоднородного уравнения (12), соответствующее заданной функции (17). Решение уравнения (19) также будем искать методом разделения переменных. Искомую функцию представим в виде произведения . Операция разделения переменных приводит к уравнениям , , (20) где - параметр разделения переменных. Общее решение первого уравнения (20) при условии имеет вид . Функция удовлетворяет однородным граничным условиям (5) и не равна нулю тождественно только при . Таким образом, получим . Решение второго уравнения (18) при будет иметь вид . Если же параметр разделения переменных равен нулю, то общие решения уравнений (19) представляют собой линейные функции от координат x и y. Однородные краевые условия (5) в этом случае будут удовлетворены, если Следовательно, решение уравнения (18), которое удовлетворяет краевым условиям (5), можно записать в виде где ; , - семейства произвольных коэффициентов. Функция Ф в уравнении (12) имеет вид (17), после некоторых преобразований частное решение (12) можно представить в виде На основании (12) и (17) получим (21) где ; , . Касательные усилия в несущих слоях (2) имеют вид (22) Краевые условия (7) разложим в ряд Фурье по собственным функциям задачи (16) на интервале , получим , , , , (23) где , . Удовлетворив краевым условиям (7), получим системы линейных уравнений относительно коэффициентов : (24) (25) Решения систем (24) и (25) запишем в виде Из приведенных формул следует, что , где - некоторые линейные выражения с постоянными коэффициентами, зависящие от (22), а . Следовательно , где . Таким образом, гиперболические функции , , и т.д. на интервале ограничены и экспоненциально стремятся к нулю с ростом n. Учитывая сказанное выше, можно сделать вывод, что внутри рассматриваемой области ряды Фурье (20) дважды дифференцируемы и удовлетворяют уравнениям (4). А ряды (21) и (22) на отрезке сходятся равномерно. Касательные напряжения в клее (3) можно представить в виде Можно отметить, что в данной формуле слагаемые, стоящие вне знака суммы, представляют собой известное одномерное решение [30] и описывают напряжения в клеевом слое при равномерной касательной нагрузке. При этом слагаемые, стоящие под знаком суммы, описывают напряжения в клее, вызванные самоуравновешенной нагрузкой, поскольку при . Можно показать, что в глубине области, при удалении от краев и , слагаемые, стоящие под знаком суммы экспоненциально убывают, что совпадает с принципом Сен-Венана. То есть если соединяемые пластины будут достаточно велики, то на удалении от края напряженное состояние соединения будет мало зависеть от конкретного распределения нагрузки по краю и будет определяться лишь суммарной величиной нагрузки . 3. Численный пример Рассмотрим клеевое соединение двух алюминиевых пластинок, имеющих размеры см, см, мм, мм. Толщина клеевой прослойки мм. Упругие характеристики материалов соединения ГПа, ГПа, ГПа. Зададим следующие краевые условия на боковых сторонах склеенных пластин: , ; , где График касательных напряжений в клее приведен на рис. 2. Напряжения на рисунке показаны в безразмерном виде как отношение действующих напряжений τ к гипотетическим напряжениям τ0, которые возникли бы при равномерном распределении нагрузки по всей плоскости клеевого шва, т.е. как отношение τ к напряжениям . Рис. 2. Касательные напряжения в клеевом слое Fig. 2. Shear stresses in the adhesive layer Для верификации предложенной методики проведен расчет напряженного состояния соединения при помощи метода конечных элементов (МКЭ) в пакете COMSOL Multiphysics 5.3. Для расчета использована трехмерная модель, генерация сетки автоматическая, характерный размер элемента 0,03 мм. На рис. 3 приведены графики касательных напряжений в середине толщины клеевого слоя вдоль оси симметрии соединения (), вычисленные при помощи предложенной методики (а), и МКЭ (б). Рис. 3. Касательные напряжения в клее вдоль прямой : а - расчет по предложенной методике; б - МКЭ Fig. 3. Shear stresses in the adhesive layer along the line : а - calculation by the proposed method; b - FEM Рис. 4. Касательные напряжения в клее вдоль прямой : а - расчет по предложенной методике; б - МКЭ Fig. 4. Shear stresses in the adhesive layer along the line : а - calculation by the proposed method; b - FEM На рис. 4 приведены графики касательных напряжений в серединной плоскости клеевого слоя вдоль края , вычисленные при помощи предложенной методики (а), и МКЭ (б). Из графиков видно, что касательные напряжения в клее, вычисленные по предложенной модели, несколько превосходят напряжения, вычисленные при помощи МКЭ. Расчет напряженного состояния данного соединения при помощи МКЭ также показал, что касательные напряжения в клее в поперечном направлении, которые обусловлены поперечными перемещениями слоев и которые в предложенной модели не учитываются, не превосходят 9 % от максимальных касательных напряжений в продольном направлении. Заключение Предложена упрощенная модель клеевого соединения, которая позволяет находить напряженное состояние соединения при произвольном нагружении соединяемых пластин касательными усилиями на боковых сторонах. Получено аналитическое решение задачи и обоснована его сходимость. Решена модельная задача. Расчеты показали, что точность предложенной приближенной методики достаточна для решения многих инженерных задач. Данный подход может быть использован для построения аналитических решений задач о напряженном состоянии клеевых соединений деталей разной ширины; соединений силовых элементов конструкции с обшивкой; соединений, которые имеют дефекты в клеевом слое; соединений ремонтных накладок с обшивкой и других задач, где требуется знать двумерное напряженное состояние клеевого соединения.

Об авторах

С С Куреннов

Национальный аэрокосмический университет им. Н.Е. Жуковского «Харьковский Авиационный Институт»

К П Барахов

Национальный аэрокосмический университет им. Н.Е. Жуковского «Харьковский Авиационный Институт»

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

  1. Analytical models of adhesively bonded joints. Part I: Literature survey / L.F.M. da Silva, P.J.C. das Neves, R.D. Adams, J.K. Spelt // Int. J. Adhes. Adhesiv. - 2009. - Vol. 29. - P. 319-330. doi: 10.1016/j.ijadhadh.2008.06.005
  2. Семин М.И. Расчеты соединений элементов конструкций из полимерных материалов на прочность и долговечность: монография. - М.: Изд-во МАДИ, 2016. - 62 с.
  3. Wang J., Zhang C. Three-parameter elastic foundation model for analysis of adhesively bonded joints // Int. J. Adhes. & Adhes. - 2009. - Vol. 29. - Р. 495-502. doi: 10.1016/j.ijadhadh.2008.10.002
  4. Amidi S., Wang J. Three-parameter viscoelastic foundation model of adhesively bonded single-lap joints with functionally graded adherends // Engineering Structures. - 2018. - Vol. 170. - Р. 118-134. doi: 10.1016/j.engstruct.2018.05.076
  5. Артюхин Ю.П. Напряжения в клеевых соединениях // Исслед. по теор. пластин и оболочек. - Казань: Изд-во Казан. ун-та, 1973. - № 10. - С. 3-27.
  6. Tsai M.Y., Morton J. An Investigation into the Stresses in Double-Lap Adhesive Joints with Laminated Composite Adherends // Int. J. Solids Struct. - 2010. - Vol. 47. - Р. 3317-3325. doi: 10.1016/j.ijsolstr.2010.08.011
  7. Gustafson P.A., Bizard A., Waas A.M. Dimensionless Parameters in Symmetric Double Lap Joints: An Orthotropic Solution for Thermomechanical Loading // Int. J. Solids Struct. - 2007. - Vol. 44. - Р. 5774-5795. doi: 10.1016/j.ijsolstr.2007.01.025
  8. Исследование напряженного состояния в составной пластинке вблизи края линии соединения в зависимости от толщины и материальных параметров соединяющей прослойки / В.М. Пестренин, И.В. Пестренина, Л.В. Ландик, Е.А. Полянина // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2014. - № 1. - C. 153-166.
  9. Васильев В.В., Боков Ю.В. Исследования напряженного состояния клеевого соединения композиционного материала с металлическим листом // Проектирование, расчет и испытания конструкций из композиционных материалов. - М.: Изд-во ЦАГИ, 1979. - Вып. 7. - С. 4-10.
  10. Liu J., Sawa T. Stress analysis and strength evaluation of single-lap band adhesive joints subjected to external bending moments // J. Adhesion Science. - 1999. - Vol. 13. - No. 6. - P. 729-749. doi: 10.1163/156856199X00965
  11. Короткова С.Е. Механика клеевых соединений. - Алчевск, 1998. - 187 с.
  12. Кутьинов В.Ф. Расчет клеевых соединений // Проектирование, расчет и испытания конструкций из композиционных материалов. - М.: Изд-во ЦАГИ, 1979. - Вып. 7. - С. 14-30.
  13. Kurennov S.S. Koshevoi A.G., Polyakov A.G. Through-Thickness Stress Distribution in the Adhesive Joint for the Multilayer Composite Material // Russian Aeronautics (Iz. VUZ) - 2015. - Vol. 58. - No. 2. - P. 145-151. doi: 10.3103/S1068799815020026
  14. Рябенков Н.Г., Артюхин Ю.П. Определение напряжений клея в соединении двух полубесконечных пластин // Исслед. по теор. пластин и оболочек. - Казань: Изд-во Казан. ун-та, 1981. - № 16. - C. 82-90.
  15. Andruet R.H., Dillard D.A., Holzer S.M. Two- and three-dimensional geometrical nonlinear finite elements for analysis of adhesive joints // International Journal of Adhesion & Adhesives. - 2001. - Vol. 21. - Р. 17-34. doi: 10.1016/S0143-7496(00)00024-5
  16. Rapp P. Mechanics of adhesive joints as a plane problem of the theory of elasticity. Part II: Displacement formulation for orthotropic adherends // Archives of Civil and Mechanical Engineering. - 2015. - Vol. 15. - Iss. 2. - P. 603-619. doi: 10.1016/j.acme.2014.06.004
  17. Analysis method for bonded patch repair of a skin with a cutout / Barut A., Hanauska J., Madenci E., Ambur D.R. // Composite Structures. - 2002. - Vol. 55. - Iss. 88. - P. 277-294. doi: 10.1016/S0263-8223(01)00158-1.
  18. Федотов А.А., Ципенко А.В., Лебедев А.И. Численное моделирование клеевого ремонтного соединения // Научный вестник МГТУ ГА. - 2018. - Т. 21(3). - С. 125-138. doi: 10.26467/2079-0619-2018-21-3-125-138
  19. Гришин В.И., Дзюба А.С., Дударьков Ю.И. Прочность и устойчивость элементов и соединений авиационных конструкций из композитов. - М.: Изд-во физ.-мат. лит., 2013. - 272 с.
  20. Design, analysis and performance of adhesively bonded composite patch repair of cracked aluminum aircraft panels / Chukwujekwu Okafor, N. Singh, U.E. Enemuoh, S.V. Rao // Composite Structures. - 2005. - Vol. 71. - P. 258-270. doi: 10.1016/j.compstruct.2005.02.023
  21. Adams R.D., Peppiat N.A. Effect of Poisson’s ratio strains in adherends on stresses of an idealized lap joint // The Journal of Strain Analysis for Engineering Design. - 1973. - Vol. 8. - No. 2. - P. 134-139.
  22. Mathias J.D., Grédiac M., Balandraud X. On the bidirectional stress distribution in rectangular bonded composite patches // Int. J. of Solids and Structures. - 2006. - Vol. 43. - P. 6921-6947. doi: 10.1016/j.ijsolstr.2006.02.016
  23. Transient hygro-thermo-mechanical stresses analysis in multi-layers bonded structure with coupled bidirectional model / R. Kessentini, O. Klinkova, I. Tawfiq, M. Haddar // International Journal of Mechanical Sciences. - 2019. - Vol. 150. - P. 188-201. doi: 10.1016/j.ijmecsci.2018.10.004
  24. Kurennov S.S. An Approximate Two-Dimensional Model of Adhesive Joints. Analytical Solution // Mechanics of Composite Materials. - 2014. - Vol. 50 - No. 1. - P. 105-114. doi: 10.1007/s11029-014-9397-z
  25. Kurennov S.S. A Simplified Two-Dimensional Model of Adhesive Joints. Nonuniform Load // Mechanics of Composite Materials. September. - 2015. - Vol. 51. - Iss. 4. - P. 479-488. doi: 10.1007/s11029-015-9519-2
  26. Васильев В.В., Лурье С.А. О сингулярности решения в плоской задаче теории упругости для консольной полосы // Изв. РАН. МТТ. - 2013. - № 4. - С. 40-49.
  27. Куреннов С.С., Танчик Е.В. Расчет напряженного состояния клеевого соединения прямоугольных пластинок разной ширины // Вестн. Моск. авиац. ин-та. - 2014. - Т. 22, № 2. - С. 162-169.
  28. Kim H., Kedward K. Stress analysis of in-plane shear loaded adhesively bonded composite joints and assemblies // Office of Aviation Research. - Washington, DC, Report No. DOT/FAA/AR017; 2001.
  29. Kim H., Kedward K.T. Stress analysis of adhesively-bonded joints under in-plane shear loading // The Journal of Adhesion. - 2001. - Vol. 76(1). - P. 1-36. doi: 10.1080/00218460108029615
  30. Карпов Я.С. Проектирование деталей и агрегатов из композитов: учебник. - Харьков: Изд-во Нац. аэрокосм. ун-та «Харьковский авиационный ин-т», 2010. - 768 с.
  31. Uslu Uysal M., Güven U. A bonded plate having orthotropic inclusion in the adhesive layer under in-plane shear loading // The Journal of Adhesion. - 2015. - Vol. 92(3). - P. 214-235. doi: 10.1080/00218464.2015.1019064

Статистика

Просмотры

Аннотация - 303

PDF (Russian) - 150

Cited-By


PlumX


© Куреннов С.С., Барахов К.П., 2019

Creative Commons License
Эта статья доступна по лицензии Creative Commons Attribution-NonCommercial 4.0 International License.

Данный сайт использует cookie-файлы

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

О куки-файлах