Отделение космического аппарата от орбитальной ступени
0 likes5,222 views
Рассматривается задача построения модели движения орбитальной ступени и космического аппарата на этапе работы средств отделения - толкателей. Приводится методика определения удлинения толкателя, записываются динамические и кинематические уравнения движения разделяемых тел.
![пружинный толкатель [1]
1 2 3 4
1 - шток; 2 - обойма внутренней пружины; 3 - корпус;
4 - технологическая заглушка.
5](https://image.slidesharecdn.com/pusher-130309045923-phpapp02/85/-6-320.jpg)
![пиротолкатель [1]
1 2
3 4
5
FT
x0
1 - пиропатрон; 2 - корпус; 3 - пассивный отделяемый блок;
4 - шток;
5 - активный отделяемый блок.
6](https://image.slidesharecdn.com/pusher-130309045923-phpapp02/85/-7-320.jpg)
![координатная форма
Координатная форма выражения для определения
перемещения штока толкателя имеет вид:
∆i =
n
(2)
p [A2T(r
(0)
2 − r
(0)
1 − A1ρ
(1)
1i ) + ρ
(2)
p ]
a
(1)T
i A1TA2n
(2)
p
Координатный столбец радиус-вектора точки контакта
относительно полюса C2 в базисе C2x2y2z2:
ρ
(2)
2i = A2T
(r
(0)
1 + A1
a
(1)
i ∆i − r
(0)
2 )
15](https://image.slidesharecdn.com/pusher-130309045923-phpapp02/85/-18-320.jpg)
More Related Content
What's hot (20)
Viewers also liked (6)
Similar to Отделение космического аппарата от орбитальной ступени (20)
More from Theoretical mechanics department (20)
Отделение космического аппарата от орбитальной ступени
- 1. отделение полезного груза от орбитальной ступени Основы синтеза механических систем Юдинцев В. В. 21 марта 2015 г. Самарский государственный аэрокосмический университет им. академика С. П. Королёва (СГАУ) 0
- 2. Средства отделения и разделения Математическая модель процесса отделения Определение перемещения штока толкателя Определение сил и моментов Уравнения движения 1
- 3. требования к системе отделения ка • Заданная скорость отделения полезного груза • Минимальная угловая скорость отделения полезного груза • Исключение взаимного соударения отработавшей ступени и полезного груза 2
- 4. средства отделения и разделе- ния
- 5. классификация средств разделения • Средства крепления и разделения: пирочеки, пирозамки, пиромеханические устройства, пироболты. • Средства отделения: пиротолкатели, пневмотолкатели, пружинные толкатели, ... • Сопутствующие элементы: электроразъемы, пневморазъемы, шпильки, ... 4
- 6. пружинный толкатель [1] 1 2 3 4 1 - шток; 2 - обойма внутренней пружины; 3 - корпус; 4 - технологическая заглушка. 5
- 7. пиротолкатель [1] 1 2 3 4 5 FT x0 1 - пиропатрон; 2 - корпус; 3 - пассивный отделяемый блок; 4 - шток; 5 - активный отделяемый блок. 6
- 8. математическая модель про- цесса отделения
- 9. схема механической системы КАОРБИТАЛЬНАЯ СТУПЕНЬ ПЛОСКОСТИ СТЫКА КА И ОРБИТАЛЬНОЙ СТУПЕНИ ЕДИНИЧНЫЙ ВЕКТОР НОРМАЛИ ПЛОСКОСТИ СТЫКА 8
- 10. допущения K K0 • КА и орбитальная ступень – абсолютно твёрдые тела. • Сумма внешних сил, действующих на систему равна нулю. • Сила i-го толкателя пропорциональная удлинению ∆i. • Направление силы действия i-го толкателя на КА и орбитальную ступень определяется единичным вектором ⃗ai. 9
- 11. исходные данные K K0 • Массы и тензоры инерции КА и орбитальной ступени, положения центров масс тел. • Линия действия силы i-го толкателя определяется единичным вектором ⃗ai: a (1) i = (a (1) ix , a (1) iy , a (1) iz )T = const. 10
- 12. исходные данные K K0 • Положение i-го толкателя на орбитальной ступени определяется вектором ⃗ρ1i: ρ (1) 1i = (ρ (1) 1ix, ρ (1) 1iy, ρ (1) 1iz)T = const. • Положение плоскости стыка КА определяется векторами ⃗np, ⃗ρp: ρ (2) p = const, n (2) p = const. 10
- 14. уравнение линии действия толкателя K K0 • Уравнение прямой (линии действия толкателя) относительно полюса C1 ⃗ρC1K = ⃗ρ1i + ⃗ai∆i • относительно полюса O ⃗ρ0K = ⃗r1 + ⃗ρ1i + ⃗ai∆i • относительно полюса C2 ⃗ρC2K = ⃗ρ2i = ⃗r1 −⃗r2 + ⃗ρ1i + ⃗ai∆i 12
- 15. уравнение плоскости поперечного стыка K K0 • Условие принадлежности точки контакта плоскости поперечного стыка (⃗ρ2i − ⃗ρp) · ⃗np = 0 • Точка контакта K должна принадлежать и плоскости и линии действия силы толкателя: { ⃗ρ2i = ⃗r1 −⃗r2 + ⃗ρ1i + ⃗ai∆i (⃗ρ2i − ⃗ρp) · ⃗np = 0 13
- 16. определение удлинения толкателя { ⃗ρ2i = ⃗r1 −⃗r2 + ⃗ρ1i + ⃗ai∆i (⃗ρ2i − ⃗ρp) · ⃗np = 0 Подставляя ⃗ρ2i из первого уравнения во второе, получим (⃗r1 −⃗r2 + ⃗ρ1i + ⃗ai∆i − ⃗ρp) · ⃗np = 0. Векторная форма выражения для определения перемещения штока толкателя: ∆i = (⃗r2 −⃗r1 − ⃗ρ1i + ⃗ρp) · ⃗np ⃗ai · ⃗np Радиус-вектор точки контакта относительно полюса C2: ⃗ρC2K = ⃗ρ2i = ⃗r1 + ⃗ai∆i −⃗r2 14
- 17. координатная форма Координатные столбцы векторов, определяющих положение и направление толкателей, заданы в базисе C1x1y1z1*: ρ (1) 1i , a (1) i Координатные столбцы векторов, определяющих положение и ориентацию плоскости стыка КА, заданы в базисе C2x2y2z2: ρ (2) p , n (2) p *Верхний индекс в скобках обозначает номер базиса, в котором записаны координатные столбцы. 15
- 18. координатная форма Координатная форма выражения для определения перемещения штока толкателя имеет вид: ∆i = n (2) p [A2T(r (0) 2 − r (0) 1 − A1ρ (1) 1i ) + ρ (2) p ] a (1)T i A1TA2n (2) p Координатный столбец радиус-вектора точки контакта относительно полюса C2 в базисе C2x2y2z2: ρ (2) 2i = A2T (r (0) 1 + A1 a (1) i ∆i − r (0) 2 ) 15
- 19. определение сил и моментов
- 20. сила пружинного толкателя Сила пружинного толкателя пропорциональна перемещению штока толкателя ∆i: Fi(∆i) = { Fi0 − Fi0−Fik hi ∆i ∆i ≤ hi 0, ∆i > hi где Fi0 – начальное усилие толкателя, Fik – конечное усилие толкателя, hi – ход толкателя. 17
- 21. сила и момент, действующие на ка Векторная форма ⃗F2i = ⃗aiFi(∆i) Координатная форма Координатный столбец силы в базисе Ox0y0z0 F (0) 2i = a (0) i Fi(∆i) = A1 a (1) i Fi(∆i) Координатный столбец силы в базисе C2x2y2z2 F (2) 2i = a (2) i Fi(∆i) = A2T A1 a (1) i Fi(∆i) 18
- 22. сила и момент, действующие на ка Векторная форма ⃗M2i(⃗F2i) = ⃗ρ2i × ⃗aiFi(∆i) Координатная форма Координатный столбец момента силы в базисе C2x2y2z2 M (2) 2i (F (2) 2i ) = ˜ρ (2) 2i a (2) i Fi(∆i) 19
- 23. матричная запись векторного произведения ˜ρ (2) 2i – кососимметричная матрица компонент координатного столбца вектора ρ (2) 2i : ˜ρ (2) 2i = 0 −ρ (2) 2iz ρ (2) 2iy ρ (2) 2iz 0 −ρ (2) 2ix −ρ (2) 2iy ρ (2) 2ix 0 Здесь и далее оператор “тильда” преобразует координатный столбец v вектора ⃗v в кососимметричную матрицу v = vx vy vz , ˜v = 0 −vz vy vz 0 −vx −vy vx 0 20
- 24. сила и момент, действующие на орбитальную ступень Векторная форма ⃗F1i = −⃗F1i = −⃗aiFi(∆i) Координатная форма Координатный столбец силы в базисе Ox0y0z0 F (0) 1i = −F (0) 2i = −a (0) i Fi(∆i) = −A1 a (1) i Fi(∆i) Координатный столбец силы в базисе C1x1y1z1 F (2) 1i = −a (1) i Fi(∆i) 21
- 25. сила и момент, действующие на орбитальную ступень Векторная форма ⃗M1i(⃗F1i) = −(⃗ρ1i + ∆i⃗ai) × ⃗aiFi(∆i) = −⃗ρ1i × ⃗aiFi(∆i) Координатная форма Координатный столбец момента силы в базисе C1x1y1z1 M (1) 1i (⃗F1i) = −˜ρ (1) 1i a (1) i Fi(∆i) ˜ρ (1) 1i – кососимметричная матрица компонент координатного столбца вектора ρ (1) 1i 22
- 27. уравнения движения центра масс ка Координатная форма m2¨x2 = ∑np i=1 F1ix, m2¨y2 = ∑np i=1 F1iy, m2¨z2 = ∑np i=1 F1iz Матричная форма m2¨r (0) 2 = np ∑ i=1 F (0) 1i 24
- 28. уравнения движения ка вокруг центра масс Координатная форма J2x ˙ω2x − (J2y − J2z)ω2yω2z = ∑np i=1 M (2) 1ix, J2y ˙ω2y − (J2z − J2x)ω2zω2x = ∑np i=1 M (2) 1ix, J2z ˙ω2z − (J2x − J2z)ω2xω2y = ∑np i=1 M (2) 1ix Матричная форма J (2) 2 · ˙ω (2) 2 + ˜ω (2) 2 J (2) 2 · ω (2) 2 = np ∑ i=1 M (2) 2i где J (2) 2 - тензор инерции КА в базисе C2x2y2z2. 25
- 29. уравнения движения ка вокруг центра масс • Для определения матрицы поворота динамические уравнения Эйлера необходимо дополнить кинематическими уравнениями, связывающими производные параметров, задающих ориентацию тела, с угловыми скоростями. • Кинематические уравнения для углов Крылова (Брайнта) ˙α21 = cos α23 cos α22 ω2x − sin α23 cos α22 ω2y, ˙α22 = sin α13ω2x + cos α13ω2y, ˙α23 = − cos α23 tan α22ω2x + sin α23 tan α22ω2y + ω2z. где ω2x, ω2y, ω2z – проекции угловой скорости КА на оси системы координат C2x2y2z2. 26
- 30. уравнения движения ка вокруг центра масс Ориентация базиса C2x2y2z2 относительно базиса Ox0y0z0 задаётся тремя последовательными поворотами на углы α21, α22 и α23. Матрица поворота: A1 = c2c3 −c2s3 s2 c1s3 + s1s2c3 c1c3 − s1s2s3 −s1c2 s1s3 − c1s2c3 s1c3 + c1s2s3 c1c2 где cj = cos α2j, sj = sin α2j, j = 1, 2, 3. 27
- 31. уравнения движения орбитальной ступени • Уравнение движения центра масс m1¨r (0) 1 = np ∑ i=1 F (0) 1i • Уравнение движения вокруг центра масс J1 · ˙ω (1) 1 + ˜ω (1) 1 J1 · ω (1) 1 = np ∑ i=1 M (1) 1i 28
- 32. уравнения движения орбитальной ступени • Ориентация базиса C1x1y1z1, связанного с орбитальной ступенью, относительно базиса Ox0y0z0 задаётся тремя последовательными поворотами на углы α11, α12 и α13. • Кинематические уравнения для углов α11, α12 и α13: ˙α11 = cos α13 cos α12 ω1x − sin α13 cos α12 ω1y, ˙α12 = sin α13ωx2 + cos α13ωy2 , ˙α13 = − cos α13 tan α12ω1x + sin α13 tan α12ω1y + ω1z. где ω1x, ω1y, ω1z – проекции угловой скорости КА 29
- 33. список использованных источников 1. Расчёт и проектирование систем разделения ступеней ракет: Учеб. пособие / К. С. Колесников, В. В. Кокушкин, С. В. Борзых, Н. В. Панкова. – Изд-во МГТУ им. Н. Э. Баумана, 2006. 2. Аналитическое проектирование механических систем: Учеб. пособие / Г. Е. Круглов. – Самараск. гос. аэрокосм. ун-т, Самара, 2001. 3. Й. Виттенбург Динамика систем твёрдых тел. М.: Мир, 1980. 30