Calculus III: Applications of Double Integrals

Double integrals allow us to move beyond calculating area and volume to solving foundational problems in physics, engineering, and probability. By integrating a function over a two-dimensional region, we can determine the mass of an object with variable composition, locate its balance point, analyze its resistance to rotation, and find average values or probabilistic outcomes. These applications are critical for modeling real-world systems, from designing machine components to assessing material stress.

Mass of a Lamina with Variable Density

The most direct physical application of the double integral is calculating the mass of a thin plate, or lamina. If the lamina occupies a region $R$ in the $xy$-plane and its density (mass per unit area) at a point $(x,y)$ is given by $\rho (x,y)$, the mass is not simply density times area because the density varies. Instead, we divide $R$ into small sub-regions, approximate the mass of each as density times area, and sum. In the limit, this Riemann sum becomes a double integral.

The mass $m$ of the lamina is: $$m={\iint }_R\rho (x,y)dA$$

For example, consider a rectangular plate defined by $0\le x\le 2$, $0\le y\le 1$, with a density function $\rho (x,y)=1+xy$ kg/m². Its mass is calculated as: $$m={\int }_0^1{\int }_0^2(1+xy)dxdy$$ First, integrate with respect to $x$: ${\int }_0^2(1+xy)dx=[x+\frac{1}{2}{x}^{2}y{]}_0^2=2+2y$. Then, integrate with respect to $y$: ${\int }_0^1(2+2y)dy=[2y+y^2{]}_0^1=3$ kg. This process forms the basis for all subsequent moment calculations.

First Moments and Center of Mass

The center of mass $(\bar{x},\bar{y})$ is the point where the lamina balances perfectly. To find it, we first compute the first moments about the coordinate axes. The moment about the $y$-axis, $M_y$, measures the tendency to rotate about that axis and depends on the distribution of mass in the $x$-direction. Similarly, $M_x$ is the moment about the $x$-axis.

They are given by: $$M_y={\iint }_{R}x\rho (x,y)dA,M_x={\iint }_{R}y\rho (x,y)dA$$

The coordinates of the center of mass are then the total moment divided by the total mass: $$\bar{x}=\frac{M_y}{m},\bar{y}=\frac{M_x}{m}$$

If the density is constant, the center of mass is the centroid of the region, depending only on its geometry. For a region with symmetry and constant density, the centroid lies on the axis of symmetry. Finding the center of mass is crucial in engineering for ensuring stability and predicting rotational dynamics.

Moments of Inertia

While the first moment uses distance to an axis, the moment of inertia (or second moment) uses the square of the distance. This quantity measures an object's resistance to changes in rotational motion about an axis. For rotation about the $x$-, $y$-, and $z$-axes (the latter being perpendicular to the $xy$-plane through the origin), the formulas are:

$$I_x={\iint }_R{y}^2\rho (x,y)dA,I_y={\iint }_R{x}^2\rho (x,y)dA,I_0={\iint }_R(x^2+y^2)\rho (x,y)dA$$

Here, $I_0=I_x+I_y$ is the polar moment of inertia about the origin. Notice the integrand is the square of the distance from the point to the axis: $y^2$ is the squared distance to the $x$-axis. Engineers use these values to design beams, shafts, and rotating machinery, as a higher moment of inertia means greater resistance to angular acceleration for a given torque.

Probability and Average Value Applications

Double integrals are powerful tools in probability theory for continuous random variables. If $X$ and $Y$ are continuous random variables with a joint probability density function $f(x,y)$ over a region $R$, then the probability that $(X,Y)$ lies in a subregion $D\subseteq R$ is: $$P((X,Y)\in D)={\iint }_{D}f(x,y)dA$$ A valid pdf must satisfy ${\iint }_{R}f(x,y)dA=1$, meaning the total probability is 1.

Closely related is the concept of the average value of a function over a region. For a function $f(x,y)$ defined over $R$, its average value is the integral of $f$ over $R$ divided by the area of $R$: $$f_{\text{avg}}=\frac{1}{\text{Area}(R)}{\iint }_{R}f(x,y)dA$$ This is analogous to adding up all possible values and dividing by the "number" of points. In probability, the expected value of a function $g(X,Y)$ is ${\iint }_{R}g(x,y)f(x,y)dA$.

Engineering Synthesis: Analyzing a Lamina

In a complete engineering analysis of a lamina, you systematically apply these concepts. Start by sketching the region $R$ and determining its bounds. Define the density function $\rho (x,y)$, which could model material corrosion (density decreasing from an edge) or a composite material. The workflow is:

1. Compute the mass $m$.
2. Compute the first moments $M_x$ and $M_y$ to find the center of mass $(\bar{x},\bar{y})$.
3. Compute the moments of inertia $I_x$, $I_y$, and $I_0$.

For instance, analyzing a quarter-circular plate of radius $a$ in the first quadrant with constant density involves using polar coordinates for simplicity. The mass is $m=\rho \cdot (\pi a^2/4)$. The moments become integrals like $M_y=\rho {\int }_0^{\pi /2}{\int }_0^a(r\cos \theta )rdrd\theta$. This structured approach allows engineers to predict how a component will behave under force and rotation.

Common Pitfalls

1. Omitting the Density Function in Moments: A frequent error is calculating a centroid (geometric center) when the problem asks for a center of mass with variable density. Remember, for center of mass, the integrand for $M_y$ is $x\rho (x,y)$, not just $x$. If density is constant, it cancels out, but you must include it in the integral initially.
2. Incorrect Order of Integration for Moments of Inertia: Confusing $I_x$ and $I_y$ is easy. Associate the variable in the integrand with the axis: $I_x$ involves $y^2$ (distance squared to the $x$-axis). Writing $I_x=\iint x^2\rho dA$ is a critical mistake that invalidates your physics.
3. Misapplying Probability Rules: In probability applications, ensure your double integral over the entire sample space $R$ equals 1 to confirm $f(x,y)$ is a valid pdf. Also, when finding a probability like $P(X>Y)$, carefully set up the region $D$: it consists of all points $(x,y)$ where $x>y$, within the bounds of $R$.
4. Neglecting Coordinate System Choice: Attempting to integrate over a circular region using rectangular coordinates often leads to extremely complex, if not impossible, integrals. Recognize symmetry and switch to polar coordinates where $dA$ becomes $rdrd\theta$ and $x^2+y^2=r^2$. This simplifies calculations for moments of inertia and mass dramatically.

Summary

• The mass of a lamina with variable density $\rho (x,y)$ over a region $R$ is $m={\iint }_R\rho (x,y)dA$.
• The center of mass $(\bar{x},\bar{y})$ balances the lamina and is found using first moments: $\bar{x}=\frac{1}{m}{\iint }_{R}x\rho dA$, $\bar{y}=\frac{1}{m}{\iint }_{R}y\rho dA$.
• Moments of inertia $I_x$ and $I_y$ quantify rotational resistance and are calculated by integrating the square of the distance to the axis times density.
• In probability, the double integral of a joint pdf over an event region gives the probability of that event.
• The average value of a function over a region is its integral divided by the area of the region.
• A complete engineering analysis synthesizes these calculations, often simplified by choosing an appropriate coordinate system like polar coordinates for circular regions.