Viscoelasticity is a material property that is relevant in a variety of nanoscale materials and interfaces in medicine and industry. Therefore, a method of mechanical quantification has become exceedingly desired. In this thesis the Atomic force microscope (AFM) is applied to accurately characterize the mechanical behavior of viscoelastic samples. The goal is to enhance viscoelastic characterization using the so-called Intermodulation AFM (ImAFM) technique by applying, adapting and improving multiple modelling and optimization methods. In ImAFM force reconstruction is performed by extracting intermodulations around resonance in the cantilever response. These intermodulations present new observables that can be used for characterization. This thesis investigates the potential of this technique in combination with an up-and-coming model describing viscoelastic interaction. A toolbox has been developed for numerical simulations of the model to resemble the experiments. The model has been evaluated in a variety of situations using sensitivity analysis in a large feasibility range, encompassing many complex dynamics. Because of the diversity in model dynamics a global optimization has been performed for experimental reconstruction.