Задача выбора оптимального закона неоднородного легирования силового стержня для заготовки анизотропного кварцевого волокна типа Панда Текст научной статьи по специальности «Физика»

УДК 681.7.068: 539.319

ЗАДАЧА ВЫБОРА ОПТИМАЛЬНОГО ЗАКОНА НЕОДНОРОДНОГО ЛЕГИРОВАНИЯ СИЛОВОГО СТЕРЖНЯ ДЛЯ ЗАГОТОВКИ АНИЗОТРОПНОГО КВАРЦЕВОГО ВОЛОКНА ТИПА ПАНДА

© 2013 Н.В. Семенов, А.Н. Труфанов, Н.А. Труфанов

Пермский национальный исследовательский политехнический университет

Поступила в редакцию 26.11.2013

Рассматривается задача поиска оптимальной зависимости легирования оксидом бора силового стержня анизотропного оптического волокна типа Панда. Задача оптимизации ставится с целью увеличения общей температурной деформации силового стержня при выполнении условий прочности на всех технологических этапах производства. При решении используется метод Нелдера-Мида. Ограничения по условиям прочности учитываются с помощью метода штрафных функций. Параметрами оптимизации являются коэффициенты функции, описывающей закон распределения по радиусу стержня оксида бора в чистом кварце. Зависимость распределения аппроксимируется двумя кубическими сплайнами. Получен оптимальный закон распределения легирующих добавок с учётом условий прочности.

Ключевые слова: силовой стержень, оптическое волокно Panda, технологические напряжения, закон легирования, оптимизация

Изготовление высокоточных приборов, таких, как волоконно-оптические датчики, накладывает серьёзные требования на уровень свето-проводящих характеристик волокна. В таких устройствах широко применяются анизотропные кварцевые оптические волокна [1-3], поддерживающие поляризацию сигнала за счет разности механических напряжений в светопроводящей жиле, например, волокна типа Панда. Особенность этих волокон заключается в наличии двух цилиндрических стержней, устанавливающихся в заготовке световода на одном диаметре по обе стороны от светопроводящей жилы. Стержень изготавливается с помощью модифицированного метода химического осаждения из газовой фазы (MCVD) [4] и содержит переменные вдоль радиуса легирующие добавки, назначение которых заключается в увеличении температурной деформации стержня по сравнению с материалом заготовки волокна для создания в светопро-водящей жиле максимальной разности главных напряжений, определяющих оптические свойства

Семенов Никита Владимирович, аспирант Труфанов Александр Николаевич, кандидат технических наук, директор Областного центра новых информационных технологий. E-mail: [email protected] Труфанов Николай Александрович, доктор технических наук, профессор, заведующий кафедрой вычислительной математики и механики. E-mail: [email protected]

волокна. Неоднородное и достаточно высокое легирование стержня приводит к возникновению больших технологических напряжений при изготовлении стержня, которые зачастую приводят к растрескиванию заготовки [5]. Задача выбора оптимального закона легирования ставится с целью подбора необходимых параметров, благодаря которым удастся увеличить общую температурную деформацию силового стержня при выполнении условий прочности на всех технологических этапах производства. При этом известно, что температурная деформация кварцевого стекла увеличивается с ростом концентрации легирующего оксида бора В20з. Задача оптимизации решается методом Нелдера-Мида [6] на основе данных, получаемых в процессе решения термомеханической задачи о поведении стержня при остывании от температуры 2100оС до комнатной температуры с учётом наличия релаксационного перехода из стеклообразного в вязкотекучее (стеклования). Ограничения по условиям прочности учитываются с помощью метода штрафных функций. Параметрами оптимизации являются коэффициенты функции, описывающей закон распределения по радиусу стержня оксида бора в чистом кварце.

Постановка задачи оптимизации. Силовой стержень представляет собой длинный круглый цилиндр радиусом г2=5 10-3 м. Стержень имеет переменную по радиусу степень

легирования ^(г) чистого кварца оксидом бора В2О3 на расстоянии от центра до радиуса Г\=Ъ 10-3 м, а остальная часть стержня состоит из чистого оксида кварца. При этом содержание легирующих добавок не превышает 10%. В качестве функции, характеризующей величину температурной деформации, используется интегральное относительное содержание оксида бора в элементе стержня единичной длины

