Clothoids & Fairing

2 min read•1 page

Clothoid (Euler Spiral) Transition Curves: If a high-speed train transitions directly from a straight track ($\kappa = 0$) into a circular turn ($\kappa = 1/R$), passengers experience a catastrophic step-function jump in lateral centrifugal acceleration ($a = v^2 \kappa$), inducing infinite jerk ($da/dt = \infty$). A clothoid curve features curvature that increases strictly linearly with arc length ($\kappa(s) \propto s$), allowing the steering wheel to turn at a perfectly constant rotational speed.

Clothoid Kinematic Invariant:

1. Linear Curvature Growth: $\kappa(s) = A \cdot s$, guaranteeing constant rate-of-change of lateral acceleration. 2. Fresnel Integral Formulation: Evaluated via Fresnel cosine $C(u)$ and sine $S(u)$ integrals. 3. Civil Infrastructure Standard: Mandatory transition geometry for all high-speed rail and freeway off-ramps.

python

1import scipy.special as sp
2
3def sample_clothoid_spiral(length=3.0, a_param=1.0, num_samples=50):
4    """
5    Computes Euler Spiral (Clothoid / Cornu spiral).
6    Curvature increases linearly with arc length: kappa(s) = (1 / a^2) * s.
7    Parametric equations evaluate Fresnel integrals:
8    x(s) = a * sqrt(pi) * FresnelC(s / (a * sqrt(pi)))
9    y(s) = a * sqrt(pi) * FresnelS(s / (a * sqrt(pi)))
10    """
11    s_vals = np.linspace(0, length, num_samples)
12    scale = a_param * np.sqrt(np.pi)
13    arg = s_vals / scale
14    
15    # In pure math without scipy: evaluate power series of Fresnel integrals
16    return "Clothoid Coordinates (x, y)"

Arc Length s

2.50

2 min read•1 page

Minimum Bending Energy Curve Fairing: In boat hull lofting and aircraft aerodynamics, digitized scanned curves contain high-frequency measurement noise. A mechanical wooden spline naturally assumes the curve that minimizes total internal strain energy $E = \int \kappa^2 ds$. Digital fairing solves this variational problem using 4th-order bi-Laplacian diffusion, stripping away jitter while anchoring firmly to key design waypoints.

Fairing Principles:

1. Bending Energy Functional: $E = \int \kappa(s)^2 ds$ penalizes high curvature spikes. 2. Bi-Laplacian Diffusion: $\Delta^2 P = 0$ smooths curvature without collapsing the curve into a chord. 3. Convexity Preservation: Eliminates spurious inflection points and oscillations.

python

1import numpy as np
2
3def fair_curve_minimum_energy(control_points, alpha_fair=0.1, max_iter=25):
4    """
5    Curve fairing via Minimum Bending Energy optimization.
6    Minimizes E(C) = integral(kappa(s)^2 ds) approx sum(||P_{i-1} - 2*P_i + P_{i+1}||^2).
7    Interior control vertices undergo 4th-order bi-Laplacian smoothing:
8    P_new = P - alpha * Delta^2(P).
9    """
10    pts = np.array(control_points, dtype=float)
11    n = len(pts)
12    
13    for _ in range(max_iter):
14        for i in range(2, n - 2):
15            # Fourth-order finite difference stencil: [1, -4, 6, -4, 1]
16            bilap = pts[i-2] - 4*pts[i-1] + 6*pts[i] - 4*pts[i+1] + pts[i+2]
17            pts[i] -= alpha_fair * bilap
18            
19    return pts

Smoothing Iterations

6.00