Abstract In this paper, a novel method is proposed to support the process of solving multi-objective nonlinear programming problems subject to strict or flexible constraints. This method assumes that the practical problems are expressed in the form of geometric programming problems. Integrating the concept of intuitionistic fuzzy sets into the solving procedure, a rich structure is provided which can include the inevitable uncertainties into the model regarding different objectives and constraints. Another important feature of the proposed method is that it continuously interacts with the decision maker. Thus, the decision maker could learn about the problem, thereby a compromise solution satisfying his/hers preferences could be obtained. Further, a new two-step geometric programming approach is introduced to determine Pareto-optimal compromise solutions for the problems defined during different iterative steps. Employing the compensatory operator of “weighted geometric mean”, the first step concentrates on finding an intuitionistic fuzzy efficient compromise solution. In the cases where one or more intuitionistic fuzzy objectives are fully achieved, a second geometric programming model is developed to improve the resulting compromise solution. Otherwise, it is concluded that the resulting solution vectors simultaneously satisfy both of the conditions of intuitionistic fuzzy efficiency and Pareto-optimality. The models forming the proposed solving method are developed in a way such that, the posynomiality of the defined problem is not affected. This property is of great importance when solving nonlinear programming problems. A numerical example of multi-objective nonlinear programming problem is also used to provide a better understanding of the proposed solving method.