Projectile & Force Motion — Projectile Motion and Motion Under Forces (Modeling) - Calculus 3

SECTION 1

Historical Context & Motivation

The study of projectile motion is one of the oldest and most consequential threads in mathematical physics, tracing its roots to ancient Greece and reaching maturity in the hands of the great natural philosophers of the Renaissance and the Enlightenment. While early thinkers debated whether a thrown stone followed a straight line before falling abruptly, the development of vector calculus eventually provided the definitive language for describing curved paths under the influence of forces. Understanding how position, velocity, and acceleration relate through differentiation and integration of vector-valued functions transformed ballistics, engineering, and orbital mechanics into exact sciences.

1638 • Galileo's Parabolic Trajectory

In Discorsi e dimostrazioni matematiche, Galileo demonstrated that projectile motion decomposes into uniform horizontal motion and uniformly accelerated vertical motion, producing a parabolic path — a radical departure from Aristotelian physics.

1687 • Newton's Principia Mathematica

Isaac Newton unified terrestrial and celestial mechanics through his three laws of motion and the law of universal gravitation, showing that forces determine acceleration and, through integration, velocity and position.

1740s • Euler's Analytical Mechanics

Leonhard Euler reformulated Newtonian mechanics in the language of differential equations and introduced component-wise analysis, laying the groundwork for the modern treatment of vector-valued position functions r(t).

1880s • Gibbs & Heaviside — Vector Notation

Josiah Willard Gibbs and Oliver Heaviside developed the compact vector notation (i, j, k basis) that allowed Newton's second law to be written as F = ma with each symbol representing a vector quantity, streamlining multivariable motion analysis.

20th c. • Modern Computational Modeling

Numerical methods — Euler's method, Runge–Kutta schemes — extended vector-valued motion modeling to situations with nonlinear drag, variable gravity, and electromagnetic forces where closed-form solutions are unavailable.

The central question this lesson addresses is: given a force field F(t) acting on an object of mass m, how do we construct the vector-valued position function r(t) that describes the object's trajectory, and how do initial conditions uniquely determine that path? This is the modeling problem at the heart of Calculus 3 motion analysis.

SECTION 2

Core Principles & Definitions

Modeling projectile motion and motion under forces in Calculus 3 rests on expressing Newton's second law in vector form and then using the calculus of vector-valued functions — differentiation and integration — to move between acceleration, velocity, and position. The following foundational ideas underpin every calculation in this domain.

Position as a Vector-Valued Function

The position of a particle at time t is r(t) = x(t) i + y(t) j + z(t) k. Each component is a real-valued function of time; together they trace a curve in ℝ² or ℝ³.

Velocity and Acceleration via Derivatives

Velocity v(t) = r′(t) and acceleration a(t) = r″(t). Differentiation acts component-wise, so each coordinate is handled independently once the vector framework is established.

Newton's Second Law in Vector Form

F(t) = m a(t) = m r″(t). Given a force field, the acceleration is a(t) = F(t)/m. This is the vector ODE from which velocity and position are recovered by integration.

Initial Conditions Fix the Solution

Two successive integrations introduce two vector constants of integration. The initial position r(0) = r₀ and initial velocity v(0) = v₀ uniquely determine the trajectory — this is the existence and uniqueness guarantee for second-order linear ODEs.

Superposition of Components

When force components along orthogonal axes are independent (e.g., gravity acts only vertically), the motion decouples into separate scalar problems. This superposition principle is what makes the ideal projectile problem tractable.

✦

KEY TAKEAWAY

Think of the force law F = m a as a recipe and integration as the cooking process. The recipe (acceleration) tells you the rate at which velocity changes. Integrating once gives you the velocity 'batter,' and integrating again gives you the finished 'cake' — the position function. The initial conditions are the starting ingredients: without knowing where and how fast the object starts, you cannot produce a unique trajectory, just as omitting a key ingredient leaves the recipe incomplete.

SECTION 3

Visual Explanation — Anatomy of a Projectile Trajectory

