AP STEM PATH
Physics C: Mechanics • Score 5 Strategy

Simple Harmonic Motion & Calculus-Driven Oscillations Guide: AP Physics C: Mechanics Score 5 for UC Berkeley

AP Physics C: Mechanics Mastery Guide

Unit: Simple Harmonic Motion & Calculus-Driven Oscillations

Target Institution: UC Berkeley | Goal: Score 5 (Physics 7A Course Exemption)

1. Introduction & AP Exam Weight

Simple Harmonic Motion (SHM) and calculus-driven oscillations represent 10–14% of the AP Physics C: Mechanics exam content. However, their strategic importance is far greater: differential equations derived from oscillatory systems frequently serve as the foundational backbone of multi-topic Free Response Questions (FRQs) combining rotational dynamics, work-energy relationships, and variable-force kinematics.

To earn a Score 5 and satisfy the prerequisite standard for UC Berkeley’s Physics 7A, you must move beyond algebraic formula substitution ($T = 2\pi\sqrt{m/k}$). You are expected to demonstrate complete fluency in setting up, solving, and interpreting linear second-order homogeneous differential equations, performing Small-Angle Approximations via Taylor Series expansions, and optimizing system parameters using multi-variable calculus.

2. Deep Concept Breakdown

2.1 The Differential Equation of SHM

Simple Harmonic Motion is fundamentally defined by a linear restoring force proportional to displacement:

$$\vec{F}_{\text{net}} = -k \vec{x}$$

Applying Newton’s Second Law in one dimension yields the governing second-order differential equation:

$$m \frac{d^2x}{dt^2} = -kx \implies \frac{d^2x}{dt^2} + \frac{k}{m}x = 0$$

Defining the angular frequency $\omega \equiv \sqrt{\frac{k}{m}}$, the equation takes its canonical form:

$$\frac{d^2x}{dt^2} + \omega^2 x = 0$$

Analytical Solution Derivation

The general solution to this linear homogeneous differential equation with constant coefficients is:

$$x(t) = A \cos(\omega t + \phi)$$

Where: * $A$ is the amplitude of oscillation. * $\omega$ is the angular frequency ($\text{rad/s}$). * $\phi$ is the phase constant determined by initial boundary conditions $x(0) = x_0$ and $v(0) = v_0$.

Differentiating with respect to time $t$:

$$v(t) = \frac{dx}{dt} = -\omega A \sin(\omega t + \phi)$$

$$a(t) = \frac{d^2x}{dt^2} = -\omega^2 A \cos(\omega t + \phi) = -\omega^2 x(t)$$

Boundary Value Optimization

Given initial position $x(0) = x_0$ and initial velocity $v(0) = v_0$:

$$x(0) = A \cos(\phi) = x_0$$

$$v(0) = -\omega A \sin(\phi) = v_0 \implies A \sin(\phi) = -\frac{v_0}{\omega}$$

Squaring and summing both boundary equations yields the exact amplitude:

$$A^2 \left(\cos^2\phi + \sin^2\phi\right) = x_0^2 + \left(\frac{v_0}{\omega}\right)^2 \implies A = \sqrt{x_0^2 + \frac{v_0^2}{\omega^2}}$$

Dividing the two yields the phase angle:

$$\tan(\phi) = -\frac{v_0}{\omega x_0} \implies \phi = \arctan\left(-\frac{v_0}{\omega x_0}\right)$$

2.2 Rotational SHM & Physical Pendulums

For a rigid body of total mass $M$ and rotational inertia $I$ pivoted about a frictionless fixed axis at distance $d$ from its center of mass ($CM$):

$$\sum \tau = I \alpha \implies -M g d \sin\theta = I \frac{d^2\theta}{dt^2}$$

$$\frac{d^2\theta}{dt^2} + \frac{Mgd}{I} \sin\theta = 0$$

Small-Angle Approximation

Using the Maclaurin Series expansion for $\sin\theta$:

$$\sin\theta = \sum_{n=0}^{\infty} \frac{(-1)^n}{(2n+1)!}\theta^{2n+1} = \theta - \frac{\theta^3}{3!} + \frac{\theta^5}{5!} - \mathcal{O}(\theta^7)$$

For small angular displacements ($\theta \ll 1 \text{ rad}$, typically $\theta \le 10^\circ$), $\sin\theta \approx \theta$. The equation simplifies to SHM form:

$$\frac{d^2\theta}{dt^2} + \left(\frac{Mgd}{I}\right)\theta = 0$$

