3.3 Differential and Integral Calculus Applied to Curves, Surfaces, and Surveying Computations

3.3 Differential and Integral Calculus Applied to Curves, Surfaces, and Surveying Computations

Calculus provides the essential mathematical framework for analyzing continuous spatial variation in geodetic engineering. From evaluating rate-of-change sensitivity in error propagation to integrating area under terrain profiles and computing earthwork volumes, differential and integral calculus are indispensable tools for professional practice.

---

Differential Calculus and Error Propagation Analysis

Differential calculus analyzes rates of change, tangent slopes of curves, and sensitivity of computed quantities to measurement errors.

Derivatives of Key Functions in Surveying

• Trigonometric Functions: $$\frac{d}{dx}(\sin x) = \cos x, \quad \frac{d}{dx}(\cos x) = -\sin x, \quad \frac{d}{dx}(\tan x) = \sec^2 x$$
• Transcendental Functions: $$\frac{d}{dx}(e^x) = e^x, \quad \frac{d}{dx}(\ln x) = \frac{1}{x}$$

General Law of Variance / Error Propagation

When a calculated quantity $Z$ is a function of independent measured variables $X_1, X_2, \dots, X_n$ with standard errors $\sigma_{X_1}, \sigma_{X_2}, \dots, \sigma_{X_n}$, the total variance $\sigma_Z^2$ is determined using partial derivatives:

$$\sigma_Z^2 = \left(\frac{\partial Z}{\partial X_1}\right)^2 \sigma_{X_1}^2 + \left(\frac{\partial Z}{\partial X_2}\right)^2 \sigma_{X_2}^2 + \dots + \left(\frac{\partial Z}{\partial X_n}\right)^2 \sigma_{X_n}^2$$

Worked Example: Error Propagation in Rectangular Area

Problem: The area of a rectangular survey parcel is $A = W \times L$. Measured dimensions are $L = 100.00 \text{ m} \pm 0.03 \text{ m}$ and $W = 50.00 \text{ m} \pm 0.02 \text{ m}$. Compute the standard error $\sigma_A$ in the calculated area.

Solution:

1. Compute partial derivatives: $$\frac{\partial A}{\partial L} = W = 50.00, \quad \frac{\partial A}{\partial W} = L = 100.00$$
2. Apply the variance propagation formula: $$\sigma_A^2 = \left(\frac{\partial A}{\partial L}\right)^2 \sigma_L^2 + \left(\frac{\partial A}{\partial W}\right)^2 \sigma_W^2 = (W \sigma_L)^2 + (L \sigma_W)^2$$
3. Substitute values: $$W \sigma_L = 50.00 \times 0.03 = 1.50 \text{ m}^2$$ $$L \sigma_W = 100.00 \times 0.02 = 2.00 \text{ m}^2$$ $$\sigma_A^2 = (1.50)^2 + (2.00)^2 = 2.25 + 4.00 = 6.25$$
4. Take square root to get standard error: $$\sigma_A = \sqrt{6.25} = 2.50 \text{ m}^2$$

---

Maxima and Minima Applications in Parabolic Vertical Curve Design

Vertical curves in highway alignment transition between grade tangents using vertical axis parabolas: $y = a x^2 + b x + c$.

Parabolic Curve Equations

If $g_1$ is initial tangent grade (%), $g_2$ is final tangent grade (%), and $L$ is total horizontal length of curve (meters):

$$y(x) = y_{\text{PVC}} + g_1 x + \left(\frac{g_2 - g_1}{2 L}\right) x^2$$

where $x$ is horizontal distance measured from the Point of Vertical Curvature (PVC).

Finding High or Low Turning Point

To locate the maximum elevation on a crest curve or minimum elevation on a sag curve, set the first derivative (slope) to zero:

$$\frac{dy}{dx} = g_1 + \left(\frac{g_2 - g_1}{L}\right) x = 0$$ $$x_{\text{turn}} = \frac{-g_1 L}{g_2 - g_1} = \frac{g_1 L}{g_1 - g_2}$$

Worked Example: High Point of a Crest Vertical Curve

Problem: A vertical crest curve has an initial grade $g_1 = +3.0\%$, final grade $g_2 = -2.0\%$, and length $L = 200 \text{ meters}$. Find the distance $x$ from PVC to the highest point on the curve.

