Адрес документа: http://law.rufox.ru/view/25/1200050146.htm


     РД 50-25645.216-90

Группа Ф40

     
     
РУКОВОДЯЩИЙ НОРМАТИВНЫЙ ДОКУМЕНТ

МЕТОДИЧЕСКИЕ УКАЗАНИЯ

Безопасность радиационная экипажа космического
аппарата в космическом полете

МЕТОД РАСЧЕТА РАСПРЕДЕЛЕНИЯ ПОГЛОЩЕННОЙ И ЭКВИВАЛЕНТНОЙ ДОЗ
КОСМИЧЕСКИХ ИЗЛУЧЕНИЙ ПО ТОЛЩИНЕ МАТЕРИАЛОВ НА ВНЕШНЕЙ
ПОВЕРХНОСТИ КОСМИЧЕСКОГО АППАРАТА НА ОРБИТАХ,
ПРОХОДЯЩИХ ЧЕРЕЗ ЕРПЗ

     
     
ОКСТУ 6968

Дата введения 1991-07-01

     
     
ИНФОРМАЦИОННЫЕ ДАННЫЕ

     
     1. РАЗРАБОТАН И ВНЕСЕН Минздравом СССР
     
     РАЗРАБОТЧИКИ
     
     Н.А.Анфимов, член-корр. АН СССР; В.В.Архангельский; В.Н.Васильев, канд. техн. наук; А.А.Волобуев; В.А.Гончарова; А.И.Григорьев, д-р мед. наук; В.Е.Дудкин, д-р физ.-мат. наук; Е.Е.Ковалев, д-р техн. наук; В.В.Козелкин, д-р техн. наук; Е.Н.Лесновский, канд. техн. наук; В.Г.Митрикас, канд. физ.-мат. наук; В.А.Панин, Е.В.Пашков, канд. техн. наук; В.М.Петров, канд. физ.-мат. наук; Ю.В.Потапов, канд. физ.-мат. наук; В.А.Шуршаков, канд. физ.-мат. наук
     
     2. УТВЕРЖДЕН И ВВЕДЕН В ДЕЙСТВИЕ Постановлением Государственного комитета СССР по управлению качеством продукции и стандартам от 28.03.90 N 660
     
     3. Срок первой проверка - 1996 г., периодичность проверки - 5 лет
     
     4. ВВЕДЕН ВПЕРВЫЕ
     
     5. ССЫЛОЧНЫЕ НОРМАТИВНО-ТЕХНИЧЕСКИЕ ДОКУМЕНТЫ
     

Обозначение НТД, на которые дана ссылка

Номер пункта, приложения

ГОСТ 25645.138-86

1.3, приложение 2

ГОСТ 25645.139-86

1.3, приложение 2

РД 50-25645.206-84

1.3

НРБ-76/87*

1.3

_______________
     * Действуют СП 2.6.1.758-99 (НРБ-99), здесь и далее по тексту. - Примечание изготовителя базы данных.

      
     
     Настоящие методические указания устанавливают метод расчета поглощенной и эквивалентной доз от протонов и электронов естественного радиационного пояса Земли (ЕРПЗ) за защитой толщиной до 1,5 г/см (в массовых единицах длины), выполненной из материалов с зарядом 15 при космических полетах на высотах от 200 до 10 км при разных наклонениях орбиты к плоскости экватора.
     
     Методические указания предназначены для расчетов дозовых нагрузок на биологические и технические объекты, в том числе на космонавтов при их работе вне космического аппарата.
     
     

