Народный космос-Z

Объявление

ВНИМАНИЕ.

Гости могут оставлять сообщения только на форуме "Поговорим".

Экскосмиты могут оставлять сообщения только в "РУГАТЬСЯ - СЮДА" и на собственном форуме, если таковой имеется.

ПОРТАЛ.

Освоение космоса.

Форум Космопорта--astrotema--Космодром "Плесецк"--"НОВОСТИ КОСМОНАВТИКИ"--ОБЩЕПЛАНЕТАРНЫЙ ФОРУМ им. Королёва.Академия космоса и космонавтики--Owalon--Путешествие в космос--Сайт "Знания-сила "

Астрономия.

Астрофорум--Y-Astro

Информация о пользователе

Привет, Гость! Войдите или зарегистрируйтесь.


Вы здесь » Народный космос-Z » Освоение околоземного космоса » Личный взгляд на создание РН для КА класса "Космический велосипед"


Личный взгляд на создание РН для КА класса "Космический велосипед"

Сообщений 61 страница 80 из 83

61

NikkyButt написал(а):

Что касается полученных по формулам результатов - поинтересуйтесь удельным импульсом современных двигателей и прикиньте - чья формула дает более похожие на правду результаты.

По памяти для высотника сегодня максимальный что-то чуть меньше 520 сек.
Мне известно, что температура в камере сгорания у современных РД достигает 4500 цельсиев. Если принять молекулярную массу продуктов сгорания 0.022 кг, а она обязана быть больше чем у водяного пара (0.018 кг), то по моей формуле получамс 5155 м/с. :)

62

По памяти для высотника сегодня максимальный что-то чуть меньше 520 сек.
Мне известно, что температура в камере сгорания у современных РД достигает 4500 цельсиев. Если принять молекулярную массу продуктов сгорания 0.022 кг, а она обязана быть больше чем у водяного пара (0.018 кг), то по моей формуле получамс 5155 м/с.

Вообще-то, УИ даже для пустотных двигателей как правило не дотягивает до 500. У РД-0120 (двигатель блока Ц "Энергии", расчитаный в основном на работу в космосе) 455сек.
У РД-180 (который мы продаём американцам) -атмосферник (что для нас более актуально) всего 337 сек.

:P

63

Suzen написал(а):

Нет, с контроллерами я к сожалению дела не имела и я программист любитель, пишу в основном разные забавные програмки для себя и для друзей, и иногда несложные игрушки для бабушки, а работаю я по другой специальности - гейм мастером.

Может быть тогда для дела освоения космоса оформите эти коды, чтобы на экране смотрелось. А то мне надо и туды и сюды. :D Руки не доходят. Правда сопротивление атмосферы так руки и не дошли сделать по Маху. Так что вы это место в исправленных кодах как-нить выделите.
Если вам, конечно, не трудно заморочиться этой программкой. :)
// Traektr0View.cpp : implementation of the CTraektr0View class
//

#include "stdafx.h"
#include "Traektr0.h"

#include "Traektr0Doc.h"
#include "Traektr0View.h"

#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif

#include "cmath"
/////////////////////////////////////////////////////////////////////////////
double gm0;
double a;
double R;//Радиус орбиты
double R0;//Радиус Земли
double vy;
double vx;
double dvy;
double dvx;
double v;//скорость ракеты
double alfa;//угол полета
double dl;
double dt;
double dR;
double dbeta;
double beta;//угловой путь ракеты
double x;
double y;
double c;
double pi;

int i0a;
int j0a;
int fl;
char znm2[20]={0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};

double mn;//масса ракеты
double m1;//масса газов, выбрасываемых за секунду
double v1;//удельный импульс
double kx1;
double kx2;
double gama;//угол положения ракеты
double vn;
double v2;
double nu;//вязкость воздуха
double r;//лобовой радиус ракеты
double vv;
double av;
// CTraektr0View

