Hockey puck with mathematical formulas and trajectory curves on ice

The Mathematics of Hockey

Discover how algebra, geometry, calculus, and physics shape every shot, pass, and save on the ice. Interactive visualizations bring the math to life.

100+ MPH Shots
200' Rink Length
Math Concepts

Algebra: The Language of Hockey

From calculating shooting percentages to predicting player statistics, algebra is fundamental to understanding hockey performance.

Shot Velocity and Distance
Middle School
High School
Advanced

Understanding Speed = Distance ÷ Time

When a hockey player shoots the puck, we can calculate how fast it's moving using a simple formula.

\(v = \frac{d}{t}\)

What each letter means:

  • v = velocity (speed) of the puck in feet per second
  • d = distance the puck travels in feet
  • t = time it takes in seconds

Example 1: Fast Shot

If a puck travels 60 feet in 0.5 seconds, its speed is 60 ÷ 0.5 = 120 feet per second!

Example 2: Slower Pass

A gentle pass covering 30 feet in 1.5 seconds has a speed of 30 ÷ 1.5 = 20 feet per second.

Example 3: Pro Shot

Professional players can shoot the puck at over 100 MPH, which is about 147 feet per second.

Example 4: Rink Length

A puck traveling the full rink length (200 feet) in 2 seconds moves at 200 ÷ 2 = 100 ft/s (68 MPH).

Linear Equations and Rate Problems

We can rearrange the velocity equation to solve for different variables and analyze more complex scenarios.

\(d = v \times t \quad\quad t = \frac{d}{v} \quad\quad v = \frac{d}{t}\)

Extended Analysis:

  • Converting units: 1 MPH = 1.467 ft/s
  • Average velocity: \(v_{avg} = \frac{v_{initial} + v_{final}}{2}\)
  • Acceleration: \(a = \frac{v_{final} - v_{initial}}{t}\)

Example 1: Race to the Goal

Player A shoots from 60 ft at 120 ft/s (t = 0.5s). Player B from 45 ft at 90 ft/s (t = 0.5s). They arrive simultaneously!

Example 2: Timing the Shot

A puck shot at 90 MPH (132 ft/s) must travel 89 feet (goal line to goal line). Time = 89/132 ≈ 0.67 seconds.

Example 3: Unit Conversion

Convert 100 MPH to ft/s: 100 × 1.467 = 146.7 ft/s. This is crucial for comparing speeds in different units.

Example 4: Acceleration Problem

Puck goes from 0 to 88 ft/s (60 MPH) in 0.1s during a slap shot. Acceleration = 88/0.1 = 880 ft/s²!

Vector Algebra and Momentum Conservation

Advanced algebraic concepts model puck collisions, deflections, and multi-directional motion.

\(m_1v_1 + m_2v_2 = m_1v_1' + m_2v_2' \quad \text{(Conservation of Momentum)}\)

Advanced Concepts:

  • Elastic collisions: Both momentum and kinetic energy conserved
  • Coefficient of restitution: \(e = \frac{v_2' - v_1'}{v_1 - v_2}\)
  • Vector decomposition: \(\vec{v} = v_x \hat{i} + v_y \hat{j}\)
  • Impulse-momentum theorem: \(\vec{F}\Delta t = m\Delta\vec{v}\)

Deflection Analysis

When a puck (6 oz) deflects off a stick blade, momentum conservation predicts the new trajectory angle and velocity.

Real Application

NHL teams use vector analysis to optimize tip-in locations where deflected shots are most likely to score.

Interactive: Calculate Puck Speed

Calculated Velocity
120.0
feet per second (81.8 MPH)
Shooting Percentage
Middle School
High School
Advanced

Understanding Percentages

Shooting percentage tells us how often a player scores when they shoot the puck.

\(\text{Shooting \%} = \frac{\text{Goals}}{\text{Shots}} \times 100\)

Example 1: Basic Calculation

If a player scores 15 goals on 100 shots, their shooting percentage is (15/100) × 100 = 15%.

Example 2: Elite Sniper

20 goals on 90 shots = (20/90) × 100 = 22.2%. This is an elite level shooting percentage!

Example 3: Average Player

An average NHL player with 8 goals on 80 shots = (8/80) × 100 = 10%.

Example 4: Career Stats

Over 5 seasons: 150 goals on 1,200 shots = (150/1200) × 100 = 12.5% career shooting percentage.

Ratios, Proportions, and Statistical Analysis

We can use algebra to predict future performance and compare players across different sample sizes.

\(\text{Expected Goals} = \text{Shots} \times \frac{\text{Shooting \%}}{100}\)

Algebraic Relationships:

  • Solving for shots: \(\text{Shots} = \frac{\text{Goals}}{\text{Shooting \%} / 100}\)
  • Combining statistics: \(\text{Total \%} = \frac{\text{Goals}_A + \text{Goals}_B}{\text{Shots}_A + \text{Shots}_B} \times 100\)
  • Weighted averages: Used when comparing multiple game performances

Prediction Problem

If a player maintains 12% shooting and takes 250 shots in a season, expected goals = 250 × 0.12 = 30 goals.

Regression to the Mean and Sample Size Effects

Advanced statistical analysis accounts for variance, confidence intervals, and probabilistic outcomes.

\(\sigma = \sqrt{\frac{p(1-p)}{n}} \quad \text{(Standard Error of Proportion)}\)

Statistical Concepts:

  • p = true shooting percentage (as a decimal)
  • n = sample size (number of shots)
  • σ = standard error (uncertainty in the estimate)
  • 95% CI: \(p \pm 1.96\sigma\) (confidence interval)

Small Sample Bias

A player with 5 goals on 10 shots (50%) likely doesn't have a true 50% shooting percentage—the sample is too small.

