An algorithm for finite deformation plasticity theory, recently proposed by Simo [5] and Simo and Meschke [6], is applied to numerical simulations of localized failure in snow. Within the algorithmic framework for finite strain plasticity the classical Drucker-Prager model and the modified Cam-Clay model are employed for the representation of snow. The respective material parameters are calibrated according to results from hydrostatic and shear-box tests of snow specimen. The effectiveness of the FE-model is demonstrated by means of re-analyses of the shear-box test, characterized by the formation of a sharp, localized failure zone. The traction mechanism of a single automobile tread moving over a snow-covered road surface is numerically simulated. The predicted failure mechanism of snow in consequence of the movement of the tread agrees relatively well with respective experimental observations. INTRODUCTORY REMARKS This paper is concerned with the numerical investigation of the traction behavior of automobile tires on snow-covered road surfaces by means of the finite element method (FEM). This type of analysis involves realistic modelling of rubber, snow and the contact at the rubber-snow interface. The traction behavior of tires for winter conditions is characterized by a complex mechanical interaction between the treads of the tires and the snow. The snow in the vicinity of the treads is subjected to large plastic deformations associated with the appearance of localized failure zones. The occurance of large non-recoverable deformations suggests using the of classical plasticity theory, embedded in the geometrically nonlinear theory of finite strain Transactions on Engineering Sciences vol 6, © 1994 WIT Press, www.witpress.com, ISSN 1743-3533 462 Localized Damage plasticity, for the modelling of snow. Relatively little research work concerning 3D constitutive models for snow is found in the open literature. In a recent work by Mohamed, Yong and Murcia [1] it is assumed that snow follows a ,/2-type of plasticity law. However, this simplified model does not account for the considerable plastification of snow when subjected to hydrostsatic pressure. Salm [4] proposed a viscoelastic constitutive law restricted to uniaxial stress states. The present numerical investigation has been accompanied by an experimental investigation of snow specimens involving hydrostatic and shear-box tests. The results of this series of tests suggest the 'critical state concept', widely used in geomechanics, to be an adequate representation of the constitutive behavior of snow. In the present analyses, the (modified) Cam-Clay model, recast within a recently proposed algorithmic framework for finite strain plasticity [5], [6] is employed. This algorithm is formulated on the basis of hyperelasticity in the current configuration, with the elastic region being defined in terms of Kirchhoff stresses. The key advantage of this method is that the structure of the return mapping algorithm of the infinitesimal theory can be taken over to the nonlinear theory without modification. A fully implicit algorithmic treatment of the nonlinear elastic constitutive law along with an implicit integration of the flow rule on the basis of the return mapping algorithm is used. The Cam-Clay model is calibrated by means of the hydrostatic and the shear-box tests, respectively. For comparison, an elastic, ideally-plastic Drucker-Prager model is adopted. The important issue of viscosity of snow is not considered in this paper. It will be addressed in a follow-up publication. Both models for snow are applied in the re-analyses of shear-box tests. As far as the Drucker-Prager model is concerned, a comparison between the new plasticity algorithm with the standard incremental hypoelastic formulation, implemented in the multipurpose program MARC is presented. The obtained results are compared with experimental data. Numerical simulations of the traction mechanism of a single automobile tread moving on a snow-covered surface, based on the Cam-Clay model for snow, are performed. REMARKS ON FINITE STRAIN PLASTICITY The extension of the classical (modified) Cam-Clay plasticity model to the finite strain regime is based on an algorithm for finite strain plasticity recently proposed by Simo [5] and Simo and Meschke [6]. The key ingredients of this algorithm are briefly outlined below: The stress response is governed by a hyperelastic law, characterized by a function of stored energy, W(F ) involving the elastic part of the deformation gradient, F*. The restriction to isotropy implies that W(F*) = M (b ), where b* is the elastic part of the left Cauchy-Green tensor. The total deTransactions on Engineering Sciences vol 6, © 1994 WIT Press, www.witpress.com, ISSN 1743-3533