Научная статья на тему 'Обработка сигналов импульсного волоконно-оптического датчика ускорения'

Обработка сигналов импульсного волоконно-оптического датчика ускорения Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
175
27
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
FIBER-OPTIC ACCELERATION SENSORS / PULSE MODULATION / FIBER OPTIC ACCELEROMETERS / ВОЛОКОННО-ОПТИЧЕСКИЕ ДАТЧИКИ УСКОРЕНИЯ / ИМПУЛЬСНАЯ МОДУЛЯЦИЯ / ВОЛОКОННО-ОПТИЧЕСКИЕ АКСЕЛЕРОМЕТРЫ / ВОЛОКОВОЛОКОННО-ОПТИЧНі ДАВАЧі ПРИСКОРЕННЯ / іМПУЛЬСНА МОДУЛЯЦіЯ / ВОЛОКОННО-ОПТИЧНі АКСЕЛЕРОМЕТРИ

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Демьяненко П.А., Зиньковский Ю.Ф.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Демьяненко П.А., Зиньковский Ю.Ф.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Signal processing of the optical fiber impulse sensor acceleration

Introduction. Feature of the pulsed fiber-optic sensor is a mutual influence of the components of the acceleration vector measured. When processing the signals of the pulse sensor should provide adequate correction value counts for correct measurement results. Main part. It consists of four sections. The first section outlines the method of determining the undistorted values of horizontal components of the acceleration vector acting on the sensor, in the absence of the vertical component of the acceleration vector. The second section outlines the method of determining the undistorted value of vertical component of the acceleration vector. Based on the results obtained in the first two chapters, the third chapter described a method of determining undistorted measured quantities of all components of the acceleration vector. In the fourth chapter provides a quantitative assessment of the effectiveness of the proposed information processing in the case of precision measurements of acceleration. Conclusions. It is pointed to the principal possibility of determining the undistorted values of samples of the vector acceleration components along all axes of sensitivity of the pulse sensor. The method of determination is proposed and a general algorithm is formulated for measuring randomly oriented components of the acceleration vector. It is evaluated the efficiency of correction of pulse sensor signals for example acceler-ometer of micro-accelerations.

Текст научной работы на тему «Обработка сигналов импульсного волоконно-оптического датчика ускорения»

Обчислювальні методи в радіоелектроніці

УДК 681.586.5; 531.768.

ОБРОБКА ІНФОРМАЦІЇ ІМПУЛЬСНОГО ВОЛОКОННО-ОПТИЧНОГО ДАВАЧА ПРИСКОРЕННЯ1

Дем’яненко П. О., к.т.н., доцент; Зіньковський Ю. Ф., д.т.н., професор

Національний технічний університет України «Київський політехнічний інститут», м.Київ, Україна pdemianenko@gmail.com

SIGNAL PROCESSING OF THE OPTICAL FIBER IMPULSE SENSOR

ACCELERATION

Dem’yanenko P. O., PhD, Associate Professor;

Zinkovsky Yu. F., Doc.of Sci(Technics), Professor

National Technical University of Ukraine «Kyiv polytechnic Institute», Kyiv, Ukraine

Вступ

Використання в сучасній техніці традиційних електричних давачів часто наштовхується на важкорозв’язувані проблеми захисту як їх самих, так і трактів каналювання їхніх сигналів від різних електромагнітних (ЕМ) завад. В цьому сенсі привабливими є волоконно-оптичні давачі (ВОД) з принциповою невразливістю їхніх сигналів збоку будь-яких ЕМ завад.

Нажаль, як показала практика, на основі звичайних аналогових ВОД неможливо побудувати прецизійних вимірювачів фізичних величин, зокрема, акселерометрів, які могли б скласти гідну конкуренцію електричним вимірювачам. Так, наприклад, описаний в [1] ВОД прискорення, який модулює інтенсивність оптичного потоку, має порогову чутливість ~ 2-10-6 g. Для ВОД прискорення, в основі роботи якого лежить ефект фотопружності [2], експериментально визначені мінімальні значення порогового приско-

