Расчет распределения зарядов пластин при наличии внешнего несимметричного поля

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

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

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

В качестве подтверждения правильности работы алгоритма приведены решения модельных задач восстановления зарядов по распределению потенциалов.

СПИСОК ЛИТЕРАТУРЫ

1. Тарновский А.С. Об определении понятий «потенциал» и «потенциальное поле» // Электричество. - 2000. - № 1. - С. 63-64.

2. Шишигин С.Л. Построение двумерной картины электростатического поля // Электричество. - 2004. - № 3. - С. 53-58.

3. Миролюбов Н.Н., Костенко М.В., Левинштейн М.Л., Тиходе-ев Н.Н. Методы расчета электростатических полей. - М.: Высшая школа, 1963. - 415 с.

4. Исаев Ю.Н., Кулешова Е.О. Расчет распределения зарядов электрического поля на поверхности плоской системы электродов,

помещенной во внешнее электростатическое поле // Энергетика: экология, надежность, безопасность: Матер. XII Всеросс. на-учно-техн. конф. - Томск: Изд-во ТПУ, 2006. - С. 69-72.

5. Тихонов А.Н., Арсенин В.Я., Тимонов А.А. Математические методы компьютерной томографии. - М.: Наука, 1985. - 160 с.

6. Тихонов А.Н., Арсенин В.Я. Методы решения некорректно поставленных задач. - М.: Наука, 1986. - 286 с.

Поступила 05.05.2008 г.

УДК 621.372.4:537.52

РАСЧЕТ РАСПРЕДЕЛЕНИЯ ЗАРЯДОВ ПЛАСТИН ПРИ НАЛИЧИИ ВНЕШНЕГО НЕСИММЕТРИЧНОГО ПОЛЯ

Е.О. Кулешова, Ю.Н. Исаев, О.В. Васильева, Д.А. Русол*

Томский политехнический университет E-mail: kuleshova@elti.tpu.ru *ОАО «НИПИ», г. Томск

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

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

При воздействии внешнего поля возникает перераспределение зарядов, результирующее поле является суперпозицией полей создаваемого заря-

дами на поверхности электродов и внешним источником. Этот факт выражается в виде интегрального уравнения Фредгольма первого рода [1, 2].

-и0(г0) + и0 = —!— [ ст(г) Жт, г = {х,у,г},

4п—о Ъ I т - то1

то = {ХрУ»2о}> г>го 6 Ъ

Принятые обозначения соответствуют [1].

Перепишем уравнение (1) в цилиндрических координатах и учтем, что уравнение должно быть записано для поверхности электрода ¿=0 (2) из [1].

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

ы(р,в) = Х Сп,тК (р)ерт,

п,т

и (р,в) = Х А„тК (р)е1вт.

п,т

где С„т и Лпт - коэффициенты разложения рядов ы=(р,в) и и=(р,в) соответственно, Дт(р) - радиальная составляющая.

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

Ы(р) = Xсп,К(р1 ит (р) = Xапрп‘р т = Ъ2’---

п п

И, если функции разложения являются собственными, то соотношения между коэффициентами имеют простой вид: сп=а„/Япт, т. е.

ыт(р)=х^гкптр ит(р)=Xа„к(р),

п п,т п

где Япт - собственные числа.

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

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

где Ч(р) - собственные функции, Я - собственные числа.

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

ы

п п

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

ир) в виде рядов (3) подставляются в интегральное уравнение (1), которое редуцируется в алгебраическое, в силу свойств собственных функций (2):

XапЧп

п п

В силу линейной независимости собственных функций получаем сп=а„/Я„. Таким образом, получаем искомое решение в виде разложения по собственным функциям

ы(р)=ХЯ ч (р)-

п Яп

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

ы(р)=Хт^~'¥ п (р),

п Яп +15

где е - бесконечно малая величина, I - мнимая единица.

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

Задача нахождения таких собственных функций Чп(р) сводится к решению интегрального уравнения Фредгольма второго рода

Ч*

здесь К(р,р) - ядро интегрального уравнения,

Як - собственные функции и собственные значения интегрального уравнения (4), соответственно,

Решение интегрального уравнения (4) будем искать в факторизованном виде:

Ч

Прежде всего определим вид азимутальной функции 0(в). Отметим, что функция 0(в) должна быть непрерывной и периодической функцией угла в с периодом 2п. Подставляя решение (5) в уравнение (4), получим

2п Я

Л К (р,в, р',в') Я(р)0(в)рёрёв = ЯЯ(р)0{в).

0 0

Произведем следующую замену переменных: %=в'=в, йС,=йв, в =С+ в, тогда

2п Я

//к (р, р\С)Я(р'Ш +в)рёрёс =

0 0

= ЯЯ

Из уравнения (6), в силу непрерывности, периодичности и единственности решения, следует, что 0(£+в)=0(£)0(в). Общее решение этого ура-

внения, имеющее период 2п, хорошо известно и имеет вид ехр(тв). Следовательно, функция 0(в) записывается следующим образом:

0

Для определения радиальной матрицы перехода подставим выражение (7) в (6) и, после преобразований, получаем однородное интегральное уравнение Фредгольма второго рода:

Я 2п

/рёрЯ

0 0

Введем обозначение:

2п

Кт(р,р) = / ёСехр(т£)К(р,р

0

Выражение (8) сводиться к:

Я

/рё рЯт

0

Для выражения (9) было рассчитано ядро при различных азимутальных составляющих т=0,1,2,... Пространственное изображение полученных собственных функций приведено на рис. 1.

В выражении (9) введем обозначение

Рис. 2. Пространственное распределение, изолинии и профили потенциальной функции и(ху) по осям: 1) х; 2) у

