We extend a recently-developed framework for isogeometric analysis of composite laminates to drive material damage evolution with a smoothed strain field. This builds on ideas from gradient-enhanced continuum damage modeling, and is intended to limit the dependence of damage predictions on the choice of discrete mesh. The resulting enhanced framework models each lamina of a composite shell structure as a Kirchhoff–Love thin shell. To account for the anisotropic damage modes of laminae, we smooth a tensor-valued strain by solving an elliptic partial differential equation (PDE) system on each lamina. This strain-smoothing PDE system is formulated to be independent of the choice of coordinates and is applicable to general manifold shell geomet...