From this, the natural angular frequency $\omega$ and period $T$ are:

$$\omega = \sqrt{\frac{Mgd}{I}} \implies T = \frac{2\pi}{\omega} = 2\pi \sqrt{\frac{I}{Mgd}}$$

2.3 Exact vs. Small-Angle Dynamics: Numerical Verification

To observe where the Small-Angle Approximation fails on non-linear systems, we implement a numerical integration of the non-linear differential equation $\frac{d^2\theta}{dt^2} + \frac{g}{L}\sin\theta = 0$ using Python (scipy.integrate.solve_ivp) compared against the linearized analytical SHM solution.

import numpy as np
from scipy.integrate import solve_ivp
import matplotlib.pyplot as plt

def pendulum_dynamics(t, y, g, L):
    """
    State vector y = [theta, omega]
    d/dt [theta, omega] = [omega, -(g/L)*sin(theta)]
    """
    theta, omega = y
    return [omega, -(g / L) * np.sin(theta)]

# Physical parameters
g = 9.81  # m/s^2
L = 1.0   # meters
omega_0 = np.sqrt(g / L)

# Time span
t_span = (0, 10)
t_eval = np.linspace(0, 10, 1000)

# Initial conditions: High amplitude (60 degrees) to demonstrate non-linearity
theta_0 = np.radians(60.0)
v_0 = 0.0
y0 = [theta_0, v_0]

# Solve non-linear ODE
sol = solve_ivp(pendulum_dynamics, t_span, y0, args=(g, L), t_eval=t_eval, rtol=1e-9, atol=1e-9)

# Linearized analytical solution (SHM)
theta_shm = theta_0 * np.cos(omega_0 * t_eval)

# Plot comparison
plt.figure(figsize=(10, 5))
plt.plot(sol.t, np.degrees(sol.y[0]), label="Exact Non-Linear Differential Eq.", color="crimson", linewidth=2)
plt.plot(t_eval, np.degrees(theta_shm), label="SHM Linearized Approximation", color="navy", linestyle="--", linewidth=1.5)
plt.title(f"Physical Pendulum Motion: Small-Angle Breakdown at $\\theta_0 = 60^\\circ$")
plt.xlabel("Time (s)")
plt.ylabel("Angle $\\theta$ (degrees)")
plt.grid(True, alpha=0.3)
plt.legend(loc="upper right")
plt.savefig("shm_comparison.png", dpi=300)
plt.show()

3. Common AP Exam Pitfalls & Score 5 Scoring Rubric Nuances

Scoring Nuance: Score 4 vs. Score 5 Performance

Problem Feature Score 4 Response Path Score 5 Response Path
Differential Equation Setup Writes $\frac{d^2x}{dt^2} = -\frac{k}{m}x$ directly from memory without free-body diagram or Newton's 2nd law explicit derivation. Derives system dynamics explicitly: $\sum F_x = m a_x \implies -k(x - x_0) = m\frac{d^2x}{dt^2}$, setting up explicit boundary variables.
Physical Pendulums Mistakes moment of inertia as $I = \frac{1}{3}ML^2$ automatically without checking pivot offset or Parallel Axis Theorem. Evaluates $I = I_{cm} + M d^2$, verifies distance $d$ relative to pivot point, and uses small-angle Taylor expansion explicitly.
Phase Angle ($\phi$) Selection Assumes $\phi = 0$ blindly, ignoring initial condition statements like $x(0) = 0$ or $v(0) = -v_{\text{max}}$. Calculates $\phi = \arctan\left(-\frac{v_0}{\omega x_0}\right)$ explicitly and determines quadrant alignment correctly.
Calculus Optimizations Treats $T$ as a fixed constant in complex geometry systems instead of differentiating $T(h)$ to find period minima. Sets $\frac{dT}{dh} = 0$, confirms local minimum using $\frac{d^2T}{dh^2} > 0$, and evaluates optimized geometric limit.

4. UC Berkeley Placement Pathway

      AP Physics C: Mechanics Exam (Score 5)
                        │
                        ▼
       Exempts: Physics 7A (4 Semester Units)
   (Mechanics, Special Relativity, Wave Motion)
                        │
                        ▼
      Enrolls Directly in: Physics 7B
   (Heat, Electricity, Magnetism & Vector Fields)
                        │
                        ▼
 Acceleration Advantage for Berkeley Engineering (EECS, ME, BioE)

Academic & Prerequisite Nuance

