Neoclassical transport simulation code (FORTEC-3D) applicable to non-axisymmetric configurations is developed. Adoption of a new hybrid simulation model, in which ion transport is solved by using the delta f Monte-Carlo method including the finite-orbit-width effects while electron transport is solved by a reduced ripple-averaged kinetic equation, makes it possible to simulate the dynamism of non-local transport phenomena with self-consistently developing radial electric field within a allowable computation time. Time evolution of radial electric field in LHD plasma is simulated in the full volume of confinement region, and the finite-orbit-width effect of neoclassical transport is found to make the negative ambipolar electric field more larger than the prediction by a local transport theory