IMPLEMENT_DYNCREATE(CTraektr0View, CView)

BEGIN_MESSAGE_MAP(CTraektr0View, CView)
//{{AFX_MSG_MAP(CTraektr0View)
    // NOTE - the ClassWizard will add and remove mapping macros here.
    //    DO NOT EDIT what you see in these blocks of generated code!
//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CTraektr0View construction/destruction

CTraektr0View::CTraektr0View()
{
// TODO: add construction code here
R=6376000;R0=6371000;
c=6371000/150;
dvx=0;dvy=0;vy=0;dt=0.01;vx=232.8;
v1=2820;//2650;
fl=100;
r=3.2;
alfa=0;pi=3.141592;
gm0=6.672*5.978*10000000000000;
v=232.8;//0;//sqrt(pow(vy,2)+pow(vx,2));
gama=0;//pi/6;
}

CTraektr0View::~CTraektr0View()
{
}

BOOL CTraektr0View::PreCreateWindow(CREATESTRUCT& cs)
{
// TODO: Modify the Window class or styles here by modifying
//  the CREATESTRUCT cs

return CView::PreCreateWindow(cs);
}

/////////////////////////////////////////////////////////////////////////////
// CTraektr0View drawing

void CTraektr0View::OnDraw(CDC* pDC)
{
CTraektr0Doc* pDoc = GetDocument();
ASSERT_VALID(pDoc);
// TODO: add draw code for native data here
CClientDC ClientDC(this);
//if(fl==10)goto a1;
a=gm0/pow(R,2);
if(alfa<pi)alfa=acos(vx/v);
else asin(vy/v);
dl=vx*dt;
dR=vy*dt;
R=sqrt(pow((R+dR),2)+pow(dl,2));
dbeta=asin(dl/R);

v_ot_m();
vy=vy-a*dt;
vx=vx+dvx;vy=vy+dvy;
v=sqrt(pow(vy,2)+pow(vx,2));

vx=vx*cos(dbeta)-vy*sin(dbeta);
vy=vy*cos(dbeta)+vx*sin(dbeta);

beta=beta+dbeta;
x=R*sin(beta)/c;
y=R*cos(beta)/c;
pDC->SetPixel(250+(int)x,200-(int)y,19);
x=R*sin(beta)*20000/6371000;
y=R*cos(beta)*20000/6371000;
pDC->SetPixel(10+(int)x,20520-(int)y,19);

ClientDC.TextOut(430,0,"Скорость:                                    ");
j0a=0;
sprintf(znm2,"%f",v);
for(i0a=0;i0a<20;i0a++){
ClientDC.TextOut(i0a*12+500,j0a*15,znm2[i0a]);
}

ClientDC.TextOut(440,15,"Радиус:                                    ");
j0a=1;
sprintf(znm2,"%f",R);
for(i0a=0;i0a<20;i0a++){
ClientDC.TextOut(i0a*12+500,j0a*15,znm2[i0a]);
}
ClientDC.TextOut(440,45,"Высота:                                    ");
j0a=3;
x=R-R0;
sprintf(znm2,"%f",x);
for(i0a=0;i0a<20;i0a++){
ClientDC.TextOut(i0a*12+500,j0a*15,znm2[i0a]);
}
    RECT Rect;
    Rect.left=0;//i00*10;
    Rect.top=0;
    Rect.right=Rect.left+500;
    Rect.bottom=Rect.top+600;
    InvalidateRect(&Rect,FALSE);
//fl=10;
//a1:;
}

/////////////////////////////////////////////////////////////////////////////
// CTraektr0View diagnostics

#ifdef _DEBUG
void CTraektr0View::AssertValid() const
{
CView::AssertValid();
}

void CTraektr0View::Dump(CDumpContext& dc) const
{
CView::Dump(dc);
}