A projectile launched from r₀ with initial velocity v₀ at angle θ follows a parabolic arc (violet curve). The horizontal component of velocity (cyan, dashed) remains constant, while the vertical component (pink) changes linearly due to gravitational acceleration −g j. At the apex, the vertical velocity is zero; the resultant velocity vector (amber) is tangent to the curve at every point.

The diagram above captures the essential geometry of ideal projectile motion. The violet parabola is the locus of the position vector r(t). At every instant the velocity v(t) = r′(t) is tangent to the path, and the acceleration a(t) = −g j is constant and directed downward. Because no horizontal force acts on the projectile, the x-component of velocity is constant; the curve's curvature is entirely due to the changing y-component. This decomposition into independent horizontal and vertical motions is the key simplification that Galileo first articulated and that vector calculus formalizes through component-wise integration.

SECTION 4

Mathematical Framework

The entire modeling pipeline can be summarized in three stages: state the acceleration from Newton's second law, integrate once to obtain velocity (applying the initial velocity condition), and integrate again to obtain position (applying the initial position condition). We develop each equation below for the general case and then specialize to the ideal projectile.

NEWTON'S SECOND LAW (VECTOR FORM)

a(t) = r″(t) = F(t) / m

Here r(t) is the position vector, F(t) is the net force, and m is the mass. For a projectile near Earth's surface with no drag, F = −mg j, so a(t) = −g j, a constant vector.

VELOCITY BY INTEGRATION

v(t) = ∫ a(t) dt = a(t) · t + C₁ ⟹ v(t) = v₀ + a · t

The constant of integration C₁ is determined by the initial velocity: v(0) = v₀, so C₁ = v₀. For constant acceleration this yields a linear function of t. For time-varying a(t), the integral may be more complex.

POSITION BY INTEGRATION

r(t) = ∫ v(t) dt = r₀ + v₀ t + ½ a t²

Applying r(0) = r₀ fixes the second constant C₂ = r₀. This is the fundamental kinematic equation in vector form — valid whenever the acceleration is constant.

COMPONENT FORM — IDEAL PROJECTILE

x(t) = x₀ + (v₀ cos θ) t , y(t) = y₀ + (v₀ sin θ) t − ½ g t²

Here θ is the launch angle above the horizontal, v₀ = |v₀| is the initial speed, and g ≈ 9.8 m/s². Eliminating t between the two equations yields the parabola y = y₀ + (tan θ)(x − x₀) − g(x − x₀)² / (2v₀² cos² θ).

ℹ

Generalization: Non-Constant Forces

When the force depends on time, position, or velocity — for instance, a drag force proportional to v(t) — the equation r″(t) = F(t, r, r′)/m becomes a system of ODEs that generally requires numerical integration (e.g., Runge–Kutta methods). However, the modeling philosophy remains the same: specify F, derive a, integrate twice, apply initial conditions.

SECTION 5

Detailed Breakdown — Common Force Models

The ideal projectile (gravity only, no air resistance) is merely the simplest member of a family of force models that arise in Calculus 3 applications. The table below catalogues the most common scenarios, their force expressions, and the nature of the resulting position function. Understanding these models equips you to set up the correct initial-value problem for a given physical situation.

Force Model	Net Force F(t)	Closed-Form r(t)?
Ideal Projectile	−mg j	Yes — quadratic in t (parabola)
Linear Drag	−mg j − β v(t)	Yes — exponential decay terms
Quadratic Drag	−mg j − c \|v\|v	No — requires numerical methods
Uniform Electric Field	qE (constant vector)	Yes — parabola (analogous to gravity)
Lorentz (Magnetic) Force	q v(t) × B	Yes — circular / helical motion

Common force models and their tractability

Three trajectories launched with identical v₀ and θ. The ideal parabola (solid violet) achieves the greatest range. Linear drag (dashed cyan) shortens the range and lowers the apex. Quadratic drag (dotted pink) produces the most dramatic reduction, with a steeper descent than ascent — a hallmark of realistic ballistic trajectories.