о

рення склали ~ 8,5-10- g.

Рекордну чутливість, що склала 10-9 g, отримали автори [3] при вимірюванні прискорення сили тяжіння, опосередкувавши результати вимірювань в інтервалі часу 100 с. Ці дані були реалізовані завдяки використанню мобільного атомного гравіметра, який розміщувався в підземній дослідницькій лабораторії з низьким рівнем шумів на глибині 500 м. Пробною (чутливою) масою гравіметра була хмаринка охолоджених до температури 2-10- K атомів Rb, які під дією сили тяжіння падали в вакуумі в тривимірну магнітно-оптичну яму, що створювалася після адіабатичного відключення лазерних променів.

Зауважимо, що всі наведені приклади стосуються лабораторних макетів вимірювачів, яким до реальних вимірювальних інструментів необхідно по-

1 http://radap .kpi.ua/radiotechnique/article/view/1075

Вісник Національного технічного університету України «КПІ»

Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №63 33

Обчислювальні методи в радіоелектроніці

долати ще не простий шлях.

Наведені вище цифри є, на перший погляд, вельми пристойними. Одначе, для розв‘язку низки нагальних задач геофізики та навігації порогову чутливість акселерометрів бажано мати від 10" g і менше. Серед таких задач можна відзначити:

• експериментальна перевірка стабільності положення центра мас Землі;

• вимірювання параметрів траєкторії зміщення полюсів Землі;

• вимірювання гравітаційної сталої;

• вивчення неоднорідностей гравітаційного поля Землі та інших космічних об‘єктів безпосередньо з орбіт їх штучних супутників за допомогою малорозмірних градієнтометрів;

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

З доступних сьогодні на ринку промислових типів вимірювачів можна відзначити, наприклад, п‘єзоелектричні акселерометри фірми ZETLAB. Серед них, наприклад, акселерометр типу АР 2050, порогова чутливість якого складає ~ 2-10"5 g [4].

Причиною незадовільних показників ВОД щодо точності є метрологічно низька якість потоку оптичного випромінювання у вимірювальному тракті ВОД, який одночасно є і носієм, і реципієнтом інформації [5]. Разом з тим, згадана вище перевага ВОД є настільки вагомою, що спонукає шукати шляхи виходу із ситуації, що склалася.

Її аналіз привів нас до усвідомлення необхідності відмови від аналогових і переходу до неаналогових (дискретних) принципів модуляції в ВОД. Реально це означає, що під час модуляції оптичного потоку в нього слід вводити нові, додаткові до оптичних але неоптичні параметри, які і будуть відтепер виконувати роль реципієнтів інформації. Такі ВОД зберігають всі переваги волоконної оптики, оскільки носієм інформації залишається оптичний потік. Разом з тим, проблема точності вимірювань переноситься із області принципово низькоточних вимірювань інтенсивності малопотужних оптичних потоків в інші, неоптичні, області, де цієї проблеми або просто немає, або вона там розв’язана на належному рівні.

Такий підхід дозволив розробити новий клас ВОД з імпульсною модуляцією інтенсивності оптичного потоку (ІВОД). Прикладом може бути да-вач лінійних прискорень, на основі якого можуть бути побудовані прецизійні акселерометри (гравіметри). Схематична будова трикоординатного ІВОД прискорення наведена на рис.1 [2]. Чутливим елементом (ЧЕ) тут є конічний маятник 1, утворений консольно закріпленим відрізком кварцового волоконного світловода (ВС), на вільному кінці якого розміщена інерційна маса (ІМ) 2. Незатухаючий рух ЧЕ підтримується почерговою сило-

34

Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №63

Обчислювальні методи в радіоелектроніці

вою дією електромагнітів 4. Розділення вхідного немодульованого і вихід

