This paper proposes a computationally technique for simulating solutions of nonlinear nonlocal stochastic dynamical systems driven by variable-order fractional Brownian motion with Hurst index. The value of the Hurst index depends on time t belong to interval ( 1/2,1). The proposed technique is adopted quadratic interpolation for fractional-order derivative. Moreover, it is exploited in the discussion of fractional stochastic financial and pendulum dynamical systems. The proficiency of the presented technique is confirmed by using of investigating statistical indicators for the stochastic approximations for various values of fractional order parameters.