автореферат диссертации по приборостроению, метрологии и информационно-измерительным приборам и системам, 05.11.16, диссертация на тему:Алгоритмическое и программное обеспечение измерительного комплекса для определения параметров распространения различных типов акустических волн в твердых средах

Король Андрей Александрович

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

05.11.16 - информационно-измерительные и управляющие системы (техника и технология)

Автореферат диссертации на соискание ученой степени кандидата технических наук

1 2 МАЙ 2011

Хабаровск -2011

4845452

Работа выполнена в государственном образовательном учреждении высшего профессионального образования «Дальневосточный государственный университет путей сообщения» (г. Хабаровск)

Научный руководитель:

доктор технических наук, профессор Кондратьев Александр Иванович

Официальные оппоненты:

доктор физико-математических наук, профессор

Касаткин Борис Анатольевич

кандидат технических наук, доцент Овчарук Валерий Николаевич

Ведущая организация:

УРАН Институт автоматики и процессов управления ДВО РАН

Защита состоится 1 июня 2011 года в 15-00 на заседании диссертационного совета ДМ 212.294.05 в государственном образовательном учреждении высшего профессионального образования «Тихоокеанский государственный университет» по адресу: 680035, г. Хабаровск, ул. Тихоокеанская, 136, ауд. 315л.

С диссертацией можно ознакомиться в библиотеке Тихоокеанского государственного университета

Автореферат разослан: 2011г.

Ученый секретарь диссертационного совета

Бурдинский И.Н.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы

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

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

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

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

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

значениях коэффициента затухания до 2000 дБ/м). Предложенные методики позволяют повысить точность, а также расширить диапазон измерений соответствующих физических величин.

Цель работы

Целью настоящей работы является разработка алгоритмического и программного обеспечения измерительного комплекса для определения параметров распространения различных типов акустических волн в твердых телах.

Решаемые задачи

1. Анализ методов измерения акустических величин и методов возбуждения и приема у.з. колебаний.

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

3. Разработка алгоритмов измерения фазовой скорости распространения и коэффициента затухания у.з. колебаний при измерениях резонансным методом.

4. Разработка программных библиотек численной обработки данных и взаимодействия с пользователем для построения программного обеспечения с учетом требований, предъявляемых к программным комплексам средств измерений в процессе их метрологической аттестации.

5. Разработка программного обеспечения измерительного комплекса, реализующего предложенные алгоритмы.

6. Исследование метрологических характеристик информационно-измерительного комплекса с учетом разработанных алгоритмов и программного обеспечения.

Методы исследования

При выполнении работы использовались: теория и методы неразрушаю-щего контроля, теория дифференциальных уравнений, теория и математический аппарат преобразования Фурье, численные методы интегрирования и дифференцирования, теория объектно-ориентированного программирования; системное программирование для операционных систем (ОС) семейства Microsoft Windows.

Научная новизна

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

2. Разработаны алгоритмы определения времени прихода сдвиговой компоненты сигнала при наличии помехи, основанные на аппроксимации фронта сдвигового импульса различными функциями.

3. Разработана методика оценки дополнительной дифракционной поправки при измерении коэффициента затухания резонансным методом на частотах ниже 5 МГц для твердых сред.

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

5. Установлено, что при величине коэффициента затухания более 1000 дБ/м происходит эффект усиления полезного сигнала, обусловленный наличием помехи от задающего генератора анализатора спектра.

Практическая значимость

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

2. Сформирован набор программных библиотек для реализации различных численных методов и взаимодействия с пользователем, позволяющих создавать на их основе программное обеспечение средств измерений, удовлетворяющее требованиям, предъявляемым при метрологической аттестации.

3. Предложенная методика определения коэффициента затухания для резонансного режима позволяет повысить верхнюю границу диапазона измерений до значения 2000 дБ/м.

Основные положения, выносимые на защиту

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

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

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

4. Структура и состав программного обеспечения измерительного комплекса для определения параметров распространения акустических волн в твердых средах, обеспечивающего заданные метрологические характеристики.

Применение результатов работы

Разработанное алгоритмическое и программное обеспечение применяется в установках высшей точности, разрабатываемых в Дальневосточном филиале ФГУП «ВНИИФТРИ», которые в 2011-2013 годах будут переведены в ранг государственных эталонов.

Внедрение результатов диссертационной работы подтверждено соответствующим актом.

На основе разработанных алгоритмов и программных библиотек на базе государственного образовательного учреждения высшего профессионального образования «Дальневосточный государственный университет путей сообщения» был реализован и внедрен программный комплекс расчета степени локального разупрочнения контактного провода.

Апробация работы

Основные положения и отдельные результаты работы докладывались и обсуждались: на межрегиональной научно-практической конференции «Информационные и коммуникационные технологии в образовании и научной деятельности» (г. Хабаровск, 2008); на пятой международной научно-практической конференции «Исследование, разработка и применение высоких технологий в промышленности» (г. Санкт-Петербург, 2008); на V Всероссийской конференции «Механика микронеоднородных материалов и разрушение» (г. Екатеринбург, 2008); на XXXIV Дальневосточной математической школе-семинаре им. академика Е.В.Золотова (г. Хабаровск, 2009); на пятой Всероссийской конференции «Математическое моделирование и краевые задачи» ММР-2008 (г. Самара, 2008); на международной научной конференции «Оптика кристаллов и наноструктур» (г. Хабаровск, 2008); на VI Всероссийском симпозиуме «Физика геосфер» (г. Владивосток, 2009).

Публикации

Основное содержание диссертационной работы отражено в 10 публикациях, из них 2 публикации в изданиях из перечня ВАК, 5 докладов на научно-практических конференциях (2 - на международных), 2 свидетельства о государственной регистрации программ для ЭВМ.

Личный вклад соискателя

Личный вклад заключается в разработке и анализе предложенных алгоритмов и их программной реализации, в разработке программно-аппаратных средств автоматизации измерительных процессов. Также при участии автора были проведены экспериментальные исследования.

Структура и объем работы

Диссертационная работа представлена введением, четырьмя главами, заключением, списком использованных литературных источников и приложениями. Работа изложена на 140 страницах машинописного текста и содержит 6 таблиц, 38 рисунков, 137 наименований литературных источников, 3 приложения.

В приложениях представлены: свидетельство о государственной регистрации программы для ЭВМ, акт внедрения результатов диссертационной работы, фрагменты исходного кода основного программного модуля.

КРАТКОЕ СОДЕРЖАНИЕ ДИССЕРТАЦИОННОИ РАБОТЫ

Во введении обосновывается актуальность темы, сформулированы цели и задачи, показана научная новизна и практическая значимость работы.

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

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

Приведен анализ методов возбуждения и приема у.з. колебаний, сделан вывод о перспективности использования бесконтактных методов. Приведены особенности оптических и емкостных методов возбуждения и приема у.з. колебаний, проанализированы основные достоинства и недостатки рассмотренных методов.

Представлены блок-схемы установок (в том числе установок высшей точности), использующихся для проведения измерений параметров распространения различных типов акустических волн, приведены их характеристики и рассмотрены особенности проведения измерений.

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

1. Групповая скорость распространения продольных у.з. колебаний.

2. Скорость распространения сдвиговых у.з. колебаний.

3. Скорость распространения и коэффициент затухания продольных у.з. колебаний в зависимости от частоты.

При построении информационно-измерительного комплекса (ИИК) для определения параметров распространения у.з. колебаний в импульсном режиме может использоваться как оптическая установка УВТ 39-А-86, так и установка ИЗУ-1, в которой в качестве излучателей и приемников у.з. колебаний используются емкостные преобразователи.

При измерении групповой скорости распространения в импульсном методе за основу принято следующее выражение:

сгр=_Ы_ > (1)

' Т2~Т\

где с/ - толщина образца; т: и т2 - времена, соответствующие особым точкам первого и второго рассматриваемых импульсов соответственно; к = (2пг-1) ~ (2п-1)\ п, т- номера «переотражений» первого и второго рассматриваемых импульсов.

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

В разработанном ПО ИИК реализованы следующие способы определения искомой величины:

По «пикам» импульсов. В данном случае за г; и т2 в выражении (1) принимаются времена, соответствующие максимумам первого и второго рассматриваемого импульсов.