Several features of the comparison are worth noting. First, the ideal trajectory is perfectly symmetric about its midpoint — a consequence of the fact that x(t) is linear and y(t) is quadratic. Once drag is introduced, the ascending branch is longer (in distance) than the descending branch because the object spends more time near the top where its speed is lowest and drag is weakest. Second, the terminal velocity concept becomes relevant for drag models: as the object falls, it may approach a speed at which drag balances gravity, limiting the impact speed. In Calculus 3 we model these effects by replacing the constant acceleration vector with a velocity-dependent one and solving the resulting differential equation system.

SECTION 6

Worked Example — Finding the Position Function and Range

A ball is launched from the origin with initial speed v₀ = 50 m/s at an angle θ = 30° above the horizontal in a uniform gravitational field (g = 9.8 m/s²). Air resistance is negligible. Determine the vector-valued position function r(t), the time of flight, the maximum height, and the horizontal range.

Ideal Projectile — Full Trajectory Analysis

1
Step 1 — State the Acceleration

The only force is gravity: F = −mg j. Dividing by m gives the acceleration vector a(t) = −g j = −9.8 j m/s². Note that there is no i-component.

a(t) = ⟨0, −9.8⟩ m/s²
2
Step 2 — Integrate to Find Velocity

v(t) = ∫ a dt = ⟨0, −9.8t⟩ + C₁. Setting t = 0: v(0) = v₀ = ⟨50 cos 30°, 50 sin 30°⟩ = ⟨25√3, 25⟩ ≈ ⟨43.30, 25⟩ m/s, so C₁ = v₀.

v(t) = ⟨25√3, 25 − 9.8t⟩ m/s
3
Step 3 — Integrate to Find Position

r(t) = ∫ v dt = ⟨25√3 · t, 25t − 4.9t²⟩ + C₂. With r(0) = ⟨0, 0⟩, we have C₂ = 0.

r(t) = ⟨25√3 · t, 25t − 4.9t²⟩ m
4
Step 4 — Time of Flight

The projectile lands when y(t) = 0 (and t > 0): 25t − 4.9t² = 0 ⟹ t(25 − 4.9t) = 0 ⟹ t = 25/4.9 ≈ 5.102 s.

T ≈ 5.10 s
5
Step 5 — Maximum Height

At the apex, $v_y$ = 0: 25 − 9.8t = 0 ⟹ t* = 25/9.8 ≈ 2.551 s. Then y(t*) = 25(2.551) − 4.9(2.551)² ≈ 63.78 − 31.87 ≈ 31.89 m. Alternatively, y_max = v₀² sin²θ / (2g) = 625 × 0.25 / 19.6 ≈ 31.89 m.

y_max ≈ 31.89 m
6
Step 6 — Horizontal Range

R = x(T) = 25√3 × 5.102 ≈ 43.30 × 5.102 ≈ 220.9 m. The closed-form expression is R = v₀² sin 2θ / g = 2500 × sin 60° / 9.8 = 2500 × 0.8660 / 9.8 ≈ 220.9 m.

R ≈ 220.9 m

SECTION 7

Strengths and Limitations of the Ideal Model

Every mathematical model involves trade-offs between analytical tractability and physical fidelity. The ideal projectile model, with its constant-acceleration assumption, is a cornerstone of introductory mechanics precisely because it yields clean closed-form solutions that illuminate fundamental relationships. However, real-world projectiles — baseballs, rockets, raindrops — operate in regimes where the model's assumptions break down.

Feature	Strength	Limitation
Analytical solutions	r(t) is a polynomial — easy to differentiate, integrate, optimize	Breaks down when force depends on v or r (drag, variable gravity)
Decoupling of axes	Horizontal and vertical motions solved independently	Drag couples the components, destroying this simplification
Symmetric trajectory	Time of flight = 2 × time to apex; simple range formula	Real trajectories are asymmetric — steeper descent, shorter range
Flat-Earth assumption	Excellent for short-range problems (lab bench, sports)	Invalid for long-range or orbital trajectories (curvature of Earth, variation of g)
Pedagogical value	Perfect vehicle for learning vector integration and initial-value problems	May instill false confidence if students mistake the model for reality

