Научная статья на тему 'Общие подходы к моделированию динамики электрофрикционного взаимодействия'

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

CC BY
136
23
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЭЛЕКТРОФРИКЦИОННОЕ ВЗАИМОДЕЙСТВИЕ / СКОЛЬЗЯЩИЙ ЭЛЕКТРИЧЕСКИЙ КОНТАКТ / ИМИТАЦИОННОЕ МОДЕЛИРОВАНИЕ / ФРАКТАЛЬНЫЕ КЛАСТЕРЫ / ELECTRO FRICTION ENGAGEMENT / SLIDE ELECTRIC CONTACT / SIMULATION SYSTEM / FRACTAL CLUSTERS

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

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

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

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

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

GENERAL APPROACHES TO SIMULATION OF THE DYNAMICS ELECTRO FRICTION ENGAGEMENT

Electro friction engagement in electric machines is considered. General approach for building of the dynamic simulation system slide electric contact is developed. Mathematical relations for general groups components of state vector of contact element are given.

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

УДК 51-74; 62-768.3

И. В. Плохов, А. В. Ильин, О. И. Козырева

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

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

Ключевые слова: электрофрикционное взаимодействие, скользящий электрический контакт, имитационное моделирование, фрактальные кластеры.

Электрофрикционным взаимодействием (ЭФВ) называют процесс передачи электрического тока через скользящий механический контакт. Такие процессы протекают в электрических машинах, снабжённым узлами скользящего токосъёма (УСТ), которые состоят из электропроводных щёток, обычно, графитовых, или металлографитовых, прижатых с помощью щёткодержателей к поверхности контактных колец или коллекторов. При ЭФВ происходит тепловыделение в зоне скользящего контакта (СК), механические вибрации, износ контактирующих поверхностей и другие явления, которые системно взаимосвязаны. Исследование ЭФВ проводят с целью создания новых более эффективных контактных пар. Экспериментально измерить большую часть динамических характеристик ЭФВ не представляется возможным, так как механические, тепловые и электрические процессы протекают в зонах контакта, закрытых для непосредственного наблюдения.

Динамические процессы, протекающие при работе УСТ можно разделить на четыре основных категории: 1) механические; 2) электромагнитные; 3) тепловые; 4) химические. Процессы взаимосвязаны непосредственно, или через промежуточные функциональные цепи (табл.).

Будем строить систему имитационного моделирования УСТ исходя из положений [1]. Разделим задачу моделирования на три части:

1. Получение механических колебаний щётки в результате кинематических возмущающих воздействий и определение неустойчивости СК;

2. Расчёт динамики токораспределения между параллельно работающими щётками;

3. Расчёт динамических процессов токопередачи в переходном слое.

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

Таблица

Процессы в УСТ Краткое описание Взаимовлияние

1. Механические Колебания щеток, вызванные кинематическими и динамическими возмущающими воздействиями. і > і

Процессы фрикционного взаимодействия микрорельефов. Ч І > і і і 1 к

2. Электромагнитные Колебательные процессы динамического токораспределения между параллельными СК. к ‘ < Г і

Эволюция контактной проводимости и токопередача в переходном слое СК. і і 1 Г і і ) к

3. Тепловые Эволюция нестационарного температурного поля в переходном слое СК. і г 1 к < Г і 1 г

Интегральный нагрев контактных пар. 1 > Г

4. Химические Окислительновосстановительные и электролитические реакции. г 1 г 1

Рис. 1. Система имитационного моделирования

Моделирование динамики электрофрикционного взаимодействия ЭФВ представляет собой многофакторную задачу исследования процессов формирования проводимости в переходном слое СК. Развитие данных процессов связано с внутренними свойствами контактного перехода, которые обусловлены характеристиками взаимодействующих тел и внешними воздействиями. Формирование внешних воздействий, таких как механические колебания щётки и ко-

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

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

Представим контакт двух микрорельефов как систему дискретных элементов. При этом будем рассматривать переходный слой как некоторое «третье тело», состоящее из множества взаимосвязанных площадок, размещённых на плоскости, равноудалённой от обоих контактирующих микрорельефов. Таким образом «третье тело» — переходный слой контакта заменяем более абстрактным объектом, позволяющим перейти к чисто математическим определениям и операциям. Назовём данный объект переходной плоскостью электрического контакта (ПП).

Разобьём переходную плоскость на множество элементарных участков одинаковой площади ДS, которые назовём контактными элементами (КЭ). Они располагаются в узлах регулярной сетки с шагом Дx по оси X и Ду — по оси Y (в практике вычислительного моделирования принимают Дx = Ду).

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

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

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

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

Таким образом, обобщённая модель скользящего электрического контакта предстаёт как регулярная совокупность КЭ на переходной плоскости, обладающих свойствами и объединённых отношениями-связями. Однотипные проводящие КЭ образуют соответствующие кластеры. Свойства КЭ проявляются при приложении внешних воздействий к контактной паре. Бинарные отношения КЭ устанавливаются через посредство физических полей.

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

