The standard material point method (MPM) encounters severe numerical difficulties in simulating shell structures. In order to overcome the shortcomings of locking effects, the discretization size of background grid should be small enough, usually smaller than 1/5 of the shell thickness, which however will lead to prohibitive computational cost. A novel solid shell material point method (SSMPM) is proposed to efficiently model the large deformation of thin structures. The SSMPM describes the material domain of shell structures by shell particles with hexahedral particle domains. The locking treatments of solid shell element are then introduced in SSMPM, which results in the correction of strain field throughout the shell thickness. Namely, the assumed natural strain (ANS) method is adopted to eliminate the shear locking and trapezoidal locking, while the enhanced assumed strain (EAS) method is employed to eliminate the thickness locking. With the precise description of bending modes, a single layer of particles and a coarse background grid could be used in shell structure simulations with the SSMPM, which dramatically increases the computational efficiency. A local multi‐mesh contact method is presented to naturally couple SSMPM and MPM for the contact situations of shells with other objects. Several numerical examples, including beam vibration, pinched cylinder with free edges and full hemispherical shell, are performed to verify and validate the SSMPM, which shows that the SSMPM considerably outperforms the standard MPM in these situations. A fluid–structure interaction problem and the penetration of a thin plate are investigated based on the contact method and the results are in good agreement with those in the literature.