Определение времени выживания космической тросовой системы

Автор: Асланов В.С., Пироженко А.В., Волошенюк О.Л., Кислов А.В., Ящук А.В.

Журнал: Известия Самарского научного центра Российской академии наук @izvestiya-ssc

Рубрика: Механика и машиностроение

Статья в выпуске: 4-1 т.12, 2010 года.

Бесплатный доступ

Предложены алгоритмы определения вероятности выживания космической тросовой системы при угрозе ее разрушения космическими частицами. Потоки частиц космического мусора рассмотрены как стационарные пуассоновские потоки. Интенсивность потоков частиц рассчитывается с помощью известных компьютерных моделей ORDEM 2000, Master 2001. На основе этого можно определены вероятностные оценки времени выживания троса и, выбраны его конструктивные характеристики.

Космические тросовые системы, космическим мусором, оценки времени выживания, модели загрязнения околоземного космического пространства

Короткий адрес: https://sciup.org/148199348

IDR: 148199348

Текст научной статьи Определение времени выживания космической тросовой системы

слоев атмосферы, проведения распределенных измерений [8-10].

Одной из основных особенностей функционирования КТС является их бьльшая уязвимость, в сравнении с КА, при столкновении с космическими частицами (космическим мусором и метеорными частицами). Маленькие размеры и большое количество таких частиц не позволяют использовать известные методы баллистики для описания их движения и оценки опасности столкновения с КТС. Имея большую скорость движения, частицы, размеры которых составляют ~ (1/3 – 1/2) толщины троса, могут перебить трос и тем самым разрушить систему. Поэтому использование КТС предполагает решение проблемы выживания тросовых соединений в условиях космического пространства. Здесь и далее, говоря о выживании КТС, мы будем подразумевать выживание только троса при угрозе его разрушения космическими частицами.

В работах [11-13] определены оценки вероятности выживания КТС при столкновении с космическими частицами. Эти оценки основываются на определении количества столкновений КТС с суммарным потоком частиц разных размеров на некотором интервале времени. Наиболее проработанным представляется алгоритм, предложенный в работе [11], позволяющий оценить влияние размера частиц на выживаемость КТС. Но данный алгоритм представляется излишне усложненным и позволяет определять вероятность выживания КТС только для монолитных конструкций троса. Вместе с тем, изменения конструкции троса позволяют существенно увеличить время его выживания [14].

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

ОПРЕДЕЛЕНИЕ ВРЕМЕНИ ВЫЖИВАНИЯ КТС С МОНОВОЛОКОННЫМ ТРОСОМ

Пусть трос состоит из одного волокна (рис. 1), которое имеет круглое поперечное сечение диаметром DT . Предполагается также, что частицы космического мусора (КМ) имеют сферич еск ую форму с некоторым диаметром d i , i = 1, p ( P — заданное число потоков, в зависимости от размера (диаметра) частиц, и каждому i -му потоку частиц соответствует среднее значение диа ме тра di ). Считаем, что d i d 2 ... dp , i = 1, p .

Не каждое столкновение частиц КМ с тросом приводит к разрушению троса. Определение условий, при которых в результате столкновения будет происходить разрушение троса, является отдельной задачей. Так же, как в [11], будем предполагать, что пр и с толкновении с частицей диаметром d i , i = 1, p трос не будет перебит, если расстояние прохождения частицы от края троса не больше, чем x . По экспериментальным данным [15] значение x может изменяться в преде-

лах

3 d t ;2 d t

, в зависимости от материала

троса. Исходя из этого условия, опасность для троса будут представлять потоки частиц КМ , средние размеры которых d i x , i = k , p , k = max { i : d i x } .

Пусть § - положение геометрического центра частицы КМ диаметра d i , i = k , p относительно левого края троса (рис. 1). В случае, когда § удовлетворяет условию

— + x § DT

2              T

d x + i

2 ,

i = k , p — трос бу-

дет перебит.

Пусть X i , i = k , p величина, равная интенсивности потока частиц КМ диаметром di , раз-

рушающих трос в единицу времени. Величина X i определяется формулой

