This paper presents a novel mechanical environmental test method that multi-output random strain responses are replicated instead of acceleration responses. In this work, the positive-definite strain reference spectral matrix and reference kurtoses for multiple control channels are defined first, and then three types of strain random signals including stationary Gaussian, stationary non-Gaussian, and non-stationary are generated. The inverse Fourier transform is used to obtain the stationary Gaussian strain random signals and the stationary non-Gaussian strain signals are obtained by the zero-memory nonlinear transformation. To synthetize non-stationary strain random signals, an amplitude modulation function is introduced with time-varying root mean square values. The correction algorithms are employed to update the drive signals for matching the strain responses compared with corresponding reference values. The simulation and experimental results prove the effectiveness and feasibility of the proposed strain response control procedure. The presented work provides new perspectives to investigate the fatigue behavior and life estimation of structure under controlled strain response environments for vibration fatigue tests.