По переднему фронту на заданном уровне. Данный случай аналогичен случаю с измерением «по пикам», за исключением того, что в формуле (1) за т, и т2 берутся времена, соответствующие приходу заранее заданного уровня [1 переднего фронта рассматриваемых импульсов. Таким образом, для каждого из двух импульсов определяется время, когда на осциллографе фиксируется значение и р=Р' мтах > где итах~ максимум соответствующего импульса.

По центрам тяжести импульсов. Поиск центров тяжести импульсов производится по формуле:

где ¡а и ц - границы интегрирования, которые соответствуют границам появления импульса на осциллограмме. Считается, что функция сигнала и(1) за пределами границ интегрирования равна нулю. Формула вычисления групповой скорости совпадает с формулой (1), при этом за и т2 принимаются значения тцт! и тцт2, полученные путем нахождения центров тяжести первого и второго рассматриваемого импульса соответственно.

По локальному максимуму автокорреляционной функции на заданном интервале. Предложенная в настоящей работе методика основывается на понятии автокорреляционной функции:

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

(2)

+ 00

Г(т)= ¡и(0-и(( + т)&,

(3)

— со

тервале [т;; ъ], при этом значения и г? выбираются таким образом, чтобы функция вычислялась для наложения двух рассматриваемых импульсов друг на друга. Значение параметра т„жс, на котором достигается искомый максимум, подставляется в выражение (1) в качестве знаменателя.

Для вычисления скорости распространения сдвиговых колебаний используется следующее выражение:

г + 21 а

где с! - толщина образца, г - время между приходом первого и второго импульсов, /у - время между приходом первого импульса и появлением на осциллограмме сдвиговой составляющей (рис. 1).

и

и

-и

■V

Рис. 1. Осциллограмма, полученная на УВТ 39-А-86. 1.1, - продольные переотраженные импульсы; Б - сдвиговый импульс

При вычислениях по формуле (4) время г между импульсами Ь] и Ь2 определяется по локальному максимуму автокорреляционной функции, получаемому при определении скорости продольных волн.

Основная задача вычисления Су сводится к нахождению точки Б прихода сдвиговой компоненты. Ввиду наличия помех, в ИИК необходимо было реализовать отказоустойчивый алгоритм поиска Б. В работе предлагается использование следующих способов.

Линейное приближение. В данном способе происходит построение линейного приближения переднего фронта сдвигового импульса. За точку Б берется точка пересечения построенной прямой с нулевым уровнем.

Квадратичное приближение. В данном случае происходит построение квадратичной функции, аппроксимирующей передний фронт сдвигового импульса. За точку Б берется точка пересечения построенной функции с нулевым уровнем.

Приближение с помощью сглаживающей функции. Данный способ подразумевает приближение фронта сдвигового импульса следующей функцией 0(1):

\,t <a

0,t>b

, G{t)-kF{t),

(5)

{t - a){t -2b + a)

(b-a)2

/

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

Среди всех сумм выбирается минимальная 1¥а/, ,„,„, а соответствующие значения атт и Ьтт используются для построения функции Отт{(). Пересечение полученной функции с нулевым уровнем считается за искомую точку 8.

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

где a, b - коэффициенты, влияющие на форму кривой графика функции. В нашем случае в выражение (7) добавляется множитель h(t) - функция Хэвисайда, так как t > 0. Поиск искомой точки S аналогичен предыдущему варианту. Происходит перебор параметров а и b с построением сумм вида (6), где в качестве функции G(t) выступает функция вида (7), умноженная на масштабирующий коэффициент и смещенная в окрестность искомой точки.

Для определения скорости распространения и коэффициента затухания продольных у.з. колебаний в зависимости от частоты были получены соотношения:

о)2.

(6)

Fah{t) = -a-t2-e-h',

(7)

(\U,{a)\-\D{aM)\\

2d

(8)

Q(®) =

2d

(9)

L-L +

| ф {со) - ф 1<о) - (co,d) + Уд (a,3d)'

где U,(си), U2(co) - спектральные характеристики двух наблюдаемых импульсов; d - толщина образца; D(oj, х) - коэффициент, учитывающий дифракционные искажения сигнала; 1^(^)1 = abs[Uiiw)},(pi((i)) = argiUx (ы)},^; |i/2(<u)| = abs[U2(cj)}, <р2(ш) = arg{U2(co)}, t2 - АЧХ, ФЧХ и времена прихода первого и второго импульсов соответственно; Ч'о(шх) - фаза D(ojx)', со = 2 л/

Для определения D(co, х) использовалось следующее соотношение (полученное в работе Гитиса М.Б., Химунина A.C.):

D{cd,X) =

1 - A-Jj^j-B-Jj^e

-К

= \D{ü),x)\erir» '

J^D . (Ю)