X,=^F,St(di), i = k,p, где AFi — среднее число частиц диаметром di, проходящих за единицу времени через единичную площадку (определяется из модели загрязнения околоземного пространства);

S T ( d i ) — приведенная площадь продольного сечения троса для i- го потока частиц КМ, равная

S t ( d i ) = L ( D T + d i - 2 x ) , i = k J ,

L — длина троса.

Тогда, интенсивность X суммарного потока частиц КМ, разрушаю щи х трос, определяется как сумма всех X i , i = k , p , т.е.

p

X = E X i , i = k , p .

i = k

Будем считать, что поток частиц КМ является стационарным пуассоновским потоком, т.е. удовлетворяет следующим условиям:

  • -    вероятность попадания того или иного числа частиц космического мусора в трос не зависит от его положения на орбите, иными словами, частицы распределены на заданной высоте с одинаковой средней плотностью;

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

  • -    вероятность попадания на малый участок троса КТС двух и более частиц пренебрежимо мала по сравнению с вероятностью попадания одной частицы (это условие означает практическую невозможность совпадения двух и более частиц).

При соблюдении данных условий, число событий на любом фиксированном интервале времени будет распределено по закону Пуассона [16]. Т.е., вероятность того, что на некотором интервале времени t произойдет столкновение с m частицами, разрушающими трос, определяется по формуле

Рис. 1. Модель взаимодействия моноволок онн ого троса и частиц КМ диаметром d i , i = k , p

P m ( t ) =

( X t ) m m !

„- X-t

e ,

t >  0 .

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

Р т ( t ) = Р ( T t ) = Р о ( t ) = e - X - t .     (1)

Функция распределения FT ( t ) времени жизни троса Т определяется равенством

FT ( t ) = Р ( T t ) = 1 - P ( T t ) .     (2)

Или, с учетом (1)

FT ( t ) = 1 - e - х .             (3)

Дифференцируя (3), находим плотность распределения Т f (t) = Xe-"'t, t > 0.

Таким образом, время жизни троса Т – показательно распределенная случайная величина.

Важной характеристикой случайной величины является ее математическое ожидание. В нашем случае математическое ожидание случайной величины Т представляет собой среднее время жизни троса и определяется по формуле mT = M [t]= j tf (t)dt = X J te-X'*dt, 0                     0

интегрируя по частям, получаем mT = —

T X '

Таким образом, — среднее время жизни

X

троса КТС в случае, когда трос является моно-волоконным. Тогда получим оценку времени t0 , на протяжении которого трос с заданной вероятностью p 0 не будет разрушен

In P 0

t 0

X

.

ОПРЕДЕЛЕНИЯ ВРЕМЕНИ ВЫЖИВАНИЯ ДЛЯ ДВУХВОЛОКОННОГО ТРОСА

Рассматривается простейшая модель многоволоконного троса. Предполагается, что трос состоит из двух идентичных волокон (рис. 2) каждый из которых имеет круглое поперечное сечение диаметром DTV . Волокна находятся на некотором расстоянии друг от друга А у . Частицы космического мусора имеют сферическую форму с некоторым диаметром di ( d i x ), i = k , p , при условии, что d 1 d 2 ... dp .

Рис. 2. Модель двухволоконного троса

Для модели двухволоконного троса предполагается, что есть потоки частиц КМ, которые могут:

  • -    перебить только одно из двух волокон;

  • -    перебить трос, т.е. два волокна одновременно.

Исходя из этого предположения, выделяем два суммарных потока частиц КМ, с размерами в диапазонах x d i < А у + 2 x и d i > А у + 2 x , представляющих опасность для троса. Интенсивности этих суммарных потоков определяются ниже.

Интенсивность суммарного потока частиц КМ, которые могут перебить только одно волокно. Как и в случае моноволоконного троса, потоки частиц КМ с размерами в диапазоне x < di < Ау + 2x могут перебить волокно диаметром DTV , если x удовлетворяет условию di                         di

- у + x < 5 < dtv - x + у,   где i = k, j , j = max{i: di < Ау + 2 x}.

Пусть X' i , i = k , j интенсивность потока частиц КМ диаметром di в единицу времени, равная

