A method combining analytical and numerical solution for material homogenization is investigated to calculate 3D transient eddy current field and losses in nonlinear and anisotropic laminated core under the condition of obvious leakage flux. With this method, the high-order harmonics of field components can be considered, and the phase difference of the field components within the thickness of the sheet is ignored. The TEAM Workshop Problem 21 is used as the numerical example to verify this method.