где А=[\-^2/{1к2г2)1 В=[А+С /(к2г2)]; C = k{jxT+4r*-x}/2, х - путь

пройденный у.з, импульсом; Л(0> «Л (О ~ функции Бесселя; к=0)/С^ /"-радиус пятна лазера (или радиус преобразователя в случае использования емкостной установки).

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

Таблица 1

Результаты измерений импульсным методом (оптическая установка, измерения проведены по уровню)

Материал образца Толщина образца, мм С/7', м/с Cs, м/с а,дБ/м*

40X13 10,84 5999+1 3261±5 20+2

Д16Т 9,89 6420+1 311Ш0 50±4

пл. кварц 19,949 5959,510,2 3749,8±1,4 1,2±0,5

* измерения проведены на установке ИЗУ, на частоте 20 МГц.

Исследования показали, что алгоритм вычисления С^ по центрам тяжести переотраженных импульсов для ряда образцов «получает» меньшее значение скорости, чем другие алгоритмы. Так, для образца из алюминия Д16Т среднее значение групповой скорости для рассматриваемого алгоритма составило С/" =6412 м/с, в то время как значения, полученные по другим алгоритмам, были не ниже 6418 м/с. Данный факт говорит о необходимости более подробного рассмотрения алгоритма и условий его применения.

Во второй главе также приводятся полученные графики зависимостей С/, и а от частоты. Графики для образца из кварца (см. табл. 1) приведены на рис. 2.

Точки на образцах, в которых проводились измерения параметров С/, и а, выбирались следующим образом. Точка «О» - центр образца. Точки «1», «2», «3» - расположены по вершинам равностороннего треугольника с центром в точке «О» на расстоянии 15 мм от нее.

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

образца из кварца получили

/■ \ ас,;*''

А С,"'

п 'Р

= 0,006%, в то время как для образца из

алюминия Д16Т

С 'V

Ч. ср

= 0,041%, что говорит, в частности, о сравнительно

V V /

высокой структурной неоднородности последнего.

Рис. 2 Зависимость С; (а) и а (б) от частоты для образца из кварца

В третьей главе рассматривается резонансный режим измерений. При построении ИИК для резонансного метода в качестве основы используется установка ИЗУ-1. Для оцифровки сигналов предлагается использование цифрового анализатора спектра С5Р-827, а также предложено использование аналогового прибора СК4-59 с последующей оцифровкой спектрограммы цифровым осциллографом.

В главе приведены методики определения скорости распространения и коэффициента затухания у.з. колебаний, приведена методика оценки дополнительной дифракционной поправки при измерениях коэффициента затухания на частотах ниже 5 МГц. Также приведена методика определения коэффициента затухания и скорости распространения ультразвуковых колебаний при наличии помехи, зависящей от частоты.

Измерение скорости продольных колебаний в ИИК может быть произведено двумя способами. Первый способ основан на следующем выражении:

2 ¿{¡„-О

п-т

(И)

где/„,/,„ - частоты акустических спектральных линий с номерами п и т соответственно.

Также в работе предложена методика, основанная на следующем выражении:

_2#„

(12)

где/, - частота максимума акустической спектральной линии (АСЛ) с номером и. При определении С/, по формуле (12) методика измерений упрощается и снижается СКО результата измерения.

При малых значениях коэффициента затухания используется следующее соотношение:

20 лАГ

* = ' (13>

где АР - ширина акустической спектральной линии по уровню 0,7071. При измерении Д/7 на уровне ¡3 от максимума:

20.

а = -—1ё а

/ЗьтАц/ Р ьтку/

+ 1

(14)

где Ау/= лАР/Су

Для резонансных методов известно, что на частотах менее 5 МГц точностные характеристики снижаются вследствие значительного влияния дифракционных эффектов и искажения спектральной картины наложением «боковых» переотражений у.з. колебаний от ограничивающих поверхностей образца. Так, дифракционная расходимость проявляется в уменьшении амплитуды и увеличении ширины акустической спектральной линии, а также ее незначительном смещении в область низких частот. Переотражения приводят к появлению дополнительных спектральных линий. В работе приведена разработанная методика оценки дополнительной дифракционной поправки при измерении коэффициента затухания ре-

зонансным методом на частотах ниже 5 МГц, позволяющая повысить точность измерений при величине коэффициента затухания менее 5 дБ/м.

Для определения дополнительной дифракционной поправки были проведены эксперименты, на основе которых было получено следующее соотношение:

(15)

где А = 0,971 ±0,001 кГц; 6 = 3,26 1(Г3 ; # = 10 ± 0,1; с1-толщина образца; а -радиус преобразователей; к =со/С[_.

При этом коэффициент затухания рассчитывается по формуле:

а

= [20/(С, 1п Ю)]ф^£п -Д/-,-Д^ор), (16)

где Л/^,,- измеренная ширина акустической спектральной линии; Л/ч, - ширина спектральной линии, определенная теоретически при «нулевом» значении коэффициента затухания.

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

При измерении коэффициента затухания и скорости распространения у.з. колебаний при наличии помехи, зависящей от частоты, используется выражение, описывающее сигнал на выходе приемного емкостного преобразователя при наличии наводки. Анализ экспериментальных данных показал, что на частотах выше 20 МГц уровень помехи (параметр А) может изменяться в зависимости от частоты:

-[л

у (17)

(Ву)2 | 2ЛуВ [С05 0//-ф)-у2СО5 (у/+ф)} \72 \+у4 -2у2 со$2ц/ -2у2соъ2у/

где А и (р— амплитуда и фаза помехи соответственно; В = Кщ.

В настоящей работе предполагается, что параметр А связан с частотой соотношением:

А(П = СгГ2 + с2/ + с3, (18)

где с,, с2, с3 - некоторые константы.

Анализ выражения (17) показывает, что при Ву >2А и АР < 0,1-А/, методика измерения и расчетные соотношения для параметров а и О. не изменяются. Если А зависит от частоты, то предлагается несколько модифицированная методика. При этом выбираются АСЛ, для которых третий член в выражении

(17) имеет знак "+", определяется положение максимума АСЛ/тах и ее номер п, далее производится вычисление С/, по формуле (12). Для определения а предлагается следующая схема (рис. 3, /0 - частота, начиная с которой производится запись спектрограммы).

Для выбранной АСЛ определяются: максимум АСЛ £/„,; величина первого минимума слева от АСЛ {/_; величина первого минимума справа от АСЛ и частота/^ (слева от АСЛ), для которой уровень сигнала меньше £/„, на величину 0,3•(£/„, - -частота/, (справа от АСЛ), для которой уровень сигнала меньше 1/т на величину 0,3■(£/,„ - и.). Далее находится ширина АСЛ Д/7 =/+ -/_ и по формуле (16) производится вычисление коэффициента затухания. Проведенные эксперименты показали, что погрешность определения а по данной схеме не превышает 0,1-а дБ/м при а б (1 + 70) дБ/м.

При Д/" > 0,1 -Д/и большом уровне наводки форма сигнала в зависимости от уровня помехи и ее фазы может изменяться в значительной степени, при этом реализация выше описанной методики может оказаться затруднительной.

Была предложена следующая методика измерения коэффициента затухания и скорости распространения резонансным методом при наличии значительного уровня помехи (а > 70 дБ/м).

Рис. 3 Схема определения Л/-4 при наличии помехи.

Шаг 1. Определяется параметр В (17), зависящий от характеристик приемного тракта установки ИЗУ, на частотах, для которых помеха отсутствует. в = и„л\-у)1у, где ит - максимум выбранной АСЛ, у = ехр(-аг-г/), а - измеренное значение коэффициента затухания.

Шаг 2. Осуществляется цифровая запись сигнала (7Э(/) таким образом, чтобы на спектрограмме наблюдалось целое число периодов. Для сокращения времени вычислений делается сдвиг на некоторую частоту /0 начала записи спектрограммы.

Шаг 3. Определяется частотная зависимость уровня помехи. Для этого, в анализируемой области, строится квадратичный тренд и определяются коэффициенты, входящие в выражение (18).

Шаг 4. Производится частичная «компенсация» помехи и «строится» модифицированный сигнал гУ,э(/) = иэ(/)-А(/) = 1/э(/) - (с1[1+с2[+с1).

Шаг 5. Определяется параметр С/ либо по особым точкам (например, по положению максимумов сигнала или по пересечению нуля), либо автокорреляционным методом. С1 = с1/(пА/}, где п, А/- количество периодов и интервал частот между выбранными точками соответственно.

Если погрешность измерения Сд не должна превышать 0,1 %, то погрешность определения частот выбранных точек должна быть не более 200 Гц при толщине образцов менее 0,05 м.

Шаг 6. Производится определение параметра а. Для этого с учетом данных, полученных на шагах 1, 4, 5, предположения, что (ро = 0, у0 = (и^тах/В {{и1Э)тах - максимальный уровень сигнала £/)Э в анализируемой частотной области) и выражения (17) строится расчетная зависимость (/(/). Далее путем подбора значений параметров у, ср (при необходимости и параметров В, С), с2, с3) добиваются того, что бы сумма квадратов разностей экспериментальной и

расчетной зависимостей сигнала 5(5, (р, у) = > найденная,

/

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

201 л 20 &Ц а =--1еу; Аа=--, (19)

й с1уЫЪ В

где А11-погрешность измерения сигнала.

Данная методика позволяет определять а при у > 0,05 и А < 30Ву при погрешности менее 10%. Для более высокого уровня помехи погрешность измерений увеличивается.

При измерении скорости по «особым точкам» вначале выделяется область, в которой находится эта точка (примерно в середине области), затем для нее строится квадратичный тренд и по его параметрам (в[, в2, вз) определяется положение точки.

Ь2

При измерении по максимуму: /т; =--—; при измерении по пересече-

2Ь\

-Ь2± |Ь|-4Ь1Ь3

нию нуля: foi = ---, где i =1 - первая точка; i =2 - вторая точка.

2 Di

Для /0, «-» берется, если dU/df в искомой точке < 0 и «+», если dU/df в искомой точке > 0. После определения положения точек скорость рассчитывается по

формуле CL --, где Af=f„,, -fm2 (или Af=f0, -foi)', n - количество периодов

между выбранными точками.

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

В третьей главе также приводятся результаты экспериментальных исследований, проведенных с целью апробации разработанных методик и подтверждения функциональных характеристик ИИК. При проведении экспериментов было установлено, что при величине коэффициента затухания более 1 ООО дБ/м происходит эффект усиления полезного сигнала, обусловленный наличием помехи от задающего генератора анализатора спектра.

В четвертой главе приведены результаты анализа точностных характеристик ИИК с учетом предложенных алгоритмов и методик, а также его метрологические характеристики (табл. 2). Приведены структура и состав ПО ИИК, обоснован выбор программных и аппаратных средств, проведен анализ нестандартных ситуаций, способных повлиять на работу ИИК в автоматическом режиме.

Таблица 2

Основные метрологические характеристики ИИК

Величина Диапазон измерений Погрешность

Импульсный режим

Групповая скорость распространения С/'', м/с 2000-15000 <0,05%

Скорость распространения сдвигового импульса СЛ, м/с 1000-8000 <0,5%

Коэффициент затухания а, дБ/м 10-1000 <20%

Резонансный режим

Скорость распространения продольных волн С/ , м/с 2000-15000 <0,1%

Коэффициент затухания а, дБ/м 0,2 < а < 1 <20%

1 <а<5 <10%

5 < а < 200 <5%

200 < а <2000 <10%

Особое внимание в работе уделено вопросу метрологической аттестации программного обеспечения. Рассмотрены государственные и международные стандарты и руководства, регламентирующие различные аспекты создания и экс-

плуатации программного обеспечения средств измерений. На всех этапах разработки ПО НИК учитывались требования рассмотренных нормативных документов и требования к различной создаваемой сопроводительной документации.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

1. Разработаны алгоритмы измерения скорости распространения и коэффициента затухания у.з. колебаний импульсным методом. Предложена методика, основанная на автокорреляционной функции.

2. Разработаны алгоритмы и методики измерения скорости сдвиговых у.з. колебаний импульсным методом.

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

4. Разработаны алгоритмы определения фазовой скорости распространения у.з. колебаний для резонансного метода.

5. Предложена методика оценки дополнительной дифракционной поправки при измерении коэффициента затухания резонансным методом на частотах ниже 5 МГц.

6. Для резонансного режима предложена и программно реализована методика определения коэффициента затухания и скорости распространения у.з. колебаний при наличии помехи, зависящей от частоты.

7. Установлено, что при величине коэффициента затухания более 1000 дБ/м происходит эффект усиления полезного сигнала, обусловленный наличием помехи от задающего генератора анализатора спектра.

8. С целью подтверждения функциональных характеристик ИИК проведены экспериментальные исследования с использованием оптической и емкостной установок. Сделаны выводы о необходимости дальнейшего изучения условий применения метода определения групповой скорости по центрам тяжести импульсов.

9. Проведена оценка составляющих погрешностей для реализованных алгоритмов и методик измерений.

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

11. Сформирован набор программных библиотек для реализации различных численных методов, позволяющих создавать на их основе программное обеспечение средств измерений, удовлетворяющее требованиям, предъявляемым при метрологической аттестации.

12. Изучены и учтены требования нормативных документов и международных руководств по программному обеспечению средств измерений. Проведена работа по подготовке программного обеспечения к метрологической аттестации.

13. Разработанное алгоритмическое и программное обеспечение применяется в установках высшей точности, разрабатываемых в Дальневосточном филиале ФГУП «ВНИИФТРИ», которые в 2011-2013 годах будут переведены в ранг государственных эталонов.

ОСНОВНЫЕ ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Король, A.A. Определение дифракционных поправок по коэффициенту затухания ультразвуковых волн для резонансного режима измерений |Текст] / A.A. Король, А.И. Кондратьев, М.С. Жукова // Измерительная техника. - 2010, №4. - С. 29-43.

2. Король, A.A. Динамика дисперсионных характеристик акустических параметров полимеризуюидихся составов (Текст) / A.A. Король,

B.И. Римлянд, В.Н. Старикова, А.И. Кондратьев // Вестник Тихоокеанского государственного университета № 3 (14), 2009. - С. 23-32.

3. Korol, A.A. Determination of the diffraction corrections to the attenuation factor of ultrasonic waves for the resonance method of measurement [Текст] / A.A. Korol, A.I. Kondratiev, M.S. Zhukova // Measurement Techniques - Volume 53, Issue4 (2010).-P. 431-438.

4. Король, A.A. Программное обеспечение установки для комплексного измерения акустических величин [Текст] / A.A. Король, А.И. Кондратьев // Межрегиональная научно-практическая конференция «Информационные и коммуникационные технологии в образовании и научной деятельности» (21-23 мая 2008 года), г. Хабарвоск: материалы конференции / под науч. ред. А.И. Мазура. - Хабаровск: изд-во Тихоокеан. гос. ун-та, 2008. - 373 с. -

C. 298-303.

5. Король, A.A. Подходы к обеспечению необходимой точности при измерениях акустических величин твердых сред [Текст] / A.A. Король, А.И. Кондратьев // Высокие технологии, фундаментальные и прикладные исследования, образование. Т. 13: Сборник трудов Пятой международной научно-практической конференции «Исследование, разработка и применение высоких технологий в промышленности». СПб.: Изд-во Политехи, ун-та, 2008. - 491 с. -С. 202-204.

6. Король, A.A. Автоматизация обработки сигналов с использованием цифровых осциллографов Lecroy [Текст] / A.A. Король // XXXIV Дальневосточная математическая школа-семинар имени академика Е.В. Золотова «Фундаментальные проблемы математики и информационных наук». - Хабаровск: Изд-во Тихоокеан. гос. ун-та, 2009 - 198 с. - С. 171-172.

7. Korol A.A. Influence of heat treatment of alloy D16 on formation acoustic spectral lines [Текст] / A.A. Korol, A.I. Kondratiev, M.S. Zhukova // Modern materials and technologies 2009: International Xth Russian-Chinese Symposium. Proceedings. - Khabarovsk: Pacific National University, 2009. - P. 461-464.

8. Король, A.A. Влияние силы прижатия электрода к образцу на форму у.з. импульсов [Текст] / A.A. Король, А.И. Кондратьев, М.С. Жукова // Физика геосфер: шестой Всероссийский симпозиум, 7-11 сентября 2009 г., Владивосток, Россия: мат-лы докл. - Владивосток: ТОЙ ДВО РАН, 2009. - 357 с. -С. 186-190.

9. Король A.A., Ли В.Н., Кондратьев А.И., Титов Е.А. Программный комплекс расчета степени локального разуплотнения контактного провода. Свидетельство о государственной регистрации программы для ЭВМ №2008615538,19.11.2008 г.

10. Король A.A., Луговой В.А., Кондратьев А.И., Базылев П.В., Ни-гай В.П. Программно-измерительный комплекс для расчета параметров распространения различных типов акустических волн в твердых средах. Свидетельство о государственной регистрации программы для ЭВМ №2011610483, 11.01.2011 г.

Король Андрей Александрович

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

Автореферат диссертации на соискание ученой степени кандидата технических наук

Подписано в печать 25.04.2011. Формат 60x84Vif„ Гарнитура Times New Roman. Уч.-изд. л. 1,3. Усл. печ. л. 1,2. Зак. 180. Тираж 100 экз.

Издательство ДВГУПС 680021, г. Хабаровск, ул. Серышева, 47.

Список сокращенных наименований.

ВВЕДЕНИЕ.

1 МЕТОДЫ ИЗМЕРЕНИЯ АКУСТИЧЕСКИХ ВЕЛИЧИН.

1.1 Методы измерения акустических величин. Источники возникновения погрешностей.

1.2 Методы измерения скорости распространения ультразвуковых колебаний.

1.3 Методы измерения коэффициента затухания ультразвуковых колебаний.

1.4 Методы возбуждения и приема ультразвуковых колебаний.

2 РАЗРАБОТКА АЛГОРИТМОВ И МЕТОДИК ИЗМЕРЕНИЯ СКОРОСТЕЙ РАСПРОСТРАНЕНИЯ И КОЭФФИЦИЕНТА ЗАТУХАНИЯ УЛЬТРАЗВУКОВЫХ КОЛЕБАНИЙ ИМПУЛЬСНЫМ МЕТОДОМ.

2.1 Измерение скорости распространения продольных ультразвуковых колебаний.

2.2 Измерение скорости распространения сдвиговых ультразвуковых колебаний.

2.3 Измерение скорости распространения и коэффициента затухания ультразвуковых колебаний в зависимости от частоты.

2.4 Проведение экспериментальных исследований.

3 РАЗРАБОТКА АЛГОРИТМОВ И МЕТОДИК ИЗМЕРЕНИЯ ФАЗОВОЙ СКОРОСТИ РАСПРОСТРАНЕНИЯ И КОЭФФИЦИЕНТА ЗАТУХАНИЯ УЛЬТРАЗВУКОВЫХ КОЛЕБАНИЙ РЕЗОНАНСНЫМ МЕТОДОМ.

3.1 Измерение скорости распространения продольных ультразвуковых колебаний.

3.2 Измерение коэффициента затухания ультразвуковых колебаний.

3.3 Методика оценки дополнительной дифракционной поправки при измерении коэффициента затухания резонансным методом на частотах ниже 5 МГц.

3.4 Методика определения коэффициента затухания и скорости распространения ультразвуковых колебаний при наличии помехи, зависящей от частоты.

3.5 Проведение экспериментальных исследований.

4 ИНФОРМАЦИОННО-ИЗМЕРИТЕЛЬНЫЙ КОМПЛЕКС.

4.1 Анализ точностных характеристик.

4.2 Основные метрологические характеристики информационно-измерительного комплекса.

4.3 Структура и состав программного обеспечения комплекса.

4.4 Определение нестандартных ситуаций в автоматическом режиме.

4.5 Выбор средств для реализации программного обеспечения.

4.6 Сопряжение с цифровым осциллографом ЬеСгоу \VS422.

4.7 Метрологическая аттестация программного обеспечения комплекса.

4.8 Качество программного обеспечения комплекса.

В настоящее время теоретические и экспериментальные исследования во многих областях науки проводятся с использованием современных информационных технологий. Применение цифровой техники позволяет значительно повысить качественные показатели измерений, такие как точность и скорость обработки данных. Другой важный аспект внедрения информационных технологий — автоматизация — позволяет упростить проведение трудоемких экспериментов, снизить влияние человеческого фактора на результаты, а также снизить сопутствующие материальные затраты [1].

Автоматизация измерений подразделяется на полную и частичную [1, 2]. Под полной автоматизацией понимается такой процесс измерений, при котором все действия — от получения первичной информации до получения конечных результатов и их анализа - происходят без участия человека. Частичная автоматизация подразумевает, что одним из звеньев в цепи получения измерительной информации является оператор, на которого могут возлагаться различные функции.

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

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

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

В рамках настоящей работы также предложен ряд методик измерений. Для импульсного режима предложены методики нахождения скорости сдвиговой составляющей импульса при наличии помехи, методика определения частотной зависимости скорости распространения и коэффициента затухания у.з. колебаний. Для резонансного режима предложена методика определения дополнительной дифракционной поправки при нахождении коэффициента затухания, а также методика определения коэффициента затухания и скорости ультразвуковых колебаний при наличии помехи, зависящей от частоты (при значениях коэффициента затухания до 2000 дБ/м). Предложенные методики позволяют повысить точность, а также расширить диапазон измерений соответствующих физических величин.

Целью настоящей работы является разработка алгоритмического и программного обеспечения измерительного комплекса для определения параметров распространения различных типов акустических волн в твердых телах.

Основные задачи настоящей работы:

1. Анализ методов измерения акустических величин и методов возбуждения и приема у.з. колебаний.

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

3. Разработка алгоритмов измерения фазовой скорости распространения и коэффициента затухания у.з. колебаний при измерениях резонансным методом.

4. Разработка программных библиотек численной обработки данных и взаимодействия с пользователем для построения программного обеспечения с учетом требований, предъявляемых к программным комплексам средств измерения в процессе их метрологической аттестации.

5. Разработка программно обеспечения измерительного комплекса, реализующего предложенные алгоритмы.

6. Исследование метрологических характеристик информационно-измерительного комплекса с учетом разработанных алгоритмов и программного обеспечения.

Методы исследования. При выполнении работы использовались: теория и методы неразрушающего контроля, теория дифференциальных уравнений, теория и математический аппарат преобразования Фурье, численные методы интегрирования и дифференцирования, теория объектно-ориентированного программирования; системное программирование для операционных систем (ОС) семейства Microsoft Windows.

Научная новизна работы заключается в следующем:

I. Предложены и реализованы алгоритмы измерения скоростей распространения и коэффициента затухания продольных и сдвиговых у.з. волн импульсным и резонансным методами.

2. Разработаны алгоритмы определения времени прихода сдвиговой компоненты сигнала при наличии помехи, основанные на аппроксимации фронта сдвигового импульса различными функциями.

3. Разработана методика оценки дополнительной дифракционной поправки при измерении коэффициента затухания резонансным методом на частотах ниже 5 МГц для твердых сред.

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

5. Установлено, что при величине коэффициента затухания более 1000 дБ/м происходит эффект усиления полезного сигнала, обусловленный наличием помехи от задающего генератора анализатора спектра.

Практическая значимость работы:

1. На основе предложенных алгоритмов разработано программное обеспечение измерительного комплекса для определения параметров распространения различных типов акустических волн в твердых средах, обеспечивающее снижение влияния субъективного фактора и, как следствие, снижение среднеквадратического отклонения (СКО) результатов измерения.

2. Сформирован набор программных библиотек для реализации различных численных методов и взаимодействия с пользователем, позволяющих создавать на их основе программное обеспечение средств измерений, удовлетворяющее требованиям, предъявляемым при метрологической аттестации.

3. Предложенная методика определения коэффициента затухания для резонансного режима позволяет повысить верхнюю границу диапазона измерений до значения 2000 дБ/м.

Применение результатов работы. Разработанное алгоритмическое и программное обеспечение применяется в установках высшей точности, разрабатываемых в Дальневосточном филиале ФГУП «ВНИИФТРИ», которые в

2011-2013 годах будут переведены в ранг государственных эталонов. 8

Внедрение результатов диссертационной работы подтверждено соответствующим актом (приложение 2).

На основе разработанных алгоритмов и программных библиотек на базе государственного образовательного учреждения высшего профессионального образования «Дальневосточный государственный университет путей сообщения» был реализован и внедрен программный комплекс расчета степени локального разупрочнения контактного провода.

Основные положения, выносимые на защиту:

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

2. Алгоритмы определения скорости распространения сдвиговых у.з. колебаний импульсным методом.

3. Алгоритмы и методики определения фазовой скорости распространения и коэффициента затухания продольных у.з. колебаний резонансным методом.

4. Структура и состав программного обеспечения измерительного комплекса для определения параметров распространения различных гипов акустических волн в твердых средах.

Апробация работы. Основные положения и отдельные результаты работы докладывались и обсуждались: на межрегиональной научно-практической конференции «Информационные и коммуникационные технологии в образовании и научной деятельности» (г. Хабаровск, 2008); на пятой международной научно-практической конференции «Исследование, разработка и применение высоких технологий в промышленности» (г. Санкт-Петербург, 2008); на V Всероссийской конференции «Механика микронеоднородных материалов и разрушение» (г. Екатеринбург, 2008); на XXXIV Дальневосточной математической школе-семинаре им. академика Е.В.Золотова (г. Хабаровск, 2009); на пятой Всероссийской конференции «Математическое моделирование и краевые задачи» ММР-2008 (г. Самара,

2008); на международной научной конференции «Оптика кристаллов и нано9 структур» (г. Хабаровск, 2008); на VI Всероссийском симпозиуме «Физика геосфер» (г. Владивосток, 2009).

Публикации. Основное содержание диссертационной работы отражено в 11 публикациях, в том числе 2 из них - в изданиях из перечня ВАК, 6 -на научно-практических конференциях (из них две - на международных), два свидетельства о государственной регистрации программ для ЭВМ.

Личный вклад автора. Личный вклад заключается в разработке и анализе предложенных алгоритмов и их программной реализации, в разработке программно-аппаратных средств автоматизации измерительных процессов. Также при участии автора были проведены экспериментальные исследования.

Структура и объем работы. Диссертационная работа представлена введением, четырьмя главами, заключением, списком использованных литературных источников и приложениями. Работа изложена на 140 страницах машинописного текста и содержит 6 таблиц, 38 рисунков, 137 наименований литературных источников, 3 приложения.

Выводы по главе.

1) Проведена оценка составляющих погрешностей для реализованных методов измерений.

2) Приведены основные метрологические характеристики ИИК с учетом предложенных алгоритмов и методик.

3) Приведена структура и состав ПО ИИК.

4) Проведен анализ нестандартных ситуаций, которые могут повлиять на работу ИИК в автоматическом режиме.

5) Приведены особенности сопряжения управляющей ЭВМ с цифровым осциллографом ЬеСгоу 422.

6) Изучены и учтены требования нормативных документов и международных руководств по программному обеспечению средств измерений. Проведена работа по подготовке ПО к метрологической аттестации, предложены инструменты для аттестации методом «черного ящика».

ЗАКЛЮЧЕНИЕ

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

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

2) Разработаны алгоритмы измерения скорости распространения и коэффициента затухания у.з. колебаний импульсным методом. Предложена методика, основанная на автокорреляционной функции.

3) Разработаны алгоритмы и методики измерения скорости сдвиговых у.з. колебаний импульсным методом.

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

5) Разработаны алгоритмы определения фазовой скорости распространения у.з. колебаний для резонансного метода.

6) Предложена методика оценки дополнительной дифракционной поправки при измерении коэффициента затухания резонансным методом на частотах ниже 5 МГц.

7) Для резонансного режима предложена и программно реализована методика определения коэффициента затухания и скорости распространения у.з. колебаний при наличии помехи, зависящей от частоты.

8) Установлено, что при величине коэффициента затухания более 1 ООО дБ/м происходит эффект усиления полезного сигнала, обусловленный наличием помехи от задающего генератора анализатора спектра.

9) С целью подтверждения функциональных характеристик ИИК проведены экспериментальные исследования с использованием оптической и емкостной установок. Сделаны выводы о необходимости дальнейшего изучения условий применения метода определения групповой скорости по центрам тяжести импульсов.

