“…In general, flat 3-node triangular shell element cannot well approximate the pure bending displacement field; thus, locking problems are serious and inevitable. As the shell thickness decreases, the convergence of the finite element solution rapidly deteriorates [8]. Researchers developed methods to overcome locking phenomena in flat triangular shell elements, such as the Mixed Interpolation of Tensorial Components method used in the MITC3 element [8], the line integration method used in TRIA3 element [42], the mixed or hybrid formulation [43][44][45], incompatible displacement methods [46][47], stabilization methods [29,[48][49][50][51][52][53][54], assumed strain methods [19,28,55], etc.…”