This paper is devoted to the theoretical and numerical study of singular perturbation problems for elliptic inhibited shells. We present a reduction of the classical membrane equations to a partial differential equation with respect to the bending displacement, which is well adapted to the study of singularities of the limit problem. For a discontinuous loading or when the boundary of the loading domain presents corners, we put in a prominent position the existence of two kinds of singularities. One of them is not classical; it reduces to a logarithmic point singularity at the corner of the loading domain. To finish numerical simulations are performed with a finite element software coupled with an anisotropic adaptive mesh generator. They enable to visualize precisely the singularities predicted by the theory with only a very small number of elements.