10) Проведена оценка составляющих погрешностей для реализованных алгоритмов и методик измерений.

11) Проведен анализ нестандартных ситуаций, которые могут повлиять на работу ИИК в автоматическом режиме. Программно реализованы механизмы, позволяющие снизить вероятность возникновения грубых промахов.

12) Разработано вспомогательное программное обеспечение, необходимое для автоматизации измерительных процессов с использованием цифровых осциллографов ЬеСгоу \¥ауе8шТег.

13) Сформирован набор программных библиотек для реализации различных численных методов и взаимодействия с пользователем, позволяющих создавать на их основе программное обеспечение средств измерений, удовлетворяющее требованиям, предъявляемым при метрологической аттестации.

14) Изучены и учтены требования нормативных документов и международных руководств по программному обеспечению средств измерений. Проведена работа по подготовке программного обеспечения к метрологической аттестации, предложены инструменты для аттестации методом «черного ящика».

15) Разработанное алгоритмическое и программное обеспечение применяется в установках высшей точности, разрабатываемых в Дальневосточном филиале ФГУП «ВНИИФТРИ», которые в 2011-2013 годах будут переведены в ранг государственных эталонов.

1. Перцовский М.И. Лабораторная автоматизация: организация современных приборных комплексов, систем проведения экспериментов и испытаний // RM MAGAZINE. 2005, № 6. - С. 46-52.

