Methods and systems using a numerical model to describe polymeric material properties are disclosed. FEM model of a product is defined. FEM model includes one or more solid elements of polymeric material. In a time-marching simulation of the product under loads, stress state of the solid elements is calculated from deformation gradient tensors. Stress state incorporates the Mullins effect and strain hardening effect, also includes elastic stress, viscoelastic stress and back stress. A yield surface is defined to determine whether the elements are under plastic deformation. Plastic strain is obtained to update the deformation gradient tensor, which is then used to recalculate the stress state. Calculations continue until updated stress state is within a tolerance of the yield surface, at which time the results of polymeric material elements are obtained. The numerical model takes into account all characteristics of a polymeric material.