STATICS • CENTROIDS AND DISTRIBUTED LOADS

Centroid via Integration — Compute centroid using integration for simple shapes (intro)

Locate the geometric balance point of planar regions using first-moment integrals.

Historical Context & Motivation

The idea of a centroid — the geometric center of a shape — has roots that reach back to antiquity. Archimedes of Syracuse was arguably the first to formalize the concept when he studied the equilibrium of planar laminae and proved that a triangle's centroid lies at the intersection of its medians. For centuries, engineers relied on geometric decomposition and symmetry arguments to locate centroids, but these techniques falter when shapes cannot be neatly subdivided into rectangles, triangles, or circles. The invention of calculus in the seventeenth century by Newton and Leibniz provided the analytical machinery to handle arbitrary curves and regions, transforming centroid calculations from clever geometry into systematic integration problems. Today, every structural engineer, mechanical designer, and civil engineer uses centroid calculations — whether computing beam deflections, analyzing distributed loads, or sizing cross-sections — making integration-based centroid determination one of the foundational skills of engineering statics.

~250 BC
Archimedes — Centers of Gravity
Archimedes establishes the principle that every planar figure has a unique center of gravity and derives the centroid of triangles, parabolic segments, and other simple shapes using the method of exhaustion.
1687
Newton's Principia
Isaac Newton publishes the Principia Mathematica, formalizing calculus and providing the integral framework that would later be applied to centroid and moment-of-inertia calculations.
1788
Lagrange — Analytical Mechanics
Joseph-Louis Lagrange publishes Mécanique analytique, unifying statics and dynamics under a single variational principle and heavily employing integration for distributed mass and area properties.
1900s
Modern Engineering Statics
Centroid calculations via integration become a standard topic in engineering curricula worldwide, directly linking to beam analysis, structural design, and finite-element methods.

The central question this lesson addresses is straightforward yet powerful: given a planar region bounded by known curves, how do we systematically compute the exact coordinates of its centroid using definite integrals? Mastering this technique provides the foundation for computing first moments of area, which in turn are essential for analyzing bending stresses, shear flows, and distributed load resultants throughout your engineering career.

Core Principles & Definitions

Before diving into the integration procedure, it is essential to establish the foundational ideas that underpin centroid calculations. The centroid of a planar area is the point (x̄, ȳ) at which the entire area could be concentrated without changing the first moment of area about any axis. For a homogeneous plate of uniform thickness, the centroid coincides with the center of mass, making it the natural balance point. The integration approach works by summing the contributions of infinitesimally thin strips — each strip's area element dA multiplied by its distance from a reference axis — and then dividing by the total area.

1

First Moment of Area

The integral Qx = ∫ y dA (or Qy = ∫ x dA) quantifies how area is distributed relative to an axis. It has units of length³ and is the numerator of the centroid formula.
2

Differential Area Element (dA)

The infinitesimal slice of area used in integration. For vertical strips, dA = f(x) dx; for horizontal strips, dA = g(y) dy. Choosing the correct element simplifies the integral dramatically.
3

Centroid Coordinates (x̄, ȳ)

Defined as x̄ = (1/A) ∫ x̃ dA and ȳ = (1/A) ∫ ỹ dA, where x̃ and ỹ are the coordinates of the centroid of each differential element, and A is the total area.
4

Symmetry Principle

If a region has an axis of symmetry, the centroid lies on that axis. This can immediately fix one coordinate, cutting the required computation in half.
KEY TAKEAWAY
Think of the centroid as the balance point of a cardboard cutout. If you cut a shape from uniform cardboard and tried to balance it on the tip of a pencil, the pencil tip would need to be placed exactly at the centroid. Integration lets us compute that balance point precisely by summing up the 'moment' contribution of every infinitesimal strip of cardboard across the entire region.

Visual Explanation — The Differential Strip Method

The diagram below illustrates the fundamental setup for computing the centroid of a region bounded by a curve y = f(x) and the x-axis, from x = a to x = b. A representative vertical strip of width dx is shown at an arbitrary position x. The strip has height f(x), so its differential area is dA = f(x) dx. The centroid of this thin rectangular strip is located at coordinates (x, f(x)/2), since the strip is centered vertically at half its own height. By integrating the product of each strip's centroid coordinate and its area over the interval [a, b], we obtain the first moments of area needed for x̄ and ȳ.

A representative vertical strip of width dx is highlighted at position x. The strip's area is dA = f(x) dx and its own centroid sits at (x, f(x)/2). The faint strips in the background illustrate how the entire region is partitioned.
💡 Choosing Your Strip Orientation
Vertical strips (integrating with respect to x) are natural when the top and bottom boundaries are given as y = f(x). Horizontal strips (integrating with respect to y) are preferable when the left and right boundaries are expressed as x = g(y). The key is to ensure that every strip spans the full width or height of the region without requiring a piecewise description.

