CC BY
50
Глаголев В.В., Маркин А.А., Фурсаев А.А. Моделирование процесса разделения композита с адгезионным слоем // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2016. - № 2. - С. 3444. DOI: 10.15593/peim.mech/ 2016.2.03
Glagolev V.V., Markin A.A., Fursaev A.A. Separation process modeling of composite with adhesive layer. PNRPU Mechanics Bulletin. 2016. No.2. Рр. 34-44. DOI: 10.15593/perm.mech/2016.2.03
■
ВЕСТНИК ПНИПУ. МЕХАНИКА № 2, 2016 PNRPU MECHANICS BULLETIN
http://vestnik.pstu.ru/mechanics/about/inf/
Б01 10.15593/регт.тесЬ/2016.2.03 УДК 539.375
МОДЕЛИРОВАНИЕ ПРОЦЕССА РАЗДЕЛЕНИЯ КОМПОЗИТА С АДГЕЗИОННЫМ СЛОЕМ
В.В. Глаголев, А.А. Маркин, А.А. Фурсаев
Тульский государственный университет, Тула, Россия
О СТАТЬЕ АННОТАЦИЯ
Получена: 16 марта 2016 г. Рассматривается модель деформирования композиционного материала с тон-
Принята: 25 мая 2016 г. ким адгезионным слоем. Напряженное состояние слоя предлагается рассматри-
Опубликована: 30 июня 2016 г. вать на основе связи средних по толщине слоя напряжений и напряжений по гра-
ницам слоя. Средние деформации слоя выражаются через его граничные перемещения. Использование средних напряжений и деформаций позволяет избежать зависимости напряженно-деформированного состояния слоя от формы его торцевых поверхностей. В рамках малых деформаций получено вариационное условие равновесия тел, соединенных посредством адгезионного слоя. Задача рассматривается в рамках линейной теории упругости. Для сопрягаемых тел законом Гука непосредственно связываются поля деформаций и напряжений. В результате свя-метод конечных элементов занная система вариационных уравнений сводится к уравнениям относительно
полей перемещений в сопрягаемых телах, в том числе и на границах контакта со слоем. В качестве параметра система вариационных уравнений в перемещениях содержит толщину адгезионного слоя. Существенно, что данная система уравнений не является дискретной, так как поля перемещений полагаются непрерывными. Для получения приближенного решения можно использовать различные аппроксимации перемещений. В частности, применялся метод конечного элемента с квадратичной аппроксимацией полей перемещений для случая плоской деформации. Исследовалось влияние характерного размера конечного элемента на сходимость решения. Установлено, что если отношение грани конечного элемента к толщине слоя равно четырем и более, то имеет место численная сходимость. В силу отсутствия сингулярности напряжений в точках сопряжения адгезионного слоя с телами предлагаемый подход позволяет использовать известные локальные критерии разрушения. Проведен анализ возможных видов разрушения композита как за счет разрушения материала слоя, так и посредством разрыва связей между слоем и смежными материалами.
© ПНИПУ
Ключевые слова: композит, подход Нейбера-Новожилова, характерный размер, линейная упругость, вариационное уравнение,
© Глаголев Вадим Вадимович - доктор физико-математических наук, профессор, e-mail: vadim@tsu.tula.ru Маркин Алексей Александрович - доктор физико-математических наук, профессор, e-mail: markin-nikram@yandex.ru Фурсаев Артем Александрович - аспирант, e-mail: artemkajs@mail.ru.
Vadim V. Glagolev - Doctor of Physical and Mathematical Sciences, Professor, e-mail: vadim@tsu.tula.ru Alexey A. Markin - Doctor of Physical and Mathematical Sciences, Professor, e-mail: markin-nikram@yandex.ru Artem A. Fursaev - Postgraduate student, e-mail: artemkajs@mail.ru
SEPARATION PROCESS MODELING OF COMPOSITE WITH ADHESIVE LAYER
V.V. Glagolev, A.A. Markin, A.A. Fursaev
Tula State University, Tula, Russian Federation
ARTICLE INFO
ABSTRACT
Received: 16 March 2016 Accepted: 25 May 2016 Published: 30 June 2016
Keywords:
composite, Neuber-Novozhilov approach, characteristic size, linear elasticity, variational equation, finite element method
The deformation model of a composite with a thin adhesive layer is examined. The consideration of a layer's stress state is based on the relationship between the average stresses by thickness and the stresses on the layer's border. The layer's medium strains are expressed in terms of its boundary displacements. The average stresses and strains are used to avoid the stress-strain state dependence on the shape of end faces. The variational condition for the equilibrium state of two bodies linked through an adhesive layer is obtained within small strains. The problem is considered in the framework of linear theory of elasticity. The Hooke's law relates the strain and stress fields in the matched bodies. As a result, the system of variational equations is reduced to the equations with respect to the displacements fields in the matched bodies including the layer's bounds. The system of variational equations with respect to displacements contains the adhesive layer thickness as a parameter. It is significant that the current equations system is not a discrete one since the displacement fields are supposed to be continuous. Various approximations for displacements may be used to obtain approximate solutions. In particular, the finite element method with a quadratic approximation for displacement fields is used for the case of plane strain. The influence of the characteristic size of a finite element on the convergence of the solution is studied. It is found that the numerical convergence is present when the ratio between the finite element faces and the layer's thickness is four or more. The proposed approach allows to use the well-known local failure criteria under the absence of stress singularity at the points of conjugation of the adhesive layer with the bodies. The analysis of the possible forms of composite destruction due to the destruction of a material layer as well as due to bonds breaking between the layer and adjacent materials is carried out.
© PNRPU
Введение
Моделирование процесса разрушения композиционных материалов, представляющих собой совокупность двух тел, соединенных посредством адгезива конечной толщины, определяется типом нарушения связи между компонентами. Рассматриваются модели адгезионного разрушения в случае, когда адгезив целиком отделяется от материала, или коге-зионного разрушения по массиву адгезива или связанных им тел. Для адгезионных моделей [1-4], как правило, пренебрегают толщиной адгезива, а его механические свойства сводятся к силам взаимодействия склеенных материалов, которые могут иметь разные механические свойства [5-7]. При когезионном разрушении [8-13] рассматривается процесс зарождения трещины в материале с конкретными механическими свойствами. Модель расслоения композита с одинаковыми физическими свойствами, но с разными прочностными характеристиками рассмотрена в работе [14]. В этом плане представляется перспективной разработка таких моделей, которые бы учитывали как механические свойства материалов композита, так и тип разрушения в зависимости от напряженно-деформированного состояния адгезива и соединенных им тел.
В настоящее время наиболее распространенными являются подходы на основе метода конечных элементов [15, 16]. Недостатком данных подходов при прямом моделировании геометрии связующего компонента является наличие точек сингулярности различных видов [17, 18]. Принимая ту или иную форму окончания связующего слоя, в решении задачи приходим к существенно различным распределениям напряжений, в том числе и сингулярным. Отметим, что в случае малости толщины связующего геометрия его окончания в зоне обрыва соединений является неопределенной и связь соответствующей области с заданной геомет-
рией является не вполне корректной. Таким образом, разработка моделей, в которых напряженно-деформированное состояние (НДС) соответствующего композита не зависит от геометрии окончания адгезионного слоя (АС), является достаточно актуальной.
Считаем, что в предлагаемой модели адгезионный слой и соединяемые им материалы могут иметь разные механические свойства. Толщина адгезива мала по сравнению с соответствующими характеристиками соединяемых тел. Для описания НДС адгезионного слоя применим концепцию слоя взаимодействия, используемую в модели трещиноподобного дефекта [19, 20]. Средние напряжения в слое выражаются через компоненты тензора напряжений на границах АС. Соответствующие граничные напряжения естественным образом формируют граничные условия для сопрягаемых со слоем материалов. Граничные напряжения в данном случае будем ассоциировать с адгезионными силами сцепления и достижение ими критического значения трактовать как адгезионное разрушение. Использование средних по толщине напряжений позволяет отказаться от конкретизации геометрии торца АС. Таким образом, задача определения предельного состояния композиционного материала под воздействием внешней нагрузки сводится к нахождению предела прочности в соответствующих материалах и по границам слоя. По достижении той или иной критериальной характеристики разрушение трактуется либо как адгезионное, либо как когезионное.
1. Постановка задачи
На рис. 1 представлено тело, состоящее из трех областей, в общем случае с различными материальными свойствами, где область 3 ассоциируется с клеевой подложкой толщиной 80, размер которой мал по сравнению с толщинами тел 1 и 2. Процесс нагру-
жения предполагаем квазистатическим и изотермическим. Условие равновесия запишем в вариационной форме:
Ц о- -5 г&: = Ц о- -5 г&: + Ц о- -5 г&: + Ц о- -5 г&: = | Р - ЬиШ, (1)
51+2+3 51 $2 53 Ь
где Р - внешняя нагрузка на контуре Ь; о - тензор напряжений; г - тензор деформаций; и - поле перемещений.
Так как система из трех тел находится в равновесии, запишем уравнения равновесия для каждого из них:
Л о- -5 = | Р - 5Ш - | Р + - 5иё1, (2)
Ь1 ь¥С
Ло--5гйУ =|Р-5иё1 Р -5гМ, (3)
$2 ЬЬ
Л о - - 5 гйУ =| Р+ - 5и+а +| Р- - 5и -М, (4)
53 Ьср
где Р+ = а+е + а+2е2, Р = -а-е - а22е2 - векторы напряжений, действующие по верхней и нижней границе АС; а^, а22, а-1, а -2 - граничные напряжения АС; и + , и- - векторы перемещения верхней и нижней границы АС; Ь1, Ь2 - внешние контуры тел 1 и 2 без учета
границ с АС. Торцы слоя считаем свободными от напряжений. При этом постулируется жесткое сцепление между границами АС с областями 1, 2. Кроме того, принимаем, что векторы напряжений на сопряженных границах АС равны и противоположны векторам напряжений сопряженных границ тела.
М
//Р
^ О е р<-►
> 80 0 ; г ' Х1(ГУ » 1
© Т 1 ^ ^ N
В
С
Рис. 1. Нагружение составного тела
Работу внутренних напряжений в АС выразим посредством средних характеристик НДС, полагая в (х1з х2) = в (х1):
■ч
Ц = 501 о--8 г ёх1 =501 о--8 г dx1,
(5)
где о, г - соответственно тензоры средних напряжений и деформаций в слое, а21 (х1 ) =
1 0,580 _ 1 0,580 _ 1 0,580
= Т- | 021 (^ Х2 )^Х2; а12 (х1 )=— | а21 (X1, Х2 )^Х2; 022 (Х1 ) = — | а22 (X1, Х2 ) ^
о 1 21 \ 1 ^ 2 / 2? \1/ О
80 -0,580 80 -0,580
0 -0,580
°11 (Х1 ) = 7Г 1 аП (X1, Х2 )^Х2.
0,580
0 -0,580
Отметим, что в силу симметрии касательных напряжений средние касательные напряжения также симметричны: а21 (х1 ) = 012 (х1 ) .
Средние деформации и перемещения определяем через их граничные значения следующим образом [19, 20]:
Л11 + ( Х- ^ Л11 ( -V
в22 (х1 ) =
^ (Х1) " М- (Х1) ^
, В11 (х1 ) = 0,5
би+( х1) би1 (х1)
V ЙХ1
бх.
1 У
бй1(х1) = и+(х1) -и- (Х1) дй2 (Х1) = 0+би-(х1р
бх2
бх1
= 0,5
V бх1
бх1
У
Ы1 (х1 ) = 0,5 (и+ (х1) + щ (х1)), и2 (х1 ) = 0,5 (и+(х1) + и2 (х1)).
(6)
(7)
(8)
е
Из выражений (7) получим представление средней сдвиговой деформации вдоль слоя:
821 (х1 ) = 812 (х1 ) = 0,5 Г
' du2 ( Xj) du1 ( x
v axi
dx2
= 0,5
ui+ (xi) " ui(xi)
2
f A..+ A.,- i
+ 0,5
^U2+(X1 ) + 5u-(X1 )
(9)
V
V dxi
dx
1 II
Из (4) и (5) с учетом (6) и (9) получаем вариационное условие равновесия АС:
0,55о Jan5
(
5и,+ (х,) (x,)
Л
V dxi
dx
i У
dXi + j°225 (U2 (Xi ) - U- (Xi )) )
+50 Ja2i5
Ui+ (Xi) - u-(Xi) + 0 5 fdu2+(xi ) | du2 (Xi )
dxi dxi
\\
dx, =
(i0)
уу
|a+25u+dxi +|a+i5ui+ dxi - Ja225w2dxi -|a2i5ui dxi.
При отсутствии торцевых нагрузок на АС, интегрируя по частям, имеем
Ja„5 dudfel dxi = 5uf(xi )dxi;
da.
v dxi у
ja2i5 dud^xi) dxi^ItT1 du^(xi^
v dxi у
dx
da 2
dx
(ii) (i2)
где I = 1,2 .
Из (10), с учетом (11) и (12), приравнивая слагаемые при одинаковых вариациях, приходим к представлениям граничных напряжений АС:
+ _ A CS daii + _ A C2 da2i
a2i = a2i - 0,550-, a22 = a22 - 0,550
n M "22
dxi
a- =a2i + 0,55,
daii dxi
a-2 =a22 + 0,55,
dxi
da2i dx
(i3)
(i4)
Подставив (13) в (2), а (14) в (3) и выполнив преобразования, обратные (11), (12), приходим к вариационным уравнениям равновесия для тела 1:
(^5 и Э5 и |
I а11-+ I а21-2^х1 1=1 Р - 5иё1 (15)
{ Эх, { Эх, I {
и тела 2:
d5u d5u
|a--5sds -|a225u2dxi -|a2i5uidxi + 0,5501 Jo,,-l-dxi + Ja2i-- dxi l = J P •5udl. (i6)
■_ d5u2
dxi
dxi
Отметим, что соотношения (15), (16) получены без ограничения на свойства материалов и форму тела.
S
2
Уравнения (15), (16) необходимо замкнуть конкретными определяющими соотношениями. Поведение материала при нагружении определяем законом Гука для тела 1:
о = 20, в; (17)
Р = 3^0, (18)
для тела 2:
о = 202 в; (19)
р = 3К20, (20)
где о - девиатор тензора напряжений; в - девиатор тензора деформаций; К{ - модуль объемного сжатия; 0j - модуль сдвига г = 1,2; р = о --Е.
В материале слоя определяющие соотношения считаем справедливыми для средних по толщине слоя характеристик НДС.
0 = 203 е; (21)
р = 3К3 0. (22)
В результате подстановки в определяющие соотношения (21) и (22) выражений компонент средних деформаций (6), (9) средние напряжения определяются через граничные перемещения и их производные. В этом случае в уравнении (15) будут присутствовать перемещения границы N8 тела 2 (см. рис. 1), а в уравнении (16) - перемещения границы БС тела 1 (см. рис. 1). Таким образом, совместное решение системы (15)-(22) сводится к определению поля перемещений и (хх, х2) в телах 1 и 2.
После определения полей перемещений в телах 1 и 2, в том числе и по границам со слоем, из (21), (22) находим средние напряжения в слое. Зная распределение средних напряжений вдоль слоя, из уравнений (13), (14) находим граничные напряжения по его границе.
2. Метод дискретного решения
Еще раз отметим, что прямое моделирование композитных материалов с явным заданием геометрии торцевой поверхности АС сопряжено с зависимостью НДС в окрестности торцевой поверхности от условий сопряжения слоя и смежных материалов. В частности, если торцевая поверхность является плоскостью, то модель будет сингулярной. Предлагаемый вариант рассмотрения модели композита за счет введения средних по толщине характеристик НДС универсален в плане независимости НДС от геометрии окончания АС.
Для численного решения системы (15)-(22) можно использовать различные дискретные методы, обеспечивающие численную сходимость. В частности, будем использовать МКЭ для дискретизации задачи. Существенно, что построение решения формально не накладывает ограничений на размер конечного элемента, который может быть меньше параметра 80, для обеспечения численной сходимости решения. Однако условия разрушения АС формулируются для 5-элементов слоя размером 80 х80 и на их верхней и нижней границах с использованием процедуры осреднения характеристик НДС, полученных в результате решения системы (15)-(22) и (13), (14):
1 Х1 + °0 1 Х1 + °0 1 Х1 +°0 а =1Г I а =^т I а- =^т I а-^. (23)
°0 х? °0 х? °0 х"
Это является следствием основного физического допущения - разрушение охватывает частицу материала с характерным размером 80 [21, 22]. В качестве критерия разрушения АС будем использовать критерий Кулона, согласно которому разрушение 5-элемента происходит при достижении максимальным главным растягивающим усредненным напряжением критического значения: ашах = ак . Для адгезионного разрушения будем сравнивать граничные напряжения а22, а^2 с адгезионной прочностью на отрыв, а граничные напряжения а 21, а^ - с адгезионной прочностью на сдвиг.
Решение задачи (15)-(22) дает распределение поля перемещений в узловых точках тела 1, в том числе и по границе с АС. После нахождения соответствующего решения возникает задача определения НДС в АС. При решении (15)-(22) используем метод конечного элемента с квадратичным распределением поля перемещений.
3. Результаты численного расчета
В качестве примера использования предлагаемой модели рассмотрим композиционный материал (см. рис. 1) в состоянии плоской деформации, состоящий из адгезионного слоя, близкого по характеристикам к эпоксидной смоле с механическими характеристиками: О = 1,3 -109 Па; К = 1,7 -109 Па; ак = 9 -107 Па - предел прочности; материал 1 и 2
соответствует характеристикам алюминиевого сплава Д16Т: О = 2,8 -1010 Па; К = 6 -1010 Па; ак = 4,2 -108 Па. Адгезионную прочность по границе смола-сплав принимаем равной аа = 2,5 -107 Па на отрыв и сдвиг. Геометрические характеристики образца выбираем следующими: ЛБ = 5-10-3 м; 50 = FN = 10-3 м; МТ = 2-10-2 м; ЛБ = 10-1 м; Б? = 5-10-2 м. На границе ЛБ задаем распределенную нагрузку Р = 1 Па, направленную под углом тс/4 к оси 0X1. Граница TR жестко закреплена от перемещений. Остальная поверхность свободна. В качестве исследуемых параметров задачи будем рассматривать граничные и средние напряжения в адгезионном слое.
На рис. 2 показана вычислительная сходимость решения для отношения максимального главного растягивающего усредненного напряжения на 5-элементе АС к соответствующему значению при грани конечного элемента А, равной параметру 80 в зависимости
от размера А в концевой области слоя относительно параметра 80.
Как видно из графика (см. рис. 2), результат расчета имеет вычислительную сходимость. Дальнейшие вычисления проводились при характерном размере конечного элемента А = 80/4.
На рис. 3 показаны распределения напряжений (23) по структурным 5-элементам слоя. Все напряжения отнесены к граничному напряжению на первом элементе а22 . Графики рис. 3, а определяют граничные напряжения, где кривая 1 задает напряжение а^2, кривая 2 - а22, кривая 3 - а21, кривая 4 - а^ . На рис. 3, б приведены средние напряжения, где кривая 1 задает напряжение а22, кривая 2 - ап, кривая 3 - а21.
а б
Рис. 3. Распределение напряжений в слое
В результате решения задачи на первом 5-элементе АС были получены следующие характеристики: а22 = 6 Па, а^2 = 6,3 Па, а21 = -1,7 Па, а^ = -2,1 Па, ашах = 6,9 Па. Следовательно, разрушение соответствующего композита произойдет при нагрузке Р = 4 • 106 Па посредством адгезионного отслоения по нижней границе слоя в силу достижения предела адгезионной прочности на отрыв. Когезионное разрушение в соответствующих композитах будет наблюдаться в случае, когда адгезионная прочность на отрыв или сдвиг будет приближаться к пределу прочности АС.
Заключение
Использование предлагаемого подхода позволяет естественным образом учитывать характер разрушения при механических воздействиях на композиты с адгезионным слоем. В этом случае задача сводится к решению вариационного уравнения в рамках тех или иных определяющих соотношений.
Работа выполнена при финансовой поддержке РФФИ (грант № 15-01-01875) и Министерства образования и науки РФ (госзадание № 467).
Библиографический список
1. Гольдштейн Р.В., Осипенко Н.М. Отслоение покрытий под действием термоупругих напряжений (Балочное приближение) // Вестник СамГУ. Естественно-научная серия. - 2007. - Т. 54, № 4.- С. 66-83.
2. Sheinman I., Kardomateas G. Energy release rate and stress intensity factors for delaminated composite laminates // International Journal of Solids and Structures. - i997. - Vol. 34. - No. 4. -P. 45i-459. DOI: i0.i0i6/S0020-7683(96)000i8-2
3. Storakers B., Andersson B. Nonlinear plate theory applied to delamination in composites // Journal of Mechanics and Physics solids. - i988. - Vol. 36. - No. 6. - P. 689-7i8. D0I:i0.i0i6/0022-5096(88)90004-X
4. Kulkarni S.V., Frederick D. Propagation of delamination in a layered cylindrical shell // International Journal of Fracture. - i973. - Vol. 9. - No. i. - P. i i3-ii5.
5. Evams A.G., Hutchinson J.W. On the mechanics of delamination and spelling on compressed films // International Journal of Solids and Structures. - i984. - Vol. 20. - No. 5. - P. 455-466. D0I:i0.i0i6/0020-7683(84)900i2-X
6. Устинов К.Б. Об отслоении слоя от полуплоскости; условия упругой заделки для пластины, эквивалентной слою // Изв. РАН. МТТ. - 20i5. - Т. 50, № i. - С. 75-95.
7. Фудзии Т., Дзако М. Механика разрушения композиционных материалов. - М.: Мир, i982. - 232 с.
8. Whitcomb J.D. Finite Element Analysis of Instability related Delamination Growth // Journal of Composite materials. - i98i. - Vol. i5. - No. 5 - P. 403-426. DOI: i0.ii77/002i99838i0i500502
9. The special issue: Cohesive models // Eng. Fract. Mech. - 2003. - Vol. 70. - No. i4. - P. i74i-i987.
10. Williams J.G. On the calculation of energy release rates for cracked laminates// International Journal of Fracture. - i988. - Vol. 36. - No. 2. - P. i0i-ii9. DOI: i0.i007/BF000i7790.
11. Bottega W.J. Structural scale decomposition of energy release rates for delamination propagation // International Journal of Fracture. - 2003. - Vol. i22. - No. i. - P. 89-i00. DOI: i0.i023/B:FRAC.0000005376.59962.f2
12. Кулиев В. Д., Борисова Н.Л. К проблеме разрушения многослойных композитных материалов // Вестн. Чуваш. гос. пед. ун-та им. И.Я. Яковлева. Серия: Механика предельного состояния. - 20i5. - Т. 26, № 4 - С. 63-7i.
13. Yang Q., Cox B.N. Cohesive models for damage evolution in laminated composites // International Journal of Fracture. - 2005. - Vol. i33. - No. 2. - P. i07-i37. DOI: i0.i007/si0704-005-4729-6
14. Астапов И.С., Астапов Н.С., Корнев В.М. Модель расслоения композита при поперечном сдвиге // Механика композиционных материалов и конструкций. - 20i5. - Т. 2i, № 2. -С. i49-i6i.
15. Davila C.G., Camanho P.P., Turon A. Effective Simulation of delamination in aeronautical structures using shells and cohesive elements // Journal of Aircraft. - 2008. - Vol. 42. - No. 2. - P. 663672. DOI: i0.25i4/i.32832
16. Schmidt P., Edlund U. A finite element method for failure analysis of adhesively bonded structures // Int. J. Adhes. Adhes. - 20ii. - Vol. 30. - No. 8. - P. 665-68i. DOI: i0.i0i6/j.ijadhadh.20i0.05.0i2
17. Hildebrand M. Non-linear analysis and optimization of adhesively bonded single lap joints between fibre-reinforced plastics and metals // Int. J. Adhes. Adhes. - i994. - Vol. i4. - No. 4. - P. 26i-267. DOI: i0.i0i6/0i43-7496(94)90039-6
18. He X. A review of finite element analysis of adhesively bonded joints // Int. J. Adhes. Adhes. -20ii. - Vol. 3i. - No. 4. - P. 248-264. DOI: i0.i0i6/j.ijadhadh.20ii.0i.006
19. Glagolev V.V., Glagolev L.V., Markin A.A. Stress-Strain State of Elastoplastic Bodies with Crack // Acta Mechanica Solida Sinica. - 20i5. - Vol. 28. - No. 4. - P. 375-383. DOI: i0.i0i6/S0894-9i66(i5)30023-9
20. Глаголев В.В., Маркин А.А. Нахождение предела упругого деформирования в концевой области физического разреза при произвольном нагружении его берегов. Прикладная механика и техническая физика. - 2012. - Т. 53, № 5. - С. 174-183.
21. Нейбер Г. Концентрация напряжений. - М.: ОГИЗ; Л.: Гостехиздат, 1947. - 204 с.
22. Новожилов В.В. О необходимом и достаточном критерии хрупкой прочности // ПММ. -1969. - Т. 33, № 2. - С. 212-222.
References
1. Goldstein R.V., Osipenko N.M. Otsloenie pokrytii pod deistviem termouprugikh napriazhenii (Balochnoe priblizhenie) [Lamination of coatings under thermoelastic stresses (beam approximation)].
Vestnik Samarskogo gosudarstvennogo universiteta. Estestvenno-nauchnaia seriia, 2007, vol. 54, no. 4, pp. 66-83.
2. Sheinman I., Kardomateas G. Energy release rate and stress intensity factors for delaminated composite laminates. International Journal of Solids and Structures, 1997, vol. 34, no. 4, pp. 451-459. DOI: 10.1016/S0020-7683(96)00018-2
3. Storakers B., Andersson B. Nonlinear plate theory applied to delamination in composites. Journal of Mechanics and Physics solids, 1988, vol. 36, no. 6, pp. 689-718. DOI: 10.1016/0022-5096(88)90004-X
4. Kulkarni S.V., Frederick D. Propagation of delamination in a layered cylindrical shell. International Journal of Fracture, 1973, vol. 9, no. 1, pp. 113-115.
5. Evams A.G., Hutchinson J.W. On the mechanics of delamination and spelling on compressed films. International Journal of Solids and Structures, 1984, vol. 20, no 5, pp. 455 - 466. DOI: 10.1016/0020-7683(84)90012-X
6. Ustinov K.B. On separation of a layer from the half-plane: Elastic fixation conditions for a plate equivalent to the layer. Mechanics of Solids, 2015, vol. 50, no. 1, pp. 62-80. DOI: 10.3103/ S0025654415010070
7. Fujii T., Dzako M. Mekhanika razrusheniia kompozitsionnykh materialov [Fracture mechanics of composite materials]. Moscow: Mir. 232 p.
8. Whitcomb J.D. Finite Element Analysis of Instability related Delamination Growth. Journal of Composite materials, 1981, vol. 15, no. 5, pp. 403-426. DOI: 10.1177/002199838101500502
9. The special issue: Cohesive models. Eng. Fract. Mech., 2003, vol. 70, no. 14, pp. 1741-1987.
10. Williams J.G. On the calculation of energy release rates for cracked laminates. International Journal of Fracture, 1988, vol. 36, no. 2, pp. 101-119. DOI: 10.1007/BF00017790
11. Bottega W.J. Structural scale decomposition of energy release rates for delamination propagation. International Journal of Fracture, 2003, vol. 122, no. 1, pp. 89-100. DOI: 10.1023/ B:FRAC.0000005376.59962.f2
12. Kuliev V.D., Borisova N.L. K probleme razrusheniia mnogosloinykh kompozitnykh materialov [In the issue of multilayer composite materials]. Vestnik Chuvashskogo gosudarstvennogo pedagogi-cheskogo universiteta im. I.Ia. Iakovleva. Mekhanika predel'nogo sostoianiia, 2015, vol. 26, no. 4, pp. 63-71.
13. Yang, Q., Cox, B.N. Cohesive Models for damage evolution in laminated composites. International Journal of Fracture, 2005, vol. 133, no. 2, pp. 107-137. DOI: 10.1007/s10704-005-4729-6
14. Astapov I.S., Astapov N.S., Kornev V.M. Model' rassloeniia kompozita pri poperechnom sdvige [Model of composite exfoliation with transverse shear]. Composite Mechanics and Design, 2015, vol. 21, no. 2, pp.149-161.
15. Davila C.G., Camanho P.P., Turon A. Effective Simulation of delamination in aeronautical structures using shells and cohesive elements. Journal of Aircraft, 2008, vol. 42, no. 2, pp. 663-672. DOI: 10.2514/1.32832
16. Schmidt P., Edlund U. A finite element method for failure analysis of adhesively bonded structures. Int. J. Adhes. Adhes, 2011, vol. 30, no. 8, pp. 665-681. DOI: 10.1016/j.ijadhadh.2010.05.012
17. Hildebrand M. Non-linear analysis and optimization of adhesively bonded single lap joints between fibre-reinforced plastics and metals. Int. J. Adhes. Adhes, 1994, vol. 14, no. 4, pp. 261-267. DOI: 10.1016/0143-7496(94)90039-6
18. He X. A review of finite element analysis of adhesively bonded joints. Int. J. Adhes. Adhes, 2011, vol. 31, no. 4, pp. 248-264. DOI: 10.1016/j.ijadhadh.2011.01.006
19. Glagolev V.V., Glagolev L.V., Markin A.A. Stress-Strain State of Elastoplastic Bodies with Crack. Acta Mechanica Solida Sinica, 2015, vol. 28, no. 4, pp. 375-383. DOI: 10.1016/S0894-9166(15)30023-9
20. Glagolev V.V., Markin A.A. Finding the elastic strain limit at the tip region of a physical cut with arbitrarily loaded faces. Journal of Applied Mechanics and Technical Physics, 2012, vol. 53, no. 5, pp. 784-792. DOI: 10.1134/S0021894412050185
21. Neuber H. Kerbspannunglehre: grunglagen fur Genaue Spannungsrechnung, Springer-Verlag, 1937. 154 p.
22. Novozhilov V.V. O neobkhodimom i dostatochnom kriterii khrupkoi prochnosti [Necessary and Sufficient Criteria of Brittle Strength]. Prikladnaia Matematika Mekhanika, 1969, vol. 33, no. 2, pp. 212-222.
