如何定义项目5位数系列方程?
10视图(30天)
显示旧的评论
我一直在试图分析定义的形状NACA系列5位数(包括反射和non-reflex)。然而,当试图制造阴谋,MATLAB产生的形状非常奇怪。对于一些翼型,它产生的形状是完美的,但对其他人来说,它不是。我有图片。我还没有找到一个模式去看哪些失败,不,还有一个原因。
几何定义函数代码:
函数[x_upper, y_upper x_lower y_lower] = NACAFoil5 (NACA5和弦)%定义的函数
x = 0:0.00001:和弦;% x坐标的定义
t = str2num ([num2str (NACA5 (4)), num2str (NACA5 (5)))) / 100;% #好< ST2NM >
y_t = 5 * t *和弦* 0.2969 *√(x /和弦)-0.126 * -0.3516 (x /和弦)* (x /和弦)^ 2 + 0.2843 * (x /和弦)。^ 3 - 0.1015 * (x /和弦)^ 4);
p = NACA5 (2) / 20;
如果长度(NACA5) = = 5
disp ([项目实施5位数系列:项目实施的num2str (NACA5 (1)) num2str (NACA5 (2) num2str (NACA5 (3) num2str (NACA5 (4) num2str (NACA5 (5))))
如果NACA5 (3) = = 0
反射= 0;
elseifNACA5 (3) = = 1
反射= 1;
结束
结束
如果反射= = 0
如果NACA5 (2) = = 1
disp (5%不反折的弯度翼型的)
r1 = 0.0580;
k1 = 361.4;
elseifNACA5 (2) = = 2
disp (10%不反折的弯度翼型的)
r1 = 0.126;
k1 = 51.64;
elseifNACA5 (2) = = 3
disp (15%不反折的弯度翼型的)
r1 = 0.2025;
k1 = 15.957;
elseifNACA5 (2) = = 4
disp (20%不反折的弯度翼型的)
r1 = 0.29;
k1 = 6.643;
elseifNACA5 (2) = = 5
disp (25%不反折的弯度翼型的)
r1 = 0.391;
k1 = 3.23;
结束
为i = 1:长度(x)
如果(x(我)/和弦)< = p
y_camber (i) = (k1/6) * (((x(我)/和弦)^ 3)——(3 * r1 * (x(我)/和弦)^ 2)+ (r1 ^ 2) * (3-r1) * (x(我)/和弦));
dy_camber (i) = (k1/6) * ((3 * (x(我)/和弦)^ 2)——(6 * r1 * (x(我)/和弦))+ (r1 ^ 2) * (3-r1));
elseif(x(我)/和弦)> p
y_camber (i) = ((k1 * r1 ^ 3) / 6) * (1 - (x(我)/和弦));
dy_camber (i) = (1 * k1 * r1 ^ 3) / 6;
结束
结束
elseif反射= = 1
如果NACA5 (2) = = 2
disp (“10%反折的弯度机翼的)
r1 = 0.13;
k1 = 51.99;
k2_k1 = 0.000764;
elseifNACA5 (2) = = 3
disp (“15%反折的弯度机翼的)
r1 = 0.217;
k1 = 15.793;
k2_k1 = 0.00677;
elseifNACA5 (2) = = 4
disp (“20%反折的弯度机翼的)
r1 = 0.318;
k1 = 6.52;
k2_k1 = 0.303;
elseifNACA5 (2) = = 5
disp (“25%反折的弯度机翼的)
r1 = 0.441;
k1 = 3.191;
k2_k1 = 0.1355;
结束
为i = 1:长度(x)
如果(x(我)/和弦)< = p
y_camber (i) = (k1/6) * ((((x (i) /和弦)r1) ^ 3)——(k2_k1) * ((1-r1) ^ 3) * x (i) - (r1 ^ 3) * ((x(我)/和弦))+ r1 ^ 3);
dy_camber (i) = (k1/6) * (3 * (((x (i) /和弦)r1) ^ 2)——(k2_k1) * ((1-r1) ^ 3) r1 ^ 3);
elseif(x(我)/和弦)> p
y_camber (i) = (k1/6) * ((k2_k1) * (((x (i) /和弦)r1) ^ 3)——(k2_k1) * ((1-r1) ^ 3) * (x(我)/和弦)- (r1 ^ 3) * ((x(我)/和弦)+ r1 ^ 3));
dy_camber (i) = (k1/6) * ((3 * k2_k1) * (((x (i) /和弦)r1) ^ 2)——(k2_k1) * ((1-r1) ^ 3) r1 ^ 3);
结束
结束
结束
θ=每股(dy_camber);
x_upper = x - (y_t。* sin(θ));
y_upper = y_camber + (y_t。* cos(θ));
x_lower = x + (y_t。* sin(θ));
y_lower = y_camber - (y_t。* cos(θ));
y_upper(结束)= 0;
y_lower(结束)= 0;
结束
解决:
clc
清晰的所有
关闭所有
% %身体方程
AirfoilTypeAsk =”为Jokowski机翼输入“J”,“它们”项目实施5系列翼型,“N”NACA翼型或缸的“C”:“;%的身体类型
类型=输入(AirfoilTypeAsk);
如果类型= =“它们”
AirfoilAsk =”进入机翼在原始矢量数[x x x x x]:“;%机翼数量
机翼=输入(AirfoilAsk);
NACA5 =机翼;
ChordAsk =“输入翼弦长:”;%机翼弦
和弦=输入(ChordAsk);
[x_upper, y_upper x_lower y_lower] = NACAFoil5 (NACA5和弦);
elseif类型= =“N”
AirfoilAsk =”进入机翼在原始矢量数[x x x]:“;%机翼数量
机翼=输入(AirfoilAsk);
项目实施=机翼;
ChordAsk =“输入翼弦长:”;%机翼弦
和弦=输入(ChordAsk);
[x_upper, y_upper x_lower y_lower] = NACAFoil(项目实施,和弦);
elseif类型= =“J”
AirfoilAsk =“进入机翼max。厚度/和弦和马克斯。曲面/和弦(t / c c / c):“;%机翼厚度和曲率
机翼=输入(AirfoilAsk);
ChordAsk =“输入翼弦长:”;%机翼弦
和弦=输入(ChordAsk);
t_c =翼(1);
C_c =翼(2);
[x_upper, y_upper x_lower y_lower] = JFoil (t_c C_c,和弦);
elseif类型= =“C”
CylinderAsk =的输入半径:;
Raduis =输入(CylinderAsk);
和弦= 2 * Raduis;
[x_upper, y_upper x_lower y_lower] =缸(和弦);
结束
图%机翼
持有在
网格在
轴平等的
情节(0:0.0001:和弦,polyval (polyfit (x_upper y_upper 10), 0:0.0001:和弦),0:0.0001:和弦,polyval (polyfit (x_lower y_lower 10), 0:0.0001:和弦),“线宽”,1.5,“颜色”,“b”)
情节(0:0.0001:和弦,(polyval (polyfit (x_upper y_upper 10), 0:0.0001:和弦)+ polyval (polyfit (x_lower y_lower 10), 0:0.0001:和弦))/ 2,“——”,“线宽”,1.5,“颜色”,“r”)
如果类型= =“N”
标题([“项目”num2str(项目实施(1))num2str(项目实施(2)num2str (NACA (3)) num2str (NACA (4)),“翼”])
elseif类型= =“J”
标题([“儒科夫斯基翼型绕流流行:t_{马克斯}/ c = '识别num2str (t_c * 100)' % C_{马克斯}/ c = 'num2str (C_c * 100)“%”])
elseif类型= =“它们”
标题([“NACA5”num2str (NACA5 (1)) num2str (NACA5 (2) num2str (NACA5 (3) num2str (NACA5 (4) num2str (NACA5 (5)),“翼”])
结束
包含(“x”)
ylabel (“y”)
传奇(的上表面,的下表面,“弧线”)
数据,包括方程使用来自以下参考:
同时,流的主要代码是为了解决在一个任意的形状,我只是想添加项目5到解算器的形状。所以请忽略其他形状类型和各自的命令。
