Wing Torsional Divergence lab report

1. Title and aim

Aeroelastic Torsional Divergence and Wing Torsional Rigidity. Measure elastic wing twist theta under aerodynamic loading in a wind tunnel and extrapolate the critical divergence dynamic pressure q_div and airspeed V_div.

2. Abstract

Six controlled aerodynamic test points were executed across standardized inputs q = 0.3 kPa to q = 1.8 kPa for Aeroelastic Torsional Divergence and Wing Torsional Rigidity. An ordinary least-squares linear regression yielded an empirical Inverse Divergence Pressure 1/q_div of 0.3125 kPa^-1 (difference from theoretical reference = 0.0046%). The empirical telemetry conclusively validates the fundamental aerospace engineering formulation with R² = 1.000000. The investigation demonstrates stable flow behavior and mechanical response across all test conditions with negligible residual variance. Data recorded through the virtual wind tunnel and telemetry bus exhibit strict adherence to established boundary layer, shock wave, flight mechanics, or structural mechanics governing equations, confirming the dimensional and similitude assertions of the curriculum.

3. Apparatus and fixed settings

• Aeroelastic Wind Tunnel Test Wing: NACA 0012 rigid aerodynamic foam shell mounted over a flexible internal steel torsion rod (elastic axis at 45% chord).
• Variable Speed Aeroelastic Wind Tunnel: Delivers controllable dynamic pressures q from 0.1 to 2.5 kPa with safety catch nets.
• Laser Tip Optical Inclinometer System: Reflects laser beam from tip wing mirror onto optical sensor to measure twist angle theta with 0.01 degree accuracy.
• Aerodynamic Center Offset Mount: Adjusts aerodynamic center-to-elastic axis offset e from 15 mm to 35 mm.
• Dynamic Pressure Pitot-Static Sensor: Continuously records test-section dynamic pressure q in kPa.
• Aeroelastic Safety Interlock System: Automatically trips wind tunnel power if wing twist angle exceeds 8.0 degrees to prevent model destruction.
• Airfoil chord c: 150 mm
• Semi-span l: 600 mm
• Aerodynamic offset e: 30 mm
• Initial incidence alpha_0: 2.0 deg
• Target q_div: 3.20 kPa

4. Theory and graph choice

Aeroelasticity investigates the dangerous mutual interaction between aerodynamic forces, elastic structural flexibility, and inertial dynamics on flight vehicles. When an aerodynamic lifting surface deflects under air loads, the structural deformation alters the aerodynamic pressure distribution, creating a closed-loop feedback mechanism. One of the most catastrophic static aeroelastic instabilities is wing torsional divergence. On an unswept wing, the aerodynamic center (where lift acts, approx quarter-chord c/4) typically lies ahead of the structural elastic axis (shear center axis, approx mid-chord c/2) by an offset distance e. When the wing experiences an initial angle of attack alpha_0, aerodynamic lift generates a nose-up twisting torque about the elastic axis. This twists the wing elastically by an angle theta, increasing the local angle of attack to (alpha_0 + theta). This additional incidence generates extra lift, creating even more torque. The structural torsional restoring stiffness resists this twisting: M_elastic = G * J * (dtheta/dz). However, aerodynamic moment scales with dynamic pressure q = 0.5 * rho * V^2, whereas structural stiffness remains constant. At the critical divergence dynamic pressure q_div, aerodynamic twisting moment completely overpowers structural torsional rigidity, causing the twist amplification factor theta / alpha_0 = 1 / (1 - q / q_div) to diverge to infinity and snapping the wing off in flight. Inverting this relation yields the linear form (1 - alpha_0 / theta) = (1 / q_div) * q. In this experiment, a wind tunnel elastic wing model is tested across dynamic pressures from q = 0.30 to 1.80 kPa, verifying q_div = 3.20 kPa.

The governing physical relationship was linearized into the standard regression form y = mx + b. Here, the independent coordinate x represents Wind Tunnel Dynamic Pressure q (kPa), and the dependent coordinate y represents Twist Amplification Metric 1 - alpha_0 / theta. The slope m directly extracts Inverse Divergence Pressure 1/q_div, while the intercept b accounts for zero-lift offsets, tare baselines, or residual reference shifts. Linear transformation allows robust parameter estimation via ordinary least squares without solver instability or local minima traps. Plotting transformed observation pairs highlights flow similitude coefficients, lift curve slopes, pressure ratios, or stiffness parameters directly as linear gradients. This linear graphical formulation provides an immediate visual check of homoscedasticity, linearity, and potential flow separation or structural non-linearities.

5. Method

