A kinematic hardening model including an associated flow rule is proposed for elastoplastic deformations in graphitic grey cast iron. Quantitatively good results are obtained when comparing with previously performed biaxial experiments. Use of a nonassociated flow rule is found to result in an undesirable weakening behavior that can be explained as a deficiency with the combination of kinematic hardening and the present choice of yield potential. The model proposed is also extended to include multilinear kinematic hardening. With this model qualitatively good agreement with experimental cyclic results from the literature is obtained. A three-dimensional FE-analysis of a cylinder head for a heavy duty Diesel engine is performed as an application. To predict initiation of thermal fatigue cracks, it is essential to use an elastoplastic material model.