Bomb Drop Equations

Arial Bombing Techniques By Crystal Pepper and Chris Wilson March 30, 2009 Abstract In this article we will explore the

Views 120 Downloads 3 File size 1MB

Report DMCA / Copyright

DOWNLOAD FILE

Recommend stories

Citation preview

Arial Bombing Techniques By Crystal Pepper and Chris Wilson March 30, 2009 Abstract In this article we will explore the bombing techniques used by the United States Army during World War II in order to produce a maximum impact on targets. We will demonstrate through calculations of differential equations how the speed of the plane, linear distance the plane is from the target, and altitude of the plane determines the time that the bombs must be deployed.

Intro During World War II the role of the bombardier was critical in the accuracy of placing the free falling bombs on target. The conditions needed were: altitude, true airspeed, bomb ballistics, trail, actual time of fall, ground speed, and drift. The altitude, which is controlled by the pilot, determines the amount of time the bomb is free falling. The pilot also controls the true airspeed which gives the bomb its initial forward velocity and in turn affects the trail (horizontal distance that the bomb is behind the plane upon impact). The bomb ballistics such as the size, shape, and density determine the air resistance the bomb will encounter as it falls. The actual time of fall is the length of time the bomb is sustained in the air from the instant of release to the instant of impact. The actual time of fall is affected by the altitude and bomb ballistics. The ground speed is the speed of the plane in relation to the Earth’s surface. This affects the range of the bomb. Lastly, the drift is determined by the direction and velocity of the wind and affects the distance the bomb will travel downwind. With the help of a mechanical analog computer, the then top secret Norden bomb sight, the bombardier was able to input the parameters altitude, airspeed, bomb ballistics, and air density. With these parameters, the accuracy of the placement of the bombs was sufficient for the technology available at the time.

1

Projectile Motion

When the bomb leaves the airplane it will not fall in a straight line. From the moment the bomb leaves the plane its trajectory will vary in curvature. Projectile motion is two-dimensional motion. However, any case of two dimensional motion can be resolved into two cases of one dimensional motion: one along the x − axis and the other along the y-axis. We will initially study the one dimensional cases separately, without considering air resistance.

1

Figure 1: Vertical Motion

1.1

Vertical Motion

We now develop a mathematical model of the bomb’s trajectory. For convenience, we will think of our bomb as an object moving in the xy plane and make use of this idea in our graphical analyses. Figure 1 shows the force on an object with a vertical motion. It is always true for an object near the Earth’s surface that the acceleration of gravity is equal to the differential change in velocity in a given time interval. dvy = a = −g dt We know that gravity g is 9.8, m/s2 . In our coordinate system, with y upward, this acceleration is a = −gy = −9.8 m/s2

(1)

This is the acceleration along the y-axis. The velocity in the y direction can now be determined by taking the anti-derivative of −g with respect to time. In order to set up our velocity equation, two initial conditions will be required. Our parameters are the initial time t = 0, and the initial vertical velocity vy = 0. Z dy = − gdt dt vyf = −gt + c, where c is a constant of integration. By letting c = vyi and vyi = 0 then it follows that vyf = vyi − gy t vyf = 0 − gt

(2)

vyf = −gt Equation (2) can be used to solve for the velocity along the y-axis at any instant t. By the same process, we can determine the displacement along the y-axis at any instant t. Integrating the velocity with respect to time gives

2

Z y=

(vyi − gt)dt

1 y = vyi t − gy t2 2 (3) 1 2 y = 0 − gy t 2 1 2 y = − gt . 2 We have now described all of the components of the bomb’s vertical drop. We can now calculate the acceleration at any time, the velocity at any time, and the position at any time.

1.2

Horizontal Motion

