A mathematical scheme is proposed here to model a damaged mechanical configuration for laminated and sandwich structures. In particular, two kinds of functions defined in the reference domain of plates and shells are introduced to weaken their mechanical properties in terms of engineering constants: a two-dimensional Gaussian function and an ellipse shaped function. By varying the geometric parameters of these distributions, several damaged configurations are analyzed and investigated through a set of parametric studies. The effect of a progressive damage is studied in terms of displacement profiles and through-the-thickness variations of stress, strain, and displacement components. To this end, a posteriori recovery procedure based on the ...