Strengths and limitations of the ideal (no-drag, constant-g) projectile model

✦

KEY TAKEAWAY

The ideal projectile model is like a frictionless ramp in engineering coursework — it strips away complication so the core mechanism (vector integration with initial conditions) stands in sharp relief. Once you have mastered the constant-acceleration pipeline, you can layer in drag, Coriolis effects, or variable gravity as perturbations, often requiring numerical solvers but the same foundational approach.

SECTION 8

Connections to Advanced Theory

The projectile motion framework is not merely an end in itself; it serves as a gateway to several deeper ideas that appear later in multivariable calculus, differential equations, and physics. Recognizing these connections early helps you see how the relatively simple act of integrating a constant acceleration vector generalizes into far richer mathematical structures.

This Lesson	Advanced Generalization
Constant acceleration a = −g j	Position-dependent acceleration a(r) in orbital mechanics (inverse-square law)
r″(t) = F/m with initial conditions	Systems of ODEs / phase-space analysis in dynamical systems theory
Component-wise integration in ℝ²	Integration along curves in ℝ³, line integrals of vector fields
Speed \|v(t)\| and arc length	Curvature κ, torsion τ, and the Frenet–Serret frame (TNB frame)
Range optimization (maximize R over θ)	Calculus of variations and optimal control theory

From ideal projectiles to advanced theory

Perhaps the most immediate extension is the tangential and normal decomposition of acceleration. In this lesson we decomposed a into i and j components, but in later sections of Calculus 3 you will learn to decompose a into a tangential component aT (changing speed) and a normal component aN (changing direction). For the ideal projectile, both components are nonzero except at the apex, where aT = 0 and all acceleration is normal, instantaneously redirecting the velocity vector from upward to downward.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

Explain why, in the ideal projectile model (no drag), the horizontal component of velocity remains constant throughout the flight. Reference Newton's second law in your explanation.

Show answer

Newton's second law states that acceleration equals net force divided by mass: a = F/m. In the ideal model, the only force is gravity, F = −mg j, which has no i-component. Therefore

a_x

= 0, meaning

v_x

(t) = v₀ cos θ for all t. A zero net force in a given direction implies zero acceleration in that direction, and hence constant velocity in that direction. This is a direct consequence of Newton's first law applied to the horizontal coordinate.

Question 2: BASIC CALCULATION

A projectile is launched from the origin with initial velocity v₀ = ⟨20, 30⟩ m/s under constant gravitational acceleration a = ⟨0, −9.8⟩ m/s². Find the vector-valued position function r(t) and determine the time at which the projectile returns to ground level (y = 0).

Show answer

Integrate a to get v(t) = ⟨20, 30 − 9.8t⟩ (using v(0) = ⟨20, 30⟩). Integrate v to get r(t) = ⟨20t, 30t − 4.9t²⟩ (using r(0) = ⟨0, 0⟩). Setting y(t) = 0: 30t − 4.9t² = 0 → t(30 − 4.9t) = 0, so t = 0 or t = 30/4.9 ≈ 6.122 s. The projectile returns to the ground at approximately t ≈ 6.12 s, at horizontal position x ≈ 20 × 6.122 ≈ 122.4 m.

Question 3: INTERMEDIATE

A particle moves under the acceleration a(t) = ⟨6t, −2, cos t⟩ with initial conditions r(0) = ⟨1, 0, 0⟩ and v(0) = ⟨0, 3, 1⟩. Find the position function r(t).

Show answer

Integrate a component-wise: v(t) = ⟨3t², −2t, sin t⟩ + C₁. Applying v(0) = ⟨0, 3, 1⟩ gives C₁ = ⟨0, 3, 1⟩, so v(t) = ⟨3t², 3 − 2t, 1 + sin t⟩. Integrate again: r(t) = ⟨t³, 3t − t², t − cos t⟩ + C₂. Applying r(0) = ⟨1, 0, 0⟩: the z-component at t = 0 is 0 − cos 0 + C₂z = 0 − 1 + C₂z, so C₂z = 1; C₂x = 1; C₂y = 0. Therefore r(t) = ⟨t³ + 1, 3t − t², t − cos t + 1⟩. You can verify by differentiating twice to recover a(t).