главному напряжению о1 и интенсивности напряжений ог):

ф(Рl,•••,РктаХ к (Рl,•••, Рк Ык к (Рl,•••, Рк Ык

'1

Ф(Р1 ^ Рк ) = 2ж\ ^(Рх ^ Рк ,'

(1)

Параметрами оптимизации р1,...,рк являются коэффициенты функции, описывающей закон распределения оксида бора в чистом кварце.

Задача оптимизации может быть сформулирована следующим образом: найти параметры распределения, обеспечивающие максимальное значение функции (1) (температурную деформацию силового элемента) при ограничениях типа неравенств (условия прочности по первому

и ограничениях типа равенств (уравнения краевой задачи по определению остаточных напряжений при изготовлении стержня с законом легирования ^(г); где [ор ]- предел прочности на растяжение, [ас ]- предел прочности материала на сжатие.

Для преобразования сформулированной задачи с ограничениями типа неравенств и равенств к задаче минимизации без ограничений типа неравенств использован метод штрафных функций [7]. Целевая функция с учетом (1) примет следующий вид:

Р ^х^^ Рк ) = -ф(Рl,•••, Рк ) +

+ 4 (к ]-кГх (р,,..., Рк ))2-(1 - ]-кГх (р,,..., Рк ))))+

(([к ]-кГ (р,,..., Рк ))2-(1 - ]-кГх (а,..., Рк

+ 4

где огтах и о1тах - максимальные значения в конструкции интенсивности остаточных напряжений и первого главного остаточного напряжения соответственно, А,2 - весовые коэффициенты для штрафных функции.

Ограничения типа равенств на минимизируемую функцию. Ограничения типа равенств на компоненты тензора остаточных напряжений формулирует математическая модель термомеханических явлений, происходящих при остывании стержня из кварцевого стекла от температур выше температуры стеклования 2100оС до комнатной температуры, которая включает в себя краевую задачу нестационарной теплопроводности и краевую задачу о напряженно-деформированном состоянии (НДС) с учетом

(2)

неоднородного поля температур и наличия релаксационного перехода (стеклования) в остывающем материале. Постановка краевой задачи нестационарной теплопроводности в области £ с границей I содержит: - уравнение теплопроводности

ят

рс (Т)— = О» (4(т)- вгаё (т))

х е 5

где с(Т), А,(Т), р - соответственно теплоемкость, теплопроводность и плотность материала. Дис-сипативный разогрев полагается пренебрежимо малым из-за низких скоростей деформаций;

- граничные условия

-4(Т ) вгаё (Т) - п = а(Т) - (Т - Т (*)) + в^ (Т )4, х е 1Ъ

где первое слагаемое правой части описывает конвективный теплоперенос, а второе - теплообмен излучением по закону Стефана-Больцмана; в - коэффициент черноты, о0 - постоянная Стефана-Больцмана, а(Т) - коэффициент теплопередачи, Тс() - температура окружающей среды, п - нормаль к границе I;

- начальные условия Несвязанная краевая задача о НДС с учетом

малости деформаций включает:

Т ( 0 х ) = Т0 ( х) х е 5

<

0

- уравнения равновесия

Шу (ст ) = 0, х е 5,

- геометрические соотношения Коши

8 * 2+ )Т), х е 5,

где и(х,£) - вектор перемещений, 8 (х, t) - тензор полных деформаций.

- граничные условия в перемещениях и напряжениях

и ( х ) = и ( х), х е 1а, ст ■ п = Р (х) х е 1ст

где 10 - части границы с заданными перемещениями и нагрузками соответственно.

Общая система уравнений задачи о напряженно-деформированном состоянии заготовки включает также определяющие соотношения. Механическое поведение кварцевого стекла в широком диапазоне температур (20-2100°С) связано с переходом материала из стекло-образного в вязкотекучее состояние. Для описания механического поведения кварцевого стекла с целью адекватного отражения указанных явлений были использованы физические соотношения максвеловского типа:

/у 4/А ^ ^ \

СТ= С ■■(8 -8Т -8в )

Т

8Т (х, t) = Е| а(х, Т(х, t))ЛТ

58 в _ 5 8 В =

I

I

58 В

5t

Ж

5t г|(Т) Л(х, Т) = Ло (хУТ{х']

упругих свойств при введении малых концентраций легирующих элементов пренебрегаем); 8 е (х, t) = 8(х, t) — 8 Т (х, t) — 8 в (х, t) - тензор упругих деформаций, 8(х, t) - тензор полных деформаций; 8 Т (х, t) - тензор температурных деформаций; 8 в (х, t) - тензор вязких деформаций; Л(х,Т), а(х,Т) - вязкость и коэффициент температурного расширения материала, зависящие от температуры и от неоднородно распределенных по объему тела легирующих примесей; ^ -начальная температура, при которой предполагается отсутствие в теле начальных напряжений, деформаций и их производных по времени (гипотеза о естественном ненапряженном и неде-формированном состоянии;

5(х, t ) = ст (х, t) —ст(x,t )Е - девиатор тензора напряжений; ст(х, t) = ст /3 - среднее напряжение; ЕЕ - единичный тензор второго ранга; R -универсальная газовая постоянная; ?у(ц(х)) -энергия активации. Девиатор тензора вязкой деформации ев (х, t) = 8 в (х, t) — 10в (х, t )Е, в силу

того, что в расплавленном состоянии материал предполагается несжимаемой жидкостью (бв = о), равен тензору вязких деформаций

=8в; ® в = 8 вм - объемная деформация.

Методика численного решения задачи определения напряженно-деформированного состояния изложена в работе [8].

Решение задачи оптимизации. Для поиска минимума функционала (2) с изложенными ограничениями использован метод Нелдера-Мида. В данной работе зависимость распределения степени легирования от радиуса была аппроксимирована двумя кубическими сплайнами. Таким образом, рассматривались функции распределения легирующего элемента ^(г) вида:

где 4( - тензор четвертого ранга упругих констант застеклованного материала (изменением

Ц (г )

N (г )Ц1+щ (г к+ыт (г к+ыт (г к, К (гК + М'т (гК + N (гк + Щ (гк,г е

-е [0,0.0015] [0.0015,0.003]

(3)

где ц1, рда, - значения концентрации легирующих элементов в центре стержня, в середине легированной зоны (г=г1/2), в конце легированной зоны (г=г1) соответственно; Ц , Цт, ЦП -значения производных d\x¡dr в этих же точках.

В результате имеем шесть независимых параметров оптимизации:

Т

0

Р1 =^1, Р2 = Ц, Р3 = Ц Р5 = Ц n Рб

m

Р4 =Ц

m,

n

о, использованы следующие [ст р ] = 35 МПа, [стс ]= 60 МПа. Последовательные результаты решения задачи

напряжении значения:

В качестве ограничении по нормальным растяги- оптимизации отражены на рис. 1. вающим напряжениям о1 и по интенсивности

Рис. 1. Изменение зависимости степени легирования силового стержня от радиуса в процессе решения задачи оптимизации, на итерациях № 1, 50, 100, 200 и 241

В результате решения поставленной задачи была определена зависимость степени легирования от радиуса (рис. 1. кривая «241»), позволяющая реализовать в силовом стержне для оптического волокна типа Панда максимальные температурные деформации при соблюдении условий прочности. Результат поиска оптимальных параметров формулы (3) имеет вид:

Р = Ц = 0,0869, р2 = Ц, = 0,0299, Р3 = Ц т = 0,0977, Р4 = = 0,0084, Р5 = Цп = 0,1037, Р6 = Цп = 0,0413.

Остаточные технологические напряжения в силовом стержне для этой зависимости распределения легирующих элементов приведены на рис. 2.

Выводы: показано, что при изготовлении методом MCVD силовых стержней с переменным содержанием легирующих добавок по радиусу для оптического волокна типа Панда, возможно увеличение величины температурной деформации стержня и, соответственно, улучшение оптических свойств волокна, с помощью использования оптимального закона легирования.

Рис. 2. Зависимость напряжении стг (радиальные напряжения), ст (окружные напряжения), <стг (осевые напряжения), <Ji (интенсивность напряжении) в силовом стержне для оптического волокна типа Panda от радиуса в случае оптимального легирования

На основе решённой задачи оптимизации можно сделать вывод, что используя полученный закон распределения легирующих добавок можно увеличить величину температурнои деформации стержня по сравнению с ранее полученными аналитическими зависимостями, сохранив его прочность. Полученное оптимальное распределение степени легирования от радиуса

существенно отличается от используемых на практике степенных зависимостей вида [8]:

('')

= a

1 -

С \

r

\ r0 У

(4)

Так, например, параметр, отражающий величину температурной деформации (1) для распределения вида (4), равен Ф^ь.р^^ 10-7 при значении степенного аргумента b=4 и a=0,1, а для полученной оптимальной зависимости этот параметр равен Ф(рь.. ,,pk)=13,82 10-7. То есть, при оптимальном легировании удается увеличить общее содержание оксида бора в стержне почти в 1,5 раза при сохранении прочности изделия.

Работа выполнена при финансовой поддержке гранта РФФИ проект №13-08-96036р_урал_а.

СПИСОК ЛИТЕРАТУРЫ:

1. Guan, R Stress birefringence analysis of polarization maintaining optical fibers / R. Guan, F. Zhu, Z. Gan et

al. // Optical Fiber Technology. 2005. Vol. 11, №3. P. 240-254.

Окоси, Т. Волоконно-оптические датчики / Т. Окоси, К. Окамото, М. Оцу и др. Пер. с япон. - Л.: Энергоатомиздат. Ленингр. отделение, 1990. 256 с. Гроднев, И.И. Оптические кабели: конструкции, характеристики, производство и применение / И.И. Гроднев, Ю.Т. Ларин, И.И. Теумин. - М.: Энергоатомиздат, 1991. 264 с.

Бурков, В.Д.Физико-технологические основы волоконно-оптической техники: учеб. пособие / В.Д. Бурков, Г.А. Иванов. - М.: ГОУ ВПО МГУЛ, 2007. 222 с.

Trufanov, A.N. 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. Vol. 16, №3. P. 156-161. Банди, Б. Методы Оптимизации. Вводный курс. -М.: Радио и связь, 1988. 42 с. Реклейтис, Г. Оптимизация в технике: в 2 кн. / Г. Реклейтис, А. Рейвиндран, К. Рэгсдел. - М.: Наука, 1980. 226 с.

Сметанников, О.Ю. Технологические и остаточные напряжения в неоднородном стеклующемся цилиндрическом стержне / О.Ю. Сметанников, НА. Труфанов // Механика композиционных материалов и конструкций. 2009. №2. С. 126-140.

THE PROBLEM OF CHOICE THE OPTIMUM LAW OF NONUNIFORM ALLOYING THE POWER CORE FOR PREPARATION THE ANISOTROPIC QUARTZ FIBRES OF PANDA TYPE

© 2013 N.V. Semenov, A.N. Trufanov, N.A. Trufanov Perm National Research Polytechnical University

The problem of searching the optimum dependence of alloying by boron oxide the power core of anisotropic optical fiber of Panda type is considered. The task of optimization is set for the purpose of increase the common temperature deformation of a power core at performing the conditions of durability at all technological production phases. At the decision the method of Nelder-Meed is used. Restrictions under the terms of durability are considered by means of penal functions method. Parameters of optimization are coefficients of the function describing the distribution law on radius of a core of boron oxide in clear quartz. Dependence of distribution is approximated by two cubic splines. The optimum distribution law of alloying additives taking into account durability conditions is received.

Key words: power core, optical fiber of Panda type, technological tension, alloying law, optimization

Nikita Semenov, Post-graduate Student

Alexander Trufanov, Candidate of Technical Sciences, Director of the Regional Center of New Informational Technologies. E-mail: [email protected]

Nikolay Trufanov, Doctor of Technical Sciences, Professor, Head of the Department of Computational Mathematics and Mechanics. E-mail: [email protected]