1. Set the root incidence angle to an initial geometric angle of attack alpha_0 = 2.0 degrees with zero wind.
2. Zero the optical laser mirror receiver, confirming theta = 0 deg relative to root incidence.
3. Run the wind tunnel to dynamic pressure q = 0.30 kPa, recording steady-state aeroelastic twist angle theta.
4. Step dynamic pressure through q = 0.60, 0.90, 1.20, 1.50, and 1.80 kPa, recording twist angle theta at each plateau.
5. Compute total incidence (alpha_0 + theta) and calculate the non-dimensional metric (1 - alpha_0 / theta).
6. Plot (1 - alpha_0 / theta) versus q, perform linear regression to determine slope 1 / q_div, and calculate divergence dynamic pressure q_div = 1 / m.

6. Observations and calculation table

All aerodynamic observations use strictly controlled 64-bit IEEE 754 precision buffers. Residuals are computed relative to the ordinary least-squares regression line.

Download the sample observations (CSV) ↓

7. Worked calculation for one observation

8. Graph, gradient and final calculation

Circular markers represent empirical aerodynamic observations. The solid line represents the ordinary least-squares linear fit.

9. Uncertainty and comparison

Percentage difference = |estimate − theoretical value| / |theoretical value| × 100 = 0.0046%.

For n = 6 observations, Σresidual² = 3.429e-09; s² = Σresidual²/(n − 2). Gradient standard error = √[s²/Σ(x − x̄)²] = 2.333e-05. Intercept standard error = √[s²(1/n + x̄²/Σ(x − x̄)²)] = 2.726e-05.

Propagating the fit parameters yields a fit standard error of 2.333e-05 kPa^-1, with an ideal coefficient of determination R² = 1.000000.

In experimental aeronautics and flight vehicle testing, empirical uncertainties originate from dynamic pressure transducer calibration drift, wind tunnel boundary layer blockage corrections, tare balance interaction tare loads, free-stream turbulence intensity, and numerical discretization (64-bit IEEE 754 precision). The ordinary least-squares regression across n = 6 standardized test points yields unbiased standard errors for both gradient and intercept. Because the governing aerodynamic and structural equations adhere strictly to conservation of mass, momentum, and energy, residual deviations are tightly bounded, yielding a coefficient of determination R² ≥ 0.999999 and confirming the theoretical scaling model. Furthermore, systematic errors—including test section solid blockage, wake blockage, buoyancy pressure gradients along the test section axis, and pitch balance sting flexure deflections—are corrected in accordance with standard aerodynamic testing protocols. Standard deviations are propagated using first-order Taylor series expansion (delta method) across transformed dimensionless parameters. Residual diagnostics demonstrate homoscedasticity across the entire operational envelope without systematic bias.

10. Result and discussion

The empirical wind tunnel and flight test observations provide rigorous quantitative verification of Aeroelastic Torsional Divergence and Wing Torsional Rigidity. In aeronautical engineering, gas dynamics, and aerospace structures, physical systems are governed by Navier-Stokes conservation laws, isentropic compressible relations, rigid-body flight dynamics equations, and continuum elasticity. Across all six observations, the empirical slope m = 0.3125 confirms the theoretical model with near-zero residual variance. Furthermore, comparison with classical aeronautical benchmarks underscores critical design principles including aerodynamic efficiency, static stability margins, control authority, thermodynamic cycle limits, and structural mass optimization. Residual analysis reveals homoscedastic dispersion without systematic non-linear patterning, validating the underlying model assumptions. Confidence intervals calculated via asymptotic standard errors align with theoretical quantiles, indicating that finite-sample distortions remain negligible under controlled test conditions. The experimental protocol confirms the robustness of the methodology across varying flight conditions.

11. Precautions and limitations

• Never operate the wind tunnel near the extrapolated divergence dynamic pressure q_div; test points must remain below 60% of q_div.
• Verify that safety optical interlocks immediately kill fan drive power if twist angle exceeds 8.0 degrees.
• Ensure the internal torsion rod is operating well within its elastic fatigue limit without plastic deformation.
• Account for wind tunnel ceiling and floor aerodynamic interference when testing high-aspect-ratio aeroelastic models.

Linear divergence theory assumes 2D strip aerodynamics without 3D wingtip vortex relief. Linear regression of the twist metric against dynamic pressure yields a slope m = 0.3125 kPa^-1 corresponding to q_div = 3.200 kPa with R^2 > 0.999999.

12. Conclusion

This laboratory investigation quantitatively verified the mathematical formulation and physical behavior of Aeroelastic Torsional Divergence and Wing Torsional Rigidity. Through rigorous tracking across six observation points, we confirmed an empirical Inverse Divergence Pressure 1/q_div of 0.3125 kPa^-1 matching the theoretical reference 0.3125 kPa^-1 within 0.0046% error. The investigation demonstrates why this aerospace method serves as a fundamental cornerstone in aircraft and spacecraft design. The verified scaling laws and parameter estimates demonstrate the practical efficacy of aerodynamic modeling and aerospace engineering principles.