Regression Analysis

NHL teams use multiple regression to model shooting % as a function of shot distance, angle, traffic, and shooter skill.

Geometry: Angles and Distances

Hockey is played on a precise geometric surface. Understanding angles, distances, and spatial relationships is crucial for strategy.

Hockey geometry visualization showing shooting angles and distances with neon overlays
Shooting Angles
Middle School
High School
Advanced

Understanding Angles to the Goal

The angle between a shooter, the puck, and the goal determines how much net the shooter can see.

Larger Angle = More Net Visible = Better Chance to Score

Example 1: Center Position

Shooting from directly in front of the goal (center ice) gives the widest angle and best scoring chance.

Example 2: Sharp Angles

Shooting from the side (near the boards) creates a sharp angle where the goalie blocks most of the net.

Example 3: The Slot

The "slot" (15-30 feet from goal, center ice) offers angles of 20-40 degrees—ideal for scoring.

Example 4: Corner Shots

From the corner (5 feet from boards, 10 feet from goal), the angle might be only 5-10 degrees—very difficult!

Example 5: Blue Line Shot

From the blue line (60 feet out, center ice), the angle is about 2×arctan(3/60) ≈ 5.7 degrees—much narrower than the slot!

Example 6: Rink Dimensions

NHL rink: 200 feet long × 85 feet wide. The goal is 6 feet wide × 4 feet tall, creating specific geometric targets.

Hockey puck with geometric shapes and distance measurements in neon colors

Trigonometry of the Goal Angle

We can calculate the exact shooting angle using trigonometric functions.

\(\theta = \arctan\left(\frac{3}{d}\right) \quad \text{(angle from side to center)}\)

Variables:

  • θ = shooting angle in radians or degrees
  • d = distance from goal in feet
  • 3 = half the goal width (6 feet wide / 2)
  • Total angle: Full viewing angle = \(2\theta\)

Example Calculation

From 30 feet away and 10 feet to the side: angle ≈ 2 × arctan(3/30) ≈ 11.3 degrees of net visible.

Optimal Position

The "slot" (center area 15-30 feet from goal) provides angles of 20-40 degrees—the sweet spot for scoring.

Example 2: Point Shot Angle

From the point (60 ft away): angle = 2×arctan(3/60) × (180/π) ≈ 5.7 degrees. Explains why point shots rarely score!

Example 3: Goalie Coverage

Goalie is ~2 feet wide. At 10 feet away, goalie blocks arctan(1/10) × (180/π) × 2 ≈ 11.4 degrees of the total angle.

Solid Angles and 3D Shooting Cones

Advanced geometry models the three-dimensional cone of possible shot trajectories.

\(\Omega = 2\pi(1 - \cos(\theta)) \quad \text{(solid angle in steradians)}\)

3D Geometric Analysis:

  • Solid angle: Measures 3D angular "area" of visible net
  • Elevation angle: Shooting high vs. low changes available targets
  • Goalie coverage: Subtracts occupied solid angle from total
  • Expected goals (xG): Models goal probability from position using multi-variable geometry

Heat Maps

NHL analytics teams create geometric heat maps showing xG from every ice position, accounting for angle, distance, and occlusion.

Defensive Geometry

Defenders position to minimize opponent shooting angles—a geometric optimization problem.

Pass Geometry
Middle School
High School
Advanced

Distance Between Two Points

To pass the puck successfully, players need to judge distances accurately.

\(\text{Distance} = \sqrt{(x_2-x_1)^2 + (y_2-y_1)^2}\)

Example 1: Simple Pass

Player at (10, 20) passes to teammate at (30, 20). Distance = √[(30-10)² + (20-20)²] = √400 = 20 feet.

Example 2: Diagonal Pass

From (0, 0) to (30, 40): Distance = √[(30)² + (40)²] = √2500 = 50 feet.

Example 3: Cross-Ice Pass

From (20, 10) to (20, 75) (across the rink width): Distance = √[(0)² + (65)²] = 65 feet straight across.

Example 4: Stretch Pass

From defensive zone (30, 40) to offensive zone (170, 50): Distance = √[(140)² + (10)²] = √19,700 ≈ 140 feet!

Example 5: Offensive Zone Triangle

Three forwards at (50,25), (50,60), and (30,42.5). Distances: √[(0)²+(35)²]=35 ft, √[(20)²+(17.5)²]≈26.5 ft, √[(20)²+(17.5)²]≈26.5 ft.

Example 6: Behind-Net Pass

From behind goal (5, 42.5) to wing (25, 70): Distance = √[(20)² + (27.5)²] = √1256.25 ≈ 35.4 feet—quick decision needed!

Distance Formula and Pythagorean Theorem

The distance formula comes directly from the Pythagorean theorem applied to coordinate geometry.

\(d = \sqrt{(x_2-x_1)^2 + (y_2-y_1)^2} \quad \text{from} \quad a^2 + b^2 = c^2\)

Breaking it Down:

  • (x₂-x₁) = horizontal distance between players
  • (y₂-y₁) = vertical distance between players
  • d = straight-line distance (hypotenuse)

Example

Player A at (20, 30) passes to Player B at (50, 50). Distance = √[(50-20)² + (50-30)²] = √[900 + 400] = √1300 ≈ 36 feet.

Pass Angle

Using \(\arctan\left[\frac{y_2-y_1}{x_2-x_1}\right]\), we can calculate the angle to aim the pass.

Vector Fields and Optimal Passing Lanes

Advanced geometric analysis models all possible passes as a vector field across the ice surface.

\(\vec{v} = (x_2-x_1, y_2-y_1) \quad ||\vec{v}|| = \sqrt{v_x^2 + v_y^2}\)

