Applied Mathematics

Chapter 6
Integral Calculus

Bachelors Program in Machine Technology
Applied Mathematics   •   Bachelor

Integral Calculus

Introduction

Integral calculus provides a systematic method for determining accumulated quantities. In machine technology, integration is used to calculate displacement from velocity, work from variable force, material mass from density, centroids of components, moments of inertia, fluid volume, and cumulative production. It is also essential when a quantity varies continuously rather than remaining constant.

Differentiation describes instantaneous change; integration reverses this process and describes accumulation. The two ideas are connected by the Fundamental Theorem of Calculus. In practical engineering, an integral may be evaluated exactly from a formula, approximated numerically from measured data, or extended to two and three dimensions for distributed areas, surfaces, and volumes.

Learning Objectives

After studying this chapter, students should be able to:

  • distinguish indefinite, definite, and improper integrals;
  • apply antiderivatives and the Fundamental Theorem of Calculus;
  • evaluate integrals using substitution and integration by parts;
  • estimate integrals using the trapezoidal and Simpson rules;
  • formulate single and multiple integrals for work, energy, mass, centroids, and moments of inertia;
  • interpret units, assumptions, accuracy, and limitations in engineering calculations;
  • use integration to analyze accumulated production and machine-system quantities.

1. Meaning of an Integral

An integral can be understood in two complementary ways. Geometrically, a definite integral represents signed area under a curve. Physically, it represents the accumulation of small contributions. If a quantity changes with position or time, the total is found by adding infinitely many small contributions.

Suppose a force F(x) acts along a straight path. Over a small displacement Δx, the approximate work is F(x)Δx. In the limiting process, the total work from x = a to x = b is

W = ∫ab F(x) dx

Here W is work in joules, F(x) is force in newtons, and x is displacement in metres. The unit check is N·m = J.

The symbol ∫ represents summation of infinitesimal contributions. The differential dx identifies the variable with respect to which accumulation occurs.

Indefinite and definite integrals

An indefinite integral is a family of antiderivatives:

∫ f(x) dx = F(x) + C, where F′(x) = f(x)

The constant C is necessary because all constant shifts have the same derivative. A definite integral has specified limits and produces a number:

∫ab f(x) dx = F(b) − F(a)

provided f is continuous on [a,b], or is otherwise integrable there. The result may be positive, negative, or zero because it is signed area. Geometric area requires attention to the sign of f.

2. Fundamental Integration Rules

For n ≠ −1,

∫ xn dx = x^(n+1)/(n+1) + C

The special logarithmic case is

∫ 1/x dx = ln|x| + C

Other frequently used rules are:

∫ ex dx = ex + C
∫ sin x dx = −cos x + C
∫ cos x dx = sin x + C

For a constant k,

∫ k f(x) dx = k∫ f(x) dx

and for integrable functions f and g,

∫ [f(x) ± g(x)] dx = ∫ f(x) dx ± ∫ g(x) dx

The power rule is not applicable when n = −1; this is why 1/x must be treated separately.

Worked Example 1: Polynomial accumulation

A conveyor speed is described by v(t) = 2t2 + 3t m/s for 0 ≤ t ≤ 4 s. Determine the distance travelled.

Distance is the integral of velocity:

s = ∫04 (2t2 + 3t) dt
s = [2t3/3 + 3t2/2]04
s = 128/3 + 24 = 66.67 m

The result is positive because the velocity is non-negative throughout the interval.

3. Substitution and Integration by Parts

Substitution

Substitution simplifies an integral by replacing a complicated inner expression with a new variable. If u = g(x), then du = g′(x)dx, and

∫ f(g(x))g′(x) dx = ∫ f(u) du

For a definite integral, either transform the limits to u-values or substitute back before applying the original limits.

Worked Example 2: Substitution

Evaluate

∫01 2x(x2 + 1)3 dx

Let u = x² + 1, so du = 2x dx. When x = 0, u = 1; when x = 1, u = 2.

∫12 u3 du = [u4/4]12 = (16 − 1)/4 = 15/4 = 3.75

Integration by parts

