Single crystal nickel base superalloy turbine blades are being utilized in rocket engine turbopumps and turbine engines because of their superior creep, stress rupture, melt resistance and thermomechanical fatigue capabilities over polycrystalline alloys. Currently the most widely used single crystal nickel base turbine blade superalloys are PWA 1480/1493 and PWA 1484. These alloys play an important role in commercial, military and space propulsion systems. High Cycle Fatigue (HCF) induced failures in aircraft gas turbine and rocket engine turbopump blades is a pervasive problem. Blade attachment regions are prone to fretting fatigue failures. Single crystal nickel base superalloy turbine blades are especially prone to fretting damage because the subsurface shear stresses induced by fretting action at the attachment regions can result in crystallographic initiation and crack growth along octahedral planes. Furthermore, crystallographic crack growth on octahedral planes under fretting induced mixed mode loading can be an order of magnitude faster than under pure mode I loading. This paper presents contact stress evaluation in the attachment region for single crystal turbine blades used in the NASA alternate Advanced High Pressure Fuel Turbo Pump (HPFTP/AT) for the Space Shuttle Main Engine (SSME). Single crystal materials have highly orthotropic properties making the position of the crystal lattice relative to the part geometry a significant factor in the overall analysis. Blades and the attachment region are modeled using a large-scale 3D finite element (FE) model capable of accounting for contact friction, material orthotrophy, and variation in primary and secondary crystal orientation. Contact stress analysis in the blade attachment regions is presented as a function of coefficient of friction and primary and secondary crystal orientation, Stress results are used to discuss fretting fatigue failure analysis of SSME blades. Attachment stresses are seen to reach peak values at locations where fretting cracks have been observed. Fretting stresses at the attachment region are seen to vary significantly as a function of crystal orientation. Attempts to adapt techniques used for estimating fatigue life in the airfoil region, for life calculations in the attachment region, are presented. An effective model for predicting crystallographic crack initiation under mixed mode loading is required for life prediction under fretting action.