Методика расчёта собственной ёмкости плоских элементов конструкции космического аппарата Текст научной статьи по специальности «Физика»

УДК 629.78

МЕТОДИКА РАСЧЁТА СОБСТВЕННОЙ ЁМКОСТИ ПЛОСКИХ ЭЛЕМЕНТОВ КОНСТРУКЦИИ КОСМИЧЕСКОГО АППАРАТА

© 2013 Н. Д. Сёмкин, В. В. Брагин, Р. А. Помельников, В. Л. Балакин

Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет)

Предложена методика расчёта распределения потенциалов на основе использования триангуляционного метода аппроксимации поверхностей, метода площадок и аналитического нахождения элементов матрицы кулоновского взаимодействия. Полученные результаты имеют хорошую сходимость с экспериментами по исследованию процессов электризации на малом космическом аппарате "АИСТ".

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

Введение

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

В настоящее время существуют несколько методик математического моделирования электризации КА [4]. Эти методики основаны на различных методах аппроксимации поверхности и таких математических методах, как метод интегральных уравнений (МИУ), называемый также методом граничных элементов, или метод конечных элементов (МКЭ). Эти математические методы имеют свои преимущества и недостатки, они применяются в зависимости от высоты орбит КА и, следовательно, в зависимости от свойств окружающей КА плазмы. Хотя для унификации математического обеспечения возможно применение МИУ для проведения анализа электризации как высокоорбитальных, так и низкоорбитальных КА.

Данный метод пространственного распределения потенциала описывается следующим соотношением [4]:

где и (г) - пространственное распределение потенциала; С (г') - относительная плотность поверхностного заряда; р(г') - распределение пространственного заряда в окружающем КА пространстве; г' -радиус-вектор центра тяжести /-го элемента поверхности КА.

Для численного решения данного интегрального уравнения относительно с(г) необходимо разбить поверхность КА на элементы, в качестве которых удобно использовать треугольники. Плотность поверхностного заряда С на каждом элементе считается постоянной, и, следовательно, интегральное уравнение сводится к системе линейных уравнений относительно плотности зарядов элементов С/:

^ АтпС/ = и/ , (2)

]

где Атп - матрица кулоновского взаи-

*

модействия; и/ - эффективный потенциал поверхности /-го элемента.

Искомым решением этой системы является матрица взаимных ёмкостей

элементов поверхности КА Стп, которая является обратной к матрице Amn, определяемой уравнением [4]:

Л

тп

=1

з/- - го

(3)

где г0 - радиус-вектор центр тяжести /-го треугольника; интегрирование по переменной г ведётся по площади треугольника. Диагональные элементы этой мат-

Г

г

}

в

ществующем аналитическом методе расчёта собственной ёмкости [5]. Данный метод нахождения собственной ёмкости сводится к аппроксимации плоской поверхности квадратной сеткой со стороной а и решению интеграла вида [5]:

¿13 а а | ^ = \ d^

dу

X - Х1)2 + (у - У1)2

(4)

рицы имеют особенность при

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

Методика расчёта собственной ёмкости пластинки

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

где х1 и у1 - координаты центра масс элемента квадратной сетки.

Для применения данного метода нахождения собственной ёмкости необходимо аппроксимировать поверхность КА треугольной сеткой.

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

Разбиваем плоскость на прямоугольники с координатами (хг, у-). Каждый прямоугольник разбиваем на две части отрезком (рис. 1) с уравнением, преобразованным к виду:

х

у = —-— (У/+1 - У/)+У/--/-— (У/+1 - У/) =Ах + в,

х/+1 х/ х/+1 х/

(5)

где

А = ^^^, В = у---х-- (У-+1 - У-).

х/+1 х/ х/+1 х/

У I У2

Ут4 УтЗ

У1

Ут2 Ут1

XI х2 х3 х4

О Хт1 Хт2 Хт3 Хт4 Хт4 Хт5 Хт5 Хт6 X

Рис. 1. Аппроксимация прямоугольной пластинки

Находим пределы интегрирования при условии:

Тогда интеграл примет вид:

