Wie zeichnet man Tragflächen der NACA 6er-Serie?

In Bezug auf die Frage zum Tragflügel NACA 64-2A015 würde ich gerne wissen, wie man diesen Tragflügel zeichnet.

Geben Sie hier die Bildbeschreibung ein

Zumindest liefern diese beiden Berichte [1 , 2] der NASA die Gleichung dafür.Geben Sie hier die Bildbeschreibung ein

Ich bin jedoch nicht in der Lage, die richtigen Parameter zu finden, um es zu zeichnen, damit es mit den in der obigen Tabelle angegebenen Koordinaten übereinstimmt.

Nacaplotter

Hier ist der Matlab-Code, den ich verwendet habe:

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

Es scheint, dass meine Leading Edge-Definition nicht richtig ist. Irgendwelche Vorschläge, wie man den vorderen Teil des Tragflächenabschnitts justiert?

Es ist schon eine Weile her, aber ich erinnere mich vage, in Vorlesungen gehört zu haben, dass die Vorderkante eines Tragflügels ein Kreis ist. Also irgendwann würde sich die Gerade aus der Gleichung in einen Kreis mit unbekanntem Radius verwandeln, vielleicht kann dir der NASA Report Server weiter helfen.
Jupp, du hast recht. Wenn ich es richtig sehe, ist dies auch in der Abbildung/Zeichnung aus dem NASA-Bericht oben angegeben. Aber ich sprach eigentlich von der Nichtübereinstimmung meiner berechneten Linie (blau) mit den roten Kreisen. Besonders der Abschnitt 1-40% Akkord.
Oh richtig, ich nahm an, dass der Startpunkt versetzt werden musste, um vom Nasenkreis statt von (0,0) abzuweichen.
Wenn Sie Matlab verwenden, kann es ein Polynom durch Ihre Datenpunkte konstruieren?

Antworten (3)

Aus dieser Referenz :

c F = K 1 T c + K 2 ( T c ) 2 + K 3 ( T c ) 3 + K 4 ( T c ) 4

Für die 64A-Serie:

  • K 1 = 8,2125018
  • K 2 = 0,7685596
  • K 3 = 1,4922345
  • K 4 = 3,6130133

c F ist der jeweilige Skalierungsfaktor für dieses Profil. Und dann geht die Referenz mit dieser kopfzerbrechenden weiteren Erklärung weiter:

Nun kann für eine spezifizierte Familie und Dicke die Dickenverteilung ohne Iteration bestimmt werden. Aus der Dicke wird der Skalierungsfaktor aus der oben gezeigten Polynomfunktion berechnet. Dann wird der Skalierungsfaktor verwendet, um die Basiswerte der Psi- und Epsilon-Funktionen für diese Tragflächenfamilie zu multiplizieren. Diese skalierten psi- und epsilon-Funktionen werden beim Abbilden der z-Ebene auf die in Fig. 1 gezeigte z'-Ebene verwendet. Die Joukowski-Funktion zeta = z' + 1/z' bildet dann die z'-Ebene auf die zeta-Ebene und ab diese Ergebnisse werden so normiert, dass die vordere Flanke bei x=0 und die hintere Flanke bei x=1 liegt.

Hoffe, das hilft, mein Kopf tut weh, nachdem ich das gelesen habe.

Die Tragflächen der 1er- und 6er-Serie waren die ersten, die eine Zieldruckverteilung verwendeten, um die Kontur der Oberfläche zu definieren. Dies war ein großer Fortschritt in der Tragflächenkonstruktion, macht es jedoch unmöglich, eine einfache Gleichung zu verwenden, um die Kontur mathematisch genau zu beschreiben. Die ältere, 4-stellige Serie verwendete einfach den Mittelwert einer Sammlung erfolgreicher früher Tragflächen .

Ich habe immer eine Sammlung von Fortran-Subroutinen verwendet, um die Kontur zu generieren, und es verwendet tabellierte Daten und Splines, um die Form zu generieren. Hier ist ein Link zu einer ähnlichen Sammlung, aber in allen Fällen erhalten Sie nur eine Annäherung. Verwenden Sie XFOIL oder einen ähnlichen Code, um die Druckverteilung wie gewünscht zu verfeinern.

Was meinen Sie damit, dass das Diagramm den Vorderkantenradius (LE) nicht angibt? Es besagt eindeutig, dass der Vorderkantenradius 1,561 % der Sehne beträgt. Nennen wir den Akkord einfach die Einheit von "1", um die Dinge einfacher zu machen. Alles, was größer oder kleiner ist, ist einfach ein Verhältnis, richtig? Daher ist der LE-Radius 1,561. Der Plot sieht symmetrisch aus, gespiegelt über die Akkordlinie. Vielleicht interpretiere ich da zu viel hinein. Unabhängig davon hoffe ich, dass dies hilft. Andere Tragflächen geben im Allgemeinen die oberen und unteren Koordinaten an, deshalb vermute ich, dass die Darstellung über die Akkordlinie symmetrisch ist.