Question 4: APPLIED

A drone drops a payload from a height of 200 m while flying horizontally at 15 m/s. The payload experiences a linear drag force F_drag = −0.5v (in addition to gravity, g = 9.8 m/s²), and the payload's mass is 2 kg. Write the system of differential equations governing x(t) and y(t), identify the initial conditions, and find the velocity function v(t) by solving the resulting first-order linear ODE for each component.

Show answer

The net force is F = ⟨−0.5 vₓ, −2(9.8) − 0.5 vᵧ⟩. Dividing by m = 2: a = ⟨−0.25 vₓ, −9.8 − 0.25 vᵧ⟩. This gives two first-order linear ODEs: vₓ′ = −0.25 vₓ and vᵧ′ = −9.8 − 0.25 vᵧ. Initial conditions: r(0) = ⟨0, 200⟩, v(0) = ⟨15, 0⟩ (dropped, so no initial vertical velocity). Solving vₓ′ + 0.25 vₓ = 0 yields vₓ(t) = 15

e^{−0.25t}

. Solving vᵧ′ + 0.25 vᵧ = −9.8 yields vᵧ(t) = −39.2(1 −

e^{−0.25t}

). The terminal vertical speed is 39.2 m/s, and the horizontal speed decays exponentially. Position is obtained by integrating each component with the given initial positions.

Question 5: CRITICAL THINKING

For the ideal projectile launched from the origin at speed v₀ and angle θ, the range is R(θ) = v₀² sin 2θ / g. (a) Prove this formula starting from the position function. (b) Use calculus to show that R is maximized at θ = π/4. (c) Discuss whether the conclusion θ* = π/4 survives the introduction of a linear drag force F_drag = −β v. Provide a qualitative argument, not a full derivation.

Show answer

(a) From r(t) = ⟨(v₀ cos θ)t, (v₀ sin θ)t − ½gt²⟩, the flight time is T = 2v₀ sin θ / g (set y = 0, t > 0). Range R = x(T) = v₀ cos θ × 2v₀ sin θ / g = v₀²(2 sin θ cos θ)/g = v₀² sin 2θ / g. (b) dR/dθ = v₀²(2 cos 2θ)/g = 0 ⟹ cos 2θ = 0 ⟹ 2θ = π/2 ⟹ θ = π/4. Since d²R/dθ² = −4v₀² sin 2θ / g < 0 at θ = π/4, this is a maximum. (c) With linear drag, the horizontal velocity decays exponentially, meaning the projectile spends less effective time traveling horizontally. To compensate, a lower launch angle (θ* < π/4) is optimal because it increases the initial horizontal speed component at the expense of height, reducing the time the object is in the air and subject to drag. The exact optimal angle depends on the drag coefficient β and can only be found numerically in general.

SUMMARY

Lesson Summary

Summary

Modeling projectile and force-driven motion in Calculus 3 centers on the vector-valued position function r(t) and the pipeline that derives it: start from Newton's second law a(t) = F(t)/m, integrate once to obtain velocity v(t) (fixing the constant via the initial velocity v₀), and integrate again to obtain position r(t) (fixing the constant via the initial position r₀). For ideal projectile motion (constant gravitational acceleration, no drag), this yields a parabolic trajectory with component equations x(t) = (v₀ cos θ)t and y(t) = (v₀ sin θ)t − ½gt², from which the time of flight, maximum height, and range are easily extracted.

Beyond the ideal case, the same modeling philosophy extends to velocity-dependent forces (drag), position-dependent forces (gravity, electric fields), and cross-product forces (magnetism). In these more complex scenarios the resulting differential equations may require numerical integration, but the underlying structure — decompose into components, apply Newton's law, integrate with initial conditions — remains invariant. Mastery of the ideal projectile problem thus equips you with the conceptual template for all of classical particle mechanics.