This study carried out fiber damage detections of laminated GFRP plate structures using a modified bi-variate Gamma function. The effects of different layup sequences of composites on the fiber damage detection are studied using the finite element commercial package and genetic algorithm. Four unknown parameters are considered to determine the shape of the damage distribution, which is a modified form of the bivariate Gamma density distribution function. The sample studies show the excellence of the proposed method from the standpoints of its computation efficiency as well as its ability to determine the complex shape of an arbitrary stiffness degradation distribution.