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. Плотность атмосферы р, кг/м3, вычисляют по формуле
р=ркад1*./с8*4, (1)
где p„= K0=l + (to+li-h+l2-h*)(F81-F0)-KF\ K1=\+(c0-\-c1-h-\-c2-h*-\’C3-h3)-co$n*+n'-h(q;(2), где coscp= -j- (2>sm$Q+cosSQ(A:'COsp+y*sinP))> |3=ctQ—5* — (*)з(^ — 10800)—{-cp!; “ 1 ‘h-{~d2'h2) * A{D)\ К2= ^Jr{b0J\~b1'hJrb2 ■ h2) -(F—F81)[F; /C4 — 1 + (e0 + ~\~езe^3) ■ ,j^p+e6 * Kp ) > t — московское декретное время, с; h — геометрическая высота над поверхностью общего земного» эллипсоида (разницей между высотами, отсчитываемыми от уровня моря и от общего земного эллипсоида, можно пренебречь), км; а©, б© — прямое восхождение и склонение Солнца, рад; (о3 =7,292115* 10—5 рад/с — угловая скорость вращения Земли;. х, tjj z — гринвичские координаты точки пространства, км; S* — звездное время в среднюю гринвичскую полночь, рад; г— x2-{-y2 + z2 — расстояние от центра Земли до точки про странства с координатами х} у, zf км; рн— плотность ночной атмосферы, кг/м3; F0 — фиксированный уровень солнечной активности, 10“22 Вт/ _ (м2-Гц); F — среднесуточное значение индекса F\0,7, 10~22 Вт/(м2*Гц); ^8i — среднее взвешенное значение индекса РюлЛО”22 Вт/(м2-Гц), осредненное за 3 оборота Солнца (81 сут) и рассчитываемое по> формуле о 2 Ft Wt р _ i=—80_ ГЫ— о 2 wb i=— 80 где Г, = 1+ Ц£- , i=—80, -79, ,0, значение индекса_г~—80 соответствует дате, отстоящей на 80 сут от даты расчета -Pst» а i = 0 соответствует дате расчета Psi, Р ю,7 — индекс солнечной активности, равный плотности потока радиоизлучения Солнца на волне длиной 10,7 см (на частоте 2800 МГц), выраженной в солнечных единицах потока 10~22 Вт/(м2-Гц); Кр — квазилогарифмический планетарный среднесуточный индекс геомагнитной активности, Kf — ЫО22 м2*Гц/Вт — нормирующий множитель, Ко — множитель, учитывающий изменение плотности атмосферы, связанное с отклонением Fqi от F0; Ki — множитель, учитывающий суточный эффект в распределении плотности, Кг ■— множитель, учитывающий полугодовой эффект, Ко, — множитель, учитывающий изменение плотности, связанное с отклонением F от К81 Ка — множитель, учитывающий зависимость между плотностью атмосферы и геомагнитной активностью, <р -— центральный угол между точкой пространства, для которой рассчитывается плотность, и точкой пространства с максимальным значением плотности в ее суточном распределении, рад, ф! ■— коэффициент модели, равный углу запаздывания максимума плотности по отношению к максимуму освещенности, рад, A(D) — множитель, характеризующий влияние полугодового эффекта на плотность атмосферы, D — количество суток, прошедших от начала года, ct2j do, du d2, bo, b\, b2, Cq, C\, c2> Сз, tio, fi\, ерь £<ъ e\> e2, £з> <64, e*» lo, h, h — коэффициенты модели, используемые для расчета плотности атмосферы при различных значениях F0 Значения коэффициентов для 120 кмсй<180 км приведены в табл. 1, для 180с/г-<600 км — в табл 2, для бООсйс 1500 км— в табл 3. Таблицы составлены для семи фиксированных уровней солнечной активности F0 = 75, 100, 125, 150, 175, 200, 250 (10~22 Вт/ <(м2*Гц)) При расчете плотности атмосферы используют значения коэффициентов модели для F0, ближайшего к значению К8ь Коэффициенты —22 Фиксированный уровень солнечной активности F0 10 Вт/(м* Гц) Обозначение Размерность 75 100 125 150 175 200 250 ах _ -18 2991 -18 1908 -18 5209 -118 6522 -18 6586 -18 6495 -18 7074 а2 KM~1/2 0 7009 0 7000 0 6419 0 6124 0 6038 0 5974 0 5772 аъ км 115 3430 114 6386 115 9569 1164154 116 3531 1162144 116 3395 d о — —5 1019 -5 1019 -5 1019 -5 1019 -5 1019 -5 1019 —51019 dx км-1 6 2581-2 6 2580-2 6 2580-2 6 2580-2 6 2580—2 6 2580-2 6 2580-2 d2 км~2 -1 6721-4 —1 6722-4 -1 6722-4 -1 6722-4 -1 6722—4 —1 6722-4 -1 6722-4 h — -0,6828 -0,7804 -0 8220 -0 7376 -0 3150 -0 5161 -0 2531 КМ-1 5 5762-3 7 1730-3 8 3300-3 7 5970-3 2 3250-3 5 3410-3 1 9290-3 ^2 км 2 9 5238—7 -5 5780-6 -1 2330—5 -1 2090-5 2 5000-6 -8 6719-6 1 4950—6 C0 —- -4 3840 -4 3840 0 9776 0 9776 -0 5632 —0 5632 0 4842 Cl КМ-1 8 0629—2 8 0629-2 -2 5700-2 -2 5700-2 5 7430-3 5 7430-3 -1 6040-2 C2 КМ-2 —4 9250—4 -4 9250-4 20270-4 2 0270-4 -9 2500—6 —9 2500-6 1 4050—4 c3 КМ-3 1 0420-6 1 0420-6 -4 7080-7 -4 7080-7 4 1670-9 4 1670-9 —3 3750-7 Щ — 15 1 5 I 5 1 5 1 5 15 15 Щ КМ-1 0,006 0 006 0006 0 006 0006 0 006 0 006 рад 0 5411 05515 05585 0 5585 0 5585 0 5585 05585 Co — -7 2380 —6 6830 -5 3520 -4 7990 -4 9030 -5 1150 —3 1370 Cl КМ-1 1203 11014 8 6150-2 7 7792-2 8 2708—2 8 5075-2 4 7742—2 C2 КМ-2 —6 4500i—4 -5 8375-4 -4 4375-4 -4 0750—4 -4 5875-4 —4 5875-4 -2 2750-4 Сз км 3 1 2080—6 1 0830-6 8 1250-7 7 7083—7 9 1670—7 8 7500-7 3 9583—7 e< — -0,1200 -0 1200 -0 1000 -0 1000 -0 1200 -0 1100 -00900 c5 — 5 0000—3 2 5000-2 2 0830-2 2 7500-2 4 1160-2 3 8100-2 3 1180—2 ce — 0 0150 0 0075 0 0063 00038 0 0014 0 0012 0 0010 к — -1 1975-2 -9 9000-3 -7 6800-3 -5 6000—3 -4 9630-3 —4 1100-3 -3 0300—3 u КМ-1 9 9830-5 8 2120-5 6 3620-5 4 6670—5 4 1360—5 3 4630-5 2 5320-5 к КМ-2 0 3 1250-9 3 1250-9 0 0 -3 1250-9 -5 5560-10 . 4 ГОСТ 25645.115—84 Коэффициенты Фиксированный уровень солнечной активности F0, 10 Вт/м2 Гц) Обозначение Размерность 75 100 125 150 175 200 250 а\ -155605 -156408 -15 2229 -16 9752 -17 3045 — 18 2660 -19 2782 KM“1/2 0 8248 0 7754 0 7569 0 6736 0 6382 0 5797 05118 аг KM 769132 67 9162 55 8165 85 4440 81 9596 100i9417 116 5792 d$ _ —0 1721 -0 1721 -0 1721 -0 1721 -0 1721 -0 1721 —0 1721 d\ KM-1 5 7560—3 5 7560-3 5 7560—3 5 7560-3 5 7560-3 5 7560-3 5 7560—3 di KM-2 -3 6350-6 -3 6350-6 -3 6350-6 —3 6350-6 —3 6350—6 —3 6350-6 —3 6350—6 h -0 8607 -0 7540 -0 5700 -0 4760 —0 2920 -0 3113 —0 3307 b{ KM-1 7 8610—3 6 8500—3 5 2500-3 4 4000—3 2 8000—3 2 8390-3 2 8780—3 b KM-2 -5 7110-6 -4 6000-6 -3 0000-6 2 4000-6 —8 0000-7 —1 0890-6 —1 3780—6 Сл 12791 1 2791 1 2903 1 2903 0 2057 0 2057 0 0015 V Ci KM-1 —1 5760—2 —I 5760—2 -1 5470-2 -1 5470—2 —2 9120-3 —2 9140—3 —2 3990—4 i C2 KM-2 6 4990—5 6 4990-5 5 9640—5 5 9640-5 1 7390—5 1 7390-5 7 0060—6 C4 KM-3 —5 1450-8 —5 1450—8 -4 5030-8 -4 5030-8 -8 5649-9 —8 5649-9 —5 9990—10 fla 15 1 5 1 5 1 5 15 15 1 5 П\ KM-1 0 006 0 006 0 006 0 006 0 006 0 006 0 006 рад 0 5411 0 5515 0 5585 0 5585 0 5585 05585 0 5585 Co -2152 — 2162 -1486 -1495 —8 1900—2 -8 2860-2 — 2048 Cl KM-1 4 1670-3 4 0860-3 3 2630-3 3 1820-3 2 3580—3 2 2780-3 3 5960—3 C2 KM-2 1 5870—6 1,2700—6 3 1430-6 2 8250—6 4 6980—6 4 3810—6 — 1 5870—6 в* KM-3 -1 6510-9 -1,5870-9 -3 4290—9 —3 3650-9 —5 2060—9 -5 1430-9 3 1750—0 1? e* „__ —01200 -0 1200 -0 1000 -0 1000 —0 1200 -0 1100 —0 0900 e$ ___, 5 0000—3 2 5000-2 2 0830—2 2 7500-2 4 4160—2 3 8100-2 3 1180—2 __ 0 0150 0 0075 0 0063 0 0038 0 0014 0 0012 0 0010 о /q _ -1 6980-2 —1 2490—2 -7 8790-3 -4 8820-3 -5 1950-3 —5 0170-3 —5 4550—3 L KM-1 1 4480—4 1 1110-4 7 2580—5 4 6920—5 4 6640-5 4 2820-5 4 2730—5 l h KM-2 -9 5350-8 -7 7060-8 -3 6580-8 -1 7420-8 -2 1640—8 -2 1320—8 -2 2730—8 ГОСТ 25645115—84 С. 5 Коэффициенты Фиксированный уровень солнечной активности FQ 10 22 Вт/(м* Гц) Обозначение Размерность 75 100 125 150 175 200 250 fli _ —33 2283 -32 7731 -316715 -29 7592 -28 8463 -262994 -23 6627 а2 КМ -1/2 01784 0 1899 02265 0 2948 0 3140 0 3817 0 4231 йг км 555 0636 5842450 571 5408 528 3389 509 7230 434 1220 336 4318 do — 1 0204 10204 1 0204 1 0204 1 0*204 10204 10204 d\ КМ-1 2 4990-3 2 4990—3 2 4990-3 2 4990-3 2 4990-3 2 4990-3 2 4990-3 d2 КМ-2 -1 5190-6 —1 5190-6 -1 5190-6 -1 5190-6 -1 5190-6 —1 5190-6 -1 5190-6 bo — 0 7833 0 7250 0 6100 0 0933 —0 3333 -0 4333 -0 1750 bi КМ"1 2 8610-3 2 6750-3 2 3430-3 3 0380-3 3 5220-3 3 5220—3 2 6420—3 b2 КМ-2 -1 9440-6 -1 7500-6 -1 4330-6 -17110-6 -1 8890-6 —1 8890-6 -1 4170—6 CQ — —4 4000 -4 4000 -8 9800 -8 9800 -15 7800 -157800 —9 7500 Cl КМ-1 30240-2 30240-2 4 0870—2 4 0870-2 5 7570-2 5 7570-2 3 3830—2 c2 КМ-2 —3 2830—5 -3 2830—5 -3 9500-5 -3 9500-5 -5 3220-5 -5 3220-5 -2 6940-5 C3 КМ-3 1 0120-8 1 0120—8 1 1230-8 1 1230-8 1 5120-8 1 5120-8 6 4810—9 По — 15 15 15 1 5 1 5 15 15 П\ КМ-1 0 006 0 006 0 006 0 006 0 006 0 006 0 006 Фг рад 0 5411 0 5515 0 5585 0 5585 0 5585 0 5585 0 5585 Со ■— -3 8000 -3 7000 -37000 -4 4000 —3 6000* -3 6000 1000 Cl КМ-1 1 9720-2 1 7830-2 1 7500-2 1 9810-2 1 6940-2 1 6530-2 26390—3 С2 КМ”2 —1 8330-5 -1 5060-5 -1 5000-5 -1 8060-5 -1 5560-5 -1 5280-5 —2 7780—7 Сз КМ-3 4 9380-9 3 5800—9 3 7040-9 49380-9 4 3210-9 4 3210-9 —6 1730—0 е4 — -01200 -0 1200 -0 1000 -0 1000 -0 1200 -01100 -00900 С5 — 5 0000-3 2 5000—2 2 0830-2 2 7500-2 4 1160-2 3 8100-2 3 1180—2 С6 — 0 0150 0 0075 0 0063 0 0038 0 0014 0 0012 0 0010 to — 1 0830-2 8 3170-3 4 6670-3 -1 3330—2 -3 5000-3 -1 5000-3 -8 3330—4 lx КМ-1 6 6940-5 48370-5 4 6060-5 7 1670-5 4 3170-5 3 2500-5 2 8170—5 h КМ-2 —4 2770—8 -3 0390-8 -27220-8 -3 5180-8 —2 0560-8 -1 3890—8 -1 1290—8 Примечание к табл 1—3 Число со знаком минус, стоящее после значения параметра, является показателем степени десяти — сомножителя зна/чения параметра С в ГОСТ 25645 115—84 Значения множителя A(D) Таблица 4 D, сут Аф) D, сут Аф) D, сут Аф) 0 —0,028 130 0013 260 0015 10 —0 045 140 —0 037 270 0 070 20 —0 047 150 —0 086 280 0 115 30 —O.Q35 160 —0 128 290 0 144 40 —0 011 170 —0 162 300 0 155 50 0 022 180 —0 185 310 0 145 60 0 057 190 —0 199 320 0 120 70 0 090 200 —0 202 330 0 084 80 0 114 210 —0 193 340 0 044 90 0 125 220 —0 173 350 0 006 100 0 118 230 —0 140 360 —0 023 ПО 0 096 240 —0 096 370 —0 040 120 0.060 250< —0 042 — — (Измененная редакция, Изм. № 1) Значения множителя A(D) с шагом 10 сут приведены в табл. 4 (на промежуточные моменты времени значения множителя A(D) находят линейной интерполяцией). 2 ТАБЛИЦЫ ПАРАМЕТРОВ ПЛОТНОСТИ АТМОСФЕРЫ Параметры плотности атмосферы для семи фиксированных уровней солнечной активности в функции высоты приведены в табл. 5—11: рн — плотность ночной атмосферы, вычисленная по формуле Ко — величина, характеризующая вековое изменение плотности ночной атмосферы в 11-летнем цикле солнечной активности * h + /2 -h2; § К i — амплитуда суточного эффекта, вычисленная по формуле К[ — cQ + Cj • /i-bca • /i2-f с3 • h3-, f / f К2, Кз, К4 —величины, характеризующие, соответственно, влияние полугодового эффекта, радиоизлучения Солнца и геомагнитной активности, вычисленные по формулам: ^С2 ■ h-\-d^ • h2\ -hl\ K'4=e0+e^h-\-e2^h2 + e^h3. Для fW'S-lO-22 Вт/(м2.Гц) И, КМ Рн» кг/м3 «о ч *3 *4 120 2 4402-08 &00000 0 00000 0 00000 0 00000 0 00000 140 3 4105—09 0 00200 011045 0 38212 0.11653 0 27675 160 1 0237-09 0 00400 017683 0 63048 0 23377 0 44597 180 3 9507-Ю 0 00599 0 24934 0 74508 0 35177 0 56306 200 1 8180-10 0 00817 031510 0 83370 048306 0 66847 250 3 3192-11 0 01326 0 59707 1 03971 0 74761 0 89994 300 7 6469—12 0 01788 1 01105 1 22755 0 98361 1 13315 350 2 0616—12 0 02202 1 51846 1 39721 1 19105 1 36687 400 6 2383—13 0 02568 2 08070 1 54870 1 36994 1 59986 450 2 0643-13 0 02887 2 65919 1 68201 1 52027 1 83087 500 7 3403—14 0 03158 3 21535 1 79715 1 64205 2 05868 550 2 7701-14 003382 3 71058 1 89411 1 73527 2 28203 600 1 0993-14 0 03557 410630 1 97290 1 79994 2 49970 650 6 3936-15 0 03627 4 16453 2 00297 1 82161 2 62967 700 4 2452—15 0 03673 415246 2 02539 1 83344 2 71603 750 3 0124—15 003698 4 08250 2 04021 1 83555 2 76259 800 2 2289-15 0 03701 3 96224 2 04744 1 82794 2 77306 850 1 6985-15 0 03683 3 79927 204707 181061 2 75112 900 1 3234-15 0 03643 360118 2 03911 1 78356 2 70050 950 1 0494—45 003582 3 37556 2 02355 1 74679 2 62489 1000 8 4405-16 0 03500 3 13000 2 00040 1 70030 2 52800 1050 6 8704-16 003396 2 87209 1 96965 1 64409 2 41353 1100 5 6495-16 0 03271 2 60942 1 93131 1 57816 2 28518 1150 46866-16 0 03125 234958 1 88537 1 50251 2 14666 1200 3 9179-16 0 02957 210016 1 83184 1 41714 2 00166 1250 3 2976-16 0 02768 1.86875 1 77071 1 32205 1 85391 1300 2 7926-16 0 02557 1 66294 1 70199 121724 1 70709 1350 2 3779-16 0 02325 1 49032 1 62567 1 10271 1 56491 1400 2 0349-16 0 02072 1 35848 1 54176 0 97846 1 43107 1450 1 7493-16 0 01797 1 27501 1 45025 0 84449 1 30929 1500 1 5100-16 0 0150] 1 24750 135115 070080 1 20325 8 ГОСТ 25645.115—84 Я, КМ Рн. кг/м3 *0 «2 *3 «4 120 2.4403—08 0 00000 0 00000 0 00000 0.00000 0.00000 но 3.6339-09 0.00166 011045 0 38179 0.11449 0.26685 160 1.1062-09 0 00332 017683 063007 0.22448 0.43137 180 4.3011-10 о. 00498 0.24934 0 74457 0 33001 0 54476 200 2.1311-10 0 00665 031510 0 83370 0.43200 0.63910 250 4 5157-11 0 01047 0.59707 1.03971 0.67100 0.85988 300 1.1712-11 0 01390 1.01105 1.22755 0.88700 1.08105 350 3 4909—12 0 01696 1 51846 1 39721 1 08000 1.30143 400 1.1538-12 0,01962 2.08076 1.54870 1.25000 1.51983 450 4.1336-13 0 02190 2.65919 1 68201 1 39700 1.73506 500 1 5808-13 0.02380 3 21535 1 79715 1 52100 1.94593 550 6 3818-14 0 02530 3 71058 1 89411 1 62200' 2 15124 600 2 6976-14 0 02643 410630 1.97290 1 70000 2 34981 650 1 2291-14 0 02692 416453 2 00297 1 72438 2 50981 700 7.4306—15 0 02728 415246 2 02539 1 74000 2.62954 750 4.9722-15 0,02750 4.08250 2.04021 1 74688 2.71156 800 3 5232-15 0 02756 3 96224 2 04744 1 чж 2 75856 850 2.5934-15 002747 379927 2.04707 1.73438 2 77322 900 1 9627—15 0 02723 3 60118 2 03911 1.71500 2 75822 950 15173-15 0 02684 3 37556 2 02355 1 68688 2 71625 1000 1.1932-15 0 02630 3 13000 2.00040 1 65000 2 65000 1050 9 5165-16 0 02560 2 87209 1 96965 1 60438 2.56215 1100 7.6802-16 0 02475 2 60942 1.93131 1 55000 2 45538 1150 6.2616—16 0 02375 2 34958 1 88537 1 48688 2.33238 1200 5.1504-16 0 02260 2 10016 1.83184 1 41500 2.19584 1250 4 2694—16 0 02130 К86875 1.77071 1 33438 2.04844 1300 3 5638-16 0 01984 1 66294 1.70199 1 24500 1.89286 1350 2.9932-16 0 01823 1 49032 1 62567 1 14688 1 73179 1400 2.5282-16 0 01647 1 35848 1 54176 1 04000 1.56792 1450 2.1463—16 0.01456 1.27501 1.45025 0 92437 1.40393 1500 1.8307-16 0 01249 1.24750 1.35115 0.80000 1.24250 ГОСТ 25645.115—84 С. 9 #, км Рн, кг/м3 *0 *2 *3 *4 120 2 4404—08 0 00000 0 00000 0 00000 000000 0 00000 140 3 8112-09 000129 0 06064 0 38179 010253 0 24100 160 1 2529-09 000258 012632 063007 019515 0 40000 180 5 2126-10 0 00387 017337 0 74457 0 27791 0 51600 200 2 7113-10 0 00517 0 22166 083370 036000 0 60229 250 6 3042-11 000798 044671 1 03971 0 55500 081001 300 1 7524-11 0 01060 0 80109 1 22755 0 73500 1 02059 350 5 5234—12 0 01304 1 25104 1 39721 0 90000 1 23145 400 1 9137—12 0 01530 1 76278 1 54870 1 05000 1 44002 450 7 1437-13 0 01737 2 30254 1468201 1 18500 1 64374 500 2 8337-13 0 01927 2 83655 1 79715 ] 30500 1 84003 550 1 1826-13 0 02097 3 33103 1 89411 1 41000 2 02631 600 5 1530-14 002250 3 75222 1 97290 1 50000 2 20002 650 2 3198-14 0 02311 3 98079 2 00297 L52751 2 35471 700 1 3239-14 002357 4 12589 2 02539 1 54793 2 47047 750 8 3692-15 0.02390 419141 204021 156119 2 55013 800 56227-15 0 02409 418576 2 04744 1 56728 2 59645 850 3 9381—15 0 02415 4 11737 2 04707 1 56621 2 61222 900 2 8445-15 0 02407 3 99467 2 03911 I 55797 2 60022 950 2 1043-15 0 02386 3 82607 2 02355 1 54257 2 56322 ЮОО 1 5872-15 0 02351 3 62000 2 00040 1 52000 2 50400 1050 1 2165—15 0 02302 3 38488 1 96965 1 49027 2 42534 1100 9 4507-16 0 02240 3 12913 1 93131 1 45337 2 33002 1150 7 4285-16 0 02164 2 86118 1 88537 1 40931 2 22082 1200 5 8989-16 002074 2 58944 1 83184 1 35808 2 10051 1250 4 7267-16 0 01971 232234 1.77071 1 29969 1 97188 1300 3 8179—16 0 01854 2 06831 1 70199 1 23413 1 83769 1350 31061-16 0 01724 1 83576 1 62567 1 16141 1 70073 1400 2 5436-16 001580 1 63312 1 54176 1 08152 1 56378 1450 2 0954-16 0 01422 1 46881 1 45025 099447 1 42961 1500 1 7356-16 001251 1 35125 135115 0 90025 1 30100 10 ГОСТ 25645.115—84 Пара* Для Fo=L50h1O' -и Вт/(м2-Гц) И, км рн, кг/ч8 К 120 2 4402-08 000000 140 3 9752—09 0 00093 160 1 3650-09 0 00187 180 5 8911—10 0 00280 200 3 0775—10 0 00381 250 7 3551-11 0 00576 300 2 1589—11 0 00763 350 7 2620—12 0 00941 400 2 6969—12 001110 450 1 0807—12 0 01270 50Q 46022—13 0 01422 550 2 0606-13 0 01565 600 9 6237—14 0 01700 650 45195—14 001839 700 2 4537-14 0 01960 750 1 4492-14 0 02063 800 9 0585-15 0 02149 850 5 9028-15 002217 900 3 9722—15 0 02268 950 2 7431-15 0 02301 1000 1 9353-15 0 02316 1050 1 3902—15 0 02314 1100 1 0142-15 0 02294 1150 7 5002-16 0 02256 1200 5 6129-16 0 02201 1250 4 2452—16 0 02129 1300 3 2415—16 002039 1350 2 4965-16 0 01931 1400 1 9378—16 0 01806 1450 1 5150-16 0 01663 1500 1 1923—16 0 01502 атмосферы К2 *3 *4 0 00000 0 00000 оооооо 035179 0 08902 0 22004 0 63007 0 16842 0 37304 0 74457 0 23814 0 49604 0 83370 0 30800 0 57298 1 03971 0 47400 0 76998 1 22755 0 62800 0 96850 1 39721 0 77000 1 16599 1 54870 0 90000 1 35994 1 68201 1 01800 1 54783 1 79715 1 12400 1 72713 1 89411 1 21800 1 89531 1 97290 1 30000 2 04986 2 00297 1 34513 2 20225 202539 1 38154 2 31133 2 04021 1 40939 2 38197 2 04744 1 42869 2 41786 204707 1 43943 2 42270 2 03911 1 44162 2 40020 2 02355 1 43525 2 35407 2 00040 1 42033 2 28800 1 96965 1 39685 2 20570 1 93131 1 36482 2 11088 1 88537 1 32423 2 00723 1 83184 1 27509 1 89846 1 77071 1 21739 1 78828 1 70199 1 15114 1 68039 1 62567 1 07633 1 57848 1 54176 0 99297 1 48627 1 45025 0 90105 1 40746 1 35115 0 80058 134575 ГОСТ 25645.115—84 С 11 И, км Рн кг/м* *о *2 *3 < 120 2 4402—08 0 00000 0 00000 0 00000 0 00000 0 00Q00 140 4 1024-09 0 00083 0 07095 0 38179 0 05950 0 20004 160 1 4314—09 0 00165 0 13595 0 63007 012100 0 34108 180 6 2542—10 000248 0 19514 0 74457 018450 0 46713 200 29171—10 0 00327 0 25038 0 83370 0 23600 0 53597 260 7 6445-11 0 00511 0 43075 1 03971 0 35800 0 71988 300 2 4187-11 0 00685 0 66595 1 22755 0 47600 0 90776 350 8 6799-12 0 00848 0 94955 1 39721 0 59000 1 09570 400 3 4148-12 0 01000 1 27514 1 54870 0 70000 1 27980 450 1 4421—12 001141 1 63629 1 68201 0 80600 1 45615 500 6 4463-13 0 01272 2 02658 1 79715 0 90800 1 62085 550 3 0198—13 001391 2 43957 1 89411 1 00600 1 77000 600 1 4717-13 0 01500 2 86886 1 97290 1 10000 1 89968 650 7 0567—14 0 01587 3 30738 2 00297 1 15790 2 02355 700 3 8252—14 0 01664 3 62736 2 02539 1 20649 211570 750 2 2382-14 001731 3 84000 2 04021 1 24564 217542 800 1 3815-14 0 01788 3 95664 2 04744 1 27534 2 20595 850 8 8756-15 0 01834 3 98862 2 04707 1 29560 2 21053 900 5 8846-15 0 01870 3 94728 2 03911 1 30641 219241 950 4 0026-15 0 01896 3 84396 2 02355 1 30778 215482 1000 2 7812-15 001911 3 69000 2 00040 1 29970 210100 1050 1 9679-15 001916 3 49674 1 96965 1 28218 2 03420 1100 1 4144—15 001911 3 27552 1 93131 1 25521 1 95765 1150 1 0306-15 0 01895 3 03768 1 88537 1 21880 1 87460 1200 7 6016-16 0 01870 2 79456 1 83184 1 17294 1 78829 1250 5 6678-16 0 01834 2 55750 1 77071 1 11764 170195 1300 4 2674-16 0 01787 2 33784 170199 1 05289 1 61884 1350 3 2415-16 001731 2 14692 1 62567 0 97870 1 54218 1400 2 4822-16 0 01664 1 99608 1 54176 0 89506 1 47522 1450 19149-16 0 01587 1 89666 1 45025 0 80198 1.42121 1500 1 4873-16 0 01500 1 86000 135115 0 69945 1.38338 Для F0=20CM 0 -и Вт/(м2Тц) И, км рн, кг/мз *0 «2 *3 < 120 2 4398-08 0 00000 0 00000 000000 0 00000 000000 140 4 2346—09 0 00068 0 07095 0 38179 0 06167 0 20500 160 1 4975-09 0 00135 013595 063007 011646 0 33700 180 6 6076-10 0 00202 019514 0 74457 016431 0 43800 200 3 5725-10 0 00269 0 25058 0 83370 0 21294 050684 250 9 6598—11 0 00436 0 43100 1 03971 0 33039 0 68009 300 32111-11 0 00591 0 66625 1 22755 044239 0 85597 350 12176-11 0 00736 0 94990 1 39721 0 54895 1.03061 400 5 0691-12 0 00870 1 27554 1 54870 0 65006 1 20015 450 2 2С48-12 000993 1 63674 1 68201 0 74573 1 36074 500 1 0701-12 0 0Н06 2 02708 1 79715 0 83595 1 50852 550 5 2923-13 0 01208 2 44012 1 89411 0 92073 1 63963 600 2 7194-13 0 01300 2 86946 1 97290 1 00006 1 75021 650 1 3621—13 0 01376 3 30738 2 00297 1 05790 1 87535 700 7 3591-14 001444 3 62736 2 02539 1 10649 1 96590 750 4 2035-14 001506 3 84000 2 04021 1 14564 2 02542 800 2 5061—34 0 01561 3 95664 2 04744 1 17534 2 05715 850 1 5462-14 0 01609 3 98862 204707 1 19560 2 06433 900 9 8134—15 001650 3 94728 2 039а 1 20641 2 05021 950 6 3783—15 0 01684 3 84396 2 02355 1 20778 2 01802 1000 4 2312-15 001711 3 69000 2 00040 1 19970 1 97100 1050 2 8572-15 0 01731 3 49674 1 96965 1 18218 1 91240 1100 1 9597-15 0 01744 3 27552 193131 1 15521 1 84545 1150 1 3630-15 0 01751 3 03768 1 88537 1 11880 1 77340 1200 9 5991-16 0 01750 2 79456 1 83184 1 07294 1 69949 1250 6 8367-16 0 01742 2 55750 1 77071 1 01764 1 62695 1300 4 9194-16 0 01728 2 33784 1 70199 0 95289 1 55904 1350 3 5732-16 0 01706 2 14692 1 62567 0 87870 1 49898 1400 2 6178-16 0 01678 1 99608 1 54176 0 79506 1 45002 1450 1 9332-16 001642 1 89666 1 45025 0 70198 1 41541 1500 1.4382—16 0 01600 1 86000 135115 0 59945 1 39838 ГОСТ 25645 115—84 С 13 Для Fo—250'Ю^22 Вт/(мЧц) Я, км Рн. кг/ма < < *2 *3 < 120 2 4400-08 000000 0 00000 0 00000 0 00000 0 00000 140 4 4429—09 000050 006630 0 38179 0 04626 017404 160 1 6242-09 000101 013220 0 63007 0 09381 0 29904 180 7 3608-10 0 00151 0 18090 0 74457 014256 039404 200 3 8814-10 000218 0 22896 0 83370 018978 0 45346 250 1 1264-10 0 00381 0 37003 1 03971 0 30268 0 59997 300 40627-11 000532 0 54387 1 22755 0 40868 0 73974 350 1 6719-11 000672 0 75005 1 39721 0 50780 0 87301 400 7 5359-12 000800 0 98811 1 54870 0,60002 1 ооооо 450 3 6346-12 000917 1 25759 1 68201 0 68535 1 12096 500 1 8483-12 0 01023 1 55806 1 79715 0 76380 J 23614 550 9 8102-13 001117 1 88906 1 89411 0,83536 1 34576 600 5 3957-13 001200 2 25014 1 97290 090002 1 45006 650 2 8915-13 001271 263719 200297 0 94362 1 52845 700 1 6267-13 001335 2 95338 2 02539 098007 1 59944 750 9 5089-14 001394 3 20292 2 04021 1 00944 1 66256 800 5 7370-14 001448 3 39067 2 04744 103172 171735 850 3 5546-14 0 01495 3 52149 2 04707 1 04692 1 76334 900 2 2532-14 001537 3 60025 2 03911 1 05503 1 80007 950 1 4568—14 001574 3 63180 2 02355 1 05606 1 82706 14 ГОСТ 25645.115—84 И» км рн, кг/м* К *1 1000 9.5846-15 0.01605 3.62100 1060 6.4045-15 0.01630 3.57272 1100 4.3394-15 0.01649 3.49181 1150 2.9773-15 0.01663 3.38314 1200 2.0662-15 0.01671 3.25157 1250 1.4489-15 0.01674 3.10195 1300 1.0259-15 0.01671 2.93916 1350 7.3277-16 0.01662 2.76804 1400 5.2773-16 0.01648 2.59346 1450 3.8298—16 0.01628 2.42029 1500 2.7991-16 0.01602 2.25338 Примечание к табл. 5—11, Число со знаком минус, степени десяти — сомножителя значения параметра. Разд. 2. (Измененная редакция, Изм. № 1). Продолжение табл. 1] «2 *3 < 2.00040 1.05000 1.84390 1.96965 1.03686 1.85007 1.93131 1.01663 1.84514 1.88537 0.98932 1.82862 1.83184 0.95492 1.80007 1.77071 0.91344 1.75902 1.70199 0.86487 1.70501 1.62567 0.80922 1.63757 1.54176 0.74648 1.55624 1.45025 0.67666 1.46056 1.35115 0.59975 1.35006 стоящее после значения параметра, является показателем ГОСТ 25645.118—84 С. 15 ПРИЛОЖЕНИЕ 1 Справочное I Значения предельных относительных отклонении плотности атмосферы m 1 п И &Р max по высоте от ее средних значений за период в годы низкого, среднего и высокого уровней солнечной активности приведены в табл 1. В соответствии с ГОСТ 25645 102—83 считают низкий уровень солнечной активности — минимум, средний — фаза роста и спада, а высокий — максимум солнечной активности Таблица I Предельные относительны^ отклсмсния плотности h, км Низкий СОЛЛГЧЮЙ j урове НЬ активности Ср дпий уровень солчсчюи активности Высокий уровень солнечной активности 8Pmin' * ^Ргпах* °° in’ % 6?шах- 0/° ®Рш in- % 5Ршах- 160 — 15 + 20 — 11 + 18 — 10 + 15 200 —17 +-25 — 17 + 22 —15 +20 240 —20 + 30 —22 +27 —20 + 30 300 —27 + 35 —30 + 40 —25 + 40 400 —45 +-50 —43 + 52 —35 +50 500 —55 + 60 —50 +57 —45 +50 600 —60 + 65 —57 + 62 —50 + 60 900 —35 +55 —50 + 60 —55 +60 1200 —40 + 50 —42 , +52 —47 +55 1500 —35 1 + 47 —37 + 47 —37 + 50 Значения 6pmin и бртах получены с использованием формулы (1) настоящего стандарта и данных об индексах солнечной активности F и К р за период 1964—1987 гг. 2. Значения предельных относительных отклонений изменения плотности атмосферы за одни сутки бр ( по высоте за период в годы низкого, среднего и высокого уровня солнечной активности приведены в табл. 2. 3. Изменения отношения плотности атмосферы при различных фиксированных уровнях солнечной активности к плотности при уровне F0=75*lO~22 Вт/ (м2-Гц) по высоте указаны в табл. 3 Нормированные автокорреляционные функции невязок модели плотности атмосферы для низкого, среднего и высокого уровней солнечной активности (СА) приведены на черт. 1. Средние квадратические отклонения невязок модели плотности атмосферы, вычисленные по торможению ИСЗ в функции высоты (до 620 км), приведены на черт. 2. Данные табл. 1—3 используют при определении диапазона изменений баллистического коэффициента, как коэффициента согласования, и как исходные данные при априорной оценке точности прогнозирования движения ИСЗ. Таблица 2 Предельные относительные отклонения плотности на 1 сут , % Низкий уровень Средний уровень Высокий уровень h, км солнечной активности солнечной активности солнечной активности брдпНп ^дгпах 3рДШ!П 3Рдтах ®Рдт in 3Рд шах 160 —15 + 15 — 12 + 17 —10 + 10 200 —15 +20 —12 + 17 —15 +15 240 —17 +25 —17 + 17 —15 +15 300 —25 +35 —22 +27 —20 +25 400 —40 +50 —32 +40 —30 +35 500 .—45 +55 —37 +57 —45 +45 600 —45 +55 —40 +60 —45 +55 900 —25 +35 —27 +27 —40 +40 1200 —35 +40 —30 +42 —25 +30 1500 —35 +45 —37 +45 —30 +35 Таблица 3 Н, км Отношения плотности атмосферы Р„. Ю-22 Вт/(м**Гц) 100 125 150 200 250 160 1.080 1.223 1.333 1.463 1.586 200 1.161 1.478 1.710 1.989 2.270 240 1.288 1.776 2.165 2.751 3.355’ 280 1.398 2.058 2.648 3.705 4.856 320 1.498 2.332 3.166 4.879 6.873 360 1.591 2.603 3.723 6.308 9,533 400 17678 2.873 4.323 8.033 12.991 500 1.883 3.555 6.031 13.943 26.442 600 2.075 4.260 8.075 22.857 50.168 900 1.923 3.336 , 5.491 12.811 26.684 1200 1.828 2.770 3.989 7.239 12.485 1500 1.755 2.374 3.044 4.375 6.156 1—высокий уровень; 2—средний уровень; 3—низкий уровень. Черт. 1 Черт. 2 ПРИЛОЖЕНИЕ 2 Справочное При определении коэффициента аэродинамического сопротивления (Сх ^ ИСЗ приняты следующие основные исходные предположения режим обтекания поверхности ИСЗ разреженным газом верхней атмосферы — свободномолекулярный, набегающий поток газа, имеющий переносную скорость и — «+[(оХгЬ состоит из молекул с температурой Т, скорости которых распределены по закону Максвелла; молекулы набегающего потока газа испытывают однократные соударения с поверхностью ИСЗ, характер отражения молекул от поверхности ИСЗ— полностью диффузный (коэффициенты аккомодации нормального и касательного импульсов равны единице), при этом после отражения молекулы имеют скорости, распределенные по закону Максвелла, с температурой поверхности Tw; затенение одних элементов ИСЗ другими учитывают по аналогии с геометрической оптикой При принятых исходных предположениях проекции на орты пит вектора аэродинамической силы Рп и Р^ , отнесенные к скоростному напору pv2/2go и действующие на элементарную площадку dF соответственно равны* P„-eos© %(z) + 1-ferfz У Ml %(z) V-All ' Pr =sin0 M, Х(г) = ехр(—г2) + |^я z(l -f erf 2), erfz= JL- ( exp(—x2)dx, (1) it 0 z= |/"^/2/Vfooco-0, где M ^—v/a^ — число Маха набегающего потока, а— скорость звука в невозмущенном потоке, м/с, у — отношение удельных теплоемкостей, п — внутренняя нормаль к элементу dF k поверхности, © — угол между вектором скорости набегающего потока и внутренней нормалью, рад, erf г — интеграл вероятности Для определения коэффициента аэродинамического сопротивления С х (а) а ИСЗ, движущегося под углом атаки а, наружную поверхность его условно разбивают на К элементов простой геометрической формы (например, пластина, цилиндрическая, коническая, сферическая поверхности и др ) Для каждого ^-го-элемента (£=1, 2, . , К) вычисляют силу аэродинамического сопротивления путем интегрирования по обтекаемой поверхности F & выражения проекции на ось ха силы, действующей на элементарную площадку dF ь . Результаты вычислений, полученных для К элементов, суммируют и относят к характерной площади ИСЗ Л* —const и аэродинамическому скоростному напору рv2/2gQ к Сха (a)=2RXa /FMpu2= 2 I (Рлсо,0+Р. sine)dFk/Fv, (2) *=i Fk где R x — проекция суммарной аэродинамической силы на ось ха (cos (ха, а *0 = 1). При определении в момент времени t коэффициента аэродинамического сопротивления Сх (ot (t)) неориентированных ИСЗ (вращающихся относительно не- а -*■ которой оси спутника) в силу малости угловой скорости w вторую составляющую скорости v не учитывают. При этом для определения С х ^ в интервале времени от t до t+T,ip усредняют по времени значения коэффициентов сопро тивления: Т СХа= -4— С CXa(a(t))dt, (3) 1ДР о где Т др — период обращения ИСЗ. Погрешность при определении коэффициентов аэродинамического сопротивления ИСЗ по приведенной методике в общем случае не более 30 %. Для ИСЗ, близких по форме к сферической, погрешность в определении С х оценивают а меньше 7 %. ПРИЛОЖЕНИЕ 3 Справочное При баллистическом обеспечении полетов ИСЗ с целью повышения точности определения и прогнозирования орбит ИСЗ рекомендуют рассчитывать плотность с учетом текущих и прогнозируемых индексов F и Ар , а также суточного и полугодового эффектов. Обеспечение данными об индексах солнечной и геомагнитной активностей F и Кр возложено в настоящее время на Службу Солнца Кроме того, данные об этих индексах публикуют в специальной литературе При использовании этих данных следует учесть запаздывание изменения плотности атмосферы по отношению к временам соответствующих солнечной и геомагнитной активностей Для солнечной активности это запаздывание — 1,7 сут, а для геомагнитной — 0,6 сут. При отсутствии текущих данных об индексах солнечной активности плотность атмосферы рекомендуют рассчитывать с использованием прогнозируемых значений этих индексов или осредненных данных за предшествующий период. В последнем случае в качестве F принимают среднее значение F10 7 за предшест вующие 30—50 сут, в качестве Fai — среднее значение Г]07за предшествующие 81 сут, а в качестве — значение, ближайшее к F8U для которого есть коэффициенты модели. При отсутствии данных об индексе Кр его значение принимают равным 2,66667. При использовании в модели 3-часовых индексов геомагнитной активности коэффициенты е±, е$ и е6 принимаются в соответствии с табл 1 В формулу (1) в выражение Ка в этом случае вместо индекса Кр подставляют индекс Кр/- , вычисляемый по формуле где Кpj — 3-часовой планетарный индекс геомагнитной активности, д_ при /=?ь0; г== (°>3 при Л>0 1о при / = 0 [0,7 при Д<0. В этом случае запаздывание для геомагнитной активности — 0,25 сут. Таблица 1 Коэффи циенты _12 Фиксированный уровень солнечюй активности Fo> Ю " Вт/(м2-Г О 75 100 125 150 175 200 250 е4 е5 е6 (Изм< — 12 .025 .0075 ?ненная ре, — 11 1 — 10 —09 —08 02625 02617 .02542 1 .02333 005625 00425j 0031251 .0025 дакция, Изм. № 1). — 07 02125 001875 0 1 ° о о 1+ 01 сл Значения индекса Кр, дят по табл 2. используемого в модели плотности атмосферы, нахо- Таблица 2 Таблица перевода Кр ЛР *Р ЛР лр *Р лр 0 0 2+ 9 5— 39 7 132 0 + 2 з- 12 5 48 7 + 154 1- 3 3 15 5+ 56 8' 17Э 1 4 3+ 18 6- 67 8 207 1 + 5 4" 22 6 80 8+ 236 2- 6 4 27 6+ 94 9- 300 2 7 4 + 32 7- 311 9 400 Индекс А р — по ГОСТ 25645 302—83 Знаки « + » и «—» справа от значения Кр означают, что к значению, заданному цифрой, нужно прибавить, соответственно, +0,3333 или —0 3333 т е например, 3+ = 3 + 0,3333 = 3,3333, а 4~ = 4+(—0,3333) =3,6667 Для нахождения промежуточных значений используют линейную интерполяцию. При отсутствии данных об индексе К р допустимо использование индекса^ осредненного по нескольким станциям, находящимся на территории СССР В качестве коэффициента согласования между плотностью атмосферы, рассчитанной в соответствии с настоящим стандартом, и фактической плотностью* атмосферы при баллистическом обеспечении ИСЗ рекомендуют использовать баллистический коэффициент ИСЗ, периодически уточняемый по результатам определения параметров орбит ИСЗ Наряду с баллистическим коэффициентом в качестве коэффициента согласования допускают использовать изменение периода обращения ИСЗ за виток, скорость изменения периода обращения, изменение большой полуоси орбиты за виток (Измененная редакция, Изм. № 1) ПРИЛОЖЕНИЕ 4 Справочное \ Для высот км плотность р рассчитывают по формуле р—ALexp(—Ku(h~hl)-\'K2i(h—hi)1) > где At Ки> Kzi — коэффициенты, приведенные в таблице, i ■— номер слоя 1 h t, км At, ki/m3 Kir 1/km K2l 1/km* I 0«Л<20 0 I2280-101 * 0 90764 10“1 —0 20452 10-2 2 20«ft<60 090130 10-1 016739 0 62669 10~3 3 60 0 31043-10-3 0 12378 —0 86999 10“3 4 100 ch 0 53675 10-6 017527 0 12870i-10-2 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 — значение индекса F8b D — время московское в сутках от начала года, PR — если PR>0, то при обращении к программе параметру АР присваивают значение Кр , если PRcO, то параметру АР присваивают значение Ар и оно переводится в Кр внутри программы (подпрограмма INTERP (АР, КР)) Выходные параметры RO ■— плотность атмосферы, кг/м5; RU .— плотность атмосферы, кгс*с£/м4 5 Используемые блоки COMMON Программа DENSIT использует три общих блока COMMON /QKOEF/Q (27,3)/YEAR/P(38) COMMON/UKOEF/U(8,4), GR, WZ, в которых должны содержаться все необходимые для вычисления константы для фиксированного уровня солнечной активности F0 6 Описания параметров блоков COMMON Q — двумерный массив размерности 27X3, содержащий коэффициенты модели плотности атмосферы для одного из фиксированных уровней F0 Первые элементы каждого столбца равны F0, далее следуют коэффициенты модели настоящего стандарта в соответствии с табл 1, 2 и 3 п 18, причем в первом столбце находятся коэффициенты для диапазона высот от 120 до 180 км, во втором — от 180 до 600 км, а в третьем — от 600 до 1500 км Р — массив, содержащий значения поправок полугодового эффекта из табл. 4 п I 8 настоящего стандарта, длина массива равна 38, П — двумерный массив размерности 8X4, содержащий коэффициенты таблицы, WZ — 0 7292115-10-4 (угловая скорость вращения Земли, рад/с), GR^= — 9 80665 кг/м3 5 6 (Измененная редакция, Изм. № 1). 7 Для присвоения нужных значений переменным из упомянутых общих блоков следует воспользоваться подпрограммой VKM(FSH) CALL VKM(FSH) Подпрограмма VKM(FSH) находит ближайшее к FSH значение F0 и засылает соответствующие этому 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, .0014D01, * —.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, .003383DO1, — .2694D—4, .648ID—3, * 1.5DO, .O06DO, .541052DO, * 1.0204DO, .00249900, —.1519D—5, * .IDO, .OQ2639DO, —.27780—6, —.6173D—9, * — .O9D0, .03118DO1, .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.D0v 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
. 12 ГОСТ 25645.115—84KW2Moo
X(f)
/=0, 1,.-. п;