7.2: Опір повітря, пропорційний швидкості
Як і в попередньому розділі, я напишуx -компонент рівняння руху, і першого і другого часу інтегралів, в лівій колонці, іy -компонент в правій колонці. x-складова опору повітря на одиницю маси єγ˙x іy -компонент єγ˙y. γОсь константа демпфування, визначена в розділі 6.3. Уx - іy -складових початкової швидкості є, відповідно,V0cosα іV0sinα. Слід легко бачити, що рівняння руху та їх тимчасові інтеграли такі:
Горизонтальний | Вертикальний | |
---|---|---|
¨x=−γ˙x | ¨y=−g−γ˙y | 7.2.1a,b |
˙x=u=V0cosα.e−γt | ˙y=v=V0sinαe−γt−ˆv(1−eγt)whereˆv=g/γ | 7.2.2a,b |
x=x∞(1−e−γt)where x∞=V0cosαγ | y=1γ(V0sinα+ˆv)(1−e−γt)−ˆvt | 7.2.3a,b |
(Якщо це не «легко видно», для горизонтального руху зверніться до розділу 3 глави 6, особливо Рівняння 6.3.2, 6.3.3 та 6.3.5, а для вертикального руху зверніться до розділу 6, розділ 3b, особливо Рівняння 6.3.24, 6.3.25 та 6.3.27.) Буде видно, щоt→∞,u→0,v→−ˆv,x→x∞. xy-рівняння до траєкторії єt -елімінант рівнянь6.2.3a і6.2.3b. Після невеликої кількості алгебри це виявляється:
y=x(V0sinα+ˆv)V0cosα+ˆvγln(1−xx∞).
Діапазон на горизонтальній площині знаходить, встановившиy = 0, щоб отримати або
x=−Aln(1−xx∞)
або
x=x∞(1−e−xA),
де
A=ˆvV0cosαγ(V0sinα+ˆv),x∞=V0cosαγ
і
ˆv=gγ.
Припустимо
- V0= 20 мс -1
- α= 50∘
- g=9.8 мс -2
- γ=1.96с -1 (∴=5 мс −1)
Тоді A = 1,613 870 65 м
а x_{\infty} = 6,55905724 м.
Спробуйте знайти діапазон на горизонтальній площині, використовуючи\ref{eq:7.2.5} або Рівняння\ref{eq:7.2.6}, або, до дев'яти значущих цифр. Яке рівняння найкраще працює? Ньютон-Рафсон може провалитися з дурною першою здогадкою - але зробити досить розумну першу здогадку не повинно бути важко. Не повинен вам розповідати, але цифра VII.2 була розрахована за даними цього прикладу.
Я роблю відповідь 6.437 584 2 м.
Ось більш складна проблема.
Загальновідомо, що при відсутності опору повітря максимальний діапазон на горизонтальній площині здійснюється шляхом вибору початкової висоти запуску, щоб бути \alpha = 45 ^{\circ}. Що робити, якщо є опір повітря, з постійною демпфуванням \gamma? Яким же тоді повинен бути кут запуску для досягнення найбільшої дальності на горизонтальній площині? Дано рівняння\ref{eq:7.2.6} x=x_{\infty}(1-e^{\frac{-x}{A}}), для якого значення \alpha є x найбільшим?
Рішення
Рівняння \ref{eq:7.2.6}, записане повністю, є
x=\frac{V_{0}cos\alpha}{\gamma}\left[1-exp\left(\frac{-\gamma(V_{0}\sin\alpha+\hat{v})x}{\hat{v}V_{0}\cos\alpha}\right)\right]. \tag{7.2.7}\label{eq:7.2.7}
Це можна написати
x=acos\alpha\left[1-exp\left(-\frac{(b \sin\alpha + 1)x}{a \cos\alpha}\right)\right], \tag{7.2.8}\label{eq:7.2.8}
де a=\frac{V_{0}}{\gamma} і b=\frac{V_{0}}{\hat{v}}=\frac{\gamma V_{0}}{g}. Ми повинні знайти, для чого значення \alpha є x найбільшим. Це здається досить простою проблемою, але на даний момент я не можу знайти хорошого способу її вирішення. Якщо у когось є підказка, дайте мені знати (jtatum@uvic.ca). Тим часом найкраще, що я можу запропонувати, для нашого конкретного числового прикладу, обчислити діапазон x, для декількох значень \alpha і побачити, де він проходить максимум. Для нашого конкретного числового прикладу, a = 10,204 081 63 м і b = 4. Ось графік діапазону проти кута запуску, для початкової швидкості 20 мс −1. Кут запуску близько 23 ^{\circ} 59′ дає дальність близько 8.4635 м Для заданого \gamma і g, оптимальний кут запуску залежить від швидкості запуску V_{0}. Це інтуїтивно очевидно?