CC BY
89
Галягин К.С., Ошивалов М.А., Савин М.А. Учет пьезооптических эффектов при моделировании теплового дрейфа волоконно-оптического гироскопа // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2015. - № 4. - С. 55-71. DOI: 10.15593/perm.mech/2015.4.04
Galyagin K.S., Oshivalov M.A., Savin M.A. Simulation of thermal drift of fiber-optic gyroscope taking into account piezo-optical effects. PNRPUMechanics Bulletin. 2015. No. 4. Рр. 55-71. DOI: 10.15593/perm.mech/2015.4.04
ВЕСТНИК ПНИПУ. МЕХАНИКА № 4,2015
PNRPU MECHANICS BULLETIN
http://vestnik.pstu.ru/mechanics/about/inf/
DOI: 10.15593/perm.mech/2015.4.04 УДК 621.372.82:621.383.8-752.4:536.7
УЧЕТ ПЬЕЗООПТИЧЕСКИХ ЭФФЕКТОВ ПРИ МОДЕЛИРОВАНИИ ТЕПЛОВОГО ДРЕЙФА ВОЛОКОННО-ОПТИЧЕСКОГО ГИРОСКОПА
К.С. Галягин, М.А. Ошивалов, М.А. Савин
Пермский национальный исследовательский политехнический университет, Пермь, Россия
о СТАТЬЕ
АННОТАЦИЯ
Получена: 1 июля 2015 г.
Принята: 13 октября 2015 г. Опубликована: 25 декабря 2015 г.
Ключевые слова: волоконно-оптический гироскоп, фазовая невзаимность, дрейф, эффект Саньяка, разностные методы, схемы укладки, пьезооптика, термоупругость, коэффициент преломления, главные напряжения, число витков
Проблема прогноза и компенсации погрешностей волоконно-оптических гироскопов, вызванных влиянием внешних факторов, существует довольно давно, и до сих пор является актуальной. Особое внимание к ней проявилось в связи с необходимостью увеличить точность приборов для установки в прецизионных устройствах. В статье подробно рассматриваются теоретические аспекты возникновения термически индуцированной фазовой невзаимности волоконно-оптического гироскопа (ВОГ). Вкратце поясняется основной принцип работы ВОГ, а также определяется проблема исследования. Содержится базовая классификация дрейфов ВОГ по причинам, их вызывающим. Приводятся основные определяющие соотношения пьезооптических эффектов, возникающих в кварцевом волокне, кроме того, демонстрируется процедура получения расчетного теплового дрейфа ВОГ. Суть методики расчета сведена к разрешению дифференциального уравнения движения разностным методом. На основании соотношений пьезооптики строится итерационный процесс для расчета времени хода лучей по оптическому контуру ВОГ. В работе также подробно излагается поэтапная методика расчета дрейфа с привлечением как сторонних программных продуктов, так и авторского программного кода. Устанавливается взаимосвязь и очередность выполнения расчетов для получения результата. Так, задача термоупругости ВОГ разрешается в инженерном пакете STAR-CCM+, а собственно расчет дрейфа ведется независимо в среде MATLAB. Приводится и анализируется качественная картина дрейфа ВОГ, полученная прямым численным моделированием без привлечения алгоритмов анализа сигнала. Преимуществом вынесенного отдельно постпроцессора является многократная обработка исходных данных, полученных из расчета термоупругости, любыми методами и способами, запрограммированными пользователем в постпроцессор.
©ПНИПУ
© Галягин Константин Спартакович - кандидат технических наук, доцент, e-mail: abins@pstu.ru Ошивалов Михаил Анатольевич - кандидат технических наук, доцент, e-mail: abins@pstu.ru Савин Максим Анатольевич - аспирант, старший преподаватель, e-mail: abins@pstu.ru
Konstantin S. Galyagin - PhD in Technical Sciences, Associate Professor, e-mail: abins@pstu.ru Michael A. Oshivalov - PhD in Technical Sciences, Associate Professor, e-mail: abins@pstu.ru Maxim A. Savin - Doctoral Student, Assistant Professor, e-mail: abins@pstu.ru
55
Galyagin K.S., Oshivalov M.A., Savin M.A. /PNRPUMechanics Bulletin 4 (2015) 55-71
SIMULATION OF THERMAL DRIFT OF FIBER-OPTIC GYROSCOPE TAKING INTO ACCOUNT PIEZO-OPTICAL EFFECTS
K.S. Galyagin, M.A. Oshivalov, M.A. Savin
Perm National Research Polytechnic University, Perm, Russian Federation
ARTICLE INFO
ABSTRACT
Received: 1 July 2015 Accepted: 13 October 2015 Published: 25 December 2015
Keywords:
fiber-optic gyroscope, phase independence, drift, Sagnak effect, difference method, packing diagram, piezo-optics, thermoelasticity, refraction coefficient, principle stresses, number of turns
The prediction and compensation of the fiber-optic gyro errors caused by the influence of the environmental factors has been the subject of many studies for a long time, but still remains an urgent problem. In the past few years an interest in this problem has been increased due to a need in improving the precision of instruments incorporated into sensitive equipment. The theoretical aspects of thermally induced phase non-reciprocity of fiber-optic gyroscope phase (FOG) are considered in detail by the authors. The operating principle of the fiber-optic gyros is outlined and the main objective of the study is formulated. The paper presents the basic classification of the FOG drifts based on the drift root causes. The main constitutive relations for piezo-optical effects occurring in silica fiber are derived and the procedure of evaluating the design value of the FOG thermal drift is described. The key idea of the computational procedure is the development of solution to the differential equation of motion by the finite difference method. The peizo-optical relations are used to elaborate the iterative scheme to calculate the time of the ray traveling through the FOG optical loop (optical path length). The paper also presents a detailed (step-wise) procedure for the temperature drift computation, which employs both the third-party programs and authors' codes. The interrelation of computational steps and the precedence of operations to be made for obtaining the required estimates have been determined. Thus, the problem of FOG thermoelasticity was solved based on the STAR-CCM+ applied package, whereas the drift computation was done separately in the MATLAB environment. A qualitative picture of the FOG drift obtained by direct numerical simulation without using the algorithms for signal analysis has been obtained and analyzed. The advantage of the externally connected processor is multiple processing of the initial data obtained from thermoelastic computations by any of the methods and algorithms which have been programmed by users in the processors.
©PNRPU
Введение
В настоящее время в гироскопических и навигационных системах широко применяются волоконно-оптические гироскопы (ВОГ). Они малы, надежны, точны, лишены недостатков, присущих инерциальным гироскопам, таких как масса движущихся частей, склонных к изнашиванию, а также существенное время готовности устройства (здесь время, затраченное на разгон маховика) [1]. Отсутствие подвижных деталей объясняется устройством ВОГ или же его минимальной конфигурацией. Волоконно-оптический гироскоп в своей минимальной конфигурации включает в себя источник света (суперлюминесцентный диод, лазер или активные эрбиевые волокна), светоделитель, модулятор, оптоволоконный контур, собиратель и детектор [2]. Импульс света, проходя светоделитель, делится на два в идеале абсолютно одинаковых луча, которые, встречно распространяясь, обходят оптоволоконный контур, затем собираются и регистрируются детектором. Далее гироскопом обрабатывается полученный сигнал, и на основе зарегистрированной интерференционной картины от световых лучей рассчитывается значение угловой скорости [3-5]. В сущности, это и есть так называемый эффект Саньяка.
Но за всеми положительными качествами ВОГ скрывается довольно существенный недостаток - дрейф показаний ввиду воздействия внешних факторов.
Так как принцип работы гироскопа основан на эффекте Саньяка, устанавливающем прямую пропорциональную зависимость между скоростью вращения контура и фазовым набегом, различные паразитные дрейфы тесно с ним связаны.
56
Галягин К.С., Ошивалов М.А., Савин М.А. /Вестник ПНИПУ. Механика 4 (2015) 55-71
Q = cXA(p/(8 nSv), (В1)
где с - скорость света; X - длина волны излучения; Дф - разность фаз встречно бегущих лучей; Sv - суммарная площадь витков [6].
Дрейфы проявляют себя в состоянии покоя ВОГ в виде кратковременного либо значительного по времени ухода полезного сигнала. Проще говоря, ВОГ регистрирует фиктивную угловую скорость. Источником дрейфа может служить множество факторов. Кратковременные уходы чаще всего бывают вызваны ударной либо вибрационной нагрузкой, длительные по времени - изменением теплового состояния оптоволоконного контура. Дрейф также могут порождать различные включения в сердцевине оптического волокна, вызывающие рэлеевское рассеяние, приводящее к непредсказуемому изменению фазы лучей. Радиационное и электромагнитное возмущения также влияют на показания ВОГ, но на сегодняшний день эти проблемы практически решены. С рэлеевским рассеянием борются, совершенствуя процесс изготовления волокна, от воздействия радиации избавляют специальные радиационно-стойкие волокна, а электромагнитные возмущения подавляются специальными экранами, установленными на изделии [6-8].
Актуальной по-прежнему остается проблема термически и механически индуцированных дрейфов [9; 10]. В основе обоих лежит зависимость коэффициента преломления кварцевой жилы оптического волокна от его напряженно-деформированного состояния (НДС).
При нестационарном неравномерном температурном возмущении оптоволоконный контур ввиду высокой жесткости кварца деформируется с крайне интенсивным ростом напряжений (порядка 0,039 МПа/мин), что приводит к значительному по амплитуде дрейфу (от 0,1 до 4 град/ч и выше в зависимости от темпа нагрева). Малые коэффициенты температуропроводности усиливают эффект, растягивая его во времени. В среднем на прогрев оптоволоконного контура от 20 °С до 60 °С и выход его на квазистационарный режим требуется не менее 40 минут при темпе нагрева изделия 1 °С/мин.
Ударные и вибрационные нагрузки не менее опасны. Здесь в первую очередь важна амплитуда дрейфа, достигающая сотен градусов в час, что совершенно неприемлемо.
Существуют как конструктивные, так и программные способы уменьшения подобных погрешностей прибора. Конструктивный путь ведет к оптимизации конструкции, геометрии изделия, подбору материалов с близкими коэффициентами температурного расширения, программные методы борьбы в основном используются для снижения термических дрейфов. Это прежде всего так называемые алгоритмы коррекции показаний гироскопа в виде прогноза и компенсации дрейфа [11; 12], основанные на анализе сигнала ВОГ в реальном времени. Этот способ довольно успешен, но он не снимает саму проблему, а лишь маскирует ее.
Таким образом, целью данной работы является исследование природы возникновения дрейфа ВОГ в условиях внешних тепловых нестационарных воздействий и разработка алгоритмов численного моделирования погрешностей прибора.
1. Математическая постановка
Одним из наименее затратных в финансовом плане путей решения возникшей проблемы является прямое численное моделирование, не привязанное к анализу сигнала. В текущей работе демонстрируется алгоритм получения сигнала ВОГ при нестационар-
57
Galyagin K.S., Oshivalov M.A., Savin M.A. /PNRPUMechanics Bulletin 4 (2015) 55-71
ном термическом возмущении непосредственно по оценке НДС оптического волокна, уложенного на каркасе гироскопа.
Вся процедура расчета и получения дрейфа с целью верификации модели разделена на два основных этапа. Первый этап не столь интересен и подразумевает обычный расчет термо-НДС конструкции ВОГ. Здесь стоит лишь отметить, что решается осесимметричная задача термоупругости (рис. 1), а оптическое волокно представлено не анизотропным усредненным материалом, а дискретизовано с учетом собственной микроструктуры [13]. Это означает наличие подробной конечно-объемной сетки на каждом витке оптического волокна. Сеткой покрываются оба защитно-упрочняющих покрытия, кварцевая оболочка и силовые нагрузочные стержни из боросиликатного стекла (рис. 2). Сам волновод, или светопроводящая жила, не моделируется ввиду ничтожно малого отличия термомеханических свойств от свойств кварцевой оболочки.
2
1
4 *
Рис. 1. Компьютерная модель: 1 - электромагнитный экран; 2 - воздух; 3 - корпус;
4 - оптическое волокно
Рис. 2. Сеточная модель оптического волокна: 1 - компаунд; 2 - вторичное защитно-упрочняющее покрытие; 3 - первичное защитно-упрочняющее покрытие; 4 - кварцевая оболочка; 5 - боросиликатные силовые стержни
Построение сетки практически полностью автоматизировано как на геометрии корпуса ВОГ, так и на геометрии оптического волокна. Расчет также автоматизирован и проходит под управлением специально разработанных макросов, своевременно экспортирующих необходимые данные в табличные файлы.
58
Галягин К.С., Ошивалов М.А., Савин М.А. /Вестник ПНИПУ. Механика 4 (2015) 55-71
Далее происходит обработка полученных файлов в среде MATLAB по специальному алгоритму. Перед непосредственным рассмотрением алгоритма стоит более подробно пояснить причины возникновения термически индуцированного дрейфа.
В оптоволоконном контуре ВОГ всегда существуют механические напряжения. Их наличие обусловлено как конструкцией чувствительного элемента, так и технологией его изготовления. В частности, напряжения в оптическом волокне возникают при изгибе и натяжении во время намотки, а также после полимеризации компаундной заливки. В большинстве современных гироскопов используется волокно с устойчивой поляризацией. Это делается для увеличения двулучепреломления и как следствие обеспечения работы световода на данной собственной поляризационной моде. Такие волокна имеют асимметричный профиль показателя преломления, который достигается созданием в области сердечника волокна поля остаточных напряжений соответствующей структуры [14, 15]. При суперпозиции их с технологическими напряжениями образуется результирующее НДС оптоволоконного контура.
В условиях внешних тепловых воздействий на оптоволоконный контур изменение температуры волокна приводит к появлению термонапряжений, изменяющих величину показателя преломления. В этом случае его следует считать функцией двух параметров -температуры и напряжений (или деформаций) n = n(t,c). Изменение показателя преломления можно представить в виде
7 dn 7 дп 7
dn ——dt н--d о —
dt до
гдп дп дол к дt до дt у
dt.
(1.1)
Наличие механических напряжений в волокне может изменить первоначальное состояние поляризации излучения. В световодах с устойчивой поляризацией две ортогональные компоненты фундаментальной моды имеют разные фазовые постоянные распространения. Это уменьшает связь по мощности между двумя поляризациями и, следовательно, связь между поляризационными модами. В таких световодах двулучепреломление создано, как правило, за счет анизотропии остаточных напряжений. Мерой этой анизотропии является модальное двулучепреломление [3]:
в — (,-Р,)/(2Я / X), (1.2)
где рх, РУ - фазовые постоянные распространения для двух направлений поляризации.
Результирующее НДС в области сердечника волокна может изменить первоначальную ориентацию площадок главных напряжений и связанную с ними ориентацию главных осей тензора диэлектрической проницаемости, в результате чего линейно поляризованное излучение становится эллиптически поляризованным и усиливается связь между поляризационными модами.
Внешние тепловые воздействия приводят к температурным деформациям каркаса и залитой компаундом катушки оптоволоконного контура. Деформации вызывают изменение площади витков, а так как от этого параметра зависит фаза Саньяка, то возникает сигнал ошибки. Экспериментальные исследования ВОГ показали, что для минимизации ошибки нужно согласовывать коэффициенты теплового расширения материала каркаса и компаунда. Тогда деформация витков контура становится симметричной относительно середины волокна, и индуцированная напряжениями невзаимность уменьшается.
59
Galyagin K.S., Oshivalov M.A., Savin M.A. /PNRPUMechanics Bulletin 4 (2015) 55-71
Практически у всех волоконно-оптических гироскопов наблюдается дрейф сигнала в течение довольно продолжительного времени после включения. Это вызвано тем, что источник излучения размещается, как правило, вблизи оптоволоконного контура и его начальный разогрев до стационарного рабочего режима аналогичен внешнему тепловому воздействию на контур. Для оценки влияния пьезооптических эффектов на дрейф ВОГ нужно установить связь между возникающим в волокне напряженно-деформированным состоянием и его поляризационно-оптическими характеристиками.
В основе метода фотоупругости лежит явление искусственного двойного лучепреломления, связанное с тем, что прозрачные материалы под воздействием напряжений приобретают оптическую анизотропию. Степень этой анизотропии зависит от НДС материала и характеризуется тензором диэлектрической проницаемости. В связи с постановкой задачи фотоупругости для оптического волокна рассмотрим некоторые основные положения метода.
Напряженное состояние в точке определяют девять компонент напряжений, которые являются составляющими тензора второго ранга
T =
а т т
x xy xz
т а т
yx y yz
т т а
zx zy z
(1.3)
Здесь а - нормальные, т - касательные напряжения. Эти напряжения на каждой грани мысленно вырезанного параллелепипеда представляют собой составляющие вектора полного напряжения на соответствующей площадке. Для касательных напряжений выполняется закон парности:
тyx Txy, Tzx Txz, Tyz Tzy, (1.4)
и, следовательно, только 6 компонент тензора напряжений независимы [16].
Напряжения на площадках, произвольно ориентированных относительно осей координат, определяются по формулам преобразования компонент тензора второго ранга при повороте координатных осей.
При любом виде напряженного состояния через каждую точку можно провести три взаимно перпендикулярные главные площадки. Нормальные напряжения на этих площадках называются главными, а касательные напряжения отсутствуют. Главные напряжения обозначаются через а1, а2, а3 и нумеруются так, что а1 > а2 > а3.
В фотоупругости важную роль играют так называемые квазиглавные напряжения. Пусть система координат x1, y1, z1 повернута на угол а вокруг оси z1, совпадающей с осью z системы x, y, z (рис. 3). Запишем формулы для напряжений по двум взаимно перпендикулярным площадкам BC и AB, которые параллельны оси z (перпендикулярной плоскости рисунка) и внешние нормали к которым параллельны осям x1 и y1. Для площадки BC
аx1 = аx cos а + аy sin а + тxy sin 2а,
< тx1 y1 = - СТx 2 СТy sin 2а + тyx cos 2а, (1.5)
тx1z =Tzy sin а + тzx cos а-
60
Галягин К.С., Ошивалов М.А., Савин М.А. /Вестник ПНИПУ. Механика 4 (2015) 55-71
Рис. 3. Напряжения по площадкам, параллельным оси z
Для площадки AB
ау1 = аx sin2 а + ау cos2 а - тxy sin 2а,
| Тylxl = Тх1 у1,
Т y1z = Т zy C0s а-т zx sin а.
(1.6)
Выберем а таким, чтобы напряжения а и а принимали экстремальное значение:
d а
xl
= 0,
d а
yi
= 0
d а d а
( -ау )п2аа + 2тху С0Б2аа = 0.
(1.7)
Если угол а = а0, то на площадках AB и BC отсутствуют касательные напряжения тху. Эти площадки называются квазиглавными. В отличие от главных площадок на квазиглавных площадках касательные напряжения Txz в общем случае не равны нулю. Нормальные напряжения на этих площадках а^) и а2z) называются квазиглавными для выбранного направления z. Угол а0 измеряется относительно а(^ (большее напряжение) и осью x.
Выбирая оси х1 и у1 так, чтобы они совпали с направлениями квазиглавных направлений, из (1.6) и (1.7) получим (полагая а = -аД
(z) 2
аx = а1 ’ C0s аа ,
(z) • 2 . (z) 2
аy = а1^т аа+аv2;cos аа,
T(z ^-а()
т = •
xy
2
^т2аа.
(1.8)
Разрешая (1.8) относительно квазиглавых напряжений, получим
61
Galyagin K.S., Oshivalov M.A., Savin M.A. /PNRPUMechanics Bulletin 4 (2015) 55-71
(z) _ ax + ay
a2 _
2 -y
+ a y _
2 д
2
_°y '
+ Txy 2,
x _ay л
+ T
xy ■
(1.9)
tg «a
(1.10)
Выражения (1.9) и (1.10) определяют квазиглавные напряжения и угол наклона площадок в системе координат x, y, z для направления z. Для данного направления эта пара квазиглавных напряжений является единственной.
Оптические свойства немагнитных анизотропных сред характеризуются тензором диэлектрической проницаемости. Аналогично тензорам напряжений и деформаций он является симметричным тензором второго ранга:
A
%11 %12 %13 %21 %22 X23
Х31 X32 Хзз
(1.11)
Аналогично квазиглавным напряжениям вводится понятие о квазиглавных компонентах диэлектрического тензора. Квазиглавные диэлектрические проницаемости определяют свойства оптической среды для заданного направления просвечивания (например, z, как для напряжений). Их величины определяются выражениями типа (1.9), а направления - формулами типа (1.10).
В оптически анизотропных средах направления плоскостей поляризации совпадают с квазиглавными направлениями тензора диэлектрической проницаемости [17, 18]. В частности, для направления просвечивания z скорости и показатели преломления двух лучей, поляризованных в квазиглавных направлениях, определяются следующим образом [19]:
u1 _ c / n1; u2 c / n ;
(1.12)
где с - скорость света.
Механически изотропная среда, находящаяся в напряженно-деформированном состоянии, приобретает оптическую анизотропию. Поэтому плоскополяризованные в квазиглавных направлениях диэлектрического тензора волны при прохождении через такую среду получают сдвиг фаз или оптическую разность хода, величина которой зависит от НДС среды.
Основы теории пьезооптического эффекта для упругих тел, подчиняющихся законам Гука, изложены в работах Неймана и Максвелла. В соответствии с ними главные значения тензора диэлектрической проницаемости линейно связаны с главными деформациями (Нейман) и главными напряжениями (Максвелл). Для напряжений эти зависимости выглядят следующим образом:
62
Галягин К.С., Ошивалов М.А., Савин М.А. /Вестник ПНИПУ. Механика 4 (2015) 55-71
Ь = Ь + c1g1 + С2 ( + g3 ) %2 = Ь + c1g2 + С2 ((3 + g1 ) Ьз = Ь + c1g3 + c2 (1 + g2 )
(1.13)
где ci, c2 - оптические (фотоупругие) постоянные материала для света с колебаниями соответственно параллельными и перпендикулярными действию главных напряжений. Они характеризуют вклад в изменение показателя преломления от соответствующих напряжений; х0 - диэлектрическая проницаемость в ненапряженном состоянии.
Разность фаз двух поляризованных волн при просвечивании напряженного элемента вдоль направления х3
Лф = 2к—,
X
где оптическая разность хода 5 при толщине элемента d
5 = d((-VxT) = d(c/u -c/U2).
Из (1.13) находим
& -&]1+c,g| +c2(2 +^ ■
V Xo
(1.14)
(1.15)
(1.16)
Для всех известных материалов второе слагаемое под радикалом <<1, поэтому
1 + c1g1 + С2 (2 + g3 )
2X 0
Ь *\Д о
(1.17)
Получая аналогичное выражение для ^д2 , из (1.16) и (1.17), находим
'(g1 -g2 )
5=/х0-
2Ьо
= СЛа1 -g2 ),
(1.18)
где
c
С = c1 c2 2n0 2n0
(1.19)
оптический коэффициент напряжений.
Как видно из (1.18), размерность у cg обратна размерности напряжений. Таким образом, по теории Максвелла оптическая разность хода пропорциональна разности главных напряжений, а направления главных диэлектрических проницаемостей соосны с главными напряжениями.
Для оптических волокон с устойчивой поляризацией показатели преломления в квазиглавных направлениях обусловлены полем остаточных напряжений. При изготовлении волоконного контура ВОТ внутри волокна возникают технологические напряжения, а в условиях внешних воздействий (температура, вибрации и т.д.) - индуцированное или «внешнее» поле напряжений. Суперпозиция всех этих полей и создает результирующие напряжения внутри волокна, которые определяют реальное значение показателей в каждый момент времени.
63
Galyagin K.S., Oshivalov M.A., Savin M.A. /PNRPUMechanics Bulletin 4 (2015) 55-71
2. Постпроцессор
В условиях нестационарного теплового режима дрейф определяется изменением показателей преломления волокна за время пробега по контуру встречно бегущих лучей [20]. Для расчета результирующей разности фаз необходимо определить оптические пути обоих лучей либо время, в течение которого они пробегают контур.
После разрешения задачи термоупругости ВОГ в пакете STAR-CCM+ имеется ряд файлов с исходной информацией для постпроцессора, рассчитывающего дрейф. Эта информация представляет собой неупорядоченный массив данных, включающий в себя зна-
чения трех главных напряжений и величину температуры в центрах каждого витка осесимметричной модели. С использованием специальной процедуры «разворачивания» осесимметричные поля искомых функций в соответствии с принятой схемой намотки волокна преобразуются в одномерные функции по длине уложенного в катушку волокна. Алгоритм «разворачивания» сортирует исходный массив данных по координатам витков в последовательности намотки для каждой ветви контура. Отметим, что существует несколько основных схем укладки оптического волокна на каркас ВОГ. Оптимальной с точки зрения минимизации теплового дрейфа является квадрупольная намотка [6, 8] или укладка (рис. 4, а), она же и используется при расчете дрейфа. Существуют также и менее удачные схемы: дипольная, шахматная и др.
а б
Рис. 4. Схемы укладки: а - квадрупольная; б - дипольная. Черная точка -начальный виток середины контура, серые точки - «левая ветвь» контура, белые точки - «правая ветвь» контура
Расчет оптического пути (или времени пробега) для каждого луча иллюстрирует рис. 5.
(2.1)
ni = %-1 +
С1°1 + С2 К +а3 )
n,
0i
где n0i - показатель преломления при начальной температуре; с1, с2 - фотоупругие постоянные материала.
По волокну в направлении х бежит луч со скоростью w. В начале и в конце каждого витка длиной l по формулам типа (2.1) рассчитываются показатели преломления
n^1, nr1, Пк, n£. Здесь верхние индексы указывают на принадлежность соответствую-
64
Галягин К.С., Ошивалов М.А., Савин М.А. /Вестник ПНИПУ. Механика 4 (2015) 55-71
щему временному слою, а нижние - началу (н) и концу (к) витка. В формуле (2.1) n0i является линейной функцией температуры и вычисляется в зависимости от температуры в данной точке и в данный момент времени.
При расчетах предполагается, что шаг интегрирования hT много больше времени пробега луча по контуру. Для произвольного витка время входа луча в виток соответствует на рисунке точке а, а выхода - точке b. Время пробега витка лучом равно Ат. Для его определения нужно проинтегрировать уравнение
dx _ с d х n( x, x) ’
(2.2)
где с - скорость света в вакууме.
Численное интегрирование уравнения выполняется в следующем порядке. Определяются скорости изменения показателей преломления в начале и в конце витка:
dnH
dnK
n„ - n„
„к-1
n„ - n„
„к-1
h
h
(2.3)
Виток разбивается на m участков длиной H = l/m, и для центра хт каждого участка вычисляется значение показателя преломления:
nm _
( к к Л /
к Пк - Пн 1
n к +—--- x
н
V
l
т
+
, dnк - d^ ^
dn„ +—к-------- x„
У
V
l
Х т ,
(2.4)
где тт - время входа луча в данный участок.
Это время для каждого последующего участка вычисляется по формуле
х _ х +
т т
nm • H /с
1 - 0,5 [ dn^ + xm (dnK - dn^) /l ]
(2.5)
После выхода из процедуры интегрирования время тт равно времени прохождения лучом данного витка Ат. Суммируя эти величины для всех витков, находим общее время
65
Galyagin K.S., Oshivalov M.A., Savin M.A. /PNRPUMechanics Bulletin 4 (2015) 55-71
пробега лучом контура ть Аналогично вычисляется время пробега встречного луча т2. Результирующую разность фаз встречно бегущих лучей можно определить следующим образом:
ЛФ = ^(Ат1; - Ат2; ), (2.6)
где N - число витков в контуре; Лть, Лт2 - время пробега i-го витка контура лучами; X -длина волны излучения.
Далее по найденной фазовой невзаимности рассчитывается «кажущаяся» угловая скорость вращения гироскопа в соответствии с известной формулой Саньяка (В1).
В настоящее время соотношение для вычисления показателя преломления витка упрощено ввиду осесимметричности задачи и неоправданности в таком случае деления витка на участки, так как нет информации об НДС в промежуточных точках. Однако эта формула потенциально применима при решении задачи о дрейфе в трехмерной постановке.
На рис. 6 можно видеть результат так называемой размотки оптического волокна в постпроцессоре. Здесь демонстрируется распределение температуры по общей длине оптического волокна. Центральные витки расположены на середине общей длины оптического волокна, от середины в обе стороны расходятся «левая» и «правая» ветви. Пульсации соответствуют переходу со слоя на слой либо через два слоя, гладкие впадины соответствуют распределению температуры по слою. Прослеживается явная асимметрия нагрева относительно центральных витков катушки, расположенных на середине общей длины оптического волокна. Подобные распределения относительно длины оптоволокна можно получить и для главных напряжений (например, рис. 7 демонстрирует распределение первого главного напряжения).
Рис. 6. Распределение температуры по длине волокна
66
Галягин К.С., Ошивалов М.А., Савин М.А. /Вестник ПНИПУ. Механика 4 (2015) 55-71
Основной задачей является получение фазовой картины и дрейфа в целом. Так, на рис. 8 изображено распределение времени пробега одного из лучей по каждому витку в зависимости от длины волокна с учетом квадрупольной схемы намотки. Отчетливо видны переходы со слоя на слой - малые скачки соответствуют переходу на соседний слой, большие скачки - переходу через два слоя. Ожидаемо и снижение времени хода луча при приближении к центру оптического волокна - центральные витки расположены непосредственно на катушке ВОТ.
При нахождении разности времен хода получена фазовая картина, изображенная на рис. 9. Симметричное изменение фазового набега из отрицательной области в положительную соответствует обгону встречно распространяющимся лучом нормального луча на первой половине оптического пути, а затем отставанию встречного луча на второй половине пути.
67
Galyagin K.S., Oshivalov M.A., Savin M.A. /PNRPUMechanics Bulletin 4 (2015) 55-71
Анализ полученной разности фаз с помощью соотношений Саньяка (В1) с учетом перевода единиц измерения приводит к получению расчетного дрейфа ВОГ (рис. 10).
Разработанная компьютерная модель позволяет оценить на качественном уровне величину дрейфа ВОГ при заданной геометрии и граничных условиях. Так как задача термоупругости разрешается в специализированном инженерном пакете, то возможно широкое варьирование граничных и начальных условий. Возможно также любое варьирование геометрических размеров конструкции и числа витков/слоев оптического волокна в волоконном контуре. Это позволяет оценить качество геометрии как катушки, так и датчика чувствительности ВОГ в целом для каждой конкретной конструкции.
Апробация разработанной модели для различных расчетных вариантов показала возможность адекватного прогноза дрейфа, вызванного изменением НДС оптоволоконного контура и конструкции датчика чувствительности ВОГ, при условии заведомо известных достоверных свойств всех материалов.
68
Галягин К.С., Ошивалов М.А., Савин М.А. /Вестник ПНИПУ. Механика 4 (2015) 55-71
Библиографический список
1. Жижин В. Волоконно-оптические датчики: перспективы промышленного применения // Электронные компоненты. - 2010. - № 12. - С. 17-23.
2. Малыкин Г.Б., Андронова И.А. Физические проблемы волоконной гироскопии на эффекте Саньяка // Успехи физических наук. - 2002. - Т. 172, № 8. - С. 849-873.
3. Логозинский В.Н., Листвин В.Н. Волоконно-оптический датчик вращения [Электронный ресурс] // Физоптика, 2011. - URL: www.fizoptika.ru/describtion/book_bind.pdf (дата обращения: 18.11.1014)
4. Czommer R. Leistungsfahigkeit fahrzeugautonomer Ortungsverfahren auf der Basis von Map-Matching-Techniken. Von der Fakultat fur Bauingenieur und Vermessungswesen der Universitat Stuttgart zur Erlangung der Wurde eines Doktors-Ingenieurs (Dr.-Ing.) Genehmigte Dissertation. - Stuttgart, 2000.
5. Neuhierl T. Eine neue Methode zur Richtungsubertragung durch Koppelung von Inertialmesstechnik und Autokollimation. Vollstandiger Abdruck der von der Fakultat fur Bauingenieur-und Vermessungswesen der Technischen Universitat Munchen zur Erlangung des akademischen Grades eines Doktor-Ingenieurs (Dr.-Ing.) Genehmigten Dissertation. - Munchen, 2005.
6. Шереметьев А.Г. Волоконный оптический гироскоп. - М.: Радио и связь, 1987. - 152 с.
7. Окоси Т. Волоконно-оптические датчики. - Л.: Энергоатомиздат, 1990. - 256 с.
8. Lefevre, Herve Fiber-optic Gyroscopes. - ARTECH HOUSE, INC., 1993.
9. Тепловой дрейф волоконного оптического гироскопа / Е.И. Вахрамеев, К.С. Галягин, А.С. Ивонин, М.А. Ошивалов, Т.А. Ульрих // Изв. вузов. Приборостроение. - 2011. - Т. 54, № 1. -С.32-37.
10. Рупасов А.В. Исследование метода локального температурного воздействия и его применение для компенсации дрейфа волоконно-оптического гироскопа: дис. ... канд. физ.-мат. наук. -СПб., 2014.
11. Шарков И.А. Исследование и пути компенсации тепловых воздействий на сигнал волоконно-оптического гироскопа: дис. ... канд. техн. наук. - СПб., 2013.
12. Прогноз и коррекция теплового дрейфа волоконно-оптического гироскопа / Е.И. Вахрамеев, К.С. Галягин, А.С. Ивонин, М.А. Ошивалов // Изв. вузов. Приборостроение. - 2013. - Т. 56, № 5. - С. 79-84.
13. Савин М.А. Дрейф волоконно-оптического гироскопа // Современные техника и технологии: сб. докл. XX Междунар. науч.-практ. конф. студентов, аспирантов и молодых ученых / НИТ-ПУ. - Томск, 2014. - Т. 1. - С. 133-134.
14. Филатов Ю.В. Волоконно-оптический гироскоп: учеб. пособие. - СПб.: Изд-во С.Петерб. гос. электротехн. ун-та «ЛЭТИ», 2003. - 52 с.
15. Ebeling K.J., Michalzik R., MahnB J. Optische Informationstechnik. - Universitat Ulm: Abtei-lung Optoelektronik, 2009. - 191 s.
16. Лехницкий С.Г. Теория упругости анизотропного тела. - М.: Наука, 1977. - 416 c.
17. Александров А.Я., Ахметзянов М.Х. Поляризационно-оптические методы механики деформируемого тела. - М.: Наука, 1973. - 576 c.
18. Сивухин Д.В. Общий курс физики. Оптика. - М.: Физматлит, 2005. - 792 c.
19. Смоленский Г.А., Писарев Р.В., Синий И.Г. Двойное лучепреломление света в магнитоупорядоченных кристаллах // Успехи физических наук. - 1975. - Т. 116, № 2. - С. 231-270.
20. Драницына Е.В., Егоров Д.А. Исследование зависимости выходного сигнала волоконнооптического гироскопа от температуры в составе бескарданного инерциального измерительного модуля // Навигация и управление движением: материалы XIV конференции молодых ученых. -СПб., 2012. - С. 447-452.
69
Galyagin K.S., Oshivalov M.A., Savin M.A. /PNRPUMechanics Bulletin 4 (2015) 55-71
References
1. Zhizhin V. Volokonno-opticheskie datchiki: perspektivy promyshlennogo primeneniia [Fiber Optic Sensors: Prospects for industrial applications]. Elektronnye komponenty, 2010, no. 12, pp. 17-23.
2. Malykin G.B., Andronova I.A. Fizicheskie problemy volokonnoi giroskopii na effekte San’iaka [Physical problems of fiber gyroscope on the Sagnac effect]. Uspekhi fizicheskikh nauk, 2002, vol. 172, no. 8, pp. 849-873.
3. Logozinskii V.N., Listvin V.N. Volokonno-opticheskii datchik vrashcheniia [The fiber optic rotation sensor]. Fizoptika, 2011, available at: www.fizoptika.ru/describtion/book_bind.pdf, svobodnyi. Iaz. rus. (accessed 18 November 1014).
4. Czommer R. Leistungsfahigkeit fahrzeugautonomer Ortungsverfahren auf der Basis von Map-Matching-Techniken. Von der Fakultat fur Bauingenieur und Vermessungswesen der Universitat Stuttgart zur Erlangung der Wurde eines Doktors-Ingenieurs (Dr.-Ing.) Genehmigte Dissertation. Stuttgart, 2000.
5. Neuhierl T. Eine neue Methode zur Richtungsubertragung durch Koppelung von Inertialmesstechnik und Autokollimation. Vollstandiger Abdruck der von der Fakultat fur Bauingenieur-und Vermessungswesen der Technischen Universitat Munchen zur Erlangung des akademischen Grades eines Doktor-Ingenieurs (Dr.-Ing.). Genehmigten Dissertation. Munchen, 2005.
6. Sheremet’ev A.G. Volokonnyi opticheskii giroskop [The fiber optic gyroscope]. Moscow: Radio i sviaz’, 1987. 152 p.
7. Okosi T. Volokonno-opticheskie datchiki [Fiber Optic Sensors]. Leningrad: Energoatomizdat, 1990. 256 p.
8. Lefevre, Herve Fiber-optic Gyroscopes. ARTECH HOUSE, INC., 1993.
9. Vakhrameev E.I., Galiagin K.S., Ivonin A.S., Oshivalov M.A., Ul’rikh T.A. Teplovoi dreif volokonnogo opticheskogo giroskopa [The thermal drift of the fiber optic gyroscope]. Izvestiia vuzov. Priborostroenie, 2011, vol. 54, no. 1, pp. 32-37.
10. Rupasov A.V. Issledovanie metoda lokal’nogo temperaturnogo vozdeistviia i ego primenenie dlia kompensatsii dreifa volokonno-opticheskogo giroskopa [The research method of local temperature effects and its application for compensation of the drift of the fiber optic gyroscope]. Thesis of doctor’s degree dissertation. Saint Petersburg, 2014.
11. Sharkov I.A. Issledovanie i puti kompensatsii teplovykh vozdeistvii na signal volokonno-opticheskogo giroskopa [Research and ways of compensating thermal effects on the signal of the fiber optic gyroscope]. Thesis of doctor’s degree dissertation. Saint Petersburg, 2013.
12. Vakhrameev E.I., Galiagin K.S., Ivonin A.S., Oshivalov M.A. Prognoz i korrektsiia teplovogo dreifa volokonno-opticheskogo giroskopa [Prediction and correction of thermal drift of the fiber optic gyroscope]. Izvestiia vuzov. Priborostroenie, 2013, vol. 56, no. 5, pp. 79-84.
13. Savin M.A. Dreif volokonno-opticheskogo giroskopa [The drift of the fiber optic gyroscope]. Sbornik dokladov XX mezhdunarodnoi nauchno-prakticheskoi konferentsii studentov, aspirantov i molodykh uchenykh ’’Sovremennye tekhnika i tekhnologii”. Tomsk, 2014, vol. 1, pp. 133-134.
14. Filatov Iu.V. Volokonno-opticheskii giroskop [Fiber Optic Gyroscope]. S. Peterburgskii gosudarstvennyi elektrotekhnicheskii universitet, 2003. 52 p.
15. Ebeling K.J., Michalzik R., MahnB J. Optische Informationstechnik. Universitat Ulm: Abteilung Optoelektronik, 2009. 191 p.
16. Lekhnitskii S.G. Teoriia uprugosti anizotropnogo tela [The theory of elasticity of an anisotropic body]. Moscow: Nauka, 1977. 416 p.
17. Aleksandrov A.Ia., Akhmetzianov M.Kh. Poliarizatsionno-opticheskie metody mekhaniki deformiruemogo tela [The polarization-optical methods for solid mechanics]. Moscow: Nauka, 1973. 576 p.
70
Галягин К.С., Ошивалов М.А., Савин М.А. /Вестник ПНИПУ. Механика 4 (2015) 55-71
18. Sivukhin D.V. Obshchii kurs fiziki. Optika [The general course of physics. Optics]. Moscow: Fizmatlit, 2005. 792 p.
19. Smolenskii G.A., Pisarev R.V., Sinii I.G. Dvoinoe lucheprelomlenie sveta v magnitouporiado-chennykh kristallakh [The birefringence of light in magnetically ordered crystals]. Uspekhi fizicheskikh nauk, 1975, vol. 116, no. 2, pp. 231-270.
20. Dranitsyna E.V., Egorov D.A. Issledovanie zavisimosti vykhodnogo signala volokonno-opticheskogo giroskopa ot temperatury v sostave beskardannogo inertsial’nogo izmeritel’nogo modulia [The dependence of the output signal of the fiber optic gyroscope with temperature as part of strapdown inertial measurement unit]. Materialy dokladov XIV konferentsii molodykh uchenykh ’’Navigatsiia i upravlenie dvizheniem”. Saint Petersburg, 2012, pp. 447-452.
71
