This article investigates the influence of random elastic modulus on beam eigenfrequencies using multiple simulation techniques: Monte Carlo simulations (employing Cholesky decomposition (MCS-CD) and Kosambi-Karhunen-Loève expansion (MCS-KKL)), Polynomial Chaos expansion (PCE), and a proposed Random Sampling method (RSM). Anomalies in Monte Carlo simulations, where normally distributed elastic modulus led to negative values and imaginary eigenfrequencies, were effectively addressed by adopting a log-normal distribution...