Substructuring based parametric reduced order modelling for large-scale dynamical systems containing viscoelasticity with application to bonded assemblies

•Parametric model order reduction strategy in context of substructuring proposed for large scale mechanical systems with viscoelasticity.•Krylov subspace based on the approximation of viscoelasticity replaces the interior and interface reduction of the Craig-Bampton method with local characteristic constraint modes.•New expression of the frequency-dependent constraint modes.•Parametric global reduction basis obtained by combination of singular value decomposition and interpolation.•Efficiency and accuracy of the proposed method demonstrated with two numerical examples.