In this study, an embedded fiber is considered in matrix under pull-out force. Then, two intensities of singular stress fields (ISSF) are discussed appearing at the fiber end A and the intersection point E of the fiber and the surface. The analysis method focuses on the FEM stress at points A and E by applying the same FEM mesh pattern to the unknown and reference 2D models. To analyze the ISSF at A, the body force method solution is used as the reference model. To analyze the ISSF at E, the reciprocal work contour integral method (RWCIM) solution is used as the reference model. Then, the two ISSFs are compared by varying the fiber bonded length lin. When lin is shorter, the singular stress at A is larger than the singular stress at E. When lin is longer, the singular stress at E is larger than the singular stress at A.
