Aerodynamic and physical foundations of Wing Spar Unsymmetrical Bending.
Aircraft structural components must achieve maximum structural efficiency with minimum weight, leading to thin-walled open and closed cross-sections (such as Z-stringers, J-sections, and asymmetric angle spars) that rarely possess two planes of structural symmetry. When an asymmetric beam section is subjected to pure vertical bending loads (e.g. aerodynamic lift acting along the y-axis), it does not simply bend in the vertical plane; instead, it undergoes simultaneous out-of-plane lateral bending deflection u along the x-axis. This cross-coupling phenomenon—known as unsymmetrical bending—arises because the geometric product of inertia I_xy = int (x * y) dA is non-zero about the geometric reference axes. According to generalized beam bending theory, the longitudinal bending stress sigma_z at any cross-sectional coordinate (x, y) is given by: sigma_z = [(M_x * I_yy - M_y * I_xy) * y + (M_y * I_xx - M_x * I_xy) * x] / (I_xx * I_yy - I_xy^2). Under pure vertical cantilever loading (M_y = 0, M_x = P_y * (L - z)), the non-zero product of inertia forces lateral deflection u = [I_xy / (E * (I_xx * I_yy - I_xy^2))] * (P_y * L^3 / 3). In this laboratory experiment, a cantilevered aluminum Z-section aircraft spar is loaded vertically with precision weights from P_y = 1.0 to 6.0 kN. Dial indicator gauges measure lateral tip deflection u, validating the linear cross-coupling compliance coefficient.
Viva voce questions & analytical theory
What is unsymmetrical bending and why does it occur in aircraft spars?
Unsymmetrical bending occurs when a beam is loaded in a plane that does not coincide with one of its principal planes of inertia. In asymmetric sections like Z-stringers or L-angles, the product of inertia I_xy is non-zero. A pure vertical load produces both vertical curvature and an out-of-plane lateral curvature, causing the beam to deflect sideways as it bends downward.
What are the principal axes of inertia of a cross section?
The principal axes (u, v) are the orthogonal axes of a cross section about which the product of inertia vanishes identically: I_uv = 0. About these axes, bending is decoupled: a load applied along a principal axis produces deflections purely in that same direction without lateral coupling.
How is the orientation of the principal axes determined mathematically?
The angle theta_p between the geometric reference axes and the principal axes is given by: tan(2 * theta_p) = -2 * I_xy / (I_xx - I_yy). The principal moments of inertia are I_u,v = (I_xx + I_yy)/2 +/- sqrt(((I_xx - I_yy)/2)^2 + I_xy^2).
Where does the neutral axis lie in unsymmetrical bending?
In unsymmetrical bending, the neutral axis (the locus of zero axial stress sigma_z = 0) is not perpendicular to the applied load vector. Its inclination angle alpha_na relative to the horizontal axis is given by tan(alpha_na) = (M_y * I_xx - M_x * I_xy) / (M_x * I_yy - M_y * I_xy).
Why do aircraft structures use Z-section stringers instead of symmetric I-sections?
Z-stringers allow easy riveting of the outer skin to one flange and internal rib attachments to the opposing flange without tool interference. Furthermore, under aerodynamic wing skin compression, Z-sections provide high local buckling resistance and convenient nesting during manufacturing.
What is shear lag in thin-walled aircraft wing structures?
Elementary beam theory assumes plane sections remain plane, implying uniform flange stress. In real wide thin-walled aircraft wing boxes, the shear flexibility of the thin skin delays the transfer of bending stresses from the spar webs to the center of the skin panels, causing peak stresses at the spar edges and reduced stress in the middle (shear lag).
Experimental precautions and wind tunnel safety
- Ensure the applied load acts purely vertically through the shear center to avoid introducing unwanted torsional twisting.
- Verify that the cantilever root clamp is completely rigid with zero rotational compliance under 6.0 kN bending moments.
- Do not exceed the elastic yield limit of the 7075-T6 aluminum alloy (sigma_yield = 500 MPa) during maximum load tests.
- Check that dial indicator contact points have low-friction rolling tips to eliminate lateral binding during deflection.
Theoretical assumptions & flow limits
- Assumes linear elastic Bernoulli-Euler beam behavior without localized flange shear lag.
- Thin-walled open sections can experience localized flange buckling if web-to-thickness ratio exceeds critical limits.
- Root clamp stress concentration induces slight localized compliance not captured in idealized beam formulas.
- Assumes constant cross-sectional geometry along the entire span (untapered spar).