At UC Berkeley, Physics 7A is the initial gateway physics course for all majors within the College of Engineering (COE) and the College of Chemistry, as well as Physics/Astrophysics majors in L&S.

5. High-Yield Practice Problem

Problem Statement

A uniform thin rod of mass $M$ and total length $L$ is pivoted frictionless about a horizontal axis located at a variable distance $h$ from its center of mass ($0 < h \le \frac{L}{2}$). A point-mass $m_0 = \frac{M}{2}$ is securely attached to the absolute bottom tip of the rod. A linear spring with spring constant $k$ is horizontally connected at the top end of the rod, anchored to a vertical wall.

       Wall
       |     k
       |───/\/\/\/\───┐  (Top End)
       |              │
       |              │
       |              ▲  Pivot (Distance h above CM)
       |              │
       |              │
       |            ( CM )
       |              │
       |              │
       |              ●  Point Mass m_0 = M/2 (Bottom End)
  1. [4 Points] Derivation of Differential Equation: Derivation of the second-order differential equation of motion governing small angular displacements $\theta(t)$ from equilibrium.
  2. [4 Points] System Angular Frequency: Expression for the natural angular frequency $\omega_0$ of the combined system in terms of $M, L, h, k,$ and $g$.
  3. [4 Points] Calculus Optimization: Assume spring force is negligible ($k \to 0$). Derive the distance $h_{\text{opt}}$ (in terms of $L$) that minimizes the period of oscillation $T$.
  4. [3 Points] Energy Kinetics: If the rod is released from rest at initial displacement $\theta_0 \ll 1 \text{ rad}$, derive the maximum angular velocity $\Omega_{\text{max}} = \left.\frac{d\theta}{dt}\right|_{\text{max}}$.

Step-by-Step Solution Checklist

Part (a): Derive the Differential Equation of Motion

Step 1: Calculate Total Rotational Inertia $I_{\text{pivot}}$ about the axis. Using the Parallel Axis Theorem for the rod ($I_{\text{cm, rod}} = \frac{1}{12}M L^2$) displaced by distance $h$:

$$I_{\text{rod}} = \frac{1}{12}ML^2 + M h^2$$

The point mass $m_0 = \frac{M}{2}$ is located at distance $\left(\frac{L}{2} + h\right)$ from the pivot point:

$$I_{\text{mass}} = m_0 r^2 = \left(\frac{M}{2}\right)\left(\frac{L}{2} + h\right)^2$$

Summing both terms:

$$I_{\text{pivot}} = \left(\frac{1}{12}ML^2 + M h^2\right) + \frac{M}{2}\left(\frac{L^2}{4} + L h + h^2\right)$$

$$I_{\text{pivot}} = M \left[ \frac{1}{12}L^2 + \frac{1}{8}L^2 + \frac{1}{2}L h + \frac{3}{2}h^2 \right] = M \left[ \frac{5}{24}L^2 + \frac{1}{2}L h + \frac{3}{2}h^2 \right]$$

Step 2: Apply Rotational Newton's Second Law ($\sum \tau_{\text{pivot}} = I_{\text{pivot}} \alpha$). Restoring torques acting for small angle $\theta$ counter-clockwise: 1. Gravity on Rod: $\tau_{g, \text{rod}} = -M g h \sin\theta$ 2. Gravity on Mass: $\tau_{g, \text{mass}} = -\left(\frac{M}{2}\right)g\left(\frac{L}{2} + h\right)\sin\theta$ 3. Spring Force: Moment arm from pivot to top end is $\left(\frac{L}{2} - h\right)$. Displacement of spring $x \approx \left(\frac{L}{2} - h\right)\theta$. $$\tau_{\text{spring}} = -k \left[\left(\frac{L}{2} - h\right)\theta\right] \left(\frac{L}{2} - h\right) = -k \left(\frac{L}{2} - h\right)^2 \theta$$

Summing restoring torques using small angle approximation $\sin\theta \approx \theta$:

$$\sum \tau = -\left[ M g h + \frac{M g}{2}\left(\frac{L}{2} + h\right) + k\left(\frac{L}{2} - h\right)^2 \right] \theta$$

$$\sum \tau = -\left[ \frac{M g L}{4} + \frac{3}{2}M g h + k\left(\frac{L}{2} - h\right)^2 \right] \theta$$

Equating $\sum \tau = I_{\text{pivot}} \frac{d^2\theta}{dt^2}$:

$$I_{\text{pivot}} \frac{d^2\theta}{dt^2} + \left[ \frac{MgL}{4} + \frac{3}{2}Mgh + k\left(\frac{L}{2} - h\right)^2 \right] \theta = 0$$

$$\frac{d^2\theta}{dt^2} + \left( \frac{\frac{MgL}{4} + \frac{3}{2}Mgh + k\left(\frac{L}{2} - h\right)^2}{M\left(\frac{5}{24}L^2 + \frac{1}{2}Lh + \frac{3}{2}h^2\right)} \right) \theta = 0$$

Part (b): Expression for Natural Angular Frequency $\omega_0$

From standard SHM differential form $\frac{d^2\theta}{dt^2} + \omega_0^2 \theta = 0$:

$$\omega_0 = \sqrt{\frac{\frac{MgL}{4} + \frac{3}{2}Mgh + k\left(\frac{L}{2} - h\right)^2}{M\left(\frac{5}{24}L^2 + \frac{1}{2}Lh + \frac{3}{2}h^2\right)}}$$

Part (c): Calculus Optimization for Minimum Period ($k \to 0$)

When $k = 0$:

$$\omega_0^2 = \frac{\frac{MgL}{4} + \frac{3}{2}Mgh}{M\left(\frac{5}{24}L^2 + \frac{1}{2}Lh + \frac{3}{2}h^2\right)} = \frac{g\left(\frac{L}{4} + \frac{3}{2}h\right)}{\frac{5}{24}L^2 + \frac{1}{2}Lh + \frac{3}{2}h^2} = \frac{g\left(L + 6h\right)}{\frac{5}{6}L^2 + 2Lh + 6h^2}$$

Since $T = \frac{2\pi}{\omega_0}$, minimizing $T$ is equivalent to maximizing $\omega_0^2$:

Let $f(h) = \frac{L + 6h}{\frac{5}{6}L^2 + 2Lh + 6h^2}$. Differentiating with respect to $h$ and setting $\frac{df}{dh} = 0$:

Applying the quotient rule $\frac{u'v - uv'}{v^2} = 0 \implies u'v - uv' = 0$:

$$u = L + 6h \implies u' = 6$$

$$v = \frac{5}{6}L^2 + 2Lh + 6h^2 \implies v' = 2L + 12h$$

$$6\left(\frac{5}{6}L^2 + 2Lh + 6h^2\right) - (L + 6h)(2L + 12h) = 0$$

$$5L^2 + 12Lh + 36h^2 - \left(2L^2 + 12Lh + 12Lh + 72h^2\right) = 0$$

$$5L^2 + 12Lh + 36h^2 - 2L^2 - 24Lh - 72h^2 = 0$$

$$3L^2 - 12Lh - 36h^2 = 0$$

Dividing entire polynomial by $3$:

$$-12h^2 - 4Lh + L^2 = 0 \implies 12h^2 + 4Lh - L^2 = 0$$

Factoring the quadratic in $h$:

$$(6h - L)(2h + L) = 0$$

Since $h > 0$ physically:

$$6h - L = 0 \implies h_{\text{opt}} = \frac{L}{6}$$

Part (d): Maximum Angular Velocity Calculation

Using Conservation of Mechanical Energy or amplitude relation of SHM:

$$\theta(t) = \theta_0 \cos(\omega_0 t)$$

$$\Omega(t) = \frac{d\theta}{dt} = -\omega_0 \theta_0 \sin(\omega_0 t)$$

$$\Omega_{\text{max}} = \omega_0 \theta_0$$

Substituting $h = h_{\text{opt}} = \frac{L}{6}$ back into $\omega_0$ (with $k=0$):

$$\omega_0^2 = \frac{g\left(L + 6\left(\frac{L}{6}\right)\right)}{\frac{5}{6}L^2 + 2L\left(\frac{L}{6}\right) + 6\left(\frac{L}{6}\right)^2} = \frac{g(2L)}{\frac{5}{6}L^2 + \frac{1}{3}L^2 + \frac{1}{6}L^2} = \frac{2gL}{\frac{4}{3}L^2} = \frac{3g}{2L}$$

$$\omega_0 = \sqrt{\frac{3g}{2L}}$$

Thus, the maximum angular velocity is:

$$\Omega_{\text{max}} = \theta_0 \sqrt{\frac{3g}{2L}}$$

Official AP Point Breakdown Rubric

Aiming for a Score 5 in Physics C: Mechanics?

Secure admission and advanced standing at top institutions like UC Berkeley with elite 1-on-1 AP STEM mentorship.

無料相談・学習プラン診断