Abstract: The iterative formulas of elastoplastic stochastic finite element method (SFEM) under cyclic loading are deduced and the random responses of local multiaxial stress and strain are calculated. According to the complication of cyclic loading the anisotropy and Bauschinger effect resulting from plastic deformation are reflected in the kinematic hardening model. The transient stress-strain relation of Jhansale model is considered in the iterative formulas. The elastoplastic finite element method (EFEM) overcomes the limitation that only the uniaxial local stress and strain could be calculated by the approximate method, establishing the base of multiaxial fatigue analysis. The component randomness is considered. The SFEM is introduced in the iterative formulas of EFEM under cyclic loading to obtain the random response of local stress and strain. The result of Mont Carlo simulation demonstrates the efficiency of the proposed method.