УДК 681.7.068:531.8
ЧИСЛЕННО-ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ ВЛИЯНИЯ ТЕХНОЛОГИЧЕСКИХ НЕСОВЕРШЕНСТВ ГЕОМЕТРИИ СИЛОВЫХ СТЕРЖНЕЙ НА НАПРЯЖЕННОЕ СОСТОЯНИЕ И ОПТИЧЕСКИЕ ХАРАКТЕРИСТИКИ
СВЕТОВОДА ТИПА "ПАНДА"
© 2014 Е.В. Праведникова, Н.В.Семенов, А.Н.Труфанов, Н.А.Труфанов Пермский национальный исследовательский политехнический университет
Поступила в редакцию 24.11.2014
Проведён статистический анализ экспериментальных данных о геометрических параметрах поперечного сечения партии оптических волокон. Показано, что при изготовлении волокон не удается получить идеальное соответствие фактических параметров проектным значениям. Выявлены наиболее типичные формы сечений и отклонения. Исследование касается таких величин как диаметры волокна, силовых стержней и светопроводящей жилы, а также положение силовых стержней и све-топроводящей жилы относительно центра волокна. Установлено, что отклонения значений этих параметров варьируются в пределах 10%. Экспериментально замеренные величины двулучепреломле-ния, важной оптической характеристики волокна, связанной с механическими напряжениями, коррелируют с такими геометрическими параметрами, как диаметры силовых стержней и расстояния от центров стержней до центра волокна. В конечно-элементном пакете ЛМ8У8 проведён численный эксперимент по нахождению остаточных напряжений в процессе изготовления кварцевого оптического волокна типа "Панда" с учётом реальной формы геометрии волокна и его конструктивных элементов. Результаты представлены в виде графиков и аналитических зависимостей, которые соответствуют полученным экспериментально и дополняют их. Установлено, что при возможных отклонениях геометрии, величина двулучепреломления световода может изменяться в пределах 30%. Ключевые слова: анизотропное оптическое волокно, статистическая обработка результатов, остаточные напряжения, двулучепреломление.
ВВЕДЕНИЕ
В настоящее время оптическое волокно считается самой совершенной физической средой для передачи информации, в том числе на значительные расстояния [1-3]. Для изготовления высокоточных датчиков таких как волоконно-оптический гироскоп используется специальный тип волокна, так называемое анизотропное одномодовое волокно [4, 5], например, волокно типа "Панда". Такое волокно сохраняет состояние поляризации вводимого в него излучения. Эффект достигается за счёт создания в светопроводящей жиле анизотропии распределения полей остаточных напряжений. Это достигается путем введения в заготовку волокна специальных силовых стержней из кварцевого стекла, легированного небольшими добавками примесей, значительно увеличивающих коэффициент линейного температурного расширения
Праведникова Екатерина Владимировна, магистрант. E-mail: [email protected].
Семенов Никита Владимирович, аспирант кафедры вычислительной математики и механики. E-mail: [email protected]
Труфанов Александр Николаевич, кандидат технических наук, директор регионального центра информатизации. E-mail: [email protected].
Труфанов Николай Александрович, доктор технических наук, профессор, заведующий кафедрой вычислительной математики и механики. E-mail: [email protected]
и уменьшающих температуру стеклования [6, 7]. В процессе изготовления волокна после охлаждении от 2000°С до температуры окружающей среды в светопроводящей жиле волокна формируется необходимый уровень разности главных механических остаточных напряжений, который в рамках действия фотоупругих эффектов способствует изменению двулучепреломления.
При проектировании конструкции волокна и выбора параметров технологического процесса его изготовления предполагается, что волокно, его заготовка и конструктивные элементы волокна и заготовки имеют идеальные геометрические формы. У волокна типа "Панда" сечения светопроводящей жилы и силовых стержней должны иметь круглую форму, центры светопрово-дящей жилы и волокна совпадать, а центры силовых элементов находится на одинаковом расстоянии и на одной прямой, проходящей через центр волокна [8-10]. Однако технологический процесс изготовления волокна связан со сложным поведением материала в широком температурном диапазоне, при котором неоднородно легированный материал из вязкотекучего состояния переходит в стеклообразное. В таких условиях, когда остаточное напряженное состояние определяется совокупностью факторов, связанных со свойствами материалов и условиями технологического процесса изготовления, имеющиеся не-
совершенства геометрии могут влиять на механические и оптические характеристики волокна. В данной работе исследуется влияние отклонений некоторых геометрических параметров силовых стержней в составе готового волокна от их проектных значений.
ИЗМЕРЕНИЕГЕОМЕТРИЧЕСКИХПАРАМЕТРОВ ПОПЕРЕЧНОГО СЕЧЕНИЯ ВОЛОКНА
Геометрические параметры поперечного сечения каждого оптического волокна типа "Панда" фиксировались с помощью измерительной системы РК-2400, всего исследовано 1400 образцов. На рис. 1 представлена схема части кругового сечения волокна с указанием характерных замеряемых геометрических параметров.
Предполагается, что волокно с границей Г, стержни и светопроводящая жила имеют идеальное круговое сечение, а разбросы возможны в величинах диаметров стержней, жилы и в расположении их центров по отношению к проектным позициям. Точки О,£1зS2,С - центры волокна, первого, второго силовых стержней и светопро-водящей жилы соответственно; [) - диаметр волокна; D1, 02, Ос - диаметры первого и второго силовых стержней, а также светопроводящей жилы соответственно; Ц1, Ь2, Ьс - расстояния от центра волокна до центров силовых стержней и жилы; (р1,(2,(с - углы между отрезками, соединяющими центры силовых стержней, светопро-водящей жилы с центром волокна и основной осью (здесь и далее под основной осью понимаем прямую, параллельную отрезку, соединяющему центры силовых стержней, проходящую через центр волокна).
СТАТИСТИЧЕСКАЯ ОБРАБОТКА РЕЗУЛЬТАТОВ ИЗМЕРЕНИЙ
На основе измеренных данных получены статистические зависимости, характеризующие реальное распределение геометрических пара-
Г.
метров волокна "Панда" и его конструктивных элементов. Полученные значения представлены на графиках в условных безразмерных единицах вычисленных по формуле: У' = У/У^ , где У -анализируемое значение, Уср - среднее арифметическое значений соответствующей геометрической характеристики. При анализе параметров геометрических характеристик силовых стержней порядок стержней (первый или второй стержень) не определяется, поэтому Ог - объединённый набор величин D1,В2, Ц - объединённый набор величин Ь1,Ъ2, а ( - объединённый набор величин .
В математическом пакете 8ТЛТ18ТЮЛ подобраны законы распределения геометрических параметров, рис. 2-7. Из результатов подбора распределения для определенной выборки присваивается тот закон распределения из исследованных, который при найденных значениях его параметров дает минимальные значения статистики Колмогорова-Смирнова и статистики хи-квадрат.
Распределения диаметров волокна, силовых стержней и светопроводящей жилы стремятся к нормальному закону. При этом диаметры волокна [) варьируется в небольших диапазонах (<1%), (рис.2), а диаметры силовых стержней Di (рис.3) и диаметры светопроводящей жилы Ос (рис.4.) варьируются в пределах 10%.
Расстояние от центров силовых стержней до центров волокна Ц изменяется в пределах 10% по нормальному закону распределения (рис. 5).
Расстояние от центра светопроводящей жилы до центра волокна Ьс подчиняется логнормаль-ному закону (рис. 6).
На основе обработки экспериментально замеренных значений двулучепреломления световода В, установлено, что данная характеристика коррелирует с такими величинами как диаметры силовых стержней Di и расстояния от центров стержней до центра волокна Ц в соответствии с формулой вычисления нормированного коэффициента корреляции:
180
Рис. 1. Схема области поперечного сечения оптического волокна типа "Панда", содержащей конструктивные элементы
Рис. 2. Гистограмма распределения диаметра
волокна
D'
0,70 0.73 0.77 0.81 0.84 0,88 0.92 0.95 0.99 1.03 1.06 1.10 1.14 1.17 1.21 1.25 0.71 0,75 0.79 0.82 0.86 0.90 0,94 0.97 1.01 1,04 1,08 1.12 1.16 1.19 1,23
Рис. 3. Гистограмма распределения диаметров силовых стержней Di'
Рис. 4. Гистограмма распределения диаметра светопроводящей жилы D '
Рис. 6. Гистограмма распределения расстояний от центра светопроводящей жилы до центра волокна Ь
=
M (R - M (R ))M (Rm - M (Rm ))
O o
nm
где Rn, Rm,
- распределение величин, м - математическое ожидание, оп,ат - средние квадратичные отклонения для Rn и Rm соответственно.
Двулучепреломление В прямо пропорционально диаметрам силовых стержней ( = 0,56) и обратно пропорционально расстоянию от
Рис. 5. Гистограмма распределения расстояний от центров силовых стержней до центров волокна Ц '
центров силовых стержней до центра волокна ( ГБь = -0,22). По результатам статистического и корреляционного анализа геометрических параметров выделены диапазоны изменения значений величин влияющих на величину двулучеп-реломления.
ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ОСТАТОЧНЫХ НАПРЯЖЕНИЙ И ДВУЛУЧЕПРЕЛОМЛЕНИЯ
Рассмотрим задачу о нахождении остаточных напряжений, возникающих в результате охлаждения волокна до температуры окружающей среды после вытяжки из заготовки под действием высокой температуры (~ 2100°С ). Процесс связан с переходом материала из вязкоупругого состояния в стеклообразное. В статье [8] показано, что при решении задачи теплопроводности об охлаждении волокна в процессе вытяжки температура на поверхности и в центре волокна отличаются несущественно. Поэтому решение квазистатической краевой задачи (считаем, что процесс зависит от времени, но при этом скорости настолько малы, что силами инерции можно пре-
небречь) об изменении полей напряжений ведём с учётом предположения об однородности температуры в поперечном сечении волокна и зависимости от времени, установленной в [8]. Тогда задача о нахождении напряжённо-деформированного состояния волокна соответствует случаю обобщенной плоской деформации и включает уравнения равновесия:
дах (х, у, t) дтху (х y, t)
ддс
ду
= 0,
дТху (с, у, t) | да у (с, у, t) = о
дх
ду
х, у £ $,
геометрические соотношения:
£х (х y, t) =
дих (х, у, 0 , диу (х, y, О ———-—1, £у (х, у, t) = —--,
дх у ду
Уху (x, У, t) =
дих(х, у, t) диу(х, у, t)
ду дх
X, у £ $.
граничные условия: ах (X у, 1 )пх + Тху (X у, 1 )пу = 0,
Тху (X у, 1 )Пх + ау (X у, 1 )Пу = 0
ху £ Г , физические соотношения:
ах (х, у, 1) = (В + 3 О) £ (х, у, 1) - 8х (X, у, 1) -
- еТ (х у, 0) +(В - 2 °) (у (х у, 0 -8у (х у, 0-ет (х у, 0)
а у (х, у, 1) = (В - з О) £ (х, у, 1) - е„ (х, у, 1) -
-£т (х у, 0) ++(В+з С) (у (х у, 0 -£,у (х у, 0 -£т (Х у, 0),
£ 1 ^ (х, у, 1) = ~(ах-а), 1 г/
(Х у,1) = ¿в (Х у,1) - 3 6в (х у,1 №,
в силу
того, что в расплавленном состоянии материал предполагается несжимаемой жидкостью (бВ = 0), равен тензору вязких деформаций еВ = £В; 6В = £Вкк - объемная деформация. В -модуль объемного сжатия; О - модуль сдвига
Вследствие неоднородного легирования силовых стержней коэффициент линейного температурного расширения а и вязкость П являются функциями концентраций легирующих элементов и координат.
и (М( X, у ))
а = а(М(х, у)), п = пМх, у))еКТ(х,у1) , где и(^(х, у)) - энергия активации, ц - универсальная газовая постоянная, /л(х, у) - концентрация легирующих элементов. Известно [6], что влияние небольших степеней легирования кварцевого стекла на упругие характеристики незначительно.
В конечно-элементном пакете Л^УЗ проводится численный эксперимент по нахождению остаточных напряжений в процессе изготовления кварцевого оптического волокна типа "Панда", позволяющий задавать реальную форму геометрии волокна и его конструктивных элементов.
В результате вычислительных экспериментов выявлено, что при отклонении формы сечения волокна от проектной величины количественно напряжения изменяются незначительно. Так, например, на рис.7. приведены распределения интенсивностей напряжений при разных отклонениях формы геометрии и разница между максимальными значениями составляет не больше 3%. Меняется качественное распределение напряжений по сечению, в том числе и в области светопроводящей жилы, что влечёт за собой возможность изменения величины двулучепрелом-ления В, которое в данном случае вычисляется по формуле:
£вТ (Ху,0 = ~(ау-а), X,у £ $ ,
где 1 - время, п - внешняя нормаль к Г ,
Т
£Т (х, у, I) = |а(Т)ёТ - температурная деформа-
То
ция, £х (х, у, 1), £ву (х, у, 1) - компоненты тензора вязких деформаций ёв. а = акк / 3 - среднее
напряжение, Е - единичный тензор второго ранга. Девиатор тензора вязкой деформации
В =
ю 2п
С(я) I'(ах ((,()- ау (г' (^{г? ( V)2 гМ( Х
0 0
ио ¿.Л
|||Е (г (V )2 тйтй(
00
где Г, (
соответствующие оси цилиндрической системы координат, а х и а у - нормальные на-
X У , ч
пряжения в поперечном сечении волокна, С (к ) -фотоупругая постоянная для данной длины волны, |Е (г, (, V ) - распределение интенсивности основной моды по поперечному сечению све-
Х
¿7. t , Па
<*Ы >Па
а) б)
Рис. 7. Интенсивность напряжений <7Ы в волокне: а - проектная форма; б - форма с отклонениями
B ■ 10-
B ■ÍO
-4
4
а) б)
Рис. 8. Зависимость модового двулучепреломления световода В от значений а - диаметров силовых стержней Д ' и £>2'; б - расстояний от центров силовых стержней до центра волокна Ь1 и Ь2
товода. В данном исследовании принималось равномерное распределение E (r, ф ,V )= E0 .
Выполнены численные исследования напряженного состояния оптического волокна Panda и определены величины модового двулучепреломления при установленных наиболее вероятных разбросах диаметров силовых стержней D1, D2 и расстояний от центров стержней до центров волокна L1 и L2 (рис.8).
Методом наименьших квадратов получены удобные в инженерных целях аппроксимирующие выражения поверхностей полными кубическими полиномами.
B(DD2') = boo + bwDi'+ boiD2'+ b2oDi '2+buDi' D2'+ b
'+ b02D2 '2 + b30D1 '3 + b21D1 '2 D2'+ b12D1' D2 '2 + b03D2 '3
b00 = 0,03249; b10 = -3022; b01 = -3022; b20 = 1,178E + 08; b11 = 1,405E + 08 ; b02 = 1,178E + 08; b30 = -1,639E + 12 ; b21
B(L1', L2') = b00 + bwL1'+ b01L2'+ b20L1 '2 + bnL1' L2'+«
+ b02 L2 1 + b30L1' + b21L1' L2'+ bxlL 1 L2 + b03 L2
b00 = 0,03924; b10 = 607,7 ; b01 = 607,7; b20 = 5,263E + 07; b11 = 7,931E + 06; b02 = -5,263E + 07; b30 = 1,28E +12;
-7,839E +10; b12
-7,839E +10 ; b03 = -1,28E + 12
-2,211E +12; b12
-2,211E +12; b03 = -1,639E + 12
Установленный численно нелинейный характер зависимости отражает установленную ранее экспериментально корреляционную зависимость, состоящую в том, что увеличение диаметра силовых стержней при фиксированном положении их центра приводит к увеличению величины двулучепреломления, а увеличение расстояний от центров силовых стержней до центра волокна Ьх и приводит к уменьшению величины двулучепреломления. При этом разброс величины двулучепреломления при разных реализациях геометрии достигает 30%.
ЗАКЛЮЧЕНИЕ
Таким образом, путем статистической обработки результатов экспериментальных измерений отклонений некоторых геометрических параметров силовых стержней в составе готового волокна от их проектных значений, а также численного исследования влияний этих отклонений на остаточные напряжения в анизотропных оптических волокнах типа "Панда", установлено существенное изменение двулучепреломления (разброс до 30%) - важной оптической характеристики волокна.
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований, гранта РФФИ проект № 13-08-96036р_ура-л_а и 14-08-31530мол_а.
СПИСОК ЛИТЕРАТУРЫ
1. Гроднев И.И., Ларин Ю.Т., Теумин И.И. Оптические кабели: конструкции, характеристики, производство и применение. М.: Энергоатомиздат, 1991. 264 с.
2. Волоконно-оптические датчики [пер. с яп.] / Т. Окоси, К. Окамото, М. Оцу, Х. Нисихара, К. Каюма, К.
Хататэ. Л.: Энергоатомиздат. Ленингр. отделение, 1990. 256 с.
3. Волоконная оптика и приборостроение / М.М. Бутусов, С.Л Галкин., С.П. Оробинский, Б.П. Пал. Л.: Машиностроение, 1987. 328 с.
4. Varnham M., Payne D., Barlow A. Analytic solution for the birefringence produced by thermal stress in polarization-maintaining optical fibers // Journal of lightware technology. 1983. V. LT-1, №2. P.332-338
5. Stress birefringence analysis of polarization maintaining optical fibers / R. Guan, F. Zhu, Z. Gan, D. Huang, S. Liu // Optical Fiber Technology. 2005. V. 11. № 3. P. 240-254.
6. Леко В.К., Мазурин О.В. Свойства кварцевого стекла. Л.: Наука, 1985 г. - 166 с.
7. Бурков В.Д., Иванов Г.А. Физико-технологические основы волоконно-оптической техники: Учеб. пособие. М.: ГОУ ВПО МГУЛ, 2007. 222 с.
8. Остаточные напряжения в заготовках силовых стержней оптического волокна / Н.А. Труфанов, О.Ю. Сметанников, А.Н. Труфанов, И.И. Крюков // Вестник ПГТУ. Технологическая механика. 2002. С.110-115.
9. Numerical analysis of residual stresses in preform of stress applying part for PANDA-type polarization maintaining optical fibers / A.N. Trufanov, O.Yu. Smetannikov, N.A. Trufanov // Optical Fiber Technology. 2010. V. 16. N 3. P. 156-161.
10. Сметанников О.Ю., Труфанов Н.А. Технологические и остаточные напряжения в неоднородном стеклующемся цилиндрическом стержне / / Механика композиционных материалов и конструкций. 2009. № 2. С. 126-140.
NUMERICAL AND EXPERIMENTAL STUDY OF INFLUENCE TECHNOLOGICAL IMPERFECTIONS OF GEOMETRIC PARAMETERS OF STRESS APPLYING PARTS ON THE STRESS STATE AND OPTICAL CHARACTERISTICS OF PANDA FIBER
© 2014 E.V. Pravednikova, N.V. Semenov, A.N. Trufanov, N.A. Trufanov
Perm National Research Polytechnic University
The experimental data of the geometry of cross sections of the party of Panda fiber were statistically analyzed. The analysis shown, that in the fiber manufacture it is not obtained to achieved a complete capability of geometry with the project design. Found the most typical forms of cross sections and deviations. The study concerned with such criteria as diameter of fiber, location of stress applying parts and fiber core relative to the center of fiber. It was found that deviation of these parameters ranges within 10%. Shown that the experimentally measured birefringence correlates to such geometrical parameters as stress rod diameter and distance between centers of rod and fiber. Was conducted a numerical experiment by finite element modeling with Ansys and were found residual stresses during the manufacturing of Panda fiber with account of the real form of fiber geometry and its constructive elements. The results are presented in graphs and analytical relationships that are consistent with the experimental data and complement them. It was found that with possible geometry deviations, the range of birefringence may vary within 30%. Key words: anisotropic optical fiber, statistical analysis of the results, residual stress, birefringence.
Ekaterina Pravednikova, Undergraduate at the Computational Mathematics and Mechanics Department. E-mail: [email protected]
Nikita Semenov, Postgraduate at the Computational Mathematics and Mechanics Department. E-mail: [email protected]
Aleksandr Trufanov, Candidate of Technics, Associate Professor at the Computational Mathematics and Mechanics Department, Director of the Regional Center of Informatization. E-mail: [email protected]
Nikolay Trufanov, Doctor of Technics, Professor, Head at the Computational Mathematics and Mechanics Department. E-mail: [email protected]