A modular Gaussian beam model is developed to simulate some ultrasonic testing configurations where multiple interfaces are involved. A general formulation is given in a modular matrix form to represent the Gaussian beam propagation with multiple interfaces. The ultrasonic transducer fields are modeled by a multi-Gaussian beam model which is formed by superposing 10 single Gaussian beams. The proposed model, referred to as "MMGB" (modular multi-Gaussian beam) model, is then applied to a typical contact and angle beam testing configuration to predict the output signal reflected from the corner of a vertical crack. The resulting expressions given in a modular matrix form are implemented in a personal computer using the MATLAB program. Simulation results are presented and compared with available experimental results.