In the previous section, we described the motion of the bomb from its acceleration to its final displacement by taking the anti-derivative of each equation. Determining the velocity and displacement of the bomb in the horizontal direction can be calculated in the same manner. First, we will need to establish an initial condition. Since the bomb is attached to a plane which is flying in a horizontal direction, it will have the same initial velocity as the plane. However, due to air resistance and gravitation, the bomb will fall in an arc like path towards the ground. In other words, there is no horizontal acceleration and the plane maintains its velocity from the time of release to impact. Therefore, the initial conditions are ax = 0 and vxi = vx . We will need to use ax to solve for our final displacement. Now the component of initial velocity along the x-axis can be written by taking the integral of the acceleration component Z dx = ax dt dt vxf = ax t + c (4) vxf = 0 + c vxf = c. The c represents an arbitrary constant which we can solve for by plugging in our initial values of vxi = vxf and ax = 0 vxi = ax t + c vxi = 0 + c vxi = c. Plugging this c value into equation (4) with the initial conditions, will give us a first order differential equation that describes the horizontal speed of our bomb. vxf = ax t + c vxf = ax t + vxi vxf = 0 + vxi vxf = vxi .

3

Figure 2: Horizontal Motion This is the component of the velocity along the x-axis at any instant t. This means that the horizontal component of velocity does not change throughout the projectile’s motion. Let vxf = vxi = vx

(5)

represent our final velocity in the horizontal direction. By incorporating a horizontal acceleration, and integrating once more, we will be able to determine the displacement of the bombs trajectory. Z xf = (ax t + vx )dt 1 ax t2 + vx t + c. 2 We will need to solve for c to determine the final displacement equation. Plugging in the initial values and letting t = 0 we will solve for the initial horizontal position xf =

1 a x t 2 + vx t + c 2 xi = 0 + 0 + c xi =

xi = c. With this c value the displacement along the x-axis at any instant t can be described by xf = ax t2 + vx + c x f = 0 + vx t + x i

(6)

x f = vx t + x i . We can now determine the velocity and displacement at any time during the flight of the bomb. With equation (6), we can plot the bomb’s trajectory on an xy coordinate system as shown in Figure 3.

2 2.1

The Real Forces of Projectile Motion Drag and the Vertical Motion

In sections 1.1 and 1.2 we showed you the idealized formulas for the motion of a projectile by neglecting the force of drag. However, in order to establish a more real world basis for our model, we will need to determine and incorporate this force. Drag is defined as the force that opposes the motion of an object in a fluid. This fluid can be anything from liquid to gas. As our bomb leaves 4

Figure 3: Horizontal and Vertical trajectory. the plane, the air that it moves through will provide this resistance. Drag R depends on the density of the air ρ, the velocity v, the air’s viscosity and compressibility D, and the reference area A. 1 DρAv 2 . (7) 2 To derive the drag equation is a tedious task that would take hours of computation. So for the sake of time we state that this equation is just the drag force without any further explanation. Since air resistance is present in all real applications, the motion of the vertical fall of the bomb will now have a drag force. Using Newtons law of force F = ma we will sum the forces in the y direction to derive our acceleration formula. R=

F = Fg − Fd . Knowing that F = ma, Fg = mg, and Fd = 12 DρAv, the equation looks like 1 ma = mg − DρAv 2 (8) DρAv . a=g− 2m Equation (8)includes all of the forces that are acting on a bomb as it falls to the ground. Before we attempt to integrate this equation, we will simplify the drag by letting k = DρA/2m. This will help keep our values straight. Therefore, with the new substitution, equation (8) becomes

5

a = g − kv (9) dv a= = g − kv. dt Equation (9) is the acceleration due to the velocity and drag. By integrating this function, we will get a formula for the velocity with respect to time. In order for us to differentiate equation (9), we will need to separate the variables as follows dv = g − kv dt dv = (g − kv) dt dv = dt g − kv Now that the function is set up, we will integrate Z Z dv = dt g − kv ln (g − kv) =t+C k ln (g − kv) = −k (t + C) eln(g−kv) = e−kt+kC g − kv = e−kt ekC Through the integration, we are left with an arbitrary constant ekc . To simplify this, let b = ekc . With the substitution we can now solve for the velocity as a function of time. g − kv = be−kt kv = g − be−kt b g v = − e−kt . k k However, this equation is not yet correct. We will need to substitute our initial values to solve for the constant b. Setting v = 0 at t = 0 produces g b − e−kt k k g b 0 = − e−k(0) k k g b 0= − k k g = b.

