A computational method for wave propagation from a point load in an anisotropic material

One approach which is employed to solve dynamic point load problems in plates and laminates is to take integral transforms to reduce the governing equations to a system of ordinary differential equations with respect to the depth variable. The solution of this system leads to expressions for the transforms of the displacement and stress components at any level in the plate and the transient response at any location may then be recovered by inversion of the multiple transforms. The formal transform inversion involves a double infinite integral but by making a change of variable this may be replaced by an infinite integral associated with a line source and a finite integral with respect to the orientation of the line. A first attempt at applying this approach to obtain the point load response of quasi-isotropic fibre composite laminate led to a non-causal predicted signal. This paper deals with an investigation of this proposed method applied to the simpler model problem of wave propagation in a two-dimensional anisotropic medium. Results are obtained for two different time histories of point loads, namely: a delta function; and a single period of a sine function. In the case of the delta function source a comparison is made with the analytic solution and the errors arising from the numerical approach are discussed. Graphs are also presented showing the non-causal contributions to the overall response which arise at individual angles of orientation of the line source.