ного модульованого оптичних потоків відбувається у воло-конно-оп-тичному розгалужу-вачі-сума-торі 5. За відсутності прискорення, вісь обертання конічного маятника співпадає з його вертикальною віссю чутливості OZ, а його горизонтальні осі чутливості ОХ i OY, визначаються осями взаємно-перпендикулярних циліндричних дзеркал 3, які за даної умови є і осями симетрії кола обігу маятника. Очевидно, за таких умов, інтервали часу перебування маятника обабіч кожної з горизонтальних осей будуть однаковими.

Особливістю таких ІВОД є вплив вертикальної складової вектора прискорення на вихідний сигнал ІВОД для горизонтальних складових цього ж вектора. Для отримання неспотворених результатів вимірювання вихідні сигнали ІВОД мають бути відповідним чином відкориговані.

Обробка інформації ІВОД при дії горизонтально орієнтованого вектора прискорення

Хай вектор прискорення а , що діє на ІВОД, співпадає з напрямком

його горизонтальної осі чутливості OY, як це показано на рис. 2, де зображено лише коло обігу, яке описується кінцем конічного маятника. Дія прискорення зсовує це коло на величину у в напрямку, протилежному напрямку дії прискорення. На рис.2 ця ситуація відображена зсувом не самого кола обігу, а адекватним йому зсувом відповідної осі чутливості: ОХ ^ ОХ. Такий зсув, очевидно, порушує симетрію інтервалів часу перебування маятника обабіч осі ОХ. Різницю цих інтервалів позначимо Аг і приймемо її

в якості міри величини прискорення а . Щоб отримати аналітичний вираз ау = f (Ат, ), пов’яжемо величини ау і у рівнянням руху ІМ конічного маятника, скориставшись для цього законами Ньютона і Гука:

тау =- Ay, (1)

або в канонічному вигляді:

у + со20у = 0.

Рис.1. Схематична будова ІВОД прискорення.

Вісник Національного технічного університету України «КПІ»

Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №63 35

Обчислювальні методи в радіоелектроніці

В цих рівняннях: m - інерційна

. 3EI j,

маса на кінці маятника; A = ——, E ,

3 ’

I =

nd4 ~64

L, d

L

відповідно: жорст-

кість, модуль Юнга, полярний момент інерції поперечного перерізу, довжина та діаметр консолі кварцового ВС;

a ,

' колова частота

Ю0 =Л — m

mL

обертального руху конічного маятника.

Користуючись рис.2, пов’яжемо величину зсуву кола обігу конічного

Рис. 2. Коло обігу маятника ІВОД. маятника-у з різницею інтервалів часу Аг перебування ІМ обабіч осі О'Х.: Аху = ij -т2 = l±- R[(л + 2ф)-(л-2ф)] =

V V

. R 4 .у 2T .у

= 4—— ф = — arcsin—= —- arcsin

(2)

ю0 R ю

R0

п

R0

де l, т і 12, т2 - частина довжини кола обігу та час перебування ІМ над віс-

лінійна тттвид-

сю О'Х і під нею, відповідно; v = ю0R, Т0 = — = 2л,.

юо V

mil

3 EI

кість та період колового руху ІМ; R, ф = arcsin— - радіус та кутова міра

R0

зсуву кола обігу ІМ.

Знак різниці Ат визначає напрямок складової вектора прискорення a : якщо Ат^ > 0, то напрямок a збігається з напрямком осі OY.

Із виразу (1):

m mid

У = ay =--------

A у

T

-L (Л

a=^т=2 ay

3EI y 4л2

Нехтуючи знаком «-» та підставивши цей вираз у (2), отримаємо:

2T

Дту = ^

л

arcsin

T2

-L с\

a

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

4л2 R

(3)

а розв’язавши (3) відносно a , маємо:

36

Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №63

Обчислювальні методи в радіоелектроніці

а =

У

ґ2кЛ 2 V T0 У

ґ

R • sin

п • Дт,

2T

V 210 У

(4)

Очевидно,

а..

0 при Аг = 0;

а

(ау)тах при (Ату)тах Т0

(для забезпечення виконання умови таутохронності руху маятника приймаємо R0 = 2 мм).

З урахуванням цього із (4) отримаємо:

(а )

\ У )п

