To build a fluid-structure interaction model of pulmonary veins with total anomalous pulmonary venous connection, which can be used to predict the result of surgical treatment and at the same time to estimate the elastic properties of pulmonary veins based on patient-specific data and clinic postoperative results. The fluid-structure interaction (FSI) model was used to simulate the anastomosis on pulmonary veins based on computed tomography angiography data collected from three children with total anomalous pulmonary venous connection (TAPVC), supra-cardiac type. The deformation and the stress of anastomosis, and also the velocity of blood flow were calculated in fluid-structure coupling algorithm. During the simulation the variable boundary conditions were applied, including the thickness of vessel wall and the vessel elasticity for which was selected a range of values. The calculation results were finally compared to postoperative results of same patients and discussed. The corresponding outcomes are given in the conclusions section. The blood flow velocity through the outlet will vary depending on the properties of vessels, including physical properties and thickness of vessel wall. The stress on vessel is lower for smaller values of Young's modulus. The calculated blood flow velocity correlates well with the postoperative results for the Young's modulus of vessels ranging from 0.5 to 1.0 MPa. The FSI model has high potential to predict the result of surgery for TAPVC and to estimate the physical properties of pulmonary vein. This model also has potential to guide the strategy for surgical treatment.