Nonlinear multi-point constraints are essential in modeling various engineering problems, for example in the context of (a) linking individual degrees of freedom of multiple nodes to model nonlinear joints, (b) coupling different element types in finite element analysis, (c) enforcing various types of rigidity in parts of the mesh and (d) considering deformation-dependent Dirichlet boundary conditions. One method for addressing constraints is the master–slave elimination, which offers the benefit of reducing the problem dimension as opposed to Lagrange multipliers and the penalty method. However, the existing master–slave elimination method is limited to linear constraints. In this paper, we introduce a new master–slave elimination method for handling arbitrary smooth nonlinear multi-point constraints in the system of equations of the discretized system. We present a rigorous mathematical derivation of the method. Within this method, new constraints can be easily considered as an item of a “constraint library”; i.e. no case-by-case-programming is required. In addition to the theoretical aspects, we also provide helpful remarks on the efficient implementation. Among others, we show that the new method results in a reduced computational complexity compared to the existing methods. The study also places emphasis on comparing the new approach with existing methods via numerical examples. We have developed innovative benchmarks which encompass all relevant computational properties, and provide analytical and reference solutions. Our findings demonstrate that our new method is as accurate, robust and flexible as the Lagrange multipliers, and more efficient due to the reduction of the total number of degrees of freedom, which is particularly advantageous when a large number of constraints have to be considered.