Accurate prediction of dynamical response of structural system depends on the correct modeling of that system. However, modeling becomes increasingly challenging when there are many candidate models available to describe the system behavior. Furthermore, uncertainties can be present even for the parameters of these model classes. The plausibility of each input-output model class of the structures with uncertain components can be determined by a Bayesian approach from measured dynamic responses to one or more input records; predictions of the structural system response to alternate input records can then be made. However, this approach may require many model simulations, even though most of those model classes are quite implausible. An approach is proposed herein to use a bound, computed from the false discovery rate, on the likelihood of measured data to falsify models considering uncertainties in the passive control devices that do not reproduce the measured data to sufficient accuracy. Response prediction is then performed using the unfalsified models in an approximate Bayesian sense by assigning weights, computed from the likelihoods, only to the unfalsified models approach incurring only a fraction of the computational cost of the standard Bayesian approach. The proposed approach for response prediction is illustrated using three structural examples: an earthquake-excited four--degree-of-freedom building model with a hysteretic isolation layer; a 1623--degree-of-freedom three-dimensional building model, with tuned mass dampers attached to its roof, subjected to wind loads; and a full-scale four-story base-isolated building tested on world's largest shake table in Japan's E-Defense lab. The results exhibit accurate response predictions and significant computational savings, thereby illustrating the potential of the proposed method.