- •Содержание
- •Определение характеристик топлива.
- •Выбор проектных параметров.
- •Проектировочный баллистический расчет.
- •Проверочный баллистический расчет.
- •Весовой расчет
- •Выбор оптимальных проектных параметров.
- •Весовой расчет ракеты при выбранных оптимальных проектных параметрах.
- •Определение основных геометрических характеристик
- •Определение тяговых характеристик
- •Объемный расчет ракеты Объемный расчет гч
- •Бак окислителя.
- •2 Ступень. Массовые расходы окислителя и горючего равны
- •Бак окислителя
- •Прикидочный расчет габаритов ду Расчет ду 1 ступени.
- •Расчет ду 2 ступени.
- •Объемный расчет приборного, хвостового и переходного отсеков
- •Определение центра тяжести ракеты на аут
- •Расчет нагрузок действующих на ракету в полете.
- •Расчет топливных баков ракеты Расчет обечаек топливных баков Бак окислителя 2-й ступени.
- •Бак горючего 2-й ступени.
- •Бак окислителя 1-й ступени.
- •Бак горючего 1-й ступени.
- •Расчет распорных шпангоутов. Распорные шпангоуты верхних днищ.
- •Распорный шпангоут нижнего днища бака горючего первой ступени.
- •Распорный шпангоут нижнего днища бака горючего второй ступени.
- •Распорный шпангоут нижнего днища бака окислителя первой ступени.
- •Распорный шпангоут нижнего днища бака окислителя второй ступени.
- •Расчет днищ топливных баков.
- •Нижнее днище бака окислителя 2-й ступени.
- •Нижнее днище бака горючего 2-й ступени.
- •Нижнее днище бака окислителя 1-й ступени.
- •Нижнее днище бака горючего 1-й ступени.
- •Верхние днища баков ракеты.
- •Расчет окантовок около круговых отверстий топливных баков.
- •Расчет окантовок бака горючего первой ступени
- •Расчет окантовок бака окислителя первой ступени
- •Расчет окантовок бака горючего второй ступени
- •Расчет окантовок бака окислителя второй ступени
- •Расчет тоннельных труб. Тоннельная труба 1-й ступени
- •Тоннельная труба 2-й ступени
- •Расчет фланцевого соединения крепления крышки люка-лаза
- •Расчет сухих отсеков ракеты
- •Расчет приборного отсека 2 ступени.
- •Расчет межбакового отсека 2 ступени.
- •Расчет хвостового отсека 2 ступени
- •Расчет приборного отсека 1 ступени.
- •Расчет хвостового отсека 1 ступени.
- •Расчет фермы переходного отсека
- •Расчет стыковочного шпангоута гч.
- •Расчет заклепок
- •Список использованных источников
Бак горючего 1-й ступени.
Определяем силы и моменты в наиболее характерных сечениях бака
Сечение 1 - 1` : Т1 = - 633218 Н, М1 = 37095 Н×м
Сечение 2 - 2` : Т2 = - 647983 Н, М2 = 18938 Н×м
Сдвигающая сила Q оказывает малое влияние на прочность бака и поэтому в дальнейшем не учитывается.
Выбираем материал обечайки бака - алюминиевый сплав АМг -6:
sв = 320 МПа , s0,2 = 280 МПа , Е = 6800 МПа
Эксплуатационное давление наддува рнад = 0,4 МПа.
Меридиональные и кольцевые усилия в отдельных точках рассмотрим в характерных точках.
Точка 1.
Меридиональное усилие:
fт ,fр - коэффициенты безопасности по усилию Т и давлению р.
Кольцевые усилия:
= 1,3 × 4×105× 0,9 = 4,68×105 Н/м
Точка 1`.
Меридиональное усилие:
Точка 2.
Меридиональное усилие:
где ризб = рнад + r×g×h×nx = 4×105 + 1000 × 9,81 × 3,632 × 6 = 6,137×105 Па
= 1,3 × 6,137×105 × 0,9 = 7,180×105 Н/м
Точка 2`.
Меридиональное усилие:
Так как на обечайку бака действуют большие усилия, то для уменьшения толщины обечайки, а как следствие массы баков, применяем вафельную конструкцию.
Будем считать, что оболочка изготовляется методом химического травления и подкреплена продольно-поперечным набором.
Определяем параметры силового набора
Толщина полотна равна
где коэффициенты равны k = 0,28 , b = 0,6 , j = 0,2, y = 5.
Примем d = 2 мм.
Толщина исходного листа
dисх = y×d = 5 × 2 = 10 мм.
Задаемся соотношением j1/j2 = 1
Шаг ребер определяем по формуле
Примем а = 0,15 м.
где d1 = d×[1 + kф×j1×(y - 1)] = 0,002×[1+0.16×0.2×(5-1)] = 0.00226 м
kф = 0,26 при r = h, k1 = 6.
h = dисх - d = 10 - 2 = 8 мм.
Шаг ребер в продольном направлении
b = a = 0,15 м.
Ширина ребер в поперечном и продольном направлении
Примем С = S = 0,005 м.
Определяем предельное напряжение
Для вафельной оболочки возможна местная потеря устойчивости. Под потерей местной устойчивости подразумевается потеря устойчивости полотна обечайки в отдельных ячейках, которая происходит хлопком.
где а0 = а - S - 2×r = 0,15 - 0,005 - 2×0,008 = 0,129 м
Расчетные напряжения в ячейках равны
Условие выполняется.
При проверке работоспособности вафельной оболочки используется понятие эквивалентной толщины, которая получается при равномерном распределении материала ребер по всей поверхности оболочки.
Проверку сохранения работоспособности будем производить в наиболее нагруженной точке обечайки (точке 2).
Меридиональные напряжения:
s1 = N1 / d1э = 2,197×105 / 0,00235 = 9,348×107 Па
где
Кольцевые напряжения
s2 = N2 / d2э = 7,180×105 / 0,00235 = 3,042×108 Па
где
Эквивалентные напряжения, действующие в обечайке
Запас прочности
h = sВ / sэкв = 320 / 269,925 = 1,18
Расчет распорных шпангоутов. Распорные шпангоуты верхних днищ.
Исходные данные:
Эксплуатационная нагрузка р = 0,4 МПа
Диаметр бака D = 1,8 м
Высота сферического днища hд = 0,3 м
Материал бака и шпангоута: АМг6: sВ = 380 МПа.
Определяем радиус сферического днища
Rсф = D2 / 8hд + hд / 2 = 1,82/8×0,3 + 0,3/2 = 1,5 м
Определяем угол b
b = arcsin R / Rсф = arcsin 0,9/1,5 = 37
Определяем усилия Т1, S1
В обечайке:
Т1 = p R / 2 = 4×105×0,9 / 2 =180 кН/м
В днище:
S1 = p Rсф / 2 =4×105×1,5 / 2 = 300 кН/м
Рис. 9 Схема распорного шпангоута
Задаемся размера шпангоутов:
А = 11,5 мм, В = 35,5 мм, С = 5,0 мм, Н = 3,0 мм, Н11 = 3,0 мм, К = 5,5 мм, О = 1,5 мм, М = 2,0 мм, N = 16,5 мм, l1 = 9,0 мм, l2 = 11,0 мм, l3 = 11,5 мм.
Определяем площади фигур
Фигура 1
F1 = A× B = 11,5 × 35,5×10-6 = 4,0825×10-4 м2
Фигура 2
F2 = l2 × H = 11,0×3 ×10-6= 3,3×10-5 м2
Фигура 3
F3 = l2 × C / 2 = 11,0×5×10-6 /2 = 2,75×10-5 м2
Фигура 4
F4 = К×N = 5,5×16,5 ×10-6= 9,075×10-5 м2
Фигура 5
F5 = l1 × H = 9,0×3×10-6 = 2,7×10-5м2
Фигура 6
F6 = l1 × M / 2 = 9,0×2,0×10-6/2 = 9,0×10-6 м2
Фигура 7
F7 =Н112× tg a /2 = 3,02×10-6 tg 53 / 2 = 5,9717×10-6 м2
Фигура 8
F8 = H11×( l3 - H11×tg a) = 3×10-3 (11,5×10-3 - 3×10-3×tg 53) = 2,256×10-5 м2
Фигура 9
F9=(K-O)×(l3-H11×tga)sina=(5,5-1,5)×10-3×(11,5×10-3-3×10-3×tg53)×sin53 =
= 1,20097×10-5 м2
Суммарная площадь равна S Fi = 6.36×10-4 м2
Определяем статические моменты фигур:
Фигура 1
Sx1 = 0
Sy1 = A B (B-H) / 2 = 11,5×35,5×10-6×(35,5-3)×10-3 / 2= 6,34×10-5 м3
Фигура 2
Sx2 = - l2 H (l2 + A) / 2 = - 11×3×10-6×(11 + 11,5)×10-3 = -3,713×10-7 м3
Sy2 = 0
Фигура 3
Sx3 = - C l2 (l3 / 3 +A / 2) / 2 = - 5×10-3×11×10-3×(11,5 / 3 + 11,5 / 2)×10-3 =
= -2,635×10-7 м3
Sy3 = C l2 (C / 3 +H / 2) / 2 = 5×10-3×11×10-3×( 5 / 3 + 3 / 2)×10-3 = 8,708×10-8 м3
Фигура 4
Sx4 = K N (K /2 + A / 2) = 5,5×16,5×10-6×(5,5 / 2 +11,5 /2)×10-3 = 7,714×10-7 м3
Sy4 = K N (N /2 - H / 2) = 5,5×16,5×10-6×(16,5 / 2 - 3 /2)×10-3 = 6,126×10-7 м3
Фигура 5
Sx5 = H l1 ( l1 / 2 +K + A / 2) = 3×9×10-6×(9. / 2 + 5,5 + 11,5 / 2)×10-3 = 4,253×10-7 м3
Sy5 = 0
Фигура 6
Sx6 = l1 M (l1 /3 + K + A / 2) / 2 = 9×2×10-6×(9 / 3 + 5,5 + 11,5 / 2 )×10-3 / 2 = 1,282×10-7 м3
Sy6 = l1 M (H /2 + M / 3) / 2 = 9×2×10-6×(3 / 2 + 2 / 3 )×10-3 / 2 = 1.95×10-8 м3
Фигура 7
Sx7 = H112 tg a ×(H11×sin a/3 + K +A/ 2) / 2 = 32×10-6 tg 57 (3×sin 57 / 3 + 5.5 + 11.5/ 2)×10-3 / 2 = = 7.195×10-8 м3
Sy7 = 7.606×10-8 м3
Фигура 8
= 3.318×10-7 м3
= 4.419×10-7 м3
Фигура 9
Sx9 = F9×(K + A/2 + S’x9 / F9) = (5.5×10-3 + 11.5×10-3 / 2 + 2.965×10-8 / 1.20097×10-5)×1.20097×10-5 = = 1.647×10-7 м3
Sy9 = F9(1/3×(l3 - H11 tg a) sin a + N - H / 2) = 1.20097×10-5(1/3×(11.5-3×tg57)×10-3×sin 57 + 16.5×10-3 - 3×10-3 / 2) = 2.042×10-7 м3
Суммарные статические моменты шпангоута равны
S Sx i = 1.259×10-6 м3
S Sy i = 8.075×10-6 м3
Определяем координаты центра тяжести шпангоута
x0 = Sy / F = 1.259×10-6 / 6.36×10-4 = 12.696 мм
y0 = Sx / F = 8.075×10-6 / 6.36×10-4 = 1,979 мм
Определяем моменты инерции фигур, входящих в шпангоут:
Фигура 1
I1 = A3×B / 12 = (11,5×10-3)3×35,5×10-3 / 12 = 4.4993×10-9 м4
Фигура 2
I2 = H×l33/12 + l2×H (l2 /2 + A/2)2 = 3×10-3×(11,5×10-3)3/12+11×10-3×3×10-3×
×(11×10-3/2 + 11,5×10-3/2)2 =4,5093×10-9 м4
Фигура 3
I3 = C×l3/12 - l2×C/2×(A/2)2 = 5×10-3×11,5×10-3/12 -11×10-3×5×10-3×
×(11,5×10-3/2)2/2 = -3,5463×10-10 м4
Фигура 4
I4 = N×K3 / 12 + K×N (K/2 + A/2)2 = 16,5×10-3×(5,5×10-3)3/12 +
+ 16,5×10-3×5,5×10-3(5,5×10-3/2 +11,5×10-3/2)2 = 6,7854×10-9 м4
Фигура 5
I5 = l13×H/12 + l1×H (l1/2 +K + A/2)2 = (9×10-3)3×3×10-3/12 +
+ 9×10-3×3×10-3-(9×10-3/2 + 5,5×10-3 + 11,5×10-3/2)2 = 6,879×10-9 м4
Фигура 6
I6 = l13×M / 12 -M×l1(A/2 + K)2 / 2 = (9×10-3)3×2×10-3/12 -
- 2×10-3×9×10-3×(11,5×10-3/2 + 5,5×10-3)2 / 2 = -1,01756×10-9 м4
Фигура 7
I7 = (Н114 tg a / 12) sin2 a + (H114×tg3 a) / 12 + H112 tg a (K + A/2)2 =
= 7,6722×10-10 м4
Фигура 8
Фигура 9
где k = (K-O + (l3 - H11 tg a) cos a) / [(l3 × H11 tg a) sin a] = 1,419
x22 = (l3 - H11 tg a) sin a = ( 11,5×10-3 - 3×10-3× tg 57) × sin 57 = 0.006 м
I9 = I’9 + F9 ( O + A/2)2 = -3,1645×10-9 + 1,20097×10-5×(1,5×10-3 + 11,5×10-3/2)2 = -2,5332×10-9 м4
Суммарный момент инерции равен S Ii = 2,610×10-8 м4
Вычисляем координаты точки m:
xm = 0 ym = K + A/2 - (N - H/2 - H11/2 cos a) / tg a = 1,825 мм
Определяем координаты точки n:
xn = S1×sin a / 107 = 3×105 sin 57 / 107 = 23,959 мм
yn = (- T1/107) + S1 cos a /107 + ym = 1,879 мм
Уравнение линии mn
Подставляем в уравнение х0, получаем y = 0,00191 м.
Определяем постоянные С1 и С2
С1 = x0 × cos a - sin a ×(y0 - ym) =
= 12,696×10-3 × cos 57 - sin 57 (1,98 - 1,825)×10-3 = 7,52 мм
С2 = х0 = 12,696 мм
Проверяем условие С2Т1 = С1S1
C2T1 = 0.0126×1,8×105 = 2268 Н×м/м
С1S1 = 0.00752 × 3×105 = 2256 Н×м/м
Так как С2Т1 ¹ С1S1, то в сечении шпангоута будет действовать крутящий момент.
Мкр = (C2T1 - С1S1) = 2268 - 2256 = 29,96 Н×м/м
Определяем напряжение, действующее в сечении шпангоута
где a0 = 90° - a = 90° - 57° = 33°
R1 = R - x0 = 0,9 - 0.0129 = 0,887 м,
xmax=A / 2+ K + l1 = 11.5 / 2 + 5.5 + 9 = 0.020 м
Определяем запас прочности шпангоута:
n = sB / s = 380 / 354,868 = 1.07