The eigenfunction expansion method is introduced into the numerical calculations of elastic plates. Based on the variational method, all the fundamental solutions of the governing equations are obtained directly. Using eigenfunction expansion method, various boundary conditions can be conveniently described by the combination of the eigenfunctions due to the completeness of the solution space. The coefficients of the combination are determined by the boundary conditions. In the numerical example, the stress concentration phenomena produced by the restriction of displacement conditions is discussed in detail.