KSR-III로켓 Gravity Turn
- 최초 등록일
- 2020.12.15
- 최종 저작일
- 2014.07
- 13페이지/
한컴오피스
- 가격 2,500원
![할인쿠폰받기](/images/v4/document/ico_det_coupon.gif)
소개글
"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