In this work the warpage of a power modules ceramic substrates due to temperature variation has been numerically calculated. It has been used a not linear finite element model, which account the experimentally characterized material properties, including the copper elastoplasticity. Model results have been compared with dedicated interferometric measurements to validate the model. As application example, it has been calculated the benefit in terms of temperature warpage realized with a different substrate design option.