All projects
MAE 5700 Finite Element Analysis, Cornell University

Aircraft Wing Structural Optimization

Shell finite-element modeling and parametric optimization under deformation and stress constraints

Role
Team of three — thickness studies, mathematical model, and optimization
Status
Completed
Aircraft Wing Structural Optimization

Headline results

7,967 kg
Final mass, Design J — 6.6% under the single-parameter optimum
0.375 m
Deformation limit — the active constraint throughout
252 MPa
Stress limit: 378 MPa yield with a 1.5 factor of safety
10
Rib and spar configurations evaluated
ANSYS MechanicalShell FEAParametric OptimizationMesh SensitivityDesign Exploration
01

The problem

Aircraft wing structures balance strength, stiffness, and weight. Excess structural mass reduces aerodynamic efficiency and increases fuel consumption; insufficient stiffness or strength leads to excessive deformation or failure under load. The objective here was to minimize the mass of a simplified thin-walled wing — skin, spars, and ribs, cantilevered and loaded by aerodynamic pressure plus gravity — subject to a maximum deformation of 0.375 m and a maximum von Mises stress of 252 MPa.

Figure 1: Wing geometry — a cantilevered thin-walled box structure of skin, spars, and ribs.
Figure 1: Wing geometry — a cantilevered thin-walled box structure of skin, spars, and ribs.
02

Why shell elements

For thin-walled components the characteristic thickness is much smaller than the in-plane dimensions, so the three-dimensional continuum can be accurately represented by its mid-surface with thickness carried through the constitutive relations rather than resolved geometrically. This drops the degree-of-freedom count sharply while retaining the dominant bending, membrane, and transverse shear behaviour.

The structural response follows from the principle of minimum total potential energy, π = U − W. For a shell, strain energy integrates membrane strains, bending curvatures, and transverse shear strains over the mid-surface, with membrane, bending, and shear stiffness matrices assembled from Young's modulus, Poisson's ratio, and thickness. Bending stiffness scales with t³ while membrane stiffness scales linearly with t — which is precisely why thickness is such a powerful design variable, and why the response surfaces later come out so strongly nonlinear.

03

Baseline and mesh sensitivity

The initial design — 1 rib, 2 spars, 10 mm skin, 10 mm ribs, 10 mm spars, 4,768 kg — failed both constraints, deforming 0.8168 m against the 0.375 m limit and reaching 309 MPa against the 252 MPa limit. That baseline is the reference point for everything that follows: the optimization is not making a working design lighter, it is finding the lightest design that works at all.

Figure 2: Initial design total deformation — 0.8168 m, well over the limit.
Figure 2: Initial design total deformation — 0.8168 m, well over the limit.
Figure 3: Initial design equivalent von Mises stress — 309 MPa against a 252 MPa allowable.
Figure 3: Initial design equivalent von Mises stress — 309 MPa against a 252 MPa allowable.

Before trusting any optimum, we ran a mesh sensitivity study to separate real stress concentrations from discretization singularities. This matters more than it sounds: a stress singularity at a re-entrant corner will keep rising as the mesh refines and never converge, and an optimizer chasing it will keep adding material to a location that is a modeling artifact.

Figure 4: Baseline mesh.
Figure 4: Baseline mesh.
Figure 5: Refined mesh used for the sensitivity comparison.
Figure 5: Refined mesh used for the sensitivity comparison.
04

Parametric optimization

With geometric mass, maximum total deformation, and maximum equivalent stress set as output parameters, the single-parameter study swept skin thickness alone. The candidate point came back at 0.01885 m skin thickness for 8,533 kg — heavier than the baseline, which is exactly right, because the baseline was infeasible. The response charts confirm the expected behaviour: mass rises linearly with skin thickness while deformation and stress fall off sharply.

Figure 6: Maximum deformation against skin thickness — the constraint that binds first.
Figure 6: Maximum deformation against skin thickness — the constraint that binds first.
Figure 7: Maximum equivalent stress against skin thickness.
Figure 7: Maximum equivalent stress against skin thickness.
Figure 8: Geometric mass rises linearly with skin thickness — the cost side of the trade.
Figure 8: Geometric mass rises linearly with skin thickness — the cost side of the trade.

Opening skin, rib, and spar thickness together brought mass down to 8,384.6 kg at 0.375 m deformation and 150 MPa stress — a better result, and the first indication that the design is deformation-controlled rather than stress-controlled.

Figure 9: Three-parameter optimum total deformation — right at the 0.375 m limit.
Figure 9: Three-parameter optimum total deformation — right at the 0.375 m limit.
Figure 10: Three-parameter optimum equivalent stress, 150 MPa — well inside the allowable.
Figure 10: Three-parameter optimum equivalent stress, 150 MPa — well inside the allowable.
05

Design exploration: rib and spar layout

Knowing the model was deformation-controlled, we bounded maximum deformation between 0.36 m and 0.37 m to sit just inside the 0.375 m limit, capped stress at 252 MPa, and minimized mass across ten evenly-spaced rib and spar configurations.

VariationRibsSparsSkin (m)Mass (kg)Max def. (m)Max stress (MPa)
Initial120.01004,7680.8168309
1-parameter120.01888,5330.368148
3-parameter120.01858,3850.375150
C110.019488,5860.3747168
D130.016358,1980.3747161
E210.018738,3290.3282149
F310.018858,4660.3750170
G220.017788,1210.3749169
H320.018458,0910.3747192
I230.018238,0060.3749196
J (recommended)330.018047,9670.37497189
Figure 11: Design J total deformation — 0.37497 m, just inside the limit.
Figure 11: Design J total deformation — 0.37497 m, just inside the limit.
Figure 12: Design J equivalent stress — 189 MPa against a 252 MPa allowable.
Figure 12: Design J equivalent stress — 189 MPa against a 252 MPa allowable.
06

Conclusion

Design J — three ribs, three spars, 18.04 mm skin, 3.094 mm ribs, 3.897 mm spars — reaches 7,967 kg while staying inside both the 0.375 m deformation and 252 MPa stress limits. Across every study the same trend held: more ribs and spars let us reduce total deformation while thinning every member. Because the design is deformation-controlled, the optimum is the configuration with the most internal structure.