The steady-state behavior of an existing plant depends on the independent input variables, process equipment and process controllers. This paper presents a method for formulating models that represent the effects of controllers when they are included within a steady-state process flowsheet. The method replaces the controller equations with the equivalent stationarity conditions representing the relationship between the controlled variables and the implemented manipulated variables at steady state. The method is demonstrated for the centralized multivariable Dynamic Matrix Control algorithm applied to two processes, binary distillation and gasoline blending. The integrated process and control system simulation is used to design controllers that improve the profitability of processes without extensive real-time calculations; this is sometimes termed self-optimizing control. For both processes, controllers were designed that yielded higher profit than standard control methods and that approached the highest possible profit achieved by frequent real-time optimization.