1. ОБЩИЕ ПОЛОЖЕНИЯ

     
     1.1. Метод расчета доз от протонов основан на предположении прямолинейного распространения и непрерывного торможения протонов в веществе. Ослабление первичного потока протонов за счет ядерных взаимодействий не учитывают.
     
     1.2. Метод расчета доз от электронов основан на рассмотрении многократного рассеяния в веществе с учетом ионизационных потерь, пренебрегая их флуктуациями. Вклад в дозу от тормозного излучения не учитывают.
     
     1.3. Входная информация, необходимая для проведения расчетов, должна содержать:
     
     - дифференциальные энергетические спектры протонов и электронов  и , вычисляемые в соответствии с ГОСТ 25645.138, ГОСТ 25645.139;
     
     - заряд  и массовое число  для вещества защиты;
     
     - значения пробегов  и ионизационных потерь  в веществе защиты и веществе-поглотителе, определяемые для протонов в соответствии с РД 50-25645.206 или приложением 1, для электронов - в соответствии с приложением 1. Разрешается использовать для определения  и  аппроксимационные выражения вида:
     
     а) протоны (для защиты из алюминия):     

    

     б) протоны (для тканеэквивалентного вещества):     

    

  

     в) электроны (для защиты из алюминия):
     

;

;

     
     г) электроны (для тканеэквивалентного вещества):
     

;    ,

     
где  в МэВ,  в г/см;

     
     - зависимость коэффициента качества  от ионизационных потерь  для протонов, определяемую аппроксимацией регламентированных в НРБ-76/87 нормативных значений, в виде:
     

     
     
2. АЛГОРИТМ РАСЧЕТА ПОГЛОЩЕННОЙ И ЭКВИВАЛЕНТНОЙ ДОЗ
ОТ ПРОТОНОВ И ЭЛЕКТРОНОВ

     
     2.1. За исходные данные принимают нормальное падение широкого пучка частиц на плоский полубесконечный слой толщиной . Защита выполнена из материала с 15, толщина защиты не превышает 1,5 г/см. Дозу определяют в точке вещества-поглотителя, располагающегося непосредственно за защитой.
     
     2.2. Алгоритм расчета доз от протонов
     
     2.2.1. Поглощенную дозу в веществе защиты от протонов, имеющих энергетический спектр , за защитой толщиной  вычисляют по формуле
     

,                                               (1)

     
где  - ионизационные потери протонов в веществе защиты, МэВ·г·см;

     
      - коэффициент перехода от поглощенной энергии к дозе, равный 1,6·10 Гр·МэВ·г;
     
      - величина поглощенной дозы, Гр;
     
      - энергия протонов на глубине защиты , связанная с энергией протонов, падающих на защиту , соотношением "пробег-энергия"
     

,                                                    (2)

     
где  и  - ионизационные пробеги протонов с энергиями  и , соответственно, в веществе защиты;
     
     0,1 МэВ.
     
     2.2.2. Поглощенную дозу в тканеэквивалентном веществе определяют по формуле
     

,                                  (3)

     
 где  - ионизационные потери протонов в тканеэквивалентном веществе.
     
     2.2.3. Эквивалентную дозу от протонов за плоским слоем толщиной  вычисляют по формуле
     

.                             (4)

     
     2.2.4. Для расчета доз от протонов по формулам (1-4) задают расчетную сетку изменения энергий, равномерную в логарифмическом масштабе: 1,122; 1,414; 1,778; 2,239; 2,818; 3,548; 4,467; 5,623; 7,079; 8,913 МэВ в каждом порядке изменения значений . При этом интервалы изменений  равны: 0,259; 0,326; 0,410; 0,517; 0,650; 0,819; 1,031; 1,298; 1,633; 2,057 МэВ.
     
     2.2.5. Вычисляют значения доз, заменяя интегрирование суммированием по .
     
     2.2.6. Переходят к другим слоям защиты и (или) тканеэквивалентного вещества.
     
     2.3. Алгоритм расчета доз от электронов
     
     2.3.1. Определяют  и , где  - дифференциальный энергетический спектр электронов ЕРПЗ в интервале энергий электронов 0,04 МэВ и 4,0 МэВ.
     
      2.3.2. Задают полное число статистических испытаний  (от 1000 до 5000 историй), присваивают начальные значения накопителей для статистической оценки дозы 0 и квадрата дозы 0, присваивают начальное значение номеру текущего испытания 0.
     
     2.3.3. Определение начальных значений величин для -го статистического испытания
     
     2.3.3.1. Присваивают начальные значения:
     
     - глубине проникновения электрона в пластину 0;
     
     - косинусу угла падения электрона  на пластину, 1 (нормальное падение).
     
     2.3.3.2. Определяют начальное значение энергии падающего электрона  (метод Ноймана):
     
     - генерируют , здесь и далее  - случайные числа, равномерно распределенные в интервале [0; 1];
     

   - вычисляют , МэВ,   

 (5)

                                                                            
