In aerospace engineering, design and optimization of mechanical structures are usually performed with respect to elastic limit. Besides causing insufficient use of the material, such design concept fails to meet the ever growing needs of the light weight design. To remedy this problem, in the present study, a shakedown theory based numerical approach for performing parametric optimization is presented. Within this approach, strength of the structure is measured by its shakedown limit calculated from the direct method. The numerical method developed for the structural optimization consists of nested loops: The inner loop adopts the interior point method to solve shakedown problems pertained to fixed design parameters, while the outer loop employs the genetic algorithm to find optimal design parameters leading to the greatest shakedown limit. The method established is first verified by the classic plate-with-a-circular-hole example, and after that it is applied to an airtight module for determining few key design parameters. By carefully analyzing results generated during the optimization process, it is convinced that the approach can become a viable means for designing similar aerospace structures.