In this work we present an approach that can be systematically used to construct nonlinear systems possessing analytical multikink profile configurations. In contrast with previous approaches to the problem, we are able to do it by using field potentials that are considerably smoother than the ones of the doubly quadratic family of potentials. This is done without losing the capacity of writing exact analytical solutions. The resulting field configurations can be applied to the study of problems from condensed matter to braneworld scenarios.Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)