ABSTRACT. Many materials used in engineering applications obey to time-dependent behaviours and the mechanical fields are affected by the time effects. As a result, the evolution of the stresses and strains in these materials appear still very complex and difficult to study. Among such cases is the situation when the material has an axisymmetric shape and when it is submitted to a complex fracture loading. In this paper, the creep loading is applied on an axisymmetric viscoelastic orthotropic material and the stress intensity factors are computed in the opening mode, in the shear mode and in the mixed mode using to a finite element approach. The uncoupling method is based on M integral, combining the virtual and real mechanical fields. In the same time, the viscoelastic effects are introduced according to the generalized Kelvin-Voigt model composed by four branches. The numerical solution is obtained with an incremental viscoelastic formulation in the time domain. Using a Compact Tension Shear (CTS) specimen, the evolutions of stress intensity factor versus time are posted in each fracture mode configuration. The obtained results demonstrate the efficiency of the proposed model.