Finite element models ranging from single to multiscale models have been widely used to gain valuable insights into the physiological delivery process and associated complication scenarios. However, the fetus descent simulation with the active uterine contraction is still challenging for validation and uncertainty quantification issues. The present study performed a fetus descent simulation using the active uterine contraction. Then, simulation outcomes were evaluated using theoretical and in vivo MRI childbirth data. Moreover, parameter uncertainty and propagation were also performed. A maternal pelvis model was developed. The active uterine contraction was modeled using a transversely isotropic Mooney–Rivlin material. Displacement trajectories were compared between simulation, theoretical and in vivo MRI childbirth data. Monte Carlo (M.C) and Polynomial Chaos Expansion (PCE) methods were applied to quantify uncertain parameters and their propagations. Obtained results showed that fetal descent behavior is consistent with the MRI-based observation as well as the theoretical trajectory (curve of Carus). The head downward vertical displacement ranges from 0 to approximately 47 mm. A reduction of 50% in uterine size was observed during the simulation. Three high-sensitive parameters ( C 1 , C 2 , Ca 0 ) were also identified. Our study suggested that the use of the active uterine contraction is essential for simulating vaginal delivery but the global parameter sensitivity, parameter uncertainty, and outcome evaluation should be carefully performed. As a perspective, the developed approach could be extrapolated for patient-specific modeling and associated delivery complication simulations to identify risks and potential therapeutic solutions.