In reliability studies, we are interested in the behaviour of a system when it interacts with its surrounding environment. To assess the system’s behaviour in a reliability sense, we can take the system’s intrinsic quality as strength and the outcome of interactions as stress. Failure is observed whenever stress exceeds strength. Taking Y as a random variable representing the stress the system experiences and random variable X as its strength, the probability of not failing can be taken as a proxy for the reliability of the component and given as P(Y<X)=1−P(X<Y). This way, in the present paper, it is considered that X and Y follow generalized extreme value distributions, which represent a family of continuous probability distributions that have been extensively applied in engineering and economic contexts. Our contribution deals with a more general scenario where stress and strength are not independent and copulas are used to model the dependence between the involved random variables. In such modelling framework, we explored the proper selection of copula models characterizing the dependence structure. The Gumbel–Hougaard, Frank, and Clayton copulas were used for modelling bivariate data sets. In each case, information criteria were considered to compare the modelling capabilities of each copula. Two economic applications, as well as an engineering one, on real data sets are discussed. Overall, an easy-to-use methodological framework is described, allowing practitioners to apply it to their own research projects.