Is there a solution to the expected value/variance for a Gaussian with random phase:

$$\cos(\omega_0 t + \phi), \qquad \phi \sim \cal{N}(0,\sigma^2) $$

?

For $t=0$, the solution is for example given here: Resultant probability distribution when taking the cosine of gaussian distributed variable