2. Дубовой Н.Д., Осокин В.И., Поротов В.Н. Автоматизация измерений и контроля электрических и неэлектрческих величин. — М: Изд-во стандартов, 1987-328 с.

3. Кондратьев А.И. Прецизионные методы и средства измерения акустических величин твердых сред. Хабаровск: изд-во ДВГУПС, 2006 - 152 е.: ил.

4. Труэл Р., Эльбаум Ч., Чик Б. Ультразвуковые методы в физике твердого тела. М.: Мир, 1978 - 544 с.

5. Ермолов И.Н. Теория и практика ультразвукового контроля. М.: Машиностроение, 1981 -240 с.

6. Бреховских Л.М. Волны в слоистых средах. М.: Наука, 1973 - 273с.

7. Викторов И.А. Звуковые поверхностные волны в твердых телах. -М.: Наука, 1981 -287 с.

8. Физическая акустика. Т. 1. Методы и приборы для ультразвуковых исследований // Под ред. У. Мазона. — М.: Мир, 1972 — 248 с.

9. Чабанов В.Е. Лазерный ультразвуковой контроль материалов. Л.: Ленинградский университет, 1986 - 242 с.

10. Hitchins D.A., Demhurst R., Palmer S.R., Scrubu C.B. Laser generation as a standard acoustic source in metals // Appl. Phys. Letters, 1981 №9., p. 677679.