Integration by parts follows from the product rule. For differentiable functions u and v,

∫ u dv = uv − ∫ v du

It is useful when the integrand is a product, such as x e^x, x sin x, or ln x. The choice of u should normally make du simpler.

Worked Example 3: Integration by parts

Evaluate ∫01 x ex dx. Choose u = x and dv = e^x dx. Then du = dx and v = e^x.

∫ x ex dx = xex − ∫ ex dx = ex(x − 1) + C

Therefore,

∫01 x ex dx = [ex(x − 1)]01 = 0 − (−1) = 1

4. The Fundamental Theorem of Calculus

The Fundamental Theorem has two related parts.

If f is continuous on an interval and

G(x) = ∫ax f(t) dt

then

G′(x) = f(x)

This states that differentiation recovers the original rate of accumulation. The evaluation form states that if F′(x) = f(x), then

∫ab f(x) dx = F(b) − F(a)

This theorem makes exact evaluation practical: find an antiderivative and evaluate it at the endpoints.

5. Improper Integrals

An integral is improper when an interval is unbounded or the integrand is unbounded at an endpoint or interior point. It is defined through a limit.

For b → ∞,

∫a^∞ f(x) dx = lim_(b→∞) ∫ab f(x) dx

The integral converges if this limit is finite; otherwise it diverges. For an interior singularity c,

∫ab f(x) dx = ∫ac f(x) dx + ∫cb f(x) dx

where each part is evaluated as a separate one-sided limit. It is not valid to cancel divergent parts merely because a symmetric expression appears finite.

A central comparison result is

∫1^∞ 1/xp dx

which converges when p > 1 and diverges when p ≤ 1. For example, ∫1^∞ 1/x2 dx = 1, whereas ∫1^∞ 1/x dx diverges.

6. Numerical Integration

Exact antiderivatives are not always available. Machine measurements may also be given only as a table. Numerical integration estimates the accumulated area from discrete values.

Composite trapezoidal rule

For n equal subintervals of width h = (b − a)/n,

∫ab f(x) dx ≈ h[f(x0)/2 + f(x1) + ⋯ + f(x_(n−1)) + f(xn)/2]

where x_i = a + ih. The method replaces each curve segment by a straight line. If f has a continuous second derivative, the global error is proportional to h².

Composite Simpson’s 1/3 rule

When n is even,

∫ab f(x) dx ≈ h/3 [f(x0) + f(xn) + 4Σ f(xi)_(odd) + 2Σ f(xi)_(even interior)]

The odd-indexed interior points receive coefficient 4 and the even-indexed interior points receive coefficient 2. For sufficiently smooth functions, the global error is proportional to h⁴, generally making Simpson’s rule more accurate than the trapezoidal rule for the same step size.

Worked Example 4: Trapezoidal estimate

A measured force is recorded as follows:

x (m)0123
F (N)10141317

Estimate work from x = 0 to x = 3 m. Here h = 1 m.

W ≈ 1[10/2 + 14 + 13 + 17/2] = 53 J

The estimate is appropriate when the measured values are sufficiently representative and the force varies smoothly between points. Measurement uncertainty and rapid changes may dominate the numerical approximation error.

Best practices for numerical integration

Use consistent units, inspect the data for outliers, choose a smaller step where the curve changes rapidly, and compare two methods when accuracy matters. A result should not be reported with more significant figures than justified by the measurements. Numerical integration is an approximation, not a substitute for understanding the physical model.

7. Applications in Machine Technology

Work and energy

For a variable force parallel to displacement,

W = ∫ab F(x) dx

If the force is conservative, the change in potential energy is related by ΔU = −W_conservative. For a linear spring obeying Hooke’s law F = kx, where k is stiffness in N/m and x is extension in m,

W = ∫0x ku du = 1/2 kx2

This is the energy stored in an ideal spring. The model assumes a linear elastic range and neglects damping and permanent deformation.

Fluid pressure force

Pressure may vary with depth. If a vertical plate has width w(y) and pressure p(y), the resultant force is

F = ∫ p(y)w(y) dy

