The aim of this study is to simulate the contraction phenomena in pregnant uterine. Differential equations are used to model the uterine contractions for three states, namely contracting, refractory, and resting. Since verification and validation are vital for simulation systems, we compared our model with uterine magnetomyogram (MMG) data which are recorded non-invasively using SQUID-Array for Reproductive Assessment (SARA) system at University of Arkansas for Medical Sciences (UAMS). It is possible to identify the contraction locations from the data via Hilbert-Huang Transform (HHT). Two simulations were done; first one was based on the mathematical model and the second one was based on the contraction data. Comparison between those simulations revealed similar characteristics.