In this study, a three-dimensional nonlinear finite element (FE) model of a female pubic symphysis was developed using automatic mesh generation techniques for tetrahedral and hexahedral elements. The geometry was based on CT scan data of a female cadaver pelvis. The model was composed of cortical and trabecular bone, a fibrocartilaginous interpubic disc, and four ligaments connecting the pubic bones. Cortical and trabecular bone were assigned linear elastic properties, while properties for the soft tissues were estimated heuristically until the overall symphysis structural behaviors agreed with the results of previously published biomechanical experiments. The nonlinear hyperelastic James-Green-Simpson constitutive material equation was used to capture axial tension and compression, while a two-term Prony series was fit to the experimental creep data. The accuracy of the estimated material properties was further studied by comparing model predictions with experimental data from slow and fast tensile tests. The results demonstrated that the two-term Prony series extension of the James-Green-Simpson material model captured the nonlinear viscoelastic behavior of the pubic symphysis to within 6% of the experimental data. Element performance tests indicated that quadratic tetrahedrons and linear hexahedrons performed roughly the same in terms of solution accuracy and computation cost.