Computational models, particularly finite-element (FE) models, are essential for interpreting experimental data and predicting system behavior, especially when direct measurements are limited. A major challenge in tuning these models is the large number of parameters involved. Traditional methods, such as one-by-one sensitivity analyses, are time-consuming, subjective, and often return only a single set of parameter values, focusing on reproducing averaged data rather than capturing the full variability of experimental measurements. In this study, we applied simulation-based inference (SBI) using neural posterior estimation (NPE) to tune an FE model of the human middle ear. The training dataset consisted of 10,000 FE simulations of stapes velocity, ear-canal (EC) input impedance, and absorbance, paired with seven FE parameter values randomly sampled within plausible ranges. The neural network learned the association between parameters and simulation outcomes, returning the probability distribution of parameter values that can reproduce experimental data. Our approach successfully identified parameter sets that reproduced three experimental datasets simultaneously. By accounting for experimental noise and variability during training, the method provided a probability distribution of parameters, representing all valid combinations that could fit the data, rather than tuning to averaged values. The network demonstrated robustness to noise and exhibited an efficient learning curve due to the large training dataset. SBI offers an objective alternative to laborious sensitivity analyses, providing probability distributions for each parameter and uncovering interactions between them. This method can be applied to any biological FE model, and we demonstrated its effectiveness using a middle-ear model. Importantly, it holds promise for objective differential diagnosis of conductive hearing loss by providing insight into the mechanical properties of the middle ear.