х,+1 Ах+Б йу

- центр масс треугольника (хт/, ут/) I йх | лежит выше гипотенузы треугольника. х, у, л/ (х - хТ1) + (у - уТ1) Тогда интеграл (3) примет вид:

(7)

х+1 у+1

Iйх I

йу

Ч Ах+Ел\(х - хТ1)2 +(у - уп )2

(6)

- центр масс треугольника (х,, у,) лежит ниже гипотенузы треугольника.

Произведём следующие замены в уравнении (7):

С=х-хт/, Э=у-уТ/, йЭ=йу

и получим с учётом изменения области интегрирования уравнения следующее уравнение:

х/+1 Ах+В - ут/

I йх I , . . (8)

х, у, - УТ,

С2 + э2

В результате интегрирования (8) получаем разницу интегралов

I

1п^ Ах + В - уТ1 + ^(х - хТ1 )2 + (Ах + В - уТ1 )2

- 1п(у, - ут, + ^(х - хт,)2 +(у, - ут,)2 )

йх.

(9)

Преобразуем натуральные логарифмы в уравнении (9) в обратные тригонометрические функции [6] и произведём замену Вут/ = В-ут,:

Arcsh

' Ах + ВуП Л V х - хт,/ ]

= 1п

Ах + Вут, + д (Ах + Вут,)2 + (х - хт, )'

х - хт.,

Arcsh

С \

у - ут,

V х - хт,/ ]

= 1п

у - ут/

+ У(у -ут/)2 +(х -хт,

х - хт/

+1п(х - хт/);

+1п(х - хт/ )

(10)

(11)

Подставим (10) и (11) в (9) и получим

х/+1 Г (Ах + В . У| х/+1 Г Г у у- ЛЛ

Arcsh

ут/

V х - хт/ ))

йх - I

Arcsh

у/ - ут/

.. V V х - хт/ ))

(12)

Для решения первого интеграла в уравнении (12) произведём замену: 1=х-хти Ж=йх, М=А-хт/+Вут/ и преобразуем пределы интегрирования. Решением интеграла будет

+

I "1/

I Arcsh

т,

М

Ах + В,.

ут/

V х хт/ )

йх = \t\Arcsh

'At + М Л

t )

л/Т+А

1п

+

2(АМ +1 + А 2t) л/1 + А2

х/+1 - хт,

(13)

2

)

х -х

х -х

\t\Arcsh М

+ '

Преобразуем (13) к виду:

