На них цифрой 1 обозначены траектории нагружения для оболочек со свободно смещающимися в продольном направлении краями, цифрой 2 отмечен случай всестороннего сжатия оболочки равномерным давлением при свободно смещающихся в продольном направлении краях, а цифрой 3 - траектории нагружения для оболочек, у которых на краях запрещены продольные перемещения. Для некоторых точек приведённых траекторий нагру-жения на рис. 2 и 3 показаны продольные сечения оболочки с профилем её деформированной срединной поверхности.
Заключение.
Полученное аналитическое решение задачи об осесимметричном деформировании цилиндрической оболочки имеет непреходящую ценность, как и любое другое точное решение в линейной и нелинейной теории стержней, пластин и оболочек, к тому же оно описывает основную ветвь траектории нагружения оболочки перед неосесимметричной потерей устойчивости. Точность определения координат точек бифуркации и построение закритических ветвей деформирования оболочки зависит от точности построения основной, докритической траектории нагружения. Поэтому представленное решение может рассматриваться как эталонное для численных подходов, применяемых к решению задач о неосесимметричной потере устойчивости и закритическом поведении цилиндрических оболочек, на этапе расчёта их предварительного деформированного состояния.
Литература
1. Вольмир А.С. Устойчивость деформируемых систем. М.: Наука. 1967, 984 с.
2. Григолюк Э.И., Кабанов В.В. Устойчивость оболочек. М.: Наука, 1978. 359 с.
3. Тимошенко С.П., Войновский-Кригер С. Пластинки и оболочки. М.: Физматгиз, 1963.
636с.
Программное и математическое обеспечение автоматизированных прецизионных лазерных измерительных систем с двумя оптически
обращёнными каналами
Сиротский А. А. МГТУ «МАМИ»
Задачи высокоточного и автоматизированного контроля актуальны для всех отраслей современного машиностроения. Значительная часть задач контроля в машиностроении может быть эффективно решена только применением компьютеризированных лазерных измерительных систем с двумя оптически обращёнными каналами (ЛИС ДОК). Рассмотрены программные средства наблюдения, сбора, записи и хранения измерительной информации, а также средства и методы автоматизированной обработки данных и визуализации объектов контроля. Описаны разработанные программно-алгоритмические решения для автоматизации многоканальных прецизионных измерительных систем, которые могут быть интегрированы в действующие АСУ ТП.
Лазерные измерительные системы (ЛИС) весьма перспективны, они свободны от большинства недостатков, присущих иным методам контроля. Прежде всего, они не подвержены влиянию электромагнитных помех и наводок, легко встраиваются в действующие технологические процессы, могут быть успешно автоматизированы и интегрированы в автоматизированные системы управления технологическими процессами и производствами (АСУ ТП). При этом важнейшей составляющей участка контроля являются средства компьютеризации и автоматизации, обеспечивающие гибкость контроля и взаимодействие со смежными системами управления.
Отличительной особенностью лазерных измерительных систем с двумя оптически об-ращёнными каналами (ЛИС ДОК) является наличие оптического блока формирования двух оптически обращённых каналов (ДОК). Оптическая часть системы ДОК обеспечивает разделение исходного лазерного луча на два луча, оптически обращённых друг относительно дру-
га путём разворота одного из лучей на 180 градусов вокруг своей оси, так что нестабильности положений энергетических центров лучей находятся в противофазе и взаимно компенсируются, а измерительные перемещения при этом находятся в одной фазе и складываются.
Рассмотрим пример функциональной схемы ЛИС ДОК для контроля прямолинейности оси отверстия в деталях типа «труба». ЛИС ДОК (рис. 1) состоит из жесткого основания 1, на котором на опорах 2 устанавливается деталь 3, привод 4 перемещения первичного измерительного преобразователя (ПИП) 5 с отражателем 6. При движении ПИП 5 вдоль трассы контроля, он в каждом текущем сечении центрируется, занимая при этом положение по центру отверстия в текущем сечении. Источником излучения служит лазер 8, а формирование ДОК происходит в блоке 7. Позиционно-чувствительные фотоприёмники (ПЧФ) 9 и 10 регистрируют перемещения, пропорциональные смещениям лазерного зонда. Сигналы с ПЧФ усиливаются по координатам X и Y усилителями 11, 12, 13, 14 и поступают в ЭВМ 15, оснащённую системой сбора, регистрации, обработки информации и управления на основе аналого-во-цифрового сигнального процессора ADSP.
i 2 3 5 6 7 8
Рис. 1. Функциональная схема ЛИС ДОК для контроля положения осей отверстий.
Измерения проводятся в M сечениях с некоторым заданным шагом перемещения привода 4. В каждом сечении Mi фиксируются напряжения Uxli, Uyli, и Ux2i, Uy2i, снимаемые с выходов ПЧФ. Величины измеренных напряжений записываются в памяти ЭВМ в цифровой форме в виде таблицы. Таким образом, необходима реализация четырёхканальной системы сбора и регистрации информации.
Основой созданного компьютеризированного измерительного комплекса является специализированный модуль PCI-1202H фирмы ICP DAS, представляющий собой так называемую «DAQ - board» (Data Acquisition Board - Плата сбора данных) - измерительный модуль, встраиваемый непосредственно в персональный компьютер (ПК).
Для обеспечения программного функционирования данного модуля существуют специализированные программные среды разработки, имеющие функции программной среды, обработки и визуального представления измерительной информации (например - среда LabVIEW фирмы National Instruments).
Средства измерений, построенные на базе ПК, встраиваемых в ПК многофункциональных и многоканальных DAQ - board, внешних программно-управляемых модулей предварительной обработки сигналов и специализированных измерительных интегрированных про-
граммных оболочек для сбора, обработки и визуального представления измерительной информации, называются виртуальными измерительными системами (ВИС).
Управление такими ВИС, как правило, осуществляется через графический пользовательский интерфейс (Graphics User Interface - GUI) при помощи технологии Drag-and-Drop («Перенёс и положил») с использованием манипулирования мышью через виртуальные элементы управления, расположенные на виртуальных приборных панелях. Разработка программного обеспечения ВИС как правило проводится на языках графического программирования, одним из которых является язык Джей.
Модуль PCI-1202H имеет 12-битный АЦП (44К измерений в секунду), программируемый усилитель с максимальным усилением 1000, аналоговые входы с диапазоном от -10 до 10 В в дифференциальном режиме, входной импеданс 10 Мом при 6 пФ. Подробнее параметры данного модуля и аналогичных ему можно найти в каталогах производителей.
Структурная схема модуля PCI-1202H показана на рис. 2. Модуль совмещается с ЭВМ семейства x86 по шине PCI, имеет встроенный формирователь шины, внутреннюю память, локальный контроллер, специализированную микросхему управления аналоговой входной частью «Magic scan controller», управляющую входным мультиплексором, программируемым усилителем, 12-битным АЦП и буфером FIFO.
XSi System
ZS
Ш BUS
ТУ
PCI Interface System
J\
V>
EPiO*
Local Syslefii Conli oHei
Siirtns Cwrtic!
ГС
10 Dits PI 1ô Lits DO
2 Channels 12Wt DA Converter
< I I —— Lttput
J 'Li 'l
Ш',
DiçLJ Üubut DA С QITÏ0
DA С GUT I
M aQic Scan Contiolioi
ЛшЬр C'Jif il
lS.IÈDil A/U Senveiisr
гт 7
-ш-
+ Mil
X
/jLLlop jipi:
Рис. 2. Структурная схема модуля РС1-1202Н.
Также в данном модуле реализованы по 16 дискретных входов и выходов и 12-битный ЦАП. Эти аппаратные возможности модуля могут быть использованы для целей реализации
управления приводом измерительного преобразователя перемещений ЛИС. Алгоритм разработанной программы представлен на рис. 3.
Рис. 3. Алгоритм программы.
Алгоритм обеспечивает непрерывное проведение измерений по четырём независимым каналам сбора данных, по два канала на каждый ПЧФ, запись результатов измерений во внутреннюю память ПК и в файлы в двух форматах: текстовом с расширением «.txt» и бинарном с расширением «.bin». Наиболее удобным является текстовый формат сохранения данных, так как он предполагает возможность простого экспорта данных для дальнейшей обработки в такие программные среды, как MS Excel и MathSoft MathCAD.
Общий вид окна программы показан на рис. 4. Окно программы состоит из трёх частей. В первой части находятся основные рычаги управления программой, здесь указывается код режима (верхний индикатор), т.е. интервал измерения с коэффициентом усиления. Чуть ниже находится элемент, который устанавливает частоту дискретизации, её можно изменить как передвижением ползунка, так и указав числовое значение в поле рядом. Ещё ниже указывается время измерения, оно вводится в миллисекундах. Также предусмотрены функциональные кнопки «Начать измерения», «STOP» и «Записать данные в файл». Если остановить измерения прежде, чем истечет указанное время, то в файле сохранится реальное время измерения, а не время, которые указывалось в настройках.
Вторая часть программы - графическая, здесь находятся четыре графика-индикатора, в которых показываются измерения в режиме реального времени. Шкалу графика можно задавать вручную или указать программе автомасштабирование._
tin Hnnii Hi; rt.ln
Рис. 4. Общий вид окна программы.
Третья часть программы состоит из графической шкалы, которая указывает процесс сохранения данных в файл. Также тут имеются четыре элемента, которые позволяют просматривать измеренные данные в массиве по четырем каналам.
Порядок работы с программой следующий:
1. Указать диапазон измеряемых данных.
2. Ввести время измерения.
3. Задать частоту дискретизации.
4. Нажать кнопку «начать измерения».
5. При необходимости остановить измерения кнопкой «STOP» или дождаться истечения указанного времени измерения.
6. После завершения операции при необходимости сохранить данные, нажав кнопку «Сохранить данные в файл».
Одно из преимуществ разработанной программы в том, что можно изменять параметры измерения (настройки программы) прямо в процессе сканирования.
Для удобства просмотра сохраненных файлов в их названии указываются текущее время, текущая дата, время измерения, частота дискретизации и коды каналов.
Разработанная программа обладает высокой гибкостью и может быть встроена в любые
компьютеризированные системы и автоматизированные комплексы.
Программа разработана специально для компьютерной системы регистрации измерительной информации с ЛИС ДОК, имеющей два двухкоординатных ПЧФ. Таким образом обеспечивается непрерывное измерение координат X и Y одновременно с двух датчиков.
Структура файла данных - построчная, с записью измерений в виде десятичных значений, разделённых пробелами или табуляторами. Пример файла данных на 10 измерений по четырём каналам измерений показан на рис. 5. Общее количество измерений программно не ограничено и определяется доступным объёмом памяти ЭВМ.
Для целей обработки и визуализации результатов измерений разработано математическое обеспечение в среде МаШСАО.
Поскольку общее количество измерений заранее не известно, то это учтено при разработке программного модуля в среде МаШСАО, который корректно работает с любым количеством входящих данных, не превышающим 1000000, что с гарантированным запасом перекрывает все возможные области применения данной системы. В среде МаШСАБ импорт данных из файла производится функцией КЕАВРКК:
А := ШАБРт(" Ба1аР11е4.1x1") (1)
где: DataFile4.txt _ имя файла данных.
Щ ОаЬаРИе4.1:н1 - Блокнот Л 1
Файл Правка Формат Справка
-0.000300 0 . 000120 0 000360 0. 000270 И [
-0.000180 0 .000180 0 ОООЗЗО 0. 000270
0.000060 0 .000240 0 ОООЗЗО 0. 000240
0.000180 0 .000270 0 000300 0. 000240
0.000420 0 .000270 0 000300 0. 000270
0.000660 0 .000300 0 000270 0. 000240
0. 000840 0 .000330 0 000270 0. 000210
0. 0009160 0 .000330 0 000240 0. 000180
0.001080 0 .000390 0 000240 0. 000150
0.001200 0 .000360 0 000270 0. 000150
1 1
<1 1
(2)
где:
Рис. 5. Структура файла данных.
Далее массив данных А разделяется на четыре массива по каналам:
их1 := С0) иу1 := С«
<
их2 := С2) иу2 := С3
их1, иу1, их2, иу2 - массивы значений напряжений по координатам с ПЧФ;
С <°> С 0> С<2) С <3> _ „ , , , - операции извлечения столбцов массива значений.
После чего проводится пересчёт измеренных напряжений в величины отклонений:
Бх1 := ■
ш
Цу1 :=
Бх2 :=
Бу 2 :=
КЦА • Ких1
иу1
КР1 • Киу1
их2 Кр2 • Ких2
иу2 К]р 2 • Киу 2
(3)
где: Kfp1, Kp2 - масштабные коэффициенты ПЧФ;
КшХ, Kux2, Kuy1, Kuy2 - коэффициенты усиления усилителей соответствующих каналов.
Математическое обеспечение без изменений корректно обрабатывает массив данных измерений любой длины, для чего производится вычисление действительного количества прочитанных строк данных REALmaxindex функцией length:
REALmaximdex := lengthA°°) (4)
Количество измерений на длине трассы контроля:
B := REAL maximdex -1 (5)
Для целей визуализации контролируемого отверстия проводится угловая аппроксимация окружности текущего сечения:
„ 2• пn
Qn:=-S--(6)
?
где: n - индекс вектора аппроксимации окружности;
S - количество точек аппроксимации; причём n := 0 • • •S. При использовании двух ПЧФ производится суммирование (усреднение) их показаний, что даёт координаты положения оси отверстия по сечениям:
Dx1 + Dx2
Dx := ■
2
(7)
X Y
m,n m,n
2
Z_
Далее, задавая матрицы " , ^, ~,п-п, в панели построения трехмерных моделей и при условии постоянства диаметра осуществляется визуализация пространственного положения и формы внутренней поверхности контролируемого отверстия, показанная на рис. 6.
'Л о
Рис. 6. Визуальная модель пространственного расположения отверстия.
Ещё больший интерес представляет воспроизведение пространственного положения самой оси отверстия, которую невозможно увидеть иными способами, отличными от математического моделирования. С этой целью создаются вектора значений отклонений положений центров отверстия и соответствующий им вектор значений трассы измерений:
XCm,п := DXm
YCm, П:= DУm (8)
ZCm,П := т
Визуализация пространственного положения оси отверстия показана на рис. 7.
На основе визуализации становится возможным оценить характер деформаций и отклонений формы деталей. Объём собранной информации позволяет получить статистические данные о средних значениях, разбросе значений, отклонениях, и провести комплексную оценку качества формообразования.
Предложенное программное и математическое обеспечение полностью отлажено и может быть применено в системах автоматизированного контроля.
Отличительной особенностью созданной системы обработки данных является её самодостаточность: в процессе использования не требуется никакой коррекции алгоритма обработки - он универсален для всех случаев. Единственным ограничением являются ресурсы ЭВМ и объём памяти.
Рис. 7. Визуальная ЭБ-модель отверстия с воспроизведением пространственного
положения оси.
Апробация разработанной методики обработки данных и визуализации показала надёжную и адекватную обработку данных из массива 30000 элементов.
Разработанная информационно-измерительная система отвечает самым высоким современным требованиям, обеспечивает высокую гибкость, и может быть гармонично интегрирована в современные АСУ ТП.
Литература
1. Сиротский А.А. Автоматизированная программная обработка результатов измерений и трехмерная визуализация объектов контроля в среде MathCAD. Материалы Международной научно-технической конференции «Системные проблемы качества, математического моделирования, информационных и электронных технологий». ISBN 5-94099-32-0. Часть 3. Секция №5 «Системные концепции в экономике, науке, технике и экологии». Москва-Сочи, 2005. М, МАМИ, 2005, 114 с., стр. 94-100.
2. Сиротский А.А. Программирование многоканального компьютеризированного комплекса сбора данных для лазерной измерительной системы. УДК 004.94. Известия Тульского государственного университета. Серия Технологическая системотехника. Выпуск 9. Труды участников Пятой международной электронной научно-технической конференции «Технологическая системотехника - 2006». Издательство ТулГУ, Тула, 2006. Стр. 139-147.
3. Сиротский А.А. Исследования точности и погрешностей лазерной измерительной системы с двумя оптически обращёнными каналами. Сборник избранных докладов 49-ой
Раздел 3. Естественные науки. МНТК ААИ России «Приоритеты развития отечественного автотракторостроения и подготовки инженерных и научных кадров» 4-го международного научного симпозиума «Современное автотракторостроение и высшая школа России», посвящённого 140-летию МГТУ МАМИ под ред. проф., д.т.н. Бахмутова С В. ISBN 5-94099-036-3. Москва, МГТУ «МАМИ», 2005. Книга 2. Стр. 459-464.
4. Сиротский А.А. Применение LABVIEW для автоматизации лазерных измерительных систем и научного эксперимента. Тезисы докладов научно-практической конференции «Образовательные, научные и инженерные приложения в среде LabVIEW и технологии National Instruments» ESEA-NI-04. Москва, РУДН, 2004.
5. Сиротский А.А. Экспериментальные исследования лазерной измерительной системы с двумя оптически обращёнными каналами. Материалы 49-й МНТК ААИ «Приоритеты развития отечественного автотракторостроения и подготовки инженерных и научных кадров». Секция 3 «Автоматизированное управление, электроника и электрооборудование автотракторной техники», часть 1. М., МГТУ «МАМИ», 2005 г. 2 стр.
6. Сиротский А.А. Технологические области применения лазерных измерительных систем для контроля прямолинейности и положения осей отверстий в различных деталях машиностроения. Материалы 49-й МНТК ААИ «Приоритеты развития отечественного автотракторостроения и подготовки инженерных и научных кадров». Секция 8 «Механообработка и сборка. Методы и средства автоматизации автотракторного производства, оптимизация технологических систем», часть 2. М., МГТУ «МАМИ», 2005 г. 3 стр.
Решение задач с плоскими треугольными элементами в граневых моделях
поверхностей
к.т.н. Суслин В.П. МГТУ «МАМИ»
Для граневых моделей поверхностей часто используются плоские треугольники в качестве аппроксимирующих элементов. К примеру, геометрические модели в формате STL используют именно такие аппроксимации. При этом с плоскими треугольными элементами решаются различные геометрические задачи, например такие, как определение принадлежности точки треугольнику, определение точки пересечения прямой и плоскости треугольника и др.
Рассмотрим первую задачу. Она является частным случаем задачи о принадлежности точки многоугольнику. Это известная задача в экранной графике. Для ее решения используется два метода, описанные в различных монографиях по машинной графике [1, 2]. Один основан на вычислении суммы углов, а второй - на подсчете числа пересечений луча, проходящего через тестируемую точку, со сторонами многоугольника. В обоих случаях требуется преобразовать систему координат так, чтобы координатная плоскость совпала с плоскостью треугольника, т.е. привести пространственную задачу к плоской.
В программном обеспечении, разработанном в лаборатории САПР МГТУ «МАМИ», используется другой метод решения задачи о принадлежности точки треугольнику, а также ряда других задач, связанных с треугольными элементами. Он не требует преобразования пространственной системы координат в плоскую, экономичен по объему вычислений и прост по логике.
Пусть имеется невырожденный треугольный элемент V1V2V3, определенный вершинами V1 =(x1 y1z1) , V2 =(x2y2z2 ) и V3 =(x3y3z3) . Введем двумерную параметризацию с
параметрами u и v так, чтобы в вершине V1 было u, v = 0, на стороне V1V2 - V = 0, 0 < u < 1 и на стороне V1V3 - u = 0, 0 < V < 1 (рис. 1). Такая параметризация связана с косоугольной системой координат [3].