The conservative elastic behavior of soft materials is characterized by a stored energy function which shape is usually specified a priori, except for some material parameters. There are hundreds of proposed stored energies in the literature for different materials. The stored energy function may change under loading due to damage effects, but it may be considered constant during unloading–reloading. The two dominant approaches in the literature to model this damage effect are based either on the Continuum Damage Mechanics framework or on the Pseudoelasticity framework. In both cases, additional assumed evolution functions, with their associated material parameters, are proposed. These proposals are semi-inverse, semi-analytical, model-driven and data-adjusted ones. We propose an alternative which may be considered a non-inverse, numerical, model-free, data-driven, approach. We call this approach WYPiWYG constitutive modeling. We do not assume global functions nor material parameters, but just solve numerically the differential equations of a set of tests that completely define the behavior of the solid under the given assumptions. In this work we extend the approach to model isotropic and anisotropic damage in soft materials. We obtain numerically the damage evolution from experimental tests. The theory can be used for both hard and soft materials, and the infinitesimal formulation is naturally recovered for infinitesimal strains. In fact, we motivate the formulation in a one-dimensional infinitesimal framework and we show that the concepts are immediately applicable to soft materials.