Delphi - Моделирование полета тела с учетом сопротивления воздуха

[W3Avenger]

Сообщение от type_Oleg

Держите .
Скорость задаю не через угол, а через проекции вектора скорости.

Что значит домножение модуля вектора скорости на некую k в mod1??

[type_Oleg]

Сообщение от W3Avenger

Что значит домножение модуля вектора скорости на некую k в mod1??

Имеете в виду в функции AxAirVect ?
Сначала в функции fAxAir вычисляем абсолютное значение ( модуль ) аэродинамического ускорения ( то есть торможения ).
Но это ускорение - величина векторная, имеющая направление навстречу вектору скорости. Для того, чтобы создать вектор, и используется AxAirVect , в которой kVect.
А вообще функция kVect используется и в других местах, например при вычислении следующих pos, vel по предыдущим, в процедуре MoveT.

kVect - это вообще функция умножения любого вектора на любое число.

Вот, закомментил часть кода:

repeat
  t:=t+dt;
  mAxAir:=fAxAir(M,S,ro,Cx,ModQuaVect(vel));  // модуль аэродин. торможения
  axAir:=AxAirVect(mAxAir,vel);  // вектор аэродин. торможения
  axSum:=SumVect(axAir,axGr);    // полный вектор ускорения = уск. св. падения + аэродин. торможения
  MoveT(dt,axSum,pos,vel);       // изменение вектора положения и вектора скорости
  Memo1.Lines.Add(FloatToStr(t));
  Memo2.Lines.Add(FloatToStr(pos.X));
  Memo3.Lines.Add(FloatToStr(pos.Y));
  Memo4.Lines.Add(FloatToStr(vel.X));
  Memo5.Lines.Add(FloatToStr(vel.Y));
 until (pos.Y<0)or(t>TMAX);    //  .. пока не упадет, или пока t<MAX
// t>TMAX - защита от дурака, если зафутболить тело с космической скоростью, то не дождешься пока оно упадет

Проверил на примере падения капель дождя.
Капля 0.5 мм у меня падает с установившейся скоростью ~ 4 м/сек, капля 5 мм - 12.5 м/сек . Сравнить - отсюда:
Скорость падения дождевых капель диаметром 0,5 мм на уровне моря и без ветра составляет от 2 до 6,6 метров в секунду, в то время как капли диаметром 5 мм имеют скорость от 9 до 30 метров в секунду
Вроде нормально..

Кстати, обнаружился в этом случае небольшой баг - нельзя задавать начальную скорость точно =0 ( и Vx=0; Vy=0 ), будет деление на 0. Пришлось для капель задать небольшую Vx=0.1.

[W3Avenger]

Сообщение от type_Oleg

Для того, чтобы создать вектор, и используется AxAirVect ...

Так, внимание, а если я, как дурак вырежу половину вашей работы и попробую задать скорость не компонентами, не векторами, а через модуль скорости и угол.... Какие могут быть варианты? Считать через массив с определенным приращением не получится?

[type_Oleg]

Сообщение от W3Avenger

.. и попробую задать скорость не компонентами, не векторами, а через модуль скорости и угол...

Только задавать скорость не через Vx, Vy , а через V и угол, а вообще векторы сохранить ?
Или вообще отказаться от типа tVector (и для координат, и для скорости, и для ускорения), и использовать отдельные переменные X,Y,Vx,Vy,Ax,Ay ?

Вобщем-то и так, и так можно. Только сами переделывайте.

Только я не понял, что значит Считать через массив с определенным приращением ...