MR tagging is considered as a valuable technique to evaluate regional myocardial function quantitatively and noninvasively, however the cumbersome and time-consuming post-processing procedures for cardiac motion tracking still hinder its application to routine clinical examination. We present a fast and semiautomatic method for tracking 3D cardiac motion from short-axis (SA) and long-axis (LA) tagged MRI images. The technique, called 3D-HARP (HARmonic Phase), is based on the HARP1 method and extends this method to track 3D motion. A material mesh model is built to represent a collection of material points inside the left ventricle (LV) wall. The phase time-invariance property of material points is used to track the mesh points. For a series of 9 timeframe MRI images, the total time required for initializing settings, building the mesh, and tracking 3D cardiac motion is approximately 10 minutes. Further analysis of Langrangian strain and twist angle demonstrates that during systole, the lateral LV wall shows a greater strain values than the septum and the SA slices from the base to the apex show a gradual change in twist pattern.