Despite of being considered a successful theory in the field of fundamental particles and interactions, the Standard Model is considered incomplete because it still leaves a number of open questions. In order to address some of these, it is necessary to study physical processes such as particle decays. In particular, the decays of heavy hadrons such as D and B mesons can be sensitive to physics beyond the Standard Model, such as new sources of charge-parity violation and new interactions. Nevertheless, at masses at the few GeV range, where perturbative calculations of quantum chromodynamics are limited, phenomenological models are needed to describe their hadronic decays. In this dissertation, a study of the 3-body decay dynamics of the Ds+ to K-K+K+ channel is presented for the first time and uses two phenomenological approaches to describe the decay amplitude: the Isobar Model and the Quasi-Model Independent Partial Wave Analysis (QMIPWA). The samples used corresponds to the data collected in Run 2 of the LHCb experiment (between the years 2016-2018) with proton-proton collisions at a center-of-mass energy of 13 TeV. The analysis is performed after a selection process to reduce background, leading to a final sample of about 100 thousand decays. To perform the amplitude analysis fit to the Dalitz plot of the decay, the remaining background and the efficiency across the Dalitz plot are parametrized. In the case of the Isobar Model, the results obtained show a dominance of the phi(1020) resonance and a composition of scalar resonances, f0(980) as the most clear one, being formed as intermediate states, and a small but not negligible contribution of spin-2 resonances. The Isobar Model offers a qualitative description of the decay dynamics but does not provide a good quality fit, even when a variety of possible intermediate states are added to the amplitude model. On the other hand, the QMIPWA as a tool to describe the scalar sector revealed itself difficult for the region of K−K+ mass above 1.3 GeV and sensitive to the higher-spin resonances added to the model, making it difficult and less reliable to interpret the results obtained. Even with the limitations found through both models, the analysis presented here represent an important step towards the understanding of the dynamics of this decay, opening the path to further studies with higher statistics in run 3 of LHCb.