Mathematical Framework

We now formalize the centroid formulas for a planar region of total area A. These expressions are direct applications of the definition: the centroid coordinate equals the first moment of area divided by the total area. The integrals below assume that the region is bounded above by y = f(x) and below by the x-axis over the interval [a, b]. Analogous forms exist for regions bounded by two curves or expressed in terms of y.

TOTAL AREA
A = ∫ₐᵇ f(x) dx
A = total area of the region; f(x) = height of the curve at position x; a, b = limits of integration along the x-axis.
CENTROID x-COORDINATE
x̄ = (1/A) ∫ₐᵇ x · f(x) dx
Each vertical strip is located at x-position x, with area dA = f(x) dx. The integrand x · f(x) is the first moment of each strip about the y-axis. Dividing by A yields the area-weighted average x-position.
CENTROID y-COORDINATE
ȳ = (1/A) ∫ₐᵇ ½ [f(x)]² dx
The centroid of each vertical strip sits at height f(x)/2. Multiplying this by the strip's area dA = f(x) dx gives ½ [f(x)]² dx. Summing and dividing by A yields ȳ.
📐 Region Between Two Curves
When the region is bounded above by y = f(x) and below by y = g(x), the formulas become: dA = [f(x) − g(x)] dx, x̃ = x, and ỹ = [f(x) + g(x)] / 2. Substitute these into the standard centroid integrals accordingly.

Notice the pattern: the x̄ integral weights each strip's area by its x-coordinate, while the ȳ integral weights each strip's area by the strip's own vertical centroid. This structure — the first moment divided by the zeroth moment (total area) — is identical to computing a weighted average. Understanding this analogy with weighted averages provides a strong conceptual anchor: the centroid is simply the area-weighted average position of all infinitesimal elements composing the region.

Vertical vs. Horizontal Strip Elements

Engineering problems frequently present regions where one strip orientation is clearly more convenient than the other. The following diagram compares vertical and horizontal strip setups for the same parabolic region, highlighting how each strip element is defined and where its centroid is located. Selecting the correct strip orientation before writing the integrals can mean the difference between a clean, one-step integration and a messy, multi-part calculation.

Left: vertical strips for the parabola y = √(4x) yield dA = √(4x) dx. Right: horizontal strips for the same curve rewritten as x = y²/4 yield dA = [h − y²/4] dy. The pink dots mark the centroid of each representative strip.
Comparison of vertical and horizontal strip formulations
PropertyVertical Strip (dx)Horizontal Strip (dy)
Integration variablexy
dAf(x) dx[x_right − x_left] dy
x̃ (moment arm for x̄)x(x_right + x_left) / 2
ỹ (moment arm for ȳ)f(x) / 2y
Best when…Top/bottom curves are functions of xLeft/right curves are functions of y

Worked Example — Quarter-Circle Region

Find the centroid (x̄, ȳ) of the region in the first quadrant bounded by the circle x² + y² = r² (a quarter-circle of radius r). We expect, by symmetry, that x̄ = ȳ.

Centroid of a Quarter-Circle
1
Step 1 — Set Up the Region and Choose StripsThe quarter-circle occupies 0 ≤ x ≤ r with the upper boundary y = √(r² − x²). Using vertical strips, dA = √(r² − x²) dx. The strip centroid is at (x, √(r² − x²)/2).
2
Step 2 — Compute Total AreaA = ∫₀ʳ √(r² − x²) dx. This integral equals the area of a quarter-circle, so A = πr²/4.
A = πr²/4
3
Step 3 — Compute x̄x̄ = (1/A) ∫₀ʳ x · √(r² − x²) dx. Use the substitution u = r² − x², du = −2x dx. When x = 0, u = r²; when x = r, u = 0. The integral becomes (1/A) · (−½) ∫_{r²}^{0} √u du = (1/A) · (½) ∫₀^{r²} u^{1/2} du = (1/A) · (½) · [⅔ u^{3/2}]₀^{r²} = (1/A) · (½) · (⅔) · r³ = (1/A) · r³/3.
x̄ = (r³/3) / (πr²/4) = 4r / (3π)
4
Step 4 — Compute ȳȳ = (1/A) ∫₀ʳ ½ (r² − x²) dx = (1/A) · ½ [r²x − x³/3]₀ʳ = (1/A) · ½ · (r³ − r³/3) = (1/A) · ½ · (2r³/3) = (1/A) · r³/3. By symmetry, this is identical to the x̄ numerator.
ȳ = 4r / (3π)
5
Step 5 — State Final Result and Verify SymmetryThe centroid of a quarter-circle of radius r is located at (4r/(3π), 4r/(3π)). For example, if r = 6 in., then x̄ = ȳ = 4(6)/(3π) ≈ 2.546 in. The result lies on the line y = x, consistent with the 45° symmetry axis of a quarter-circle.
(x̄, ȳ) = (4r/(3π), 4r/(3π))