Solution:

1. Identify grades: $g_1 = +0.03$, $g_2 = -0.02$, length $L = 200 \text{ m}$.
2. Substitute into the turning point formula: $$x = \frac{g_1 L}{g_1 - g_2} = \frac{3.0 \times 200}{3.0 - (-2.0)} = \frac{600}{5.0} = 120 \text{ meters}$$

---

Arc Length Integration of Curves and Ellipsoidal Meridian Arcs

Arc length integration calculates exact spatial distances along curved paths.

Standard Arc Length Integral

The arc length $S$ of a curve $y = f(x)$ between $x = a$ and $x = b$ is given by:

$$S = \int_{a}^{b} \sqrt{1 + \left(\frac{dy}{dx}\right)^2} dx$$

Meridian Arc Length on a Reference Ellipsoid

The distance along a meridian on an ellipsoid from equator to latitude $\phi$ requires integrating the meridian radius of curvature $M(\phi)$:

$$M(\phi) = \frac{a(1 - e^2)}{(1 - e^2 \sin^2 \phi)^{3/2}}$$ $$S(\phi) = \int_{0}^{\phi} M(\phi) d\phi = a(1 - e^2) \int_{0}^{\phi} (1 - e^2 \sin^2 \phi)^{-3/2} d\phi$$

Using binomial expansion, this integral is evaluated as a series:

$$S(\phi) = A_0 \phi - A_2 \sin(2\phi) + A_4 \sin(4\phi) - A_6 \sin(6\phi) + \dots$$

---

Area Integration under Curves and Numerical Integration

Determining cross-sectional area under irregular terrain boundaries uses numerical integration techniques.

Numerical Integration Profiles:
Trapezoidal Rule: Linear segments between offset ordinates
Simpson's 1/3 Rule: Parabolic arcs through triplets of ordinates

Trapezoidal Rule

Assumes linear variation between offset ordinates spaced at uniform interval $d$:

$$\text{Area} = \frac{d}{2} \left[ y_0 + y_n + 2 \sum_{i=1}^{n-1} y_i \right]$$

Simpson's 1/3 Rule

Assumes parabolic arcs through sets of three consecutive ordinates (requires an even number of intervals $n$):

$$\text{Area} = \frac{d}{3} \left[ y_0 + y_n + 4 \sum_{i=\text{odd}} y_i + 2 \sum_{i=\text{even}} y_i \right]$$

Comparison of Area Integration Methods

---

Volume of Earthwork: End-Area and Prismoidal Formulas

Volume calculations determine cut and fill quantities for highways, railways, and site grading.

Average End-Area Formula

$$\text{Volume}_{\text{ea}} = \left(\frac{A_1 + A_2}{2}\right) L$$

where $A_1$ and $A_2$ are cross-sectional areas at two stations separated by distance $L$.

Prismoidal Formula

Provides exact volume for regular prismoids:

$$\text{Volume}_p = \frac{L}{6} (A_1 + 4 A_m + A_2)$$

where $A_m$ is the cross-sectional area at mid-length ($L/2$), calculated from averaged mid-dimensions.

Prismoidal Correction ($C_p$)

$$\text{Volume}_p = \text{Volume}_{\text{ea}} - C_p$$ $$C_p = \frac{L}{12} (c_1 - c_2)(w_1 - w_2)$$

where $c_1, c_2$ are center cuts/fills and $w_1, w_2$ are total cross-section widths.

Worked Example: Prismoidal Volume Calculation

Problem: Two highway cross sections separated by $L = 100 \text{ m}$ have end areas $A_1 = 450 \text{ m}^2$ and $A_2 = 650 \text{ m}^2$. The mid-section area computed from average dimensions is $A_m = 540 \text{ m}^2$. Compute the prismoidal volume $V_p$.

Solution:

1. Apply the prismoidal formula: $$V_p = \frac{L}{6} (A_1 + 4 A_m + A_2)$$
2. Substitute values: $$V_p = \frac{100}{6} [450 + 4(540) + 650]$$ $$V_p = \frac{100}{6} [450 + 2160 + 650] = \frac{100}{6} [3260] = \frac{326000}{6} \approx 54,333.33 \text{ m}^3$$