Numerous powerful methods exist for developing reduced-order models (ROMs) from finite element (FE) models. Ensuring the accuracy of these ROMs is essential; however, the validation using dynamic responses is expensive. In this work, we propose a method to ensure the accuracy of ROMs without extra dynamic FE simulations. It has been shown that the well-established implicit condensation and expansion (ICE) method can produce an accurate ROM when the FE model’s static behaviour are captured accurately. However, this is achieved via a fitting procedure, which may be sensitive to the selection of load cases and ROM’s order, especially in the multi-mode case. To alleviate this difficulty, we define an error metric that can evaluate the ROM’s fitting error efficiently within the displacement range, specified by a given energy level. Based on the fitting result, the proposed method provides a strategy to enrich the static dataset, i.e. additional load cases are found until the ROM’s accuracy reaches the required level. Extending this to the higher-order and multi-mode cases, some extra constraints are incorporated into the standard fitting procedure to make the proposed method more robust. A curved beam is utilised to validate the proposed method, and the results show that the method can robustly ensure the accuracy of the static fitting of ROMs.