Рис. 3. Пространственное распределение, изолинии и профили 50

Рис. 4. Пространственное распределение, изолинии и профили распределения зарядов Б(хгу) по осям: 1) х; 2) у

и

Рис. 5. Пространственное распределение, изолинии и профили потенциальной функции и(хгу) по осям: 1) х; 2) у

Тт (р, р') = К (р, р')р ёр. Они имеют N соответствующих линейно независи-

т мых собственных векторов, удовлетворяющих

Тт = кт р Ар, у = 0... N, Ар = N, р - ъес~ условию [6]:

ТтЧт (р) = ЯтЧт (р),

Матрицы Т не симметричны и имеют N различных собственных значений Ят, у=0,1,2,...,Ж где Ym(P)=[Чom(P),Чlm(P),Ч2m(P),...ЧNm(P)]

Для определения коэффициентов разложения необходимо использовать собственные функции ФДг) транспонированной матрицы [ Тт]Т

(Тт)ТФт (р) = ЯтФт (р),

где Ф"(р)=[ Фот(р), Ф/Чр), Ф2и(р),... Ф/(р)]Т.

Известно [6], что собственные векторы матриц Тт и [Тт]Т, соответствующие различным собственным значениям, ортогональны и для них справедливы условия:

Я

/рё рф т

т(Р)^7 (Р) = Ntf, j = 0,1,2,..., N,

где 6 у ^ У, Щт - условия нормировки.

Продемонстрируем работу алгоритма на нескольких азимутальных зависимостях и покажем совпадение результатов расчета по выше описанному алгоритму с непосредственным вычислением интеграла Фредгольма (1).

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

и(р,в) = 0,1-0,2р2 -0,5р2совв + 0,5роо$2в.

Рассматривать данное распределение потенциала будем как сумму трех составляющих при различных значениях азимутального индекса т. В нашем случае т=0,1,2.

При т=0 составляющая распределения потенциала имеет вид ио(р)=0,1-0,2р2, при т=1 и1(р,в)=-0,5р2ео8в, при т=2 и2(р,в)=-0,5рсо82в.

Модельное распределение потенциала и(р,в), являющееся суммой всех выше перечисленных составляющих, представлено на рис. 2.

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

Для т=0 и к=0 составляющая распределения потенциала ио(р)=0,1=]СЛп°Фп°(р). Умножим правую и левую часть равенства на Ч0(р):

Я Я

/ рё рЧ к (р)и 0=х а / рё рЧк (р)ф: (р).

0 п 0

Учитывая условия ортогональности, получаем:

]рё рЧ Кр)и 0=х а: ^т.

0 п

Отсюда определяем коэффициенты Лп, используя выражение

А: ^ 1Р'1Н / и (Р) Чп (Р) РёР.

Используя (3), определяем коэффициенты разложения распределения зарядов

С 0 = А

я0 '

Распределение зарядов определим как

Ы0(р) = ХФ °(р) С °.

п=0

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

Аналогично определяем коэффициенты разложения ы(р, в) при т=0 для полинома второй степени ио(р)=-0,2р2.

Результирующее распределение зарядов для нулевой азимутальной составляющей изображено на рис. 3.

Аналогично были рассмотрены составляющие распределения потенциала при т=1 и т=2.

На рис. 4 представлено пространственное распределение зарядов, соответствующее модельному распределению потенциала.

Критерием правильности расчетов является непосредственное вычисление интеграла (2). В качестве подынтегральной функции ы(р',в) подставляем полученное распределение зарядов. В результате получаем распределение потенциала, представленное на рис. 5, которое совпадает с модельным распределением потенциала, рис. 6. Данный пример наглядно демонстрирует эффективность работы алгоритма. Переход от интегральных уравнений к системе алгебраических уравнений существенно упрощает решение и требует меньших временных затрат.

Рис. 6. Профили потенциальной модельной функции U(x,y) и восстановленного распределения потенциала соответственно по осям: 1,2) х; 3,4) y

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

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

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

СПИСОК ЛИТЕРАТУРЫ

1. Исаев Ю.Н., Кулешова Е.О., Васильева О.В., Русол Д.А. Расчет распределения зарядов пластин при наличии внешнего несимметричного поля // Известия Томского политехнического университета. - 2008. - Т. 312. - № 4. - С. 70-75.

2. Миролюбов Н.Н., Костенко М.В., Левинштейн М.Л., Тиходе-ев Н.Н. Методы расчета электростатических полей. - М.: Высшая школа, 1963. - 415 с.

3. Тихонов А.Н., Арсенин В.Я., Тимонов А.А. Математические методы компьютерной томографии. - М.: Наука, 1985. - 160 с.

4. Преображенский Н.Г., Пикалов В.В. Неустойчивые задачи диагностики плазмы. - Новосибирск: Наука, 1982. - 238 с.

5. Тихонов А.Н., Арсенин В.Я. Методы решения некорректно поставленных задач. - М.: Наука, 1986. - 286 с.

6. Баринов В.А., Совалов С.А. Режимы энергосистем: методы анализа и управления. - М.: Энергоатомиздат, 1990. - 440 с.

Поступила 05.05.2008 г.

УДК 621.311.001

ЗАКОНОМЕРНОСТИ ФОРМИРОВАНИЯ СПЕКТРОВ И СТРУКТУР СОБСТВЕННЫХ ЭЛЕКТРОМЕХАНИЧЕСКИХ КОЛЕБАНИЙ

С.В. Шиловский

Институт автоматизации энергетических систем, г. Новосибирск E-mail: shi@iaes.ru

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

Введение

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

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

В основе электромеханических переходных процессов энергообъединений лежат свободные колебательные движения синхронных машин.

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

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