Vector Analysis:

  • Pass vector: Magnitude (speed needed) and direction
  • Intercept geometry: Defenders' positions create exclusion zones
  • Voronoi diagrams: Partition ice into control zones by player
  • Optimization: Find shortest unobstructed path between teammates

Computational Geometry

AI systems calculate thousands of potential passes per second, evaluating which lanes are open using line-segment intersection algorithms.

Triangle Offense

Three players forming a triangle maximize passing options—any two players have a direct passing lane.

Trigonometry: Sine, Cosine, and Tangent in Hockey

Trigonometric functions are essential for analyzing angles, trajectories, and the relationship between forces in hockey.

Trigonometry functions visualized through hockey puck trajectory with sine and cosine waves
Sine, Cosine, and Tangent Basics
Middle School
High School
Advanced

Understanding Right Triangles

When a player shoots at an angle, we can use right triangles to understand the geometry of the shot.

SOH-CAH-TOA: A memory trick for sin, cos, and tan!

The Basics:

  • Sine: Opposite side ÷ Hypotenuse (SOH)
  • Cosine: Adjacent side ÷ Hypotenuse (CAH)
  • Tangent: Opposite side ÷ Adjacent side (TOA)

Example 1: Shot Angle

If a player is 30 feet from goal and 15 feet to the side, tan(angle) = 15/30 = 0.5, so angle ≈ 27°.

Example 2: Height and Distance

A puck shot at 30° travels horizontally 60 feet. Its maximum height involves sin(30°) = 0.5.

Example 3: Corner Angle

From the corner (40 ft away, 28 ft to side): tan(θ) = 28/40 = 0.7, giving angle of about 35°.

Example 4: Slot Position

From the slot (25 ft away, 5 ft to side): tan(θ) = 5/25 = 0.2, giving angle of about 11° from center.

Example 5: Ramp Shot

Shooting upward at 10°: if speed is 100 ft/s, vertical component = 100×sin(10°) ≈ 17.4 ft/s, horizontal = 100×cos(10°) ≈ 98.5 ft/s.

Example 6: Bank Pass Off Boards

Puck hits boards at 30° angle. Reflection angle also 30° (angle in = angle out). Use tan(30°) ≈ 0.577 to calculate position.

Trigonometric Functions and Unit Circle

Trigonometry extends beyond right triangles to analyze any angle and rotational motion.

\(\sin(\theta) = \frac{y}{r} \quad \cos(\theta) = \frac{x}{r} \quad \tan(\theta) = \frac{y}{x} = \frac{\sin(\theta)}{\cos(\theta)}\)

Key Trigonometric Concepts:

  • sin²(θ) + cos²(θ) = 1 (Pythagorean identity)
  • Period: sin and cos repeat every 360° (2π radians)
  • Special angles: sin(30°) = 0.5, cos(60°) = 0.5, tan(45°) = 1
  • Inverse functions: arcsin, arccos, arctan find angles from ratios

Example 1: Launch Angle

A shot with initial velocity 100 ft/s at 20° has horizontal component: v_x = 100 × cos(20°) ≈ 94 ft/s.

Example 2: Vertical Component

Same shot has vertical component: v_y = 100 × sin(20°) ≈ 34 ft/s upward.

Example 3: Finding the Angle

If horizontal velocity is 80 ft/s and vertical is 60 ft/s, angle = arctan(60/80) ≈ 37°.

Example 4: Pass Direction

To pass from (20,30) to (50,60), the angle is arctan[(60-30)/(50-20)] = arctan(1) = 45°.

Example 5: Crossbar Height

From 20 ft away, to hit crossbar (4 ft high): angle = arctan(4/20) = arctan(0.2) ≈ 11.3° above horizontal.

Example 6: Deflection Angle

Puck approaches at 25°, stick redirects at 60° to ice. Change in angle: 60° - 25° = 35° deflection using trig identities.

Advanced Trigonometry and Fourier Analysis

Complex trigonometric analysis models periodic motion, oscillations, and wave-like patterns in player movement.

