The results indicate that viscoelastic-GK calculation with relaxation of shear modulus and of bulk modulus using the multilayer shell element is the most appropriate method for predicting the warpage.