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). Для иллюстрации алгоритма решена «прямая» задача по нахождению из интегрального уравнения Фредгольма производной интенсивности света на выходе оптоволокна от управляющего напряжения для случая заданного равномерного закона функции плотности вероятностей распределения напряжения и далее для найденной функции решена «обратная» задача по нахождению функции плотности вероятностей ; численное решение обратной задачи осуществлено в различных приближениях, распределение узловых точек , внутри соответствующих интервалов и искомые узловые значения находили из условия минимизации суммарных невязок по заданным и рассчитываемым значениям функции .

About the authors

A A Pan’kov

Perm National Research Polytechnic University

References

Statistics

Views

Abstract - 150

PDF (Russian) - 80