где  и  - минимальное и максимальное значения энергии в спектре электронов;
     
     - если , то переходят к генерации  и .
     
     Начальной энергии электрона  присваивают значение ; .
     
     2.3.3.3. Полагают 1.
     
     2.3.4. Движение электрона через слои вещества в пластине
     
     2.3.4.1. Вычисляют параметр экранирования Мольер
     

,                        (6)

     
где .

     
     2.3.4.2. Определяют средний пробег электрона до рассеяния:
     

, г/см,                                       (7)

     
где 6,022·10, 1/моль;

     
,                            (8)

     
где 2,818·10 см.

     
     2.3.4.3. Определяют пробег электрона до рассеяния
     

, г/см.                                                   (9)

     
     2.3.4.4. Рассчитывают пробег электрона с энергией
     

, г/см,                                                     (10)

     
где  - аппроксимирующая функция, связывающая энергию электрона с его пробегом в веществе защиты.
     
     2.3.4.5. Рассчитывают координату точки рассеяния :
     

, г/см.                                        (11)

     
      2.3.4.6. Рассчитывают энергию электрона в точке рассеяния
     

, МэВ,                                           (12)

     
где  - функция, обратная .

     
     Если , перейти к п.2.3.3.2.
     
     2.3.4.7. Рассчитывают косинус угла рассеяния:
     

.                                          (13)

     
     2.3.4.8. Рассчитывают косинус угла между импульсом электрона после рассеяния и нормалью к пластине
     

,                                (14)

     
где .

     
     2.3.4.9. Если 0, то переходят к п.2.3.3.2.
     
     2.3.4.10. Если , то переходят к п.2.3.3.2.
     
     2.3.4.11. Если 0, то присваивают ;
     
     ;  и переходят к п.2.3.4.
     
     2.3.5. Определение вклада в поглощенную дозу от -го испытания
     
     2.3.5.1. Вычисляют
     

, Гр,                                            (15)

     
где  - аппроксимирующая функция, связывающая энергию электрона и ионизационные потери в веществе-поглотителе.
     

;                                                           (16)

     
.

     
     2.3.5.2. Если , то переходят к п.2.3.3.
     
     2.3.6. Определяют средние значения поглощенных доз
     

, Гр                                                       (17)

     
и средние квадратичные отклонения

     
.                                             (18)

     
     2.3.7. Поглощенную (эквивалентную) дозу от спектра электронов, распределенных по закону  в интервале  вычисляют по формуле
     

, Гр.                                             (19)

     
     2.4. Примеры расчета поглощенной и эквивалентной доз от протонов и электронов ЕРПЗ приведены в приложении 2.
     
     

ПРИЛОЖЕНИЕ 1
Справочное

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

     
1. Протоны

     
Таблица 1

     

Энергия, МэВ

Алюминий

Тканеэквивалентное вещество

,
МэВ·г·см

,
г·см

,
МэВ·г·см

,
г·см

1.120-01

4.377+02

2.941-04

8.793+02

1.514-04

1.410-01

4.132+02

3.623-04

8.267+02

1.854-04

1.770-01

3.826+02

4.530-04

7.631+02

2.308-04

2.240-01

3.468+02

5.823-04

6.897+02

2.958-04

2.810-01

3.104+02

7.565-04

6.159+02

3.835-04

3.540-01

2.809+02

1.004-03

5.317+02

5.115-04

4.460-01

2.614+02

1.344-03

4.622+02

6.977-04

5.610-01

