Dynamic substructuring methods encompass a range of techniques, which allow for the decomposing of large structural systems into multiple coupled subsystems. This decomposition of structures into smaller domains has the principle benefit of reducing computational time for dynamic simulation of the system by considering multiple smaller problems rather than a single global one. In this context, dynamic substructuring methods may form an essential component of hybrid simulation, wherein they can be used to couple physical and numerical substructures at reduced computational cost. Since most engineered systems are inherently nonlinear in nature, particular potential lies in incorporating nonlinear treatment in existing sub-structring schemes, which are mostly developed within a linear problem scope.The most widely used and studied dynamic substructuring methods are classical linear techniques such as the Craig-Bampton and MacNeal-Rubin methods. These methods are widely studied and implemented in many commercial FE packages. However, as linear methods they naturally break down in the presence of nonlinearities. Recent advancements in substructuring have involved the development of enrichments to the linear substructuring methods, which allow for some degree of nonlinearity to be captured. The use of mode shape derivatives has been shown to be able to capture geometrically non-linear effects as an extension to the Craig-Bampton method. Other candidates include the method of Finite Element Tearing and Interconnecting. Linear substructuring methods have been used in several cases for hybrid simulation, rendering additional benefits in removing non-physical high frequency modes which may be excited due to controller tracking errors in hybrid simulations.In this work, a virtual hybrid simulation is presented in which a linear elastic vehicle frame supported on four nonlinear spring damper isolators is decomposed into separate domains. One domain consisting of the finite element model of the vehicle frame, which is reduced using the Craig-Bampton method to only include modes relevant to the structural response. The second domain consists of the nonlinear isolators whose restoring forces are characterised by nonlinear spring and damper forces. Coupling between the models is carried out using a Lagrange multiplier method and time series simulations of the system are conducted and compared to the full global system with regards to simulation time and accuracy.