AP Physics C: Mechanics Mastery Guide

Simple Harmonic Motion & Calculus-Driven Oscillations

---

1. Introduction & AP Exam Weight

Simple Harmonic Motion (SHM) and differential oscillations represent one of the most mathematically rigorous units on the AP Physics C: Mechanics exam, accounting for 10%–14% of the total score. While introductory physics treats harmonic motion through algebraic plug-and-chug formulas ($T = 2\pi \sqrt{m/k}$), the AP Physics C exam demands a calculus-first approach.

To earn a 5, you must be prepared to: 1. Derive non-standard second-order differential equations from Newton’s Second Law ($\sum F = ma$) or Rotational Dynamics ($\sum \tau = I\alpha$). 2. Apply Taylor series expansions (such as the small-angle approximation $\sin\theta \approx \theta$) to linearize non-linear restoring forces. 3. Solve oscillatory mechanics using energy conservation methods with definite integrals. 4. Calculate exact angular frequencies ($\omega$), periods ($T$), and state equations for complex systems such as physical pendulums, torsional oscillators, and variable-mass spring systems.

---

2. Deep Concept Breakdown

2.1 The Differential Equation of SHM

Simple Harmonic Motion is uniquely defined by a dynamic system where the acceleration of a body is directly proportional to its displacement from equilibrium and directed oppositely:

$$a(t) = -\omega^2 x(t)$$

Expressed as a homogeneous, second-order linear ordinary differential equation (ODE):

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

Analytical Solution Derivation

To solve this differential equation, we assume a continuous trial solution (ansatz) of the form:

$$x(t) = C e^{rt}$$

Substituting $x(t)$ 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$$

Thus, the general solution in complex exponential form is:

$$x(t) = c_1 e^{i\omega t} + c_2 e^{-i\omega t}$$

Applying Euler's formula ($e^{i\theta} = \cos\theta + i\sin\theta$) and applying boundary conditions $x(0) = A$ and $v(0) = 0$, we arrive at the standard kinematically real sinusoidal form:

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

Differentiating with respect to time $t$:

$$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)$$

Where: * $A$ = Amplitude (maximum linear or angular displacement) * $\omega$ = Angular frequency ($\text{rad/s}$), where $\omega = 2\pi f = \frac{2\pi}{T}$ * $\phi$ = Phase constant (determined by initial conditions at $t = 0$)

---

2.2 Rotational Oscillations: The Physical Pendulum

Unlike a simple pendulum with a point mass $m$ suspended from a massless string of length $L$, a physical pendulum consists of an arbitrary rigid body of mass $M$ pivoted about a frictionless axis passing through point $P$, located a distance $d$ from its Center of Mass ($\text{CM}$).

       Pivot (P)
         o
        /| 
       / | 
      /  | d
     /   v
    o--- (CM)
   /   \
  /     \  Rigid Body
 (_______)

Analytical Derivation of Angular Frequency $\omega$

Applying Newton's Second Law for Rotation:

$$\sum \tau_P = I_P \alpha$$

The restoring torque is produced exclusively by gravity acting at the Center of Mass:

$$\tau_P = -M g d \sin\theta$$

Equating torque to rotational inertia:

$$I_P \frac{d^2\theta}{dt^2} = -M g d \sin\theta \implies \frac{d^2\theta}{dt^2} + \frac{M g d}{I_P} \sin\theta = 0$$

For small angular displacements ($\theta \ll 1 \text{ rad}$), we apply the first-order Taylor polynomial expansion:

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

Substituting this approximation linearizes the differential equation into standard SHM form:

$$\frac{d^2\theta}{dt^2} + \left(\frac{M g d}{I_P}\right)\theta = 0$$

Comparing this directly to $\frac{d^2\theta}{dt^2} + \omega^2 \theta = 0$, we read off the angular frequency and period:

$$\omega = \sqrt{\frac{M g d}{I_P}}$$

$$T = \frac{2\pi}{\omega} = 2\pi \sqrt{\frac{I_P}{M g d}}$$

Note: $I_P$ must be calculated relative to the pivot point using the Parallel Axis Theorem if the rotational inertia about the center of mass $I_{\text{cm}}$ is known:

$$I_P = I_{\text{cm}} + M d^2$$

---

2.3 Energy Conservation Proof via Calculus

The total mechanical energy $E$ of a linear simple harmonic oscillator is:

$$E = U + K = \frac{1}{2} k x^2 + \frac{1}{2} m v^2$$

To prove energy conservation analytically, take the time derivative of $E$:

$$\frac{dE}{dt} = \frac{d}{dt} \left( \frac{1}{2} k x^2 + \frac{1}{2} m v^2 \right) = k x \frac{dx}{dt} + m v \frac{dv}{dt}$$

Since $v = \frac{dx}{dt}$ and $a = \frac{dv}{dt}$:

$$\frac{dE}{dt} = k x v + m v a = v (k x + m a)$$

Because $F_{\text{net}} = m a = -k x$ for an ideal oscillator, $k x + m a = 0$. Therefore:

$$\frac{dE}{dt} = v (0) = 0 \implies E = \text{constant}$$

Substituting $x(t) = A \cos(\omega t + \phi)$ and $v(t) = -A\omega \sin(\omega t + \phi)$ with $\omega^2 = \frac{k}{m}$:

$$E = \frac{1}{2} k A^2 \cos^2(\omega t + \phi) + \frac{1}{2} m \left( -A \omega \sin(\omega t + \phi) \right)^2$$

$$E = \frac{1}{2} k A^2 \cos^2(\omega t + \phi) + \frac{1}{2} m A^2 \left(\frac{k}{m}\right) \sin^2(\omega t + \phi)$$

$$E = \frac{1}{2} k A^2 \left( \cos^2(\omega t + \phi) + \sin^2(\omega t + \phi) \right) = \frac{1}{2} k A^2$$

---

2.4 Numerical Simulation: Linear vs. Non-Linear Pendulum Oscillation

When angles grow large, the small-angle approximation fails, and period lengthening occurs. The exact non-linear equation $\frac{d^2\theta}{dt^2} + \frac{g}{L}\sin\theta = 0$ requires numerical integration.

Below is a production-grade Python computational script using the Euler-Cromer integration scheme to demonstrate phase space trajectories $(x, v)$ and compare exact non-linear oscillations versus the linearized SHM approximation.

import numpy as np
import matplotlib.pyplot as plt

def simulate_pendulum(theta0_deg, L=1.0, g=9.81, t_max=10.0, dt=0.001):
    """
    Simulates linear and non-linear pendulum dynamics.

    Parameters:
        theta0_deg (float): Initial angular position in degrees.
        L (float): Length of the pendulum (m).
        g (float): Gravitational acceleration (m/s^2).
        t_max (float): Total time for simulation (s).
        dt (float): Time step resolution (s).
    """
    theta0 = np.radians(theta0_deg)
    omega0 = 0.0

    N = int(t_max / dt)
    time = np.linspace(0, t_max, N)

    # Non-linear arrays
    theta_nl = np.zeros(N)
    omega_nl = np.zeros(N)
    theta_nl[0] = theta0
    omega_nl[0] = omega0

    # Linearized SHM arrays
    theta_lin = np.zeros(N)
    omega_lin = np.zeros(N)
    theta_lin[0] = theta0
    omega_lin[0] = omega0

    omega_sq = g / L

    # Euler-Cromer Integration Scheme
    for i in range(1, N):
        # Non-Linear Execution: d^2(theta)/dt^2 = - (g/L)*sin(theta)
        alpha_nl = -omega_sq * np.sin(theta_nl[i-1])
        omega_nl[i] = omega_nl[i-1] + alpha_nl * dt
        theta_nl[i] = theta_nl[i-1] + omega_nl[i] * dt

        # Linear Execution: d^2(theta)/dt^2 = - (g/L)*theta
        alpha_lin = -omega_sq * theta_lin[i-1]
        omega_lin[i] = omega_lin[i-1] + alpha_lin * dt
        theta_lin[i] = theta_lin[i-1] + omega_lin[i] * dt

    return time, theta_nl, omega_nl, theta_lin, omega_lin

# Execute simulation at large initial angle (60 degrees) to display non-linear departure
time, theta_nl, omega_nl, theta_lin, omega_lin = simulate_pendulum(theta0_deg=60.0)

# Phase Space Plotting
plt.figure(figsize=(10, 5))
plt.plot(theta_nl, omega_nl, label="Exact Non-Linear ($\sin\\theta$)", color="red", linewidth=1.5)
plt.plot(theta_lin, omega_lin, label="Linearized SHM ($\\theta$)", color="blue", linestyle="--", linewidth=1.5)
plt.title("Phase Space Trajectory $(\\theta \\text{ vs } \\omega)$ at $\\theta_0 = 60^\\circ$")
plt.xlabel("Angular Displacement $\\theta$ [rad]")
plt.ylabel("Angular Velocity $\\omega$ [rad/s]")
plt.grid(True)
plt.legend(loc="upper right")
plt.tight_layout()
plt.show()

---

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

Critical Exam Pitfalls