X 1’ = А FS <"( d , ), i = kJ , где S Х'! d i ) — приведенная площадь продольного сечения волокна заданного диаметра DTV для i- го потока частиц КМ, равная

S i(1) ( d i ) = L ( D tv + d i - 2 x ) , i = kJ .

Потоки частиц КМ, размеры которых находятся в диапазоне di > Ау + 2x , могут перебить также только одно волокно, если x лежит di                                di

на участке — — + x < 5 < DTV +Ау + x —  , i = j +1, p , j = max{i: di < Ау + 2 x}.

Интенсивность X 2 , i = j + 1, p в этом случае определяется формулой

X =^2( di), где

S 1(2) ( d i ) = L ( d tv у ) .

Тогда интенсивность Xi суммарного потока частиц КМ, которые могут перебить только одно волокно, равна jp

X 1 = 2 хц + 2 X 12) . i = k           i = j + 1

Интенсивность суммарного потока частиц, которые перебивают два волокна одновременно. Потоки частиц КМ с размерами в диапазоне d i у + 2 x могут перебить трос диаметром DT = 2 D TV + А у , в том случае, если x лежит на участке

D TV

dd

+ Ay - .2 + x < £ < DTV - x + _L , i — j +1, p , j — max{i: d < Ay + 2 x}.

Интенси вность каждого такого потока частиц 2 2, i , i j + 1, p определяется формулой

2,—AFS2( di), i — j +1, p , где

S 2( d i ) L ( d i -A y - 2 x ) .

Тогда интенсивность 2 2 суммарного потока частиц, которые могут пер ебить т рос, определяется суммой всех 2 2, i , i j + 1, p

p

2 ^ 2 i , i j + 1, p .

i = j + 1

Вероятностные оценки времени выживания для двухволоконного троса. Как и в случае моноволоконного троса считаем, что суммарный поток частиц КМ, которые могут перебить два волокна одновременно и суммарный поток частиц, которые перебивают только одно волокно, представляют собой стационарные пуассоновские потоки [16], с интенсивностями 2 1 и 2 2 соответственно.

Пусть Т — случайная величина, описывающая время жизни двухволоконного троса, P T ( t ) — вероятность того, что время жизни троса больше t . Тогда, исходя из (1) справедливо равенство

Р т ( t ) = P o( t ) ,                (4)

где Po(t ) - вероятность неразрушения троса на интервале времени t. Величина Po(t ) определяется как сумма вероятностей того, что за время t не будет перебито ни одного волокна и, что за время t будет перебито только одно волокно, т.е.

P o ( t ) = ( р ( t ) ) 2 P 2 ( t ) + 2 P ( t )(1 - P , (t )) P 2 ( t ) ,(5) где P1(t ) - вероятность того, что частицами первого суммарного потока (потока частиц с интенсивностью X i ) не будет перебито ни одного волокна;

P 2 ( t ) - вероятность того, что частицами второго суммарного потока (потока с интенсивностью X 2 ) не будет перебит трос.

С учетом (1), P 1 (t ) , P 2( t ) примут вид

P( t) = e -21,(6)

P2( t) = e -22 t.(7

Заметим, что

2 — 21+ 22,(8)

где 2 - это интенсивность суммарного потока всех частиц КМ диаметром d i x , i = p , k , которые могут разрушить одно волокно.

Подставив выражение (5) в (4), и учитывая (6)-(8), получаем

P o( t ) = 2 e - 2 t - e4 2 1 + 2 ) t .            (9)

Как и в случае с моноволоконным тросом, аналогично (2), функция распределения времени жизни троса Т равна

F ( t ) = 1 + e "( 2 1 + 2 ) t - 2 e - 2 t .       (10)

Дифференцируя (10) находим плотность распределения Т f (t) = -(21 + 2) e -(21 + 2) t + 22e -2t.

Среднее время жизни троса, равное математическому ожиданию величины Т , распределенной по показательному закону, определяется по формуле

X

mT

o

-

2 + 2 1

2 2 1 + 2    1 ^      2 1 л

