ГОСТ 25645.115-84 Атмосфера Земли верхняя. Модель плотности для баллистического обеспечения полетов искусственных спутников Земли

95 коп.

государственный стандарт

СОЮЗА ССР

АТМОСФЕРА ЗЕМЛИ ВЕРХНЯЯ

МОДЕЛЬ плотности для

БАЛЛИСТИЧЕСКОГО ОБЕСПЕЧЕНИЯ ПОЛЕТОВ ИСКУССТВЕННЫХ СПУТНИКОВ ЗЕМЛИ

ГОСТ 25645.115—84

Издание официальное

ГОСУДАРСТВЕННЫЙ КОМИТЕТ СССР ПО УПРАВЛЕНИЮ КАЧГСГВОМ ПРОДУКЦИИ И СТАНДАРТАМ

MotKBd

УДК 551 51*001.572:006.354 Группа Т27

государственный стандарт союза сер

Атмосфера земли верхняя

МОДЕЛЬ ПЛОТНОСТИ ДЛЯ БАЛЛИСТИЧЕСКОГО ОБЕСПЕЧЕНИЯ ПОЛЕТОВ ИСКУССТВЕННЫХ СПУТНИКОВ ЗЕМЛИ

Earth upper atmosphere.

Density model for ballistics support of flights of artificial Earth satellites

ГОСТ

25645.115—84

ОКСТУ 0080

Дата введения 01.07*85

Настоящий стандарт устанавливает модель плотности, методику расчета и значения параметров плотности атмосферы Земли в диапазоне высот от 120 до 1500 км для различных уровней солнечной активности при известных дате и времени полета искусственного спутника Земли (ИСЗ).

Стандарт предназначен для баллистического обеспечения полетов искусственных спутников Земли и приведения результатов расчетов к одинаковым условиям.

Требования настоящего стандарта являются рекомендуемыми.

(Измененная редакция, Изм. № 1).

1. ОСНОВНЫЕ ПОЛОЖЕНИЯ И ФОРМУЛЫ ДЛЯ РАСЧЕТА

ПЛОТНОСТИ АТМОСФЕРЫ

1.1. Модель плотности верхней атмосферы Земли построена по данным о торможении ИСЗ за период с 1964 по 1987 г.

1.2. Модель плотности атмосферы представлена в аналитическом виде. Коэффициенты модели плотности атмосферы приведены к фиксированным значениям индексов солнечной активности (далее — фиксированным уровням солнечной активности).

1.3. Изменения плотности атмосферы относительно среднего значения в фазе цикла солнечной активности приведены в справочном приложении 1.

1.4. Методика расчета коэффициента аэродинамического сопротивления ИСЗ приведена в справочном приложении 2.

Издание официальное

★ © Издательство стандартов, 1984

© Издательство стандартов, 1991 Переиздание с Изменениями Настоящий стандарт не может быть полностью или частично воспроизведен, тиражирован и распространен без разрешения Госстандарта СССР

1*5. Рекомендации по использованию модели плотности верхней атмосферы для баллистического обеспечения полетов ИСЗ> представлены в справочном приложении 3.

1.6. Описание и текст программы расчета плотности атмосферы. Земли на языке Фортран представлены в справочном приложении 4.

1.7. Угловые скорости вращения атмосферы и Земли приняты равными. Плотность атмосферы для высоты 120 км определяют по ГОСТ 4401—81.

1.8. Плотность атмосферы р, кг/м³, вычисляют по формуле

р=р_кад₁*./с₈*4, (1)

где p„=<V^exPl^ai—^1a2(^h—<?₃)^iy2 L где a_Q= 9,80665 кг/м³;

K₀=l + (to+li-h+l₂-h*)(F₈₁-F₀)-K_F\ K₁=\+(c₀-\-c₁-h-\-c₂-h*-\’C₃-h³)-co$ⁿ*+ⁿ'-^h(q;⁽2),

