목차
◎ 실험 개요
◎ 실험 과정
◎ 실험 고찰 및 분석
1) 추력 곡선 그래프
2) 프로그래밍을 이용한 비추력
[MATLAB]
◎ 실험 과정
◎ 실험 고찰 및 분석
1) 추력 곡선 그래프
2) 프로그래밍을 이용한 비추력
[MATLAB]
본문내용
해당하는 힘을 무차원 계수 단위당 힘으로 환산한다.
a=zeros(1,length(y1)); % 초기치 51을 0으로 만들기 위해 임의의 배열을 설정하고 값을 대입함.
a(1,:)=52; % 보정치 이외에 로켓 및 테이프 무게 포함하여 51에서 52.
Y1=(y1-a)*N; % 그래프를 그리기 위해 초기치를 모두 0으로 설정하도록 하며 무차원 계수를 힘으로 변환함.
Y2=(y2-a)*N;
% 그래프를 두개 같이 프린팅함.
plot(t,Y1,'r--',t,Y2,'b-.');
grid;
legend('First Test','Second Test');
title('B6-4 Type Rocket Engine Thrust Test');
xlabel('Time[sec]with unit interval 25ms'); ylabel('Measured Force[N]');
% 비추력 구하기.
%수치적분 기법 복합 사다리꼴 법칙[composite trapezoidal rule]을 이용
It1=0;
It2=0;
k=1.5/0.025; % 사출 지연 시간을 제외한 추력발생 시간/time interval
for i=1:k
It1=It1+(0.5*h)*(Y1(1,i)+Y1(1,i+1)); %total impulse 계산
It2=It2+(0.5*h)*(Y2(1,i)+Y2(1,i+1));
end
T1a=It1/1; %평균 추력을 구함.
T2a=It2/1; %평균 추력을 구함.
mp=0.02; % 연소한 추진제의 총 질량
Is1=It1/(mp*g); % 비추력1을 구함
Is2=It2/(mp*g); % 비추력2을 구함
Is1
Is2
a=zeros(1,length(y1)); % 초기치 51을 0으로 만들기 위해 임의의 배열을 설정하고 값을 대입함.
a(1,:)=52; % 보정치 이외에 로켓 및 테이프 무게 포함하여 51에서 52.
Y1=(y1-a)*N; % 그래프를 그리기 위해 초기치를 모두 0으로 설정하도록 하며 무차원 계수를 힘으로 변환함.
Y2=(y2-a)*N;
% 그래프를 두개 같이 프린팅함.
plot(t,Y1,'r--',t,Y2,'b-.');
grid;
legend('First Test','Second Test');
title('B6-4 Type Rocket Engine Thrust Test');
xlabel('Time[sec]with unit interval 25ms'); ylabel('Measured Force[N]');
% 비추력 구하기.
%수치적분 기법 복합 사다리꼴 법칙[composite trapezoidal rule]을 이용
It1=0;
It2=0;
k=1.5/0.025; % 사출 지연 시간을 제외한 추력발생 시간/time interval
for i=1:k
It1=It1+(0.5*h)*(Y1(1,i)+Y1(1,i+1)); %total impulse 계산
It2=It2+(0.5*h)*(Y2(1,i)+Y2(1,i+1));
end
T1a=It1/1; %평균 추력을 구함.
T2a=It2/1; %평균 추력을 구함.
mp=0.02; % 연소한 추진제의 총 질량
Is1=It1/(mp*g); % 비추력1을 구함
Is2=It2/(mp*g); % 비추력2을 구함
Is1
Is2
소개글