We discuss how the mixed model and the asymmetric mixed model family of bivariate extreme value can be extended to bivariate extreme value copulas with polynomial dependence function of arbitrary degree. An algorithm for fitting extreme value copulas with polynomial dependence functions to data will be presented and various practical issues that arise when fitting bivariate extreme value copula models will be discussed.