This preview has intentionally blurred sections. Sign up to view the full version.
View Full DocumentThis preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: A2 + dzdx*cos2t*dtheta end B0 = B0/pi A1 = A1*2/pi A2 = A2*2/pi % alpha loop for i = 1:31 alphadeg(i) = i11 alpha = alphadeg(i)*pi/180 A0 = alpha  B0 Cl(i) = 2*pi*(A0 + A1/2) Cm(i) = pi/4*(A2  A1) Xcp(i) = 1/4*(1 + pi/Cl(i)*(A1  A2)) end Xac = .25 alpha0 = B0  A1/2 alpha0 = alpha0*180/pi % load experimental data data2 = load('naca4415 alpha vs Cl.txt') [m,n] = size(data2) for i = 1:m Clexp(i) = data2(i,2) end data3 = load('naca4415 alpha vs Cm.txt') [m,n] = size(data2) for i = 1:m Cmexp(i) = data3(i,2) end % graph figure(1) xlabel('Angle of Attack (degrees)') title('Moment Coefficient Data') figure(2) xlabel('Angle of Attack (degrees)') title('Lift Coefficient Data') figure(3) xlabel('Angle of Attack (degrees)') title('Location of Center of Pressure') Page 2 Sheet1 clear all clc Page 3 Sheet1 plot(alphadeg,Cm,alphadeg,Cmexp,'ro') ylabel('Cm') plot(alphadeg,Cl,alphadeg,Clexp,'ro') ylabel('Cl') plot(alphadeg,Xcp) ylabel('Xcp') Page 4...
View
Full Document
 Spring '08
 Hindman
 NJ, upper layer fxup, load experimental data, Moment Coefficient Data, Lift Coefficient Data, A1*2/pi A2

Click to edit the document details