Developments in rapid manufacturing techniques such as electron beam freeform fabrication have made it easier to manufacture objects with complex shapes (like panel with curvilinear stiffeners). In this paper, design optimization of stiffened panels is carried out using different types of response surfaces. Placement and sizing variables are optimized using integrated high-performance, commercially available software. The proposed approach has two optimization loops. The sizing optimization is carried out in the first loop through finite element analysis using gradient-based optimization method for different placement variables. In the second loop, optimal placement of stiffeners is carried out using response surfaces. Because of response surface construction using whole design space, the proposed approach does not suffer from getting trapped in a local minima. The proposed approach is applied to square panels subjected to in-plane biaxial and uniform surface pressure loadings. Panel weight is minimized with buckling, displacements, and stress constraints. An approach to perform both placement and sizing optimization of panels with curvilinear stiffeners is also discussed and a surrogate model is developed for placement and shape optimization of stiffeners.Nomenclature a = length of the stiffened plate b = width of the stiffened plate c = stiffener shape design variable d = direction vector perpendicular to stiffener reference curve h = height of the stiffener t = thickness of the plate w = width of the stiffener w i = weights of the control point P i of nonuniform rational basis-splines curve w i;j = weights of the control point P i;j of nonuniform rational basis-splines surface y 1 = right end y coordinate of the first stiffener y 2 = left end y coordinate of the first stiffener y 3 = right end y coordinate of the second stiffener y 4 = left end y coordinate of the second stiffener Ct = nonuniform rational basis-splines parametric curve L = length of the stiffener N i;p = B-spline basis function N st = number of the stiffeners P i = control point's coordinates of nonuniform rational basis-splines curve P i;j = control point's coordinates of nonuniform rational basis-splines surface Su; v = nonuniform rational basis-splines parametric surface W = mass of the stiffened panel expressed in terms of y 1 , y 2 , y 3 and y 4 , kg W 1 = mass of the stiffened panel in terms of y 1 , kg W 2 = mass of the stiffened panel in terms of y 2 , kg W 3 = mass of the stiffened panel in terms of y 3 , kg W 4 = mass of the stiffened panel in terms of y 4 , kg W = mass of the stiffened panel expressed in terms of y 1 , y 2 , y 3 , y 4 and , kg = shape design parameter for the curvilinear stiffeners = buckling eigenvalue = density of the material of the stiffened plate = von Mises stress in the stiffened plate