Directional anisotropic hardening under non-proportional loading is modeled. A recently proposed coupled quadratic-nonquadratic (CQN) yield function (Lee et al., 2017b) successfully captured the directional hardening behavior under the proportional loadings. This paper shows a constitutive equation to capture both directional hardening response and Bauschinger effect through an extended model of the CQN yield function with the homogeneous anisotropic hardening (HAH) model (Barlat et al., 2011). In the present study, the role of the CQN yield model is to capture the directional hardening behavior, and the HAH model is to follow non proportional loading including the Bauschinger effect. The validation, with some material data, shows that the present model can follow the directional hardening under non-proportional loading in the stress-strain data. It is also shown that the present model can describe the yield surface distortion with Bauschinger and anisotropic hardening response. Finally, it is compared with other pre-existing models in order to analyze the utility of the proposed model.