A numerical method based on integration along bicharacteristics is presente
d for analysing two-dimensional stress wave propagation problem in layered
orthotropic plates. This method is an extension of the method proposed by C
lifton for isotropic plane problem. The numerical stability is discussed. T
he method is applied to simulate stress wave propagation in a semi-infinite
plate composed of two different orthotropic materials.