где coscp= -j- (2>sm$Q+cosSQ(A:'COsp+y*^siⁿP))>

|3=ctQ—5* — ⁽*⁾з(^ — 10800)—{-cp!;

“ 1 ‘h-{~d₂'h²) * A{D)\

К2= ^^Jr{b₀^J\~b₁'h^Jrb₂ ■ h²) -(F—F₈₁)[F;

/C₄ — 1 + (e₀ + ~\~^ез^e^³) ■ ^,j^p+^e6 * Kp ) >

t — московское декретное время, с;

h — геометрическая высота над поверхностью общего земного» эллипсоида (разницей между высотами, отсчитываемыми от уровня моря и от общего земного эллипсоида, можно пренебречь), км; а©, б© — прямое восхождение и склонение Солнца, рад;

(о₃ =7,292115* 10^—5 рад/с — угловая скорость вращения Земли;. х, tjj z — гринвичские координаты точки пространства, км;

S* — звездное время в среднюю гринвичскую полночь, рад;

г— x²-{-y² + z² — расстояние от центра Земли до точки про

странства с координатами х_} у, z_f км;

р_н— плотность ночной атмосферы, кг/м³;

F₀ — фиксированный уровень солнечной активности, 10“²² Вт/ _ (м²-Гц);

F — среднесуточное значение индекса F\0,7, 10~²² Вт/(м²*Гц);

^8i — среднее взвешенное значение индекса РюлЛО”²² Вт/(м²-Гц), осредненное за 3 оборота Солнца (81 сут) и рассчитываемое по> формуле

о

² F_t W_t

р _ i=—80_

^ГЫ— о

² w_b

i=— 80

где Г, = 1+ Ц£- , i=—80, -79, ,0,

значение индекса_г~—80 соответствует дате, отстоящей на 80 сут

от даты расчета -Pst» а i = 0 соответствует дате расчета Psi,

Р ю,7 — индекс солнечной активности, равный плотности потока радиоизлучения Солнца на волне длиной 10,7 см (на частоте 2800 МГц), выраженной в солнечных единицах потока 10~²² Вт/(м²-Гц);

К_р — квазилогарифмический планетарный среднесуточный индекс геомагнитной активности,

Kf — ЫО²² м²*Гц/Вт — нормирующий множитель,

Ко — множитель, учитывающий изменение плотности атмосферы, связанное с отклонением Fqi от F₀;

Ki — множитель, учитывающий суточный эффект в распределении плотности,

Кг ■— множитель, учитывающий полугодовой эффект,

Ко, — множитель, учитывающий изменение плотности, связанное с отклонением F от К₈₁

Ка — множитель, учитывающий зависимость между плотностью атмосферы и геомагнитной активностью,

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

ф! ■— коэффициент модели, равный углу запаздывания максимума плотности по отношению к максимуму освещенности, рад,

A(D) — множитель, характеризующий влияние полугодового эффекта на плотность атмосферы,

D — количество суток, прошедших от начала года,

ct_2j do, d_u d₂, bo, b\, b₂, Cq, C\, c_2> Сз, tio, fi\, ерь £<ъ ^e\> ^e2, £з> <64, e*» lo, h, h — коэффициенты модели, используемые для расчета плотности атмосферы при различных значениях F₀ Значения коэффициентов для 120 кмсй<180 км приведены в табл. 1, для 180с/г-<600 км — в табл 2, для бООсйс 1500 км— в табл 3.

Таблицы составлены для семи фиксированных уровней солнечной активности F₀ = 75, 100, 125, 150, 175, 200, 250 (10~²² Вт/

<(м²*Гц)) При расчете плотности атмосферы используют значения коэффициентов модели для F₀, ближайшего к значению К₈ь

. 4 ГОСТ 25645.115—84

ГОСТ 25645115—84 С. 5

Примечание к табл 1—3 Число со знаком минус, стоящее после значения параметра, является показателем степени десяти — сомножителя зна/чения параметра

С в ГОСТ 25645 115—84

Значения множителя A(D)

Таблица 4

(Измененная редакция, Изм. № 1)

Значения множителя A(D) с шагом 10 сут приведены в табл. 4 (на промежуточные моменты времени значения множителя A(D) находят линейной интерполяцией).

2 ТАБЛИЦЫ ПАРАМЕТРОВ ПЛОТНОСТИ АТМОСФЕРЫ

Параметры плотности атмосферы для семи фиксированных уровней солнечной активности в функции высоты приведены в табл. 5—11:

р_н — плотность ночной атмосферы, вычисленная по формуле

Ко — величина, характеризующая вековое изменение плотности ночной атмосферы в 11-летнем цикле солнечной активности

* h + /₂ -h²;

§

К i — амплитуда суточного эффекта, вычисленная по формуле

К[ — c_Q + Cj • /i-bc_a • /i²-f с₃ • h³-,

f / f

К2, Кз, К4 —величины, характеризующие, соответственно, влияние полугодового эффекта, радиоизлучения Солнца и геомагнитной активности, вычисленные по формулам:

^С₂ ■ h-\-d^ • h²\

-h^l\

K'₄=e₀+e^h-\-e₂^h² + e^h³.

Для fW'S-lO-²² Вт/(м².Гц)

8 ГОСТ 25645.115—84

ГОСТ 25645.115—84 С. 9

10 ГОСТ 25645.115—84

Пара*

атмосферы

ГОСТ 25645.115—84 С 11

. 12 ГОСТ 25645.115—84

ГОСТ 25645 115—84 С 13

Для Fo—250'Ю^²² Вт/(мЧц)

14 ГОСТ 25645.115—84

Примечание к табл. 5—11, Число со знаком минус, степени десяти — сомножителя значения параметра.

Разд. 2. (Измененная редакция, Изм. № 1).

Продолжение табл. 1]

стоящее после значения параметра, является показателем

ГОСТ 25645.118—84 С. 15

ПРИЛОЖЕНИЕ 1

Справочное

ИЗМЕНЕНИЯ ПЛОТНОСТИ АТМОСФЕРЫ В ФАЗЕ ЦИКЛА СОЛНЕЧНОЙ АКТИВНОСТИ

I Значения предельных относительных отклонении плотности атмосферы m 1 п И &Р max по высоте от ее средних значений за период в годы низкого, среднего и высокого уровней солнечной активности приведены в табл 1.

В соответствии с ГОСТ 25645 102—83 считают низкий уровень солнечной активности — минимум, средний — фаза роста и спада, а высокий — максимум солнечной активности

Таблица I

Предельные относительны^ отклсмсния плотности

Значения 6p_min и бртах получены с использованием формулы (1) настоящего стандарта и данных об индексах солнечной активности F и К р за период 1964—1987 гг.

2. Значения предельных относительных отклонений изменения плотности атмосферы за одни сутки бр ₍ по высоте за период в годы низкого, среднего и высокого уровня солнечной активности приведены в табл. 2.

3. Изменения отношения плотности атмосферы при различных фиксированных уровнях солнечной активности к плотности при уровне F₀=75*lO~²² Вт/ (м²-Гц) по высоте указаны в табл. 3

Нормированные автокорреляционные функции невязок модели плотности атмосферы для низкого, среднего и высокого уровней солнечной активности (СА) приведены на черт. 1.

Средние квадратические отклонения невязок модели плотности атмосферы, вычисленные по торможению ИСЗ в функции высоты (до 620 км), приведены на черт. 2.

Данные табл. 1—3 используют при определении диапазона изменений баллистического коэффициента, как коэффициента согласования, и как исходные данные при априорной оценке точности прогнозирования движения ИСЗ.

Таблица 2

Таблица 3

1—высокий уровень; 2—средний уровень; 3—низкий уровень.

Черт. 1

Черт. 2

ПРИЛОЖЕНИЕ 2

Справочное

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

При определении коэффициента аэродинамического сопротивления (С_х ^

ИСЗ приняты следующие основные исходные предположения режим обтекания поверхности ИСЗ разреженным газом верхней атмосферы — свободномолекулярный, набегающий поток газа, имеющий переносную скорость и — «+[(оХ^гЬ состоит из молекул с температурой Т, скорости которых распределены по закону Максвелла; молекулы набегающего потока газа испытывают однократные соударения с поверхностью ИСЗ, характер отражения молекул от поверхности ИСЗ— полностью диффузный (коэффициенты аккомодации нормального и касательного импульсов равны единице), при этом после отражения молекулы имеют скорости, распределенные по закону Максвелла, с температурой поверхности T_w; затенение одних элементов ИСЗ другими учитывают по аналогии с геометрической оптикой

При принятых исходных предположениях проекции на орты пит вектора аэродинамической силы Р_п и Р^ , отнесенные к скоростному напору pv²/2go и

действующие на элементарную площадку dF соответственно равны*

P„-eos©

%(z)

KW2Moo

+

1-ferfz У Ml

%(z) V-All '

P_r =sin0

X(f)

M,

Х(г) = ехр(—г²) + |^я z(l -f erf 2), erfz= JL- ( exp(—x²)dx, (1)

it 0

z= |^/"^/2/Vf_ooco-0,

где M ^—v/a^ — число Маха набегающего потока,

а— скорость звука в невозмущенном потоке, м/с,

у — отношение удельных теплоемкостей, п — внутренняя нормаль к элементу dF k поверхности,

© — угол между вектором скорости набегающего потока и внутренней нормалью, рад, erf г — интеграл вероятности

Для определения коэффициента аэродинамического сопротивления С _х (а)

а

ИСЗ, движущегося под углом атаки а, наружную поверхность его условно разбивают на К элементов простой геометрической формы (например, пластина, цилиндрическая, коническая, сферическая поверхности и др ) Для каждого ^-го-элемента (£=1, 2, . , К) вычисляют силу аэродинамического сопротивления путем интегрирования по обтекаемой поверхности F & выражения проекции на ось х_а силы, действующей на элементарную площадку dF ь . Результаты вычислений, полученных для К элементов, суммируют и относят к характерной площади ИСЗ Л* —const и аэродинамическому скоростному напору рv²/2g_Q

к

Сх_а (a)=2R_Xa /F_Mpu²= 2 I (Р_лсо,0+Р. sine)dF_k/F_v, (2)

*=i ^Fk

где R _x — проекция суммарной аэродинамической силы на ось х_а (cos (х_а, а

*0 = 1).

При определении в момент времени t коэффициента аэродинамического сопротивления С_х (ot (t)) неориентированных ИСЗ (вращающихся относительно не-

а -*■

которой оси спутника) в силу малости угловой скорости w вторую составляющую скорости v не учитывают. При этом для определения С х ^ в интервале

времени от t до t+T,_ip усредняют по времени значения коэффициентов сопро тивления:

Т

С_Ха= -4— С C_Xa(a(t))dt, (3)

¹ДР о

где Т д_р — период обращения ИСЗ.

Погрешность при определении коэффициентов аэродинамического сопротивления ИСЗ по приведенной методике в общем случае не более 30 %. Для ИСЗ,

близких по форме к сферической, погрешность в определении С _х оценивают

а

меньше 7 %.

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

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

ПОЛЕТОВ ИСЗ

При баллистическом обеспечении полетов ИСЗ с целью повышения точности определения и прогнозирования орбит ИСЗ рекомендуют рассчитывать плотность

с учетом текущих и прогнозируемых индексов F и Ар , а также суточного и полугодового эффектов.

Обеспечение данными об индексах солнечной и геомагнитной активностей F и Кр возложено в настоящее время на Службу Солнца Кроме того, данные об этих индексах публикуют в специальной литературе При использовании этих данных следует учесть запаздывание изменения плотности атмосферы по отношению к временам соответствующих солнечной и геомагнитной активностей Для солнечной активности это запаздывание — 1,7 сут, а для геомагнитной — 0,6 сут.

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

В последнем случае в качестве F принимают среднее значение F_{10 7} за предшест

вующие 30—50 сут, в качестве Fai — среднее значение Г_]07за предшествующие

81 сут, а в качестве — значение, ближайшее к F_8U для которого есть коэффициенты модели. При отсутствии данных об индексе Кр его значение принимают равным 2,66667.

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

в выражение Ка в этом случае вместо индекса Кр подставляют индекс К_р/- , вычисляемый по формуле

где Кpj — 3-часовой планетарный индекс геомагнитной активности,

/=0, 1,.-. п;

д_ ^{при /=?ь0;} _г== (°>^{3 при Л>0}

1о при / = 0 [0,7 при Д<0.

В этом случае запаздывание для геомагнитной активности — 0,25 сут.

Таблица 1

Коэффициенты множителя Ка для 3-часовых индексов

Значения индекса Кр, дят по табл 2.

используемого в модели плотности атмосферы, нахо-

Таблица 2

Таблица перевода

Индекс А р — по ГОСТ 25645 302—83

Знаки « + » и «—» справа от значения Кр означают, что к значению, заданному цифрой, нужно прибавить, соответственно, +0,3333 или —0 3333 т е например, 3⁺ = 3 + 0,3333 = 3,3333, а 4~ = 4+(—0,3333) =3,6667

Для нахождения промежуточных значений используют линейную интерполяцию.

При отсутствии данных об индексе К _р допустимо использование индекса^ осредненного по нескольким станциям, находящимся на территории СССР

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

(Измененная редакция, Изм. № 1)

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

ОПИСАНИЕ И ТЕКСТ ПРОГРАММЫ РАСЧЕТА ПЛОТНОСТИ АТМОСФЕРЫ

ЗЕМЛИ НА ЯЗЫКЕ ФОРТРАН

\ Для высот км плотность р рассчитывают по формуле

р—A_Lexp(—Ku(h~h_l)-\'K2i(h—hi)¹) >

где A_t Ки> Kzi — коэффициенты, приведенные в таблице, i ■— номер слоя

2 Для высот й>120 км плотность атмосферы р рассчитывают по формуле 1 настоящего стандарта

3 Обращение к программе

CALL DENSIT (ХК, ТН, SZ, НКМ, SUN, АР, FI, FSH, D, PR, RO, RU)

4 Описание параметров Входные параметры

ХК — массив, содержащий координаты точки в гринвичской системе координат, длина массива равна 3,

ТН — московское декретное время, с,

SZ — звездное время в гринвичскую полночь, рад,

НКМ — высота ИСЗ над поверхностью общего земного эллипсоида, км SUN — массив, содержащий прямое восхождение (a q ) и склонение (6q )

Солнца, рад, длина массива равна 2, _

АР — значение индекса А_р или Кр ,

F1 —значение индекса F,

FSH — значение индекса F_8b D — время московское в сутках от начала года,

PR — если PR>0, то при обращении к программе параметру АР присваивают значение Кр , если PRcO, то параметру АР присваивают значение А_р и оно

переводится в Кр внутри программы (подпрограмма INTERP (АР, КР)) Выходные параметры RO ■— плотность атмосферы, кг/м⁵;

RU .— плотность атмосферы, кгс*с^£/м⁴ 5 Используемые блоки COMMON Программа DENSIT использует три общих блока COMMON /QKOEF/Q (27,3)/YEAR/P(38)

COMMON/UKOEF/U(8,4), GR, WZ, в которых должны содержаться все необходимые для вычисления константы для фиксированного уровня солнечной активности F₀

6 Описания параметров блоков COMMON Q — двумерный массив размерности 27X3, содержащий коэффициенты модели плотности атмосферы для одного из фиксированных уровней F₀ Первые элементы каждого столбца равны F₀, далее следуют коэффициенты модели настоящего стандарта в соответствии с табл 1, 2 и 3 п 18, причем в первом столбце находятся коэффициенты для диапазона высот от 120 до 180 км, во втором

— от 180 до 600 км, а в третьем — от 600 до 1500 км

Р — массив, содержащий значения поправок полугодового эффекта из табл. 4 п I 8 настоящего стандарта, длина массива равна 38,

П — двумерный массив размерности 8X4, содержащий коэффициенты таблицы, WZ — 0 7292115-10-⁴ (угловая скорость вращения Земли, рад/с), GR^=

— 9 80665 кг/м³

5 6 (Измененная редакция, Изм. № 1).

7 Для присвоения нужных значений переменным из упомянутых общих блоков следует воспользоваться подпрограммой VKM(FSH)

CALL VKM(FSH)

Подпрограмма VKM(FSH) находит ближайшее к FSH значение F₀ и засылает соответствующие этому Fo значения констант в COMMON блоки программы DENSIT

ТЕКСТ ПРОГРАММЫ

SUBROUTINE DENSIT (ХК, ТН, SZ, НК.М, SUN,

* АР, FI, FSH, Г, PR, RO, RU)

IMPLICIT REAL*8(A—Н, О—Z)

REAL*® Kl, К2, КЗ, К4, КО, КР DIMENSION ХК(3), SUN (2)

COMMON/UKOEF/U(8, 4), GR, WZ COMMON/QKOEF/Q (27,3) /YEAR/P (88)

H = HKM

RAD = DSQRT(XK (1)*XK(1) +XK(2) *XK(2) +XK(3) *XK(3))

X=XK(1)/RAD

Y=XK(2)/RAD

Z=XK(3)/RAD

IF (H LT.120.DQ) GOTO 5

J=1

IF(H.GE.l,8OD0) J = 2'

IF (H GE 600 DO) J = 3 HD = H*H

RO = DEXP(Q(2, J)—Q(3, J)*DSQRT(H—Q(4, J))) RO=RO*GR

Kl =Q(8, J)+Q(9, J)*H+Q(10, J)*HD+Q(M, J)*HD*H SD = DSIN(SUN(2))

CD = DCOS(SUN (2))

G = S UN (1) +Q (14, J)—SZ—WZ* (TH—10600.DO)

SG = DSFN(G)

CG=DCOS (G)

CF=Z*SD+|CD*(X*CG+Y*SG)

CP=((1D0+CF)«.6D0)**((Q(I2, J)+Q(13, J)*H)*.5D0)

Kl'= 1.D0+K1*CP RO = Kl*RO I = INT(D/10iD0) + l!

SI = P(1) + (P(I+1)—P(I))*(D—(I—1)* 10 DO)/10.D0 K2=Q(16, J)+Q(16, J)*H+Q(17, J)*HD K2=1.D0+K2*SI

RO = K2*RO

K3 = Q(5, J)+Q(6, J)*H+Q(7, J)*HD A= (FI—FSH)/FSH K3 = 1.D0+K3*A RO = K3*RO A=FSH—Q( 1, J)

KO=Q(25, J)+Q(26, J)*H+Q(27, J)+HD K0= l.DO+KOi*A RO = KO*RO

IF (PR GT G.DG) GOTO 3 CALL INTERP(AP, KP)

GOTO 4

3 КР = АР

4 DKP = KP*KP

A = Q(22, J)+Q(2'3, J)*KP+Q(24, J)*DKP

K4 = Q(18, J)+Q(10, J)*H+Q(20, J)*HD+Q(21, J)*HD*H

K4=1.D0+K4*A

RO = K4*RO

GOTO 8

5 CONTINUE DO 6 1=1, 9

IF (H.GE.U(K)—I, 1)) GOTO 7

6 CONTINUE PRINT 120

20 FORMAT (5X, 'H<0')

RETURN

7 1=10—1

DH = H—U (I, 1)

DHD = DH*DH

RO=U(I, 2)*DEXP(—U(I, 3)*DH-f-U(I, 4)*DHD)

8 RU = RO/GR RETURN END

SUBROUTINE VKM(FSH)

IMPLICIT REAL*8(A—H, O—Z)

DIMENSION XF(7)

DIMENSION Rl,(27,7), R2(27,7), R3{27,7)

DIMENSION Rll (54), R12(81), R13(54)

DIMENSION R211(541), R2'2(81), R23i(64)

DIMENSION R31>(54), R32(8fl), R33(54)

EQUIVALENCE (R1 (1,1), R11/(1))

EQUIVALENCE (Rt(l,3), RL2(1))

EQUIVALENCE (RI(1,6), R13(1))

EQUIVALENCE (R|2l(Ll), R21(l))

EQUIVALENCE (R2(1,3), R22(l))

EQUIVALENCE (R2(L,6), R23(l))

EQUIVALENCE (R3i(l,l), R31/(l))

EQUIVALENCE (R3(l,3), R32f(l))

EQUIVALENCE (R3(l,6), R33(l))

COMMON/UKOEF/U (8,4), GR, WZ COMMON/QKOEF/Q (27,8) /YEAR/P (38)

DATA XF/87.5D0, 1[I2‘5|D0I, 187.5D0, 162.5DQ, 187.5D0, 225.D0, 40C.DO/ DATA R11/75.DO, —18.2991(00, .70(0900, 115.3430D0,

* —.682800, .5576D—2, .95238D—6,

* —4.384D0, .08062900, —.4925D— 3, .1I042D—5,

* L5D0, .00600, .541052D0,

* —5.10L9DOL .0625800, —.16720—3,

* — 7.238DO, .1(20300, —.645D—3, —.12080—5,

* —.1'2D0, .0O6DO, .01БОО,

* — J01188HX), .9983D—4, 0.001,

* lOO.DO, — 13.1909D0, .700, lil4.63e&DO,

* —.7804DO, 71730—2, —.5578D—5,

* —4.384D0, .080629DO, —.4025D—3, .1042D—5,

* 1.5D0, .006DO, .541J052D0,

* —5.1O19D0, .0625600, —.16720—3,

* —6.683D0, . 1101(400, —.583750—3, .H083D—5,

* —.1200, .02600, .0O75D0,

* —.'009900, .8212D—4, .31050—8/

DATA R12/1(26.00, —18.5209DO, .6419D0, 115.9569D0,

* —.822D0, .833D—2, —.12330—4,

* .977600, —J02S7DO, .2027D—3, — .4708D—6,

* 1.500, .006D0, .541O52D0I,

* — 5.1019D0, ,06t258D0, —.115720—3,

* —5.36200, .8615D—1, — .443760—3, .8126D—6,

* —.IDO, .0208800, .006300,

* — .0076800, .63620—4, .3H26D—8,

* 15O.D0, — 18.652200, .612400, 116.4154D0,

* — 7376D0, .7597D—2, —.1,209D—4,

* .9776D0, —.'0257DO, .20270—3, —.47060—6,

* 1.5D0, .00600, .54)1051200,

* —5.101900, .0625800, —.1672D—3,

* —4.799D0, .07779200, —.40750—3, 77083D—6,

* —.100, .027500, .0038DO,

* —.O066DO, .46670—4, 0 DO,

* 17500, — 18.6586D0,. 6038DO, 116 353100,

* —.31500, .0232600, 25D—5,

* — .563200, .5743D—2, —9250—5, .4167D—8,

* 1.5D0, .00600, .541052D0,

* —5.101900, .0625800, — .1672D—3,

%—•4.903D0, .08270800, —.46875D—3, .9167D—6,

* —.1200, .0411600, .0014D0¹,

* —.4963D—2, .41360—4, 0 DO/

DATA R 13/200.00, — 18.64(9500, .597400, 116.2144D0,

* —j516100, .005341 DO, —.867H9D—5,

* —,.563200, .5743D—2, —.9)250—53, .4167D—8,

* 1.500, .OQ6DO, .541.05200,

* —5.1019D0, /06258D0, —.1672D—3,

* —5.11500. .08607500, — .45675D—3, .875D—6,

* — .1100, .0381001, 001200,

* —.0041,100, .3463D—4, —.3125(0—8,

* 250.00, —18.707400, .5771200, 116.3395D0,

* — .26(3100, .1929D—2, —.14950—5,

* .4842D4, —.01,60400, .1405D—3, —.3375D-6,

* 1.500, /00600, .54105(200,

* —5 Ю19DO, 06258DO, — 161720—3,

^5,1 —3 13700, 047742D0, — 2275D—3, 39583D—б,

* — O0DO, 0314800, 001 DO,

* — ООЗОЗDO, 2532D—4, — 5556D—9/

DATA R01/75DO, —15 56Q5D0, 8248DO, 76 9132 00,

* — 8607D0, 00786100, — 571 ID—5,

* 1 279100, — 01576D0, 64990—4, — 5145D—7,

* 1 5D0, '006D0, 54105200,

* — 1721 DO, .0KK756D0, — 3605D—5

* — 2152D0, 4167D—2, 1587D— 5, — 1651D—8,

* — 12D0, 006DO 01500,

* — 01698D0, 1448D—3, — 9535D—7,

* 100 DO,—15 640800, 7754D0, 67 916200,

* — 754D0 00686D0, — 460—5,

* 1 2791D0, — 01576D0, 6499D—4, — 5145D—7,

* 1 500, 00600, 541062D0,

* — 1721D0, 005756D0, — 3636D—5,

* — 2162D0, 004086DO, 127D—5, — 15870—8,

* — 12D0, 025D0, 0IO75DO,

* — 01249D—2, 111 ID—3, — 7706D—7/

DATA К22/125 DO —16 222900, 7569D0, 55 816500,

* — 57D0, .0052500, — 3D—5,

* 1 29O3D0, — 0166700, 5964D—4, — 4503D—7,

* 1500, O06D0, 541052DO,

* — 1721D0, 0O5756D0, — 3635D—5,

* — 1486D0, 003263DO, 3143D—5, — 3429D—8,

* — IDO, 0208300, 006300,

* — 78790 —2, 7258D—4, — 3658D—7,

* 150 DO, —16 9752D 0, 673600, 85 444000,

* — 476D0, 004400, — 24D—5,

* 1 2903D0 —0154700, 69640—4, —45030—7,

* 1 5D0, 00600, 54105200,

* — 1721D0, 00575600,-36350—5,

* — 1495D0, — 3182D—2, 2826D—5, — 3365D—8,

* — IDO, 027500 003600,

* — 4882D—2, 4692D—4, — 1742D—7,

* 175 DO, —17 3045D0, 6782D0, 81 9596D0,

* — 292D0, 002800, — 8D—6,

* 205700, — 2912D—2, .17390—4, — 85649D—8,

* li 500, 00600, 54106200,

* — 172IDO, 005756DO, — 3635D—5,

* —081900, 2368D—2, 4698D—5, — 52I06D—8,

* — 12D0, O4116D0, 001400,

— 5195D—2, 4664D—4, —21640—7/

DATA R23/200 DO, — 18 266D0, 5797DO, HOC 941700,

* — 3113D0, 002839DO, — 1089D—5,

* .2067D0, —.00291 IDO, .K739D—4, — .86649D—8,

* 1.5D0, .006DO, .541062D0,

* —.172 IDO, .00575600, —.3635D—5,

* — .082186ID0, .002278(00, .43810—5, — .5143D—8,

* —.11D0, .038IDO, 0012D0,

* —.5017D—2, .42820—4, —.2132D—7,

* 250.DO, —19.278200, .5M&D0, 116.579I2D0,

* —.3307D0, .00287800, —.13780—5,

* ,O0‘15DO, —.2399D—3., .7O06D—5, — .5909D—9,

* 1 500, .006DO, .54106200,

* —.1,72100, .C05756D0, —.368&D—5,

* —.2048D0, .OOB696DO, —.L587D—5, —.3175D—9,

* — .09D0, .034 18D0, .001 DO,

* —.5455D—2, .42730—4, — .2273D—7/

DATA R31/75.DO, —33.2283D0, .1784D0, 555.06(3600,

* .7833D0, .002861 DO, — .19440—5,

* —4.4D0, .03024DO, — .3288D—4, .1012D—7,

* 1.5DO, .00600, .54105200,

* 1.0204DO, .0O2499D0, —.16190—5,

* —3.8DO, .01972D0, —.16830—4, .4938D—8,

* — .1200, .00600, .016D0,

* .1083D0, .66940—4, — 4277D—7,

* 100.DO,—32.773100, .1899D0, 584.24500,

* .725D0, .002675DO, —. 1750—S,

* —4.400, 0302400^ —.3283D—4, .10120—7,

* 1 5D0, .0O6DO, .541O52D0,

* 1.0204DO, .00249900, —.1510D—5,

* —3.7D0, .01783DO, —.15060—4, .358D—8,

* — .12D0, .025D0, .0076DO,

* .8317D—2, .4837D—4, —.30390—7/

DATA R32/125.DO, —31.6716D0, .2265DO, 571.5408DO,

* .6IDO, .234(30—2, —.14330—5,

* —8.98D0, J04087D0, —.395D—4, .1123D—7,

* I.5D0, .00600, .541O52D0,

* l'.02O4D0, .00249900, —.1519D—5,

* —3.7D0, .0175D0, —.15D—4, .37040—8i,

* —.IDO, .02083DO', .0063DO,

* .4667D—2, .4606D—4, —.2722D—7,

* ISO. DO, —29.7592DO, .294800, 528.3389DO,

* —.O933D0, .008038DO, —.171 ID—5,

* —8.98D0, .04087D0, —.3950—4, .11230—7,

* 1.500, .006DO, .541052DO,

* 1.02O4D0, .OOE499DO, — .161,9D—5,

* —4.4D0, . 1981D—1, —.1806D—4, .49380—8,

* —.IDO, .0275DO, .008800,

* —.0133300, .71I&7D—4, —.3618D—7,

* 175.DO, — 28.8463D0, .314D0, 509.72000,

* — .3333DG, .3522D—2, —.13890—5,

* —15.78DO, .05757DO, —.6322D—4^ .15P2D—7,

* 1.5D0, .096D0, .54105200^

* 1.02O4D0I, .002499DO, —. 1510D—5,

* —3.6D0, ,0169400, —.1556iD—4, .40121D—

* —.12D0, .04 Г16DO, .001400,

* —.0035D0, .4317—4, —.2056D—7/

DATA R33/20IO.DO, —26.2994D0, .381700, 434.122D0,

* —.4333D0, .60352200, —. 18890—5,

* — 15.78D0, .05757DO, — 532t2D—4, .1612D—7,

* 1.5DO, .00600, .541062D0,

* 1.0I2O4DO, .002499 00, —. 1619D—6,

* — 3.6DO, .O1653D0, — .15380—4, .4321D—8,

* —.11D0, .036 IDO* .0011200,

* —..GOI15DO, .3250—4, —.1389D—7,

* 250.DO, — 28.6627D0, .423IDO, 336.4318D0,

* — .175D0, .0O2642DO, — .1417D—5,

* —9.75D0, .003383DO¹, — .2694D—4, .648ID—3,

* 1.5DO, .O06DO, .541052DO,

* 1.0204DO, .00249900, —.1519D—5,

* .IDO, .OQ2639DO, —.27780—6, —.6173D—9,

* — .O9D0, .03118DO¹, .001 DO,

* — .8333D—3, .21817D—4, — .M29D—7/

WZ = .72921H5O0D—4

GR = 9.8O655D0 DO 1 J — 1,7

IF (FSH.LE.XF(J)) GOTO 5 1 CONTINUE

J = 6

5 DO 4 N — 1,27 Q(N, 1) — R1 ('N, J)

Q(N, 2) — R2(N, J)

Q(N, 3)=R3(N, J)

4 CONTINUE RETURN END

SUBROUTINE INTERP(AP, KP)

IMPLICIT REAL*8(A—H, O—Z)

DIMENSION PAP (28)

REAL*& KP

DATA PAP/O.DO, 2.DO, 3.D0, 4.DO, 5.D0, 6.D0, 7.D0, 9.D0, 12.D0, 15.D0,

* 18.D0, 22.D0, 27.D0, 32.DO, 39.DO, 4(8.00, 56.D0, 67.D0, 8aD0, 94.D0,

* 111.DO, 132.DO, 154.DO, 179.D0, 2O7.D0_v 236.D0, 300.DO, 40O.DO/

D = .333333333333D0

DKP—O.DO

DO 10 1=1,28

IF (АР LE.PAP (1+1)) GOTO 1 DKP = DKP+D 10 CONTINUE 1 A=PAP(I+1)—PAP (I)

B = AP—PAP (I)

C = B/A

KP = DKP+D*C

RETURN

END

BLOCK DATA

IMPLICIT REAL*8(A—H, O—Z)

COMMON/QKOEF/Q (26,3) /YEAR/P (38)

COMMON/UKOEF/U (8,4), GR, WZ

DATA U/O.DO, 20.DO, 60.DO, 100.D0, 150.DO, 300.DO, 60O.D0, 900.DO,

* 1.238D0, .09013DO, .31043D—3, .53675D—6, .20078D—8, .106610—10,

* 11273D—1>2, .56916D—14, .90764D—1, .16739D0, .1237600, .17527D0,

* .46®25D—1, .19885D—1, .14474D—1, .39247D—2, — .20452D—2,

* .62689D—3, —.68999D—3, .1287D—2, .10167D—3, .97266D—5, .15127D—4,

* O.DO/"

DATA

R/—.02800, —.045D0, —.04700, —.O35D0, —.OllODO, 0.02200, .057D0, .09DO,

* .114D0, , 125D0, .lli&DO, .09600, .06DO, .O13D0, -4137D0, —.O86D0, —.12800,

* —.16200, —.185D0, _ 199D0, —.202DO, —.193D0, —.173D0; — 1400, — —.O96D0,

* — 042D0, 0.015DO, O7D0, .115D0, .144D0, .153D0, 146£>0, I2D0, 084DO,

* .044D0, .O06D0, —,Oi23DO, —.04D0/

END

(Измененная редакция, Из м. Ml).

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

1. РАЗРАБОТЧИКИ

И. К Бажинов, д-р техн. наук; П. А. Буцко, канд. техн. наук; И. И. Волков, канд. техн. наук; В. Н. Жуков; Ю. А. Климов, канд. техн. наук; Е. И. Князева; Н. Н. Князев; В. Я. Пономарев, канд. техн. наук; В. Г. Фарафонов, канд. техн. наук; Н. М. Щербаков, канд. техн. наук; В. В. Ястребов, В. Д. Ястребов, канд. техн. наук

2. УТВЕРЖДЕН И ВВЕДЕН В ДЕЙСТВИЕ Постановлением Государственного комитета СССР по стандартам от 24.08.84 № 2997

3. ВВЕДЕН взамен ГОСТ 22721—77

4. ССЫЛОЧНЫЕ НОРМАТИВНО-ТЕХНИЧЕСКИЕ ДОКУМЕНТЫ

5. Переиздание (май 1991 г.) с Изменением № 1, утвержденным в декабре 1990 г. (ИУС 3—91)

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

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

ГОСТ 4401—81 ГОСТ 25645.102—83 ГОСТ 25645.302—83

п 1 7

Приложение 1 Приложение 3

Редактор В. С. Бабкина Технический редактор В. //. Малькова Корректор О. Я. Чернецова

Сдано в наб. 21.02.91 Подп,. к печ. 07.06.91 2,0 уел. п. л. 2,13 уел. кр.-отт. 2,32 уч.-изд. л.

Тираж 3000 экз. Цена 95 к.

Ордена «Знак Почета» Издательство стандартов, 123557, Москва, ГСП,

Новопресненск1ий пер., 3.

Калужская типография стандартов, ул. Московская, 256. Зак. 424