Two closely-related fictitious domain methods for solving problems involving multiple interfaces are introduced. Like other fictitious domain methods, the proposed methods simplify the task of finite element mesh generation and provide access to solvers that can take advantage of uniform structured grids. The proposed methods do not involve the Lagrange multipliers, which makes them quite different from existing fictitious domain methods. This difference leads to an advantageous form of the inf–sup condition, and allows one to avoid time-consuming integration over curvilinear surfaces. In principle, the proposed methods have the same rate of convergence as existing fictitious domain methods. Nevertheless it is shown that, at the cost of introducing additional unknowns, one can improve the quality of the solution near the interfaces. The methods are presented using a two-dimensional model problem formulated in the context of linearized theory of elasticity. The model problem is sufficient for presenting method details and mathematical foundations. Although the model problem is formulated in two dimensions and involves only one interface, there are no apparent conceptual difficulties to extending the methods to three dimensions and multiple interfaces. Further, it is possible to extend the methods to nonlinear problems involving multiple interfaces.