2.379+02

1.806-03

3.996+02

9.663-04

7.100-01

2.121+02

2.472-03

3.423+02

1.371-03

8.800-01

1.885+02

3.324-03

2.923+02

1.910-03

1.770+00

1.211+02

9.383-03

1.830+02

5.911-03

2.240+00

1.029+02

1.361-02

1.555+02

8.709-03

2.810+00

8.751+01

1.946-02

1.312+02

1.272-02

3.540+00

7.384+01

2.877-02

1.092+02

1.885-02

4.460+00

6.236+01

4.238-02

9.065+01

2.814-02

5.610+00

5.256+01

6.256-02

7.488+01

4.217-02

7.100+00

4.398+01

9.370-02

6.119+01

6.431-02

8.800+00

3.731+01

1.358-01

5.045+01

9.506-02

1.120+01

3.105+01

2.067-01

4.176+01

1.476-01

1.410+01

2.600+01

3.093-01

3.476+01

2.241-01

1.770+01

2.178+01

4.621-01

2.894+01

3.382-01

2.240+01

1.811+01

6.990-01

2.392+01

5.177-01

2.810+01

1.516+01

1.045+00

1.991+01

7.802-01

3.540+01

1.264+01

1.575+00

1.652+01

1.185+00

4.460+01

1.055+01

2.375+00

1.373+01

1.799+00

5.610+01

8.813+00

3.572+00

1.144+01

2.721+00

7.100+01

7.375+00

5.427+00

9.516+00

4.156+00

8.800+01

6.278+00

7.936+00

8.073+00

6.104+00

1.120+02

5.264+00

1.213+01

6.746+00

9.372+00

1.410+02

4.476+00

1.813+01

5.718+00

1.406+01

1.770+02

3.842+00

2.684+01

4.894+00

2.089+01

2.240+02

3.310+00

4.008+01

4.205+00

3.129+01

2.810+02

2.898+00

5.854+01

3.627+00

4.585+01

3.540+02

2.562+00

8.541+01

2.239+00

6.708+01

4.460+02

2.297+00

1.234+02

2.897+00

9.719+01

5.610+02

2.093+00

1.760+02

2.633+00

1.389+02

7.100+02

1.933+00

2.502+02

2.427+00

1.980+02

8.800+02

1.826+00

3.408+02

2.287+00

2.702+02

     
     
 2. Электроны

     
Таблица 2

     

Энергия, МэВ

Алюминий

Тканеэквивалентное вещество

,
МэВ·г·см

,
г·см

,
МэВ·г·см

,
г·см

0,010

16,49

3,539·10

21,89

2,592·10

0,015

12,20

7,111·10

15,98

5,304·10

0,020

9,844

1,170·10

12,79

8,827·10

0,030

7,287

2,367·10

9,368

1,810·10

0,040

5,909

3,900·10

7,547

3,008·10

0,050

5,039

5,738·10

6,408

4,451·10

0,060

4,439

7,855·10

5,626

6,121·10

0,080

3,661

1,284·10

4,617

1,007·10

0,10

3,117

1,872·10

3,994

1,474·10

0,15

2,513

3,659·10

3,142

2,903·10

0,20

2,174

5,804·10

2,711

4,624·10

0,30

1,839

1,083·10

2,285

8,677·10

0,40

1,680

1,652·10

2,085

1,327·10

0,50

1,592

2,260·10

1,972

1,820·10

0,60

1,540

2,894·10

1,901

2,335·10

0,80

1,486

4,206·10

1,825

3,408·10

1,0

1,465

5,546·10

1,789

4,508·10

1,5

1,460

8,912·10

1,764

7,306·10

2,0

1,475

1,224·10

1,767

1,011·10

3,0

1,510

1,869·10

1,792

1,563·10

4,0

1,540

2,491·10

1,818

2,102·10

     
     Примечание. При пользовании табл.1 и 2 для промежуточных значений ,  и  применяют линейную интерполяцию между двумя соседними значениями.
     
     

