Abstract
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 three-dimensional equilibrium equations for shell structures in orthogonal curvilinear coordinates is introduced. The theoretical framework for the two-dimensional shell model is based on a unified formulation able to study and compare several Higher-order Shear Deformation Theories (HSDTs), including Murakami's function for the so-called zig-zag effect. Thus, various higher-order models are used and compared also to investigate the differences which can arise from the choice of the order of the kinematic expansion. Their ability to deal with several damaged configurations is analyzed as well. The paper can be placed also in the field of numerical analysis, since the solution to the static problem at issue is achieved by means of the Generalized Differential Quadrature (GDQ) method, whose accuracy and stability are proven by a set of convergence analyses and by the comparison with the results obtained through a commercial finite element software.