The development of the liquid element that is capable of modeling the liquid behavior in flexible liquid storage tanks is achieved throughout several stages. First, the classical isoparametric liquid element, which is also used to model the liquid motion inside the liquid domain, is developed. Then, the element is enhanced with the nonlinear sloshing and liquid-structure interaction in order to model the liquid boundaries.