The present invention relates to a finite element analysis method to implement the surface tension, and more particularly to node (node) ID number to N from 1 for the N nodes constituting the boundary surface of the constituting the analysis target medium and a first step of assigning values ​​in sequence, and storage by extracting the coordinate values, the unique number from the node of the unique number n n-1, and a unique number claim calculating each of the unit direction vectors for the node of the n + 1 2 phase with, the nodes of a unique number n and, integral with the third step of calculating respective distances between a unique number n-1, and a unique number n + 1 node of the surface tension with respect to the calculated distance in a third step value the fourth step and, by multiplying the Gaussian integral of said fourth step to a unit direction vector of the second stage Co., for implementing the surface tension to a fifth step of storing in the external value of the node of the serial number n to It relates to an element analysis method.