Основные внешние воздействия, поступающие на входы системы «скользящий электрический контакт», можно подразделить на следующие категории: 1) контролируемые и целенаправленно регулируемые воздействия; к ним относятся постоянные составляющие электрического напряжения на контактной паре и усилия нажатия, скорость относительного перемещения элементов контактной пары; 2) контролируемые, но нерегулируемые воздействия (внешние условия); к ним относятся переменные составляющие электрического напряжения на контактной паре и усилия нажатия, температура и химический состав среды, рельеф контактирующих поверхностей; 3) неконтролируемые воздействия (возмущения).

Переходный слой представим совокупностью КЭ, размещённых в узлах

регулярной сетки на 1111 [2]. Каждый КЭ опишем вектором состояния V, имеющим п компонент. Таким образом, определим контакт как ограниченную область W ^мерного векторного пространства состояний, которую можно описать соответствующей матрицей п-го порядка. Определяющей операцией на W является преобразование координат данного пространства, которое задаётся гомоморфизмом класса Е (силовой вектор) на W с определяющей операцией f и автоморфизмом области W с определяющими операциями и w2. Операция f представляет собой абстрактное умножение силового вектора Е на W. Базисом векторного пространства является орт единичных параметров КЭ. Т.е. класс W с определяющими операциями и w2 представляет собой кольцо конечной размерности с единицей. Заданные внутри W бинарные операции можно охарактеризовать в общем виде следующим образом: — абстрактное сложение, реализующее полевую суперпозицию в пространстве состояний; — абстрактное

умножение, реализующее эволюционные отношения элементов класса W.

Компонентами силового вектора Е являются независимые входные воздействия: усилие нажатия ЕН, скорость относительного перемещения в кон-

тактной паре vk, электрическое напряжение и, приложенное к контакту извне, температура внешней среды 0СР, вектор концентрации химических реагентов среды Z.

Основные компоненты вектора состояния контактного элемента разобьём на следующие классы:

1. Механический VM:

- контактная жёсткость СМ;

- контактное демпфирование Кк

- контактное сближение Yk;

- контактное усилие Fk;

- микротвёрдость q;

- плотность р;

- коэффициент граничного трения ктр.

2. Электрический ^:

- напряжение U (одинаково для всех контактных элементов);

- напряжение фриттинга ^;

- электропроводимость рэл или электросопротивление R;

- контактная ёмкость СЭЛ;

- ток I через КЭ.

3. Тепловой VI:

- теплоёмкость С0;

- теплопроводность к;

- температура плавления 0П;

- температура испарения 0И;

- мгновенная температура 0 .

4. Химический ^:

- вектор химического состава % ;

- интенсивность основных химических реакций Vh.

Определим вектор состояния суммой четырёх составляющих V=V + VE + V, + V. Описанные классы компонентов вектора состояния содержат не только параметрические, но и энергетические составляющие Fk , I , 0. Поэтому классифицируем компоненты по их роли в определяющих операциях №1 и -№2 на '. В результате получим два подкласса — параметрический и энергетический, которые определяют составляющие вектора состояния V = Ув + Ур.

Воздействие силового вектора Р на вектор состояния V передаётся через энергетическую составляющую VF следующим образом. Для подкласса Vм усилие нажатия Рн влияет на величину контактного усилия РК и сближения микрорельефов Yk в каждом КЭ. Т.к. все компоненты подкласса Vв в разной степени зависят от контактного усилия, то VB=f(FH).

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

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

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

По отношению к основной функции — токопередаче — векторное поле состояний W представим как множество ограниченных зон, которые можно разделим на два основных типа: проводящих элементов Wп и непроводящих элементов WН. Термин «непроводящая зона» является в определённой мере условным, т.к. проводимость существует, но она существенно ниже, чем для «проводящих зон».

Зоны контактной проводимости представляют собой синергетические структуры — кластеры, развитие которых напоминает перколяцию с вытеснением [3-6]. Растущие кластеры проводимости объединяются во фрактальные образования — кластерные агрегаты, которые также являются кластерами. Их размерность в заданный момент времени определяет мгновенное значение контактной проводимости [5, 6].

Абстрактная операция реали-

зует эволюционные отношения (ЭО) элементов класса W, которые в процессе развития кластеров сохраняют постоянную структуру.

Прямые ЭО устанавливают законы развития системы внутри каждого цикла генерации кластеров, наследственные — определяют правила передачи эволюционных свойств из поколения в поколение.

Преобразование пространства состояний связано с наличием автоморфизма в области W и инерционности компонентов вектора состояния. Структуру характеризует количественное и качественное соотношение классов Wп и WН, представляющих кластерные агрегаты. Визуально представим данные структуры, используя сечения пространства состояний по определённым характеристическим координатам. Условимся базовыми сечения по электропроводимости рЭЛ, температуре 0 и контактному сближению Ук, которые при вычислениях представим матрицами натуральных чисел, а при графическом отображении — оттеночными портретами. На них интенсивности цвета каждого эле-

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

