# Abstract

The algorithm of finding the distribution function of pressure along the three-phase fiber optic sensor based on the results of light intensity proceeding from a fiber optic phase measured on the edge section of the sensor is developed for a case of nonlinear "function of a luminescence" which is a dependence between the intensity of light and voltage acting on the electrophosphor. The problem is reduced to the solution of the Fredholm integral equation of the 1st kind with the differential kernel depending on the calculated effective parameters of the sensor and on the derivative set function of a luminescence of an electrophosphor. The analytical solution has been obtained for the function of probabilities density of pressure distribution for a special case when the kernel is expressed by the delta-function; and the Fredholm integral equation is reduced to an algebraic one. "Direct" and "reverse" problems of the Fredholm integral equation for a case of nonlinear function of a luminescence of an electrophosphor are solved as an illustration of the algorithm. The light intensity derivative at the optical fiber output of control voltage for the set uniform law of probability-density function with regard to pressure distribution is found in a direct problem. As for the reverse problem, the probability-density function is determined in comparison with the known exact solution using the direct problem solution for the derivative of light intensity. The numerical solution of the reverse problem is carried out in different approximations in which the distribution of nodal points in intervals and required nodal function values of probability-density pressure are found from a condition of minimizing summary discrepancies based on the values of the light derivative intensity which have been set and calculated on each step based on control voltage at the output of the optical fiber.

# Full Text

