Страница 32: ОСТ 134-1044-2007. Аппаратура, приборы, устройства и оборудование космических аппаратов. Методы расчета радиационных условий на борту космических аппаратов и установления требований по стойкости радиоэлектронной аппаратуры космических аппаратов к воздействию заряженных частиц космического пространства естественного происхождения (81742)

Изменение крупномасштабных магнитных полей гелиосферы при­нимают пропорциональным изменению полярного магнитного поля Солнца, величина и полярность которого принимаются зависящими от уровня СА и четности-нечетности цикла СА( ВЛ)

где S=1 при t-t* >0 (иначе S=-l);

tj есть момент перемены знака полярного магнитного поля в ом цикле солнечной активности, отождествляемый с моментом максиму­ма СА;

/7

. W^mm- минимальные уровни СА, окаймляющие н-цикл СА.

3. Динамика крупномасштабной модуляции ГКЛ характеризуется эф­фективным модуляционным потенциалом гелиосферы Я_о(/гд) (для частиц с жесткостью R в данный момент времени t), величину которого вычис­ляют по формуле:

- дДо, R, 0]} .= 0.37 + 3 • 10^ ■ [t - дД«, R,t)], (В.2)

где дДи,аО - время запаздывания изменения потоков частиц по сравне­нию с изменением солнечной активности.

4. Время запаздывания изменения потоков частиц ГКЛ по сравнению с изменением солнечной активности - дДи.Лд) принимается зависимым от величины магнитной жесткости частиц (R), от четности или нечетности цикла (п) и от фазы СА:

дДтг, л, г) = о,5[т; +т_ (/г)]+о.5[т₊ - г_(л)]- ^[месяцев] (В.З)

где амплитуда задержки в четных циклах солнечной активности (т₊) не зависит от жесткости частиц:

Т₊ = 15[месяцев], (В .4)

а в нечетных циклах (Т_):

Г_(я) - 15-/Г⁰’⁴⁵[месяцев]. (В.5)

Временная функция изменения величины задержки из (В.З) принима­ется в виде

(B.6)

■jwfn п

где 6_wt= 16 месяцев.

5. В модели вычисляются дифференциальные жесткостные и энерге­тические спектры потоков частиц.

Связь между значениями жесткости R (ГВ) и значениями кинетиче­ской энергии частиц Е (ГэВ/нуклон) вычисляют по формулам - для заданных значений жесткости R частиц

Гт

» (В.7)

Л| J

где - массовое число и заряд частицы;

т, - масса покоя частицы - m_t= 0.51 МэВ - для электронов, т_р= 0.938 ГэВ - для протонов, m_z>₂ = 0.939 ГэВ/нуклон для ТЗЧ;

- для заданных значений кинетической энергии частиц Е

r = ^-^e(e+2^) (В.8)

1. Дифференциальные жесткостные спектры потока частиц ГКЛ Ф,(^д) (см2 _{с С}р ГВ)-1 ддя частиц жесткости R в момент времени t вычис­ляют по формуле

ФДМ-^' > (В.9)

где C_ira_it^ - параметры немодулированного жесткостного спектра частиц сорта ”Г, приведенные в табл. В.1;

“ безразмерный параметр, определяемый по формулеZ /Й?

А, (Я, 0 = 55+ШгА И”'.") ’7^) “1

где р - отношение скорости частицы к скорости света, определяемое формулой

2. Средние квадратичные отклонения расчетных значений потока час- ТИЦ вычисляют по формуле

0

-її

Л

1+ т~

1. Дифференциальные энергетические спектры частиц ГКЛ E(E.t)

(см2 _Ср _с ГэВ)~1 для частиц энергии Е в момент времени t вычисляют по формуле

(В.13)

F, (£,/)=Ф₍(й,г)Д—

(см² ср с МэВ)^_1 - для частиц с (протонов);

(см² ср с МэВ/нуклон)"1 для частиц с л, >2 (ионов).

2. Средние квадратичные отклонения расчетных значении рассчитывают по формуле

r

(ВЛ4)

V'c.f

Таблица В. 1 - Параметры жесткостных

спектров 34 ГКЛ

Приложение Г
(рекомендуемое)
Методы расчета геоцентрических и L-B координат

Г, 1 Метод расчета геоцентрических координат

Массу КА принимают ничтожной по сравнению с массой притяги­вающего тела. Движение КА рассматривают как невозмущенное, кепле­ровское движение материальной точки, которое происходит под действи­ем только одной центральной силы притяжения.

Связь абсолютной геоцентрической (х,у,^) и орбитальной систе­мами координат иллюстрирует рис. Г.1 из [4].

Вводится понятие эксцентрической аномалии Е, которая связана с ис­тинной аномалией .9 следующим соотношением где е - эксцентриситет орбиты.

В этом случае невозмущенное движение КА будет полностью опи­сываться уравнением Кеплера

£ - е sin Д - , (Г.2)

а

где д - гравитационная постоянная;

М - масса Земли;

і ~ время движения КА, отсчитываемое с момента прохождения его точки перигея;

г_п,г_а— расстояния от центра Земли до перигея и апогея орбиты, со­ответственно;

2a

x,y,z - оси абсолютной геоцентрической системы координат (ось X направлена в точку весеннего равноденствия Ч^у);

линия узлов - линия пересечения плоскости орбиты с плоскостью эк­ватора;

£2 -долгота восходящего узла~угол между положительным направле­нием оси х и направлением линии узлов в восходящий узел, при этом восходящий узел определяется как узел орбиты, который проходит КА, двигаясь из области отрицательных аппликат в область положительных. Угол £2 изменяется от 0° до 60°;