1. Failure to Prove SHM via Differential Equations
2. The Mistake: Writing down $T = 2\pi\sqrt{m/k}$ or stating "the system oscillates, so $\omega = \sqrt{\dots}$" without establishing the second-order differential equation.

3. The Fix: On AP Physics C Free Response Questions (FRQs), you must set up $\sum F = m \frac{d^2x}{dt^2}$ or $\sum \tau = I \frac{d^2\theta}{dt^2}$ and rearrange it explicitly into the standard template $\frac{d^2x}{dt^2} + C x = 0$ before identifying $\omega = \sqrt{C}$.

4. Miscalculating the Pivot Inertia ($I_P$) in Physical Pendulums

5. The Mistake: Substituting $I_{\text{cm}}$ into $T = 2\pi\sqrt{I/(mgd)}$ instead of computing $I_P = I_{\text{cm}} + md^2$.

6. The Fix: Always define $d$ explicitly as the distance from the pivot point to the object's center of mass, and evaluate $I_P$ about the pivot point using the Parallel Axis Theorem.

7. Neglecting Torsional and Spring Constants Signage

8. The Mistake: Forgetting the restoring negative sign when setting up torque/force equations, leading to incorrect differential equation signs ($\frac{d^2x}{dt^2} - \omega^2 x = 0$), which describes exponential runaway growth rather than oscillation.

---

Score 4 vs. Score 5 Solution Nuance Comparison

Example Problem:

A uniform rod of mass $M$ and length $L$ is pivoted at one end. A spring of constant $k$ is attached at the opposite end, perpendicular to the rod. Derive the angular frequency $\omega$ for small oscillations.

Pivot
  o-------------------------===== Spring (k)
  |<------- Length L ------->|

================================================================================
SCORE 4 STUDENT RESPONSE (Sloppy & Inferential)
================================================================================
Torque = Force * distance
Torque of spring = -k * x * L
x = L * theta
Torque = -k * L^2 * theta
I for rod = (1/3) * M * L^2
(1/3) M L^2 * alpha = -k L^2 theta
alpha = -3 k theta / M
T = 2*pi * sqrt(M / 3k)
omega = sqrt(3k / M)

Critique: The student gets the final answer correct, but skips the explicit differential 
form statement required by standard College Board rubrics. On a tough FRQ grading 
scale, this loses the point awarded for "stating a differential equation that 
models the physical system."
================================================================================

================================================================================
SCORE 5 STUDENT RESPONSE (Rigorous, Formal Derivation)
================================================================================
Step 1: Set up the rotational second law relative to the pivot axis P:
\sum \tau_P = I_P \alpha

Step 2: Calculate rotational inertia about the end pivot:
I_P = I_{\text{cm}} + M d^2 = \frac{1}{12} M L^2 + M \left(\frac{L}{2}\right)^2 = \frac{1}{3} M L^2

Step 3: Relate linear spring stretch x to small angular displacement \theta:
For small \theta: x \approx L \theta
Restoring force of spring: F_s = -k x = -k L \theta
Restoring torque: \tau_s = -F_s \cdot L = -k L^2 \theta

Step 4: Formulate the differential equation:
I_P \frac{d^2\theta}{dt^2} = -k L^2 \theta
\left(\frac{1}{3} M L^2\right) \frac{d^2\theta}{dt^2} + k L^2 \theta = 0

Step 5: Simplify into standard homogeneous SHM ODE form:
\frac{d^2\theta}{dt^2} + \left( \frac{3k}{M} \right) \theta = 0

Step 6: Identify \omega^2 from standard form \frac{d^2\theta}{dt^2} + \omega^2 \theta = 0:
\omega^2 = \frac{3k}{M} \implies \omega = \sqrt{\frac{3k}{M}}

Critique: Explicit differential setup, rigorous geometric substitution, full point 
retention under strict AP rubric criteria.
================================================================================

---

4. Carnegie Mellon University Placement Pathway

At Carnegie Mellon University (CMU), students in the Carnegie Institute of Technology (CIT) and Mellon College of Science (MCS) who achieve a Score of 5 on the AP Physics C: Mechanics exam receive credit for:

