Dam operation and management have become more complex recently because of the need for considering hydraulic structure sustainability and environmental protect on. An Earthfill dam that includes a powerhouse system is considered as a significant multipurpose hydraulic structure. Understanding the effects of running hydropower plant turbines on the dam body is one of the major safety concerns for earthfill dams. In this research, dynamic analysis of earthfill dam, integrated with a hydropower plant system containing six vertical Kaplan turbines (i.e., Haditha dam), is investigated.In the first stage of the study, ANSYS-CFX was used to represent one vertical Kaplan turbine unit by designing a three-dimensional (3-D) finite element (FE) model. This model was used to differentiate between the effect of turbine units’ operation on dam stability in accordance to maximum and minimum reservoir upstream water levels, and the varying flowrates in a fully open gate condition. In the second stage of the analysis, an ANSYS-static modeling approach was used to develop a 3-D FE earthfill dam model. The water pressure pattern determined on the boundary of the running turbine model is transformed into the pressure at the common area of the dam body with turbines. The model is inspected for maximum and minimum upstream water levels. Findings indicate that the water stress fluctuations on the dam body are proportional to the inverse distance from the turbine region. Also, it was found that the cone and outlet of the hydropower turbine system are the most affected regions when turbine is running. Based on the attained results, a systematic operation program was proposed in order to control the running hydropower plant with minimized principal stress atselected nodes on the dam model and the six turbines.