Dielectric elastomer (DE) actuators are deformable capacitors capable of a muscle-like actuation when charged. When subjected to voltage, DE membranes coated with compliant electrodes may form wrinkles due to the Maxwell stress. Here, we develop a numerical approach based on the finite element method (FEM) to predict the morphology of wrinkled DE membranes mounted on a rigid frame. The approach includes two steps: (I) pre-buckling and (II) post-buckling. In step I, the first buckling mode of the DE membrane is investigated by substituting the Maxwell stress with thermal stress in the built-in function of the FEM platform simulia abaqus. In step II, we use this first buckling mode as an artificial geometric imperfection to conduct the post-buckling analysis. For this purpose, we develop an equivalent model to simulate the mechanical behavior of DEs. Based on our approach, the thickness distribution and the thinnest site of the wrinkled DE membranes subjected to voltage are investigated. The simulations reveal that the crests/troughs of the wrinkles are the thinnest sites around the center of the membrane and corroborate these findings experimentally. Finally, we successfully predict the wrinkles of DE membranes mounted on an isosceles right triangle frame with various sizes of wrinkles generated simultaneously. These results shed light on the fundamental understanding of wrinkled dielectric elastomers but may also trigger new applications such as programmable wrinkles for optical devices or their prevention in DE actuators.