v=

Our calculations reveal that b is equal to g. This will make the velocity equation a little bit easier to look at.

6

g g − e−kt k k  g 1 − e−kt v= (10) k  g v= 1 − e−kt . k This formula tells us that at any time t we can evaluate the velocity of the bomb. We would like to note that as t → ∞, v → g/k, called the terminal velocity. Up to this point, we have formulated equations for the acceleration and velocity of the bomb in the vertical direction. Neither of these equations, however, will tell us where the bomb will land when a given height and time are known. We will integrate once more to find the displacement of the bomb v=

g k dy g = k Z dt g dy = k g y= k v=

1 − e−kt



 1 − e−kt Z  1 − e−kt dt   e−kt t+ + c. k

Yet again, we are left with a constant of integration c. Providing some initial conditions will eliminate the constant. With the initial condition y = y0 at t = 0 gives   g e0 y0 = 0+ +c k k   g 1 +c y0 = k k g y0 = 2 + c k g c = y0 − 2 . k Plugging in the c value will give the final equation for displacement   g e−kt g y= t+ + y0 − 2 k k k g g g (11) y = t + 2 e−kt + y0 − 2 k k k  g g y = y0 + t − 2 1 − e−kt . k k This final equation represents the general solution to our system in the vertical direction. We can plot this, as shown in figure 4, given any initial height of the plane. We have now established formulas for the acceleration, velocity, and position for our trajectory in the vertical direction. We will then follow up with the derivation for the equations with drag in the horizontal direction.

7

Figure 4: Vertical Motion with Drag

2.2

Drag and the Horizontal Motion

In section 2.1, we derived the the general solution to the vertical displacement by integration. In the same way, we can derive the general solution to the horizontal dislpacement. By taking the same steps of integration we can go through this process without too much explanation. Starting with the initial forces in the x direction F = −Fd .

(12)

We can see from equation (12) that the only force in the horizontal direction is the force of drag. Solving this equation will provide a formula for the velocity. For convienance, we will let k = DρA/2m for our drag force. dv = −kv dt Seperate variables then integrate to solve for the final velocity equation a=

8

dv dt dv Z v dv v Z dv v lnv

= −kv = −kdt Z = −kdt Z = −k dt

(13)

= −kt + c

v = e−kt+c v = c1 e−kt . The integration constant c1 can be evaluated with the initial conditions t = 0 and v = v0 . v = c1 e−kt v 0 = c1 e 0 v 0 = c1 . Then plugging this value back into the equation (13) will give v = v0 e−kt

(14)

This is the final eqaution for the velocity with drag. We will then integrate, just as we did with the vertical equation, to solve for the displacement. v = v0 e−kt dx = v0 e−kt dt Z Z dx = v0 e−kt dt Z x = v0 e−kt dt e−kt + c. −k Yet again we are left with the integral constant that can be solved with the initial conditions t = 0 and x = 0 x = v0

e−kt +c −k e0 0 = v0 +c −k v0 0= +c −k

x = v0

v0 =c k 9

Figure 5: Horizontal Motion with Drag Therefore x = v0

e−kt v0 + −k k

(15)  v0 1 − e−kt k Equation (15) is our final formula needed to solve for the discplacement of the bomb on the x-axis. Setting an altitude and initial speed of the plane, we can plot this trajectory as shown in figure 5. x=

3

Conclusion

With equations (11) and (15), y = y0 +

 g g t − 2 1 − e−kt k k

 v0 1 − e−kt k we can plot the trajectory of the bomb in an xy coordinate system. These equations can be used to solve the position of the bomb when dropped from any height and distance. Setting an altitude and velocity of the plane, we have plotted these equations as demonstrated in figure 6. x=

10

Figure 6: Horizontal and Vertical Motion with Drag

References [1] Probable Bomb Impact Areas (PJR) I04FL040501 [2] Research and Education Association. 1996 Differential Equations.

11