{"id":250,"date":"2012-09-07T17:12:19","date_gmt":"2012-09-07T23:12:19","guid":{"rendered":"http:\/\/www.ghyzmo.com\/?p=250"},"modified":"2017-11-08T16:58:51","modified_gmt":"2017-11-08T23:58:51","slug":"ballistics-of-galileo","status":"publish","type":"post","link":"http:\/\/www.ghyzmo.com\/ballistics-of-galileo\/","title":{"rendered":"Ballistics of Galileo"},"content":{"rendered":"

\u00a9Fernando Caracena 31 August 2012<\/em><\/p>\n

Independent motions<\/em><\/strong><\/h2>\n

One of Galileo's great discoveries is that motion in three dimensions can be decomposed into three separate motions along each of three mutually orthogonal directions. Specifically, up and down motions are governed by the acceleration of the Earth's gravity, whereas the two other, horizontal motions can be unaccelerated to the extent that friction becomes negligible. For these reasons, a projectile launched into a free trajectory, moves in a two-dimensional plane. Once launched, no acceleration acts across to the right or left of the launch plane, so that the projectile always remains in the same plane.<\/p>\n

Vertical motion<\/em><\/strong><\/h2>\n

Free fall<\/em><\/h3>\n

One of the first questions one can solve in Galileo's physics is: how long does it take (T) for an object, released from rest (horizontal speed =0), to fall though a given distance (h)?\u00a0 Galileo is said to have demonstrated that objects of different sizes fall at the same rate under gravity. This is true as long as air resistance is not a problem. On the Moon, an astronaut can drop a hammer and a feather together, and they will hit the regolith<\/a> at the same time; however, on the Earth, air resistance is not negligible for the feather. Air resistance is negligible when you drop a baseball and bowling ball together. Having the bottoms lined up horizontally when released, they will hit the ground at the same time.<\/p>\n

The acceleration of gravity<\/a> is a \"constant\" that has been measured many times. Actually, it varies by very small departures from its standard value (g) because of variations in the Earth's local density and, the Earth's spin. For most practical purposes, we can take the value of the acceleration of the Earth's gravity as the following:<\/p>\n

g=9.8 m\/s\/s (or 32 ft\/s\/s).<\/p>\n

Locally, we define the coordinate axis as pointing up and the value of height as z. In this case, the acceleration of gravity, which acts downward, has a negative value,<\/p>\n

a=-g\u00a0\u00a0\u00a0\u00a0 (constant).\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (1)<\/p>\n

The equations of motion referred to this axis are the following:<\/p>\n

w= w0<\/sub> -g * t\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (2a)<\/p>\n

z=z0<\/sub>+w0<\/sub>*t-1\/2 * g *t2<\/sup> , \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (3b)<\/p>\n

where the initial conditions require the following substitutions:<\/p>\n

w<\/sup>0 <\/sub>= 0 and z0<\/sub> =h.\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (4)<\/p>\n

In this case, the equations of motion are reduced to the following:<\/p>\n

w=-g * t\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (2b)<\/p>\n

z=h- 1\/2 * g *t2<\/sup>,\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (3b)<\/p>\n

where the bottom of the trajectory is taken as z=0.<\/p>\n

from (3b) we can solve for the fall time by setting z=0,<\/p>\n

h- 1\/2 * g *T2<\/sup>=0, \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (5)<\/p>\n

which is readily solved to be,<\/p>\n

 <\/p>\n

T= (2 *h \/g)<\/span><\/sub>1\/2<\/sup><\/span>.\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (6)<\/p>\n

N.B., that the exponent \"1\/2\" in the above means square root.<\/p>\n

The value of solving a problem, such as the above, algebraically, is the economy of a general solution, both in terms of units and values for various astronomical bodies. For example, the lunar acceleration of gravity is a about 1.6 m\/s\/s; on Mars, it is about 3.8 m\/s\/s. The algebraic solution, (6), applies to the surface of any planet and a release point of any height.<\/p>\n

Choose Earth's value of g and plug in a height of 16 feet into (6) and the value of g in corresponding units (32 m\/s\/s), and you will get the result of 1 second for the time of fall. Also note that we can answer the question of how fast the object impacts the ground by inserting the value, 1 s, in (2b), which results in 32 ft\/s. This shows that we can use intuition only so far before a proficiency with algebra adds a whole lot more.<\/p>\n

Upward shot<\/em><\/h3>\n

The problem that is opposite to the above one is where an object is propelled straight up from the surface. How high will it go and how long will it take to get there?<\/p>\n

Instead of the the initial conditions in (4), those that apply here are<\/p>\n

z0<\/sub> =0 and w0<\/sub> is not zero.\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (7)<\/p>\n

Under these initial conditions, the equations for a projectile shot straight up are:<\/p>\n

w=w0<\/sub> - g * t\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (2c)<\/p>\n

and<\/p>\n

z=\u00a0w0<\/sub>* t - 1\/2 * g *t2<\/sup>. \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (3c)<\/p>\n

To solve for the time of rise, note that at the highest point, the upward velocity falls to zero, which in (2c) gives the following:<\/p>\n

0=w0<\/sub> - g * T,<\/p>\n

or<\/p>\n

T=w0<\/sub>\/g.\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (8)<\/p>\n

Plug this value of time into (2c) and it results in the maximum height at time T,<\/p>\n

h= * T - 1\/2 *g *T2<\/sup>,<\/p>\n

h=w0<\/sub>2<\/sup>\/g- 1\/2 *g *(w0<\/sub>\/g)2<\/sup>,<\/p>\n

h=1\/2 *w0<\/sub>2<\/sup>\/g . \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (9)<\/p>\n

Plugging the value of the rise time (T) in (8) into (9), and solving for T, we get<\/p>\n

T= (2 *h \/g)1\/2<\/sup>, \u00a0which is exactly the same as the drop time (6).<\/p>\n

\"\"<\/a>

Fig. 1 Projectile launched straight with upward component of velocity of 70 ft\/s.<\/p><\/div>\n

<\/h3>\n

<\/h3>\n

<\/h3>\n

<\/h3>\n

<\/h3>\n

<\/h3>\n

<\/h3>\n

<\/h3>\n

<\/h3>\n

python code for Fig. 1.<\/em><\/p>\n

#---------------------------------------------------------------------------------------------------<\/p>\n

Listing 1<\/em><\/h4>\n

ipython --pylab
\nfrom pylab import *
\nfrom matplotlib.font_manager import FontProperties
\nfrom matplotlib.backends.backend_agg import RendererAgg<\/p>\n

w0 = 70. #\u00a0 ft\/s
\nz0= 0\u00a0\u00a0 #\u00a0 ft
\ng = 32.\u00a0 #\u00a0 ft\/s\/s
\na=-g
\ntmax=w0\/g
\ndef zr(t):
\nreturn z0 + w0 * t + 0.5 * a* t*t<\/p>\n

zmax=zr(tmax)
\nt=arange(0,2*tmax+.005, 0.01)
\nplot(t,zr(t), color='k')
\nplot([0,4.5],[0,0],color='r')
\nplot([tmax,tmax],[0,zmax], color='b')
\ntitle('Projectile shot straight up', color='k')
\nxlabel('Time (s)', color='k')
\nylabel('z (ft)', color='k')
\n#--------------------------------------------------------------------------------------------------------<\/p>\n

Ballistics<\/em><\/strong><\/h2>\n

If a projectile is launched from the surface with a vertical component of velocity, w0<\/sub>,\u00a0 and horizontal component u0<\/sub>, then the vertical part of the motion is described by the above equations, but the horizontal motion is described independently as unaccelerated motion as follows:<\/p>\n

x = x0<\/sub>\u00a0 + u0<\/sub>*t. \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (10a)<\/p>\n

Label the coordinates as zero at the launch point, in which case x0<\/sub>=0, and<\/p>\n

x = u0<\/sub>*t. \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (10b)<\/p>\n

This (10b) represents a rescaling of the horizontal axis. You can compute x and z for any value of t during the projectile's flight, and plot a map of the path of the projectile's flight in physical space. Or, you can eliminate t from the equations of motion in favor of x as an independent variable. Fig. 2 displays the first alternative, which is plotted by the additional python code given below.<\/p>\n

 <\/p>\n

\"\"<\/a>

Fig. 2 A plot of a projectile motion, when fired at an angle, having an upward velocity component of 70. ft\/s and a horizontal one of 30. ft\/s.<\/p><\/div>\n

 <\/p>\n

 <\/p>\n

 <\/p>\n

 <\/p>\n

 <\/p>\n

 <\/p>\n

 <\/p>\n

 <\/p>\n

 <\/p>\n

python code (follows listing 1)<\/em><\/h3>\n

listing 2<\/em><\/h4>\n

#Horizontal component of initial velocity
\nu0 = 30\u00a0\u00a0 #ft\/s
\ndef xr(t):
\nreturn\u00a0 u0 * t\u00a0\u00a0 # horizontal distance traveled
\nfigure()
\nplot(xr(t),zr(t), color='k')
\nplot([0,140.],[0.,0.], color='r')
\nplot ([xr(tmax),xr(tmax)],[0,zmax], color='b')
\ntext(xr(tmax),2.5,r'65.6 ft')
\ntext(xr(tmax),70,r'76.6 ft')
\ntitle('Projectile shot at angle', color='k')
\nxlabel('x (ft)', color='k')
\nylabel('z (ft)', color='k')<\/p>\n

#-------------------------------------------------END----------------------------------------------------<\/p>\n","protected":false},"excerpt":{"rendered":"

\u00a9Fernando Caracena 31 August 2012 Independent motions One of Galileo's great discoveries is that motion in three dimensions can be decomposed into three separate motions along each of three mutually orthogonal directions. Specifically, up and down motions are governed by … Continue reading →<\/span><\/a><\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[33,10,12,34,30],"tags":[39,18,40],"class_list":["post-250","post","type-post","status-publish","format-standard","hentry","category-frame-of-reference","category-mathematics","category-physics","category-relativity","category-vectors","tag-ballistics","tag-galileo","tag-projectile"],"_links":{"self":[{"href":"http:\/\/www.ghyzmo.com\/wp-json\/wp\/v2\/posts\/250","targetHints":{"allow":["GET"]}}],"collection":[{"href":"http:\/\/www.ghyzmo.com\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"http:\/\/www.ghyzmo.com\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"http:\/\/www.ghyzmo.com\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"http:\/\/www.ghyzmo.com\/wp-json\/wp\/v2\/comments?post=250"}],"version-history":[{"count":18,"href":"http:\/\/www.ghyzmo.com\/wp-json\/wp\/v2\/posts\/250\/revisions"}],"predecessor-version":[{"id":287,"href":"http:\/\/www.ghyzmo.com\/wp-json\/wp\/v2\/posts\/250\/revisions\/287"}],"wp:attachment":[{"href":"http:\/\/www.ghyzmo.com\/wp-json\/wp\/v2\/media?parent=250"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"http:\/\/www.ghyzmo.com\/wp-json\/wp\/v2\/categories?post=250"},{"taxonomy":"post_tag","embeddable":true,"href":"http:\/\/www.ghyzmo.com\/wp-json\/wp\/v2\/tags?post=250"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}