Congenital heart malformations occur to substantial number of pregnancies. Studies showed that abnormal flow biomechanical environments could lead to malformations, making it important to understand the biomechanical environment of the developing heart. We performed 4D high-frequency ultrasound scans of chick embryonic hearts at HH25 to study the biomechanics of the whole heart (atria and ventricle). A novel and high-fidelity motion estimation technique, based on temporal motion model and non-rigid image registration algorithm, allowed automatic tracking of fluid-structure boundaries from scan images, and supported flow simulations. Results demonstrated that atrial appendages were the most contractile portion of the atria, having disproportionately high contribution to atrial blood pumping for its volume in the atria. However, the atria played a small role in blood pumping compared to the ventricle, as it had much lower ejection energy expenditure, and as the ventricle appeared to be able to draw inflow from the veins directly during late diastole. Spatially and temporally averaged wall shear stresses (WSS) for various cardiac structures were 0.062-0.068 Pa, but spatial-averaged WSS could be as high as 0.54 Pa in the RV. WSS was especially elevated at the atrial inlet, atrioventricular junction, regions near to the outflow tract, and at dividing lines between the left and right atrium and left and right side of the ventricle, where septation had begun and the lumen had narrowed. Elevated WSS could serve as biomechanics stimulation for proper growth and development.