Паньков А.А. Диагностирование дефектов пропитки армирующих нитей полимерного композита встроенным оптоволоконным датчиком с распределенной брэгговской решеткой // Вестник Пермского национального исследовательского политехнического университета. Механика. 2020. № 3. С. 60-72. DOI: 10.15593/perm.mech/2020.3.07
Pan'kov A.A. Diagnostics of impregnation defects of reinforcing filaments of polymer composite with built-in fibre-optic sensor with distributed Bragg grating. PNRPU Mechanics Bulletin, 2020, no. 3, pp. 60-72. DOI: 10.15593/perm.mech/2020.3.07
ВЕСТНИК ПНИПУ. МЕХАНИКА
№ 3,2020 PNRPU MECHANICS BULLETIN
http s://ered.p stu .ru/index.php/mechanics/index
Б01: 10.15593/регш.шесЬ/2020.3.07 УДК 539.3; 531.787.5
ДИАГНОСТИРОВАНИЕ ДЕФЕКТОВ ПРОПИТКИ АРМИРУЮЩИХ НИТЕЙ ПОЛИМЕРНОГО КОМПОЗИТА ВСТРОЕННЫМ ОПТОВОЛОКОННЫМ ДАТЧИКОМ С РАСПРЕДЕЛЕННОЙ БРЭГГОВСКОЙ РЕШЕТКОЙ
А.А. Паньков
Пермский национальный исследовательский политехнический университет, Пермь, Россия
О СТАТЬЕ АННОТАЦИЯ
Разработана математическая модель однонаправленного волокнистого полимерного композиционного материала со встроенным в армирующее волокно (нить из элементарных волокон) оптоволоконным датчиком с распределенной брэгговской решеткой с целью диагностирования дефектов пропитки нити - нахождения вероятности дефекта пропитки как относительной длины локальных участков нити без пропитки, т.е. без заполнения связующим пространства между ее элементарными волокнами. Использована методика цифровой обработки спектра отражения по решению интегрального уравнения Фредгольма 1-го рода с целью нахождения искомой информативной функции плотности распределения осевых деформаций по длине чувствительного участка оптоволоконного датчика. Подход предполагает, что оптоволоконный датчик внедрен в композиционный материал на этапе его изготовления, при этом слабоотражательный характер чувствительного участка оптоволокна позволяет осуществить линейное суммирование коэффициентов отражения от его различных локальных участков вне зависимости от их взаимных расположений. Разработан алгоритм численной обработки функции плотности распределения деформаций для нахождения искомой вероятности наличия дефектов пропитки по длине нити. Выявлено, что функция плотности распределения имеет ярко выраженные информативные импульсы, по расположению и величине которых могут быть найдены искомые значения вероятности наличия дефектов пропитки по длине нити. Представлены результаты диагностирования различных значений искомой вероятности наличия дефекта пропитки нити по результатам численного моделирования измеряемых спектров отражения и искомой функции плотности распределения деформаций по длине чувствительного участка оптоволоконного датчика при различных значениях объемной доли нитей, сочетаниях поперечной и продольной нагрузок представительной области однонаправленного волокнистого композиционного материала в сравнении с графиками для случая без нагрузки.
©ПНИПУ
Получена: 09 июня 2020 г. Принята: 31 августа 2020 г. Опубликована: 09 ноября 2020 г.
Ключевые слова:
оптоволокно, брэгговская решетка, распределенный датчик, диагностика деформаций, дефектоскопия, численное моделирование.
© Паньков Андрей Анатольевич - д.ф.-м.н., проф., e-mail:a_a_pankov@mail.ru, : 0000-0001-8477-5206
Andrey A. Pan'kov - Doctor of Physical and Mathematical Sciences, Professor, e-mail: a_a_pankov@mail.ru, : 0000-0001-8477-5206
Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
DIAGNOSTICS OF IMPREGNATION DEFECTS OF REINFORCING FILAMENTS OF POLYMER COMPOSITE WITH BUILT-IN FIBRE-OPTIC SENSOR WITH DISTRIBUTED BRAGG GRATING
A.A. Pan'kov
Perm National Research Polytechnic University, Perm, Russian Federation
ARTICLE INFO ABSTRACT
Mathematical model of unidirectional fibrous polymer composite material with optical fiber sensor built into reinforcing fiber (filament of elementary fibers) with distributed Bragg grating is developed in order to diagnoste defects of filament impregnation - finding probability of impregnation defect as relative length of local sections of filament without impregnation, i.e. without filling binder of space between its elementary fibers. The technique of digital processing of reflection spectrum according to the solution of the integral Fredholm equation of the 1st kind is used in order to find the desired informative function of density of distribution of axial strains along the length of the sensitive section of the fibre-optic sensor. The approach assumes that the optical fiber sensor is embedded in the composite material at the stage of its manufacture, wherein the low-reflective nature of the sensitive portion of the optical fiber allows linear summation of reflection coefficients from its various local portions regardless of their mutual positions. Algorithm of numerical processing of strain distribution density function is developed for finding of sought probability of presence of impregnation defects along filament length. It has been revealed that the distribution density function has pronounced informative pulses, from the location and value of which the sought-after values of probability of presence of impregnation defects along the length of the filament can be found. The results of diagnostics of different values of the sought probability of the filament impregnation defect are presented based on the results of numerical simulation of the measured reflection spectra and the sought function of strain distribution density along the length of the sensitive section of the optical fiber sensor at different values of the volume fraction of the filaments, combinations of transverse and longitudinal loads of the representative domain of the unidirectional fibrous composite material in comparison with graphs for the case without load.
© PNRPU
Received: 09 June 2020 Accepted: 31 August 2020 Published: 09 November 2020
Keywords:
optical fiber, Bragg grating, distributed sensor, strain diagnostics, flaw detection, numerical simulation.
Введение
Сложность многоуровневых иерархических структур современных полимерных композиционных «материал-конструкций» и многообразие технологий их изготовления, различные виды используемых при этом пре-прегов и методов их формования обусловливают многообразие и сложность диагностирования возникающих внутри конструкций дефектов в виде расслоений, непроклеев, локальной пористости, трещин и зон неполного отверждения связующего, подмятий, искривлений, разрывов и дефектов пропитки армирующих волокон при формовании [1-4]. Обзор основных технологий формования изделий авиационного назначения из полимерных композиционных материалов дан в [5, 6], одна из них - «ЯТМ-технология» [7] формования пропиткой под давлением - основана на инжекции связующего (смолы) под давлением в закрытую полость формы с заготовкой из «сухого» армирующего волокна (нитей), при этом смола растекается по всему объему полости формы, вытесняя находящийся там воздух, и пропитывает волокна, заполняя внутриволоконное пространство между элементарными волокнами. При этом существенное значение имеют величины коэффициентов вязкости связующего и проницаемости армирующего волокна, что влияет на скорость технологического про-
цесса, качество пропитки и значение давления, необходимое для получения качественных полимерных композиционных конструкций. Также необходимо учитывать, что связующее может начать отверждение, не успев равномерно распределиться по армирующим волокнам. Вязкость связующего можно варьировать в широком диапазоне через изменения количества вводимого растворителя и/или температуры. В [8] проведены исследования влияния объемной доли дисперсных микросфер на реологические свойства композитного связующего в виде сферопластика с регулируемой вязкостью для заполнения участков многослойных сотовых конструкций на основе эпоксидного олигомера.
Структура имеет определяющее значение при формировании физико-механических, термоупругих, прочностных, усталостных свойств композиционных материалов, и, как следствие, технологические и эксплуатационные дефекты структуры приводят к снижению прочности, надежности, ресурса конструкций. Выявление дефектов структуры особенно актуально для современных высоконагруженных полимерных композиционных конструкций авиационного назначения. Для достижения высоких и стабильных значений характеристик полимерных композиционных материалов и конструкций особое внимание уделяется выявлению, локации, идентификации и уточненному диагностированию раз-
нообразных технологических дефектов, в частности пор, существенным образом влияющих (даже при их малой объемной доле) на стабильность свойств материалов. По месту расположения пор внутри композиционного материала различают поры в связующем, на межфазных границах связующее/армирующие элементы и внутри армирующей нити между элементарными волокнами. Поры в связующем появляются в результате попадания пузырьков воздуха в смолу при ее перемешивании с различными ингредиентами связующего: инициаторами, замедлителями, красителями - и зависит от вязкости и скорости полимеризации связующего. Появление пор внутри армирующей нити обусловливается этими же причинами и, дополнительно, трудностью проникновения связующего в занятое воздухом пространство между элементарными волокнами, наличием или отсутствием соответствующего аппретирующего состава на поверхностях элементарных волокон, оптимального для данного связующего. Поры в композиционном материале появляются также в результате использования армирующего материала повышенной влажности, когда при отверждении в процессе термохимических реакций излишняя влага в армирующем материале испаряется и пары воды, попадая в смолу, образуют в материале пустоты [9, 10]. В обзоре [11] рассмотрены основные современные способы определения объемной доли пор в изделиях из полимерных композиционных материалов с использованием ультразвуковых методов не-разрушающего контроля через измерение коэффициента затухания, скорости продольных ультразвуковых волн, по энергии шумов обратного рассеяния и по прошедшему сигналу в зеркально-теневом режиме. Применение компьютерной томографии при контроле полимерных композиционных конструкций позволяет с высокой точностью и достоверностью определить большинство внутренних дефектов, например микропористости. Совмещение компьютерной томографии и ультразвукового сканирования на некоторых образцах дает расширенную картину наличия дефектов. Методика выявления технологических дефектов и микроповреждений структуры полимерного композиционного материала с использованием рентгенографического подхода даны в [12], где также приведены результаты численного моделирования механического поведения композитной конструкции с дефектом.
Актуальность решения проблемы непрерывного мониторинга состояния ответственных высоконагру-женных полимерных композиционных конструкций обусловила появление и развитие встроенных систем неразрушающего контроля. «Интелектуальные» композиционные материалы и конструкции имеют встроенные функциональные компоненты систем управления и адаптации - сенсоры, определяющие и реагирующие на изменения состояния (характеристик) конструкции под воздействием внешних факторов, и исполнительные (актюаторные) элементы для изменения состояния системы в нужном направлении, которые образуют систе-
му управления с обратной связью и интегрируются в основную конструкцию [13], где рассмотрены концепции создания интеллектуальных полимерных композиционных материалов с оптоволоконными системами встроенного контроля напряженно-деформированного состояния конструкций. «Информкомпозиты» со встроенной системой неразрушающего самоконтроля посредством функциональных сенсорных элементов для регистрации и контроля параметров своего состояния рассмотрены в [14] и являются основой перехода авиационной техники на качественно новый уровень безопасности и надежности ее эксплуатации. Перспективными объектами для использования интеллектуальных материалов являются лопасти несущего винта вертолета для контроля деформированного состояния несущих элементов и управления формой аэродинамической поверхности винта, шумопоглощающие панели для корпуса авиационного двигателя [15-17] и внутренней обшивки салона пассажирского авиалайнера со встроенными в материал исполнительными пьезоэлектрическими элементами и оптоволоконными датчиками с решеткой Брэгга [13]. Оптоволоконные датчики с брэгговскими решетками используют для контроля скорости течения фронта связующего при формовании полимерной композиционной конструкции инфузион-ными методами и для последующего мониторинга этими же датчиками напряженно-деформированного состояния изготовленной конструкции в процессе эксплуатации. Высокотемпературные оптоволоконные брэгговские решетки используют в оптоволоконных сенсорных системах, интегрированных в керамические композиционные материалы аэрокосмического назначения. Интеллектуальные сенсорные системы могут быть встроены в материал для управления процессом «самозалечивания» повреждений материала через контролируемое раскрытие распределенных в нем капсул с лечащими агентами [18]. В [19] исследованы направления развития распределенных волоконно-оптических сенсорных систем и способы обработки информации от них.
Актуальность и подходы к использованию оптоволоконных датчиков в качестве составной части систем мониторинга состояния авиационных конструкций из полимерных композиционных материалов рассмотрены в [20]. Результаты локации зоны внешнего ударного воздействия и индикации возникающих повреждений (расслоений) в углепластиках посредством встроенной волоконно-оптической системы, подтверждающиеся результатами ультразвукового неразрушающего контроля, даны в [21, 22]. В [23, 24] рассмотрены классификация различных дефектов, методы их обнаружения, локации и особенности формы информативных откликов оптоволоконных датчиков, встроенных в полимерные композиционные материалы при сканировании волнами Лэмба. Распределенные оптоволоконные датчики встроенных систем контроля основаны на регистрации вынужденного комбинированного рассеяния:
на основе рассеяния Рэлея, оптовременная рефлекто-метрия на основе рассеяния Римана, оптовременная рефлектометрия на основе рассеяния Бриллюэна и на использовании протяженных слабоотражающих брэг-говских решеток. Использование «чирпированных» брэгговских решеток с линейно изменяющимся периодом и, как следствие, трапециевидным спектром отражения для диагностики неоднородностей и повреждений материала рассмотрено в [25]. Экспериментальные исследования по мониторингу состояния конструкции через измерение информационных откликов от встроенных оптоволоконных брэгговских решеток, вызванных внешними ультразвуковыми воздействиями, даны в [26]. Отметим другие обзорные работы [27-30] по интегрированию оптоволоконных брэгговских решеток в структуру полимерного композиционного материала. Для совместного контроля температуры и напряжений в [31] использованы специально подобранные различные покрытия на чувствительных участках оптоволокна с брэгговскими решетками для разделения упругой и температурной составляющих, измеряемых датчиком «полных» термоупругих продольных деформаций. Результаты математического моделирования и экспериментальных исследований функционирования встроенных в полимерные композиционные материалы оптоволоконных брэгговских решеток для решения различных задач, в частности, диагностирования градиентного поля деформаций в пластине из стеклопластика даны в [32], влияния размеров и упругих свойств оптоволокна и клеевой прослойки на передаточную функцию (показания) закрепленного на поверхности пластины датчика - в [33], обнаружения и идентификации различных типов повреждений в сегменте композитного шпангоута авиационного двигателя - в [34]. Современные технологии производства оптоволокон позволяют варьировать их свойства в широких пределах, удовлетворяя требованиям в зависимости от области применения [35].
Научным заделом решения многих «обратных задач» обработки изображений [36, 37] и косвенных измерений распределенных параметров для случаев, когда непосредственное измерение значений этих параметров затруднительно или невозможно, является теория интегральных уравнений 1-го рода - хорошо развитый раздел функционального анализа [38]. К обратным задачам, в частности, относятся задачи удаленного измерения параметров астрофизических объектов по результатам измерений на поверхности Земли или на искусственных спутниках, геофизические задачи исследования глубоко залегающих объектов по измерениям на земной поверхности; задачи компьютерной томографии внутренних органов по измерениям информативных интегральных ослаблений внешнего рентгеновского излучения при последовательном прохождении луча по различным по плотности тканям сводятся к решению интегральных уравнений. Например, к решению интегральных уравнений Фредгольма 1-го рода сводится
задача удаленного измерения намагниченности объекта по измерениям внешнего магнитного поля [39]. Методы неразрушающего контроля элементов конструкций, во многих случаях, также используют цифровую обработку измеряемых интегральных акустических, электрических или оптических информативных сигналов, что, как следствие, приводит к необходимости решения соответствующих интегральных уравнений. Так, к решению интегральных уравнений Фредгольма 1-го рода приводят задачи нахождения спектров распределений диагностируемых величин, например, давления или температуры по длине оптоволоконных датчиков [40-43] по результатам измерений интегральных спектров для интенсивностей «механолюминесцентного» светового потока [40] или для коэффициентов отражения от «длинной» протяженной брэгговской решетки [41-43] на выходе из оптического волокна. Актуальной остается задача разработки методик цифровой обработки спектров отражения встроенного в полимерный композиционный материал оптического волокна с длинной брэгговской решеткой для мониторинга и диагностирования технологических, эксплуатационных дефектов на структурном уровне протяженного участка материала.
Цель - разработка математической модели диагностирования технологических дефектов пропитки полимерным связующим пространства между элементарными волокнами внутри армирующих нитей в представительном объеме однонаправленного волокнистого композиционного материала посредством встроенного в нить оптоволоконного датчика с распределенной брэг-говской решеткой с использованием методик [41-43] цифровой обработки спектра отражения по решению интегрального уравнения Фредгольма 1-го рода.
1. Модель случайной структуры композита с дефектами пропитки нитей
Рассмотрим математическую модель трехфазного волокнистого композиционного материала V (рис. 1), в котором однонаправленные армирующие волокна в виде цилиндрических нитей из большого числа элементарных однонаправленных волокон (1-я фаза) с круговым поперечным сечением расположены в непрерывном связующем (2-я фаза), при этом связующее пропитывает нить и заполняет пространство между ее элементарными волокнами (с дефектами пропитки в виде пустот - 3-я фаза, т.е. с отсутствием связующего и, как следствие, поперечных связей между элементарными волокнами на локальных цилиндрических дефектных участках V0 нитей). В результате армирующие волокна представляют собой последовательно соединенные в случайном порядке локальные участки нитей с пропиткой и без пропитки элементарных волокон связующим. Рассматриваем случай осесимметричного на-гружения напряжениями с, сгг представительной области композиционного материала V.
й
r \...........4
,/wi-
ЫМ
............N
im ^--Л. I
4—y/
y-f............\v
ущ
D
а б
Рис. 1. Трехуровневая полидисперсная структура (а) и распределение области дефектов пропитки V0 вдоль нити (б) с оптоволоконным датчиком (красная пунктирная линия) в области V под нагрузкой ст, ст=
Fig. 1. Three-level polydisperse structure (a) and distribution of impregnation defect area V0 along thread (b) with fiber optic sensor (red dashed line) in area V under load ст, ст=
Для вычисления эффективных трансверсально-изотропных упругих свойств нити с пропиткой и композита в целом используем известные решения осесим-метричной задачи теории упругости для элементарной двухфазной цилиндрической ячейки полидисперсной структуры [44, 45], при этом для участков нити без пропитки упругие свойства связующего в нити приравниваем нулю. Считаем, что оптоволоконный датчик с распределенной брэгговской решеткой длиной I внедрен в армирующую нить конструкции из полимерного композиционного материала на этапе ее изготовления, при этом часть чувствительного участка длиной 10 расположена внутри дефектной области нити У0 (см. рис. 1), где интенсивности отраженной 1Я и проходящей 11} составляющих входящего /0 в оптоволокно света, коэффициент отражения Я = 1Я / 10. Искомая вероятность p0 наличия дефекта пропитки по длине нити обусловлена особенностями технологии отверждения полимерного связующего.
Для рассматриваемого случая осесимметричного деформирования трансверсально-изотропного упругого материала фаз и композита гидростатическое ст = (стгг + ст00) / 2 и осевое ст22 напряжения
Szz =CTzz / E3 - SCT :
(1)
ст= ki2s+ CU33Szz:
CT zz C1133S + C3333S zz
выражаются через величины относительного изменения объема е = егг + е00 и осевую деформацию е22 и компоненты тензора упругих свойств С , объемный модуль при плоской деформации к12 = (СШ1 + С1122)/2, ось г совпадает с координатной осью г3 (см. рис. 1). При этом для осевой деформации е22 имеем
где величина 5 = 2у31 /Е1 = С1133 /(к12С3333 -С2^), коэффициент Пуассона у31 , модули Юнга Е13 вдоль осей г13 с учетом
Е = ЕЕ = С - С2 / к (2)
-^2 ' 3 3333 1133 12 • У^-1
Решения для константы к*2 и компонент С*133, С3*333 тензора эффективных трансверсально-изотропных упругих свойств однонаправленного волокнистого композита С* для случаев (рис. 1, б), когда композит - волокна нити в связующем (с пропиткой) С*' и волокна нити без связующего (без пропитки) С*', и случая, когда композит - нити с эффективными свойствами (с пропиткой или без пропитки) в связующем (матрице) С* , найдены из решения осесимметричной задачи теории упругости для элементарной двухфазной цилиндрической ячейки полидисперсной структуры [44, 45] в виде
k12 k(2)12 + V1k12ai1, С1133 = C(2)1133 + V1C1133а11
C3333 =< C3333 > +V1C1133a12
(3)
с использованием общего решения для радиальных перемещений в волокне и(Г)г = Сг при г е (0; а), где константа С = (а11е* + а12е*) / 2, осевая деформация е* = е33, коэффициенты
а11 = (к(2)12 + °(2)12 ) / О) , а12 = -(1 - У1)С1133 / а0 ,
а0 = к(1)12 - V1к12 + С(2)12 , (4)
а
с учетом обозначений разностей модулей фаз:
к12 = к(№ - ^(2)12 , С1133 = С(1)1133 - С(2)1133, модуль сдвига
( = С
(2)12 (2)1212
Для случая отсутствия пропитки нити (в области V, на рис. 1, б) единственная ненулевая компонента С*33 трансверсально-изотропного тензора C*" эффективных упругих свойств нити (как системы невзаимодействующих между собой однонаправленных по оси г3 волокон нити) С3*3'33 = Е*'' = ^£(1)3 равна эффективному модулю Юнга нити Е*' (2), который рассчитывается через продольный модуль Юнга Е(1)3 и относительное объемное содержание у = у ф 1 волокон в нити и известен в механике композитов как «решение Фойгта».
2. Цифровая обработка информативного спектра коэффициента отражения
Диагностирование дефектов пропитки армирующей нити полимерного композиционного материала (см. рис. 1) посредством встроенного в нить оптоволоконного датчика с распределенной брэгговской решеткой осуществим по результатам цифровой обработки информативного спектра коэффициента отражения Я(Х) оптоволокна через нахождение и анализ «спектра деформаций» / (е) - функции плотности вероятностей случайной величины осевой деформации е, на произвольном
элементарном участке брэгговской решетки оптоволокна. Считаем, что случайная структура и деформационные поля внутри полимерного композиционного материала (см. рис. 1) в целом и осевые деформации е, на представительном (длиной I) чувствительном участке оптоволокна в частности обладают свойством статистической однородности и эргодичности [46], что позволяет при осреднении, например, осевых деформаций е, по длине , чувствительного участка оптоволокна применить подходы «геометрической вероятности». В общем случае, когда пространственную ориентацию чувствительного прямолинейного участка оптоволокна задает некоторый единичный направляющий вектор d с координатами С1, значения продольных (по направлению d) деформаций е, на элементарном чувствительном участке
е, =
связаны с компонентами ев тензора деформаций г на
этом участке оптоволокна. Для частного случая (см. рис. 1), когда ориентация чувствительного участка оптоволокна вдоль координатной оси 7, имеем равенство е, = еш (1).
Искомая функция плотности распределения /е (е) осевых деформаций е, по длине чувствительного уча-
стка , оптоволоконного датчика с распределенной сла-боотражающей брэгговской решеткой может быть найдена [41] из решения интегрального уравнения Фред-гольма 1-го рода
ш а) = } к а - 8)/(^
(5)
по результатам измеряемых значений спектра коэффициентов отражения ш(X) = Я("к) на «входе-выходе» оптоволокна для входящих в оптоволокно световых волн длиной X, где искомая функция /(8) - / (е)е=х, ядро
Фредгольма
к (X) - ^ «0(Х, I')
Х0 ^ д1'
(6)
\1'=0
рассчитывается через считающийся известным спектр коэффициентов отражения Я0(Х, I') для начальной однородной распределенной брэгговской решетки с некоторой малой длиной , . Интегральное уравнение (5) получено на основе представления результирующего коэффициента отражения Я для неоднородно деформированной по длине распределенной слабоотражаю-щей брэгговской решетки в виде интеграла
Я(Х) = | с1Я =| К0(Х - X0e)С/e
(7)
элементарных коэффициентов отражения сЯ по всей длине чувствительного участка оптоволокна I для некоторой световой волны с длиной X, где функция коэффициентов
вд =
дЯ0(Х,,') д,
с учетом равенства нулю значения Я (X, 0) = 0, при этом считаем, что для случая малых значений , имеем прямую пропорциональность Я^,I' ) = К^)!' для зависимости значений Я^^, I' ) от длины сенсорного участка , , и при однородной осевой деформации е сенсорного участка , функция коэффициентов отражения - это результат «смещения» функции Я^^, I') на величину X0e по оси абсцисс X. В (7) элементарную длину С1е чувствительного участка (или суммарную длину участков) оптоволокна с осевой деформацией ее (е,е + Се) представляем в виде
С1е = / (е)Се
(8)
на основе двойного равенства для «геометрической вероятности» СР = С1е /1 = /е (е)Се, где искомая функция плотности распределения / (е) неоднородных осевых деформаций е по всей длине I чувствительного участка оптоволокна, при этом выполняется равенство
le (e) = l J fe (e ')de ',
—да
область значений ls e (0; l). В результате подстановки (8) в (7) получим уравнение Фредгольма 1-го рода
да
R(k) = l J Ko(X —e)f (e)de
—да
или в канонической форме (5). Ранее в [43, 44] методика цифровой обработки спектра отражения R(X) по решению интегрального уравнения Фредгольма 1-го рода (5)-(8) [41] была применена для решения задач диагностирования объемного деформированного состояния [44] и локальной микропористости [43] внутри полимерного материала в окрестности внедренного оптоволоконного датчика с распределенной брэгговской решеткой.
Функция плотности распределения f (e) продольных деформаций е, вдоль чувствительного участка оптоволокна может быть найдена из решения интегрального уравнения Фредгольма 1-го рода (5), (6). При этом функция f (e) будет иметь два характерных максимума в виде локальных «левого» и «правого» импульсов при значениях аргумента e1,0 соответственно. В предположении об однородности осевых деформаций e1,0 по длине каждого из двух типов участков нити: без дефекта e1 и с дефектом e0 пропитки (рис. 1, б) - решение для функции f (e) имеет вид
fe (e) = Ê Pf S(ef — e), (9)
f =0
где дельта-функция Дирака 5(e).
Искомое численное значение вероятности наличия дефекта пропитки нити P0 как величины относительного объема области нити без пропитки V0 в объеме нити (см. рис. 1, б) находим интегрированием
Po =J fe (e)de (10)
eo
по значениям осевых деформаций пор e0 , т.е. по абсциссам ширины локального правого импульса функции fe (e). Выполняются равенства для нахождения вероятности пропитки нити
P = 1 — P0 =J fe (e)de (11)
e1
через значение p0 (10) и/или через интеграл по абсциссам e1 ширины локального левого импульса функции
да
f (e) с учетом J fe (e)de = 1.
—да
3. Результаты численного моделирования
Численное моделирование процесса диагностирования дефектов пропитки армирующих нитей осуществим
для математической модели представительного объема трехфазного однонаправленного волокнистого композиционного материала (рис. 1, а) с полной корреляцией случайных событий появлений дефектов пропитки в сечениях всех нитей в каждом поперечном сечении композита с целью более «яркого» проявления анализируемых информативных эффектов. При этом возможная корреляция взаимного расположения и размеров дефектов пропитки по длине нитей в этой модели не учитывалась. В разложении (9) осевые деформации е10 на цилиндрических участках нити с полной пропиткой (без дефекта) е1 и с полным отсутствием пропитки (с дефектом) е0 рассчитываем по формуле е0 =е22 (1) через значения величин к*2 , С*133, С3*333 (3), (4), где численные значения относительной объемной доли нитей у = 0,9; 0,95; 0,98; 0,99. Индекс «2» у
величин к(1,2)12 , С(1,2)1133 , С(1,2)3333 , ((2) в (3), (4) указывает на принадлежность к эпоксидному связующему, при этом индекс «1» для случая нахождения деформации е1 - принадлежность к материалу с тензором упругих свойств С*', а для случая нахождения деформации е0 - принадлежность к материалу с тензором
упругих свойств С*'', относительное объемное содержание элементарных волокон в нити у' = у''= 0,6. Упругие свойства изотропных фаз: полиэфирных волокон Е = 20 МПа, у1 = 0,31 и эпоксидного связующего Е2 = 3 ГПа, у2 = 0,4 заданы через модуль Юнга Е и коэффициент Пуассона V. Значения независимых компонентов тензора эффективных трансверсально-изотропных упругих свойств нити с пропиткой -полиэфирные волокна/эпоксидное связующее: С*; 11 = 0,806 ГПа, С*122 = 0,323 ГПа, С*133 =
= 0,448 ГПа, С3*333 = 1,568 ГПа, С*323 = 0,275 ГПа - были рассчитаны по решению (3); параметры используемой брэгговской решетки описаны в [47].
Результаты численного моделирования даны на рис. 2, 3 в виде графиков спектров коэффициента отражения Я и функций плотности распределения деформаций /е (е) при различных значениях искомой вероятности наличия дефекта пропитки нитей р0 в однонаправленном волокнистом композиционном материале под равномерной поперечной нагрузкой а с комбинацией различных значений осевой нагрузки а22 при объемной доле нитей у = 0,99 (см. рис. 2) и при различных значениях объемной доли нитей у = 0,9, 0,95, 0,98, 0,99 под фиксированной нагрузкой а, а22 и искомым значением вероятности наличия дефекта пропитки р0 (см. рис. 3) в сравнении с графиками для случая без нагрузки.
0,02 -
0,015 -
0,01 -
0,005
0,80
0,90 1,00 1,10 1,20
1
0,8 0,6 -0,4 -0,2 -
0 +0 0,00
f ,103
0,02
0,04 б
0,06
0,08 e„
0,02
0,015
0,01 -
0,005
0
0,80 0,90 1,00 1,10 1,20 1,30 в
1
0,8 0,6 0,4 -0,2 -
0 4® 0,00 e,
f ,103
0,05
0,10
0,15 e„
0
R л ±
¡Л f i
b
i* -J. W c
--Jffm.............mVi g
1
0,8 0,6 0,4 -0,2 -
0,80 0,90 1,00 1,10 1,20 1,30 1,40 1,50 д
0 -ю-0,00 e1
f ,103
0,10
0,20
0,30
e
Рис. 2. Спектры коэффициента отражения R и плотности распределения деформаций f (e) при вероятности наличия дефекта пропитки p0 = 0,2 (о), 0,7 (◊), 0,99 (□) под нагрузкой ст = 2 МПа, CTzz = 5 МПа (а), (б), 8 МПа (в), (г), 15 МПа (д), (е), объемная доля нитей v1 = 0,99, без нагрузки ст = ст= =0 (а)
Fig. 2. Spectra of coefficient of reflection R and distribution density of strains f (e) at the probability of existence of defect of impregnation p0 = 0.2 (о), 0.7 (0),0.99 (□) under loading ст = 2MPa , ст= = 5MPa (a), (b), 8МПа (с), (d), 15MPa (e), f), a volume fraction of threads v1 = 0.99 , without loading ст = ст= =0 ( <=i)
e
0
e
а
e
г
e
е
0,02 -
0,015 -
0,01 -
0,005
0
0,8 -,
0,6 -
0,4 -
0,2
f ,103
0,80 0,90
1,00
1,10
1,20
1,30
0,00
0,05
0,10
0,15
0,20
а б
Рис. 3. Спектры коэффициента отражения R (а) и плотности распределения деформаций f (e) (б) при объемной доле нитей v1 = 0,9 (о), 0,95 (◊), 0,98 (А), 0,99 (□) под нагрузкой ст = 2 МПа, ст= = 10МПа, вероятность наличия дефекта
пропитки p0 = 0,7, без нагрузки ст = ст= =0 ( i=)
Fig. 3. Spectra of reflection coefficient R (a) and density of strains distribution f (e) (b) at volume fraction of threads v1 = 0.9 (о), 0.95 (◊), 0.98 (А), 0.99 (□) under load ст = 2MPa , ctzz = 10MPa , probability of presence of impregnation defect p0 = 0.7 ,
without load ст = ст= =0 (a)
По графикам на рис. 3, б, г, е находим численные значения деформаций е1 = 2,77 10-3, е0 = 8,12 10-2 (рис. 3, б), 5,21-10-3, 1,53-10-1 (см. рис. 3, г), 1,09-10-2, 3,20 -10-1 (рис. 3, е) для различных значений нагрузки ст 22 = 5 МПа, 8 МПа, 15 МПа. По графикам на рис. 3, б находим численные значения деформации е1 = 6,05-10-3, 6,46-10-3, 6,74-10-3, 6,84-10-3 для различных значений v1 = 0,9...0,99. Искомые численные значения вероятности наличия дефекта пропитки нити р0 находим интегрированием (10) по абсциссам локального правого импульса функции /е (е) осевых деформаций пор е0 на рис. 2, б, г, е, рис. 3, б.
Заключение
Разработана математическая модель однонаправленного волокнистого полимерного композиционного материала со встроенным в армирующее волокно (нить из элементарных волокон) оптоволоконным датчиком с распределенной брэгговской решеткой (см. рис. 1) с целью диагностирования дефектов пропитки нити -нахождения вероятности дефекта пропитки р0 как относительной длины 10 /1 локальных участков нити без пропитки 10, т.е. без заполнения связующим пространства между ее элементарными волокнами. Использована методика [41-43] цифровой обработки спектра отражения Я(Х) по решению интегрального уравнения
Фредгольма 1-го рода с целью нахождения искомой информативной функции плотности распределения /е (е) осевых деформаций е по длине чувствительного
участка I оптоволоконного датчика. Подход предполагает, что оптоволоконный датчик внедрен в композиционный материал на этапе его изготовления, при этом слабоотражательный характер чувствительного участка оптоволокна позволяет осуществить линейное суммирование коэффициентов отражения от его различных локальных участков вне зависимости от их взаимных расположений. Разработан алгоритм численной обработки функции плотности распределения деформаций /е (е) для нахождения искомой вероятности р0 наличия дефектов пропитки по длине нити (9)-(11). Выявлено (9), что функция плотности распределения /е (е) имеет ярко выраженные информативные импульсы (см. рис. 2, б, г, е, рис. 3, б), по расположению и величине которых могут быть найдены искомые значения вероятности р0 наличия дефектов пропитки (10) по длине нити. Представлены результаты диагностирования различных значений искомой вероятности наличия дефекта пропитки нити р0 по результатам численного моделирования измеряемых спектров отражения Я (Я,) и искомой функции плотности распределения деформаций /е (е) по длине чувствительного участка оптоволоконного датчика (см. рис. 2, 3) при различных значениях объемной доли нитей V , сочетаниях поперечной ст и продольной ст22 нагрузок представительной области однонаправленного волокнистого композиционного материала в сравнении с графиками для случая без нагрузки.
Математическая модель (см. рис. 1) и методика диагностирования дефектов пропитки армирующих нитей (5)—(11) внутри полимерного волокнистого композиционного материала для случая его осесимметричного
0
e
0
e
e
0
деформирования могут быть адаптированы как на другие «простые» случаи упругого деформирования, например, изгиб и/или кручение стержней с волокнистой структурой с использованием известных аналитических решений [48, 49] для различных вариантов (осевой, окружной, спиральной) ориентации чувствительного участка встроенного оптоволокна относительно оси стержня, так и на случаи «сложных» нагружений элементов конструкций из полимерных композиционных материалов, например, композитного шпангоута авиационного двигателя [34] с использованием известных численных методов и пакетов прикладных программ механики де-
Библиографический список
1. Воробей В.В., Маркин В.Б. Контроль качества изготовления и технология ремонта композитных конструкций. -Барнаул: Изд-во МЦ ЭОР, 2015. - 181 с.
2. Бормотов А.Н. Декомпозиция систем и иерархические структуры показателей качества композитов // Современные наукоемкие технологии. - 2016. - № 9-2. - С. 196-203.
3. Карташова Е.Д., Муйземнек А.Ю. Технологические дефекты полимерных слоистых композиционных материалов // Известия вузов. Поволжский регион. Технические науки. -2017. - № 2. - С. 79-89.
4. Мурашов В.В., Румянцев А.Ф. Дефекты монолитных деталей и многослойных конструкций из полимерных композиционных материалов и методы их выявления. Часть 1. Дефекты монолитных деталей и многослойных конструкций из полимерных композиционных материалов // Контроль. Диагностика. - 2007. - № 4. - С. 23-32.
5. Методы формования изделий авиационного назначения из ПКМ (обзор) / Е. Д. Колпачков, А. П. Петрова, А. О. Курносов, И.И. Соколов // Труды ВИАМ. - 2019. - № 11. - С. 22-36.
6. Куприянова О.А., Сержантова М.В., Банщикова М.Н. Параметры контроля качества производства композитных рефлекторов антенн // Журнал Сибирского федерального университета. Серия: Техника и технологии. - 2015. - № 3. - С. 324-330.
7. Коваленко А.В. Исследование свойств связующего для формования изделий методом пропитки под давлением // Труды ВИАМ. - 2015. - № 1. - С. 30-37.
8. Сферопластик с регулируемой вязкостью для заполнения участков сотовых конструкций / А.В. Коваленко, Н. К. Сидельников, И. И. Соколов, К. О. Тундайкин // Труды ВИАМ. - 2019. - № 11. - С. 37-43.
9. Гайдачук О.В., Кондратьев А.В., Омельченко К.В. Анализ технологических дефектов, возникающих в серийном производстве интегральных авиаконструкций из полимерных композиционных материалов // Авиационно-космическая техника и технология. - 2010. - № 3. - С. 11-20.
10. Гайдачук А.В., Сливинский М.В. Технологические дефекты сотовых конструкций из полимерных композиционных материалов // Авиационно-космическая техника и технология. - 2005. - № 4. - С. 5-8.
11. Диков И.А., Бойчук А.С. Способы определения объемной доли пор в полимерных композиционных материалах с помощью ультразвуковых методов неразрушающего контроля (обзор) // Труды ВИАМ. - 2017. - № 2. - С. 80-95.
12. Исследование механического поведения композитной конструкции при разработке методики достоверного выявления дефектов рентгенографическим методом неразрушающего кон-
формируемого твердого тела. Планируется проведение экспериментальных исследований для отработки предложенной методики диагностирования дефектов пропитки армирующих нитей внутри полимерного композиционного материала.
Благодарность
Результаты получены при выполнении государственного задания Министерства науки и высшего образования Российской Федерации на выполнение фундаментальных научных исследований (проект № Б8КМ-2020-0026).
троля / А.Н. Аношкин, В.Ю. Зуйко, В.М. Осокин, К.А. Пеленев, А.А. Третьяков // Математическое моделирование в естественных науках. - 2017. - № 1. - С. 169-172.
13. Федотов М. Ю. Концепции создания и тенденции развития интеллектуальных материалов (обзор) // Авиационные материалы и технологии. - 2015. - № 1. - С. 71-80.
14. Встроенный контроль: от датчиков до информкомпо-зитов / Е.Н. Каблов, Г.А. Морозов, Д.В. Сиваков, И.Н. Гуляев // Авиационные материалы 75 лет. Избранные труды. - 2007. -С. 331-342.
15. Патент ЯИ № 2643692. Волоконно-оптический датчик объемного напряженного состояния / Паньков А.А., опубл.: 05.02.2018, Бюл. № 4; заявка № 2017111405 от 04.04.2017 г.
16. Патент ЯИ № 2604174. Система гашения спектра акустических шумов / Паньков А.А., Аношкин А.Н., Писарев П.В., опубл. 10.12.2016 Бюл.№ 34; заявка ЯИ № 2015144968 от 19.10.2015 г.
17. Патент ЯИ № 2697168. Лопасть воздушного винта с управляемой геометрией профиля / Паньков А.А., Аношкин А.Н., Писарев П.В., опубл.: 12.08.2019 Бюл. № 23, заявка № 2018140318 от 14.11.2018 г.
18. Колобков А.С., Малаховский С.С. Самозалечивающиеся композиционные материалы (обзор) // Труды ВИАМ. -2019. - № 1. - С. 47-54.
19. Сорокин К.В., Мурашов В.В. Мировые тенденции развития распределенных волоконно-оптических сенсорных систем (обзор) //Авиационные материалы и технологии. -2015. - № 3. - С. 90-94.
20. Качура С.М., Постнов В.И. Перспективные оптоволоконные датчики и их применение (обзор) // Труды ВИАМ. -2019. - № 5. - С. 52-61.
21. Возможности контроля внешних механических воздействий волоконно-оптической системой диагностики, встроенной в углепластики различных типов / М. Ю. Федотов, О.Н. Будадин, С.А. Васильев, О.И. Медведков, С.О. Козельская // Контроль. Диагностика. - 2019. - № 3. - С. 38-47.
22. Федотов М.Ю. Развитие технологии оптического контроля конструкций из ПКМ волоконно-оптическими датчиками: автореф. дис. ... канд. техн. наук. - М., 2019. - 24 с.
23. Дифференцирование дефектов в ПКМ по отклику волоконно-оптических сенсоров (обзор). Часть I / Л.А. Каша-рина, В.В. Махсидов, О.И. Смирнов, И.А. Рузаков // Труды ВИАМ. - 2019. - № 2. - С. 97-104.
24. Кашарина Л.А., Махсидов В.В. Дифференцирование дефектов в ПКМ по отклику волоконно-оптических сенсоров (обзор). Часть II // Труды ВИАМ. - 2019. - № 6. - С. 43-50.
25. Рузаков И.А. Мониторинг деформационного состояния элементов конструкции из ПКМ на основе волоконно-оптических датчиков (обзор) // Труды ВИАМ. - 2019. - № 4. -С. 88-98.
26. Обнаружение ультразвуковых воздействий с применением волоконных решеток Брэгга / А.А. Власов, А.С. Алей-ник, Ф.А. Шуклин, А.Н. Никитенко, Е.А. Моторин, А.Ю. Ки-реенков // Научно-технический вестник информационных технологий, механики и оптики. - 2019. - № 5. - С. 809-817.
27. Волоконно-оптические датчики для мониторинга коррозионных процессов в узлах авиационной техники (обзор) / Е. Н. Каблов, О. В. Старцев, И. М. Медведев, И. С. Шелемба // Авиационные материалы и технологии. - 2017. - № 3. - С. 26-34.
28. Рузаков И.А. Мониторинг деформационного состояния элементов конструкции из ПКМ на основе волоконно-оптических датчиков (обзор) // Труды ВИАМ. - 2019. - № 4. -С. 88-98.
29. Применение оптоволоконных технологий при создании встроенных систем самодиагностики авиационных конструкций / А.Н. Серьезнов, А.Б. Кузнецов, А.В. Лукьянов, А.А. Брагин // Научный вестник НГТУ. Авиационная и ракетно-космическая техника. - 2016. - № 3. - С. 95-105.
30. Будадин О.Н., Кульков А.А., Кутюрин В.Ю. Волоконно-оптические датчики с решетками Брэгга для мониторинга напряженно-деформированного состояния изделий из композиционных материалов // Конструкции из композиционных материалов. - 2018. - № 2. - С. 60-67.
31. Разработка чувствительного элемента волоконно-оптического тензометрического датчика на основе решеток Брэгга / А.С. Мунько, С.В. Варжель, С.В. Архипов, К.А. Кон-нов, А.Б. Петров // Известия высших учебных заведений. Приборостроение. - 2017. - № 4. - С. 340-346.
32. Экспериментальные исследования по определению деформаций образцов из полимерного композиционного материала с применением волоконно-оптических датчиков / Н.А. Кошелева, Г.С. Шипунов, А.А. Воронков, Н.П. Мерку-шева, А. А. Тихонова // Вестник Пермского национального исследовательского политехнического университета. Аэрокосмическая техника. - 2017. - № 50. - С. 26-35.
33. Созонов Н.С., Шардаков И.Н. Экспериментально-теоретическое исследование влияния конструктивных элементов волоконно-оптического датчика деформации на его показания и методика корректировки передаточной функции // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2017. - № 4. - С. 173-188.
34. Численное моделирование напряженно-деформированного состояния композитного шпангоута авиационного назначения для разработки методики контроля с применением волоконно-оптических датчиков / А.Н. Аношкин, В.Ю. Зуйко, К.А. Пеленев, П.В. Писарев, Г.С. Шипунов // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2018. - № 4. - С. 47-57.
References
1. Vorobej V.V., Markin V.B. Kontrol' kachestva izgotovleniya i tekhnologiya remonta kompozitnyh konstrukcij [Quality Control and Technology of Repair of Composite Structures]. Barnaul, Izd-vo MC EOR, 2015, 181 p.
2. Bormotov A.N. Dekompoziciya sistem i ierarhicheskie struktury pokazatelej kachestva kompozitov [Decomposition of systems and hierarchical structures of quality indicators of compo-
35. Качура С.М., Постнов В.И. Перспективные оптоволоконные датчики и их применение (обзор) // Труды ВИАМ. -
2019. - № 5. - С. 52-61.
36. Изображающий спектрометр на основе дискретного интерференционного фильтра / Р.В. Скиданов, А.А. Морозов, А.П. Порфирьев, В.А. Бланк // Компьютерная оптика. -2015. - № 5. - С. 716-720.
37. Скиданов Р.В., Морозов А.А., Гоголева С.Ю. Решение уравнений Фредгольма первого рода с использованием метода расширенных регуляризованных нормальных уравнений в задаче восстановления гиперспектральных изображений // Информационные технологии и нанотехнологии: материалы междунар. конф. и молодежной школы, Самара, 17-19 мая 2016 г. - Самара: Изд-во Сам. гос. аэрокосм. ун-та им. акад. С.П. Королева. - С. 273-278.
38. Численные методы решения некорректных задач / А.Н. Тихонов, А.В. Гончарский, В.В. Степанов, А.Г. Ягола. -М.: Наука, 1990. - 232 с.
39. Лукьяненко Д.В., Ягола А.Г. Использование многопроцессорных систем для решения обратных задач, сводящихся к интегральным уравнениям Фредгольма 1-го рода // Труды Института математики и механики УрО РАН. - 2012. -№ 1. - С. 222-234.
40. Pan'kov A.A. Piezoelectroluminescent fiber-optic sensors for temperature and deformation fields // Sensors and Actuators A: Physical. - 2019. - Vol. 288. - P. 171-176.
41. Патент RU № 2723921. Способ измерения деформаций / Паньков А.А., опубл.: 18.06.2020 Бюл. № 17, заявка № 2019136251 от 11.11.2019 г.
42. Паньков А.А. Математическая модель диагностирования микропористости материалов оптоволоконным датчиком с распределенной брэгговской решеткой // Оптический журнал. - 2020. - № 4. - С. 3-10.
43. Паньков А.А. Криволинейный оптоволоконный датчик деформаций с распределенной брэгговской решеткой в структуре полимерного композита // Оптический журнал. -
2020. - № 8. - С. 3-11.
44. Кристенсен Р. Введение в механику композитов. -М.: Мир, 1982. - 334 с.
45. Паньков А.А. Статистическая механика пьезокомпо-зитов. - Пермь: Изд-во Перм. гос. техн. ун-та, 2009. - 480 с.
46. Волков С.Д., Ставров В.П. Статистическая механика композитных материалов. - Мн.: Изд-во Белорус. гос. ун-та, 1978. - 208 с.
47. Pan'kov A.A. Mathematical model for diagnosing strains by an optical fiber sensor with a distributed Bragg grating according to the solution of a Fredholm integral equation // Mechanics of Composite Materials. - 2018. - No. 4. - P. 513-522.
48. Лехницкий С.Г. Кручение анизотропных и неоднородных стержней. - М.: Наука, 1971. - 240 с.
49. Лехницкий С.Г. Теория упругости анизотропного тела. - М.: Наука, 1977. - 416 с.
sites]. Sovremennye naukoemkie tekhnologii, 2016, no. 9-2, pp. 196-203.
3. Kartashova E.D., Mujzemnek A.YU. Tekhnologicheskie defekty polimernyh sloistyh kompozicionnyh materialov [Technological defects of polymer layered composite materials]. Izvestiya vysshih uchebnyh zavedenij. Povolzhskij region. Tekhnicheskie nauki, 2017, no. 2, pp. 79-89.
4. Murashov V.V., Rumyancev A.F. Defekty monolitnyh detalej i mnogoslojnyh konstrukcij iz polimernyh kompozicionnyh materialov i metody ih vyyavleniya. CHast' 1. Defekty monolitnyh detalej i mnogoslojnyh konstrukcij iz polimernyh kompozicionnyh materialov [Defects of monolithic parts and multilayer structures made of polymer composite materials and methods of their detection. Part 1. Defects of monolithic parts and multilayer structures made of polymer composite materials]. Kontrol'. Diagnostika. 2007, no. 4, pp. 23-32.
5. Kolpachkov E.D., Petrova A.P., Kurnosov A.O., Soko-lov I.I. Metody formovaniya izdelij aviacionnogo naznacheniya iz PKM (obzor) [Methods of forming aircraft products from PKM (review)]. Trudy VIAM, 2019, no. 11, pp. 22-36.
6. Kupriyanova O.A., Serzhantova M.V., Banshchikova M.N. Parametry kontrolya kachestva proizvodstva kompozitnyh reflektorov antenn [Parameters of quality control of production of composite antenna reflectors]. ZHurnal Sibirskogo federal'nogo universiteta. Seriya: tekhnika i tekhnologii, 2015, no. 3, pp. 324-330.
7. Kovalenko A.V. Issledovanie svojstv svyazuyushchego dlya formovaniya izdelij metodom propitki pod davleniem [Study of the properties of the binder for forming articles by pressure impregnation]. Trudy VIAM, 2015, no. 1, pp. 30-37.
8. Kovalenko A.V., Sidel'nikov N.K., Sokolov I.I., Tundajkin K.O. Sferoplastik s reguliruemoj vyazkost'yu dlya zapolneniya uchastkov sotovyh konstrukcij [Spheroplastic with adjustable viscosity to fill sections of cellular structures]. Trudy VIAM, 2019, no. 11, pp. 37-43.
9. Gajdachuk O.V., Kondrat'ev A.V., Omel'chenko K.V. Analiz tekhnologicheskih defektov, voznikayushchih v serijnom proizvodstve integral'nyh aviakonstrukcij iz polimernyh kompozicionnyh materialov [Analysis of technological defects arising in serial production of integral aircraft structures from polymer composite materials]. Aviacionno-kosmicheskaya tekhnika i tekhnologiya, 2010, no. 3, pp. 11-20.
10. Gajdachuk A.V., Slivinskij M.V. Tekhnologicheskie defekty sotovyh konstrukcij iz polimernyh kompozicionnyh materialov [Technological defects of cellular structures made of polymer composite materials]. Aviacionno-kosmicheskaya tekhnika i tekhnologiya, 2005, no. 4, pp. 5-8.
11. Dikov I.A., Bojchuk A.S. Sposoby opredeleniya ob"emnoj doli por v polimernyh kompozicionnyh materialah s pomoshch'yu ul'trazvukovyh metodov nerazrushayushchego kontrolya (obzor) [Methods of determining the volume fraction of pores in polymer composite materials using ultrasonic methods of non-destructive testing (review)]. Trudy VIAM, 2017, no. 2, pp. 80-95.
12. Anoshkin A.N., Zujko V.YU., Osokin V.M., Pelenev K.A., Tret'yakov A.A. Issledovanie mekhanicheskogo povedeniya kompozitnoj konstrukcii pri razrabotke metodiki dostovernogo vyyavleniya defektov rentgenograficheskim metodom nerazru-shayushchego kontrolya [Study of mechanical behavior of composite structure during development of the method of reliable detection of defects by X-ray method of non-destructive control]. Matematicheskoe modelirovanie v estestvennyh naukah, 2017, no. 1, pp. 169-172.
13. Fedotov M.Yu. Koncepcii sozdaniya i tendencii razvitiya intellektual'nyh materialov (obzor) [Concepts of creation and trends of development of intellectual materials (review)]. Aviacionnye materialy i tekhnologii, 2015, no. 1, pp. 71-80.
14. Kablov E.N., Morozov G.A., Sivakov D.V., Gulyaev I.N. Vstroennyj kontrol': ot datchikov do informkompozitov [Built-in control: from sensors to information complexes]. Aviacionnye materialy 75 let. Izbrannye trudy, 2007, pp. 331-342.
15. Patent RU no. 2643692. Volokonno-opticheskij datchik ob"emnogo napryazhennogo sostoyaniya [Fiber-optic sensor of volumetric stressed state] / Pan'kov A.A., opubl.: 05.02.2018, Byul. no. 4; zayavka no. 2017111405 ot 04.04.2017 g.
16. Patent RU no. 2604174. Sistema gasheniya spektra akusticheskih shumov [Acoustic noise spectrum damping system] / Pan'kov A.A., Anoshkin A.N., Pisarev P.V., opubl. 10.12.2016 Byul. no. 34; zayavka RU no. 2015144968 ot 19.10.2015 g.
17. Patent RU no. 2697168. Lopast' vozdushnogo vinta s upravlyaemoj geometriej profilya [Propeller blade with controlled profile geometry] / Pan'kov A.A., Anoshkin A.N., Pisarev P.V., opubl.: 12.08.2019 Byul. no. 23, zayavka no. 2018140318 ot 14.11.2018 g.
18. Kolobkov A.S., Malahovskij S.S. Samozalechivayu-shchiesya kompozicionnye materialy (obzor) [Self-killing composite materials (review)]. Trudy VIAM, 2019, no. 1, pp. 47-54.
19. Sorokin K.V., Murashov V.V. Mirovye tendencii razvitiya raspredelennyh volokonno-opticheskih sensornyh sistem (obzor) [World trends in the development of distributed fibre-optic sensor systems (review)]. Aviacionnye materialy i tekhnologii, 2015, no. 3, pp. 90-94.
20. Kachura S.M., Postnov V.I. Perspektivnye optovolo-konnye datchiki i ih primenenie (obzor) [Promising fiber-optic sensors and their use (review)]. Trudy VIAM, 2019, no. 5, pp. 52-61.
21. Fedotov M.YU., Budadin O.N., Vasil'ev S.A., Med-vedkov O.I., Kozel'skaya S.O. Vozmozhnosti kontrolya vneshnih mekhanicheskih vozdejstvij volokonno-opticheskoj sistemoj diagnostiki, vstroennoj v ugleplastiki razlichnyh tipov [Capabilities of control of external mechanical effects by fibre-optic diagnostics system built into carbon plastics of various types]. Kontrol'. Diagnostika, 2019, no. 3, pp. 38-47.
22. Fedotov M.Yu. Razvitie tekhnologii opticheskogo kontrolya konstrukcij iz PKM volokonno-opticheskimi datchikami [Development of technology of optical control of structures from PKM with fibre-optic sensors]: avtoreferat dis. ... kandidata tekhnicheskih nauk; [Mesto zashchity: Nauch.-issled. in-t introskopii MNPO "Spektr"]. Moskva, 2019, 24 p.
23. Kasharina L.A., Mahsidov V.V., Smirnov O.I., Ruzakov I.A. Differencirovanie defektov v PKM po otkliku volokonno-opticheskih sensorov (obzor). CHast' I [Differentiation of defects in PKM by response of fibre-optic sensors (overview). Part I]. Trudy VIAM, 2019, no. 2, pp. 97-104.
24. Kasharina L.A., Mahsidov V.V. Differencirovanie defektov v PKM po otkliku volokonno-opticheskih sensorov (obzor). CHast' II [Differentiation of defects in PKM by response of fibre-optic sensors (overview). Part II]. Trudy VIAM, 2019, no. 6, pp. 43-50.
25. Ruzakov I. A. Monitoring deformacionnogo sostoyaniya elementov konstrukcii iz PKM na osnove volokonno-opticheskih datchikov (obzor) [Monitoring of deformation state of structural elements from PKM based on fibre-optic sensors (overview)]. Trudy VIAM, 2019, no. 4, pp. 88-98.
26. Vlasov A.A., Alejnik A.S., SHuklin F.A., Nikitenko A.N., Motorin E.A., Kireenkov A.YU. Obnaruzhenie ul'trazvukovyh vozdejstvij s primeneniem volokonnyh reshetok Bregga [Detection of ultrasonic effects using Bragg fiber grating]. Nauchno-tekhnicheskij vestnik informacionnyh tekhnologij, mekhaniki i optiki, 2019, no. 5, pp. 809-817.
27. Kablov E.N., Starcev O.V., Medvedev I.M., Shelemba I.S. Volokonno-opticheskie datchiki dlya monitoringa korrozionnyh processov v uzlah aviacionnoj tekhniki (obzor) [Fibre-optic sensors for monitoring corrosion processes in aviation equipment nodes (review)]. Aviacionnye materialy i tekhnologii, 2017, no. 3, pp. 26-34.
28. Ruzakov I.A. Monitoring deformacionnogo sostoyaniya elementov konstrukcii iz PKM na osnove volokonno-opticheskih datchikov (obzor) [Monitoring of deformation state of structural elements from PKM on the basis of fibre-optic sensors (overview)]. Trudy VIAM, 2019, no. 4, pp. 88-98.
29. Ser'yoznov A.N., Kuznecov A.B., Luk'yanov A.V., Bragin A.A. Primenenie optovolokonnyh tekhnologij pri sozdanii vstroennyh sistem samodiagnostiki aviacionnyh konstrukcij [Application of fiber-optic technologies in creation of built-in systems of self-diagnostics of aviation structures]. Nauchnyj vestnik NGTU. Aviacionnaya i raketno-kosmicheskaya tekhnika, 2016, no. 3, pp. 95-105.
30. Budadin O.N., Kul'kov A.A., Kutyurin V.Yu. Volokon-no-opticheskie datchiki s reshyotkami Bregga dlya monitoringa napryazhenno-deformirovannogo sostoyaniya izdelij iz kompo-zicionnyh materialov [Fibre-optic sensors with Bragg screens for monitoring of stressed-deformed state of products from composite materials]. Konstrukcii iz kompozicionnyh materialov, 2018, no. 2, pp. 60-67.
31. Mun'ko A.S., Varzhel' S.V., Arhipov S.V., Konnov K.A., Petrov A.B. Razrabotka chuvstvitel'nogo elementa volokonno-opticheskogo tenzometricheskogo datchika na osnove reshetok Bregga [Development of a sensitive element of a fibre-optic strain gauge based on Bragg gratings]. Izvestiya vysshih uchebnyh zavedenij. Priborostroenie, 2017. no. 4, pp. 340-346.
32. Kosheleva N.A., SHipunov G.S., Voronkov A.A., Merkusheva N.P., Tihonova A.A. Eksperimental'nye issledovaniya po opredeleniyu deformacij obrazcov iz polimernogo kompo-zicionnogo materiala s primeneniem volokonno-opticheskih datchikov [Experimental studies on determination of deformations of samples from polymer composite material using fibre-optic sensors]. Vestnik PNIPU. Aerokosmicheskaya tekhnika, 2017, no. 50, pp. 26-35.
33. Sozonov N.S., SHardakov I.N. Experimental-theoretical study of influence of structural elements of fiber-optic deformation sensor on its readings and method of correction of transfer function. PNRPU Mechanics Bulletin, 2017, no. 4, pp. 173-188.
34. Anoshkin A.N., Zujko V.YU., Pelenev K.A., Pisarev P.V., Shipunov G.S. Simulation of stress-strain state of composite frame of aviation purpose for development of control method using fibre-optic sensors. PNRPU Mechanics Bulletin, 2018, no. 4, pp. 47-57.
35. Kachura S.M., Postnov V.I. Perspektivnye optovolokon-nye datchiki i ih primenenie (obzor) [Promising fiber-optic sensors and their use (review)]. Trudy VIAM, 2019, no. 5, pp. 52-61.
36. Skidanov R.V., Morozov A.A., Porfir'ev A.P., Blank V.A. Izobrazhayushchij spektrometr na osnove diskretnogo interferen-cionnogo fil'tra [Depicting spectrometer based on discrete interference filter]. Komp'yuternaya optika, 2015, no. 5, pp. 716-720.
37. Skidanov R.V., Morozov A.A., Gogoleva S.Yu. Reshenie uravnenij Fredgol'ma pervogo roda s ispol'zovaniem metoda rasshirennyh regulyarizovannyh normal'nyh uravnenij v zadache vosstanovleniya giperspektral'nyh izobrazhenij [Solution of Fredgolm Equations of the First Kind Using the Method of Ex-
tended Regulated Normal Equations in the Problem of Restoring Hyperspectration Images]. Materialy mezhdunarodnoj konferencii i molodyozhnoj shkoly «Informacionnye tekhnologii i nano-tekhnologii», Samara, 17 19 maya 2016 g., Samara: Izd.-vo Sam. gos. aerokosm. un-ta im. akad. S.P. Korolyova, pp. 273-278.
38. Tihonov A.N., Goncharskij A.V., Stepanov V.V., Yagola A.G. CHislennye metody resheniya nekorrektnyh zadach [Numerical methods of solving incorrect problems]. Moscow, Nauka, 1990, 232 p.
39. Luk'yanenko D.V., YAgola A.G. Ispol'zovanie mnogoprocessornyh sistem dlya resheniya obratnyh zadach, svodyashchihsya k integral'nym uravneniyam Fredgol'ma 1-go roda [Use of multiprocessor systems to solve reverse problems, reduced to integral equations of Fredholm of the 1st kind]. Trudy Instituta matematiki i mekhaniki UrO RAN, 2012, no. 1, pp. 222-234.
40. Pan'kov A.A. Piezoelectroluminescent fiber-optic sensors for temperature and deformation fields. Sensors and Actuators A: Physical, 2019, vol. 288, pp. 171-176.
41. Patent RU no. 2723921. Sposob izmereniya deformacij [Deformation Measurement Method] / Pan'kov A.A., opubl.: 18.06.2020 Byul. no. 17, zayavka no. 2019136251 ot 11.11.2019 g.
42. Pan'kov A.A. Matematicheskaya model' diagnostiro-vaniya mikroporistosti materialov optovolokonnym datchikom s raspredelyonnoj breggovskoj reshyotkoj [Mathematical model of diagnostics of microporosity of materials with optical fiber sensor with distributed Bragg grating]. Opticheskij zhurnal, 2020, no. 4, pp. 3-10.
43. Pan'kov A.A. Krivolinejnyj optovolokonnyj datchik deformacij s raspredelyonnoj breggovskoj reshyotkoj v strukture polimernogo kompozita [Curved fiber-optic deformation sensor with distributed Bragg grating in the structure of polymer composite]. Opticheskij zhurnal, 2020, no. 8, pp. 3-11.
44. Kristensen R. Vvedenie v mekhaniku kompozitov [Introduction to composites mechanics]. Moscow, Mir, 1982, 334 p.
45. Pan'kov A.A. Statisticheskaya mekhanika p'ezokompo-zitov [Statistical mechanics of piezocomposites]. Perm', Izd-vo Perm. gos. tekhn. un-ta, 2009, 480 p.
46. Volkov S.D., Stavrov V.P. Statisticheskaya mekhanika kompozitnyh materialov [Statistical mechanics of composite materials]. Minsk, Izd-vo Belorus. gos. un-ta, 1978, 208 p.
47. Pan'kov A.A. Mathematical model for diagnosing strains by an optical fiber sensor with a distributed Bragg grating according to the solution of a Fredholm integral equation [Mathematical model for diagnosing strains by an optical fiber sensor with a distributed Bragg grating according to the solution of a Fredholm integral equation]. Mechanics of Composite Materials, 2018, no. 4, pp. 513-522.
48. Lehnickii S.G. Kruchenie anizotropnih i neodnorodnih sterjnei [Twisting of anisotropic and inhomogeneous rods]. Moscow, Nauka, 1971, 240 p.
49. Lehnickii S.G. Teoriya uprugosti anizotropnogo tela [The theory of elasticity of an anisotropic body]. Moscow, Nauka, 1977, 416 p.