For a stationary liquid, p = ρg h, where ρ is fluid density, g is gravitational acceleration, and h is depth below the free surface. The pressure unit is pascal, N/m².

Mass of a variable-density component

For a slender rod with linear density λ(x) in kg/m,

m = ∫ab λ(x) dx

For a three-dimensional body with density ρ(x,y,z),

m = ∭V ρ(x,y,z) dV

If density is constant, this reduces to m = ρV.

Centroid and centre of mass

For a plane lamina with area A,

x̄ = 1/A ∫A x dA, ȳ = 1/A ∫A y dA

For a body of mass m,

x̄ = 1/m ∫V xρ dV, ȳ = 1/m ∫V yρ dV, z̄ = 1/m ∫V zρ dV

For uniform density, the centroid and centre of mass coincide. Symmetry should be used before integrating: a symmetric uniform plate has its centroid on the symmetry axis.

For a region between y = f(x) and y = g(x), with f(x) ≥ g(x) on [a,b], vertical strips give

dA = [f(x) − g(x)]dx

and

A = ∫ab [f(x) − g(x)] dx

The first moments are

∫A x dA = ∫ab x[f(x) − g(x)]dx
∫A y dA = ∫ab [f(x)2 − g(x)2]/2 dx

Moments of inertia

The area moment of inertia about the x-axis is

Ix = ∫A y2 dA

and about the y-axis is

Iy = ∫A x2 dA

Their units are m⁴. These quantities measure resistance to bending and are geometric properties, not mass properties.

The mass moment of inertia about an axis is

I = ∫ r2 dm

where r is the perpendicular distance from the mass element dm to the axis. Its units are kg·m² and it measures resistance to angular acceleration.

The parallel-axis theorem for area is

Ix = Ix,c + Ad2

where I_x,c is the centroidal area moment, A is area, and d is the perpendicular distance between parallel axes. For mass moments,

I = Ic + md2

These formulas require parallel axes, with one axis passing through the centroid or centre of mass.

8. Multiple Integration

A double integral accumulates a quantity over an area:

∫∫R f(x,y) dA

A triple integral accumulates it over a volume:

∭V f(x,y,z) dV

For a rectangular region R = [a,b] × [c,d],

∫∫R f(x,y)dA = ∫ab ∫cd f(x,y) dy dx

Under suitable continuity conditions, the order may be reversed. For circular or rotational geometry, polar coordinates may simplify the calculation:

dA = r dr dθ

with x = r cos θ and y = r sin θ. The factor r is the Jacobian and must not be omitted.

Worked Example 5: Centroid of a triangular plate

Consider the uniform triangular region bounded by x = 0, y = 0, and y = 2 − x. Its vertices are (0,0), (2,0), and (0,2). The area is

A = ∫02 (2 − x)dx = 2

The x-coordinate of the centroid is

x̄ = 1/A ∫02 x(2 − x)dx
= 1/2 [x2 − x3/3]02
= 1/2(4 − 8/3) = 2/3 m

Using horizontal or vertical strips, or the known triangle result, the y-coordinate is also 2/3 m. Thus the centroid is (2/3, 2/3), as expected from the average of the three vertex coordinates for a uniform triangle.

9. Accumulated Production Quantities

If production rate varies with time, cumulative production is the integral of the rate. If q(t) is measured in components per hour, then

Q(T) = Q(0) + ∫0T q(t)dt

where Q(0) is initial inventory and Q(T) is inventory or total completed output at time T. If q(t) is a net rate, it may include production minus scrapped units, downtime losses, or rework.

Similarly, if energy consumption rate P(t) is measured in watts,

E = ∫0T P(t)dt

The result is joules when time is in seconds. If power is in kW and time is in hours, the result is kWh.

Illustrative case: Variable-rate machining cell

Background. A machining cell produces parts at a rate that changes during a six-hour shift because of setup, steady operation, tool changes, and inspection. The idealized rate is

q(t) = 20 + 5t − 0.5t2

where t is in hours and q is parts per hour, for 0 ≤ t ≤ 6.

Problem. Estimate total production and determine whether a target of 240 parts is reached, ignoring scrap.

Analysis and solution.

