Мне нужно самое простое уравнение подъемной силы, которое когда-то было решено с помощью Mathcad, дает реальную вертикальную скорость самолета и неявно высоту h(t).
Также сопротивление экв. (1) приводит к хорошему решению, Vh(t), которое увеличивается и, наконец, достигает предела (максимальная горизонтальная скорость) и остается на этом уровне в течение всего полета, подъемная сила экв. (2) стабилизируется при подъеме (t) - m * g = ct. > 0 и, как следствие, Vh(t) продолжает расти бесконечно, поскольку уравнение типа m * dVv(t)/dt = ct. приводит к решению Vv(t), линейно растущему со временем.
Вопрос: Совершенно очевидно, что вертикальная скорость самолета Vh(t) не может расти бесконечно. Как я могу стабилизировать его до постоянного значения. Что я должен добавить в экв. лифта?
Перетащите экв. : m * dVh(t)/dt = T - Drag(t) (1),
Перетащите (t) = 0,5 * Cd * r * S * (Vh (t) + Vw (t)) ^ 2,
Поднимите экв. : m * dVv(t)/dt = Lift(t) - m * g (2),
Lift(t) = 0,5 * Cl * r * S * (Vh(t) + Vw(t))^2,
где:
Во-первых, нужно добавить эффект плотности : плотность воздуха меняется с высотой, и это влияет как на подъемную силу, так и на сопротивление (по крайней мере, если вы моделируете воздушно-реактивные двигатели).
Для лифта , использовать
Для перетаскивания , вы должны только максимально упростить уравнение. Ваше уравнение еще проще и не моделирует увеличение сопротивления при увеличении подъемной силы. Простейшее практическое уравнение будет выглядеть так:
Больший подъем потребует большей тяги и ограничит ускорение вверх. Подъем выше уменьшит тягу пропорционально плотности и ограничит возможный диапазон скоростей.
Номенклатура:
коэффициент подъемной силы (обычно от 0 до 1,5)
3.14159
удлинение крыла (отношение размаха к средней хорде)
коэффициент Освальда крыла (используйте 0,8, если сомневаетесь)
коэффициент аэродинамического сопротивления при нулевой подъемной силе (используйте 0,02, если сомневаетесь)
Фриман