In this paper, two techniques are proposed for determining the large displacement statistics of random exciting continua endowed with fractional derivative elements: Boundary Element Method (BEM) based Monte Carlo simulation; and Statistical Linearization (SL). The techniques are applied to the problem of nonlinear beam and plate random response determination in the case of colored random external load. The BEM is implemented in conjunction with a Newmark scheme for estimating the system response in the time domain in conjunction with repeated simulations, while SL is used for estimating efficiently and directly, albeit iteratively, the response statistics.