The main research goal is the development of a numerical methodology for solving parametric optimization problems of steel structures with orientation to software implementation in a computer-aided design system. The paper introduces a new mathematical model for parametric optimization problems of steel structures. The design variable vector includes geometrical parameters of the structure (node coordinates), cross-sectional dimensions of the structural members, as well as initial pre-stressing forces introduced into the specified redundant members of the structure. The system of constraints covers load-carrying capacities constraints formulated for all design sections of structural members of the steel structure subjected to all ultimate load case combinations. The displacements constraints formulated for the specified nodes of the steel structure subjected to all serviceability load case combinations have been also included into the system of constraints. The method of the objective function gradient projection onto the active constraints surface with simultaneous correction of the constraints violations has been used for solving the parametric optimization problem. A numerical algorithm for solving the formulated parametric optimization problems of steel structures has been developed in the paper. The comparison of the optimization results of truss structures presented by the paper confirms the validity of the optimum solutions obtained using the proposed numerical methodology.