--------— — 1 +-- 2 ( 2 + 2 1 ) 2 (    2 + 2 1 )

где – среднее время жизни троса, если он со- стоит из одного волокна диаметром DTV ;

– увеличение среднего времени 2 ( 2 + 2 1 )

жизни троса, при добавлении еще одного волокна.

В предель но м случае, т.е. при условии, что A y d i , i k , p значения 2 1 и 2 равны, и среднее время жизни двухволоконного троса в 1,5 раза больше в сравнении с моноволоконным тросом (11).

Время t0 , на протяжении которого трос с заданной вероятностью P o не будет разрушен, определяется из следующего равенства

2 e - 2 t o - e ~( 2 1 + 2 ) t o

Р o .

Отсюда to

1     ln( P °)

- 2 " ln(2 - e-2 1)"

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

Оценки интенсивностей потоков частиц КМ на конкретных высотах проводятся на основе расчетных моделей, описывающих метеорные условия и распределение космического мусора. Наиболее известные и используемые компьютерные модели– ORDEM 2000 и Master 2001, разработанные в NASA и ESA соответственно [17,18]. Эти модели созданы на основе данных, полученных при исследовании космического пространства (приборами, устанавливаемыми на метеорологических ракетах, спутниках, и т.д.) и в результате наземных измерений (радиолокационными и оптическими средствами).

Модель ORDEM 2000 применима для исследования околоземного пространства на высотах от 200 км до 2000 км. Интервал прогнозирования: 1991-2030 гг. [17].

Данная модель позволяет определять пространственную плотность космических частиц (число частиц в единице объема, 1/км3), их среднюю орбитальную скорость (км/с) и значения потоков частиц (число столкновений за единицу времени с единичной площадкой, 1/м2) на каждом сегменте заданной орбиты, а также усредненные значения полного потока для каждого заданного размера частиц на всем диапазоне высот. Исходными данными являются параметры орбиты (высота апогея и перигея, наклонение, аргумент перигея), количество сегментов Nseg , ( Neg e [1 , 100] ), на которое будет разделена плоскость орбиты, год наблюдения и, если наблюдение проводится с Земли географическая широта, угол азимута и высоты. Все характеристики потоков определяются для частиц размером от 10 мкм до 1 м (для шести фиксированных размеров: 10 мкм, 100 мкм, 1 мм, 1 см, 10 см, 1 м).

Модель Master 2001 (ESA Meteoroid and Space debris Terrestrial Environment Reference) применима для исследования околоземного пространства вплоть до высот геостационарной орбиты [18]. Интервал прогнозирования: 19602050 гг. В Master 2001, в отличие от ORDEM 2000, плотность потока — это поток частиц через поверхность некоторого объекта сферической формы с единичной площадью поперечного сечения (1 м2).

В сравнении с ORDEM 2000 модель Master 2001 позволяет более точно:

  • -    задавать положение орбиты КА (дополнительными параметрами являются долгота восходящего узла и истинная аномалия);

  • -    задавать интервал времени (с точностью до дня);

  • -    учитывать более широкий диапазон частиц (диаметр (м) и масса (кг) частиц могут изменяться в пределах [10-6, 100]).

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

ВЫВОДЫ

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

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

В силу отмеченного, представляется, что предложенные алгоритмы дают несколько завышенные оценки вероятности разрушения троса частицами космического мусора. Тем не менее, они позволяют оценивать вероятность времени выживания троса, и, выбирать конструктивные параметры троса в соответствии с задачей создания КТС. Например, во многих задачах экспериментальных исследований КТС предполагается использовать трос длиной до 1000 м [7,9,21,22]. Продолжительность экспериментов — до месяца. Полученные оценки, показывают, что для неразрушения троса частицами космического мусора с вероятностью большей 0.98, достаточно взять моноволоконный трос диаметром не менее 1 мм.

Работа выполнена при финансовой поддержке РФФИ (проект 09-01-00384а).