ПРИЛОЖЕНИЕ 2
Справочное

     

     1. Пример расчета поглощенной и эквивалентной дозы от протонов ЕРПЗ  

     
     1.1. Выбираем точку околоземного космического пространства с координатами 0,2 Гс и 2,0. Для эпохи максимума солнечной активности по ГОСТ 25645.138 находим коэффициенты для расчета энергетического спектра плотности потока протонов с энергией больше :
     

,                                                  (20)

     
     5,9910·10, -8,6988·10, 5,2703·10, -1,5394·10, 2,1688, -1,1921.
     
     Дифференциальный энергетический спектр  вычисляем из выражения (20).
     
     1.2. Из ГОСТ 25645.138 для выбранных ,  находим 0,1 МэВ, 200 МэВ.
     
     1.3. Задаем заряд вещества защиты 13 и массовое число 27.
     
     1.4. Задаем толщину защиты 0,01; 0,05; 0,1; 0,5; 1 г/см.
     
     1.5. Вычисляем энергетический спектр плотности потока протонов в прот/(см·с). Для значений энергий из п.2.2.4 получаем:
     
     

Таблица 3

     

, МэВ

0,1122

2,73·10

0,1414

2,40·10

0,1778

2,28·10

0,2239

2,29·10

0,2818

2,35·10

0,3548

2,44·10

0,4467

2,50·10

0,5623

2,51·10

0,7079

2,44·10

0,8913

2,30·10

1,122

2,07·10

1,414

1,78·10

1,778

1,48·10

2,239

1,17·10

2,818

9,01·10

3,548

6,69·10

4,467

4,88·10

5,623

3,43·10

7,079

2,37·10

8,913

1,67·10

11,22

1,12·10

14,14

7,65·10

17,78

5,22·10

22,39

3,51·10

28,18

2,39·10

35,48

1,60·10

44,67

1,05·10

56,23

6,80·10

70,79

4,17·10

89,13

2,54·10

112,2

1,35·10

141,4

6,63

177,8

2,89

112,9

-

     
     
     1.6. Вычисляя по полученным результатам дифференциальный энергетический спектр протонов и подставляя полученные значения в формулы (1-4), получаем мощности доз , рад/с, и , бэр/с.
     
     

Таблица 4

     

, г/см

0,01

0,05

0,1

0,5

1,0

4,22·10

7,64·10

3,04·10

3,09·10

1,14·10

7,64·10

8,26·10

2,50·10

1,32·10

3,75·10

     
     
      2. Пример расчета поглощенной (эквивалентной) дозы от электронов ЕРПЗ
     
     1. Выбираем точку околоземного космического пространства с координатами 0,2 Гс и 2,0. Для эпохи максимума солнечной активности по ГОСТ 25645.139 находим коэффициенты для расчета энергетического спектра плотности потока электронов с энергией больше
     

,                                                (21)

     
где  - энергия электрона, кэВ;

     -8.2008+00; 2.1010+01; -9.5091+00; 1.6550+00; -1.3243-01.
     
     Дифференциальный энергетический спектр  вычисляем из выражения (21).
     
     2. Из ГОСТ 25645.139 для выбранных  и  находим 0,04 МэВ, 4,0 МэВ.
     
     3. Задаем заряд вещества защиты 13 и массовое число 27 г/моль.
     
     4. Задаем толщины защиты 0,005; 0,01; 0,02; 0,05; 0,1; 0,2; 0,5; 1,0; 1,5.
     
     5. Задаем полное число статистических испытаний 1000.
     
     6. Вычисляя по изложенному выше алгоритму, получаем значения мощности доз в ткани  за защитой из алюминия.
     
     

Таблица 5

     

, г/см

0,005

0,01

0,02

0,05

0,1

, сГр/с

3,46·10

2,65·10

1,40·10

1,83·10

2,26·10

, г/см

0,2

0,5

1,0

1,5

-

, сГр/с

1,99·10

5,34·10

6,11·10

6,45·10

-