Разработка библиотеки конечных элементов для САПР упругих конструкций герконов
21
13-08-01250); программы фундаментальных исследований Отделения нанотехнологий и информационных технологий РАН (проект № 2.11); программы ЕБТЬДТЯиЗ: проекты 1.2/ЕЬМ-121/2011/13, 2.1/ЕЬМ/184/2011/14 (2012—2013 гг.).
ЛИТЕРАТУРА
Беллман Р., Дрейфус С. Прикладные задачи динамического программирования / Пер. с англ.; Под ред.
А. А. Первозванского. М.: Наука, 1965.
Игорь Васильевич Лысенко
Анатолий Иванович Птушкин
Борис Владимирович Соколов
Сведения об авторах д-р техн. наук, профессор; Санкт-Петербургский институт информатики и автоматизации РАН, лаборатория информационно-аналитических технологий в экономике; E-mail: ilys@iias.spb.su
канд. техн. наук, профессор; Военно-космическая академия им. А. Ф. Можайского, кафедра организации эксплуатации вооружения и военной техники, Санкт-Петербург; E-mail: anatoly.ptushkin2011@yandex.ru д-р техн. наук, профессор; Санкт-Петербургский институт информатики и автоматизации РАН, зам. директора по научной работе; E-mail: sokol@iias.spb.su
Рекомендована кафедрой СПИИРАН
Поступила в редакцию 15.03.10 г.
УДК 621.01
В. Л. Ткалич, Р. Я. Лабковская, О. И. Пирожникова
РАЗРАБОТКА БИБЛИОТЕКИ КОНЕЧНЫХ ЭЛЕМЕНТОВ ДЛЯ САПР УПРУГИХ КОНСТРУКЦИЙ ГЕРКОНОВ
Обсуждается вопрос формирования библиотеки конечных элементов для САПР современных конструкций герконов в целях создания защитных сигнализационных систем. Приведены данные по истории развития разработок и производства герконов.
Ключевые слова: герконы, конечные элементы, упругие чувствительные элементы, САПР, защитные сигнализационные системы.
Несмотря на почти 40-летнюю историю применения значение магнитоуправляемых герметизированных контактов (герконов) как элементов коммутационной техники не снижается. Мировой рынок герконов достаточно стабилен, им практически нет замены во многих отраслях промышленности и народного хозяйства.
Достоинства и простота конструкции герконов, способность управлять внешним магнитным полем (созданным электрическим током или постоянным магнитом) позволяют применять их в качестве контактных пар всевозможных реле, концевых выключателей, тумблеров, кнопок клавиатуры, датчиков положения, скорости и ускорения, индикаторов перемещения и распределения сигналов, преобразователей неэлектрических величин в электрические аналоги, элементов высокочастотных и измерительных систем спецтехники, защитных сигнализационных систем и т. п.
Ведущими производителями герконов на мировом рынке являются фирмы OKI (Япония), Hamlin (США), Standex (Великобритания), Philips, C. P. Clare, Fujitsu (Япония), Gunther (Германия), Aleph Nippon (Япония), Рязанский завод металлокерамических приборов (Россия).
22 В. Л. Ткалич, Р. Я. Лабковская, О. И. Пирожникова
Статичность объемов производства герконов в мире связана с рядом их неоспоримых качеств, таких как: полностью герметизированный металлический контакт, в связи с чем герко-ны могут работать в условиях повышенной влажности и запыленности, в агрессивных средах, при температурах от -60 до +150 °С; малая мощность управления (50—200 мВт); низкое электрическое сопротивление (0,05—0,2 Ом); высокое сопротивление изоляции (1010—1012 Ом); быстродействие (0,5—1,5 мс); полная гальваническая развязка цепей управления и нагрузок; большой срок службы (106—108 переключений) [1].
Производство герконов, как правило, является высокоавтоматизированным, прецизионным, требующим высокой квалификации обслуживающего персонала и специального технологического оборудования, что связано с необходимостью разработки САПР герконовой язычковой и мембранной элементной базы.
Формирование библиотеки конечных элементов в приложении к плоским и мембранным упругим чувствительным элементам (УЧЭ) датчиков систем управления позволит применять данную элементную базу при создании высоконадежных и конкурентоспособных защитных сигнализационных систем.
В настоящее время наиболее широко применяются модифицированный элемент с пятью степенями свободы в узлах и изопараметрический элемент с 36 степенями свободы в узлах, относящиеся к четырехугольным элементам пластин и оболочкам.
Рассмотрим элемент первого типа: в системе криволинейных ортогональных координат в каждом узле вектор его обобщенных узловых перемещений имеет пять компонент — три компоненты вектора перемещений и два угла поворота, т. е.
{U}T = [W!^^^].
Углы поворота определяются по формулам
1 dw3 U2 ; . = dw3 + W1
ttl =----; =--1 ,
1 r дф R2 2 & V
где s, ф — меридиональная и окружная криволинейные координаты срединной поверхности оболочки; R1, R2 — главные радиусы кривизны; r — радиус оболочки в сечении, перпендикулярном оси вращения.
В общем случае элемент представляет собой произвольный четырехугольник на срединной поверхности оболочки. Для формирования определяющих характеристик элемента используется нормализованная система координат -, п, на которую отображается четырехугольник.
Криволинейные ортогональные координаты s, ф и произвольная осесимметричная поверхность задаются в системе координат -, п формулами:
4
ф = Z Li(- п )ф i;
i=1
4
s = Z Li(- n) si;
i=1
r
4 4 fbr ^ 4 „ J дгЛ
= Z Hi (-, n)r- +Z Hu (-, n) d- +Z H2i (-, n)
i=1 i=1
i=1
dnJ/
4 4 дг 4 дг
7 = X И, (П)7, + X Нц (П ) — + X Н2; (П ) —. ¿=1 ,=1 £ ,=1 Здесь п) — билинейные полиномы Лагранжа; п), Н1г(^, п), п) — кубические функции Эрмита; si, ф,, г,, г, — узловые значения аппроксимируемых функций.
Разработка библиотеки конечных элементов для САПР упругих конструкций герконов 23 Аппроксимирующие функции имеют следующий вид:
L = 4(1+^i)(i+nni);
Hi = Ф1 (^)ф1 (п); Нц = ф1 (п)Ф2 (0; H2i = Ф1 (^)Ф2 (п); Ф1 (К) = -4(2 + 3К0 - К30 ); Ф2 ) = 4 (К30 + К20 - К0 -1);
К0 ^0, П0; ^ ; П0 = пп; К , n,
где n — узловые значения нормализованных координат.
Для оболочки с произвольной формой меридиана главные радиусы кривизны R1, R2 определяются согласно известным соотношениям дифференциальной геометрии с учетом того,
dr dz
что r = r(5), z = z(5), — = cos 0, — = sin 0 ( 0 — угол между осью вращения z и нормалью к
ds ds
срединной поверхности оболочки).
Использование аппроксимации обеспечивает плавные изменения на границах сопрягаемых элементов радиусов r, R1, R2. Производные по нормализованным координатам в узловых точках определяются по формуле
dq Л ( dq Л ( ds
П, i = 1, •••, 4.
dXJj yds J
Для оболочек канонической формы можно определять радиусы кривизны и их производные по известным зависимостям.
Для приведения общей матрицы жесткости к одной размерности с матрицей преобразования координат необходимо вводить фиктивный угол поворота узлов элемента вокруг нормали. Для этого вводится шестое узловое перемещение — угол поворота uzi, тогда расширение блоков матрицы жесткости [Ky] осуществляется по закону:
>£М, i = j;
4 =
- 3 yEhA, i * j,
где Е — модуль упругости материала, к — толщина слоя многослойной оболочки, так что элемент остается самоуравновешенным в локальной системе координат; изменение коэффи-
—3 —1
циента у в диапазоне 3-10 <у<3 • 10 практически не оказывает влияния на конечный результат.
Недостатками описанного элемента являются необходимость введения фиктивного угла поворота и2г при сборке ансамбля конечных элементов, а также низкий порядок аппроксимации тангенциальных перемещений и, как следствие, сравнительно низкая точность определения внутренних силовых факторов [2].
В этом смысле более перспективными представляются конечные элементы изопарамет-рического типа, простейшим из которых является криволинейный четырехугольник с девятью степенями свободы в узле. В общем случае границы элемента могут не совпадать с линиями кривизны. В качестве узловых неизвестных конечного элемента принимаются девять величин — три компоненты вектора перемещения и их первые производные по координатам произвольной криволинейной системы координат:
{U }T =
дц дц ди2 ди2 дыз дыз
U1--U2--U3 ■
5а1 да2 5а1 5а 2 5а1 да 2 где иг, г = 1, 2, 3, — составляющие вектора перемещений в локальной системе координат, связанной с координатными линиями; аг, г = 1, 2, 3, — координатные линии.
При сборке ансамбля в узловых точках компоненты вектора перемещений одного из элементов выбираются в качестве основных. Компоненты узлового вектора перемещений
24 В. Л. Ткалич, Р. Я. Лабковская, О. И. Пирожникова
примыкающей оболочки выражаются через основные исходя из условий сопряжения оболочек. С повышением порядка производных в векторе обобщенных узловых перемещений выполняется последовательное повышение порядка точности аппроксимации поля перемещений. Соответственно может быть повышена точность аппроксимации геометрии оболочки.
При выборе конечных элементов стержней, составляющих единый ансамбль с базовыми конечными элементами пластин и оболочек, необходимо обеспечить идентичный выбор вектора обобщенных узловых перемещений и аппроксимирующих функций [3], с тем чтобы обеспечить выполнение общих требований к конечным элементам.
Таким образом, сформирована библиотека конечных элементов УЧЭ герконов для САПР магнитоуправляемых устройств. Показано, что при использовании моментной схемы конечных элементов повышается точность расчетов деформаций и напряжений магнито-управляемых УЧЭ герконов на 5—7 %.
СПИСОК ЛИТЕРАТУРЫ
1. Карабанов С. М., Майзельс Р. М., Шоффа В. Н. Магнитоуправляемые контакты (герконы) и изделия на их основе. Справочное руководство. М.: Интеллект, 2011. 432 с.
2. Ткалич В. Л. Надежность магнитоуправляемых контактов в системах управления. Монография. СПб: СПб ГИТМО(ТУ) 2000. С. 100.
3. Ткалич В. Л., Лабковская Р. Я., Пирожникова О. И. Метод повышения надежности упругих чувствительных элементов систем управления и автоматики // Науч.-техн. вестн. СПбГУ ИТМО. 2011. № 1(71). С. 136—138.
Вера Леонидовна Ткалич
Римма Яновна Лабковская
Ольга Игоревна Пирожникова
Рекомендована кафедрой проектирования и безопасности компьютерных систем
Сведения об авторах д-р техн. наук, профессор; Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кафедра проектирования и безопасности компьютерных систем; E-mail: Vera_Leonidovna_Tkalich@mail.ru
аспирант; Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кафедра проектирования и безопасности компьютерных систем; E-mail: studsovet_itmo@mail.ru
аспирант; Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кафедра проектирования и безопасности компьютерных систем; E-mail: studsovet_itmo@mail.ru
Поступила в редакцию 06.03.12 г.