CTraektr0Doc* CTraektr0View::GetDocument() // non-debug version is inline
{
ASSERT(m_pDocument->IsKindOf(RUNTIME_CLASS(CTraektr0Doc)));
return (CTraektr0Doc*)m_pDocument;
}
#endif //_DEBUG

/////////////////////////////////////////////////////////////////////////////
// CTraektr0View message handlers
void CTraektr0View::v_ot_m()
{
CClientDC ClientDC(this);
if(fl==100){mn=26570;m1=144;fl=0;kx10();}
if(fl==0){if(mn<=14906){mn=mn-696;m1=144;fl=1;kx10();}}
//отстрел первой ступени
                if(fl==1){if(mn<=7730){v1=3200;m1=18;mn=mn-1164;fl=2;kx10();}}
//отстрел второй, 1164 - это масса шести пустых блоков

//третьей
if(fl==2){if(mn<=2930){mn=mn-180;fl=3;kx10();}}

/////////регулировка угла положения ракеты
//if(v>=464)gama=0.5+log(((v)/4715)+1);//3.8;//4152
if(v>=232)gama=0.275+((v+330)*pi/(10000*2));//0.27//0.005+log(((v)/1715)+1);//3.8;//4152
if(gama>=(pi/2))gama=pi/2;
////////
if(fl<=3){
    if(mn>932.1){
    vn=kx1/mn;
    mn=mn-kx2;
    nu=1.72/pow(2.718,(0.000125*(R-R0)));
    v2=6*pi*nu*r*(v-232.8)*dt/mn;
//СОПРТИВЛЕНИЕ АТМОСФЕРЫ
    vn=vn-v2;
/////
    if(v>2490)//2804)//2854)//2740)//2750)
    {
        gama=pi/2;
        av=gm0/pow(R,2)-pow(vx,2)/R;
        vv=av*dt;
        if(vn>2*vv){if(vy<1.8*vv)gama=acos(1.57*vv/vn);}
    }
/////
    dvx=vn*sin(gama);dvy=vn*cos(gama);
    }
    else {dvx=0;dvy=0;}
}
ClientDC.TextOut(440,30,"Масса:                                    ");
j0a=2;
memset(znm2,32,sizeof(znm2));
sprintf(znm2,"%f",mn);
for(i0a=0;i0a<20;i0a++){
ClientDC.TextOut(i0a*12+500,j0a*15,znm2[i0a]);}
}

void CTraektr0View::kx10()
{
kx1=m1*v1*dt;kx2=m1*dt;
}

64

NikkyButt написал(а):

Вообще-то, УИ даже для пустотных двигателей как правило не дотягивает до 500. У РД-0120 (двигатель блока Ц "Энергии", расчитаный в основном на работу в космосе) 455сек.

Значит у меня в голове засел какой-то планируемый РД. :) Которого еще нет.

65

Учтите ка народ, что к нам шпиона запросто подбросить могут, а то еще и хуже, поскольку есть такие силы, которые за ракетчиками ведут охоту, так что лучше язык держать ведь за зубами и не светиться собственными фотками и прочими делами...
А девочка умна уж очень, подбросила не фотку, а блесну...

66

серж я валяюсь под столом....ГыГы :lol:

67

ничего когда чеченцы схватят за одно место и заставят делать снаряды ракетные, будешь сидеть в такой яме!!!
Как кавказкий пленник.
Мальчик, ты реальной жизни ни черта не знаешь...
И лучше и не знать...

68

NikkyButt написал(а):

