Menambahkan seretan udara ke persamaan lintasan bola golf

Saya sedang mengembangkan permainan golf 2D di VB.NET 2005, tapi saya terjebak pada bagaimana menerapkan hambatan udara atau angin yang harus memengaruhi bola.

Saya sudah memiliki persamaan ini untuk proyektil:

• untuk kecepatan awal bola golf saat dipukul atau ditembakkan$v_0$

• Komponen vertikal dan horizontal kecepatan bola golf:

  $$\begin{align} v_x &= v_0 cos(\theta) \\ v_y &= v_0 sin(\theta) - gt* \end{align}$$

• Jarak bola golf vertikal dan horizontal:

  $$\begin{align} x &= v_0cos(\theta)t \\ y &= v_0sin(\theta) t - (0.5)gt^2 \end{align}$$

Bagaimana cara menambahkan hambatan udara ke persamaan ini untuk memengaruhi kecepatan bola golf dengan benar? Saya tidak tahu bagaimana melakukannya, apakah ada yang bekerja dengan persamaan yang sama?

Saya tidak yakin apakah bahkan ada bentuk tertutup untuk drag atau angin, tetapi cukup mudah untuk mensimulasikan secara bertahap (seperti semua perpustakaan fisika lakukan):

1. atur kondisi awal Anda:

   $$x, y, v_x, v_y \; (\text{for }t=0)$$

2. perbarui posisi:

   $$x = x + (v_x \times dt) \\ y = x + (v_y \times dt)$$

   (di mana dt adalah waktu berlalu sejak pembaruan terakhir, alias waktu delta)

3. hitung pembantu kecepatan ini:

   $$\begin{align} v^2 &= (v_x)^2 + (v_y)^2 \\ \lvert v \rvert &= \sqrt{v^2} \end{align}$$

   $\lvert v \rvert$$v$

4. menghitung gaya seret:

   $$f_{drag} = c \times v^2$$

   (di mana c adalah koefisien gesek kecil! )

5. mengumpulkan kekuatan:

   $$\begin{align} f_x &= \left(-f_{drag} \times {v_x \over \lvert v \rvert}\right) \\ f_y &= \left(-f_{drag} \times {v_y \over \lvert v \rvert}\right) + (-g \times mass) \end{align}$$

   $mass$

6. perbarui kecepatan:

   $$v_x = v_x + f_x \times \frac{dt}{mass} \\ v_y = v_y + f_y \times \frac{dt}{mass}$$

Itu pada dasarnya Metode Euler untuk mendekati fisika itu.

Sedikit lebih banyak tentang bagaimana simulasi seperti yang diminta dalam komentar:

• $(t = 0)$

Ini pada dasarnya sama dengan rumus lintasan dasar Anda di mana setiap kemunculan t diganti dengan 0.

• $KE = 0.5m(V^2)$$t$$v^2$

• $PE = m \times g \times y$

• $(x,y)$$t_1$$t = 0$$t = t_1$

• $(x,y)$$t_1$$t_2$$t_1 \lt t_2$$t_1$$t_2$

Pseudo-Code:

simulate(v0, theta, t1)
  dt = 0.1
  x = 0
  y = 0
  vx = v0 * cos(theta)
  vy = v0 * sin(theta)
  for (t = 0; t < t1; t += dt)
    x += vx * dt
    y += vy * dt
    v_squared = vx * vx + vy * vy
    v_length = sqrt(v_squared)
    f_drag = c * v_squared
    f_grav = g * mass
    f_x = (-f_drag * vx / v_length)
    f_y = (-f_drag * vy / v_length) + (-f_grav)
    v_x += f_x * dt / mass
    v_y += f_y * dt / mass
  end for
  return x, y
end simulate