Green's functions for elastic solids with random properties are usually derived by means of the perturbation method. This paper deals with a new approach that has the potential to deal with a large variability of random shear modulus based on the transformation of a polynomial chaos. The deterministic Green's functions for stresses and displacements and the principal values in the boundary integrals caused by a pressure load on the surface are nonlinear transformations of the random variables. A series of transformations of the polynomial chaos is used to transform significant parts of the equation. The first operation is a projection of the log normal distributed shear modulus to a series of Hermite's polynomials based on a Gaussian variable. The second operation is the determination of an arbitrary potential of the wave velocity. The last operation, similar to the first one, consists in the determination of an exponential function depending on the inverse of the wave velocity. These operations, together with multiplications and summations, transform the complete relation from the random shear modulus to Green's functions and principal values. The inversion of the system matrix is already derived for the random finite element approach. The operations are independent of the specific problem and can be applied to almost all acoustic media and similar nonlinear problems.