[Roger Wilco уже писал формулу для скорости истечения продуктов сгорания из сопла, приведу ее в более понятном виде (см. рисунок в конце сообщения).

Эта формула точно определяет скорость потока газа на выходе из сопла в зависимости от параметров в КС и давления на выходе (известное) для идеального случая.

Если R в этой формуле - молярная газовая постоянная - то "ваша" формула - чистейший бред, потому что дает скорость истечения при 4500 цельсиев 501,3 м/с. :lol:
Как енто понимать? :mrgreen:

69

Если R в этой формуле - молярная газовая постоянная - то "ваша" формула - чистейший бред, потому что дает скорость истечения при 4500 цельсиев 501,3 м/с. 
Как енто понимать?

Ну, давайте разбираться по порядку.

Во-первых, формула совсем не моя, а заимствованная из общей теории ракетных двигателей.

Во-вторых, R - это удельная газовая постоянная, R=R0/мю, где R0 - универсальная газовая постоянная, R0=8314 [кДж/(кмоль К)]; мю - молярная масса.

В-третьих, каким образом у вас вышло 501,3 м/с - ума не приложу. Отпишите, если не трудно, параметры, что вы задавали, и будем разбираться.

70

R=R0/mu

R - удельная газовая постоянная
R0 - универсальная газовая постоянная = 8314
mu - молекулярная масса

Для
R=350
n=k=1.2
Tk=4500K
Pk=5МПа
Pa=0.05МПа
Скорость получается 3182м/с.

Попробуйте ещё раз. Если не получится - воспользуйтесь калькулятором. Говорят - помогает :)

Отредактировано Roger Wilco (2006-06-15 17:10:51)

71

Roger Wilco написал(а):

R=R0/mu

R - удельная газовая постоянная
R0 - универсальная газовая постоянная = 8314
mu - молекулярная масса

О, нас уже двое.  :)

72

Roger Wilco написал(а):

R0 - универсальная газовая постоянная = 8314

В какой системе единиц живем? Я ващето привык, что R0=8,314 дж/(К*моль) Ладно, проехали.
У меня получилось 3695. Как это согласуется с вашими 455 сек? Моя формула ближе. Завтра ее уточню. Оказывается зависимость температуры от давления интегрируется хорошо. Но уже башка не варит. :)

73

Мне непонятно одно: какой смысл изобретать велосипед (к названию темы фраза не относится), если еще задолго до нашего рождения были выведены теории расчета и проектирования, десятилетиями проверенные на практике?

Похвально, что вы пытаетесь все сделать сами "с нуля", но не проще использовать опыт поколений и сосредоточить усилия на более важных вопросах (таких как расчет ракеты и ее элементов по опробованным методикам, поиск средств, станков, заинтересованных людей)?

74

По вашей формуле:
Tk=4500K
Pk=5МПа
Pa=0.05МПа
ro0=1кг/м3 (условно - плотность воздуха)

v^2=v0^2+(2*T*P0*ln(P/P0))/(ro0*T0)
v=sqrt(v0^2+((2*T*P0*ln(P/P0))/(ro0*To)))
v0 - скорость на входе в сопло (условно 100м/с)
T   - темпер-ра в КС (4500)
P0 - давление на срезе сопла (0.05МПа)
P   - давление в КС (5МПа)
T0 - температура на срезе сопла (2000К)

Получается, около 1023 м/с
Э?
Укажите, где я не прав.

Я ващето привык, что R0=8,314 дж/(К*моль)

R0=8314     Дж/(кмоль*к)
R=8314/mu Дж/(кг*к)

Отредактировано Roger Wilco (2006-06-15 18:27:49)

75

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

76

серж написал(а):

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

Интересно узнать, а что за конкуренты, может ссылочку дадите?

77

Берт Рутан,
проект Фалкон, и многие другие, я ведь предлагал с ними договариваться и совместное предприятие открывать, но никто и не чешется, а занимаются школьными игрушками...

78

серж написал(а):

Берт Рутан,
проект Фалкон, и многие другие, я ведь предлагал с ними договариваться и совместное предприятие открывать, но никто и не чешется, а занимаются школьными игрушками...

Ну, в любой ситуации к верному решению можно прийти только посредством дискуссий и споров, так что это отнюдь не "школьные игрушки".

А производство своей ракеты начнется, думаю, не раньше чем через 10 лет. Можете звать меня пессимистом.

