This study is based on the numerical simulation of the mechanical response of yttrium-stabilized zirconia ceramic (3Y-TZP) dental implants as a function of their intrinsic geometry and masticatory loads. Samples (n=20) of 3Y-TZP ceramics were compacted, sintered at 1500 °C -2h, and characterized by relative density, X-Ray diffraction (XRD), and scanning electron microscopy (SEM). The elastic parameters (modulus of elasticity and Poisson ratio), used in the numerical simulations, were measured by the Impulse Excitation Technique, and the bending strength was obtained using piston-on-three-balls testing. An authorial implant design and, comparatively, commercial implant CAD models were used in this study as an initial geometry of dental implant in a typical adult mandible anatomy. From CAD and CAE techniques, finite element models were generated for all implant geometries. Loading cases were considered based on different intensities (100N to 500N) and orientation angles (45° or 90°) to reproduce the human masticatory efforts. The numerical predictions were compared with finite element simulations of gold-standard titanium-based implants. The investigated 3Y-TZP sintered ceramics presented high densification (> 99%), with a microstructure formed by submicron equiaxed tetragonal zirconia grains. The 3Y-TZP average bending strength obtained from piston-on-three-balls testing is 1192 ± 99 MPa. For both dental implant geometries, the zirconia implants showed average strength of less than 550 MPa, which, in turn, is independent of the masticatory load value or orientation angle. All finite element predictions are 50% inferior to the corresponding measured flexural strength values and preliminarily enable the 3Y-TZP ceramics for dental implant applications without fracture risk.