Abstract: A three dimensional penalty equilibrium hybrid element method is presented in this paper. Suggestions on the selection of penalty number are made. Numerical experiment shows that the method is stable, reliable and efficient without incurring additional computational cost. In particular, it is effective in improving the accuracy of numerical results in distorted grids.