( At + М Л

+

V t у (

л/1 + А2

+ А ■ М +1 + А2Х)+ 1п

лЛ

л/1 + А7

(14)

00

Преобразуем подкоренное выражение из уравнения (14):

((1 + А2 )■ t2 + М2 + 2 А ■ Mt )■ (1 + А2) = А2М2 + 2 AMt + 2 АЪМХ +12 + + 2 A2t2 + А4Х2 + М2.

(15)

Часть выражения (15) можно преобразовать в квадрат суммы: А 2М 2 + 2AMt + 2 А3 МХ +12 + 2 А 212 + А 4t2 = (аМ +1 + А 21 )2. Тогда подкоренное выражение в уравнении (14) будет иметь вид:

((1 + А2 )■ Х 2 + М2 + 2 А ■ М )■(! + А2 )=(АМ + Х + А2Х )2 + М2

(16)

(17)

Уравнение (14) преобразуем к виду:

( (

\t\Arcsh\ АХ + М | +

М

1 + А2

1п

(ам + Х + А 2 Х )2

¡АМ + Х + А 2 ХI + М 2 + АМ + Х + А 2 Х

2

л Л

( 2 Л

+

+ 1п

V V v1 + А2 0

(18)

Для преобразования натурального логарифма к обратной гиперболической функции приведём уравнение (18) к следующему виду:

( ( Г/-^--Л Л

и, , ( АХ + М \t\Arcsh\-I +

М

1 + А2

1п

(ам + Х + А 2 Х) + М 2 + АМ + Х + А 2 Х

М

V

+

+ 1п

V V

( 2М Л

1 + А2

(19)

1п

Заменим часть уравнения (19) на

( \2

АМ + Х + А2Х " + М2 + АМ + Х + А2Х

__:_

М

= Arcsh

( АМ + Х + А2ХЛ М

2

Х

Х

Окончательно после подстановки (20) в (19) с учётом обратных преобразований получим

( ( I ^ 2 Л Л х'+1- хТ,

+

\t\Arcsh

( At + АхТ, + Бш Л t

+

АхТ + БуТ1

1 + А

2

Arcsh

(

а(АХп + Бут.)+1 + At

АхТ + БуТ1

+ 1п

V V

2(Ахт1 + Бут.^

1 + А

2

хг хТг

(21)

Решение второго интеграла уравнения (12) будет иметь вид:

х,+1

х,+1 | Arcsh

х,

Г \

У, - Ут,

V х - хТ, 0

х — хт, \Arcsh

г \

У, - Ут,

+

+ 1 У, - Ут,^™^

V х - хт, ) х — хт, V У, - Ут, )

(22)

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

С,1 = Iх,+1- ^М™^

Ах,+1 + БУТ! Л Ахт, + БУ,

Arcsh

, I а(АхТ, + БУТ1 )+ х,+1 - хТ, + А 2 (х,+1 - хТ, )

1 + А2

Ахт, + БУт

- х, -Xтj\Arcsh

Ах, + Б

Ут,

Ахт, + БУТ1

1 + А2

(

Arcsh

а(ахт, + БУт. )+ х, - хт, + А2 (х, - хт,)

Ахт, + БУТ1

х,+1 хТ,

-I |х,+1 - хт, \Arcsh |+| У, - Ут,1А^ I х+1 - хт' 11+1 |х, - хт,\А^ I ^^ |+| У, - УТг\А^ I х -

У, - Ут,

х, хт,

VУ, - Ут, )) .

(23)

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

нений, позволяющих получить плотность

В качестве примера проведён расчёт

поверхностного заряда (потенциал) по- потенциалов для плоского элемента кон-

верхности. При другом методе триангуля- струкции КА (солнечной батареи). Полу-

ции поверхности (аппроксимации поверх- чена следующая матрица собственных

ности треугольниками другой формы) [7] ёмкостей:

2.407 0.960 0.519 0.365 0.508 0.724 0.367 0.347

0.960 2.407 0.724 0.508 0.365 0.519 0.347 0.367

0.508 0.724 2.407 0.960 0.353 0.532 0.508 0.724

0.365 0.519 0.960 2.407 0.264 0.350 0.365 0.519

0.519 0.365 0.347 0.264 2.407 0.960 0.519 0.365

0.724 0.508 0.532 0.353 0.960 2.407 0.724 0.508

0.367 0.347 0.519 0.365 0.508 0.724 2.407 0.960

0.347 0.367 0.724 0.508 0.365 0.519 0.960 2.407

х

t

t

В результате получаем следующее распределение потенциала поверхности плоского элемента КА (рис. 2).

2М/

/5.32 ~2 /

/',.22

Рис. 2. Распределение собственных потенциалов элементов поверхности КА

В примере при расчёте шаг сетки составлял 0,05 м и суммарный потенциал на пластинке составил 32,8 В. При пересчёте на реальную солнечную батарею, например с площадью 2 м2, потенциал составит 3280В. Расчёт был проведён для условий функционирования КА в диапазоне высот 550-600 км.

Выводы

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

2. Данная методика позволяет определить потенциал поверхности не только КА, но и рассчитать собственную ёмкость бесконтактных ёмкостных датчиков или определить потенциал исследуемых в ускорители частиц образцов материалов.

3. Измеренные потенциалы поверхности МКА «АИСТ» имеют хорошее совпадение с результатами, полученными в ходе моделирования.

Библиографический список

1. Графодатский, О.С. Взаимодействие спутников связи с окружающей средой [Текст] / О.С. Графодатский, Ш.Н. Исляев. - Томск: МГП «РАСКО», 1993. -208 с.

2. Моделирование распределения электромагнитного поля при электростатическом разряде на поверхности космического аппарата [Текст] / Н.Д. Семкин, В.Л. Балакин, Р.А. Помельников [и др.] // Вестн. СГАУ, - Самара. - 2012. - №2(33).

- С.112-119.

3. Методы защиты негерметичных отсеков космического аппарата от проникновения электромагнитного поля [Текст ] / Н.Д. Семкин, В.Л. Балакин, Р.А. Помельников [и др.] // Вестн. СГАУ, Самара. - 2012. - №2(33). - С.120-125.

4. Модель космоса [Текст]: в 2 т. / под ред. М.И. Панасюка, Л.С. Новикова. -Т. 2: Воздействие космической среды на материалы и оборудование космических аппаратов. - М.: КДУ, 2007. - 1143 с.

5. Иоссель, Ю.Я. Расчет электрической ёмкости [Текст] / Ю.Я. Иоссель, Э.С. Кочанов, М.Г. Струнский. - Л.: ЭНЕРГОИЗДАТ. - 1981. - 285 с.

6. Янпольский, А.Р. Гиперболические функции [Текст] / А.Р. Янпольский. -М.: Гос. изд-во физ.-математ. лит. - 1960.

- 195 с.

7. Галанин, М.П. Разработка и реализация алгоритмов трехмерной триангуляции сложных пространственных областей: прямые методы [Электронный ресурс]/ www.keldysh.ru/papers/2006/prep10/ ргер2006 10.html

METHOD OF CALCULATING SELF-CAPACITANCE IN PLANAR STRUCTURAL ELEMENTS OF A SPACECRAFT

© 2013 N. D. Syomkin, V. V. Bragin, R. A. Pomelnikov, V. L. Balakin

Samara State Aerospace University named after academician S.P. Korolyov (National Research University)

The paper presents a methodology for calculating distribution of potentials based on the triangulation method of plane approximation, surface method and analytical determination of the matrix elements of Coulomb interaction. The results obtained demonstrate good agreement with the experiments on the electrization process in small spacecraft «AIST».

Triangulation, self-capacitance, method of integral equations, spacecraft, electrization.

Информация об авторах

Сёмкин Николай Данилович, доктор технических наук, профессор кафедры радиотехники и медицинских диагностических систем, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). E-mail: [email protected]. Область научных интересов: системы управления движением, физика высокоскоростного удара, космические исследования.

Брагин Василий Владимирович, аспирант кафедры радиотехники и медицинских диагностических систем, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). E-mail: bragin [email protected]. Область научных интересов: электромагнитная совместимость, физика космической плазмы.

Помельников Роман Александрович, младший научный сотрудник, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). E-mail: r [email protected]. Область научных интересов: космическая физика, масс-спектрометрия, электромагнитная совместимость.

Балакин Виктор Леонидович, доктор технических наук, профессор кафедры космического машиностроения, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). E-mail: [email protected]. Область научных интересов: динамика полёта и управление движением летательных аппаратов, космическое машиностроение.

Syomkin Nikolay Danilovfch, doctor of technical science, professor, the department of radio engineering and medical diagnostic systems, Samara State Aerospace University named after academician S.P. Korolyov (National Research University). E-mail: [email protected]. Area of research: high-speed impact, mass spectrometry, space device engineering.

Bragin Vasily Vladimirovich, postgraduate student, the department of radio engineering and medical diagnostic systems, Samara State Aerospace University named after academician S.P. Korolyov (National Research University). E-mail: bragin [email protected]. Area of research: electromagnetic compatibility, space plasma physics.

Pomelnikov Roman Alexandrovich, junior research associate, Samara State Aerospace University named after academician S.P. Korolyov (National Research University). E-mail: r [email protected]. Area of research: space physics, mass-spectrometry, electromagnetic compatibility.

Balakin Victor Leonidovich, doctor of technical science, professor, the department of space engineering, Samara State Aerospace University named after academician S.P. Koroly-ov (National Research University). E-mail: [email protected]. Area of research: flight dynamics, flying vehicle motion control, space engineering.