You can also use Plaxis but you have to select the Advanced calculation mode where you should provide the SWCC. In this mode the suction stress is considered in the definition of effective stresses.
I recommend Geo-Slope in version 2007 or later, to observe in the module SIGMA/W: in page 119 of the SIGMA/W Engineering Book, especially how is introduced a volumetric water content funtion.