Состояние КЭ определим положением изображающей точки в п-мерном пространстве ', а изменение состояния — фазовой траекторией. Состояние системы представим изображающей областью, охватывающей все изображающие точки КЭ, а изменение состояния системы — семейством фазовых траекторий КЭ, представляющих след изображающей области в п-мерном пространстве состояний. Векторное подпространство ' опишем матрицей ||Уу|| векторов состояния.

% ^2 . К

VI = сч > 2 2 >

Vm2 . V тп

(1)

Гомоморфизм силового вектора Е с компонентами Рн, ук, и, Z в пространство состояний ' зададим операцией абстрактного умножения матрицы элементов ||у„| на элемент Е.

Р ®Ж=Р ®| к| =

V11 ^12 . . ^п Р ® v11 Р ®^2 . .. Р ® Vln

р ® "^21 ^22 . = Р ® V21 Р ® V22 .

Vт2 . . V тп Р ® Vтl Р ® Vт2 . .. Р ® Vтn

Автоморфизм на множестве ' определим операцией абстрактного суммирования его элементов по следующему правилу:

=

\ 1 ® \ 1+1, если 3 \ 1+1 еЖ V, 1 ® V, 1 _1, если 3 V, 1 _1 е Ж

\ 1 ® V+и, если3V+и V, 1 ® V,_11, если3V,_11 еЖ

(3)

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

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

- механический класс:

Ск = I (,7,, Н); Кк = / (9,7,, Н); 7, = I (к, С^, Кк, V,);

Рк = I(,7,, I); р=119, Н);

- электрический класс:

Рэл = 1 (,7к,Н, Х ),

где X - фрактальная размерность кластера, содержащего КЭ; иф = 1(9, Н,7к);

Сэл = I (к, Н, 9);

1 = 1 (и, рэл );

- тепловой класс:

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

С9 = I (р,9, Н); к = I(, Н,9);

9 = УI1, рЭЛ , С9, Ч, 7К , ^ );

9п = I (Н);

9я = I (Н);

- химический класс:

Н = I ( ,9);

^ =I (н,9,Р).

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

Выводы

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

Литература

1. Плохов И. В. Модель динамики токопередачи через скользящий контакт Электротехника. №2. 2005. С. 28-33.

2. Плохов И. В. Комплексная диагностика и прогнозирование технического состояния узлов скользящего токосъёма турбогенераторов. Диссертация д.т.н. С.-Петербург, СПбГПУ, 2002.

3. Федер Е. Фракталы. М. :Мир. 1991.

4. Смирнов Б.М. Физика фрактальных кластеров. М. : Наука. 1991.

5. Фракталы в физике / Труды VI международного симпозиума по фракталам в физике. М. : Мир. 1988.

6. Плохов И. В., Савраев И. Е. Вычислительная модель развития перколяцион-ных кластеров контактной проводимости // Труды Псковского политехнического института. С.-Петербург/Псков : СПбГТУ, 1997. Вып.1. С.51-54.

I. V. Plohov, A. V. Iliin, O. I. Kozyreva

GENERAL APPROACHES TO SIMULATION OF THE DYNAMICS ELECTRO FRICTION ENGAGEMENT

Electro friction engagement in electric machines is considered. General approach for building of the dynamic simulation system slide electric contact is developed. Mathematical relations for general groups components of state vector of contact element are given.

Keywords: electro friction engagement, slide electric contact, simulation system, fractal clusters.

Плохов Игорь Владимирович — заведующий кафедрой «Электропривод и системы автоматизации» ФГБОУ ВПО ПсковГУ, д.т.н., профессор, igor_plohov@list.ru.

Ильин Александр Викторович — старший преподаватель кафедры «Электропривод и системы автоматизации» ФГБОУ ВПО ПсковГУ, аспирант, al.ilyin@yandex.ru.

Козырева Оксана Игоревна — старший лаборант кафедры «Электропривод и системы автоматизации» ФГБОУ ВПО ПсковГУ, аспирант, ks_33n@rambler.ru.

УДК 51-74; 62-768.3

И. В. Плохов, А. В. Ильин, О. И. Козырева

ОПРЕДЕЛЕНИЕ КОМПОНЕНТ ВЕКТОРА СОСТОЯНИЯ ДИНАМИЧЕСКОЙ МОДЕЛИ ЭЛЕКТРОФРИКЦИОННОГО ВЗАИМОДЕЙСТВИЯ

Рассматривается динамическая модель электрофрикционного взаимодействия. Заданы три основных класса компонент вектора состояния: механический, тепловой и химический, для которых приведены основные математические соотношения.

Ключевые слова: электрофрикционное взаимодействие, компоненты вектора состояния, скользящий электрический контакт, контактный элемент, динамическая модель.

Для реализации динамической модели электрофрикционного взаимодействия (ЭФВ) необходимо определить компоненты вектора состояния [1, 2]. В данной статье уделим внимание заданию трёх основных классов компонент:

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