^2 к''

V То У

Ro

3 EI R

(4a)

При визначенні порогової чутливості акселерометра на основі такого ІВОД. будемо виходити з того, що вимірювання часових інтервалів Т і т2 здійснюється шляхом підрахунку чисел імпульсів високочастотного (/твЧ = 100 МГц) генератора щ i n2, якими, відповідно, заповнюються ці інтервали, тобто: т = ПТГВЧ, а т2 = п1ТГВЧ. Таким чином, різниця тривалостей вимірюваних інтервалів визначиться як

Ап

Ату = АпТГВЧ =

f

ГВЧ

де Аn = n - n. Очевидно, умовою вимірювання порогової величини прискорення буде: А n = 1.

Як і раніше, будемо вважати, що вимірювання періодів Т0 і Tz здійснюється шляхом підрахунку числа імпульсів високочастотного генератора m0 i mz, якими, відповідно, заповнюються ці часові інтервали, тобто:

То = тоТГВЧ > а Tz

m Т

m2T ГВЧ .

Тоді величина різниці тривалостей вимірюваних інтервалів буде:

Ат

Ат z АтТГВЧ

f

ГВЧ

де Ат = т - т7

Враховуючи, що (Ату )пор = Тгеч << То, а отже, із виразу (4) отримаємо:

Sin