Введение Существуют две группы волоконно-оптических датчиков: первая - с внешним чувствительным элементом, когда измерение параметров производится в области за пределами волокна и оптические волокна используются лишь для переноса света к датчикам и в обратном направлении, вторая - чистоволоконные, в основном интерферометрические датчики, в которых оптическое волокно используется в качестве чувствительного элемента. Оптические модуляторы являются ключевыми составляющими элементами волоконно-оптических систем, выполняющими различные функции, в том числе модуляцию амплитуды, фазы, частоты и поляризации. В большинстве случаев свет модулируется путем изменения оптических свойств материала модулятора при воздействии управляющего сигнала. Механизм связи управляющего сигнала со свойствами материала может быть электрооптическим, акустооптическим или магнитооптическим. В интегрально-оптических модуляторах волноводы встраиваются непосредственно в материал модулятора. Это существенно снижает требования к электрической мощности и исключает необходимость подсоединения к оптическим волокнам внешних чувствительных элементов. Один из наиболее популярных методов модуляции интенсивности света основан на изгибании волокна, вызывающем потери излучения. Такие датчики микроизгибов используются в тех приложениях, где измеряемые параметры: деформацию, давление, силу, положение, ускорение можно механически преобразовать в перемещение устройства, которое деформирует волокно, при этом потери излучения возрастают и количество прошедшего света уменьшается. Оптоволоконные датчики могут эффективно использоваться для диагностики напряженного состояния и дефектоскопии на структурном уровне внутри композитных элементов конструкций, в частности полимерных волокнистых изделий авиационного назначения. Изготовление композитных конструкций происходит в несколько этапов, на каждом из которых возможно образование дефектов, характерных для данной технологической стадии, поэтому необходимо проводить контроль качества на всех этапах создания и эксплуатации композитной конструкции [1]. По своей природе композиционный материал является сложным объектом для диагностики и контроля механических свойств на структурном уровне, он обладает анизотропией, существенной зависимостью свойств от температуры и технологии изготовления, но в самом принципе создания композиционного материала и конструкции одновременно из ряда отдельных компонентов, в частности: волокон, ткани, связующего, фактически создания материала-детали, можно найти решение проблемы обеспечения его надежности. Это создание материала с возможностью диагностики и в будущем с возможностью управления своими характеристиками. Для этого в структуру материала необходимо добавить информационные или интеллектуальные компоненты, обеспечивающие регистрацию и передачу информации о состоянии материала [1-5], а при необходимости и физическое воздействие с целью изменения его характеристик. Для дефектоскопии полимерных композитных конструкций в структуру композита на этапе изготовления может быть [2] введено оптоволокно, и о наличии дефектов в композите судят опосредованно через появление дефектов или разрыва оптоволокна, которые приводят к снижению яркости или отсутствию исходящего из оптоволокна светового потока при заданном входящем световом потоке. В [3] предложен волоконно-оптический датчик давления на основе использования «туннельного эффекта», в [4] рассматривается возможность создания системы распределенных волоконно-оптических датчиков давления на основе учета влияния микроизгиба оптоволокна на эффект обратного рассеяния светового потока. Наиболее близким к рассматриваемому в работе оптоволоконному датчику давления является механолюминесцентный датчик давления [5], который преобразует механическую энергию воздействия в оптическое излучение с использованием явления механолюминесценции: триболюминесценции или деформационной люминесценции [6-9]. Явление механолюминесценции наблюдается у различных кристаллофосфоров, при пластической деформации которых происходит движение электрически заряженных дислокаций [10, 11]. Цель - разработка математической модели диагностики распределения давления по длине трехфазного пьезоэлектролюминесцентного оптоволоконного датчика [12] по результатам замеряемой на торцевом сечении датчика интенсивности исходящего из оптоволоконной фазы света с учетом пьезоэлектролюминесцентного эффекта для случая нелинейной зависимости интенсивности света от действующего на электролюминофор пьезо- и управляющего электронапряжения. 1. Оптоволоконный пьезоэлектролюминесцентный датчик давления Рассмотрим датчик давления (рис. 1), который представляет собой составное трехфазное волокно: оптоволоконная фаза 1, электролюминесцентная 2 и пьезоэлектрическая с радиальной поляризацией 3 фазы, с внутренним фотопрозрачным (перфорированным) 4 и внешним 5 управляющими электродами. Действие напряжений на внешней боковой поверхности трехфазного волокна приводит к деформациям пьезоэлектрического слоя 3 и появлению в нем и в смежном электролюминесцентном слое 2 электрического поля, приводящего к свечению электролюминофора в оптоволокно 1; напряжение на управляющих электродах 4, Рис. 1. Составное трехфазное волокно с управляющим напряжением и интенсивностью света на выходе из оптоволокна под действием нагрузки Fig. 1. Compound three-phase fiber with the control voltage and intensity of light at the output of the optical fiber under loading 5 изменяет напряжение на электролюминесцентном слое 2 и величину интенсивности света на выходе оптоволокна. Такие волокнистые датчики давления могут быть использованы для диагностики напряженного состояния и дефектоскопии композитов, гидро- и аэродинамики, в медико-биологических исследованиях; датчики могут вплетаться в структуру полимерных волокнистых композитных элементов конструкций на этапе их создания в процессе намотки для осуществления контроля за напряженным состоянием и возможным появлением дефектов на структурном уровне композита во время их эксплуатации. 2. Математическая модель Осесимметричная математическая модель датчика давления (см. рис. 1) представляет собой трехфазный коаксиальный цилиндр, состоящий из расположенных соосно центрального оптоволокна (1-я фаза, ), электролюминесцентного (2-я фаза, ) и пьезоэлектрического (3-я фаза, ) слоев, разделенных внутренним перфорированным фотопрозрачным (между 1-й и 2-й фазами, ) и внешним ( ) управляющими электродами; толщины внутренних электродов считаем пренебрежимо малыми по отношению к характерным размерам сечения модели в цилиндрической системе координат . На границе с радиусом между 1-й и 2-й фазами (внутреннем электроде) задан потенциал , на внешней границе составного волокна с радиусом заданы потенциал и радиальное напряжение , напряжение на управляющих электродах , (1) задано значение осевой деформации для всех трех фаз и составного цилиндра в целом. Считаем, что на всех двух межфазных цилиндрических концентрических границах с радиусами , выполняются условия непрерывности радиальных перемещений , напряжений и дополнительно на границе 2-й и 3-й фаз выполняется условие непрерывности радиальной индукции электрического поля. Эти условия позволяют определить все константы интегрирования , , , , , , , , в общих решениях для каждой из фаз: для 1-й фазы при , , (2) для 2-й фазы при , , (3) для 3-й пьезоэлектрической фазы при (4) где коэффициенты , , , , , , компоненты трансверсально изотропных тензоров упругих , , , , пьезомеханических , и диэлектрических свойств 3-й фазы в цилиндрической системе координат , , z с учетом изотропии в плоскости и симметрии , , , [13, 14]. В результате по найденным константам интегрирования , (3) найдем искомое решение для электрического напряжения (5) на электролюминесцентном слое (см. рис. 1), которое рассчитывается через действующее на внешней границе напряжение и заданное управляющее варьируемое напряжение (1) между внутренним и наружным электродами составного волокна, потенциал (3) на межфазной границе в (5) находим через найденные из решения системы уравнений значения констант (3), для которых справедливы разложения , аналогичные (5). Константы , в (5) являются эффективными параметрами датчика и зависят от электроупругих свойств и геометрических параметров: радиуса оптоволокна и толщин электролюминесцентного и пьезоэлектрического слоев составного трехфазного волокна (3) и определяются на основе решения соответствующей системы линейных алгебраических уравнений относительно всех , …, неизвестных констант интегрирования. Аналогично (2)-(5) может быть получено решение для электроупругих полей в трехфазном слоистом волокне с внутренним и внешним управляющими электродами (см. рис. 1) при динамическом нагружении с использованием известных решений стационарной [15-17] и нестационарной [18] связанных задач электроупругости. 3. Алгоритм расчета функции распределения напряжений по интенсивности света на выходе оптоволокна Рассмотрим вероятностную модель, в которой распределение действующего на боковую цилиндрическую поверхность волоконного датчика внешнего напряжения (давления) по координате описываем случайным полем напряжения со свойствами статистической однородности и эргодичности на некоторой представительной длине; электрическое напряжение на электролюминофоре датчика также является случайной величиной, связанной с зависимостью (5). Найдем функцию плотности вероятностей распределения напряжения по анализу интенсивности света на выходе из оптоволокна датчика. Измеряемую величину интенсивности света на выходе из оптоволоконной фазы составного трехфазного волокна выразим через вероятности происхождения взаимообратных событий: , в виде неравенств: , соответственно. Так как из условия следует с учетом связи (5) для заданного значения управляющего напряжения (1) и исследуемых случайных величин и , то величина интенсивности света (6) связана с искомой функцией распределения для , (7) где напряжение , (8) максимальное значение интенсивности света на выходе из оптоволокна с насыщенной электролюминесценцией по всей длине при (рис. 2). После дифференцирования (6) по управляющему напряжению получим или (9) с учетом равенств (8), , , функция плотности вероятностей распределения напряжения Далее по заданной функциональной зависимости (см. рис. 2) (10) найдем производную , Рис. 2. Заданная «S-образная» (сплошная линия) и кусочно-постоянная аппроксимация (пунктирная линия) зависимости интенсивности светоизлучения электролюминофора от действующего на него напряжения Fig. 2. Given "S-shaped" (continuous line) and piecewise constant approximation (dotted line) of dependence of intensity of light-emission of an electrophosphor from voltage operating on it в результате имеем , (11) где (12) с учетом (8). После подстановки дифференциала (11) в (9) получим интегральное уравнение Фредгольма 1-го рода , (13) границы области интегрирования , . Интегральное уравнение Фредгольма (13) преобразуем к каноническому виду , (14) где разностное ядро , (15) известная (16) и искомая (17) функции с учетом новых переменных , , (18) границы интегрирования , или , (19) с учетом , (рис. 2). Различные аналитические и численные методы решения интегральных уравнений Фредгольма (14) даны в [19-23]. В частном случае, когда зависимость - ступенчатая функция (см. рис. 2), тогда производная в ядре уравнения (13) выражается через дельта-функцию Дирака следующим образом: , и из (12) получим или , (20) где с учетом (12), (13); искомая функция плотности вероятностей , (21) константа в (21) является эффективной характеристикой составного трехфазного волокна. 3. Решение тестовых прямой и обратной задач 3.1. Прямая задача. Для случая равномерного распределения (рис. 3, в) напряжения по длине внешней поверхности волокна с функцией плотности вероятностей (22) из уравнения (13) получим или (23) Рис. 3. Заданная функция (а) и найденные зависимости интенсивности света и производной на выходе из оптоволокна от управляющего напряжения (б); искомая функция плотности равномерного распределения напряжения по длине волокна (в) Fig. 3. Given function (a) and found light intensity dependences and derivative at the output of the optical fiber of the control voltage (b), the required probability-density function for uniform distribution of stress along the fiber (с) с параметром интегрирования , связью дифференциалов ; в результате из (23) следует решение , (24) где , находим по заданной функции на рис. 3, а. На рис. 3, б построены зависимости интенсивности света и производной (24) на выходе из оптоволокна от управляющего напряжения с параметрами датчика: , , границами действующего напряжения: , при заданном равномерном законе распределения (22). Пунктирные графики на рис. 3 приведены для случая, когда функция задана в виде ступенчатой функции (рис. 3, а) с параметрами , ; на рис. 3, б построены соответствующие этому случаю зависимости интенсивности света и производной на выходе из оптоволокна от управляющего напряжения при использовании для напряжения по длине волокна (рис. 3, в) равномерного закона распределения на интервале , где рассчитываются по значениям , управляющего напряжения (см. рис. 3, б) с учетом (21). 3.2. Обратная задача. При численном решении интегрального уравнения Фредгольма (14) в классе кусочно-постоянных функций определяем в узловых точках координаты «невязки» между заданными значениями функции и рассчитанными по известным фиксированным значениям ядра и варьируемым значениям искомой функции в узловых точках координаты ; граничные значения , , шаг по координате , . Распределение узловых точек , внутри соответствующих интервалов и искомые узловые значения находим из условия минимизации суммарных невязок . Решения обратной задачи в различных приближениях: начальном (Δ) и уточненном (○) при увеличенном числе узловых точек обозначены на рис.3,в пунктирными линиями в сравнении с «точным» решением (□), полученным ранее из решения прямой задачи; решения (Δ), (○) (см. рис. 3, в) найдены по заданной функции (см. рис. 3, а) электролюминофора, зависимости производной интенсивности света на выходе из оптоволокна от управляющего напряжения (□) (см. рис. 3, б). Заключение Разработана математическая модель пьезоэлектролюминесцентного оптоволоконного датчика давления, в котором механолюминесцентный эффект возникает в результате взаимодействия электролюминесцетного и пьезоэлектрического покрытий поверх оптоволокна, наличие дополнительных управляющих электродов расширяет возможности анализа распределения давления по длине волокна. Разработан алгоритм нахождения функции распределения по длине трехфазного оптоволоконного датчика значений напряжения , действующего на внешнюю цилиндрическую боковую поверхность, по результатам замеряемой на торцевом сечении датчика интенсивности света из оптоволоконной фазы для случая нелинейной зависимости интенсивности света от действующего на электролюминофор напряжения задача сведена к решению интегрального уравнения Фредгольма 1-го рода (14) с разностным ядром (15), зависящим от рассчитываемых эффективных параметров датчика и от производной заданной функции свечения электролюминофора. Для частного случая, когда зависимость - ступенчатая функция, получено аналитическое решение для функции плотности вероятностей напряжений (21); здесь ядро выражается через дельта-функцию и интегральное уравнение Фредгольма (13) сводится к алгебраическому (20). Для иллюстрации алгоритма решена «прямая» задача по нахождению из интегрального уравнения Фредгольма производной интенсивности света на выходе оптоволокна от управляющего напряжения для случая заданного равномерного закона функции плотности вероятностей распределения напряжения и далее для найденной функции решена «обратная» задача по нахождению функции плотности вероятностей ; численное решение обратной задачи осуществлено в различных приближениях, распределение узловых точек , внутри соответствующих интервалов и искомые узловые значения находили из условия минимизации суммарных невязок по заданным и рассчитываемым значениям функции .

