For the last three decades, bistable composite laminates have gained publicity because of their outstanding features, including having two stable shapes and the ability to change these states. A common challenge regarding the analysis of these structures is the high computational cost of existing analytical methods to estimate their natural frequencies. In the current paper, a new methodology combining the Finite Element Method (FEM) and Multi-Objective Genetic Programming (MOGP) is proposed for the analysis of bistable composite structures, leading to some analytical relations derived to obtain the modal parameters of the shells. To achieve this aim, the data extracted from FEM, consisting of the ratio of the length to width (a/b) and the thickness (t) of the laminate, is split into Train and Validation, and Test, subsets. The former is used in MOGP, and four formulas are proposed for the prediction of the free vibration parameters of bistable laminates. The formulas are checked against the Test subset, and the statistical indices are calculated. An excellent performance is observed for all GP formulas, which indicates the reliability and accuracy of the predictions of these models. Parametric studies and sensitivity analyses are conducted to interpret the trend of input parameters in the GP models and the level of sensitivity of each natural frequency formula to the input parameters. These explicit mathematical expressions can be extended to the other bistable laminates to obtain their natural frequencies on the basis of their geometrical dimensions. The results are validated against the experimental data and verified against FEM outcomes. @en