Numerical Solution of the Wave Propagation Problem in a Plate

In this work, we use the phase velocity method in combination with finite element method to compute the dispersion curve for phase velocity of an ultrasonic pulse traveling in a thin isotropic plate. This method is based on the numerical solution of the wave propagation equations for several selected frequencies. To solve these equations, a second order difference scheme is used to discretize the temporal variable, while spatial variables are discretized using the finite element method. The variational formulation of the problem corresponding to a fixed value of time is obtained and the existence and uniqueness of the solution is proved. A priori error estimates in the energy norm and in the L2 norm are also obtained. The open software FreeFem++ is used with quadratic triangular elements to compute the displacements. Numerical experiments show that the velocities computed from the approximated displacements for different frequency values are in good agreement with analytical dispersion curve. This confirms that the proposed symbiosis between finite element and phase velocity method is suitable for computing dispersion curves in more general wave propagation problems, where the geometry is complex and the material is anisotropic.