enter image description here

For a given launch velocity $v$ and launch angle $\theta$, the trajectory of a projectile may be described by the standard formula $$y=x\tan\theta-\frac {gx^2}{2v^2}\sec^2\theta$$

For different values of $\theta$ what is the envelope of the different trajectories? Is it a parabola itself?

The standard solution to this "envelope of safety" problem is to state the formula as a quadratic in $\tan\theta$ and set the discriminant to zero. The resulting relationship between $x,y$ is the envelope.

This question is posted to see if there are other approaches to the solution.

Edit 1

Thanks for the nice solutions from Jack and Blue, received so far. From the solution of the envelope it can be worked out that the envelope itself corresponds to the right half of the trajectory of a projectile launched at $(-\frac{v^2}{g^2},0)$ at a launch angle $\alpha=\frac{\pi}4$ and a launch velocity $V=v\sqrt2$. This means that both vertical and horizontal components of the launch velociy are equal to $v$. It would be interesting to see if these conclusions can be inferred from the problem itself by inspection and without first solving it. If so, then this would form another solution.

See also this other question posted subsequently.

  • $\begingroup$ I think it is more like an ellipse. $\endgroup$ – Archis Welankar Oct 24 '15 at 10:05
  • $\begingroup$ @ArchisWelankar - What makes you think that might be the case? $\endgroup$ – hypergeometric Oct 24 '15 at 10:13
  • $\begingroup$ When the angle is very small like 10 then the curve becomes elongated and more like an ellipse than a parabola. $\endgroup$ – Archis Welankar Oct 24 '15 at 10:25
  • $\begingroup$ @ArchisWelankar - 1. The curve or trajectory of the projectile is always a parabola (even when elongated, at small launch angles) as given by the formula. 2. An ellipse doesn't have to be elongated. 3. The question refers to the envelope of different trajectories and not any one particular trajectory. $\endgroup$ – hypergeometric Oct 24 '15 at 10:34
  • $\begingroup$ "The standard solution to this "envelope of safety" problem is to state the formula as a quadratic in tanθ and set the discriminant to zero. " I don't understand this approach, can someone explain ? $\endgroup$ – Dat Sep 11 at 9:38

Yes it is. To find the envelope, we just have to find the intersections between two trajectories associated to two slightly different angles. If we solve $$ x\tan\theta -\frac{gx^2}{2v^2}\sec^2(\theta) = x\tan(\theta+\varepsilon) -\frac{gx^2}{2v^2}\sec^2(\theta+\varepsilon) $$ we get $x=0$ or $$ x=\frac{2v^2}{g}\cdot\frac{\tan(\theta)-\tan(\theta+\varepsilon)}{\sec^2(\theta)-\sec^2(\theta+\varepsilon)}$$ and by letting $\varepsilon\to 0$ we get $x=\frac{v^2}{g}\cot(\theta)$, from which $y=\frac{v^2}{2g}\left(2-\frac{1}{\sin^2(\theta)}\right) $.

It follows that the equation of the envelope is given by: $$ y = \frac{v^2}{2g}\left(1-\left(\frac{gx}{v^2}\right)^2\right)=\frac{v^2}{2g}-\frac{g}{2 v^2}\,x^2$$ that clearly is a parabola with vertex in $\left(0,\frac{v^2}{2g}\right)$ through the points $\left(\pm\frac{v^2}{g} ,0\right)$.

Geometric sketch of the projectile motion

We may notice that the envelope and the trajectory with $\theta=\frac{\pi}{4}$ are homothetic, and the dilation ratio is just $2$. The vertices of the trajectories lie on an ellipse that is tangent to the envelope parabola, with centre at $\left(0,\frac{v^2}{4g}\right)$, a vertex in the origin and a vertex at $\left(\frac{v^2}{2g},\frac{v^2}{4g}\right)$.

  • 3
    $\begingroup$ Thank you for a beautiful and creative solution! And nice graphics! A big +1! $\endgroup$ – hypergeometric Oct 24 '15 at 13:29
  • $\begingroup$ Please check out comments under Edit 1 in the question and see if you have any further thoughts. $\endgroup$ – hypergeometric Oct 24 '15 at 14:21

@Jack provides a very nice and intuitive derivation of the envelope as the points of intersection of infinitely-close members of the curve family. The Wikipedia "Envelope" entry provides this less-illuminating abstraction:

The envelope of the family [of curves parameterized by $t$ is] the set of points for which $$F(t, x, y) = \frac{\partial F}{\partial t}(t,x,y) = 0 \tag{$\star$}$$ for some value of $t$ [...].

In $(\star)$, $F$ is the function that, when set equal to $0$, defines each curve in the family. For the question at hand, we have (with parameter $\theta$ instead of $t$) $$F(\theta,x,y) = -y + x\tan\theta -\frac{g x^2}{2v^2}\sec^2\theta \tag{1}$$ Therefore, differentiating with respect to $\theta$ gives $$\frac{\partial F}{\partial \theta}(\theta,x,y) =x\sec^2\theta -\frac{g x^2}{2v^2}\cdot 2 \sec^2\theta\tan\theta = \frac{x\sec^2\theta}{v^2} ( v^2 - g x \tan\theta) \tag{2}$$ Solving $\partial F/\partial \theta = 0$ for $\theta$ (noting that $\sec\theta$ never vanishes) gives $$\tan\theta = \frac{v^2}{g x} \qquad\text{so that}\qquad \sec^2\theta = 1 + \tan^2\theta = \frac{g^2 x^2 + v^4}{g^2 x^2}$$

Substituting into $(1)$, and setting $F=0$, we have

$$y = x\;\frac{v^2}{gx} - \frac{g x^2}{2 v^2}\;\frac{g^2 x^2 + v^4}{g^2 x^2} = \frac{v^2}{2g} - \frac{g x^2}{2v^2} = \frac{v^4-g^2x^2}{2gv^2}$$

which agrees with Jack's answer.

  • $\begingroup$ Thanks for working out the solution from the article. Nice find! +1. $\endgroup$ – hypergeometric Oct 24 '15 at 13:55
  • $\begingroup$ Please check out comments under Edit 1 in the question and see if you have any further thoughts. $\endgroup$ – hypergeometric Oct 24 '15 at 14:21

Here's another approach. Not as elegant as Jack's though. Perhaps similar to Blue's. (Edit: Revised slightly to use substitution of $H=v^2/2g$ instead of $\lambda=g/v^2$ for more convenient reference, per solution by Narasimham)

This approach maximises the value of $y$ for any given value of $x$.

For a given value of $x$, say $x=k$, $$\begin{align} y&=kT-\frac {k^2}{4H}(1+T^2)\qquad\text{where $T=\tan\theta, H=\frac {v^2}{2g}$}\\ \frac{dy}{dT}&=k-\frac {k^2}{4H}\cdot 2T=0\qquad\text{when $T=\frac {2H}k$}\\ \text{At $T=\frac {2H}k$}:\qquad\qquad y&=k\cdot\frac{2H}k-\frac{k^2} {4H}\left(1+\frac {4H^2}{k^2}\right) \color{lightgrey}{=2H-\frac {k^2}{4H}-H}\\ &=H-\frac {k^2}{4H} \end{align}$$

Hence the equation of the envelope is $$y=\frac {v^2}{2g}-\frac {gx^2}{2v^2}\quad\blacksquare$$

  • $\begingroup$ This nicely exploits the fact that, in this particular case, we know that the envelope sits on top of the family of curves, and therefore consists of the points that maximize $y$ for a given $x$, over all $t$. (Of course don't always have such knowledge, but no matter.) Re-writing in terms of parameter $T$ is also a helpful simplification. And, to be explicit: the connection with the process in my answer is that the criterion for maximizing $y$ ---that is, setting $dy/dT = 0$--- corresponds directly to $\partial F/\partial \theta = 0$, because $y$ drops out of the latter equation. $\endgroup$ – Blue Oct 24 '15 at 21:38
  • $\begingroup$ Thanks for your kind comment! $\endgroup$ – hypergeometric Oct 25 '15 at 5:57

I shall outline the method how we get the parabola. Usual notation

$$ \ddot y = -g, \dot y = - g t + v \sin \theta , y =- g t^2/2 + v t \sin \theta\, +0 \tag {1} $$ $$\ddot x = 0, \dot x = v \cos \theta = const , x = v t \cos \theta +0 \tag{2}$$

Eliminating time $t$ between (1),(2) you got this parabola equation already.

Let $\tan \alpha = T$; Parabola equation in other words

$$ y = x T - g/2 * ( x/ v \cos\theta)^2 = x T - ( g x^2/2 v^2) ( 1+T^2) \tag{3} $$

Differentiate partially with respect to $T$ and simplify, $ T = v^2 / gx \tag{4} $

Eliminate $T$ between (3) and (4)

$$ y = v^2/2g - g x^2 /(2 v^2) = H - x^2/(4H) \tag{5}, $$

same as what was obtained before by Jack and Blue.

If you denote height reached by projectile on vertical firing $ H = v^2/(2g) \tag{6} $ you would notice that the envelope is profiled exactly as a parabolic mirror with focus at gun delivery point, focal length is exactly H. Vertical force of gravity is acting like light :)..