11. Бондаренко A.H., Дробот Ю.Б., Кондратьев А.И. Прецизионные акустические измерения оптическими и емкостными методами. — Владивосток: ТОЙ ДВО АН СССР, 1990 240 с.

12. Кондратьев А.И. Исследование бесконтактных методов возбуждения ультразвуковых колебаний // Автореф. дис канд. физ.- мат. наук. Владивосток, 1983 - 22 с.

13. Пябус Г.В., Мельканович А.Ф., Кукшулей Л.М. Установка для измерения коэффициента затухания ультразвука в твердых телах // Дефектоскопия. 1987, №2. - С. 57-62.

14. Dewhurst R.J., Hutchins D.A., Palmer S.B. Quantitative measurements of laser generated acoustic waveforms // J. Appl. Phys. 1982. V.53. N 6. - p. 4064-4071.

15. Королев M.B., Карпельсон A.E. Широкополосные ультразвуковые преобразователи. M.: Машиностроение, 1982 -157 с.

16. Кажис Р.И. Ультразвуковые контрольно-измерительные системы. -Вильнюс.: Мокслас, 1986-216 с.

17. Ультразвуковые преобразователи для неразрушающего контроля / Под ред. И.Н. Ермолова. М.: Машиностроение, 1986 - 277 с.

18. Бражников Н.И. Ультразвуковая фазометрия. М.: Энергия, 1968.272 с.

19. Rose I. L., Meger P. A. Ultrasonic Signal-Processing Concepts for Measurement the Thickness of Thin Lasers // Mater evaluation. 1974. 32. N 12. -p. 249-255.

20. Цветянский B.Jl. О прохождении ультразвуковых колебаний через контактный слой при акустических исследованиях твердых тел // Акуст. Журнал-1981, №4.-С. 610-615.

21. Архипов В.И., Кондратьев А.И. Исследование прохождения ультразвукового импульса через слой жидкости // Дефектоскопия. — 1994, № 1. — С. 21-25.

22. Алешин Н.П., Белый В.Е., Вопилкин А.Х., Ермолов И.Н., А.К. Гурвич Методы акустического контроля металлов — М.: Машиностроение, 1989-455 с.

23. Архипов В.И., Бондаренко A.H., Дробот Ю.Б., Троценко В.П. Образцовая лазерная установка для аттестации акустических мер по скорости ультразвука // Измерительная техника — 1984, №2. — С. 60-62.

24. Недбай А.И., Федоровский Г.Д. Импульсно-ингерференционный метод измерения скорости ультразвука // Дефектоскопия — 1985, №1. — С. 4954.

25. Hsu T.S., Marston J.B. Measurement of velocity in the spin-glass CuMn // J. Appl. Phys, 1987, N5. p. 2074-2077.

26. Тулуз Д., Лоней К. Автоматизированная система для измерения ослабления ультразвука и относительных изменений его скорости // Приборы для научных исследований 1988, №3. - С. 98-102.

27. Wang Y., Wei М. An apparaturs for Ultrasonic Velocity and Attenuation Measurements // Appl. Acoust 1989, N2. - p. 1-5.

28. Баранов B.M. Акустические измерения в ядерной энергетике. М.: Энергоатомиздат, 1990. - 320 с.

29. Жилинскас Р.П., Илгунас В.И., Яронис Э.П. Цифровой интерферометр для измерения скорости ультразвука // Акуст. Журнал 1967. Т. 10. №1 -С. 118-119.

30. Kinra V.K., Dayal V. A new Technique for Ultrasonic // Nondestructive Evaluation on Thin Specimens // Experimental Mechanics 1988. 28. N3. - p. 288-297.

31. Willams A.O. The Piston Source at High Frequencies // J. Acoust. Soc. Amer- 1951. V. 23. N1 p. 1-6.

32. Меркулова B.M. Расчет характеристик направленности поршневогоизлучателя в импульсном режиме // Дефектоскопия 1967, №1. — С. 7-12.117

33. Ермолов И.Н. Методы расчета акустического тракта ультразвукового дефектоскопа. 1. Акустическое поле нормального контактного искателя. 1967. №3. С. 41-50.

34. Гитис М.Б. О дифракционных поправках в импульсном режиме // Акуст. Журнал 1972, Т. 18. №1. - С. 42-48.

35. Гитис М.Б., Серегин Е.И. О дифракционных поправках при ультразвуковых измерениях в импульсном режиме // Дефектоскопия — 1978. №3. -С. 90-93.

36. Серегин Е.И., Зайцев Б.Л., Коржик И.Г. Способ учета дифракционных явлений при измерении скорости распространения ультразвука в металлах // Тез. докл. IX Всес конф. «Неразрушающие физические методы».: Минск, 1981. А-127. С. 234-236.

37. Архипов В.И., Бондаренко А.Н., Кондратьев А.И. Оптический метод измерения скорости сдвиговых волн // Измерительная техника 1984. № З.-С. 27-28.

38. Бондаренко А.Н., Дробот Ю.Б., Вологдин В.К. Оптическое устройство для измерения групповой скорости ультразвука // Измерительная техника- 1980. № З.-С. 68-69.

39. Бондаренко А.Н., Дробот Ю.Б., Кондратьев А.И. Ultrasound velocity laser meter // Poster session. Tent world conf. nondestructive Testing. 26 th August. 1989.