\(e^{i\theta} = \cos(\theta) + i\sin(\theta) \quad \text{(Euler's formula)}\)

Advanced Concepts:

  • Complex exponentials: Simplify rotational calculations
  • Trigonometric identities: sin(A+B) = sin(A)cos(B) + cos(A)sin(B)
  • Law of cosines: c² = a² + b² - 2ab·cos(C) for non-right triangles
  • Fourier series: Decompose periodic motion into sin/cos components

Example 1: Triangle Solving

Three players form a triangle. Sides 20ft, 30ft, angle between is 60°. Third side = √(400+900-2×20×30×0.5) ≈ 26.5ft.

Example 2: Skating Pattern

Player's x-position oscillates: x(t) = 40 + 20cos(πt). This periodic motion has amplitude 20ft, period 2 seconds.

Example 3: Force Resolution

A 150 lb force at 25° to ice: horizontal component = 150cos(25°) ≈ 136 lb, vertical = 150sin(25°) ≈ 63 lb.

Example 4: Stick Blade Angle

Blade at 35° to ice applies force F. Puck acceleration uses F×cos(35°) for forward motion component.

Projectile Motion with Trigonometry
Middle School
High School
Advanced

Shooting at an Angle

When a player shoots the puck at an angle (not straight ahead), the puck's path splits into horizontal and vertical motion.

Example 1: Flat Shot

A shot straight ahead (0° angle) has all its speed going forward, none going up or down.

Example 2: Rising Shot

A shot at 45° splits the speed equally between forward motion and upward motion.

Example 3: High Clear

A 60° shot goes more up than forward—useful for clearing the puck out of the defensive zone.

Example 4: Low Hard Shot

A 15° shot keeps most speed going forward with just a slight rise—perfect for top-shelf goals.

Example 5: Saucer Pass

A "saucer pass" lifts at 20° to clear sticks. If shot at 60 ft/s, it rises 60×sin(20°) ≈ 20.5 ft/s vertically.

Example 6: Wrap-Around Angle

Coming around the net at 80° angle: nearly all motion is sideways (sin(80°)≈0.98), very little forward (cos(80°)≈0.17).

Multiple hockey shot trajectories showing different launch angles with trigonometric analysis

Velocity Components and Range

Breaking velocity into components using sine and cosine lets us predict where the puck will land.

\(v_x = v_0 \cos(\theta) \quad v_y = v_0 \sin(\theta)\)
\(\text{Range} = \frac{v_0^2 \sin(2\theta)}{g} \quad \text{(for projectile from ground)}\)

Projectile Variables:

  • v₀ = initial velocity (total shot speed)
  • θ = launch angle above horizontal
  • v_x = horizontal velocity component
  • v_y = vertical velocity component
  • g = 32.2 ft/s² (gravity)

Example 1: 30° Launch

Shot at 80 ft/s at 30°: v_x = 80cos(30°) ≈ 69 ft/s, v_y = 80sin(30°) = 40 ft/s.

Example 2: Maximum Range

For a 100 ft/s shot: Range = (100)²sin(90°)/32.2 = 10,000/32.2 ≈ 310 feet at 45° angle!

Example 3: Time in Air

Shot at 60 ft/s at 40°: time = 2v_y/g = 2×(60sin40°)/32.2 ≈ 2.4 seconds in the air!

Example 4: Shallow Angle

A 10° shot at 120 ft/s: v_x = 120cos(10°) ≈ 118 ft/s forward, v_y = 120sin(10°) ≈ 21 ft/s up.

Example 5: Optimal Range Angle

To maximize distance, shoot at 45°! At 90 ft/s: Range = (90)²×sin(90°)/32.2 = 8100/32.2 ≈ 252 feet maximum possible!

Example 6: Clearing the Crossbar

Shot from 15 ft at 25°, 80 ft/s. At goal: time = 15/(80cos25°) ≈ 0.207s. Height = 80sin(25°)×0.207 - 16.1×(0.207)² ≈ 6.3 ft. Clears 4 ft bar!

Optimization and Trajectory Analysis

Calculus combined with trigonometry optimizes shot angles for specific situations.

\(\frac{d(\text{Range})}{d\theta} = \frac{2v_0^2\cos(2\theta)}{g} = 0 \implies \theta = 45°\)

Optimization Concepts:

  • Maximum range: 45° angle for projectile from level ground
  • Target optimization: Different optimal angle when shooting from height
  • Wind resistance: Changes optimal angle (typically lower than 45°)
  • Goal clearance: Minimum angle to clear goalie at given distance

Example 1: Clearing Crossbar

From 20 ft away, 4 ft height, clear 4 ft crossbar: minimum angle θ where 20tan(θ) ≥ 4, so θ ≥ arctan(0.2) ≈ 11°.

Example 2: Goalie Screen

Goalie at 6 ft away, 5 ft tall, shooter 25 ft out: angle to clear = arctan[(5-3)/(25-6)] ≈ 6° minimum.

Example 3: Dump and Chase

Optimize flip angle from 100 ft: accounting for drag, optimal angle is ≈ 35° for maximum distance with elevation.

Example 4: Top Shelf Shot

From 30 ft, hit top corner at (30, 0, 4): requires sin(θ) = 4/√(30²+4²) where total distance is hypotenuse.

Angular Velocity and Rotation
Middle School
High School
Advanced

Spinning and Rotating

When the puck spins or players turn, we measure rotation in degrees or fractions of a full circle.

Example 1: Quarter Turn

A player makes a 90° turn (quarter circle) to change direction from skating forward to skating sideways.

Example 2: Half Spin

A puck spinning at 180° (half rotation) per second completes a full rotation every 2 seconds.

Example 3: Full Rotation

A figure skater does a complete 360° spin—that's one full rotation around their axis.

Example 4: Fast Spin

A slap shot puck can spin at 1,000+ RPM (revolutions per minute)—over 16 complete rotations per second!

Example 5: Sharp Turn

A player making a 180° turn (π radians) in 0.8 seconds has angular velocity ω = π/0.8 ≈ 3.93 rad/s.

Example 6: Skating Circle

Skating a full circle (360° = 2π radians) in 4 seconds: ω = 2π/4 ≈ 1.57 rad/s (one quarter turn per second).

Angular Velocity and Radians

Angular velocity measures how fast something rotates, typically in radians per second.

\(\omega = \frac{\Delta\theta}{\Delta t} \quad \text{(angular velocity)}\)
\(v = r\omega \quad \text{(linear velocity from angular velocity)}\)

Rotation Concepts:

  • ω (omega) = angular velocity in radians/second
  • 1 revolution = 2π radians ≈ 6.28 radians
  • 1 RPM = 2π/60 ≈ 0.105 rad/s
  • r = radius from center of rotation

Example 1: Puck Spin

Puck spinning at 600 RPM = 600 × 0.105 ≈ 63 rad/s angular velocity.

Example 2: Edge Speed

Puck (radius 1.5 inches = 0.125 ft) spinning at 60 rad/s: edge speed = 0.125 × 60 = 7.5 ft/s!

Example 3: Player Turn

Player turns 90° (π/2 radians) in 0.5s: angular velocity = (π/2)/0.5 ≈ 3.14 rad/s.

Example 4: Curved Skating

Skating in a 10 ft radius circle at 20 ft/s: angular velocity = v/r = 20/10 = 2 rad/s.

Example 5: Stick Blade Rotation

During a wrist shot, blade rotates 45° (π/4 rad) in 0.05s: ω = (π/4)/0.05 ≈ 15.7 rad/s—very fast motion!

Example 6: Converting RPM

Puck spinning at 720 RPM: ω = 720 × (2π/60) = 720 × 0.1047 ≈ 75.4 rad/s angular velocity.

Moment of Inertia and Rotational Dynamics

Advanced rotational physics uses trigonometry to analyze spinning objects and torque.

\(L = I\omega \quad \tau = I\alpha \quad KE_{rot} = \frac{1}{2}I\omega^2\)

Rotational Dynamics:

  • L = angular momentum (conserved in isolated systems)
  • I = moment of inertia (resistance to rotation)
  • τ (tau) = torque (rotational force)
  • α (alpha) = angular acceleration
  • Puck: I = ½mr² for solid disk

Example 1: Spinning Puck

6 oz puck (0.17 kg), radius 0.038 m: I = ½(0.17)(0.038)² ≈ 0.00012 kg·m². At 60 rad/s: L = 0.0072 kg·m²/s.

Example 2: Friction Torque

Ice friction creates torque τ = F×r. If F = 0.01 N at edge, τ = 0.01 × 0.038 ≈ 0.00038 N·m decelerates spin.

Example 3: Magnus Effect

Spinning puck at ω = 50 rad/s creates lift force perpendicular to velocity. Curve depends on sin(θ) where θ is spin axis angle.

Example 4: Rotational Energy

Puck spinning at 60 rad/s: KE_rot = ½(0.00012)(60)² ≈ 0.22 J—small compared to translational KE but affects stability!

Calculus: Motion and Change

Calculus helps us understand how things change over time—essential for analyzing acceleration, trajectory curves, and optimization.

Velocity and Acceleration
Middle School
High School
Advanced

Speed vs. Acceleration

Speed tells us how fast something is moving. Acceleration tells us how fast the speed is changing.

Acceleration = Change in Speed / Time

Example 1: Constant Speed

If a puck slides at 50 ft/s the whole way, acceleration is zero—no change in speed.

Example 2: Rapid Acceleration

A slap shot accelerates the puck from 0 to 147 ft/s (100 MPH) in 0.1 seconds: a = 1,470 ft/s²!

Example 3: Ice Friction

A puck at 100 ft/s decelerates at -0.3 ft/s² due to friction. After 10s: v = 100 - (0.3 × 10) = 97 ft/s.

Example 4: Wrist Shot

Wrist shots have lower acceleration than slap shots: 0 to 88 ft/s (60 MPH) in 0.05s gives a = 1,760 ft/s².

Derivatives: Rates of Change

Calculus introduces the derivative, which gives us instantaneous rates of change.

\(v(t) = \frac{dx}{dt} \quad\quad a(t) = \frac{dv}{dt} = \frac{d^2x}{dt^2}\)

Understanding Derivatives:

  • x(t) = position of puck at time t
  • \(v(t) = \frac{dx}{dt}\) = velocity (first derivative of position)
  • \(a(t) = \frac{dv}{dt}\) = acceleration (first derivative of velocity, second derivative of position)
\(\text{Example: } x(t) = 5t^2 + 2t \rightarrow v(t) = 10t + 2 \rightarrow a(t) = 10\)

Constant Acceleration

If a = 10 ft/s², the puck gains 10 ft/s of velocity every second.

Finding Maximum Velocity

When acceleration becomes zero (a = 0), velocity reaches its maximum or minimum.

Vector Calculus and Trajectory Optimization

Multi-dimensional calculus models the puck's motion in both horizontal directions and vertical (when airborne).

\(\vec{r}(t) = x(t)\hat{i} + y(t)\hat{j} + z(t)\hat{k} \quad \vec{v}(t) = \frac{d\vec{r}}{dt} \quad \vec{a}(t) = \frac{d\vec{v}}{dt}\)

Multi-variable Calculus:

  • Position vector: r⃗(t) tracks puck in 3D space
  • Velocity components: v_x, v_y, v_z can change independently
  • Acceleration due to gravity: a_z = -g = -32.2 ft/s²
  • Air resistance: F_drag = -½ρv²C_dA (quadratic drag model)
\(\text{Differential equation: } m\frac{d^2\vec{r}}{dt^2} = \vec{F}_{gravity} + \vec{F}_{drag} + \vec{F}_{friction}\)

Numerical Integration

NHL tracking systems use Runge-Kutta methods to solve these differential equations and predict puck trajectories in real-time.

Optimization

Calculus of variations finds the optimal launch angle and velocity to maximize shooting accuracy given constraints.

Puck Trajectory
Middle School
High School
Advanced

Step-by-Step: How the Puck Moves

Let's break down a puck's flight into small time steps to see exactly how it moves. We'll shoot at 30° with initial speed of 80 ft/s.

Initial Conditions:

  • Total speed: 80 ft/s
  • Angle: 30° above horizontal
  • Horizontal speed (v_x): 80 × cos(30°) = 80 × 0.866 ≈ 69.3 ft/s
  • Vertical speed (v_y): 80 × sin(30°) = 80 × 0.5 = 40 ft/s
  • Starting height: 3 feet (stick height)

📊 Part 1: Iteration in X-Dimension (Horizontal Motion)

Simple pattern: Horizontal speed stays constant (no gravity sideways!)

Each 0.1 second:
• t = 0.0s: x = 0 ft
• t = 0.1s: x = 0 + 69.3×0.1 = 6.93 ft
• t = 0.2s: x = 6.93 + 69.3×0.1 = 13.86 ft
• t = 0.3s: x = 13.86 + 69.3×0.1 = 20.79 ft
• t = 0.4s: x = 20.79 + 69.3×0.1 = 27.72 ft
• t = 0.5s: x = 27.72 + 69.3×0.1 = 34.65 ft
→ Each step adds the same distance: 6.93 ft

Key Insight:

The horizontal position increases linearly with time because there's no horizontal force acting on the puck. Each time step adds exactly the same distance!

📉 Part 2: Iteration in Y-Dimension (Vertical Motion)

Complex pattern: Gravity pulls down, so vertical speed decreases each step!

Gravity pulls down at 32.2 ft/s²
Each 0.1 second, speed changes by: -32.2×0.1 = -3.22 ft/s

• t = 0.0s: y = 3 ft, v_y = 40 ft/s ↑
• t = 0.1s: y = 3 + 40×0.1 - ½×32.2×(0.1)² = 6.84 ft, v_y = 40 - 3.22 = 36.78 ft/s ↑
• t = 0.2s: y = 6.84 + 36.78×0.1 - ½×32.2×(0.1)² = 10.36 ft, v_y = 36.78 - 3.22 = 33.56 ft/s ↑
• t = 0.3s: y = 10.36 + 33.56×0.1 - ½×32.2×(0.1)² = 13.57 ft, v_y = 33.56 - 3.22 = 30.34 ft/s ↑
• t = 0.4s: y = 13.57 + 30.34×0.1 - ½×32.2×(0.1)² = 16.45 ft, v_y = 30.34 - 3.22 = 27.12 ft/s ↑
• t = 0.5s: y = 16.45 + 27.12×0.1 - ½×32.2×(0.1)² = 19.02 ft, v_y = 27.12 - 3.22 = 23.90 ft/s ↑
→ Height increases but by LESS each step (slowing down!)

Key Insight:

The vertical position follows a parabolic curve because gravity continuously accelerates the puck downward. The puck rises quickly at first, slows down, reaches a peak, then falls!

🎯 Part 3: Combining X and Y - The Complete Path

The puck's position at each moment combines both dimensions:

Time (s) X Position (ft) Y Height (ft) Position (x, y)
0.0 0.00 3.00 (0, 3)
0.1 6.93 6.84 (6.93, 6.84)
0.2 13.86 10.36 (13.86, 10.36)
0.3 20.79 13.57 (20.79, 13.57)
0.4 27.72 16.45 (27.72, 16.45)
0.5 34.65 19.02 (34.65, 19.02)

✨ These points trace out a curved path—the puck's trajectory!

Key Insight:

When we plot the X and Y positions together, we see the classic parabolic trajectory! The puck moves forward at a constant rate while simultaneously rising and falling due to gravity.

🔬 Part 4: From Iterations to Calculus

What happens if we make the time steps smaller and smaller?

  • ✓ With 0.1s steps: we calculated 5 positions
  • ✓ With 0.01s steps: we'd calculate 50 positions (more accurate!)
  • ✓ With 0.001s steps: we'd calculate 500 positions (even better!)
  • With INFINITELY small steps: That's CALCULUS!

🎓 Calculus finds the exact position at ANY instant by using a mathematical limit as the time step approaches zero. Instead of adding up tiny changes, we use derivatives and integrals to get perfect accuracy!

From Iterations to Continuous Functions with Calculus

Instead of calculating position at individual time steps, calculus gives us formulas that work for ANY time!

🔄 Review: What We Learned from Iterations

  • X-direction: Position increases by v_x × Δt each step
  • Y-direction: Position increases by v_y × Δt, but v_y decreases by g × Δt each step
  • The pattern: As steps get smaller, accuracy improves
\(x(t) = x_0 + v_x \cdot t \quad \text{(horizontal: constant velocity)}\)
\(y(t) = y_0 + v_y \cdot t - \frac{1}{2}gt^2 \quad \text{(vertical: constant acceleration)}\)

Understanding the Formulas:

  • x(t) = x₀ + v_x·t — This is just "distance = speed × time" because horizontal speed is constant!
  • y(t) = y₀ + v_y·t - ½gt² — The -½gt² term accounts for gravity pulling down continuously
  • v₀ (or v_y) = initial vertical velocity component
  • g = 32.2 ft/s² (acceleration due to gravity)

📊 Comparing Iteration vs. Calculus Formula

Same example: 80 ft/s at 30°, starting at 3 ft height

At t = 0.5 seconds:

Iteration method: x = 34.65 ft, y = 19.02 ft (from adding 5 steps)

Calculus formula:
x(0.5) = 0 + 69.3×0.5 = 34.65 ft
y(0.5) = 3 + 40×0.5 - ½×32.2×(0.5)² = 3 + 20 - 4.025 = 18.975 ft

✓ Same answer! But calculus gives it instantly, no iteration needed.

🎯 Calculus Power: Velocity from Position

Taking the derivative gives us velocity at any moment:

Position: x(t) = 69.3t
Velocity: v_x(t) = dx/dt = 69.3 ft/s (constant!)

Position: y(t) = 3 + 40t - 16.1t²
Velocity: v_y(t) = dy/dt = 40 - 32.2t ft/s (decreasing!)

At t = 0.5s: v_y = 40 - 32.2×0.5 = 23.9 ft/s
At t = 1.24s: v_y = 40 - 32.2×1.24 = 0 ft/s (peak height!)
At t = 2.0s: v_y = 40 - 32.2×2.0 = -24.4 ft/s ↓ (falling!)

📈 Calculus Power: Acceleration from Velocity

Taking the derivative again gives us acceleration:

Velocity: v_x(t) = 69.3
Acceleration: a_x(t) = dv_x/dt = 0 ft/s² (no horizontal force!)

Velocity: v_y(t) = 40 - 32.2t
Acceleration: a_y(t) = dv_y/dt = -32.2 ft/s² (gravity!)

This acceleration is constant—gravity always pulls down with the same force!

🎓 Finding Max Height with Calculus

The puck reaches maximum height when vertical velocity = 0:

Set v_y(t) = 0:

40 - 32.2t = 0

t = 40/32.2 ≈ 1.24 seconds

Substitute into y(t):

y(1.24) = 3 + 40(1.24) - 16.1(1.24)² = 3 + 49.6 - 24.77 ≈ 27.8 feet

The puck reaches 27.8 feet at its peak!

🌟 The Big Idea

Iterations show us HOW calculus works by adding up tiny changes. Calculus formulas give us the EXACT answer by taking the mathematical limit as those changes become infinitely small. This is the power of derivatives and integrals—they turn step-by-step arithmetic into instant, precise calculations!

Advanced Calculus: Differential Equations and Numerical Methods

When forces change continuously (like drag that depends on velocity), we need differential equations and return to iteration methods—but much more sophisticated ones!

🔄 Full Circle: Back to Iterations (But Smarter!)

Simple projectile motion (constant gravity) has exact formulas. But add air resistance, and suddenly we need numerical methods—fancy iterations!

\(m\frac{dv_x}{dt} = -\frac{1}{2}\rho v^2 C_d A \cdot \frac{v_x}{v} \quad m\frac{dv_y}{dt} = -mg - \frac{1}{2}\rho v^2 C_d A \cdot \frac{v_y}{v}\)

Advanced Trajectory Modeling:

  • ρ = air density (~0.075 lb/ft³ in ice rinks)
  • C_d = drag coefficient (≈0.5 for hockey puck)
  • A = cross-sectional area (puck: π × 1.5² ≈ 7 in²)
  • v = √(v_x² + v_y²) = total velocity magnitude
  • Problem: Drag force depends on v, which depends on position, which depends on v... it's circular!

🔬 Euler's Method: Simple Numerical Integration

The basic idea: Use derivatives to predict the next small step

Algorithm:
1. Start with initial position (x₀, y₀) and velocity (v_x₀, v_y₀)
2. Calculate forces: F_x = -drag_x, F_y = -mg - drag_y
3. Calculate acceleration: a_x = F_x/m, a_y = F_y/m
4. Update velocity: v_x₁ = v_x₀ + a_x×Δt, v_y₁ = v_y₀ + a_y×Δt
5. Update position: x₁ = x₀ + v_x₀×Δt, y₁ = y₀ + v_y₀×Δt
6. Repeat from step 2 with new values!

→ This is iteration with physics—calculus tells us the rate of change at each instant!

🎯 Runge-Kutta Method (RK4): The Gold Standard

More accurate than Euler—samples the derivative at multiple points per step

The RK4 approach:

  1. Calculate slope k₁ at the beginning of the interval
  2. Calculate slope k₂ at the midpoint using k₁
  3. Calculate slope k₃ at the midpoint using k₂
  4. Calculate slope k₄ at the end using k₃
  5. Weighted average: Δy = (k₁ + 2k₂ + 2k₃ + k₄)/6 × Δt

Result: With the same step size, RK4 is ~100× more accurate than Euler!

📊 Example: Comparing Methods for Real Puck Flight

80 ft/s shot at 30° with air resistance included:

Method Step Size Max Height Accuracy
No drag formula N/A 27.8 ft Exact (but unrealistic)
Euler 0.1s 26.2 ft ±5% error
Euler 0.01s 26.45 ft ±0.5% error
RK4 0.1s 26.48 ft ±0.01% error
NHL tracking systems use RK4 with Δt = 0.001s for real-time trajectory prediction!

🌪️ Adding Complexity: Magnus Effect (Spinning Puck)

When the puck spins, it curves! This adds another force to our differential equations:

Magnus lift force: F_L = ½ρv²C_L A
Direction: perpendicular to velocity and spin axis

Enhanced differential equations:
m(dv_x/dt) = -drag_x + Magnus_x
m(dv_y/dt) = -mg - drag_y + Magnus_y

→ Now we need 3D vector calculus and numerical integration!

🎓 The Beautiful Connection

We started with simple iterations to understand motion step-by-step.
Calculus gave us exact formulas for simple cases (constant forces).
But real physics brings us back to iteration—sophisticated numerical methods that use calculus (derivatives) at each tiny step to handle complex, changing forces.

This is how modern physics works: analytical calculus provides insight, numerical calculus provides solutions!

💻 Real-World Application

NHL player tracking systems (like NHL Edge IQ) use these exact methods to:

  • ✓ Predict puck trajectory in real-time (RK4 integration, 1000 Hz)
  • ✓ Calculate expected goals (xG) by simulating thousands of shots
  • ✓ Optimize goalie positioning based on trajectory predictions
  • ✓ Detect when shots are tipped/deflected by analyzing acceleration discontinuities

Physics: Forces and Energy

Physics principles govern every aspect of hockey, from collisions to friction to energy transfer.

Kinetic Energy
Middle School
High School
Advanced

Energy of Motion

Kinetic energy is the energy an object has because it's moving. Faster = more energy!

\(KE = \frac{1}{2} \times \text{mass} \times \text{speed}^2\)

Example 1: Puck Energy

A 6-oz puck at 100 ft/s: KE = ½ × (0.375 lb) × (100)² ÷ 32.2 ≈ 58 ft-lb of energy!

Example 2: Doubling Speed

Same puck at 200 ft/s has KE = ½ × 0.375 × (200)² ÷ 32.2 ≈ 233 ft-lb—four times more energy!

Example 3: Player Collision

A 200-lb player at 20 MPH (29.3 ft/s): KE = ½ × 200 × (29.3)² ÷ 32.2 ≈ 2,665 ft-lb—huge impact!

Example 4: Goalie Save

A goalie stopping a 100 MPH shot (147 ft/s) absorbs about 125 ft-lb of energy in their catching glove!

Work-Energy Theorem

The work done on an object equals its change in kinetic energy.

KE = ½mv²     W = ΔKE = ½m(v_f² - v_i²)

Key Concepts:

  • m = mass (puck: 6 oz = 0.375 lb = 0.17 kg)
  • v = velocity in ft/s or m/s
  • W = work done (force × distance)
  • Energy units: Joules (metric) or foot-pounds (imperial)

Example Calculation

A 6-oz puck at 100 MPH: KE = ½(0.17 kg)(44.7 m/s)² ≈ 170 Joules—enough to seriously injure someone!

Power Output

If a player accelerates a puck to 100 MPH in 0.1 seconds, power = ΔKE/Δt = 170J/0.1s = 1,700 Watts!

Energy Conservation and Loss Mechanisms

Real systems lose energy to friction, deformation, and other dissipative forces.

E_total = KE + PE + U_internal     dE_total/dt = -P_friction

Energy Analysis:

  • KE = ½mv² (translational kinetic energy)
  • KE_rot = ½Iω² (rotational kinetic energy, I = moment of inertia)
  • PE = mgh (gravitational potential energy)
  • U_internal = heat, deformation, vibration
  • P_friction = F_friction × v (power loss to friction)

Coefficient of Restitution

e = √(KE_after/KE_before) for collisions. Hockey pucks on ice: e ≈ 0.6 (60% energy retained in bounce).

Stick Flex Energy

Elastic potential energy stored in stick flex: U = ½kx². When released, converts to puck KE—adding 10-15 MPH to shot velocity.

Friction and Momentum
Middle School
High School
Advanced

Why the Puck Slows Down

Friction is a force that opposes motion. Even though ice is slippery, friction still slows the puck.

Example 1: Puck Momentum

A 6-oz (0.375 lb) puck at 100 ft/s has momentum = 0.375 × 100 = 37.5 lb·ft/s.

Example 2: Player Momentum

A 180-lb player skating at 20 ft/s has momentum = 180 × 20 = 3,600 lb·ft/s—nearly 100× more than the puck!

Example 3: Head-on Collision

Player A (200 lb, 15 ft/s) hits Player B (180 lb, -10 ft/s). Total momentum before = 3,000 + (-1,800) = 1,200 lb·ft/s.

Example 4: Friction Slowing

Ice friction force of 0.01 lb on a 0.375 lb puck creates deceleration of 0.01/0.375 × 32.2 ≈ 0.86 ft/s².

Newton's Laws and Momentum

Newton's laws govern all motion in hockey.

F = ma     p = mv     F = dp/dt

Newton's Laws Applied:

  • First Law: Puck continues moving until friction/collision stops it
  • Second Law: F = ma relates force, mass, and acceleration
  • Third Law: When stick hits puck, puck pushes back on stick with equal force
  • Friction force: F_f = μ × N (μ ≈ 0.01-0.02 for ice)

Deceleration from Friction

a = F_f/m = μg ≈ 0.32 ft/s². A puck at 100 ft/s takes about 300 feet to stop!

Impulse

J = FΔt = Δp. A stick in contact with puck for 0.02s delivering 100 lb force: Δp = 2 lb·s, changing velocity dramatically.

Tribology and Momentum Transfer

Advanced physics of ice friction involves thin water layers, pressure melting, and viscoelastic deformation.

F_f = μ(v,T,P) × N     (velocity, temperature, pressure dependent)

Complex Friction Models:

  • Pressure melting: Puck pressure creates thin water layer, reducing friction
  • Temperature dependence: Optimal ice temperature: 15-20°F for minimum friction
  • Velocity dependence: μ increases slightly with speed due to viscous drag in water layer
  • Surface roughness: Ice resurfacing (Zamboni) minimizes asperity contact
Momentum conservation: Σp⃗_before = Σp⃗_after    (for all collisions)

Two-Body Collisions

For player-player collision: m₁v₁ + m₂v₂ = m₁v₁' + m₂v₂'. Combined with energy equation gives post-collision velocities.

Center of Mass

In multi-player collisions, center of mass velocity: v_cm = (m₁v₁ + m₂v₂)/(m₁ + m₂) remains constant.

Interactive: Momentum Collision

Before Collision:
Total Momentum = 1200 lb·ft/s
After Collision (elastic):
Player 1 velocity = -5.8 ft/s
Player 2 velocity = 13.2 ft/s
After Collision:
Total Momentum = 1200 lb·ft/s
Momentum is conserved!

Putting It All Together

Real hockey scenarios involve multiple mathematical concepts working simultaneously.

Complete Shot Analysis

1. Geometry: Shooting Angle

Player evaluates angle to goal using trigonometry: θ = 2·arctan(3/d). Finds optimal shooting position.

2. Physics: Force Application

Player flexes stick (elastic PE stored). Releases: F = ma accelerates puck. Kinetic energy = ½mv².

3. Calculus: Trajectory

Differential equations model path: x(t), y(t), z(t). Gravity and drag affect motion continuously.

4. Algebra: Outcome Probability

Statistical model: P(goal) = f(angle, distance, velocity, goalie position). Uses regression from thousands of shots.

Mathematics Makes Hockey Possible

Every aspect of hockey—from the design of the rink to the curve of a slap shot—relies on mathematical principles. Players who understand these concepts, even intuitively, gain a competitive advantage. And for fans and analysts, mathematics reveals the hidden beauty and complexity of the game.