With regards to the composition of natural rocks including voids or pores, deformation behavior is strongly affected by variation in porosity. By using a statistical damage-based approach, the characteristics of strain softening and hardening under the influence of voids and volume changes are investigated in the present paper. Suppose that a rock consists of three parts: voids, a damaged part, and an undamaged part. The effects of voids and volume changes on rock behavior are first analyzed through determination of the porosity and an associated damage model is then developed. Later, a statistical evolution equation describing the influence of the damage threshold on the propagation condition of rock damage is formulated based on measurement of the mesoscopic element strength. A statistical damage constitutive model reflecting strain softening and hardening behavior for rocks loaded in conventional triaxial compression is further developed and a corresponding method for determining the model parameters is also provided. Theoretical results of this proposed model are then compared with those observed experimentally. Finally, several aspects of the present constitutive model, which affect the relevant behavior of rocks, are particularly discussed.