40. Кондратьев А.И., Криницин Ю. M. Измерение параметров распространения акустических волн резонансным методом на межпиковых частотах // Автометрия 2000. № 1. - С.115.

41. Зонов И.В. Определение дальней зоны преобразователей работающих в импульсном режиме // Дефектоскопия 1981. №7. - С. 70-77.

42. Бондаренко А.Н., Кондратьев А.И. Измерение дисперсии скорости и затухания упругих волн // Акуст. журнал — 1981. т. 27. № 1. — С. 51 -57.

43. Физическая акустика. / Под ред. У. Мэзона Т.5 // М.: Мир, 1973.332 с.

44. Совершенствование установки высшей точности для воспроизведения единицы скорости распространения продольных у.з. волн в твердых средах / Научн. руководитель Архипов В.И. // Отчет НПО «Дальстандарт». гос per. №76047269.: Хабаровск, 1989. 75 с.

45. Базылев П.В., Бондаренко А.Н., Луговой В.А. Лазерная установка для измерения скорости распространения поверхностных волн Рэлея. // Дефектоскопия- 1990. №Ю С. 91-93.

46. Гакенхеймер Д. Численные результаты в задаче Лэмба о действии сосредоточенной динамической нагрузки // Прикладная механика — 1970. №2. -С. 2.

47. Колесников А.Е. Ультразвуковые измерения // М.: Изд-во стандартов, 1974.-288 с.

48. McSkimin Н. J. Notes and References for the Measurement of Elastic Module Means of Ultrasonic Waves // J. Acoust. Soc. Amer 1961. v.33. N.5. - p. 606-616.

49. Kinra V.K., Daual V. A new Technique for Ultrasonic // Nondestructive Evaluation on Thin Specimens // Experimental Mechanics 1988. V. 28. N 3 - p. 288-297.

50. Batra N. K., Delsanto P.P. Simultaneous measurements of ultrasonic phase velocity and attenuation in solids // Rev Progr. Quant. Nondistr. Eval. Vol 6 A lat Haef Proc. 13 th Amtr Rev. Progr., 1987. p. 491-499.

51. Sorhau S.K., Kline R.A., Migonogna R. Phase and group velocity considerations for dynamic modulus measurement in anisotropic media // Ultrasonics 1992. V. 30. N6. - p. 373-382.

52. Ditchi Т., Alquie C., Lewiner J. Broadband determination of ultrasonic attenuation and phase velocity in insulating materials // J Acoust. Soc. Amer. -1993. V. 94. N6.-p. 3061-3066.

53. Болтарь K.O., Мансфельд Г.Д. Возбуждение ультразвуковых импульсов в твердых телах // ПТЭ — 1977. №1. —С. 128-131.

54. Reynolds W.N. The analysis of ultrasonic wave attenuation spectra in metals // 3-rd European Conference on Nondestructive Testing. Florence, 1984. -p.355-361.

55. Приходько B.H. Ультразвуковой контроль межкристаллитной коррозии аппаратуры в производственных условиях. — М.: ГОСИНТИ, 1964. -244 с.

56. Зайцев Б.Л. Серегин Е.И. Методы и средства аттестации стандартных образцов, используемых в ультразвуковом неразрушающем контроле // Измерительная техника 1988. №10. — С. 58-60.

57. Воробьев В.В., Зубков Л.А. Затухание ультразвука и рассеяние света в гетерогенных фотохромных стеклах // Акуст. журнал 1990. Т. 36. №2. -С. 20-26.

58. Breseale М.А., Cantrell John A., Heuman J.S. Ultrasonic wave velocity and attenuation measurement // Meth. ExP. Phys. V. 19. N.Y. e.a., 1981. p. 67-70.

59. Wong P.H., Garmi C.W. Ultrasonic velocity and attenuation measurement WHR computer controling pioni sensetive detection technique // Rev. Sci. Ingening- 1986. V. 50. N 12. p. 3085-3090.

60. Canela G., Tadei M. Ultrasonic attenuation in direct contact, in the far field and at high frequencies in thin samples / 3-rd European conference on Nondestructive Testing // Florence, 1984. p. 234-245.

61. Generazio E.R. The Role of the Reflection Coefficient in Precision Measurement of Ultrasonic Attenuation // Material Evaluation 1985. V. 43. N 7. - p. 995-1004.

62. Tang H.M., Toksoz M.,N., Tarif P., Wikens R.H. A metod for measuring acoustic wave attenuation in the laboratory // J. Acoust. Soc. Ainer. — 1988. V.83. N 2. p. 453-462.

63. Кукшулей Л.М., Мельканович А.Ф. Многократное отражение ультразвукового импульса в плоскопараллельном образце / Тез. Докл. IX Всес конф. «Неразрушающие физические методы и средства контроля» // Минск, 1981. Секция А. С. 28-30.

64. Ермолов И.Н. Методики измерения затухания продольных волн // Дефектоскопия — 1995. №7 С. 3-13.

65. Меркулов Л.Г., Токарев В.А. Физические основы спектрального метода измерения затухания ультразвуковых волн в материалах // Дефектоскопия. 1970. №4. - С. 3-11.

66. Кононенко B.C. Исследование погрешности измерения в ультразвуковом резонаторе, связанной со спектром его собственных частот // Акуст. журнал 1984. Т. 30. №6. - С. 785-789.

67. Кононенко B.C. Прецизионный метод измерения коэффициента поглощения ультразвука в жидкостях на частотах 0,1 20 МГц // Акуст. журнал - 1987. Т. 33. №4. - С. 688-694.

68. Кондратьев А.И., Король A.A., Жукова М.С. Определение дифракционных поправок по коэффициенту затухания ультразвуковых волн для резонансного режима измерений. // Измерительная техника 2010, №4 — С. 2439.

69. Бондаренко А.Н., Кондратьев А.И., Луговой В.А. Тонкопленочный самоустанавливающийся емкостный преобразователь ультразвуковых сигналов // Приборы и техника эксперимента. — 1988, № 2. С. 197.

70. Гусаков С. А., Кондратьев А. И. Образцовая установка для комплексного измерения акустических параметров материалов // Измерительная техника 1989, № 7. - С. 50.

71. Кривошеев И. А., Кондратьев А. И. Исследование работы емкостного преобразователя в низкочастотном диапазоне // Дефектоскопия 1989. №7.-С. 13.

72. Bass R. Diffraction effects in the ultrasonic field of a piston source // J. Acoust. Soc. America. 1958. V. 36. № 7. - p. 602.

73. Кондратьев А. И. Прецизионные измерения скорости и затухания ультразвука в твердых телах // Акустический журнал — 1990. Т. 36. № 3. С. 470.

74. Бункин Ф.В., Комиссаров В.М. Оптическое возбуждение звуковых волн (обзор) // Акуст. журнал 1973. т. 19. № 3. - С. 305-320

75. Яворский Б.М., Детлаф А.А. Справочник по физике. — М.: Наука, 1968.-558 с.

76. Новацкий В. Вопросы термоупругости. М.: Изд-во АН СССР, 1962-380 с.

77. Анисимов С.И., Имас Я.Д., Романов Г.С., Ходыко Ю.В. Действие излучений большой мощности на металы М.: Наука, 1970. - 272 с.

78. Агеев В.П. Исследование механического действия импульсного излучения С02-лазера на твердые мишени в газовой среде // Квантовая электроника- 1977. т. 4. №3 С. 310-319.

79. Броуд Г. Расчеты взрывов на ЭВМ. Гидродинамика взрывов. — М.: Мир. 1976.-270 с.

80. Шкарлет Ю.М. Бесконтактные методы ультразвукового контроля. — М.: Машиностроение, 1974. — 56 с.

81. Гитис Н.Б., Добромыслов В.М. Определение некоторых параметров датчиков ультразвуковых дефектоскопов // Дефектоскопия 1971, №8. - С. 44-49.

82. Legros D., Lewinder I., Biguard P. Geniration of Ultrasound by a dielectric Transduced // J. Acoust. Soc. Amer 1972. V. 52. N 1. - p. 196-198.

83. Shourbagy Salan A.M., Okeda Shiro, Matuttis Hans-Georg J. Acoustic of sound propagation in granular materials in one, two and three dimensions // Phys. Soc. Japan 2008, 77 №3 — пер. с яп., англ.

