Lieber Jondel, deine 30° interessieren niemanden. Du kannst es dir aber selber ausrechnen, wenn du deine Zahlen in meinen Code eingibst, mir ist meine Zeit dafür zu schade. Zudem hast du keinen ballistischen Koeffizienten oder Cw-Wert angegeben. Ich illustriere derweil lieber meine eigene Aussage:
Dateianhang:
Dateikommentar: t = 0..0.63 sek
Huelse-Projektil.gif [ 9.99 KiB | 9590-mal betrachtet ]
Code:
kg = 1; m = 1; sek = 1; grad = Pi/180;(*Dimensionen*)
d = 45*10^-4;(*Diameter*)
M = 547/1000000 kg;(*Masse Projektil*)
g = 981/100 m sek^-2;(*Erdbeschleunigung*)
p = 1293/1000 kg m^-3;(*Luftdicihte*)
J = 48 kg m^2 sek^-2;(*Energie Projektil in Joule*)
j = 1 kg m^2 sek^-2;(*Energie Hülse in Joule*)
V0 = Sqrt[2 J/M];(*Mündungsgeschwindigkeit*)
v0 = Sqrt[2 j/M];(*Hülsenauswurfsgeschwindigkeit*)
a = 0 grad;(*Winkel*)
h = 2 m;(*Höhe*)
ch = -2/100; cv = ch;(*Widerstandsbeiwert*)
Ah = Pi*d^2/4 m^2; Av = Ah;(*Fläche*)
tm = 1.3 sek; (* Flugzeit *)
X1 = NDSolve[{ch*Ah*p*rx'[t]^2/2 == M*rx''[t], rx'[0] == V0*Cos[a],
rx[0] == 0}, rx[t], {t, 0, tm}];(*horizontale Komponente Projektil*)
x1 = NDSolve[{ch*Ah*p*rx'[t]^2/2 == M*rx''[t], rx'[0] == v0*Cos[a],
rx[0] == 0}, rx[t], {t, 0, tm}];(*horizontale Komponente Hülse*)
Y1 = NDSolve[{-M*g + Sign[ry'[t]]*cv*Av*p*0.5*ry'[t]^2 == M*ry''[t],
ry'[0] == V0*Sin[a], ry[0] == h},
ry[t], {t, 0, tm}];(*vertikale Komponente Projektil*)
y1 = NDSolve[{-M*g + Sign[ry'[t]]*cv*Av*p*0.5*ry'[t]^2 == M*ry''[t],
ry'[0] == v0*Sin[a], ry[0] == h},
ry[t], {t, 0, tm}];(*vertikale Komponente Hülse*)
Show[ListPlot[
Table[Flatten[Evaluate[{rx[t] /. X1, ry[t] /. Y1}]], {t, 0, tm,
0.01 sek}], AspectRatio -> 1/3, Frame -> True,
PlotRange -> {{0, 260}, {0, 2}}, ImageSize -> 500],
ListPlot[Table[
Flatten[Evaluate[{rx[t] /. x1, ry[t] /. y1}]], {t, 0, tm,
0.01 sek}], AspectRatio -> 1/3, Frame -> True,
PlotRange -> {{0, 300}, {0, 2}}, ImageSize -> 500]] (*Plot*)
(*Yukterez*)
Illustrierend, Yukterez