Integration vs. Composite-Shape Methods

Engineering practice offers two principal approaches to finding centroids: the integration method presented in this lesson, and the composite-shape (tabular) method, which breaks a complex area into simple sub-areas whose centroids are known from standard tables. Each has distinct strengths and limitations, and a practicing engineer needs to be proficient in both.

Integration vs. composite-shape centroid methods
CriterionIntegration MethodComposite-Shape Method
ApplicabilityAny region with integrable boundaries — curves, polynomials, trig functionsRegions decomposable into rectangles, triangles, semicircles, etc.
AccuracyExact (closed-form solution when antiderivative exists)Exact for combinations of standard shapes; approximate otherwise
SpeedSlower — requires setting up and evaluating integralsFaster for common cross-sections; look up tabulated centroids
Prerequisite knowledgeIntegral calculus, substitution, partial fractionsCentroid table values, algebraic bookkeeping
Handles cutouts / holesYes — adjust limits or subtract integralsYes — subtract negative area terms
⚖️ WHEN TO USE WHICH
Think of the composite-shape method as assembling LEGO bricks — fast and effective when your shape is built from standard pieces. The integration method is like sculpting clay — it handles any arbitrary contour but requires more mathematical effort. In practice, many real cross-sections (I-beams, channels, angles) are best handled by the composite approach, while curved or continuously varying profiles demand integration.

Connection to Advanced Theory

The centroid integration you are learning now is the first rung on a ladder of increasingly sophisticated area properties. The second moment of area (moment of inertia) replaces the linear moment arm (x or y) with a squared moment arm (x² or y²), and it governs bending stiffness and stress distribution in beams. The product of inertia uses the product xy as the moment arm and is essential for analyzing unsymmetric bending. Mastering centroid integrals now ensures you have the procedural and conceptual foundation for all of these higher-order calculations.

Hierarchy of area-moment integrals
QuantityIntegral FormPhysical Significance
Area (zeroth moment)A = ∫ dATotal size of the region
First moment (Q)Q = ∫ x̃ dA or ∫ ỹ dALocates the centroid; used in shear stress (VQ/It)
Second moment (I)I = ∫ y² dAGoverns bending rigidity (EI) and normal stress (My/I)
Product of inertia (Ixy)I_xy = ∫ xy dARequired for principal axes and unsymmetric bending

Additionally, the centroid concept extends beyond two dimensions. In three-dimensional solid mechanics and finite-element analysis, centroids of volumes, surface areas, and line elements are computed using triple, surface, and line integrals respectively. The parallel-axis (Steiner's) theorem, which you will encounter when computing moments of inertia of composite sections, explicitly requires the centroid location — another reason why accurate centroid calculation is so fundamental.

Practice Problems

PROBLEM 1CONCEPTUAL
A region is symmetric about both the x-axis and the y-axis. Without performing any integration, what can you conclude about its centroid? Explain why the first-moment integrals confirm your answer.
PROBLEM 2BASIC CALCULATION
Find the centroid (x̄, ȳ) of the triangular region bounded by y = 0, x = 0, and y = 3 − (3/4)x (a right triangle with base 4 and height 3). Use vertical strips.
PROBLEM 3INTERMEDIATE
Determine x̄ for the region bounded by y = x² and y = 4 (i.e., the area between the parabola and the horizontal line). Use horizontal strips.
PROBLEM 4APPLIED
A retaining wall cross-section is modeled as the region under y = 2√x from x = 0 to x = 9 ft. A distributed soil pressure acts horizontally at the centroid height. Find ȳ to determine the line of action of this resultant force.
PROBLEM 5CRITICAL THINKING
Consider the family of curves y = xⁿ (n > 0) over the region bounded by y = xⁿ, y = 0, and x = 1. Derive general expressions for x̄(n) and ȳ(n), and analyze their behavior as n → 0⁺ and n → ∞. What do the limiting cases represent geometrically?

Lesson Summary

The centroid of a planar region is computed by dividing the first moment of area by the total area. Using vertical strips, x̄ = (1/A) ∫ x f(x) dx and ȳ = (1/A) ∫ ½ [f(x)]² dx; using horizontal strips, analogous forms with g(y) apply. The choice of strip orientation should match the form in which the boundary curves are given, and symmetry should always be exploited to eliminate unnecessary integrals.

This integration method handles any region with integrable boundaries and provides the exact centroid location, in contrast to the composite-shape method which is limited to combinations of standard shapes. The centroid calculation is the foundation for higher-order area properties — the moment of inertia and product of inertia — and is essential for locating resultant forces from distributed loads in beam and structural analysis.

Varsity Tutors • Statics • Centroid via Integration — Compute centroid using integration for simple shapes (intro)