CC BY
55
УДК 620.171.2
А.В. Ильичев1, А.М. Губин1, А.Р. Акмеев1, Н.В. Иванов1 ОПРЕДЕЛЕНИЕ ОБЛАСТИ
МАКСИМАЛЬНЫХ СДВИГОВЫХ ДЕФОРМАЦИЙ
ДЛЯ ОБРАЗЦОВ УГЛЕПЛАСТИКА ПО МЕТОДУ ИОСИПЕСКУ,
С ИСПОЛЬЗОВАНИЕМ ОПТИЧЕСКОЙ СИСТЕМЫ ИЗМЕРЕНИЙ
DOI: 10.18577/2307-6046-2018-0-6-99-109
Представлены исследования характера изменения области локализации максимальной сдвиговой деформации при испытании образцов углепластика в четырех плоскостях армирования по стандарту ASTMD5379 и ГОСТР 56799. Рассматриваются две марки углепластика на основе углеродной равнопрочной ткани и углеродного ровинга (с ориентацией волокон только в одном направлении). Измерение деформационного поля образца проводили с помощью оптической системы корреляции цифровых изображений. Рассматривается также вопрос влияния различной длины участка измерения деформации как для цифровых экстензометров (метод корреляции цифровых изображений), так и для наклеиваемых тензодатчиков, на достоверность полученных значений.
Ключевые слова: метод корреляции цифровых изображений, сдвиг в плоскости, сдвиг по Иосипеску, углепластик, композиционный материал, ASTMD5379, ГОСТР 56799.
A.V. Ilichev1, A.M. Gubin1, A.R. Akmeev1, N.V. Ivanov1
DEFINITION OF AREA
OF THE MAXIMUM SHEAR DEFORMATIONS
FOR CFRP SAMPLES ON IOSIPESCU METHOD,
WITH USE OF OPTICAL SYSTEM OF MEASUREMENTS
In this work researches nature of change of area of localization of the maximum shear deformation have been provided when testing carbon plastic specimens in four planes on the ASTM D5379 standard (GOSTR 56799). In article there are two kinds of carbon fiber-reinforced polymer. First based on carbon fabric and second based on carbon roving (with orientation offibers only in one direction). Measurements of deformation field of sample were carried out by means of optical system of digital image correlation (DIC). The question of influence of different length of site of measurement of deformation as for digital extensometers (DIC), and for pasted strain gages, on reliability of the received values also is considered.
Keywords: digital image correlation, shear in planes, shear according to Iosipescu, carbon fiber reinforced polymer, composite material, ASTM D5379, GOST R 56799.
^Федеральное государственное унитарное предприятие «Всероссийский научно-исследовательский институт авиационных материалов» Государственный научный центр Российской Федерации [Federal State Unitary Enterprise «All-Russian Scientific Research Institute of Aviation Materials» State Research Center of the Russian Federation]; e-mail: admin@viam.ru
Введение
В условиях современного проектирования и конструирования деталей, узлов и прочих элементов воздушных судов необходимо обладать обширными и достоверными комплексными данными о упруго-прочностных характеристиках применяемых материалов. Для анизотропных слоистых полимерных композиционных материалов (ПКМ)
данные по упруго-прочностным характеристикам, полученным после приложения определенного вида нагрузки и в определенном направлении к различным плоскостям образца, могут на порядок отличаться между собой [1-3].
Одними из важнейших характеристик, используемых в расчете конструкции, является сдвиговая прочность (т) и модуль сдвига (О) полимерного композиционного материала. Прочность при сдвиге характеризует влияние касательных напряжений в материале. Модуль сдвига описывает характер деформационных процессов при воздействии касательных напряжений в зоне упругих деформаций.
Существует ряд работ, связанных с изучением локализации деформационных процессов в рабочей зоне образцов при нагружении - в частности, с применением метода конечных элементов. В работе [4] описывается введение дополнительного расчетного коэффициента для вычисления скорректированного модуля сдвига. Авторами статьи рекомендуется учитывать при расчете модуля сдвига не целое сечение между надрезами, а конкретный объем образца в области аппликации тензодатчиков [4, 5]. Однако тензодатчики могут не захватывать всю область сдвиговой деформации или выходить за ее пределы. Расчет же модели с применением метода конечных элементов требует ввода большого массива данных, уже известных и полученных опытным или теоретическим путем. В работе не описываются критерии определения границ такой области.
В данном исследовании вместо методики конечно-элементного анализа применяется метод оптического измерения деформации образца в процессе его непосредственного испытания. Этот метод основан на корреляции цифровых изображений и позволяет измерять все деформационное поле на поверхности образца. Для обеспечения точного измерения механических характеристик материала на сдвиг, основной задачей является определение локализации области максимальных сдвиговых деформаций и ее объема, а также критерия выбора границы этой области. Это необходимо для получения корректных данных с устройств измерения деформации - экстензометров, тензодатчиков, бесконтактных средств измерения. Такой подход необходим для проверки данных, полученных ранее на материале двух марок углепластика с использованием тензометрических датчиков.
Работа выполнена в рамках реализации комплексных научных проблем 2.2. «Квалификация и исследования материалов» и 13.2. «Конструкционные ПКМ» («Стратегические направления развития материалов и технологий их переработки на период до 2030 года») [6].
Материалы и методы
Наиболее универсальным методом испытаний на плоскостной сдвиг является стандарт ЛБТМ Б5379, именуемый также методом Иосипеску (с 01.01.2017 г. действует ГОСТ Р 56799, гармонизированный со стандартом ЛБТМ Б5379). Данный метод позволяет определить сдвиговые характеристики в шести плоскостях армирования материала (рис. 1). Существует ряд аналогичных методов проведения испытаний для определения прочности при сдвиге. Эти методы подробно изложены и оценены в работе [7] авторов из университета Оттавы. В этой работе метод испытаний по Иосипеску в рейтинге методов на сдвиг в плоскости занимает первое место наряду со стандартом ЛБТМ Б3518 (данный стандарт распространяется только на испытание образцов в плоскостях ху и ух).
В стандарте на испытания для определения деформационных характеристик рекомендуется производить аппликацию тензорезисторов под углами ±45° к оси приложения нагрузки (рекомендуемая схема аппликации, согласно ЛБТМ Б5379 и ГОСТ Р 56799, изображена на рис. 2). Рекомендуемая длина измерительной решетки тензорезистора 1,5 мм, номинальное сопротивление 350 Ом, предельная деформация 3% [8-10].
В ЛБТМ Б5379 не описывается область локализации максимальных сдвиговых деформаций и характер ее изменения. Следует учесть, что для разных плоскостей армирования локализация зоны максимальных сдвиговых деформаций может существенно различаться.
Рис. 1. Ориентация слоев в образцах для определения свойств при сдвиге в шести различных направлениях [9]
Рис. 2. Схема аппликации тензометров на поверхность образца для определения модуля сдвига [9]
Более наглядно изучить область сложнонапряженного состояния и вызванных ею деформаций можно с применением метода корреляции цифровых изображений (КЦИ или DIC - Digital Image Correlation). Принцип работы КЦИ основан на анализе
цифровых изображений и нахождении коэффициента корреляции, который определяется по смещению точек спекл-структуры, предварительно нанесенных на образец [11]. Область распознавания должна быть контрастной, изотропной и нерегулярной. Поверхность образца должна быть однородной черного или белого цвета. Далее на образцы наносят черные метки (точки) на белый фон или белые - на черный фон. В процессе приложения нагрузки и смещения нанесенных точек на поверхности образца данная система позволяет определять поля деформации объекта исследования вплоть до разрушения, выделять и измерять любые участки цифровыми экстензометрами и тензо-датчиками, которые можно выбирать и назначать на любом участке образца, непосредственно в программе обработки скоррелированных изображений [12-14].
Точность измерения перемещения оптической системы Ую-ЗБ:
- по осям х и у 1/50000 от поля обзора;
- по оси г 1/25000 поля.
Локальная разрешающая способность по перемещению составляет 0,01 пикселя (1 мкм - для области размером 100*100 мм). Диапазон измерения деформации - от 0,005 до 2200%.
Для проведения исследований области максимальных сдвиговых деформаций в автоклаве по стандартному режиму изготовили плиты двух марок углепластика: на основе равнопрочного углеродного тканого наполнителя и на основе углеродного ровин-га, ориентированного в одном направлении. В качестве матрицы в углепластиках использовали связующее эпоксидного типа ВСЭ-1212 [15-19]. Толщина плит составляла 2 и 20 мм, что обусловлено схемой вырезки образцов. Из плит затем изготовляли образцы, размеры которых соответствуют стандарту ЛБТМ Б5379 и ГОСТ Р 56799 (рис. 2 и 3). Образцы вырезали в двух плоскостях армирования и для каждой плоскости в двух перпендикулярных направлениях (ху, хг, ух, уг) [20].
Рис. 3. Исходные образцы углепластика на основе углеродной равнопрочной ткани (а) и углеродного ровинга (б), изготовленные по стандарту Л8ТМ Б5379
Данные направления измерения сдвиговой прочности принято обозначать как т12, т13, т21, т23. В данной работе не отражены плоскости (гх, гу), что связано с высокой трудоемкостью изготовления плит из ПКМ. Однако значения сдвиговой прочности и модуля сдвига для направлений гх и гу теоретически должны соответствовать направлению ух в случае материала на основе углеродного ровинга, ориентированного
в одном направлении (для материала на основе равнопрочной ткани значения характеристик будут отличаться).
Испытания ПКМ проводили в соответствии со стандартом ASTM D5379 и ГОСТ Р 56799 c применением оптической системы измерения полей перемещений и деформации Vic-3D, с условием ограничения максимального нагружения образца до 5%-ной деформации сдвига (у) - формулы (1) и (2), согласно стандарту ASTM D5379:
Yi = N+451 + N-451; (1) Г 5%
уа = min < ^ , (2) [у при полной нагрузке
где Yi - деформация сдвига в заданной точке; s+45 - нормальная деформация (угол +45°) в заданной точке; s-45 - нормальная деформация (угол -45°) в заданной точке; yа - предельное значение деформации сдвига.
Рис. 4. Образцы углепластика на основе углеродной равнопрочной ткани (а) и углеродного ровинга (б) после нанесения спекл-структуры для оптической системы измерения поля деформации
Испытания на асимметричный изгиб проводили на аттестованной универсальной статической испытательной машине Walter+Bai LFM-100 с классом точности 0,5 при нагрузке - от 200 Н до 100 кН по ISO 7500-1. Скорость нагружения образца при проведении испытания составляла 2 мм/мин с контролем по деформации 0,01 мин-1. На рис. 4 приведены снимки подготовленных образцов с нанесенной спекл-структурой. Обработку данных проводили с использованием программного обеспечения DION-Pro, VicSnap, Vic-3D.
Результаты и обсуждение
В ходе проведения испытаний образцы доводили до максимальной нагрузки, соответствующей разрушению или достижению предельной деформации (рис. 5). Поле деформации при максимальном напряжении, предшествующее разрушению или достижению 5%-ной деформации сдвига, показано на рис. 6. Поля деформации лежат в плоскости касательных напряжений xy относительно оптического вида образца.
Рис. 5. Разрушение образца в испытательном приспособлении по стандарту Л8ТМ Б5379 и ГОСТ Р 56799
Направление
а)
о)
с 2.0 мм
12
<*У)
. 2.3 мм
Г-
13 (дг) й*1.
п г
■м а 21 (>*) Щ*1
с 2.0мм (2.7 мм
Щ- ■■ ■1 23 (уг) й
2,5 мм
Рис. 6. Общий вид деформационного поля образца в плоскости ху для углепластика на основе углеродного ровинга (а) и углеродной равнопрочной ткани (б)
В ходе анализа полученных изображений определено предельно-напряженное состояние образцов двух марок углепластика. На областях поля деформации образцов определяли равномерную зону максимальной сдвиговой деформации (область, включающую максимальное значение деформации с ограничением по коэффициенту вариации 5%), в которой назначались цифровые экстензометры и тензодатчики. На основе измерений, полученных цифровыми датчиками деформации, определяли максимальное напряженное состояние образца - вплоть до 5%-ной деформации сдвига (согласно стандарту ЛБТМ Б5379 и ГОСТ Р 56799). Последующее возрастание нагрузки и деформации из расчета исключали. На рис. 6 измерена ширина однородной области максимальных сдвиговых деформаций вдоль оси симметрии образца. Области максимальных сдвиговых деформаций имеют ряд различий по геометрической форме и однородности в зависимости от ориентации направления армирующего наполнителя (направления и плоскости вырезки образца). Наименьшая область деформации обнаружена для образца на основе углеродного ровинга, ориентированного в одном направлении с направлением
плоскости 13, что частично объясняется обнаружением в образце дефектов в виде не-сплошностей, которые могли сформироваться на этапе изготовления материала.
Для каждого образца назначали цифровые экстензометры: тип 1 - длиной 6 мм; тип 2 - длиной 3 мм. Центры отрезков (длины) этих цифровых экстензометров совпадают с плоскостью рабочего сечения образца (площадь поперечного сечения образца между надрезами), а также данные экстензометры расположены под углами ±45° к этой плоскости (рис. 7). Результаты измерения представлены на рис. 8, а данные приведены в табл. 1. Образцы углепластика на основе углеродного ровинга со схемами армирования 13, 21 и 23 доводили до разрушения (деформация сдвига менее 5%). Образцы углепластика на основе углеродного ровинга с направлением плоскости 12 и углепластика на основе равнопрочной ткани с направлением плоскости 12, 13, 21 и 23 доведены до нагрузки, ограниченной достижением значения сдвиговой деформации 5%.
Рис. 7. Схематичное изображение назначенных цифровых экстензометров для образцов углепластика на основе углеродного ровинга (а) и углеродной равнопрочной ткани (б)
Направление
а)
о)
X
21
О*)
23 0*)
Рис. 8. Назначенные цифровые экстензометры и поля сдвиговой деформации на поверхности образцов углепластика на основе углеродного ровинга (а) и углеродной равнопрочной ткани (б)
Таблица 1
Результаты измеренной деформации цифровыми экстензометрами
Схема армирования материала и направление
Показания цифровых экстензометров
Тип 1
(длина цифрового экстензометра 6 мм)
£+45, %
£.45, %
Тип 2
(длина цифрового экстензометра 3 мм)
£+45, %
£.45, %
Углепластик на основе углеродного ровинга
12 (xy)
13 (xz) 21 (yx) 23 (yz)
2,269 1,080 0,662 0,554
1,788 0,784 0,775 0,522
2,374 1,025 0,821 0,567
2,185 1,003 0,947 0,558
Углепластик на основе углеродной равнопрочной ткани
12 (xy)
13 ( xz) 21 (yx) 23 (yz)
1,921 2,523 2,464 2,445
1,993 2,113 2,272 2,022
2,404 2,800 2,666 2,717
2,464 2,220 2,457 2,346
На полученных изображениях (рис. 6) определены области максимальных сдвиговых деформаций, соответствующие участку определенного спектра, которые могут быть посчитаны и выражены в точках растрового изображения либо, при введении коэффициента масштаба, в мм .
Общий алгоритм действий выглядит следующим образом:
1. В любом растровом редакторе изображений есть функция выделения области пикселей одного оттенка с заданным допуском отклонения по тону. Воспользовавшись таким инструментом, можно получить зону максимальной сдвиговой деформации одного цвета и использовать ее далее как самостоятельное изображение.
2. Воспользовавшись гистограммой с информацией по каналам RGB, по которой также можно посчитать количество точек в картинке, получим численное выражение площади искомой области в пикселях.
3. Зная исходное разрешение цифрового изображения, можно перевести найденные пиксели в мм2 изображения на экране.
4. Имеющиеся значения площадей зон в мм можно привести к реальному размеру путем отношения характерных линейных размеров образца к отображаемым и введения масштабного коэффициента.
Измеренная площадь и рассчитанный из нее объем области максимальных сдвиговых деформации представлены в табл. 2. Напряжение сдвига в объеме можно описать соотношением максимальной нагрузки, приложенной к образцу (которая соответствует максимальной нагрузке при разрушении образца или достижения 5%-ной деформации сдвига), и общего объема, с наиболее равномерной зоной локализации максимальных сдвиговых деформаций в образце. В результате можно определить и рассчитать всю область с учетом локализации сдвиговой деформации в образце для различных материалов, типов ориентации армирующего наполнителя и плоскости вырезки образцов.
Для оценки полученных данных сопоставляли результаты, приведенные в работе, с ранее полученными значениями. Средние значения получены на шести образцах для каждого направления и плоскости (табл. 3). Определение модуля сдвига проводили на участке в диапазоне от 0,15 до 0,4% деформации сдвига, согласно рекомендациям, описанным в ASTM D5379 и ГОСТ Р 56799. Измерения в ранних исследованиях проводили с использованием тензометрических датчиков с длиной измерительной решетки 5 мм. В табл. 3 также приведены значения, полученные с применением метода КЦИ.
Таблица 2
Объем локализации максимальных сдвиговых деформаций и напряжение сдвига _ в объеме по методу Иосипеску__
Схема армирования Максимальная 2 площадь, мм Объем, мм3 Максимальная Напряжение сдвига в объеме, Н/мм3
материала и направление нагрузка, Н
Углепластик на основе углеродного ровинга
12 (ху) 127,04 476,39 3507 7,4
13 X) 79,52 218,68 1007 4,6
21 (ух) 100,99 369,63 2425 6,6
23 У) 127,61 444,08 1436 3,2
Углепластик на основе углеродной равнопрочной ткани
12 (ху) 132,38 432,90 4090 9,4
13 х) 146,93 424,64 2802 6,6
21 (ух) 111,75 368,77 4277 11,6
23 у) 127,89 337,63 2508 7,4
Результаты определения прочности при сдвиге (т) и модуля сдвига (О
Таблица 3
Схема армирования материала и направление
Данные, полученные с цифровых экстензометров длиной 3 мм (Р1С)
т, МПа
О, ГПа
Данные, полученные с тензодатчиков (длина измерительной решетки 5 мм)
Т0(тенз), МПа
в
(тенз);
ГПа
Углепластик на основе углеродного ровинга
12 (ху)
13 X) 21 (ух) 23 У)
77 30* 53
33
3,6 1,0* 2,1
2,3
90 65 52 30
3,5 3,3 3,3 2,1
Углепластик на основе углеродной равнопрочной ткани
12 (ху)
13 X) 21 (ух) 23 У)
100 78 106 75
4.0 2,7
3,9
3.1
114 85
115
83
4.1
3.2 3,6 3,1
* Структурный дефект образца.
Заключения
По данным представленным в табл. 1, расхождения значений цифровых экстензометров типа 1 и 2 достигают 20%. Поскольку база цифрового экстензометра типа 1 (6 мм) выходит за область уже определенной максимальной локализации сдвиговой деформации, то измеряемые с его помощью значения деформации ниже и в направлении ±45°. В большинстве случаев это может привести к получению завышенного значения прочности при сдвиге, так как достижение ограничения 5%-ной деформации сдвига будет фиксироваться с задержкой. Данные, приведенные в табл. 3, показывают следующее. Так, значения, полученные ранее с применением тензодатчиков с длиной измерительной решетки 5 мм, оказались в некоторых случаях выше значений, полученных с применением методики КЦИ. Это подтверждает рекомендации стандарта ASTM D5379 и ГОСТ Р 56799 в части использования датчика с длиной измерительной решетки 1,5 мм и предельной деформацией не менее 3%. Такие датчики являются менее распространенными и имеют более высокую стоимость. Поэтому иногда этим условием пренебрегают и используют датчики с большей базой тензометрической решетки. По результатам проведенного исследования можно утверждать, что для образцов из углепластика на основе однонаправленного наполнителя и равнопрочной ткани в направлении плоскостей 12, 13, 21 и 23 допустимо использование тензодатчиков с измерительной решеткой 3 мм, что установлено при определении локализации области максимальных сдвиговых деформаций.
Таким образом, по результатам проведенных исследований определена локальная зона максимальных деформаций, вызванная сложнонапряженным состоянием. Установлено, что область локализации сдвиговых деформаций имеет различный характер в зависимости от направления плоскости, подверженной сдвигу к плоскости армирования материала. Рекомендуемые датчики (согласно стандарту ASTM D5379 - длиной 1,5 мм) с учетом правильной аппликации будут находится в области максимальных сдвиговых деформаций для большинства ПКМ. Использование датчиков с большей длиной измерительной решетки или нарушение требований их размещения может привести к включению в измерения области с меньшей долей сдвиговой деформации. Это, в свою очередь, может привести к завышению значений физико-механических характеристик, полученных в ходе испытаний. Определены также критерии нахождения области максимальных сдвиговых деформаций по максимальной деформации, с установлением коэффициента вариации 5%. Рассчитана величина напряжения в объеме локализации области максимальных сдвиговых деформаций, учитывающая действительную область деформирования, а не только общее рабочее сечение образца.
ЛИТЕРАТУРА
1. Каблов Е.Н. Композиты: сегодня и завтра // Металлы Евразии. 2015. №1. С. 36-39.
2. Каблов Е.Н. Контроль качества материалов - гарантия безопасности эксплуатации авиационной техники // Авиационные материалы и технологии. 2001. №1. С. 3-8.
3. Каблов Е.Н. Материалы нового поколения // Защита и безопасность. 2014. №4. С. 28-29.
4. Xavier J.C., Garrido N.M., Oliveira M. et al. A comparison between the Iosipescu and off-axis shear test methods for the characterization of Pinus Pinaster Ait // Composites. Part A: Applied Science and Manufacturing. 2004. Vol. 35 (7-8). P. 827-840.
5. He J., Chiang M.Y., Hunston D.L., Han C.C. Application of the V-notch Shear Test for Unidirectional Hybrid Composites // Journal of Composite Materials. 2002. Vol. 36. No. 23. P. 2653-2666.
6. Каблов Е.Н. Инновационные разработки ФГУП «ВИАМ» ГНЦ РФ по реализации «Стратегических направлений развития материалов и технологий их переработки на период до 2030 года» // Авиационные материалы и технологии. 2015. №1 (34). С. 3-33. DOI: 10.18577/20719140-2015-0-1-3-33.
7. Lee S., Munro М. Evaluation of in-plane shear test methods for advanced composite materials by the decision analysis technique // Composite. 1986. Vol. 17. No. 1. P. 13-22.
8. Ильичев А.В., Раскутин А.Е., Гуляев И.Н. Сравнение геометрических размеров образцов ПКМ, используемых в международных стандартах ASTM и отечественных ГОСТ // Новости материаловедения. Наука и техника: электрон. науч.-технич. журн. 2015. №4. Ст. 05. URL: http://www.materialsnews.ru (дата обращения: 05.06.2018).
9. ГОСТ Р 56799-2015. Композиты полимерные. Метод определения механических характеристик при сдвиге на образцах с V-образными надрезами. М.: Стандартинформ, 2016. 23 с.
10. ASTM D5379/D5379M-05. Standard Test Method for Shear Properties of Composite Materials by the V-Notched Beam Method. West Conshohocken: ASTM International, 2012. DOI: 10.1520/D5379_D5379-12. 13 p.
11. Sutton M.A., Orteu J.-J., Schreier H. Image Correlation for Shape, Motion and Deformation Measurements. Columbia: University of South Carolina, 2009. 364 p.
12. Ильичев А.В., Раскутин А.Е. Исследование влияния концентратора напряжений на напряженно-деформационное состояние углепластика методом корреляции цифровых изображений // Авиационные материалы и технологии. 2014. №3 (32). С. 62-66. DOI:10.18577/2071-9140-2014-0-3-62-66.
13. Луценко А.Н., Одинцев И.Н., Гриневич А.В. и др. Исследование процесса деформации материала оптико-корреляционными методами // Авиационные материалы и технологии. 2014. №S4. С. 70-86. DOI: 10.18577/2071-9140-2014-0-s4-70-86.
14. Третьякова Т.В., Третьяков М.П., Вильдеман В.Э. Оценка точности измерений с использованием видеосистемы анализа полей перемещений и деформаций // Вестник ПНИПУ. Сер.: Механика. 2011. №4. С. 15-28.
15. Раскутин А.Е. Стратегия развития полимерных композиционных материалов // Авиационные материалы и технологии. 2017. №8. С. 344-348. Б01: 10.18577/2071-9140-2017-0-8-344-348.
16. Раскутин А.Е. Российские полимерные композиционные материалы нового поколения, их освоение и внедрение в перспективных разрабатываемых конструкциях // Авиационные материалы и технологии. 2017. №8. С. 349-367. Б01: 10.18577/2071-9140-2017-0-8-349-367.
17. Мухаметов Р.Р., Ахмадиева К.Р., Чурсова Л.В., Коган Д.И. Новые полимерные связующие для перспективных методов изготовления конструкционных волокнистых ПКМ // Авиационные материалы и технологии. 2011. №2. С. 38-42.
18. Мухаметов Р.Р., Ахмадиева К.Р., Ким М.А., Бабин А.Н. Расплавные связующие для перспективных методов изготовления ПКМ нового поколения // Авиационные материалы и технологии. 2012. №8. С. 260-265.
19. Мухаметов Р.Р., Меркулова Ю.И., Чурсова Л.В. Термореактивные полимерные связующие с прогнозируемым уровнем реологических деформационных свойств // Клеи. Герметики. Технологии. 2012. №5. С. 19-21.
20. Вильдеман В.Э., Третьякова Т.В., Третьяков М.П. Экспериментальное исследование закономерностей деформирования и разрушения материалов при плоском напряженном состоянии // Проблемы машиностроения и надежности машин. 2010. №5. С. 106-111.
