KSR-III로켓 Gravity Turn

소개글

"KSR-III로켓 Gravity Turn"에 대한 내용입니다.

목차

1. 개요

2. 이론정리
가. KSR-III 로켓과 유사한 로켓의 제원
나. The Rocket Equation

3. MATLAB 프로그램 사용
가. 고도에 따른 밀도
나. 속도에 따른 항력계수

4. 결론
가. 초기 발사각 70°
나. 초기 발사각 75°
다. 초기 발사각 80°

5. 고찰

6. 참고문헌

본문내용

1.개요
표준대기표와 KSR-III로켓의 자료를 통해 초기 발사각에 따른 로켓의 궤적을 분석할 수 있다.

<중 략>

%% Gravity Turn Trajectory 미분방정식
function dfdt=GTT1(t,y)
global m0 g0 T A Re hm mc
v=y(1); %% 속도
gr=y(2); %% 비행경로각
x=y(3); %% 지상거리
h=y(4); %% 고도
if t<57
m=m0-mc*t; %% 연소종료 직전까지의 발사체 총 중량 변화
else
m=m0-mc*57;T=0;A=0.8; %% 1단 연소 종료 후 발사체 총 중량 및 발사체 면적
end

g=g0/(1+h/Re)^2; %%고도에 따른 중력가속도 변화
rh=density(h); %% 고도에 따른 밀도변화
Cd=cd(v,h); %%속도에 따른 항력계수 변화
D=1/2*rh*v^2*A*Cd; %% 항력

if h <= hm %% 미분방적식(속도,비행경로각,지상거리,고도)
dv_dt=T/m-D/m-g;
dgr_dt=0;
dx_dt=0
dh_dt=v;
else
dv_dt=T/m-D/m-g*sin(gr);
dgr_dt=-1/v*(g-v^2/(Re+h))*cos(gr);
dx_dt=v*cos(gr);
dh_dt=v*sin(gr);
end
dfdt=[ dv_dt,dgr_dt, dx_dt,dh_dt]';

참고 자료

<표-3 표준 대기표,고도에 따른 밀도와 음속의 변화> 출처-유체역학 교재 Frank M.White