Q = ∫06 (20 + 5t − 0.5t2)dt
Q = [20t + 2.5t2 − t3/6]06
Q = 120 + 90 − 36 = 174 parts

The target of 240 parts is not reached under the stated model. If the actual process includes a startup delay or scrap, the effective total will be lower. If the rate is based on samples rather than a validated model, numerical integration of shift records may be more defensible.

Outcome and learning points. Integration converts a time-dependent rate into a total. The units must be checked: (parts/hour) × hour = parts. Production decisions should also consider capacity, quality, downtime, maintenance, and uncertainty rather than relying on a single idealized rate.

10. Verification, Assumptions, and Limitations

A professional integral calculation should include the following checks:

  1. Model check: Confirm that the integrand represents a rate, density, force, or other elemental quantity.
  2. Limits check: Confirm that the interval or region matches the physical system.
  3. Unit check: The differential contributes units. For example, N integrated over m gives J.
  4. Sign check: Determine whether negative values represent direction, opposing force, or an invalid model.
  5. Magnitude check: Compare the result with a simple upper or lower bound.
  6. Singularity check: Identify discontinuities and improper limits before applying an antiderivative.
  7. Numerical check: Refine the step size or compare methods when using measured data.

Integration can produce a mathematically precise answer to an inaccurate model. Real machines involve friction, wear, backlash, thermal effects, measurement error, and nonlinear material behaviour. A linear spring formula, for example, is unsuitable beyond the elastic range. Numerical methods also inherit errors from sensors, interpolation, rounding, and insufficient sampling.

Computer algebra, spreadsheets, simulation software, and AI-assisted coding can help evaluate integrals and visualize results. However, the engineer remains responsible for selecting the model, checking assumptions, validating units, and detecting plausible-looking but physically impossible answers. Confidential production data should not be transferred to unapproved external systems, and automated results should be traceable and reviewable.

Contemporary Practice and Best Practices

Modern manufacturing systems collect high-frequency data from sensors, programmable logic controllers, and manufacturing execution systems. Integration is used to estimate cumulative energy, tool wear exposure, material usage, and production output. Digital twins may numerically integrate differential models in real time. Adaptive quadrature and finite-element methods extend integration to complex geometries and fields.

Good practice includes defining the physical system before choosing an integration method, using dimensional analysis, preserving raw data, recording numerical tolerances, and reporting uncertainty. When a design decision depends on an integral, an independent calculation or benchmark case should be used for verification.

Chapter Summary

Integral calculus measures accumulation. An indefinite integral gives a family of antiderivatives, while a definite integral gives signed accumulation over specified limits. The Fundamental Theorem of Calculus connects differentiation and integration and allows definite integrals to be evaluated from antiderivatives.

Substitution and integration by parts provide systematic techniques for more complex expressions. Improper integrals are defined through limits and may converge or diverge. When formulas are unavailable or data are measured, the trapezoidal and Simpson rules provide numerical estimates, with accuracy depending on smoothness, step size, and data quality.

In machine technology, integration calculates work, energy, mass, fluid forces, centroids, area moments of inertia, mass moments of inertia, and accumulated production. Double and triple integrals extend the method to areas and volumes. Reliable engineering use requires correct limits, consistent units, physically appropriate assumptions, numerical verification, and responsible interpretation of uncertainty.

Key Terms

Integral: A mathematical operation representing accumulation or an antiderivative.

Indefinite integral: A family of antiderivatives differing by a constant.

Definite integral: The signed accumulation of a function over an interval.

Antiderivative: A function whose derivative is the given integrand.

Integrand: The function being integrated.

Fundamental Theorem of Calculus: The theorem connecting differentiation with definite integration.

Improper integral: An integral involving an infinite interval or an unbounded integrand.

Convergence: Existence of a finite limiting value for an improper integral or approximation process.

Substitution: A change of variable used to simplify integration.

Integration by parts: A technique based on the product rule.

Numerical integration: Approximation of an integral using discrete values.

Trapezoidal rule: A numerical method that replaces curve segments with straight lines.

Simpson’s rule: A numerical method based on quadratic interpolation.

