面元法计算超音速机身的气动参数的编程实验
机身参数如下图所示：
面元法计算机身的气动参数的matlab程序代码（仅供参考）
function [] = surface(a,b,D,n)
P=0;
Q=0;
R=0;
L=D*4.2;
A=16.5*pi/180;B=20*pi/180;
S=pi*(L/4*tan(A)+D/2)^2
for i=1:11
X(i)=-a*(i-1)/10;
r1=sqrt(1-(X(i)+a)*(X(i)+a)/(a*a))*b;
for j=1:(n+1)
Y(i,j)=r1*cos(2*pi*(j-1)/n);
Z(i,j)=r1*sin(2*pi*(j-1)/n);
end
end
for i=12:21
X(i)=-a-(i-11)*(L/4-a)/10;
r2=-(X(i)+a)*tan(B)+b;
for j=1:(n+1)
Y(i,j)=r2*cos(2*pi*(j-1)/n);
Z(i,j)=r2*sin(2*pi*(j-1)/n);
end
end
for i=22:31
X(i)=-L/4-(i-21)*(L/2)/10;
r3=D/2;
for j=1:(n+1)
Y(i,j)=r3*cos(2*pi*(j-1)/n);
资源评论
