A novel procedure for large deformation nonstationary random response computation of shell structures with spatial uncertainty is presented. The procedure is free from the limitations associated with those employing perturbation approximation techniques, such as the so-called stochastic finite element method and probabilistic finite element method, for systems with spatial uncertainties. In addition, the procedure has several important and excellent features. Chief among these are: (a) ability to deal with large deformation problems of finite strain and finite rotation; (b) application of explicit linear and nonlinear element stiffness matrices, mass matrix, and load vectors reduces computation time drastically; (c) application of the averaged deterministic central difference scheme for the updating of co-ordinates and element matrices at every time step makes it extremely efficient compared with those employing the Monte Carlo simulation and the conventional central difference algorithm; and (d) application of the time co-ordinate transformation enables one to study highly stiff structural systems.