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. Comparative analyses focused on covariance variation of the first three eigenfrequencies with correlation length and standard deviation of the random field, highlighting nuanced differences between normal and log-normal distributions. PCE exhibited distinct responses, showcasing variations in covariance with different distributions. The study culminates in eigenfrequency estimation using the proposed RSM, wherein the beam is discretised into n elements with randomly assigned elastic moduli. The mean and variance of eigenfrequencies are compared with existing methods, which represent an alternative method for achieving similar outcomes.These comparative studies provide a comprehensive understanding of how different statistical treatments and simulation methods impact the reliability and accuracy of eigenfrequency predictions in beams with random elastic properties, thus contributing valuable insights for structural analysis and design under uncertainty.