Берт Рутан все-таки обитает в США, где отношение к гражданскому ракетостроению не такое недоверчивое как в нашей стране. Да и со средствами у них получше будет.

Отредактировано NikkyButt (2006-06-16 19:22:03)

79

Roger Wilco написал(а):

По вашей формуле:
Tk=4500K
Pk=5МПа
Pa=0.05МПа
ro0=1кг/м3 (условно - плотность воздуха)

v^2=v0^2+(2*T*P0*ln(P/P0))/(ro0*T0)
v=sqrt(v0^2+((2*T*P0*ln(P/P0))/(ro0*To)))
v0 - скорость на входе в сопло (условно 100м/с)
T   - темпер-ра в КС (4500)
P0 - давление на срезе сопла (0.05МПа)
P   - давление в КС (5МПа)
T0 - температура на срезе сопла (2000К)

Получается, около 1023 м/с
Э?
Укажите, где я не прав.

Неправ я. P0 = 100000 Па- нормальное атмосферное давление. А в логарифме давление на срезе нужно обозначить иначе. P_c например. Ро нулевое для водяного пара равно молярная масса 0.018 кг разделить на 0.0224 метра в кубе. Я чтобы не путаться все считаю в одной системе единиц. T0=273,15 кельвинов.
Я настолько привык использовать эту формулу для давления на срезе в 1 атм., что невольно ... Извиняюсь
Формула строго выглядит так
v=sqrt(v0^2+((2*T*P0*ln(P/P_с))/(ro0*T0)))
T- температура в сопле. В данном случае я ее считаю постоянной
P0=100000 Па
P- давление на входе в сопло
P-c- давление на срезе
ro0- молекулярная масса в килограммах деленная на 0.0224 метра в кубе
T0- 273.15 кельвина
Извиняюсь еще раз. Я думал, что все эти величины очевидны.
NikkyButt

А производство своей ракеты начнется, думаю, не раньше чем через 10 лет. Можете звать меня пессимистом.

Вы оптимист.
Я долго атаковал высокие инстанции. Потом со мной произошли такие страшные чудеса, что вам о них лучше не знать. Так что сидите в своих конторах и не рыпайтесь, пока ветры наверху не подуют на нас. А до тех пор, подготовим все детали, чтоб осталось только собрать. Чтобы проверить формулы не обязательно делать полномасштабный двигатель. Если формула будет давать надежные результаты для модельки РД тягой пять килограмм и температурой горения 1000 Цельсиев, то и для восьмитонника она даст надежные результаты. :)

80

Roger Wilco написал(а):

Пусть будет технотемка для непродвинутых пользователей, не разделяющих ваши перспективные взгляды на методику проектирования РД.

Я, конечно, надеюсь, что у вас всё получится, но, если вдруг каким-то чудом проверенная 70ю годами практики классическая теория окажется, ткскзть, правее, - неплохо бы иметь запасной проект. Чтобы времени не терять.

А я не против классической методики проектирования РД. Если вы заглядываете в новостную ленту космонавтики, то должны знать, что после установки новых РД семерка сможет выводить на орбиту до 16 тонн. Расход горючки на космонавта снизится до 50 тонн. А вы не задумывались, почему бы не сделать полностью транспортный двухтонный корабль, вместо того, чтобы всякий раз таскать на Союзе 4 тонны ненужного орбитального отсека?
А потому что никто не хочет заморачиваться на создание маленькой ракеты, в которой с новыми двиглами как раз и тратилось бы 25тонн горючей пары, как у нас. Ну ладно, за комфортный полет на Клипере можно потратиться и на 50 тонн горючей пары.
Кстати, семерку можно легко сделать полностью многоразовым ракетоносителем. Нужно всего лишь отказаться от стереотипов мышления. :D


Вы здесь » Народный космос-Z » Освоение околоземного космоса » Личный взгляд на создание РН для КА класса "Космический велосипед"