We present a model of fracture formation on surfaces of bilayered materials. The model makes it possible to synthesize patterns of fractures induced by growth or shrinkage of one layer with respect to another. We use the finite element method (FEM) to obtain numerical solutions. This paper improves the standard FEM with techniques needed to efficiently capture growth and fractures. Reference P. Federl and P. Prusinkiewicz: Finite element model of fracture formation on growing surfaces. In M. Bubak, G. van Albada, P. Sloot and J. Dongarra (Eds.): Proceedings of Computational Science. ICCS 2004 (Krakow, Poland, June 6−9, 2004), Part II, Lecture Notes in Computer Science 3037, Springer, Berlin, pp. 138-145.