### A A Pan’kov

Perm National Research Polytechnic University

# References

1. Бростилов С.А., Мурашкина Т.И., Пивкин А.Г. Перспективы использования волоконно-оптических датчиков давления для систем контроля и испытаний космической и авиационной техники // Математическое моделирование в машино- и приборостроении: сб. науч. тр. - Пенза: Изд-во Пенз. гос. ун-та, 2013. - № 10 (спец. вып.). - С. 55-66.
2. Пат. РФ № 2488772 C2. Способ неразрушающего контроля деталей из полимерных композиционных материалов / Ю.Ф. Огнев, О.Ш. Бердиев; опубл. 27.07.2013.
3. Бадеева Е.А., Бростилов С.А., Юрова О.В. Волоконно-оптический датчик давления на основе туннельного эффекта // Современная электроника. - 2011. - № 2. - С. 26-27.
4. Томышев К.А., Баган В.А., Астапенко В.А. Распределённые волоконно-оптические датчики давления для применения в нефтегазовой промышленности // Труды МФТИ. - 2012. - Т. 4, № 2. - С. 64-72.
5. Макарова Н.Ю. Моделирование выходного сигнала механолюминесцентного датчика динамического давления // Наука и образование науч. изд. - 2015. - № 6. - С. 187-200.
6. Tатмышевский К.В. Научные основы расчета и проектирования механолюминесцентных чувствительных элементов датчиков импульсного давления: автореф. дис. … докт. техн. наук. - М., 2010. - 33 с.
7. Real-time mechanoluminescence sensing of the amplitude and duration of impact stress / B.P. Chandra, V.K. Chandra, S.K. Mahobia, P. Jha, R. Tiwari, B. Haldar // Sensors and Actuators A: Physical. - 2012. - Vol. 173. - P. 9-16.
8. Detecting Mechanoluminescence From ZnS:Mn Powder Using a High Speed Camera / L. Kobakhidze, C.J. Guidry, W.A. Hollerman, R.S. Fontenot // IEEE Sensor Journal. - 2013. - Vol. 13. - No. 8. - Р. 3053-3059.
9. Chandra B.P., Chandra V.K. Dynamics of the mechanoluminescence induced by elastic deformation of persistent luminescent crystals // Journal of Luminescent. - 2012. - Vol. 132. - No. 3. - P. 858-869.
10. Банишев А.Ф., Банишев А.А., Лотин А.А. Исследование люминесценции и механолюминесценции мелкодисперсного порошка SrAl2O4:(Eu2+, Dy3+) в матрице полимера // Физика и химия обработки материалов. - 2012. - № 5. - С. 89-92.
11. Писаревский А.И., Татмышевский К.В., Голубев А.М. Сравнение особенностей механолюминесценции в кристаллах ZnS и (BaSr)Al2O4 // Наука и образование. МГТУ им. Н.Э. Бауман а. Электрон. журн. - 2012. - № 1. - URL: http://technomag.bmstu.ru/doc/297102.html (дата обращения: 01.05.2015).
12. Пат. РФ Волоконно-оптический датчик давления / А.А. Паньков; заявл. RU № 2016136058 от 06.09.2016.
13. Паньков А.А. Статистическая механика пьезокомпозитов. - Пермь: Изд-во Перм. гос. техн. ун-та, 2009. - 480 с.
14. Паньков А.А. Электромагнитоупругость и эффективные свойства пьезокомпозитов // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2012. - № 4. - С. 79-98.
15. Гринченко В.Т., Улитко А.Ф., Шульга Н.А. Электроупругость. - Киев: Наукова думка, 1989. - 280 с. - Механика связанных полей в элементах конструкций: в 5 т., Т. 5.
16. Партон В.З., Кудрявцев Б.А. Электромагнитоупругость пьезоэлектрических и электропроводных тел. - М.: Наука, 1988. - 472 с.
17. Дианов Д.Б., Кузьменко А.Г. Расчет цилиндрического пьезокерамического преобразователя, совершающего радиально-симметричные колебания // Акуст. журн. - 1970. - Т. 16, № 1. - С. 42-48.
18. Шляхин Д.А. Нестационарная осесимметричная задача электроупругости для анизотропного пьезокерамического радиально-поляризованного цилиндра // Изв. РАН. МТТ. - 2009. - № 1. - С. 73-81.
19. Михлин С.Г., Смолицкий X.Л. Приближенные методы решения дифференциальных и интегральных уравнений. - М.: Наука, 1965. - 384 с.
20. Верлань А.Ф. Сизиков В.С. Интегральные уравнения: методы, алгоритмы, программы. - Киев: Наукова думка, 1986. - 543 с.
21. Васильева А.Б., Тихонов А.Н. Интегральные уравнения. - 2-е, стереотип. изд. - М.: Физматлит, 2002. - 360 с.
22. Полянин А.Д., Манжиров А.В. Справочник по интегральным уравнениям. - М.: Физматлит, 2003. - 608 с.
23. Латыпов А.Ф. Численные методы решения линейных интегральных уравнений Фредгольма и Вольтерра 1-го рода // Обратные и некорректные задачи математической физики: материалы междунар. конф., посвященной 75-летию академика М.М. Лаврентьева, 20-25 августа 2007 г. - Новосибирск, 2007.

# Statistics

#### Views

Abstract - 121

PDF (Russian) - 63