Я запитав це питання на Aviation.stackexchange, але після цього я подумав, що, можливо, краще розмістити його тут. Тим більше, що надійшло хороших відповідей на покоління NACA з 5 цифр серії .
Я хотів би обчислити профіль NACA 64-2A015. За допомогою Aviation.stackexchange я дізнався, що A-Версія профілю створена для полегшення виготовлення шляхом потовщення заднього краю перерізу (прямим контуром від 80% хорди назад).
Однак я намагаюся побудувати профіль на основі рівнянь, наведених тут і тут .
Реалізація в MATLAB виглядає приблизно так: Ось код Matlab, який я використав:
a = 0.4;
b = 1.0; % caution for NON-unity entries change the equation for h
% c = 1; to simplifiy the equation the chord is set to 1
cl = 1;
g = -1/(b-a) * (a^2 * (1/2 * log(a) -1/4) - b^2 * (1/2 * log(b) -1/4)); % g = -1/(1-a) * (a^2 * (1/2 * log(a) -1/4) + 1/4)
h = 1/(1-a) * (1/2*(1-a)^2 * log(1-a) -1/4*(1-a)^2)+g; % simplified version for b = 1: h = 1/(b-a) * (1/2*(1-a)^2 * log(1-a) - 1/2 * (1-b)^2 * log(1-b) + 1/4*(1-b)^2 - 1/4*(1-a)^2) + g
x = (0:0.001:1);
y = cl/(2*pi*(a+b)) * ( 1/(b-a) .* (1/2 * (a-x).^2 .* log(abs(a-x)) - 1/2 .* (b-x).^2 .* log(abs(b-x)) + 1/4 .* (b-x).^2 - 1/4 .* (a-x).^2) - x.*log(x) + g - h.*x); % y = cl/(2*pi*(a+1)) * ( 1/(1-a) .* (1/2 * (a-x).^2 .* log(abs(a-x)) - 1/2 .* (1-x).^2 .* log(abs(1-x)) + 1/4 .* (1-x).^2 + 1/4 .* (a-x).^2) - x.*log(x) + g - h.*x);
L6j01_x = [0, 0.005, 0.0075, 0.0125, 0.025, 0.05, 0.075, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1];
L6j01_y = [0, 0.01193, 0.01436, 0.01815, 0.02508, 0.03477, 0.04202, 0.04799, 0.05732, 0.06423, 0.06926, 0.0727, 0.07463, 0.07487, 0.07313, 0.06978, 0.06517, 0.05956, 0.05311, 0.046, 0.03847, 0.03084, 0.02321, 0.01558, 0.00795, 0.00032];
plot(x,y), axis equal, hold on
plot(L6j01_x, L6j01_y,'ro'), hold off
Але результати не відповідають точкам даних з Nasa-Reports
.
Що я пропустив, щоб реалізувати в рівнянні, чи є інше зображення крила (за винятком радіуса переднього та заднього краю)?
