An improved numerical (finite element) procedure is developed to investigate the
stress field in symmetrically laminated composites of finite dimensions. Emphasis is placed
on assessing the singular behavior of stresses in regions close to ply-interfaces