In this paper, the exponentiated Kumaraswamy-power function distribution is introduced. Some structural properties of the proposed distribution including the shapes of the density, hazard and quantile functions are explored. Besides, skewness and kurtosis measures are investigated. The method of maximum likelihood is used for estimating model parameters. For different parameter settings and sample sizes, a simulation study is performed and the performance of the new distribution is compared with some well-known distributions. Then, an application is presented with a real data set to illustrate the usefulness of the proposed distribution.