The above procedure method is indicated by Blue in Wiki, is referred to as C-discriminant method to obtain envelopes and singular solutions.

Like what you said in your edit and I about mirror, they are ploys for remembering curves using similarities..

enter image description here

For the image I took values of $ g=9.8 m/s^2 , v = 2 m/s $

The partial differentiation way and elimination is the right way to look at, perhaps not as what you said (standard solution.. way).


The answer to your second question, i.e., to determine if it is going to be a parabola envelope without going through all of analysis... I can only reply with extended C-discriminant, strengthening the same result by another path.

$p$ discriminant method is also relevant, but I defer it, but best is to refer to differential calculus books of authors e.g., A.R. Forsythe.

I shall expand on the C-discriminant, wherein a two parameter system of equations variation of any of the two parameters leads to the same parabola envelope. This is in reply to your question , Why do you also plot partial concentric circles?

Well, they are circles alright, but not concentric circles. They expand and come down slowly with time. What you can see as plotted are the peripheries traced out descending with time .

But first quickly watch some fireworks to see what I am discussing about:

Expanding & Descending Fire-Works Circle Periphery

It is common experience to see fire ball bright splinters expanding to bigger circles as the entire cluster comes down slowly with time. The periphery is a portion of a circle whose radius increases. Center of circle is always descending by gravity.

The two parameters are $ \theta, t $ angle of elevation at first burst or fire, and time $t$.

By C-discriminant method the parabola envelope is the eliminant of either $\theta$ variable or

$$ F(x,y,\theta) =0 ,\, F_{\theta } (x,y,\theta) =0 \tag{6} $$

or $t$ time variable.

$$ F(x,y, t) =0 ,\, F_{t } (x,y, t) =0 \tag{7} $$

The first one is already discussed, the second one is expanding/descending fireworks circles as already stated.

In the latter case working is:

$$ x = v t \cos \theta , y = v t \sin \theta - g t^2/2 \tag{8}$$

$$ (\frac{x}{vt})^2 + (\frac{y+ gt^2/2}{vt}) ^2 = 1 \tag{9}$$

$$ x^2 + ( y + g t^2/2)^2 = v^2t^2 \tag{10} $$

which is a Circle.

To find its envelope, as before partially differentiate with respect to time $t$ and cancel $ 2t$ on either side of equation , bring $ v^2/g$ to right side :

$$ y + gt^2/2 = v^2/g \tag {11} $$

$$ x^2 + (v^2/g)^2 = 2 v^2/g * ( v^2/g-y) \tag {12} $$

$$ x^2 + ( 2 H)^2 += 4 H ( v^2/g -y )\tag{13} $$

$$ x^2 = 4 H ( H-y) \tag{14} $$

which is the same parabola envelope obtained earlier with $\theta $ as parameter. End points $ ( x=0, y=H ; x= 2 H, y= 0 ) $

Although circle end traces are visible in a fireworks display it needs imagination as before with variable gun barrel angle to see that each circle is tangent to a fixed envelope. Hope you enjoyed it. FireWorksEnvelope

  • $\begingroup$ Thanks for your solution and insight on the parabolic mirror! (+1). Nice graphics too! In addition to different parabolas, why do you also plot partial concentric circles? $\endgroup$ – hypergeometric Oct 25 '15 at 5:56
  • $\begingroup$ Also, is there an argument which would directly conclude that the envelope is a parabola with focus at launch point and tip at the peak of the vertical launch, without first going through the analysis? $\endgroup$ – hypergeometric Oct 25 '15 at 6:03
  • $\begingroup$ The question in my comment above has been answered here. $\endgroup$ – hypergeometric Sep 30 '18 at 17:14

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.