Integral CalculusIntroductionIntegral 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 ObjectivesAfter 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 IntegralAn 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 integralsAn 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 RulesFor 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 accumulationA 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 PartsSubstitutionSubstitution 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: SubstitutionEvaluate ∫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 partsIntegration 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 partsEvaluate ∫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 CalculusThe 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 IntegralsAn 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 IntegrationExact 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 ruleFor 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 ruleWhen 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 estimateA measured force is recorded as follows: 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 integrationUse 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 TechnologyWork and energyFor 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 forcePressure 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 componentFor 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 massFor 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 inertiaThe 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 IntegrationA 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 plateConsider 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 QuantitiesIf 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 cellBackground. 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 LimitationsA professional integral calculation should include the following checks: - Model check: Confirm that the integrand represents a rate, density, force, or other elemental quantity.
- Limits check: Confirm that the interval or region matches the physical system.
- Unit check: The differential contributes units. For example, N integrated over m gives J.
- Sign check: Determine whether negative values represent direction, opposing force, or an invalid model.
- Magnitude check: Compare the result with a simple upper or lower bound.
- Singularity check: Identify discontinuities and improper limits before applying an antiderivative.
- 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 PracticesModern 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 SummaryIntegral 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 TermsIntegral: 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 QuestionsShort Answer Questions- Distinguish between an indefinite integral and a definite integral.
- State the Fundamental Theorem of Calculus in evaluation form.
- Why is ∫1/x dx equal to ln|x| + C rather than a power-rule expression?
- State the condition on n for the improper integral ∫1^∞ 1/xn dx to converge.
- Give the units of area moment of inertia and mass moment of inertia.
Long Answer Questions- Explain how integration is used to calculate work done by a variable force, including assumptions and unit verification.
- Derive the centroid formulas for a region bounded by two functions using vertical strips.
- Compare the composite trapezoidal rule with Simpson’s 1/3 rule, including their data requirements and expected accuracy.
- Explain the difference between geometric centroid, centre of mass, area moment of inertia, and mass moment of inertia.
- Describe a professional procedure for validating an integral-based engineering calculation.
Analytical / Application Questions- 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.
- 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.
- 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 AnswersQuestion 1Which 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 2Evaluate ∫02 3x2 dx. A. 4 B. 6 C. 8 D. 12 Correct Answer: C. 8 Question 3For 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 4The 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 5If 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 6The 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 7For 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 8Which 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 9Using 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 10In 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
|