AP Physics C: Mechanics Mastery Guide
Unit 5: Simple Harmonic Motion & Calculus-Driven Oscillations
1. Introduction & AP Exam Weight
Simple Harmonic Motion (SHM) and continuous oscillations represent approximately 10%–14% of the multiple-choice questions on the AP Physics C: Mechanics exam, and appear as a major component of Free-Response Questions (FRQs) nearly every year.
Unlike AP Physics 1, which primarily evaluates algebraic formula application ($T = 2\pi\sqrt{\frac{m}{k}}$), AP Physics C evaluates SHM through the lens of differential calculus, rotational dynamics, and non-linear system perturbations. You are expected to derive oscillatory equations of motion directly from fundamental principles ($\sum \mathbf{F} = m\mathbf{a}$ or $\sum \boldsymbol{\tau} = I\boldsymbol{\alpha}$), establish boundary conditions, and solve second-order ordinary differential equations (ODEs).
Core Conceptual Scope
- The SHM Differential Criterion: Proving that a system executes SHM by demonstrating that its acceleration is directly proportional and opposite in direction to its displacement: $$\frac{d^2x}{dt^2} + \omega^2 x = 0$$
- Kinematics & Energy Interdependence: Solving $x(t)$, $v(t)$, and $a(t)$ given initial conditions $x(0) = x_0$ and $v(0) = v_0$, while integrating potential energy functions $U(x) = -\int F \, dx$ to analyze phase-space dynamics.
- Complex Systems:
- Physical Pendulums: Extended rigid bodies oscillating under gravity.
- Torsional Oscillators: Rotational elasticity dynamics ($\tau = -\kappa \theta$).
- Systems with Variable Springs/Masses: Equivalent spring constants ($k_{\text{eff}}$) for series and parallel configurations, as well as fluid/buoyancy-driven oscillations.
2. Deep Concept Breakdown
A. Strict Calculus Derivation of the SHM Equations
1. Translational Oscillator (Mass-Spring)
Consider a particle of mass $m$ attached to an ideal spring with spring constant $k$ on a frictionless horizontal surface. Applying Newton’s Second Law:
$$\sum F_x = m a_x \implies -k x = m \frac{d^2x}{dt^2}$$
Rearranging into standard second-order linear homogeneous differential equation form:
$$\frac{d^2x}{dt^2} + \left(\frac{k}{m}\right) x = 0$$
Define the angular frequency $\omega \equiv \sqrt{\frac{k}{m}}$. The differential equation becomes:
$$\frac{d^2x}{dt^2} + \omega^2 x = 0$$
To solve this equation, assume an Ansatz solution of the form $x(t) = C e^{rt}$. Substituting this into the ODE yields the characteristic equation:
$$r^2 e^{rt} + \omega^2 e^{rt} = 0 \implies r^2 + \omega^2 = 0 \implies r = \pm i\omega$$
The general complex solution is:
$$x(t) = C_1 e^{i\omega t} + C_2 e^{-i\omega t}$$
Using Euler's formula ($e^{i\theta} = \cos\theta + i\sin\theta$) and enforcing real boundary conditions yields the real physical form:
$$x(t) = A \cos(\omega t + \phi)$$
Where: * $A$ is the amplitude of oscillation ($m$). * $\omega$ is the angular frequency ($\text{rad/s}$), related to period $T$ by $\omega = \frac{2\pi}{T} = 2\pi f$. * $\phi$ is the phase constant determined by initial boundary conditions $x(0)$ and $v(0)$.
Differentiating $x(t)$ with respect to time yields velocity and acceleration:
$$v(t) = \frac{dx}{dt} = -A\omega \sin(\omega t + \phi)$$
$$a(t) = \frac{d^2x}{dt^2} = -A\omega^2 \cos(\omega t + \phi) = -\omega^2 x(t)$$
2. The Physical Pendulum Derivation
Consider an arbitrarily shaped rigid body of mass $M$ pivoted about a frictionless horizontal axis at distance $d$ from its center of mass ($\text{CM}$), with rotational inertia $I$ about the pivot point.
Pivot (O)
*
/ \
/ \ d
/ θ \
/_______* Center of Mass (CM)
| |
| Rigid | Mg (Vector downward)
| Body |
\_______/
Applying Newton’s Second Law for Rotation:
$$\sum \tau_O = I \alpha \implies -M g d \sin\theta = I \frac{d^2\theta}{dt^2}$$
$$\frac{d^2\theta}{dt^2} + \frac{M g d}{I} \sin\theta = 0$$
For small angular displacements ($\theta \ll 1\text{ rad}$), we apply the Maclaurin series expansion $\sin\theta = \theta - \frac{\theta^3}{3!} + \dots \approx \theta$. The equation reduces to:
$$\frac{d^2\theta}{dt^2} + \left(\frac{M g d}{I}\right)\theta = 0$$
This matches the SHM canonical form $\frac{d^2\theta}{dt^2} + \omega^2 \theta = 0$, yielding:
$$\omega = \sqrt{\frac{M g d}{I}} \implies T = \frac{2\pi}{\omega} = 2\pi \sqrt{\frac{I}{M g d}}$$
B. Python Numerical Simulation: Non-Linear Pendulum Integration
When the small-angle approximation breaks down ($\theta_0 > 15^\circ$), the linear ODE becomes non-linear. Below is a Python program utilizing the 4th-Order Runge-Kutta (RK4) integration method to simulate large-amplitude oscillatory motion and contrast it with ideal SHM.
import numpy as np
import matplotlib.pyplot as plt
# Physical Parameters
g = 9.81 # Acceleration due to gravity (m/s^2)
L = 1.0 # Length of simple pendulum (m)
omega_0 = np.sqrt(g / L) # Small-angle angular frequency
# Time discretization
dt = 0.001 # Time step (s)
t_max = 10.0 # Total duration (s)
t = np.arange(0, t_max, dt)
# System of First-Order ODEs:
# d(theta)/dt = omega
# d(omega)/dt = -(g/L)*sin(theta)
def derivatives(state):
theta, omega = state
dtheta_dt = omega
domega_dt = -(g / L) * np.sin(theta)
return np.array([dtheta_dt, domega_dt])
def solve_rk4(theta_init, omega_init):
num_steps = len(t)
state = np.zeros((num_steps, 2))
state[0] = [theta_init, omega_init]
for i in range(num_steps - 1):
y_n = state[i]
k1 = derivatives(y_n)
k2 = derivatives(y_n + 0.5 * dt * k1)
k3 = derivatives(y_n + 0.5 * dt * k2)
k4 = derivatives(y_n + dt * k3)
state[i+1] = y_n + (dt / 6.0) * (k1 + 2*k2 + 2*k3 + k4)
return state[:, 0], state[:, 1]
# Simulate small angle (10 deg) vs large angle (80 deg)
theta_small_rad = np.radians(10.0)
theta_large_rad = np.radians(80.0)
theta_small, _ = solve_rk4(theta_small_rad, 0.0)
theta_large, _ = solve_rk4(theta_large_rad, 0.0)
# Ideal Harmonic Model for Large Angle Comparison
shm_large = theta_large_rad * np.cos(omega_0 * t)
# Plotting Results
plt.figure(figsize=(10, 5))
plt.plot(t, np.degrees(theta_large), 'r-', label=r'Exact RK4 Solution ($\theta_0 = 80^\circ$)')
plt.plot(t, np.degrees(shm_large), 'b--', label=r'Linearized SHM Approximation ($\theta_0 = 80^\circ$)')
plt.title("Non-Linear Pendulum Dynamics vs. Linearized SHM", fontsize=12)
plt.xlabel("Time (s)", fontsize=10)
plt.ylabel("Angular Position (Degrees)", fontsize=10)
plt.grid(True)
plt.legend(loc='upper right')
plt.savefig("pendulum_comparison.png", dpi=300)
plt.show()
3. Common AP Exam Pitfalls & Score 5 Scoring Rubric Nuances
To secure a Score 5, your mathematical derivations must be rigorous and free of unstated assumptions. Below are common student errors contrasted with high-scoring solutions.
+----------------------------------------------------------------------------------------------------+
| SCORE 4 vs. SCORE 5 COMPARISON |
+--------------------------------------------------+-------------------------------------------------+
| SCORE 4 RESPONSE (Common Pitfalls) | SCORE 5 RESPONSE (Exemplary Physics C Solution) |
+--------------------------------------------------+-------------------------------------------------+
| States T = 2π√(I / mgd) directly for a physical | Proves SHM first by calculating net torque |
| pendulum without deriving it from τ_net = I*α. | τ_net = -mgd sinθ, stating sinθ ≈ θ for θ << 1, |
| | then showing d²θ/dt² = -(mgd/I)θ = -ω²θ. |
+--------------------------------------------------+-------------------------------------------------+
| Evaluates maximum speed in non-quadratic potential| Uses U(x) = -∫ F dx to find exact potential, |
| energy functions using v_max = ωA (only valid | then sets E_total = K(x) + U(x) = U(A) and |
| for quadratic potential U = 1/2 k x²). | solves v_max = √((2/m)(U(A) - U_min)). |
+--------------------------------------------------+-------------------------------------------------+
| Forgets parallel-axis theorem when calculating | Explicitly calculates I_pivot = I_cm + md² |
| rotational inertia I for pivoted extended bodies.| before substituting I into the period equation. |
+--------------------------------------------------+-------------------------------------------------+
| Treats equivalent spring constants incorrectly | Derives k_eff for series (1/k_eq = 1/k1 + 1/k2) |
| without showing force-balance analysis. | or parallel (k_eq = k1 + k2) using displacement |
| | and force consistency equations. |
+--------------------------------------------------+-------------------------------------------------+
Critical Scoring Rubric Nuances
- The "Prove SHM" Mandate: If an FRQ asks you to "Show that the system undergoes simple harmonic motion," you will receive 0 points for simply quoting an equation for period $T$. You must:
- Set up Newton's Second Law ($\sum F = ma$ or $\sum \tau = I\alpha$).
- Perform algebraic manipulations to reach the form $\frac{d^2x}{dt^2} = -C x$ (where $C$ is a constant composed of system properties).
- Conclude explicitly: "Because acceleration is proportional to negative displacement, the motion is simple harmonic with $\omega = \sqrt{C}$."
- Phase Angle Determination: Always check $x(0)$ and $v(0)$. If $x(0) = 0$ and $v(0) = +v_{\max}$, the solution is best expressed as $x(t) = A\sin(\omega t)$, or $x(t) = A\cos(\omega t - \pi/2)$. Writing $x(t) = A\cos(\omega t)$ in this scenario results in loss of phase points.
4. Stanford University Placement Pathway
Institutional Exemption Benchmark
- Exempted Course: PHYSICS 41 (Mechanics, 4 Units).
- Target Score: 5 on AP Physics C: Mechanics.
- Placement Trajectory: Direct eligibility to enroll in PHYSICS 43 (Electricity & Magnetism) or move directly to advanced mechanical engineering prerequisites such as ENGR 14 (Applied Mechanics) and PHYSICS 110 (Advanced Mechanics).
Strategic Advantage for Engineering Science Track
Mastering calculus-driven oscillations provides a direct structural advantage in higher-level coursework at Stanford:
[AP Physics C: Mechanics (Score 5)]
│
▼
[PHYSICS 41 Exemption (4 Units)]
│
┌──────────┴────────────────────────┐
▼ ▼
[PHYSICS 43: E&M] [ENGR 14: Applied Mechanics]
(Oscillatory AC circuits, (Vibrational analysis, modal
LC/RLC differential ODEs) frequencies, damping matrices)
│ │
▼ ▼
[CME 102 / ENGR 155] [PHYSICS 110 / ME 161]
(Differential Equations for (Lagrangian Dynamics, Advanced
Engineers: Forced Oscillations) Vibrations & Signal Processing)
- Signal Processing & Fourier Analysis: The differential operator mechanics mastered in SHM form the backbone of spectral decomposition, Fourier series, and frequency-domain dynamics used in Stanford’s Electrical Engineering (
EE 102A) and Mechanical Engineering (ME 161) curricula. - Complex ODE Foundations: Linear constant-coefficient differential equations learned in calculus-driven SHM are directly applied in
CME 102(Ordinary Differential Equations for Engineers).
5. High-Yield Practice Problem & Step-by-Step Solution Checklist
Problem
A uniform thin rod of length $L$ and total mass $M$ is pivoted smoothly about a frictionless horizontal hinge located at a distance $x$ above its center of mass ($\text{CM}$). A small block of mass $m = \frac{M}{2}$ is rigidly attached to the absolute bottom tip of the rod.
Pivot Point O
o------------| (Distance x above CM)
| |
| |
| Rod (M) | Distance (L/2 - x)
| |
|----------- CM
| |
| | Distance (L/2)
| |
[m = M/2]----+ Bottom Tip
- Part A: Calculate the total rotational inertia $I_O$ of the combined system about the pivot point $O$ as a function of $M$, $L$, and $x$.
- Part B: Derive the differential equation of motion governing small-angle oscillations $\theta(t)$ of the system about its equilibrium position.
- Part C: Express the angular frequency $\omega$ of small oscillations in terms of $x$, $L$, and $g$.
- Part D: Determine the optimal pivot distance $x$ (expressed as a fraction of $L$) that minimizes the period of oscillation $T$.
Step-by-Step Solution & Scoring Checklist
Part A Solution
- Rotational inertia of a uniform rod about its $\text{CM}$ is $I_{\text{rod, CM}} = \frac{1}{12} M L^2$.
- Applying the Parallel Axis Theorem to find $I_{\text{rod}, O}$: $$I_{\text{rod}, O} = I_{\text{rod, CM}} + M x^2 = \frac{1}{12} M L^2 + M x^2$$
- Rotational inertia of the point-mass $m = \frac{M}{2}$ located at a distance $d_m = \left(\frac{L}{2} + x\right)$ from pivot $O$: $$I_{\text{mass}, O} = m (d_m)^2 = \left(\frac{M}{2}\right) \left(\frac{L}{2} + x\right)^2$$
- Total rotational inertia $I_O$: $$I_O = \frac{1}{12} M L^2 + M x^2 + \frac{M}{2} \left(\frac{L^2}{4} + L x + x^2\right)$$ $$I_O = M \left( \frac{1}{12} L^2 + x^2 + \frac{1}{8} L^2 + \frac{1}{2} L x + \frac{1}{2} x^2 \right) = M \left( \frac{5}{24} L^2 + \frac{1}{2} L x + \frac{3}{2} x^2 \right)$$
Part B Solution
- Find the distance $d_{\text{CM, sys}}$ from pivot $O$ to the system's combined center of mass: $$d_{\text{CM, sys}} = \frac{M(x) + m\left(\frac{L}{2} + x\right)}{M + m} = \frac{M x + \frac{M}{2}\left(\frac{L}{2} + x\right)}{\frac{3}{2} M} = \frac{\frac{3}{2} M x + \frac{1}{4} M L}{\frac{3}{2} M} = x + \frac{1}{6} L$$
- Total system mass $M_{\text{sys}} = M + \frac{M}{2} = \frac{3}{2} M$.
- Set up rotational Newton's Second Law ($\sum \tau_O = I_O \alpha$): $$\tau_{\text{net}} = - M_{\text{sys}} g d_{\text{CM, sys}} \sin\theta = I_O \frac{d^2\theta}{dt^2}$$ $$-\left(\frac{3}{2} M\right) g \left(x + \frac{1}{6} L\right) \sin\theta = I_O \frac{d^2\theta}{dt^2}$$
- Apply small-angle approximation $\sin\theta \approx \theta$: $$\frac{d^2\theta}{dt^2} + \left[ \frac{\frac{3}{2} M g \left(x + \frac{1}{6} L\right)}{I_O} \right] \theta = 0$$
Part C Solution
- From the standardized differential equation $\frac{d^2\theta}{dt^2} + \omega^2 \theta = 0$: $$\omega = \sqrt{\frac{\frac{3}{2} M g \left(x + \frac{1}{6} L\right)}{I_O}}$$
- Substitute $I_O = M \left( \frac{5}{24} L^2 + \frac{1}{2} L x + \frac{3}{2} x^2 \right)$: $$\omega = \sqrt{\frac{\frac{3}{2} g \left(x + \frac{1}{6} L\right)}{\frac{5}{24} L^2 + \frac{1}{2} L x + \frac{3}{2} x^2}} = \sqrt{\frac{g \left(\frac{3}{2} x + \frac{1}{4} L\right)}{\frac{5}{24} L^2 + \frac{1}{2} L x + \frac{3}{2} x^2}}$$
Part D Solution
- Minimizing period $T = \frac{2\pi}{\omega}$ is mathematically equivalent to maximizing $\omega^2(x)$: $$f(x) \equiv \omega^2 = g \cdot \frac{\frac{3}{2} x + \frac{1}{4} L}{\frac{3}{2} x^2 + \frac{1}{2} L x + \frac{5}{24} L^2}$$ Multiply numerator and denominator by $24$ to clean up fractions: $$\omega^2 = g \cdot \frac{36 x + 6 L}{36 x^2 + 12 L x + 5 L^2} = 6 g \cdot \frac{6x + L}{36x^2 + 12Lx + 5L^2}$$
- Take derivative $\frac{d(\omega^2)}{dx}$ and set to 0 using quotient rule: Let $u = 6x + L \implies u' = 6$. Let $v = 36x^2 + 12Lx + 5L^2 \implies v' = 72x + 12L$. $$\frac{d(\omega^2)}{dx} = 0 \implies u' v - u v' = 0$$ $$6(36x^2 + 12Lx + 5L^2) - (6x + L)(72x + 12L) = 0$$ $$216x^2 + 72Lx + 30L^2 - (432x^2 + 72Lx + 72Lx + 12L^2) = 0$$ $$-216x^2 - 72Lx + 18L^2 = 0$$
- Divide by $-18$: $$12x^2 + 4Lx - L^2 = 0$$
- Factor the quadratic equation: $$(6x - L)(2x + L) = 0$$
- Since $x > 0$ (pivot location above $\text{CM}$): $$6x = L \implies x = \frac{1}{6} L$$
AP Scoring Rubric Checklist (15-Point FRQ Style)
| Criteria | Allocated Points |
|---|---|
| Part A: Correct application of parallel axis theorem for rod | 1 Pt |
| Part A: Correct calculation of attached point mass inertia | 1 Pt |
| Part A: Correct total $I_O$ algebraic expression | 1 Pt |
| Part B: Correct expression for system center of mass location | 1 Pt |
| Part B: Setting up torque differential equation $\tau = I \alpha$ | 1 Pt |
| Part B: Correct application of small-angle approximation ($\sin\theta \approx \theta$) | 1 Pt |
| Part B: Final differential equation formatted as $\frac{d^2\theta}{dt^2} + \omega^2\theta = 0$ | 1 Pt |
| Part C: Substitution of $I_O$ into angular frequency expression | 2 Pts |
| Part C: Correct final simplified expression for $\omega$ | 1 Pt |
| Part D: Setting up condition for minimum $T$ / maximum $\omega^2$ via calculus | 2 Pts |
| Part D: Correct derivative computation using Quotient/Chain Rule | 1 Pt |
| Part D: Correct physical root extraction yielding $x = \frac{1}{6} L$ | 2 Pts |
| TOTAL | 15 Pts |