Mathematical model of pulse scanning of pressure along a piezo-electro-luminescent fiber-optical sensor
- Authors: Pan’kov A.A1
- Affiliations:
- Perm National Research Polytechnic University
- Issue: No 1 (2018)
- Pages: 73-82
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/mechanics/article/view/108
- DOI: https://doi.org/10.15593/perm.mech/2018.1.06
- Cite item
Abstract
The mathematical model is developed aimed at locating pressure non-uniformity along a fiber optic-piezo-electro-luminescent sensor using the location scanning electrical video pulse with a step by step change of its value. The algorithm of finding the function of pressure distribution along a local section and the whole length of the sensor is developed based on the results of the intensity of light proceeding from a fiber optic phase measured on the edge section of the sensor for a case of nonlinear “luminescence function”, which is a dependence between the intensity of light and voltage acting on the luminescent element. The problem is reduced to the solution of the Fredholm integral equation of the 1st kind with the differential kernel depending on the operating and informative transfer coefficients of the sensor and on the set luminescence function of the element. Analytical solutions for the pressure distribution functions along the sensor are obtained for special cases, when the kernel or the function of distribution density is expressed via delta function, and the Fredholm integral equation becomes algebraic. Domains of admissible values of the operating voltage for various modes of diagnostics of pressure distribution are defined. Results of numerical solutions of direct and reverse problems for non-uniform distribution of pressure by means of “pointed” scanning of this pressure are presented by an extremely narrow impulse of the operating voltage. Luminescence functions at the exit from the optical fiber at various time points and values of the impulse of the operating voltage taking into account the set luminescence function of an electroluminescent element are found in the direct problem. The distribution of pressure depending on the values of the luminescence intensity function at various time points is found in the reverse problem.
Full Text
Введение Сенсорика за последние годы сформировалась в самостоятельную высокотехнологичную ветвь измерительной техники [1], одной из центральных задач которой является разработка все более совершенных датчиков для диагностирования различных физико-механических параметров исследуемых объектов. К новым датчикам предъявляются все более высокие требования. В частности, наряду с высокими метрологическими характеристиками особое значение придается их высокой надежности, долговечности, стабильности, малым габаритам, массе и энергопотреблению, возможности встраивания в объекты и системы мониторинга, совместимости с микроэлектронными и компьютерными устройствами обработки информации, низкой себестоимости, что в наибольшей степени характеризует именно оптико-электронные волоконные датчики. В [1] рассмотрены классификация, принципы действия, конструкции, параметры и характеристики оптико-электронных и волоконно-оптических датчиков, дан анализ влияния различных факторов на характеристики датчиков и приведены принципиальные схемы построения электронных и оптических датчиков для решения различных промышленных, технологических и медицинских задач. Различие между активными и пассивными датчиками обусловлено различием принципов их функционирования, отражающих, в свою очередь, фундаментальные отличия: активные датчики сами являются генераторами информативного выходного сигнала различной физической природы, пассивные - имеют лишь информативное воздействие на процесс, в частности, через изменение удельного сопротивления изменяется сила электрического тока в проводнике или через изменение показателя преломления оптоволокна (световода) изменяются выходные характеристики светового потока. В активном пьезоэлектрическом датчике параметр измерения, в частности внешнее давление, воздействует на пьезоэлектрический элемент, что приводит к деформации и появлению на его гранях электрических зарядов противоположного знака, величина которых связана с измеряемым давлением. Световод может использоваться в качестве линии передачи информативного сигнала и/или чувствительного элемента, в частности, в виде участка световода с дифракционной решеткой Брэгга [2-6]; механолюминесцентный эффект использован в конструкциях датчиков [7]. В [8-11] представлены результаты теоретических и экспериментальных исследований по оценке механических свойств и деформаций полимерных композиционных материалов с помощью интегрированных в структуры материалов оптических волокон с решетками Брэгга. В [12] найдено распределение деформации по длине оптоволоконного датчика, и установлена минимально допустимая длина датчика для различных случаев нагружений в результате численного моделирования деформационных процессов в системе «подложка-клей-оптоволоконный датчик» в рамках линейной теории термовязкоупругости. Возможность измерения градиентных полей деформаций волоконно-оптическими датчиками с решетками Брэгга продемонстрирована в [13] на примере пластины с вырезами в сравнении с результатами, полученными при использовании трехмерной цифровой оптической системы, и с результатами численного моделирования методом конечных элементов. Оценки технологических деформаций в полимерных композиционных образцах с использованием внедренных волоконно-оптических датчиков с решетками Брэгга даны в [14]. Другим перспективным решением проблемы диагностирования неоднородных полей деформаций и наличия дефектов в элементах композитных конструкций являются сенсорные пьезоэлектрические сети [15-17], внедренные в структуру конструкции. Моделирование процессов в пьезоэлектрических чувствительных элементах датчиков проводится на основе решений соответствующих связанных краевых задач электроупругости аналитическими [18-21] или численными [22-26] методами механики сплошной среды. В [27-30] предложены новые конструкции пьезоэлектролюминесцентных оптоволоконных датчиков для уточненного диагностирования давления [27-29] и сложного объемного напряженно-деформированного состояния композитных конструкций [30] с использованием алгоритмов [28, 29] обработки приемником-анализатором интегральных оптических сигналов на выходе из световода датчика. Информативный световой сигнал возникает на локальном участке датчика в силу механолюминесцентного эффекта, обусловленного взаимодействием пьезоэлектрического и электролюминесцентного слоев датчика, и передается по световоду к приемнику-анализатору; наличие управляющих электродов позволяет диагностировать локации неоднородностей давления по длине датчика. Цель - разработка математической модели импульсной локации неоднородностей давления по длине оптоволоконного пьезоэлектролюминесцентного датчика [27] с использованием локационного сканирующего электрического видеоимпульса на управляющих электродах по результатам замеряемой на торцевом сечении волоконного датчика интенсивности исходящего из оптоволокна света для различных моментов времени; для частного случая - сканирования предельно узким импульсом - получить аналитическое решение для диагностируемого неоднородного давления, что значительно упростит известный алгоритм [29] обработки светового сигнала на выходе из оптоволокна датчика. 1. Оптоволоконный пьезоэлектролюминесцентный датчик Датчик давления (рис. 1) [27-29] представляет собой составное слоистое волокно, состоящее из оптоволокна 1, электролюминесцентного 2 и с радиальной поляризацией пьезоэлектрического 3 цилиндрических концентрических слоев с внутренним фотопрозрачным, в частности, перфорированным или сеточным 4 и внешним 5 непрерывными управляющими электродами. Действие радиального напряжения (давления) на некотором локальном участке внешней боковой поверхности датчика приводит к соответствующим деформациям участка пьезоэлектрического слоя 3 и появлению в нем и в смежном с ним локальном участке электролюминесцентного слоя 2 электрического поля, приводящего к свечению локального участка электролюминесцентного слоя внутрь (через фотопрозрачный электрод 4) оптоволокна 1; управляющие электроды 4, 5 с напряжением дают возможность изменять, в частности, однородно или неоднородно по длине датчика электрическое напряжение на электролюминесцентном слое 2 и как результат величину интенсивности света на выходе из оптоволокна. 5 1 3 2 4 Рис. 1. Датчик давления с управляющим напряжением и интенсивностью света на выходе из световода Fig. 1. Sensor of pressure with the operating voltage and intensity of light at the exit from the light guide 2. Математическая модель сканирования Рассмотрим процесс сканирования неоднородного давления по оси датчика локационным электрическим управляющим напряжением на электродах, которое представим в виде суммы некоторого постоянного по длине значения и локационного электрического видеоимпульса прямоугольной формы , распространяющегося в положительном направлении оси . (1) В произвольный момент времени импульсная функция отлична от нуля и равна значению лишь на локальном участке с координатами и длиной с центром в точке , где координаты начального и конечного положений центра импульса и соответственно, координата длина электрода время прохождения центром одиночного импульса от начального до конечного положений скорость распространения импульса . Таким образом, для диагностики давления на участке с центром имеем (после проведенных в циклах диагностики замеров интенсивностей света на выходе оптоволокна) функцию , заданную через узловые значения интенсивности света на выходе оптоволокна и управляющего напряжения на участке , где моменты времени прохождения импульсом этого участка Однородная составляющая и величина импульса постоянны в течение каждого цикла диагностики и варьируются, в частности увеличиваются на некоторую малую величину лишь при начале следующего цикла для случая , где минимальное и максимальное значения амплитуда импульса - число циклов диагностики. Рассмотрим вероятностную модель, в которой распределение действующего на боковую цилиндрическую поверхность волоконного датчика внешнего давления по координате описываем непрерывной случайной величиной , реализациями которой являются различные значения действующего давления , где ; аналогично для случайной величины реализациями являются различные значения на соответствующем локальном участке датчика. Электрическое напряжение на электролюминесцентном локальном участке датчика (2) также является случайной величиной, связанной с давлением параметрами датчика , . Найдем функции плотности распределения этого давления по анализу интенсивности света на выходе из оптоволокна датчика. Считаем, что свойства электролюминесцентного элемента заданы «S-образной» кривой зависимости интенсивности свечения от приложенного напряжения с характерными точками заданных пороговых напряжений: для начала свечения и для начала насыщенного свечения электролюминесцентного элемента (рис. 2). Алгоритм расчета функции плотности распределения давления на локальном участке оптоволоконного датчика рассмотрен в [29], согласно которому искомая функция является решением интегрального уравнения Фредгольма 1-го рода , (3) Рис. 2. Заданная функция и производная интенсивности свечения электролюминесцентного элемента от значений приложенного напряжения Fig. 2. The set function and derivative of luminescence intensity of the electroluminescent element on the values of the enclosed voltage где границы интегрирования (4) с учетом (2). В каноническом виде интегральное уравнение Фредгольма (3) имеет вид , (5) где разностное ядро (6) с учетом новых переменных , , (7) известная с учетом (7) и искомая функции, границы интегрирования: , или , с учетом (4), , . Различные аналитические и численные методы решения интегральных уравнений Фредгольма (5) даны в [31-35]. В частном случае, когда зависимость - ступенчатая функция (пунктирная линия на рис. 2), производная в ядре (6) уравнения (5) выражается через дельта-функцию Дирака и из (3) получим решение для искомой функции плотности распределения , (8) где константа максимальная интенсивность свечения участка электролюминесцентного элемента оптоволоконного датчика длиной . Для случая, когда на рассматриваемом участке диагностируемое давление или когда ширина сканирующего импульса является предельно малой величиной и давление для функция плотности распределения вырождается к виду (9) и из интегрального (3), (5) получим алгебраическое уравнение , (10) функция интенсивности свечения на выходе из оптоволокна , (11) где обусловлено возможным свечением участка . В частности, значение в (1) может быть подобрано из условия отсутствия свечения на выходе из световода в отсутствие импульса. В результате получим искомое решение для величины диагностируемого давления (12) где - функция, обратная известной (см. рис. 2) величина (13) выражается через экспериментально измеренную интенсивность для соответствующего значения управляющего напряжения . Отметим, что в (12) - это давление в малой окрестности внешней границы сечения датчика с координатой ; величина рассчитывается через интенсивность света на выходе из световода в момент времени где - скорость распространения электрического импульса (света), начальный момент времени соответствует моменту входа импульса в датчик. 3. Области допустимых значений параметров управляющего напряжения Для диагностики оптоволоконным датчиком значений давления из всего рабочего диапазона необходимо выполнение условий: во-первых, минимальное значение управляющего сигнала должно приводить к прекращению свечения ( ) даже на участке, где действует максимальное давление , и, во-вторых, максимальное значение управляющего сигнала должно приводить к максимальному «насыщенному» свечению ( ) даже на участке, где действует минимальное давление . Эти условия с учетом того, что напряжение на электролюминесцентном элементе, например на участке действия импульса, имеет вид (2), запишем как систему неравенств (14) с использованием заданных пороговых напряжений для начала свечения и для начала насыщенного свечения электролюминесцентного элемента. Для случая, когда величина импульса , систему (14) преобразуем к виду (15) с учетом равенств (1), амплитуды импульса и допуская для общего случая вариации однородной составляющей управляющего напряжения на интервале значений . Таким образом, область допустимых значений параметров управляющего сигнала найдем из решения системы (16) где коэффициенты (17) с учетом того, что . Для случая, когда в (15), (16) полагаем равенства (18) В частности, при амплитуда импульса выбирается из области допустимых значений (16) (19) Для случая, когда импульс отсутствует, область допустимых значений границ для варьирования величины однородной составляющей управляющего напряжения находим из системы (20) 4. Результаты численного моделирования Численный расчет области допустимых значений параметров: управляющего напряжения на рис. 3 проведен (18) для значений границ: действующего на участке давления , пороговых напряжений для начала и насыщенного свечений электролюминесцентного элемента (см. рис. 2), параметров датчика , [28, 29] и рассчитанных значений коэффициентов (17). Вертикальным красным отрезком на рис. 3 обозначены допустимые значения амплитуды импульса для случая (19) отсутствия однородной составляющей в управляющем напряжении (1). На рис. 4, а построены зависимости интенсивности света (с учетом , где интенсивность свечения участка , который светится лишь при значении и производной на выходе из оптоволокна от управляющего напряжения при заданном равномерном законе распределения действующего давления на участке в границах , с параметрами датчика относительной длиной участка Значение амплитуды импульса было выбрано с учетом области допустимых значений параметров управляющего напряжения (см. рис. 3). Пунктирные графики (см. рис. 4, а) зависимостей интенсивности света и производной от управляющего напряжения на участке построены для случая, когда функция задана в виде ступенчатой функции (см. рис. 2) с параметрами Рис. 3. Область допустимых значений параметров , управляющего напряжения Fig. 3. Domain of admissible values of parameters: , of operating voltage На рис. 4, б для случая наличия лишь однородной по длине датчика составляющей варьируемого управляющего напряжения в отсутствие импульса построены зависимости интенсивности света и производной (5) на выходе из оптоволокна от управляющего напряжения при заданном равномерном законе распределения действующего давления по всей длине в границах Допустимые значения границ варьирования величины управляющего напряжения найдены из решения (20). Пунктирные графики на рис. 4, б рассчитаны для случая, когда функция задана в виде ступенчатой функции (см. рис. 2). Граничные значения давлений связаны с соответствующими характерными значениями управляющего напряжения зависимостью с учетом (2); для рассматриваемого случая (см. рис. 4, б) имеем значения Рис. 4. Рассчитанная функция свечения из оптоволокна и производная для диагностирования давления на локальном участке (а) и по всей длине (б) датчика Fig. 4. The calculated function of luminescence from optical fiber and derivative for diagnosing of pressure on the local segment (а) and along the whole length (b) of the sensor Рис. 5. Распределение давления (а) и функция свечения из оптоволокна для различных величин электрического импульса сканирования (б) Fig. 5. Distribution of pressure (а) and luminescence function from the optical fiber for various values of an electric impulse of scanning (b) Для неоднородного распределения давления (рис. 5, а) были решены прямая (11) и обратная (12), (13) задачи для заданной на рис. 2 функции свечения, где время прохождения импульсом всей длины датчика со скоростью ширина импульса начальный момент времени соответствует моменту входа импульса в датчик, момент выхода света из световода. В прямой задаче найдены (рис. 5, б) функции свечения на выходе из оптоволокна для различных численных значений: (○), (Δ), (□), (◊), ( ), ( ) электрического импульса сканирования; предел насыщенного свечения наступает при значении (пунктирная линия). В обратной задаче найдено распределение давления (см. рис. 5, а) по значениям функции интенсивности свечения , в частности для значения управляющего напряжения (◊), для которого отсутствуют «неинформативные зоны» из участков с «нулевыми» ( ) и «насыщенными» значениями интенсивности свечения (см. рис. 5, б), где относительная ширина импульса . Отметим, что если для некоторого значения есть неинформативные зоны для функции интенсивности свечения , то необходимо дополнительно провести измерения для другого значения уже с информативными участками в этих зонах, в частности взаимодополняющие графики (□), ( ) на рис. 5, б. Таким образом, сканирование неоднородного давления предельно узким видеоимпульсом позволяет получить аналитическое решение (12), что значительно упрощает известный алгоритм [29] обработки светового сигнала на выходе из оптоволокна. Заключение Разработаны алгоритм и математическая модель локации неоднородностей давления по длине оптоволоконного пьезоэлектролюминесцентного датчика с использованием локационного сканирующего электрического видеоимпульса с пошаговым изменением величины импульса на каждом цикле прохождения исследуемого локального участка . Показано, что расчет функции плотности распределения давления на локальном участке внешней боковой цилиндрической поверхности оптоволоконного датчика по результатам замеряемой на торцевом сечении датчика интенсивности света из оптоволоконной фазы для случая нелинейной зависимости интенсивности света от действующего на электролюминесцентный элемент напряжения (см. рис. 2) сводится к решению интегрального уравнения Фредгольма 1-го рода (5) с разностным ядром (6), зависящим от передаточных коэффициентов , датчика, относительной длины исследуемого локального участка и производной (см. рис. 2) заданной функции свечения электролюминесцентного элемента. Определены области допустимых значений управляющего напряжения датчика для различных режимов диагностики распределения давления (см. рис. 3). Для случая, когда ширина ( ) сканирующего импульса является предельно малой величиной, функция плотности распределения выражается через обобщенную дельта-функцию Дирака (9) и из интегрального уравнения Фредгольма (3), (5) получено аналитическое решение для величины диагностируемого давления (12) по результатам измерений интенсивности света на выходе из оптоволокна в различные моменты времени и управляющего напряжения . Представлены результаты численных решений прямой (11) и обратной (12), (13) задач для неоднородного распределения давления (см. рис. 5, а) посредством «точечного» сканирования этого давления предельно узким импульсом управляющего электрического напряжения В прямой задаче найдены (см. рис. 5, б) функции свечения на выходе из оптоволокна для различных моментов времени и значений с учетом заданной на рис. 2 функции свечения электролюминесцентного элемента; в обратной задаче найдено распределение давления (см. рис. 5, а) по значениям функции интенсивности свечения .About the authors
Andrey A Pan’kov
Perm National Research Polytechnic University
References
- Козлов В.Л. Оптоэлектронные датчики. - Минск: Изд-во Белорус. гос. ун-та, 2005. - 116 с.
- Окоси Т. Волоконно-оптические датчики. - Л.: Энергоатомиздат, 1990. - 256 с.
- Guemes A., Fernandez-Lopez A., Soller B. Optical fiber distributed sensing - physical principles and applications // Structural Health Monitoring. - 2010. - Vol. 9. - No. 3. - P. 233-245.
- Suresh R., Tjin S.C., Hao J. Fiber Bragg Grating // Smart Materials in Structural Health Monitoring, Control and Biomechanics. - Berlin: Springer, Heidelberg, 2012. - P. 413-439.
- Prabhugoud M., Peters K. Efficient simulation of Bragg grating sensors for implementation to damage identification in composites // Smart Materials & Structures. - 2003. - Vol. 12. - No. 6. - P. 914-924.
- Методы исследования свойств материалов при интенсивных динамических нагрузках / под общ. ред. д-ра физ.-мат. наук М.В. Жерноклетова; ФГУП РФЯЦ-ВНИИЭФ. - Саров, 2003. - 403 с.
- Татмышевский К.В. Механолюминесцентный чувствительный элемент: математическая модель и динамические свойства // Приборы и системы. Управление, контроль, диагностика. - 2005. - № 4. - С. 35-39.
- Матвеенко В.П., Федорова В.А., Шардаков И.Н. Теоретическое обоснование возможности построения волоконно-оптической системы мониторинга деформаций земной поверхности // Изв. РАН. МТТ. - 2013. - № 5. - С. 46-52.
- Шардаков И.Н., Созонов Н.С., Цветков Р.В. Экспериментально-теоретические основы автоматизированных систем деформационного мониторинга с использованием волоконно-оптических элементов // Вестн. Перм. науч. центра. - 2016. - Окт.-дек. - С. 91-95.
- К вопросу об интеграции оптоволокна в ПКМ и измерении деформации материала с помощью волоконных брэгговских решеток / М.А. Зуев, В.В. Махсидов, М.Ю. Федотов, А.М. Шиенок // Механика композиционных материалов и конструкций. - 2014. - Т. 20, № 4. - С. 568-574.
- Применение оптического волокна в качестве датчиков деформации в полимерных композиционных материалах / Е.Н. Каблов, Д.В. Сиваков, И.Н. Гуляев, К.В. Сорокин, М.Ю. Федотов, Е.М. Дианов, С.А. Васильев, О.И. Медведков // Все материалы. Энцикл. справ. - 2010. - № 3. - С. 10-15.
- Наймушин И.Г., Труфанов Н.А., Шардаков И.Н. Численный анализ деформационных процессов в оптоволоконном датчике // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2012. - № 1. - С. 104-116.
- Измерение неоднородных полей деформаций встроенными в полимерный композиционный материал волоконно-оптическими датчиками / А.Н. Аношкин, А.А. Воронков, Н.А. Кошелева, В.П. Матвеенко, Г.С. Сероваев, Е.М. Спаскова, И.Н. Шардаков, Г.С. Шипунов // Изв. РАН. МТТ. - 2016. - № 5. - С. 42-51.
- Матвеенко В.П., Шардаков И.Н., Кошелева Н.А. Оценка технологических деформаций в полимерных композиционных образцах на основе использования внедренных волоконно-оптических датчиков деформаций // XХ Зимняя школа по механике сплошных сред: тез. докл. / РИО УрО РАН. - Пермь, 13-16 февраля 2017 г. - Екатеринбург: 2017. - С. 219.
- Patent US 20060254366 A1. Sensor and sensor array for monitoring a structure / Caroline Williamson, Lisa Fixter (nee Humberstone), Andrew Clarke; published on November 16, 2006.
- Patent US 7458266 B2. Method and apparatus for detecting a load change upon a structure and analyzing characteristics of resulting damage / Shawn J. Beard, Xinlin Qing, Hian Leng Chan, Chang Zhang, Fuo-Kuo Chang; published on December 2, 2008.
- Patent US 6399939 B1. Sensor array system / Mannur J. Sundaresan, Anindya Ghoshal, Mark J. Schulz; published on June 4, 2002.
- Гринченко В.Т., Улитко А.Ф., Шульга Н.А. Электроупругость. - Киев: Наукова думка, 1989. - 280 с.
- Партон В.З., Кудрявцев Б.А. Электромагнитоупругость пьезоэлектрических и электропроводных тел. - М.: Наука, 1988. - 472 с.
- Дианов Д.Б., Кузьменко А.Г. Расчет цилиндрического пьезокерамического преобразователя, совершающего радиально-симметричные колебания // Акуст. журн. - 1970. - Т. 16, № 1. - С. 42-48.
- Шляхин Д.А. Нестационарная осесимметричная задача электроупругости для анизотропного пьезокерамического радиально поляризованного цилиндра // Изв. РАН. МТТ. - 2009. - № 1. - С. 73-81.
- Finite Element Analysis of Functionally Graded Piezoelectric Spheres / A. Ghorbanpour Arani, R. Kolahchi, A.A. Mosalaei Barzoki, A. Loghman, F. Ebrahimi // Finite Element Analysis - Applications in Mechanical Engineering. Ed. Farzad Ebrahimi. - InTech, 2012. - 380 p.
- Numerical determination of effective properties of voided piezoelectric materials using BNM / H. Wang, G. Tan, S. Cen, Z. Yao // Engineering Analysis with Boundary Elements. - 2005. - Vol. 29. - P. 636-646.
- Li Z.H., Wang C., Chen C.Y. Effective electromechanical properties of transversely isotropic piezoelectric ceramics with microvoids // Comput. Mater. Sci. - 2003. - Vol. 27. - No. 3. - P. 381-392.
- Bishay P.L., Dong L., Atluri S.N. Multi-physics computational grains (MPCGs) for direct numerical simulation (DNS) of piezoelectric composite/porous materials and structures // Computational Mechanics. - 2014. - Vol. 54. - No. 5. - P. 1129-1139.
- Kondaiah P., Shankar K., Ganesan N. Pyroeffects on multiphase magneto-electro-elastic sensor patch bonded on mild steel plate // International Journal on Smart Sensing and Intelligent Systems. - 2014 September. - Vol. 7. - No. 3. - P. 1134-1155.
- Пат. РФ № 2630537. Волоконно-оптический датчик давления / Паньков А.А., опубл. 11.09.2017 г.; заявка RU № 2016136058 от 06.09.2016 г.
- Паньков А.А. Математическое моделирование пьезоэлектролюминесцентного эффекта и диагностика распределения давления по длине оптоволоконного датчика // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2016. - № 4. - С. 259-272.
- Pan’kov A.A. Piezoelectroluminescent optical fiber sensor for diagnostics of the stress state and defectoscopy of composites // Mechanics of Composite Materials. - 2017. - Vol. 53. - No. 2. - P. 229-242.
- Пат. РФ № 2643692. Волоконно-оптический датчик объемного напряженного состояния / Паньков А.А., опубл. 05.02.2018 г.; заявка № 2017111405 от 04.04.2017 г.
- Михлин С.Г., Смолицкий X.Л. Приближенные методы решения дифференциальных и интегральных уравнений. - М.: Наука, 1965. - 384 с.
- Верлань А.Ф. Сизиков В.С. Интегральные уравнения: методы, алгоритмы, программы. - Киев: Наук. думка, 1986. - 543 с.
- Васильева А.Б., Тихонов А.Н. Интегральные уравнения. - М.: Физматлит, 2002. - 360 с.
- Полянин А. Д., Манжиров А. В. Справочник по интегральным уравнениям. - М.: Физматлит, 2003. - 608 с.
- Латыпов А.Ф. Численные методы решения линейных интегральных уравнений Фредгольма и Вольтерра 1-го рода // Обратные и некорректные задачи математической физики: материалы междунар. конф., посвященной 75-летию академика М.М.Лаврентьева, 20-25 августа 2007 г. - Новосибирск, 2007.