• Exempted Course: 33-141 Physics I for Engineering (12 Units)

                              AP Physics C: Mechanics Score 5
                                            │
                                            ▼
                              Waive: 33-141 Physics I (12 Units)
                                            │
                    ┌───────────────────────┴───────────────────────┐
                    ▼                                               ▼
     Direct Acceleration Track                      Engineering Core Prerequisites
  33-142: Physics II for Engineering               18-100: Intro to ECE
  (Calculus-driven E&M, Maxwell's Eq.)              24-101: Intro to MechE
                    │                               16-299: Intro to Robotics
                    │                                               │
                    └───────────────────────┬───────────────────────┘
                                            ▼
                           Early Advancement into Advanced Topics:
                       • State-Space System Dynamics
                       • Linear Feedback Controls
                       • Structural Vibration Analysis

Strategic Academic & Engineering Advantages

1. Immediate Acceleration to 33-142 (Physics II for Engineering): Waiving 33-141 allows students to register for 33-142 in their first semester. 33-142 applies multivariate calculus to field theory and electrodynamics. Understanding the differential equation of mechanical SHM ($\ddot{x} + \omega^2 x = 0$) directly mirrors second-order electrical circuits ($LC$ and $RLC$ oscillators):

$$\frac{d^2 q}{dt^2} + \frac{R}{L}\frac{dq}{dt} + \frac{1}{LC} q = 0$$

1. Vital Foundation for Dynamic Systems and Hardware Engineering:
2. Electrical & Computer Engineering (ECE): Second-order differential equations introduced in AP Physics C Oscillations serve as the direct prerequisite foundation for 18-202 (Mathematical Foundations of Computer Engineering) and 18-220 (Electronic Devices & Circuits).
3. Robotics Institute (16-xxx) & Mechanical Engineering (24-xxx): Modern robotics relies on modeling dynamic systems via modern feedback controls. Mass-spring-damper equations dictate robotic joint impedance, actuator compliance, and structural resonance prevention in high-frequency manipulators.

---

5. High-Yield Practice Problem & Step-by-Step Solution Checklist

The Problem

A physical pendulum consists of a non-uniform rigid rod of total mass $M$ and length $L$, pivoted at a frictionless hinge at top end $O$. The rod's mass per unit length $\lambda(x)$ varies linearly from the pivot point ($x=0$) to the bottom end ($x=L$) according to:

$$\lambda(x) = b x$$

where $b$ is a known positive constant. A light spring with force constant $k$ is attached horizontally to the rod at a distance $h = \frac{3}{4} L$ below the pivot $O$, anchored to a wall.

 Wall
  |  Pivot O
  |== o 
  |   | \
  |   |  \
  |---|---vvvvv--- Spring (k)  [at x = 3/4 L]
  |   |    \
  |   |     \  Non-uniform rod: lambda(x) = b x
  |   v      \

(a) Determine the constant $b$ in terms of total mass $M$ and length $L$.
(b) Derive the position of the rod's Center of Mass ($x_{\text{cm}}$) relative to pivot $O$.
(c) Calculate the rotational inertia $I_O$ of the rod about the pivot $O$ using integration.
(d) Set up the second-order differential equation governing the angular displacement $\theta(t)$ for small oscillations ($\theta \ll 1$).
(e) Express the period of small oscillations $T$ in terms of $M, L, k, g,$ and fundamental constants.

---

Step-by-Step Solution & Scoring Checklist

Part (a): Determine Constant $b$

To find total mass $M$, integrate the linear density $\lambda(x) = bx$ over the length $L$:

$$M = \int_0^L \lambda(x) \, dx = \int_0^L b x \, dx = \left[ \frac{1}{2} b x^2 \right]_0^L = \frac{1}{2} b L^2$$

Solve for $b$:

$$b = \frac{2M}{L^2}$$

• Scoring Rubric - Part (a) [2 Points]:
• +1 Point: Correct setup of integral $M = \int \lambda(x) dx$ with proper limits ($0$ to $L$).
• +1 Point: Correct integration yielding $b = \frac{2M}{L^2}$.

---

Part (b): Calculate Center of Mass $x_{\text{cm}}$

Apply the continuous center of mass definition:

$$x_{\text{cm}} = \frac{1}{M} \int_0^L x \lambda(x) \, dx = \frac{1}{M} \int_0^L x (b x) \, dx = \frac{b}{M} \int_0^L x^2 \, dx$$

Substitute $b = \frac{2M}{L^2}$:

$$x_{\text{cm}} = \left(\frac{2M/L^2}{M}\right) \left[ \frac{x^3}{3} \right]_0^L = \frac{2}{L^2} \left( \frac{L^3}{3} \right) = \frac{2}{3} L$$

• Scoring Rubric - Part (b) [2 Points]:
• +1 Point: Setting up the integral $x_{\text{cm}} = \frac{1}{M}\int x \lambda(x) dx$.
• +1 Point: Correct integration yielding $x_{\text{cm}} = \frac{2}{3} L$.

---

Part (c): Calculate Rotational Inertia $I_O$

By definition of moment of inertia for a continuous rigid body:

$$I_O = \int_0^L x^2 dm = \int_0^L x^2 \lambda(x) \, dx = \int_0^L x^2 (b x) \, dx = b \int_0^L x^3 \, dx$$

$$I_O = \left(\frac{2M}{L^2}\right) \left[ \frac{x^4}{4} \right]_0^L = \left(\frac{2M}{L^2}\right) \left( \frac{L^4}{4} \right) = \frac{1}{2} M L^2$$

• Scoring Rubric - Part (c) [2 Points]:
• +1 Point: Correct calculus setup $I_O = \int x^2 \lambda(x) dx$.
• +1 Point: Correct analytical result $I_O = \frac{1}{2} M L^2$.

---

Part (d): Derive Differential Equation

Identify restoring torques acting about pivot point $O$: 1. Gravitational Torque (acting at $x_{\text{cm}}$):

$$\tau_g = -M g x_{\text{cm}} \sin\theta = -M g \left(\frac{2}{3} L\right) \sin\theta$$

1. Spring Restoring Torque (acting at $h = \frac{3}{4} L$): For small angular displacement $\theta$, linear stretch is $y = h \theta = \frac{3}{4} L \theta$. Spring force $F_s = -k y = -k \left(\frac{3}{4} L \theta\right)$.

$$\tau_s = -F_s \cdot h = -k \left(\frac{3}{4} L \theta\right) \left(\frac{3}{4} L\right) = -\frac{9}{16} k L^2 \theta$$

Apply Newton's Second Law for Rotation $\sum \tau_O = I_O \alpha$:

$$I_O \frac{d^2\theta}{dt^2} = \tau_g + \tau_s$$

$$\left(\frac{1}{2} M L^2\right) \frac{d^2\theta}{dt^2} = -M g \left(\frac{2}{3} L\right) \sin\theta - \frac{9}{16} k L^2 \theta$$

Apply small angle approximation $\sin\theta \approx \theta$:

$$\left(\frac{1}{2} M L^2\right) \frac{d^2\theta}{dt^2} + \left( \frac{2}{3} M g L + \frac{9}{16} k L^2 \right) \theta = 0$$

Divide the entire equation by $I_O = \frac{1}{2} M L^2$:

$$\frac{d^2\theta}{dt^2} + \left( \frac{\frac{2}{3} M g L + \frac{9}{16} k L^2}{\frac{1}{2} M L^2} \right) \theta = 0$$

$$\frac{d^2\theta}{dt^2} + \left( \frac{4 g}{3 L} + \frac{9 k}{8 M} \right) \theta = 0$$

• Scoring Rubric - Part (d) [4 Points]:
• +1 Point: Statement of $\sum \tau = I \alpha$ with two torque terms having consistent signs.
• +1 Point: Correct application of distance to center of mass for gravity torque and correct displacement application for spring torque.
• +1 Point: Explicit application of small-angle approximation $\sin\theta \approx \theta$.
• +1 Point: Algebraic reduction into standard SHM differential form $\frac{d^2\theta}{dt^2} + \omega^2 \theta = 0$.

---

Part (e): Compute Oscillatory Period $T$

From the standard form differential equation $\frac{d^2\theta}{dt^2} + \omega^2 \theta = 0$, we extract:

$$\omega^2 = \frac{4 g}{3 L} + \frac{9 k}{8 M} = \frac{32 M g + 27 k L}{24 M L}$$

$$\omega = \sqrt{\frac{32 M g + 27 k L}{24 M L}}$$

Since period $T = \frac{2\pi}{\omega}$:

$$T = 2\pi \sqrt{\frac{24 M L}{32 M g + 27 k L}}$$

• Scoring Rubric - Part (e) [2 Points]:
• +1 Point: Equating coefficient of $\theta$ to $\omega^2$.
• +1 Point: Correct substitution into $T = \frac{2\pi}{\omega}$ yielding exact simplified expression.

---

Final Master Checklist for Exam Day

• [ ] Always set up Newton's Second Law ($\sum F = m a$ or $\sum \tau = I \alpha$) using explicit derivative notation ($\frac{d^2x}{dt^2}$ or $\frac{d^2\theta}{dt^2}$).
• [ ] Verify that pivot torque uses moment of inertia evaluated at that specific pivot point ($I_P = I_{\text{cm}} + M d^2$).
• [ ] Apply linear/small-angle approximations explicitly ($\sin\theta \approx \theta$, $\tan\theta \approx \theta$, $(1+x)^n \approx 1+nx$).
• [ ] Put second-order ODEs into canonical form $\frac{d^2x}{dt^2} + \omega^2 x = 0$ to identify angular frequency without hand-waving.