The current industrial practice used at the preliminary design stage of complex structures involves the use of multifidelity submodelling simulations to predict failure behaviour around geometric and structural design features of interest, such as bolts, fillets, and ply drops. A simplified global model without the design features is first run and the resulting displacement fields are transferred to multiple local models containing the design features of interest. The creation of these high-fidelity local feature models is highly expert dependent, and their subsequent simulation is highly time-consuming. These issues compound as these design features are typically repetitive in complex structures. This leads to long design and development cycles. Application of machine learning to this framework has the potential to capture a structural designer’s modelling knowledge and quickly suggest improved design feature parameters, thereby addressing the current challenges. In this work, we provide a proof of concept for a machine learning assisted preliminary design workflow, see Figure 1, whereby feature-specific surrogate models may be trained offline and used for faster and simpler design iterations. The key challenge is to maximise the prediction accuracy of failure metrics whilst managing the high dimensions required to represent design feature simulation parameters in a minimum training dataset size. These challenges are addressed using: a modified Latin Hypercube Sampling scheme adjusted to improve design of experiment in composite materials; a bi-linear work-equivalent homogenisation scheme to reduce the number of nodal degrees of freedom; a non-local volume-averaged stress-based approach to reduce the number of target features; and linear superposition of stacked bi-directional LSTM neural network models. This methodology is demonstrated in a case study of predicting the stresses of open hole composite laminates in an aerospace C-spar structure. Results highlight the high accuracy (>90%) and time saving benefit (>15x) of this new approach. This methodology may be used to faster correct and iterate the preliminary design of any large or complex structure where there are repetitive localised design features that may contribute to failure, such as in Formula 1 or wind turbines. Combined with exascale computing this methodology may also be applied for predictive virtual testing of digital twins.