84. Petculesau P., Prodan G., Zugan R. An automated algorithm for simultaneously determining ultrasonic velocity and attenuation // Rom. J. Physics -2005, № 7-8, -C. 691-700.

85. Terao Michihito. Measurement techniques of acoustic properties of materials by using impedance tubes // Nihon onkyo gakkaishi, Acoustic Soc. Japan -2007, 63 №10 p. 606-611, пер. с яп.

86. Cantrel J.H.Jr. Breseall M.A. Capacative driver ting finite amplitude ultrasonic waves in Solids // Abstr. PaP. 7 th Int. SymP. Nonlin Acoust. Bracksburg 1976. V. l.-P. 92-94.

87. Бондаренко A.H., Маслов Б.Я., Рудая Б.Б., Троценко В.П. Оптическая установка для измерения сверхмалых акустических колебаний // ПТЭ -1975. №6. -С. 211-213.

88. Троценко В.П. Разработка и исследование оптических методов и средств измерения малых акустических сигналов на основе стабильных газовых лазеров // Автореф. дис канд. физ.-мат. наук. Ленинград: 1978 — 21 с.

89. Исследование методов генерации и измерения малых ультразвуковых перемещений поверхности твердых тел (диапазон частот 1 кГц //3014. 8

90. МГц. диапазон амплитуд 10" // 5ТО" м) / Бондаренко А.Н., Троценко В.П. (отв. исполнитель). Кондратьев А.И., Криницин Ю.М., Луговой В.А. // Отчет НПО «Дальстандарт». roc per. №76047269. Хабаровск. 1979 - 131 с.

91. Айвазян С.А. Прикладная статистика. Основы эконометрики. Том 2. — М.: Юнити-Дана, 2001. —- 432 с.

92. А.Г.Зюко Д. Д.Кловский, В.И.Коржик, М.В.Назаров. Теория электрической связи: Учебник для вузов, под ред. Д.Д.Кловского. — М.: Радио и связь, 1999. 432 с: 204 ил.

93. Дрейпер Н., Смит Г. Прикладной регрессионный анализ (книга 2). -М.: Финансы и статистика, 1986. — 351 с.

94. Римлянд В.И., Старикова B.H., Кондратьев А.И., Король A.A. Динамика дисперсионных характеристик акустических параметров полимери-зующихся составов. Вестник ТОГУ — 2009, № 3 (14) — С. 23-32 .

95. Бахвалов, Н. С. Численные методы / Н. С. Бахвалов, Н. П. Жидков, Г. Н. Кобельков. -М.: БИНОМ. Лаб. знаний, 2003 632 с.

96. Вержбицкий, В. М. Численные методы. Математический- анализ и обыкновенные дифференциальные уравнения. — М.: Высш.шк., 2001 — 383 с.

97. Волков, Е. А. Численные методы. СПб.: Лань, 2004 — 248 с.

98. Э. Таненбаум. Архитектура компьютера.5-е изд. СПб.: Питер, 2007. - 844 с.

99. Эхтер Ш., Роберте Дж. Многоядерное программирование. СПб.: Питер, 2010-316 е.: ил.

100. Richard Gerber, Andrew Binstock. Programming with Hyper-Threading Technology: How to Write Multithreaded Software for Intel IA-32 Processors Intel Press, 2004.

101. Richter J., Nasarre Ch. Windows via C/C++ Microsoft Press, 2007.

102. Дж. Рихтер. Windows для профессионалов СПб.: Питер, 2000752 с.

103. Richard Gerber, Aart J.С. Bik, Kevin В. Smith, Xinmin Tian. The Software Optimization Cookbook, Second Edition Intel Press, 2006.

104. Осциллографы LECROY WaveSurfer 422/424/432/434/452/454 Руководство по эксплуатации — M.: 2007 148 c.

105. Дедюхин А.А. Основные характеристики современных осциллографов. ЗАО "ПриСТ" 2004 эл. ресурс., режим доступа: http://www.prist.ru/info.php/articles/generalchmodernoscilloscopes.htm

106. Гмурман В. Е. Теория вероятностей и математическая статистика М.: Высш.шк., 2003 - 479 с.

107. ГОСТ Р ИСО/МЭК 7498-1-99. Информационная технология. Взаимосвязь открытых систем. Базовая эталонная модель. Часть 1. Базовая модель.

108. МИ 2083-90 ГСИ. Измерения косвенные. Определение результатов измерений и оценивание их погрешностей.

109. ГОСТ P 8.596-2002 ГСИ. Метрологическое обеспечение измерительных систем. Основные положения.

110. ГОСТ Р 8.654-2009 Требования к программному обеспечению средств измерений. Основные положения.

111. МИ 2891-2004 ГСИ. Общие требования к программному обеспечению средств измерений.

112. Ли В.Н., Кондратьев А.И., Король А.А., Титов Е.А. Программный комплекс расчета степени локального разуплотнения контактного провода -свидетельство государствершой регистрации программы для ЭВМ №2008615538, 2008.

113. Randy Abernethy, Randy Morin, Jesus Chahin. COM/Dcom Unleashed: SAMS, USA, 1999 667 p.

114. Nathan Wallace. COM/DCOM Blue Book: Fast Paced beaming: Cor-iolis Group Books, USA, 1999 747 p.

115. Бобровский С. EntireX технология ближайшего будущего. — PCWeek.ru Электронный ресурс. - М.: 1998. - режим доступа: http://www.pcweek.ru/themes/detail.php?ID=46806

116. Marco Cantu. Delphi 2010 Handbook: Self-Publishing Lulu.com, 2010. -318 p.

117. Marco Cantu. Mastering Borland Delphi 2005: Self-Publishing Lu-lu.coin, 2005.-345 p.

118. Фрэнк Хэйес. Distributed Component Object Model. перев. с англ. // COMPUTER WORLD РОССИЯ - M.: 1999 (№ 24).

119. Ямпольский Д.А. К вопросу о метрологической аттестации программного обеспечения вагонных весов. // Приборы+Автоматизация 2008 №1(91).

120. Кудеяров Ю.А. Аттестация программного обеспечения средств измерений. Учебное по-собие. -М.: ФГУП «ВНИИМС», 2006.

121. Чуновкина А.Г., Слаев В.А., Степанов А.В., Звягин Н.Д. Оценивание неопределенности измерений при использовании программ обработки данных // Измерительная техника. 2008, №7.

122. WELMEC 7.1. Software Requirements on the Basis of the Measuring Instruments Directive (MID).

123. WELMEC 7.2. Issue 1. Software Guide (Measuring Instruments Directive 2004/22/EC).

124. МИ 2174-91. ГСИ. Аттестация алгорнитмов и программ обработки данных при измерениях. Основные положения.

125. МИ 2955-2005. ГСИ. Типовая методика аттестации программного обеспечения средств измерений и порядок ее проведения.

126. Cook H.R., Сох M.G., Dainton М.Р., Harris P.M. Methodology for Testing Spreadsheets and Other Packages Used in Metrology. A Case Study. Report to National Measurements System Policy Unit, 1999.

127. Система добровольной сертификации программного обеспечения средств измерения и информационно-измерительных систем и аппаратно-программных комплексов. Правила функционирования М.: ВНИИМС, 2009

128. Левин С.Ф. Статистические методы и метрологическая аттестация программного обеспечения измерительных систем // Измерительная техника. 2008, № 11.

129. Дудыкин А.А., Кудеяров Ю.А., Паньков А.Н. Проблемы аттестации встроенного программного обеспечения средств измерений // Законодательная и прикладная метрология. — 2007, № 1.

130. Марченко Е. Что такое качество программного обеспечения? -Электронный ресурс., режим доступа: http://software-testing.ru/library/testing/ general-testing/124.

131. Попов, А. Метрики качества программного обеспечения Электронный ресурс., режим доступа: http://www.pmprofy.ru/content/ rus/67/672-article.asp

132. Терехов А. Современные модели качества программного обеспечения // BYTE/Россия №12, 1999.

133. Клюев В.В., Соснин Ф.Р., Ковалев А.В. и др. Неразрушающий контроль и диагностика / Ред. Клюев В.В.- 3-е изд., перераб. и доп. М.: Машиностроение, 2005 - 656 е.: ил.