In the present paper it is aimed to perform the stochastic dynamic analysis of fluid and fluidstructure systems by using the Lagrangian approach. For that reason, variable-number-nodes twodimensional isoparametric fluid finite elements are programmed in Fortran language by the authors and incorporated into a general-purpose computer program for stochastic dynamic analysis of structure systems, STOCAL. Formulation of the fluid elements includes the effects of compressible wave propagation and surface sloshing motion. For numerical example a rigid fluid tank and a dam-reservoir interaction system are selected and modeled by finite element method. Results obtained from the modal analysis are compared with the results of the analytical and numerical solutions. The Pacoima Dam record S16E component recorded during the San Fernando Earthquake in 1971 is used as a ground motion. The mean of maximum values of displacements and hydrodynamic pressures are compared with the deterministic analysis results.