Contact mechanics based solution to predict modulus of asphalt materials with high porosities