A gaussian beam based recursive stiffness matrix model to simulate ultrasonic array signals from multi-layered media