Спосіб визначення різниці фаз

Реферат: Спосіб визначення різниці фаз полягає в тому, що промінь від джерела когерентного випромінювання перетворюють на два пучки. Оптичною системою формують зображення контрольованого об'єкта на матриці цифрового реєстратора, підключеного до пристрою зберігання і обробки даних. Реєструють послідовно за допомогою матриці розподіл інтенсивності першого пучка, розподіл інтенсивності другого пучка, розподіл інтенсивності при одночасному використанні першого і другого пучків для початкового і зміненого стану об'єкта контролю. На основі даних вимірювань інтенсивностей для кожного елемента матриці реєстратора розраховують значення фаз і різниць фаз, відповідних різним станам об'єкту. В процесі послідовної реєстрації зображень за допомогою масок частково перекривають кожний з пучків так, щоб області перекриття першого і другого пучків відповідали різним ділянкам матриці цифрового реєстратора, порівнюють інтенсивності для елементів матриці реєстратора, що знаходяться в області кожної з масок, на зображеннях, зареєстрованих з одним пучком, і зображеннях, зареєстрованих з двома пучками. Враховують зміну інтенсивності для даних ділянок при розрахунку фаз. UA 72975 U (54) СПОСІБ ВИЗНАЧЕННЯ РІЗНИЦІ ФАЗ UA 72975 U UA 72975 U 5 10 15 20 25 30 35 40 45 50 55 60 Корисна модель належить до оптичних методів контролю і може бути використана для неруйнівного контролю і технічної діагностики виробів машинобудування. Зокрема, даний спосіб може бути використаний для досліджень на основі методів голографічної і спеклінтерферометрії, а також ширографії, при умові реєстрації за допомогою матричних систем цифрового запису зображень, які відображають зміни, що відбуваються в об'єкті контролю або в зоні контролю. Дані методи контролю базуються на визначенні величини різниці фаз світлових пучків, що формують зображення контрольованого об'єкта в площині матриці світлочутливої системи реєстрації, і зміни цієї різниці фаз внаслідок деформування об'єкта контролю. В даний час існують різноманітні способи визначення різниці фаз пучків, що формують зображення. Зокрема, широко відомий метод фазових зсувів, який передбачає використання пристроїв, заснованих на цифровій реєстрації зображення об'єкта (ПЗЗ-матриць), і високошвидкісних пристроїв введення, зберігання і обробки відеоінформації і дозволяє розраховувати різницю фаз пучків, що інтерферують, використовуючи серію послідовних кадрів, одержаних при контрольованій зміні фази одного із світових пучків, задіяних у формуванні зображень [1-5]. Зміна фази може здійснюватися різними способами. Найбільш поширеним способом зміни фази світлового пучка є керований зсув або обертання елементів оптичної схеми (дзеркал, лінз, інших оптичних елементів), що формують пучок [6, 7]. Крім того, керована зміна фази пучка може бути здійснена шляхом розташування в каналі його розповсюдження багатоступінчатої матриці, що складається з прозорих пластин із змінною за заданим законом товщиною. При цьому зміна оптичної довжини шляху для даного пучка приводить до зміни різниці фаз світових пучків, що реєструються матрицею [8]. Одним із різновидів методу зсуву фази є використання при формуванні зображення світлових пучків різної поляризації [9]. Способи визначення фази, засновані на застосуванні керованого зсуву фази або на використанні поляризаційних властивостей пучків, що беруть участь у формуванні зображення на матриці реєстратора, вимагають дуже точного юстування і якості виготовлення елементів оптичної схеми установки, ретельного контролю процесів управління величиною зсуву оптичних елементів і зміною поляризації світла. Будь-яка зміна величини зсуву фази, що обумовлена нестаціонарністю процесів в ланцюгах управління механізмами зсуву і повороту, механічним зносом елементів управління зсувом, або зміною властивостей матеріалів ступінчастих матриць, використовуваних для створення фазового зсуву, призводить до зниження точності вимірювань і, отже, до зниження достовірності результатів контролю. Крім того, метод фазового зсуву не може бути використаний у разі дослідження процесів, які не забезпечують точне відтворення стану контрольованого об'єкту, зокрема, при дослідженні процесів повзучості і текучості. Відомий також спосіб визначення фази, обраний як прототип, згідно з яким промінь від джерела когерентного випромінювання перетворюють у два пучки, за допомогою оптичної системи формують зображення контрольованого об'єкта на матриці цифрового реєстратора, підключеного до пристрою зберігання і обробки даних, реєструють послідовно за допомогою матриці розподіл інтенсивності першого пучка, розподіл інтенсивності другого пучка, розподіл інтенсивності при одночасному використанні першого і другого пучків для початкового і зміненого стану об'єкта контролю, на основі даних вимірювань інтенсивностей для кожного елемента матриці реєстратора розраховують значення фаз і різниць фаз, відповідних різним станам об'єкту [10]. Алгоритм обробки дозволяє обчислити значення різниці фаз, що характеризує деформацію об'єкта під впливом навантаження, без використання методу зсуву фази, що дозволяє усунути всі недоліки, пов'язані з точністю управління величиною зсуву. Крім того, оскільки для визначення різниці фаз використовується всього одне зображення, відповідне даному стану об'єкта, усувається проблема відтворюваності стану об'єкта, а отже розширюється область застосування методу. Перевагою даного методу у порівнянні із методом зсуву фаз є також скорочення кількості зображень, необхідних для отримання даних, які використовуються для розрахунку фази. Це приводить до істотного прискорення процесу контролю без зниження точності результатів контролю. Описаний спосіб є найбільш близьким по технічній суті до корисної моделі, що заявляється. Проте він, як і методи зсуву фази, не дозволяє вирішити проблему усунення погрішності визначення фази, обумовленої зміною інтенсивності пучків, що реєструються ПЗЗ-матрицею в різні моменти часу, тобто при роздільній реєстрації кожного із пучків та при їх одночасній реєстрації для різних станів навантаження об'єкта контролю. Зміна інтенсивності пучків можлива внаслідок зміни рівня загальної освітленості зони контролю, а також внаслідок флуктуацій інтенсивності лазерного випромінювання. Крім того, не передбачений механізм 1 UA 72975 U 5 10 15 20 25 30 35 40 45 50 контролю зміни часу експозиції при реєстрації зображень, пов'язаний з особливостями процесу реєстрації сигналу в ПЗЗ-матриці. В основу корисної моделі поставлена задача удосконалення процесу реєстрації інтерферограм з метою виключення похибок, пов'язаних з флуктуаціями інтенсивності лазерного випромінювання і зміною умов освітленості контрольованого об'єкта. Вирішення цієї задачі дозволить підвищити достовірність і точність результатів контролю. Поставлена задача вирішується тим, що в способі, який полягає в тому, що промінь від джерела когерентного випромінювання за допомогою системи оптичних елементів перетворюють у два пучки, за допомогою оптичної системи формують зображення контрольованого об'єкта на матриці цифрового реєстратора, підключеного до пристрою зберігання і обробки даних, реєструють послідовно за допомогою матриці розподіл інтенсивності першого пучка, розподіл інтенсивності другого пучка, розподіл інтенсивності при одночасному використанні першого і другого пучків для початкового і зміненого стану об'єкта контролю, на основі даних вимірювань інтенсивностей для кожного елемента матриці реєстратора розраховують значення фаз і різниць фаз, відповідних різним станам об'єкту, згідно з корисною моделлю в процесі послідовної реєстрації зображень при одночасному освітленні матриці реєстратора двома пучками за допомогою масок частково перекривають кожний з пучків так, щоб області перекриття першого і другого пучків відповідали різним ділянкам матриці цифрового реєстратора, порівнюють інтенсивності для елементів матриці реєстратора, що знаходяться в області кожної з масок, на зображеннях, зареєстрованих з одним пучком, і зображеннях, зареєстрованих з двома пучками, і враховують зміну інтенсивності для даних ділянок при визначенні коригуючих коефіцієнтів і подальшому розрахунку фаз. Використання алгоритму порівняння інтенсивності в області масок на знімках, одержаних в різні моменти часу, дозволяє врахувати всі зміни інтенсивності пучків, що реєструються ПЗЗматрицею в різні моменти часу, обумовлені як флуктуацією інтенсивності когерентного випромінювання, так і зміною загальної освітленості зони контролю. Крім того, даний алгоритм дозволяє врахувати також зміни, обумовлені можливою різницею в часі експозиції при послідовній реєстрації знімків. Спосіб дозволяє здійснювати розрахунок різниці фаз незалежно від того, який з методів буде використаний в процесі експериментальних досліджень, тобто він може бути застосований і в голографічній, і в спекл-інтерферометрії, і в ширографії. На фіг. 1, 2 наведені креслення, які пояснюють умови здійснення реєстрації зображень в ході контролю, де на фіг. 1 наведена експериментальна схема для контролю із застосуванням методу голографічної інтерферометри, а на фіг. 2 представлена картина розподілу інтенсивності в площині матриці реєстратора на етапі реєстрації зображення вихідного стану об'єкта контролю із застосуванням двох пучків одночасно. При здійсненні контролю з використанням методу голографічної інтерферометри заявлений спосіб реалізовується таким чином. Промінь від джерела випромінювання 1 за допомогою світлодільника 2 розбивається на два напрями - опорний і предметний (див. фіг. 1). За допомогою оптичних елементів 3, 4 формують предметний пучок, який освітлює об'єкт контролю 5, і опорний пучок. Розсіяний об'єктом контролю 5 предметний пучок і опорний пучок за допомогою системи формування зображення 6 направляють на реєстратор 7. Інформація про розподіл інтенсивності в площині матриці цифрового реєстратора зображень вводиться в пристрій обробки інформації, де шляхом аналізу одержаних зображень за певними алгоритмами перетворюється в картину розподілу різниці фаз і потім у поле переміщень поверхні об'єкта. На світлодільник 8, розташований в системі формування зображення, наносять дві маски 9 і 10, що не перекриваються. Одна з них розташовується з боку предметного пучка (предметна маска 10), а інша - з боку опорного пучка (опорна маска 9). Перекривають за допомогою екрана опорний пучок повністю і реєструють розподіл інтенсивності по перетину предметного пучка в області реєстратора. Розподіл інтенсивності   пр одержують у вигляді матриці дійсних чисел Iij розміром n×m, де n і m є характеристиками ПЗЗ-матриці. Перекривають повністю за допомогою екрана предметний пучок і реєструють 55   оп розподіл інтенсивності у вигляді матриці Iij по перетину опорного пучка в області реєстратора. Прибирають екрани і реєструють зображення контрольованого об'єкта в початковому стані    деф (матриця Iij ), а потім в деформованому (матриця Iij ). Для обчислення фазового зсуву слід скористатися співвідношеннями 2 UA 72975 U оп оп Iij  Iпр  Iij  2 Iпр Iij cos  ij ; ij ij   ( деф оп оп Iij  Iпр  Iij  2 IпрIij cos  ij  ij 1) ij ij , 5 де γ - мультиплікативний чинник, який описує вплив некогерентної складової і обчислюється шляхом аналізу отриманих величин [10]; αij - фаза, яка являє собою величину, що випадковим чином змінюється від піксела до піксела і є проміжною величиною при обчисленнях; φij - різниця фаз, яка являє собою величину, пов'язану з переміщеннями об'єкта контролю між останніми експозиціями і є метою обчислень.    деф Формули (1) справедливі в тому випадку, якщо до моменту отримання матриць Iij і Iij     елементи матриць Iпр і Iпр не змінилися. Флуктуації лазерного випромінювання і зміна ij ij загального освітлення зони контролю та часу експонування матриць, призводять до зміни 10      оп елементів цих матриць, їх нові елементи при отриманні Iij позначені відповідно як Iпр 1 і Iij 1 , а ij   деф Iij при отриманні - відповідно мають наступний вигляд:   і I  . Кориговані з урахуванням цих змін формули (1) Iпр 2 ij оп 2 ij оп оп Iij  Iпр1  Iij 1  2 Iпр1Iij 1 cos  ij ; ij ij   деф оп оп Iij  Iпр 2  Iij 2  2 Iпр 2Iij 2 cos  ij  ij , ij ij звідки можна обчислити різниці фаз за формулою 15 ij  arccos оп Iij  Iпр 2  Iij 2 ij оп 2 Iпр 2Iij 2 ij  arccos оп Iij  Iпр1  Iij 1 ij оп 2 Iпр1Iij 1 ij .(2)     пр пр Перехід від старих елементів матриць Iij і Iij до нових можна здійснити на основі математичної моделі тих змін, які з ними відбулися. У загальному випадку зміни можуть бути представлені функціональними залежностями Iпр1  f1 Iпр , p1, p 2 ,...,p k ; ij ij Iпр 2 ij оп Iij 1 оп Iij 2 20 25   f I  f I  f I  , p , p ,...,p ; , p , p ,...,p  ; , p , p , ...,p  , пр 2 ij 1 2 k оп 3 ij 1 2 k оп 4 ij 1 2 k де f1, f2, f3, f4 - відомі функції, які визначаються шляхом експериментальних досліджень; p1, p2,…, рk - наперед невідомі коригуючі коефіцієнти. Оскільки під час реєстрації світлові пучки частково перекриваються за допомогою масок, то на зображенні, яке фіксується матрицею реєстратора, будуть присутні три області. Позначимо літерою А ділянку матриці реєстратора, яка внаслідок застосування опорної маски реєструє лише інтенсивність предметного пучка, а літерою В - ділянку матриці реєстратора, яка внаслідок застосування предметної маски реєструє лише інтенсивність опорного пучка (див. фіг. 2). Для позначення цих вікон на зображенні об'єкта контролю після його навантаження використовуються літери С і D відповідно. Підсумовуючи квадрати відхилень передбачуваних елементів матриці від отриманих за всіма відміченими ділянками, одержимо   J   Iij  f1 Iпр , p1, p 2 ,...,p k ij A  B 30 (3)     2 оп Iij  f3 Iij , p1, p 2 ,...,p k 2   деф   Iij  f 2 Iпр , p1, p 2 ,...,p k ij C  D  деф Iij    2 jg  f 4 Iij , p1, p 2 ,...,p k 2  ,(4)  де E - сума за точками ділянки Е, де Е приймає значення А, В, С і D. Диференціюючи величину J за коригуючими коефіцієнтами р1, р2,…, pk і прирівнюючи результат нулю, отримаємо систему алгебраїчних рівнянь, розв'язання якої дає значення коригуючих коефіцієнтів р1, р2,…, рk. Обчислені значення коефіцієнтів використовуються для         оп оп знаходження за формулою (3) елементів матриць Iпр 1 , Iij 1 , Iпр 2 і Iij 2 . Виходячи з отриманих ij ij 3 UA 72975 U 5 внаслідок цього розрахунку величин, за формулою (2) обчислюють значення різниці фаз. Сама величина J при визначених таким чином значеннях коригуючих коефіцієнтів свідчить про ступінь відповідності використовуваної математичної моделі реальним вимірюванням і може бути використана для відбраковування результатів вимірювань шляхом порівняння із пороговим значенням. Найпростіша математична модель, яка відповідає пропорційній зміні інтенсивності освітлення матриці реєстратора, може бути представлена рівняннями: Iпр1  Iпр p1 ij ij ; Iпр 2  Iпр p 2 ij ij ; ( оп Iij 1 5) оп  Iij p 3 ; оп оп Iij 2  Iij p k . Використовуючи рівняння (5) для розрахунку J і використовуючи умову рівності нулю похідної, отримують наступні формули для визначення коригуючих коефіцієнтів р1, р2, р 3, р4: p1   Iij  Iпр ij A A ; деф p 2   Iij  Iпр ij C C ; ( p3   Iij  B B оп Iij 6) ; деф p 4   Iij D оп  Iij D . 10 Визначені за допомогою (6) значення р1, р2, р3, р4 використовують для розрахунку різниці фаз, що може бути знайдена наступним чином:  ij  arccos 15 20 25 30 оп Iij  p 2Iпр  p 4Iij ij оп 2 p 2 p 4Iпр Iij ij  arccos оп Iij  p1Iпр  p 3Iij 1 ij . оп 2 p1p 3Iпр1Iij 1 ij Джерела інформації: 1. Карпюк Б.В., Солодкин Ю.Н. Анализ погрешностей измерения фазы интерферометром с управляемым фазовым сдвигом // Автометрия.-1992. - № 6. - c. 16-21. 2. Ильиных С.П., Гужов В.И. Обобщенный алгоритм расшифровки интерферограмм с пошаговым сдвигом // Автометрия.-2002. - Т. 38, № 3. - c. 123-126. 3. Гужов В.И. Практические аспекты измерения фазы в интерферометрии // Автометрия.1995, - № 5. - c. 25-31. 4. Гуров И.П. Компьютерная обработка интерференционных сигналов на основе алгоритма управляемого фазового сдвига // Оптический журнал. 1998. N10. - c. 38-42. 11 5. Патент України 93335 МПК G01В 9/021 Імерсійний метод отримання топограм поверхонь дифузно розсіюючих об'єктів / Сминтина В.А., Тюрин О.В., Попов А.Ю. та ін. - № а201006137, Заявлено 20.05.2010; Опубл. 25.01.2011 // Бюл. № 2.-2011. 6 6. US Patent 7274467 B1 Int. Cl. G01B 11/02 Phase shifting interferometric method, interferometer apparatus and method of manufacturing an optical element / Bernd Doerband, Stefan Schulte. - Appl. No. 11/027989. Filed: 04.01.2005. Date of Patent: 25.09.2007. 6 7. US Patent 5473434 Int. Cl. G01B 9/02 Phase shifting interferometer and method for surface topography measurement / Peter de Groot. - Appl.No.242789. Filed: 16.05.1994. Date of Patent: 05.12.1995. 6 8. US Patent 7119910 B1 Int. Cl. G01B 11/02 Phase shifting wavefront interference method / Uve Schellhorn. - Appl. No. 10/697239. Filed: 31.10.2003. Date of Patent: 10.10.2006. 6 9. US Patent 7379188 B2 Int. Cl. G01В 9/02 Phase shift interferometer / Maarten Jansen. - Appl. No. 11/487430. Filed: 17.06.2007. Date of Patent: 27.05.2008. 4 UA 72975 U 6 10. US Patent 7280187 B1 Int. Cl. G01В 9/02 Method for resolving phase in electronic speckle interferometry / Sanichiro Yoshida. - Appl. No. 11/248856. Filed: 12.10.2005. Date of Patent: 09.10.2007. 5 10 15 20 ФОРМУЛА КОРИСНОЇ МОДЕЛІ Спосіб визначення різниці фаз, який полягає в тому, що промінь від джерела когерентного випромінювання перетворюють на два пучки, за допомогою оптичної системи формують зображення контрольованого об'єкта на матриці цифрового реєстратора, підключеного до пристрою зберігання і обробки даних, реєструють послідовно за допомогою матриці розподіл інтенсивності першого пучка, розподіл інтенсивності другого пучка, розподіл інтенсивності при одночасному використанні першого і другого пучків для початкового і зміненого стану об'єкта контролю, на основі даних вимірювань інтенсивностей для кожного елемента матриці реєстратора розраховують значення фаз і різниць фаз, відповідних різним станам об'єкта, який відрізняється тим, що в процесі послідовної реєстрації зображень за допомогою масок частково перекривають кожний з пучків так, щоб області перекриття першого і другого пучків відповідали різним ділянкам матриці цифрового реєстратора, порівнюють інтенсивності для елементів матриці реєстратора, що знаходяться в області кожної з масок, на зображеннях, зареєстрованих з одним пучком, і зображеннях, зареєстрованих з двома пучками, і враховують зміну інтенсивності для даних ділянок при розрахунку фаз. 5 UA 72975 U Комп’ютерна верстка Л.Литвиненко Державна служба інтелектуальної власності України, вул. Урицького, 45, м. Київ, МСП, 03680, Україна ДП “Український інститут промислової власності”, вул. Глазунова, 1, м. Київ – 42, 01601