п (Ат у Ір

2Та

П (АТ У Гр

2Та

(ау )пор =

2 n3R Т3 f

Т 0 J ГВЧ

(4б)

Динамічний діапазон прискорень, які можуть бути виміряними вздовж

2

Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №63

37

Обчислювальні методи в радіоелектроніці

осі OY за допомогою ІВОД, визначиться як:

(ay ).

N..

y /max

(ay)

у У пор

2 2 Т

2т f — 2 1 о

1 0J ГВЧ ■ гті

п п Т

ГВЧ

Таким чином, при дії на ІВОД прискорення, вектор якого співпадає з напрямком горизонтальної осі OY, розрахунок значення прискорення ау

потребує: вимірювання тривалостей інтервалів часу перебування маятника обабіч осі ОХ, відповідно, т і т2; обрахування різниці часових інтервалів Лт = т - т2; обрахування періоду обертання маятника Т0, як суми цих інтервалів: T = ^ + т2; та згідно з (4) - знаходження значення а .

Все мовлене стосовно складової прискорення a повною мірою стосується і складової прискорення a ;

У випадку довільної орієнтації вектора прискорення в горизонтальній площині його модуль визначиться як геометрична сума складових а і ах:

І 2 2"

а = +Іа + а .

гор у х у

Обробка інформації ІВОД при дії вертикально орієнтованого вектора прискорення

У випадку дії на ІВОД прискорення вздовж осі OZ, бічних зсувів кола обігу маятника не відбувається і симетрія інтервалів часу перебування маятника обабіч горизонтальних осей не порушується. Одначе, при цьому буде змінюватися сама тривалість періоду обертання конічного маятника.

Тривалість періоду, за умови коли на ІВОД діє прискорення а2, позначимо як T (на відміну від Т0, коли az = 0). Тоді в якості інформаційного параметра для визначення прискорення а2, можна використовувати різницю тривалостей згаданих періодів Лт2 = Т0 - Т2.

Період обертання маятника при дії прискорення а2 визначає вираз [6]:

Т =

2п

\ш±а,

mb3 L

Т

То

ґ т Л2 То

V 2п V

1

az L z

(5)

тут знак «+» відповідає випадкові, коли напрямок складової вектора a2 збігається з напрямком осі OZ.

Використовуючи (5), запишемо вираз для Лт2, а розв’язавши його відносно складової a отримаємо:

a.

L

Ґ 2п Л 2 1 - (т Л ±оі 2

V Т0 У Т V Lz V

(6)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

38

Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №63

Обчислювальні методи в радіоелектроніці

Як і раніше, будемо вважати, що вимірювання періодів Т0 і Т, здійснюється шляхом підрахунку числа імпульсів високочастотного генератора n0 i nz, якими заповнюються ці часові інтервали: Т0 = n0T^4, a Tz = nzT^4. Різниця тривалостей вимірюваних періодів визначиться як:

Ах z =АпТгвч = -Ап-,

J ГВЧ

де А n = n0 - nz.

Порогова величина прискорення буде виміряна за умови: А n = 1. Враховуючи, що (Ахz) = ТГВЧ << Т0, із виразу (6) отримаємо:

( az) = L

V z / пор

Л2Ґ 2п

Т

V Т0 у

Ах,

Т

0 у

8ж2 L

Т

ГВЧ

(6а)

При визначенні максимального значення прискорення, яке може бути виміряним вздовж осі OZ, маємо певну невизначеність, оскільки немає чітких фізичних обмежень ні на саму величину прискорення, ні на можливість її визначення. Приймемо величину () такою, за якої період обертання маятника зменшиться, наприклад, вдвічі, тобто:

(Ах ) = Т-(Т) = Т - 0,5Т = 0,5Т.

V z /max 0 V z / max 0 z 0 5 0

Враховуючи сказане, із виразу (6) отримаємо:

( az )

max

.2 L

12n2^

r\

(6б)

Динамічний діапазон прискорень, які можуть бути виміряними вздовж осі OZ за допомогою ІВОД:

(a) 3

Дт _ У z /max _ f

N z = / \ = ~ Т 0 J ГВЧ

(az) 2

v z / пор

3 Т 2 Т

2 1 ГВЧ

Як бачимо, динамічні діапазони вимірюваних прискорень для горизонтальних та вертикальної осей ІВОД виявилися вельми близькими.

При дії на ІВОД прискорення, напрямок вектора якого співпадає з напрямком вертикальної осі OZ, визначення значення складової прискорення a, потребує: вимірювання тривалості періоду обертання маятника Т2; обрахування різниці періодів обертання: Ат2 = Т0 - Т; обрахування за виразом (6) значення a .

Обробка інформації ІВОД при дії довільно орієнтованого вектора прискорення

Ситуація, коли вектор прискорення має довільну орієнтацію означає, що в загальному випадку ненульовими є всі три проекції цього вектора на

Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №63

39

Обчислювальні методи в радіоелектроніці

осі координат. За таких умов при вимірюванні горизонтальних складових прискорення отримувані значення ах та а будуть залежати від того, яке

значення має складова прискорення а2. Дійсно, поява вертикальної складової вектора прискорення а2 обумовлює згідно з (5), зміну періоду обертання маятника: T0 ^ Tz. І якщо тепер при вимірюваннях горизонтальної складової вектора прискорення, наприклад, а користуватися алгоритмом дій, згідно з п.1, то при цьому отримаємо вже не те значення Дт , яке відповідало умові az = 0, а деяке інше його значення (Дтy) , підстановка якого у вираз (4) обумовить спотворення значення а ■

Таким чином, в процесі вимірювання горизонтальних складових вектора прискорення необхідно коригувати їх значення відповідно зі значеннями вертикальної складової вектора прискорення

Використовуючи вирази (3) і (5), запишемо вираз для різниці інтервалів часу Дт за умови дії на ІВОД вертикальної складової прискорення а2:

( Дт У )- =

2T

0

ПІ ±

ґ Т Л2

TSL

V 2п j

шrcsm <

L

а

T2

-і с\

4 п 2

( T ^ 2 і 1

І ± 0 ■ а

V 2п J L z

а

(7)

Аналіз рівняння (7) показує, що при дії на ІВОД складової прискорення а, неспотворене значення складової прискорення а може бути знайденим, як розв’язок цього рівняння:

1 ±

ау = 4П R

ґ т \2

TL

V 2п J

1

L

а

T

■ sm

0

(T) 21

І ± 0 ■ ~ а

V V 2п J L z

2T

КУ

(8)

Як і раніше, все мовлене стосовно складової прискорення а повною мірою стосується і складової прискорення а ■

Сформулюємо загальний алгоритм проведення вимірювань всіх складових довільно орієнтованого вектора прискорення, що діє на ІВОД. Для цього необхідно виконати такі операції:

• вимірювання тривалості періоду обертання маятника Tz;

• обрахування різниці періодів обертання: Дт2 = T0 - T;

• обрахування значення а2 через значення T, і Дт2;

40

Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №63

Обчислювальні методи в радіоелектроніці

вимірювання тривалості інтервалів часу (т \ і (т2) перебування

маятника окремо з кожного боку осі ОХ;

обрахування різниці часових інтервалів (Лт^ ) = (т )z - (т2 )z;

через а2, (Лт^ ) , T і вираз (8), обрахування значення ау. аналогічним чином обрахувати значення ах.

Оцінка ефективності обробки інформації ІВОД

Дамо оцінку ефективності запропонованої корекції сигналів ІВОД. Для цього порівняємо значення горизонтальних складових вектора прискорення, отриманих без обробки інформації та за умови наявності такої обробки.

Для цього досить скласти відношення значень а , обрахованих за виразами (4) і (8), і дослідити поведінку цього відношення при зміні величини а2. Запишемо це відношення і після елементарних перетворень отримаємо:

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

а(уУ')

K = у

п (8)

ау

sin

ҐП-^Ту ^

V 2T0 У

(T ^ 2 1 1

1 ± 0 •—а

V 2п у L z

■sin

п

і

1 ±

ґ т \2

То_

V 2п У

L

а

2T

'(Ч.) z

0

(9)

де а.) та а{у] - горизонтальні складові прискорення, обраховані, відповідно, за виразами (4) і (8).

Отриманий результат розглядаємо як проміжний, тому, що він не придатний для аналізу, оскільки в ньому не відображено в явному вигляді

зв’язок між величинами Лт^ і (Лт^ ) . Пов’яжемо їх між собою.

Як було показано (див. вирази (5) і (7)), поява вертикальної складової прискорення а, обумовлює зміну періоду обертання маятника (T0 ^ Tz),

внаслідок чого Лт^ трансформується в (Лт^) . Очевидно, різниця інтервалів часу перебування маятника обабіч осі ОХ є пропорційною періоду його обертання за будь-якої величини цього періоду. Виходячи з цього, можна записати:

Лт у = Та =

КІ t \

1 ±

ґ т \2 т0

V 2п у

L

■а.

Користуючись цим виразом, можна суттєво спростити вираз (9):

1

Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №63

41

Обчислювальні методи в радіоелектроніці

а Г K = V п(8) ay

1

1 ±

С т \ To

V 2п J

(9а)

L

-а„

Одначе, для випадку вимірювача надмалих прискорень цей вираз буде малоінформативним, адже абсолютна величина коефіцієнта K буде завжди дуже близькою до одиниці. З огляду на це, доцільно ввести новий коефіцієнт & , який буде відображати тільки ту частину коефіцієнта K , на

яку він, власне, і буде відрізнятися від одиниці:

1 (10)

кУ = 1 - Ky =

1 ±

3EI 1

mL2 а.

Оцінка числового значення коефіцієнта

3EI

mL

у виразі (10) для акселе-

-2 2 -1

рометра надмалих прискорень дає: 2,640 c •м . Підрахунки проводилися для таких даних: E = 7,3 1010 Н-м-2; I = 1,2-10-17 м4; m = 0,01 кг; L = 0,1 м, що є параметрами маятника на основі стандартного кварцового ВС [2]. За цих параметрів період вільного обертання маятника Г0 = 12,2 с.

Практично у всьому діапазоні змін вертикальної складової вектора при, . , ч 8л2L 1

скорення (від (az ) =

= 4 -10 11 м - с-2 до az = 2,6 -10 3м - с 2)

-2

т3 f

1 о J Г

ГВЧ

другий доданок, що стоїть в знаменнику виразу (10), значно перевищує одиницю. З огляду на це, у вказаному діапазоні змін az можна користуватися спрощеним виразом для обрахунку & :

, mL

kVZ =--az

yz 3EI

(10а)

Рис. 3. Залежність коригуючого доданку kyz від вертикальної складової aZ.

який, власне, відображає пряму пропорційну залежність к (az), що добре видно

на рис. 3.

З першого погляду на рис.3 здавалося б, що мала величина ко-

ригуючого доданку цілком виправдовує використання спрощених вира-

2

42

Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №63

Обчислювальні методи в радіоелектроніці

зів (як наприклад, (10а)) для полегшення процедури проведення обрахунків. Такий підхід дійсно може бути прийнятним, але тільки тоді, коли не висувається високих вимог до прецизійності вимірювань.

В правій частині графіка (де вертикальна складова прискорення змінюється від а = 2,6• 1C3м• с“2 до (а ) = 12л2= 8• 10“2м • с“2) і стає поміт-

z \ z ' max гл2 у

T0

ним відхилення залежності k (az) від прямої пропорційності (пунктирна

крива лінія на рис.3), при розрахунках коригуючого доданку k необхідно

користуватися повним виразом (10). Слід мати на увазі і те, що динамічний діапазон значень прискорень, вимірюваних за допомогою ІВОД, може сягати 10-ти порядків, а тому при виконанні прецизійних вимірювань підрахунки мають вестися за точними формулами і з використанням чисел відповідної розрядності.

Висновки

1. Очікувані параметри вимірювача надмалих лінійних прискорень на основі ІВОД є унікальними: згідно з виразами (4б) і (6а) розрахункові порогові чутливості складають: ~ 10- м-с- та ~ 10- м^с- , відповідно, вздовж горизонтальної та вертикальної осей.

Максимальні значення прискорень, які можна виміряти вздовж тих же осей, згідно з виразами (4а) і (6б), складають: ~ 10- м^с- та ~ 10- м^с- .

2. Показано принципову можливість визначення неспотворених значень складових вектора прискорення вздовж усіх осей чутливості ІВОД.

Для цього розроблено метод математичної корекції вихідного сигналу ІВОД та загальний алгоритм опрацювання складових довільно орієнтованого вектора прискорення.

3. На прикладі вимірювача надмалих лінійних прискорень проведена оцінка ефективності корекції сигналів ІВОД у всьому діапазоні вимірювань. Показано, що коли значення вертикальної складової прискорення az не перевищують 2,6-10- м-с- , опрацювання сигналу ІВОД можна проводити за спрощеними формулами.

Перелік посилань

1. Soref R. Tilting-mirror Fiber Optic Accelerometer / R. Soref, D. H. McMahon // Appl. Opt. - 1984. - Vol. 23. - pp. 486-491.

2. Spillman W. D. Multimode Fiber Optic Accelerometer Based on the Photoelastic Effect / W. D. Spillman // Appl. Opt. - 1982. - Vol. 21. - pp. 2653-2658.

3. Фара Т. Транспортируемый гравиметр на холодных атомах разработки LNE-SYRTE: работа в подземных условиях в режиме наилучшей чувствительности / Т. Фара, К. Герлен, А. Ландражен, Ф. Буйе, С. Гаффе, Ф. Перейра дос Сантос, С. Мерле // Гироскопия и навигация. - 2014. - T. 86, №3. - с. 3-14.

4. Высокочувствительные акселерометры [Электронный ресурс]. - Режим доступа: http://zetview.com/catalog/vibrodats/icp/sense.php

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

5. Демьяненко П.А. Предельные возможности аналоговых оптических датчиков в

Вісник Національного технічного університету України «КПІ»

Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №63 43

Обчислювальні методи в радіоелектроніці

составе ВОС / П. А. Демьяненко // Радиотехника. - 1988. - №2, с.88-90.

6. Демьяненко П. А. Прецизионный цифровой акселерометр с волоконнооптическим датчиком / П. А. Демьяненко, Ю. Ф. Зиньковский, М. И. Прокофьев // Известия вузов. Радиоэлектроника. - 1997. - Том 40, №1. - с.39-47.

References

1. Soref R. and McMahon D. H. (1984) Tilting-mirror Fiber Optic Accelerometer. Appl. Opt., Vol. 23, pp. 486-491.

2. Spillman W. D. (1982) Multimode Fiber Optic Accelerometer Based on the Photoelastic Effect. Appl. Opt., Vol. 21, pp. 2653-2658.

3. Farah T., Guerlin C., Landragin A., Bouyer Ph., Gaffet S., Pereira Dos Santos F. and Merlet S. (2014) Underground Operation at Best Sensitivity Of The Mobile Lne-Syrte Cold Atom Gravimeter. Giroskopiya i navigatsiya, Vol. 86, No 3, pp. 3-14.(in Russian)

4. Highly sensitive accelerometers. Available at: http://zetview.com/catalog/

vibrodats/icp/sense.php

5. Dem'yanenko P.A. (1988) Predel'nye vozmozhnosti analogovykh opticheskikh datchikov v sostave VOS [Extreme performance of analog optical sensors as a part of a fiber optical system]. Radiotekhnika, No 2, pp. 88-90.

6. Dem'yanenko P. A., Zin'kovskii Yu. F. and Prokofev M. I. (1997) Pretsizionnyi tsifrovoi akselerometr s volokonno-opticheskim datchikom [Precision digital accelerometer with a fiber-optic sensor]. Izvestiya vuzov. Radioelektronika, Vol 40, No 1, pp.39-47.

Дем’яненко П.О., Зіньковський Ю.Ф. Опрацювання сигналів імпульсного волоконно-оптичного давача прискорення. Розроблено алгоритм опрацювання сигналів тривимірних волоконно-оптичних давачів прискорення з імпульсною модуляцією інтенсивності оптичного потоку, яким притаманне явище взаємного впливу складових вектора прискорення, котрі спрямовані вздовж різних осей чутливості давача. Застосування розробленого алгоритму опрацювання в акселерометрах, побудованих на основі таких давачів, забезпечує отримання неспотворених значень всіх складових вектора прискорення спрямованих вздовж різних осей чутливості давача. Дана кількісна оцінка ефективності запропонованої обробки інформації для випадку проведення прецизійних вимірювань прискорення.

Ключові слова: волоконно-оптичні давачі прискорення, імпульсна модуляція, волоконно-оптичні акселерометри.

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

Ключевые слова: волоконно-оптические датчики ускорения, импульсная модуляция, волоконно-оптические акселерометры.

Вісник Національного технічного університету України «КПІ» 44 Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №63

Обчислювальні методи в радіоелектроніці

Dem'yanenko P.O., Zinkovskiy Y.F. Signal processing of the optical fiber impulse sensor acceleration.

Introduction. Feature of the pulsed fiber-optic sensor is a mutual influence of the components of the acceleration vector measured. When processing the signals of the pulse sensor should provide adequate correction value counts for correct measurement results.

Main part. It consists of four sections.

The first section outlines the method of determining the undistorted values of horizontal components of the acceleration vector acting on the sensor, in the absence of the vertical component of the acceleration vector.

The second section outlines the method of determining the undistorted value of vertical component of the acceleration vector.

Based on the results obtained in the first two chapters, the third chapter described a method of determining undistorted measured quantities of all components of the acceleration vector.

In the fourth chapter provides a quantitative assessment of the effectiveness of the proposed information processing in the case of precision measurements of acceleration.

Conclusions. It is pointed to the principal possibility of determining the undistorted values of samples of the vector acceleration components along all axes of sensitivity of the pulse sensor.

The method of determination is proposed and a general algorithm is formulated for measuring randomly oriented components of the acceleration vector.

It is evaluated the efficiency of correction of pulse sensor signals for example accelerometer of micro-accelerations.

Keywords: fiber-optic acceleration sensors, pulse modulation, fiber optic accelerometers.

Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2015. — №63

45

i Надоели баннеры? Вы всегда можете отключить рекламу.