AP Physics C: Mechanics Mastery Guide
Work-Energy Theorem with Velocity-Dependent Drag Forces
1. Introduction & AP Exam Weight
The integration of velocity-dependent non-conservative forces into the Work-Energy Theorem represents one of the most mathematically demanding topics in AP Physics C: Mechanics. While basic energy conservation problems appear routinely, velocity-dependent drag forces—governed by linear ($F_d = -bv$) or quadratic ($F_d = -cv^2$) models—frequently anchor the calculus-heavy Free-Response Questions (FRQs), typically accounting for 14% to 20% of total exam points in Free Response Section testing.
Conceptual Scope
Standard Newtonian mechanics treats work as a spatial line integral of a position-dependent force:
$$W = \int_{\vec{r}_1}^{\vec{r}_2} \vec{F} \cdot d\vec{r}$$
When force depends explicitly on velocity $\vec{F}(\vec{v})$, work becomes inherently path- and time-dependent. To solve these systems, you must bridge differential equations with dynamic limits of integration, converting time derivatives into spatial derivatives using the chain rule identity:
$$a = \frac{dv}{dt} = v\frac{dv}{dx}$$
UC Berkeley Context
For prospective STEM majors targeting the University of California, Berkeley, securing a Score 5 on the AP Physics C: Mechanics exam satisfies the prerequisite for Physics 7A (4 semester units), placing you directly into Physics 7B (Heat, Electricity, and Magnetism).
Mastery of velocity-dependent mechanics demonstrates the mathematical fluency in differential equations and path integrals required by Berkeley’s College of Engineering (EECS, Mechanical, Civil) and the College of College of Chemistry.
2. Deep Concept Breakdown
Differential Formulation of the Work-Energy Theorem
The Work-Energy Theorem states that the net work done on a particle equals its change in kinetic energy:
$$W_{\text{net}} = \Delta K = K_f - K_i = \frac{1}{2}m v_f^2 - \frac{1}{2}m v_i^2$$
When non-conservative forces such as velocity-dependent drag are present:
$$W_{\text{net}} = W_{\text{cons}} + W_{\text{nc}} = -\Delta U + W_{\text{drag}}$$
$$\implies \Delta K + \Delta U = W_{\text{drag}} = \int_{x_0}^{x_f} F_d(v) \, dx$$
Because $F_d(v)$ depends on velocity rather than explicit position $x$, the spatial integral $\int F_d(v) \, dx$ cannot be integrated directly without expressing $v$ as a function of $x$.
Derivation 1: Spatial Velocity Profile via Spatial Acceleration Substitution
Consider a particle of mass $m$ moving horizontally subject only to a linear drag force $F_d(v) = -bv$.
1. Set up Newton's Second Law:
$$m a = -bv$$
2. Apply the Differential Chain Rule:
Substitute $a = v \frac{dv}{dx}$:
$$m v \frac{dv}{dx} = -bv$$
3. Separate Variables:
Assuming $v \neq 0$:
$$m \, dv = -b \, dx$$
4. Integrate with Definite Limits:
Integrate from initial state $(x_0 = 0, v_0)$ to final state $(x, v(x))$:
$$\int_{v_0}^{v(x)} m \, dv' = \int_{0}^{x} -b \, dx'$$
$$m (v(x) - v_0) = -bx$$
$$v(x) = v_0 - \frac{b}{m}x$$
Derivation 2: Work Done by Linear Drag Force
We evaluate $W_{\text{drag}}$ using two complementary methodologies to verify mathematical consistency.
Method A: Direct Spatial Integration
$$W_{\text{drag}} = \int_{0}^{x} F_d(v) \, dx' = \int_{0}^{x} -b \left(v_0 - \frac{b}{m}x'\right) dx'$$
$$W_{\text{drag}} = -b \left[ v_0 x' - \frac{b}{2m}(x')^2 \right]_{0}^{x} = -bv_0 x + \frac{b^2}{2m}x^2$$
Method B: Kinetic Energy Verification
$$\Delta K = \frac{1}{2}m (v(x))^2 - \frac{1}{2}m v_0^2$$
$$\Delta K = \frac{1}{2}m \left(v_0 - \frac{b}{m}x\right)^2 - \frac{1}{2}m v_0^2$$
$$\Delta K = \frac{1}{2}m \left(v_0^2 - \frac{2bv_0 x}{m} + \frac{b^2 x^2}{m^2}\right) - \frac{1}{2}m v_0^2$$
$$\Delta K = -bv_0 x + \frac{b^2}{2m}x^2$$
Both methods yield identical results, proving that the non-conservative work done by linear drag is explicitly path-dependent and bounded by the terminal stopping distance $x_{\text{stop}} = \frac{m v_0}{b}$.
Computational Simulation: Linear vs. Quadratic Drag Energetics
The following Python script computes state variable trajectories, kinetic energy dissipation, and cumulative non-conservative work performed under both linear ($F_d = -bv$) and quadratic ($F_d = -cv^2$) drag regimes.
import numpy as np
def simulate_drag_energetics(
m: float,
v0: float,
b: float,
c: float,
dt: float = 0.001,
t_max: float = 10.0
):
"""
Simulates horizontal motion and energetic dynamics under Linear vs. Quadratic drag.
Parameters:
m (float): Mass of the object (kg)
v0 (float): Initial velocity (m/s)
b (float): Linear drag coefficient (kg/s)
c (float): Quadratic drag coefficient (kg/m)
dt (float): Time step for numerical integration (s)
t_max (float): Total simulation time (s)
"""
steps = int(t_max / dt)
# Pre-allocate arrays for Linear Drag simulation
t = np.linspace(0, t_max, steps)
x_lin = np.zeros(steps)
v_lin = np.zeros(steps)
W_drag_lin = np.zeros(steps)
# Pre-allocate arrays for Quadratic Drag simulation
x_quad = np.zeros(steps)
v_quad = np.zeros(steps)
W_drag_quad = np.zeros(steps)
# Initial Conditions
v_lin[0] = v0
v_quad[0] = v0
# Euler-Cromer Integration Loop
for i in range(1, steps):
# --- Linear Drag Dynamics ---
F_lin = -b * v_lin[i-1]
a_lin = F_lin / m
v_lin[i] = v_lin[i-1] + a_lin * dt
x_lin[i] = x_lin[i-1] + v_lin[i] * dt
# Incremental work: dW = F * dx
W_drag_lin[i] = W_drag_lin[i-1] + F_lin * (v_lin[i] * dt)
# --- Quadratic Drag Dynamics ---
F_quad = -c * (v_quad[i-1]**2)
a_quad = F_quad / m
v_quad[i] = v_quad[i-1] + a_quad * dt
x_quad[i] = x_quad[i-1] + v_quad[i] * dt
W_drag_quad[i] = W_drag_quad[i-1] + F_quad * (v_quad[i] * dt)
# Compute Kinetic Energy Changes
K_initial = 0.5 * m * (v0**2)
K_final_lin = 0.5 * m * (v_lin[-1]**2)
K_final_quad = 0.5 * m * (v_quad[-1]**2)
delta_K_lin = K_final_lin - K_initial
delta_K_quad = K_final_quad - K_initial
print(f"=== Energetics Summary at t = {t_max}s ===")
print(f"Linear Drag : Delta K = {delta_K_lin:.4f} J | W_drag = {W_drag_lin[-1]:.4f} J")
print(f"Quadratic Drag : Delta K = {delta_K_quad:.4f} J | W_drag = {W_drag_quad[-1]:.4f} J")
return {
"t": t,
"x_lin": x_lin, "v_lin": v_lin, "W_lin": W_drag_lin,
"x_quad": x_quad, "v_quad": v_quad, "W_quad": W_drag_quad
}
if __name__ == "__main__":
# Test Parameters: m = 1.5 kg, v0 = 20.0 m/s, b = 0.5 kg/s, c = 0.05 kg/m
results = simulate_drag_energetics(m=1.5, v0=20.0, b=0.5, c=0.05)
3. Common AP Exam Pitfalls & Score 5 Scoring Rubric Nuances
Contrast: Score 4 vs. Score 5 Performance
| Problem Element | Score 4 Student Approach | Score 5 Student Approach |
|---|---|---|
| Work Integral Setup | Attempts to calculate $W = F \cdot d$ using algebraic multiplication ($W = -bv \cdot x$), incorrectly assuming constant force over distance. | Defines work through differential integration: $W = \int F(v) \, dx$, applying $dx = v \, dt$ or $dx = \frac{v}{a} \, dv$ to change variables. |
| Variable Substitution | Gets stuck when integrating $\int -bv \, dx$ because $v$ varies with $x$. Often treats $v$ as a constant $v_{\text{avg}}$. | Replaces acceleration $a$ with $v \frac{dv}{dx}$ in Newton's 2nd Law to immediately establish an integrable differential equation relating $v$ and $x$. |
| Integration Limits | Uses indefinite integrals and forgets the constant of integration $C$, or applies incorrect, non-matching limits on opposite sides of equation. | Consistently applies matching definite limits: $\int_{v_0}^{v(x)} m \, dv' = \int_{0}^{x} -b \, dx'$, evaluating boundary conditions cleanly. |
| Sign & Vector Errors | Fails to account for direction in dot products, writing $\vec{F}_d \cdot d\vec{r} = +bv \, dx$, leading to unphysical kinetic energy gains. | Correctly models resistive forces as negative work ($\vec{F}_d \cdot d\vec{r} = -F_d \, dx$), verifying that mechanical energy strictly dissipates. |
Scoring Rubric Criteria Breakdown (AP Free-Response Standards)
When AP Readers grade non-conservative work problems with differential drag, points are allocated according to strict, objective criteria:
- Newton’s Second Law Setup (+1 point): Writing $\Sigma F = m a$ with the velocity-dependent term correctly expressed ($m \frac{dv}{dt}$ or $m v \frac{dv}{dx}$).
- Chain Rule / Variable Separation (+1 point): Correctly substituting $a = v \frac{dv}{dx}$ and separating variables such that all terms containing $v$ are on one side and $x$ terms are on the other.
- Integration & Boundary Execution (+1 point): Demonstrating correct integration (e.g., $\int \frac{v}{a(v)} dv$) with proper definite limits corresponding to physical boundary conditions.
- Final Algebraic Isolation (+1 point): Algebraically isolating the requested variable ($v(x)$, $x_{\text{max}}$, or $W_{\text{drag}}$) without algebraic errors.
4. UC Berkeley Placement Pathway
Exemption Mechanics & Academic Credit
Achieving a Score 5 on the AP Physics C: Mechanics exam unlocks specific administrative and academic placement benefits at UC Berkeley:
AP Physics C: Mechanics (Score 5)
│
├─► Satisfies: Physics 7A Requirement
├─► Credits Awarded: 4 Units (UC Berkeley Transcript)
└─► Placement: Direct Entry into Physics 7B (Heat, Electricity & Magnetism)
Academic Impact Across Berkeley Majors
- Mechanical Engineering & Civil Engineering: Skipping Physics 7A accelerates your track toward core upper-division prerequisites, such as Engineering 7 (Numerical Methods) and ME 104 (Engineering Mechanics II / Dynamics).
- Electrical Engineering & Computer Sciences (EECS): Exemption clears fundamental physics requirements, allowing earlier entry into EECS 16A/16B (Designing Information Devices and Systems) or CS 61A.
- Physics & Astrophysics Majors: Accelerates entry into Physics 7C (Optics and Modern Physics) and Physics 105 (Analytic Mechanics), where drag force differential equations are revisited using Lagrangian and Hamiltonian formalisms.
Transitioning to Physics 7B: Mathematical Readiness
Physics 7B assumes complete mastery of: 1. Multivariable path integrals ($\int \vec{E} \cdot d\vec{\ell}$) 2. Non-conservative field dissipation (e.g., inductive drag, eddy currents) 3. Differential equations governing transient energy transfers (RC, RL, and RLC circuit dynamics)
Mastering drag force derivations using spatial substitutions ($v \frac{dv}{dx}$) prepares you for the dynamic equations encountered in Berkeley's Physics 7B curriculum.
5. High-Yield Practice Problem
Problem Statement
A projectile of mass $m$ is launched vertically upward from ground level ($y = 0$) with an initial launch speed $v_0$. As the projectile travels upward through the air, it experiences a quadratic drag force given by:
$$\vec{F}_d = -c v^2 \hat{j}$$
where $c$ is a positive constant and $v$ is the instantaneous speed. Acceleration due to gravity $g$ acts vertically downward.
▲ +y
│
│ [ Mass m ] ▲ v(y)
│ │
│ ├── F_drag = -c v^2
│ └── F_gravity = -m g
│
y = 0 ──┴───────────────────────
(a) Write the differential equation governing the upward motion of the projectile using Newton's Second Law in terms of $v$, $dy$, $dv$, $m$, $g$, and $c$.
(b) Derive an explicit expression for the maximum height $y_{\text{max}}$ attained by the projectile in terms of $m$, $g$, $c$, and $v_0$.
(c) Using the Work-Energy Theorem, derive an expression for the total work done by the drag force $W_{\text{drag}}$ during the upward flight from $y = 0$ to $y = y_{\text{max}}$.
(d) Derive an expression for the speed $v_f$ of the object when it returns to ground level ($y = 0$). Show that in the limit $c \to 0$, $v_f \to v_0$.
Step-by-Step Solution & Rubric Checklist
Part (a): Differential Equation Setup
-
Apply Newton's Second Law for the upward motion ($+y$ defined as upward): $$\Sigma F_y = -m g - c v^2 = m a$$
-
Replace $a$ with the spatial differential form $a = v \frac{dv}{dy}$: $$-m g - c v^2 = m v \frac{dv}{dy}$$
-
Rearrange into separable form: $$m v \frac{dv}{dy} = -(m g + c v^2)$$
Part (b): Derive Maximum Height $y_{\text{max}}$
-
Separate variables: $$\frac{m v}{m g + c v^2} \, dv = -dy$$
-
Set up definite integrals. At $y = 0$, $v = v_0$; at $y = y_{\text{max}}$, $v = 0$: $$\int_{v_0}^{0} \frac{m v}{m g + c v^2} \, dv = \int_{0}^{y_{\text{max}}} -dy$$
-
Reverse limits on the left side to eliminate the negative sign: $$\int_{0}^{v_0} \frac{m v}{m g + c v^2} \, dv = y_{\text{max}}$$
-
Solve using $u$-substitution ($u = m g + c v^2 \implies du = 2 c v \, dv \implies v \, dv = \frac{du}{2c}$): $$y_{\text{max}} = \frac{m}{2c} \int_{m g}^{m g + c v_0^2} \frac{du}{u} = \frac{m}{2c} \left[ \ln u \right]_{m g}^{m g + c v_0^2}$$
$$y_{\text{max}} = \frac{m}{2c} \ln \left( \frac{m g + c v_0^2}{m g} \right) = \frac{m}{2c} \ln \left( 1 + \frac{c v_0^2}{m g} \right)$$
Part (c): Work Done by Drag Force $W_{\text{drag}}$
-
Apply the Work-Energy Theorem for the upward path: $$W_{\text{net}} = \Delta K = K_f - K_i = 0 - \frac{1}{2} m v_0^2 = -\frac{1}{2} m v_0^2$$
-
Express net work as the sum of conservative and non-conservative components: $$W_{\text{net}} = W_{\text{gravity}} + W_{\text{drag}}$$
-
Calculate gravity's work over the displacement $y_{\text{max}}$: $$W_{\text{gravity}} = -m g y_{\text{max}} = -m g \left[ \frac{m}{2c} \ln \left( 1 + \frac{c v_0^2}{m g} \right) \right]$$
-
Substitute into the Work-Energy equation and solve for $W_{\text{drag}}$: $$-\frac{1}{2} m v_0^2 = -m g y_{\text{max}} + W_{\text{drag}}$$
$$W_{\text{drag}} = \frac{1}{2} m v_0^2 - m g y_{\text{max}} = -\frac{1}{2} m v_0^2 + \frac{m^2 g}{2c} \ln \left( 1 + \frac{c v_0^2}{m g} \right)$$
Part (d): Downward Motion, Return Speed $v_f$, and Asymptotic Limit
-
For downward motion, gravity acts downward ($-y$), while drag acts upward ($+y$): $$\Sigma F_{\text{down}} = m g - c v^2 = m v \frac{dv}{dy'}$$ (where $y'$ is measured downward from $y_{\text{max}}$, so $y' = 0 \to y_{\text{max}}$, starting from $v = 0$ to $v_f$).
-
Separate variables and integrate: $$\int_{0}^{v_f} \frac{m v}{m g - c v^2} \, dv = \int_{0}^{y_{\text{max}}} dy' = y_{\text{max}}$$
-
Substitute $u = m g - c v^2 \implies du = -2 c v \, dv$: $$-\frac{m}{2c} \ln \left( \frac{m g - c v_f^2}{m g} \right) = y_{\text{max}}$$
-
Equate to the expression for $y_{\text{max}}$ from Part (b): $$-\frac{m}{2c} \ln \left( 1 - \frac{c v_f^2}{m g} \right) = \frac{m}{2c} \ln \left( 1 + \frac{c v_0^2}{m g} \right)$$
$$\ln \left( 1 - \frac{c v_f^2}{m g} \right)^{-1} = \ln \left( 1 + \frac{c v_0^2}{m g} \right)$$
$$\frac{1}{1 - \frac{c v_f^2}{m g}} = 1 + \frac{c v_0^2}{m g}$$
- Solve algebraically for $v_f^2$: $$1 - \frac{c v_f^2}{m g} = \frac{1}{1 + \frac{c v_0^2}{m g}} = \frac{m g}{m g + c v_0^2}$$
$$\frac{c v_f^2}{m g} = 1 - \frac{m g}{m g + c v_0^2} = \frac{c v_0^2}{m g + c v_0^2}$$
$$v_f = \frac{v_0}{\sqrt{1 + \frac{c v_0^2}{m g}}}$$
- Evaluate Limit as $c \to 0$: $$\lim_{c \to 0} v_f = \lim_{c \to 0} \frac{v_0}{\sqrt{1 + \frac{c v_0^2}{m g}}} = \frac{v_0}{\sqrt{1 + 0}} = v_0$$
This verifies energy conservation in the limit of zero air resistance.
Official AP Scoring Rubric Checklist
Part (a)
[1 Point] Correct application of Newton's 2nd Law with -mg and -cv^2.
[1 Point] Correct substitution of spatial acceleration a = v(dv/dy).
Part (b)
[1 Point] Separation of variables with appropriate integral limits.
[1 Point] Execution of logarithmic integration step.
[1 Point] Correct expression for y_max.
Part (c)
[1 Point] Application of Work-Energy Theorem linking W_net to Delta K.
[1 Point] Correct final algebraic expression for W_drag.
Part (d)
[1 Point] Differential setup for downward flight with inverted drag sign.
[1 Point] Correct expression for v_f.
[1 Point] Validated limit proof showing v_f -> v_0 as c -> 0.
TOTAL: 10/10 Points (Score 5 Standard)