Most of the previous research has been focused on developing and validating constitutive models to predict response of unsaturated soils in low-medium suction range. However, there is a scarcity of efforts in developing soil models to simulate its mechanical response in high suction range, particularly above the residual suction. This article presents a new constitutive model introducing net stress and suction as two independent variables. Furthermore, non-associative flow rule incorporating modified stress-dilatancy relationship to take unsaturated state into account is introduced to improve the model results in low-medium to high suction range. The essential soil model parameters are calibrated using suction-controlled triaxial test results for predictions of compacted silty sand response at high values of total suction above residual suction. Preliminary simulations show that proposed model can reasonably simulate the post-peak strain softening response obtained from suction-controlled CTC tests above residual suction value with reasonable accuracy. Although, the proposed model captures initial compression followed by dilation volumetric response with reasonable accuracy, it needs some improvements to be able to capture volumetric response accurately over entire suction range.