This paper presents a fuzzy set method to enforce complementarity constraints (CCs) in a nonlinear interior point method (NIPM)-based optimization. NIPM is a Newton-type approach to nonlinear programming problems, but it adopts log-barrier functions to deal with the obstacle of managing inequality constraints. The fuzzy-enforcement method has been implemented for CCs, which can be incorporated in optimization problems for real-world applications. In this paper, numerical simulations that apply this method to power system optimal power flow problems are included.