[기체동역학]Gas Dynamics - 경사충격파 식 유도 및 메틀랩 파일

소개글

기체동역학 (Gas Dynamics) 과제구요
경사충격파 식(쎄타-델타-M)유도 하였구요
손으로 유도해서 스캔 떴습니다.
그 결과를 메틀랩으로 그래프로 나타내었습니다.
메틀랩 소스 파일도 함께 첨부하니 참고하시구요. 수정해서 쓰세요~
매틀랩 6.5로 작업했습니다.

목차

1. 경사충격파(Oblique shock wave) 식 유도
2. matlab 소스
3. 그래프

본문내용

Ⅱ. matlab m file
clc; clear all;

r=1.4;
m0=[100:20:400 450 500 600 800 1000 2000];
for i=1:1:22
m=m0(i)/100;
j=1;
for b=0.1:0.1:90
a(i,j)=(atan(2*cot(b*pi/180)*(m^2*(sin(b*pi/180))^2-1)/(m^2*(r+cos(2*b*pi/180))+2)))*180/pi;
j=j+1;
end
end

b=0.1:0.1:90;
plot(a,b);

AXIS([0 50 0 90])
grid on;
title(`Oblique-shock wave (\gamma = 1.4)`)
xlabel(`Deflection angle(degree) \delta`)
ylabel(`Shock-wave angle(degree) \theta`)
text(43.88, 74.69,`\rightarrowM=1`);
text(28.61, 67.16,`\rightarrowM=4`);
text(19.45, 52.13,`\rightarrowM=5`);
text(8.05, 48.12,`\rightarrowM=8`);
text(2.44, 48.51,`\rightarrowM=10`);
text(0.39, 56.91,`\rightarrowM1=20`);
text(38.35, 83.61, `strong shocks`);
text(17.08, 5.18, `weak shocks`);

참고 자료

없음