Centroid: The geometric centre of an area or volume.

Centre of mass: The mass-weighted average position of a body.

First moment: An integral used in determining centroids.

Area moment of inertia: A geometric measure used in bending analysis.

Mass moment of inertia: A measure of resistance to angular acceleration.

Jacobian: The coordinate-transformation factor, such as r in polar area integration.

Accumulation rate: A quantity per unit time whose integral gives a total.

Review Questions

Short Answer Questions

  1. Distinguish between an indefinite integral and a definite integral.
  2. State the Fundamental Theorem of Calculus in evaluation form.
  3. Why is ∫1/x dx equal to ln|x| + C rather than a power-rule expression?
  4. State the condition on n for the improper integral ∫1^∞ 1/xn dx to converge.
  5. Give the units of area moment of inertia and mass moment of inertia.

Long Answer Questions

  1. Explain how integration is used to calculate work done by a variable force, including assumptions and unit verification.
  2. Derive the centroid formulas for a region bounded by two functions using vertical strips.
  3. Compare the composite trapezoidal rule with Simpson’s 1/3 rule, including their data requirements and expected accuracy.
  4. Explain the difference between geometric centroid, centre of mass, area moment of inertia, and mass moment of inertia.
  5. Describe a professional procedure for validating an integral-based engineering calculation.

Analytical / Application Questions

  1. A spring has stiffness 800 N/m and is stretched from 0 to 0.06 m. Calculate the stored energy using integration and state the modelling assumption.
  2. A pump delivers water at a rate q(t) = 12 + 2t litres per minute for 0 ≤ t ≤ 10 minutes. Determine the total volume delivered and express it in litres.
  3. A measured force varies linearly from 50 N at x = 0 m to 110 N at x = 2 m. Estimate the work using the trapezoidal rule and explain why the result is exact for this idealized linear force.

Multiple Choice Questions with Answers

Question 1

Which expression represents the work done by a variable force F(x) acting along a path from x = a to x = b?

A. W = ∫ab F(x)dx

B. W = F(b) − F(a)

C. W = ∫ab xF(x)dx

D. W = dF/dx

Correct Answer: A. W = ∫ab F(x)dx

Question 2

Evaluate ∫02 3x2 dx.

A. 4

B. 6

C. 8

D. 12

Correct Answer: C. 8

Question 3

For the composite Simpson’s 1/3 rule, the number of equal subintervals must be:

A. Odd

B. Even

C. Prime

D. Infinite

Correct Answer: B. Even

Question 4

The improper integral ∫1^∞ 1/x2 dx is:

A. Equal to zero

B. Divergent

C. Convergent and equal to 1

D. Undefined because the lower limit is finite

Correct Answer: C. Convergent and equal to 1

Question 5

If q(t) is measured in parts per hour, the units of ∫0T q(t)dt are:

A. Parts per hour

B. Hours per part

C. Parts

D. Parts per hour squared

Correct Answer: C. Parts

Question 6

The polar-coordinate area element is:

A. dA = dr dθ

B. dA = r dr dθ

C. dA = r2 dr dθ

D. dA = 2πr dr

Correct Answer: B. dA = r dr dθ

Question 7

For a uniform triangular lamina, the centroid is located at:

A. The midpoint of its longest side

B. The intersection of its angle bisectors only

C. The average of the coordinates of its three vertices

D. The vertex with the greatest angle

Correct Answer: C. The average of the coordinates of its three vertices

Question 8

Which quantity has units of m⁴?

A. Mass moment of inertia

B. Area moment of inertia

C. Force

D. Work

Correct Answer: B. Area moment of inertia

Question 9

Using the linear spring law F = kx, the work required to stretch an ideal spring from 0 to x is:

A. kx

B. 1/2 kx2

C. k/x

D. 2kx2

Correct Answer: B. 1/2 kx2

Question 10

In the trapezoidal rule, interior data points generally receive:

A. Half the endpoint weight

B. Zero weight

C. The full step-size weight

D. A weight of four times the endpoint weight

Correct Answer: C. The full step-size weight

Applied Mathematics
Chapter 6 of 15