Enhanced models for the nonlinear bending of planar rods: localization phenomena and multistability