Список литературы Определение времени выживания космической тросовой системы

  • Lorenzini E. C., Cosmo M.L. Mission analysis of spinning systems for transfers from low orbits to geostationary//Journal of Spacecraft and Rockets. 2000. V.37. № 2. P.165-172.
  • Hoyt R.P., Uphoff C.W. Cislunar Tether Transport System//Journal of Spacecraft and Rockets. 2000. V. 37. № 2. P. 177-186.
  • Nordley G.D., Forward R.L. Mars Earth Rapid Interplanetary Tether Transport System: Initial Feasibility Analysis//Journal of Propulsion and Power. 2001. V. 17. P.499-507.
  • Rapid Interplanetary Tether Transport/R.P. Hoyt, R.L. Forward, G.D. Nordley, C.W. Uphoff//50th International Astronautical Congress, Netherlands, Amsterdam, Oct 1999. IAC 99 A.5.10.
  • Terminator Tether: A Spacecraft Deorbit Device/R. L. Forward, R. P. Hoyt, C. W. Uphoff//Journal of Spacecraft and Rockets. 2000. V. 37. № 2. P. 187-196.
  • Orbital Maneuvering with Spinning Electrodynamics Tethers/J. Pearson [et al.]//2nd International Energy Conversion Engineering Conference, Rhode Island, 16 -19 August 2004. AIAA 2004 5715.
  • Электродинамическая тросовая система увода космических аппаратов с орбит: исследование на наноспутниках/А. П. Алпатов, Ф. Н. Гребенкин, А. В. Мищенко, А. В. Пироженко//Вiсник днiпропетровського унiверситету. 2006. № 2/2. С.5-10.
  • Белецкий В.В., Левин Е.М. Динамика космических тросовых систем. М.: Наука, 1990. 329 с.
  • Ротационное движение комических тросовых систем/А. П. Алпатов [и др.] Днепропетровск: Институт технической механики НАН Украины и НКА Украины, 2001. 404 с.
  • Pradeep S., Kumar К. Extension of tethered satellites in the atmosphere//Acta Astronautica. 2003. V. 52. P. 1-10.
  • Anz-Meador P. D. Tether Debris Interactions in Low Earth Orbit//Proceedings of the American Institute of Physics Conference. 2001. V. 552. № 1. P. 525-531.
  • Benefits and risks of using electrodynamics tethers to de orbit spacecraft/C. Pardini, T. Hanada P. Krisko//57th International Astronautical Congress, 2006. IAC 06 B6.2.10.
  • Forward R.L. Retrive tether survival probability//38th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, Indiana, Indianapolis, 7 10 July 2002. AIAA 2002 4047.
  • Волошенюк О.Л., Храмов Д.A. Оценка вероятности выживания космических тросовых систем при столкновении с частицами космического мусора//Техническая механика. 2008. № 1. C. 3-12.
  • Lorenzini E.C., Cosmo M.L. Tethers in Space Handbook. 3rd edition. Smithsonian Astrophysical Observatory, 1997. 241 p.
  • Вентцель Е.С. Теория вероятностей. М.: Высш.шк., 1998. 576 с.
  • The New NASA Orbital Debris Engineering Model ORDEM2000/Lion J. C., Matney M. J., Anz-Meador P. D., Kessler D., [et al.]. Johnson Space center, Houston, Texas, USA, 2002. NASA/TP 2002 210780.
  • The MASTER 2001/P. Wegener, J. Bendisch, [et al.]//J. Advances in Space Research. 2004. V. 34, № 5. P. 959-968.
  • Wallace B.K. SEDS tether deployment ground test. Washington, April, 1995. P. 653-668.
  • Attitude and Orbit Determination of a Tethered Satellite System/K. T. Alfriend, W. J. Barnds, S. L. Coffey, L. M. Stuhrenberg//AAS/AIAA Astrodynamics Specialist Conference. Halifax, Nova Scotia, Canada, 14-17 August 1995.
  • Technology Demonstrator of a Standardized Deorbit Module Designed for CubeSat and RocketPod Applications/Voronka N.R. [et al.]//19th Annual AIAA/USU Conference on Small Satellites, India, Logan, August 2005. Logan, 2005.
  • DTUSat 1 Homepage. URL: http://dtusat1.dtusat.dtu.dk (дата обращения 10.09.2009).
Еще
Статья научная