The Fourier transform of $y\left(x,t\right)$ from the time to the frequency domain is given by $Y\left(x,\omega\right)=\int_{-\infty}^{\infty}y\left(x,t\right)e^{i\omega t}dt$ and satisfies the differential equation: $$ EI\frac{\partial^{4}Y\left(x,\omega\right)}{\partial x^{4}}+\mu\left(i\omega\right)^{2}Y\left(x,\omega\right)=\int_{-\infty}^{\infty}P\delta\left(x-u\left(t\right)\right)e^{i\omega t}dt $$
Now, use this property of the delta function: $$ \delta\left(g\left(s\right)\right)=\sum_{i}\frac{\delta\left(s-s_{i}\right)}{\left|g'\left(s_{i}\right)\right|} $$
where $s_{i}$ are the roots of the function, so that $g\left(s_{i}\right)=0$.
$$ \int_{-\infty}^{\infty}P\delta\left(g\left(t\right)\right)e^{i\omega t}dt=\int_{-\infty}^{\infty}P\sum_{i}\frac{\delta\left(t-t_{i}\right)}{\left|g'\left(t_{i}\right)\right|}e^{i\omega t}dt=P\sum_{i}\frac{e^{i\omega t_{i}}}{\left|g'\left(t_{i}\right)\right|} $$
Then we have: $$ \frac{\partial^{4}Y\left(x,\omega\right)}{\partial x^{4}}-\frac{\mu\omega^{2}}{EI}Y\left(x,\omega\right)=\frac{P}{EI}\sum_{i}\frac{e^{i\omega t_{i}}}{\left|g'\left(t_{i}\right)\right|} $$
The Green's function solution to this equation is: $$ Y\left(x,\omega\right)=\int_{-\infty}^{\infty}\frac{P}{EI}\sum_{i}\frac{e^{i\omega t_{i}}}{\left|g'\left(t_{i}\right)\right|}G\left(x,x'\right)dx' $$ The Green's function satisfies the equation: $$ \frac{\partial^{4}G\left(x,x'\right)}{\partial x^{4}}-\psi^{4}G\left(x,x'\right)=\delta\left(x-x'\right) $$
Now we have to invert the Fourier transform. 818