Как объяснить отклонение траектории свободно падающего тела вблизи больших вертикальных скоростей воздуха (эффект сопротивления), и как это влияет на расчёт времени падения с большой высоты
Коротко — причина отклонения и изменение времени падения в том, что сила сопротивления зависит от относительной скорости тела и воздуха: она направлена против вектора относительной скорости v⃗−w⃗\vec v-\vec wv−w, а не обязательно строго вертикально. Из этого вытекают два эффекта: изменение вертикальной скорости (время падения) и появление горизонтальной компоненты ускорения (смещение). Основные формулы и выводы - Сила сопротивления (квадратичная модель): F⃗D=−12ρCDA ∣v⃗−w⃗∣(v⃗−w⃗),
\vec F_D = -\tfrac12\rho C_D A\,|\vec v-\vec w|(\vec v-\vec w), FD=−21ρCDA∣v−w∣(v−w),
где v⃗\vec vv — скорость тела, w⃗\vec ww — скорость воздуха, ρ\rhoρ — плотность воздуха, CDC_DCD — коэффициент сопротивления, AAA — поперечная площадь. - Вертикальное уравнение движения (ось вниз положительна), при произвольной вертикальной компоненте ветра w(t,z)w(t,z)w(t,z): mdvdt=−mg+mdwdt−12ρCDA ∣v−w∣(v−w).
m\frac{d v}{dt} = -mg + m\frac{dw}{dt} - \tfrac12\rho C_D A\,|v-w|(v-w). mdtdv=−mg+mdtdw−21ρCDA∣v−w∣(v−w).
Если w=w=w= const, то dwdt=0\frac{dw}{dt}=0dtdw=0 и удобно перейти к относительной скорости u=v−wu=v-wu=v−w: mdudt=mg−12ρCDA u∣u∣.
m\frac{du}{dt} = mg - \tfrac12\rho C_D A\,u|u|. mdtdu=mg−21ρCDAu∣u∣. - Терминальная (установившаяся) относительная скорость: mg=12ρCDA vt,rel2⇒vt,rel=2mgρCDA.
mg=\tfrac12\rho C_D A\,v_{t,rel}^2\quad\Rightarrow\quad v_{t,rel}=\sqrt{\frac{2mg}{\rho C_D A}}. mg=21ρCDAvt,rel2⇒vt,rel=ρCDA2mg.
Тогда скорость тела относительно земли при установившемся движении: v=w±vt,relv = w \pm v_{t,rel}v=w±vt,rel (знак по направлению движения). - В случае отсутствия ветра (w=0w=0w=0) аналитическое решение при нулевой начальной скорости (вниз положительно): v(t)=vt tanh (gtvt),y(t)=y0−vt2gln cosh (gtvt),
v(t)=v_t\:\tanh\!\Big(\frac{g t}{v_t}\Big),\qquad y(t)=y_0-\frac{v_t^2}{g}\ln\!\cosh\!\Big(\frac{g t}{v_t}\Big), v(t)=vttanh(vtgt),y(t)=y0−gvt2lncosh(vtgt),
где vt=vt,relv_t=v_{t,rel}vt=vt,rel. Время падения с высоты HHH находится из уравнения H=vt2gln cosh (gTvt).
H=\frac{v_t^2}{g}\ln\!\cosh\!\Big(\frac{gT}{v_t}\Big). H=gvt2lncosh(vtgT). - При постоянном вертикальном ветре www переход u=v−wu=v-wu=v−w даёт те же формы для u(t)u(t)u(t), а для координаты: y(t)=y0−wt−vt,rel2gln cosh (gtvt,rel),
y(t)=y_0 - w t - \frac{v_{t,rel}^2}{g}\ln\!\cosh\!\Big(\frac{g t}{v_{t,rel}}\Big), y(t)=y0−wt−gvt,rel2lncosh(vt,relgt),
и время падения TTT решается из H=wT+vt,rel2gln cosh (gTvt,rel).
H = wT + \frac{v_{t,rel}^2}{g}\ln\!\cosh\!\Big(\frac{gT}{v_{t,rel}}\Big). H=wT+gvt,rel2lncosh(vt,relgT). Интерпретация и практические следствия - Если воздух поднимается (восходящий поток, w>0w>0w>0), относительная скорость вниз уменьшается, сопротивление слабее и время падения увеличивается; при www близком к vt,relv_{t,rel}vt,rel тело может фактически «зависнуть» или даже подняться. При сильном нисходящем потоке время падения уменьшается. - Горизонтальные скорости воздуха дают горизонтальную компоненту силы сопротивления; тело постепенно приобретает горизонтальную скорость, стремящуюся к скорости воздуха, поэтому появится значительное горизонтальное смещение ~воздух·время падения. При градиентах ветра и турбулентности траектория дополнительно усложняется. - Для больших высот и сильного сопротивления реальная траектория и время падения обычно требуют численного решения системы ОДУ (или моделирования с турбулентностью). Для простых оценок можно применять приведённые аналитические формулы при предположении постоянных параметров и квадратичного сопротивления. Если нужно, могу дать пример численного расчёта времени падения для конкретных m, CD, A, ρ, w, Hm,\;C_D,\;A,\;\rho,\;w,\;Hm,CD,A,ρ,w,H.
Основные формулы и выводы
- Сила сопротивления (квадратичная модель):
F⃗D=−12ρCDA ∣v⃗−w⃗∣(v⃗−w⃗), \vec F_D = -\tfrac12\rho C_D A\,|\vec v-\vec w|(\vec v-\vec w),
FD =−21 ρCD A∣v−w∣(v−w), где v⃗\vec vv — скорость тела, w⃗\vec ww — скорость воздуха, ρ\rhoρ — плотность воздуха, CDC_DCD — коэффициент сопротивления, AAA — поперечная площадь.
- Вертикальное уравнение движения (ось вниз положительна), при произвольной вертикальной компоненте ветра w(t,z)w(t,z)w(t,z):
mdvdt=−mg+mdwdt−12ρCDA ∣v−w∣(v−w). m\frac{d v}{dt} = -mg + m\frac{dw}{dt} - \tfrac12\rho C_D A\,|v-w|(v-w).
mdtdv =−mg+mdtdw −21 ρCD A∣v−w∣(v−w). Если w=w=w= const, то dwdt=0\frac{dw}{dt}=0dtdw =0 и удобно перейти к относительной скорости u=v−wu=v-wu=v−w:
mdudt=mg−12ρCDA u∣u∣. m\frac{du}{dt} = mg - \tfrac12\rho C_D A\,u|u|.
mdtdu =mg−21 ρCD Au∣u∣.
- Терминальная (установившаяся) относительная скорость:
mg=12ρCDA vt,rel2⇒vt,rel=2mgρCDA. mg=\tfrac12\rho C_D A\,v_{t,rel}^2\quad\Rightarrow\quad
v_{t,rel}=\sqrt{\frac{2mg}{\rho C_D A}}.
mg=21 ρCD Avt,rel2 ⇒vt,rel =ρCD A2mg . Тогда скорость тела относительно земли при установившемся движении: v=w±vt,relv = w \pm v_{t,rel}v=w±vt,rel (знак по направлению движения).
- В случае отсутствия ветра (w=0w=0w=0) аналитическое решение при нулевой начальной скорости (вниз положительно):
v(t)=vt tanh (gtvt),y(t)=y0−vt2gln cosh (gtvt), v(t)=v_t\:\tanh\!\Big(\frac{g t}{v_t}\Big),\qquad
y(t)=y_0-\frac{v_t^2}{g}\ln\!\cosh\!\Big(\frac{g t}{v_t}\Big),
v(t)=vt tanh(vt gt ),y(t)=y0 −gvt2 lncosh(vt gt ), где vt=vt,relv_t=v_{t,rel}vt =vt,rel . Время падения с высоты HHH находится из уравнения
H=vt2gln cosh (gTvt). H=\frac{v_t^2}{g}\ln\!\cosh\!\Big(\frac{gT}{v_t}\Big).
H=gvt2 lncosh(vt gT ).
- При постоянном вертикальном ветре www переход u=v−wu=v-wu=v−w даёт те же формы для u(t)u(t)u(t), а для координаты:
y(t)=y0−wt−vt,rel2gln cosh (gtvt,rel), y(t)=y_0 - w t - \frac{v_{t,rel}^2}{g}\ln\!\cosh\!\Big(\frac{g t}{v_{t,rel}}\Big),
y(t)=y0 −wt−gvt,rel2 lncosh(vt,rel gt ), и время падения TTT решается из
H=wT+vt,rel2gln cosh (gTvt,rel). H = wT + \frac{v_{t,rel}^2}{g}\ln\!\cosh\!\Big(\frac{gT}{v_{t,rel}}\Big).
H=wT+gvt,rel2 lncosh(vt,rel gT ).
Интерпретация и практические следствия
- Если воздух поднимается (восходящий поток, w>0w>0w>0), относительная скорость вниз уменьшается, сопротивление слабее и время падения увеличивается; при www близком к vt,relv_{t,rel}vt,rel тело может фактически «зависнуть» или даже подняться. При сильном нисходящем потоке время падения уменьшается.
- Горизонтальные скорости воздуха дают горизонтальную компоненту силы сопротивления; тело постепенно приобретает горизонтальную скорость, стремящуюся к скорости воздуха, поэтому появится значительное горизонтальное смещение ~воздух·время падения. При градиентах ветра и турбулентности траектория дополнительно усложняется.
- Для больших высот и сильного сопротивления реальная траектория и время падения обычно требуют численного решения системы ОДУ (или моделирования с турбулентностью). Для простых оценок можно применять приведённые аналитические формулы при предположении постоянных параметров и квадратичного сопротивления.
Если нужно, могу дать пример численного расчёта времени падения для конкретных m, CD, A, ρ, w, Hm,\;C_D,\;A,\;\rho,\;w,\;Hm,CD ,A,ρ,w,H.