A calculation method of 
sound field distribution of a 
phased array ultrasonic probe in flaw detection in a multilayer medium. The method comprises the steps of: S1, modifying a Rayleigh integral expression to obtain an expression of 
phased array sound field transmission; S2, establishing a 
fast algorithm of path in an arbitrary multilayer medium; and S3, calculating 
chip delay time and excitation 
delay time: S3-1, establishing a three-dimensional 
system of coordinates comprising a probe 
wafer, a wedge and a work piece; S3-2, conducting 
Fourier transform on a drive 
signal; S3-3, forming a path matrix; S3-4, calculating a 
point source delay time, 
diffusion attenuation, and reflection and transmission through each section; S3-5, calculating 
sound field amplitude and phase transformation generated from a 
target signal, and conducting inverse transformation to obtain sound field from each 
point source drive 
signal transmitting to the target calculation point; S3-6, superimposing to obtain sound field value at the target point; and S3-7, dividing units and sequentially calculating to obtain sound field distribution in the target area. The invention can effectively assist 
process design and defect evaluation 
phased array detection of complex multilayer structure.