This paper describes a fuzzy-based system for analyzing the stress intensity factors (SIFs) of three-dimensional (3D) cracks. A geometry model, i.e. a solid containing one or several 3D cracks is defined. Several distributions of local node density are chosen, and then automatically superposed on one another over the geometry model by using the fuzzy knowledge processing. Nodes are generated by the bucketing method, and ten-coded quadratic tetrahedral solid elements are generated by the Delaunay triangulation techniques. The singular elements such that the mid-point nodes near crack front are shifted at the quarter-points, and these are automatically placed along the 3D crack front. The complete finite element(FE) model is generated, and a stress analysis is performed. The SIFs are calculated using the displacement extrapolation method. To demonstrate practical performances of the present system, semi-elliptical surface cracks in a inhomogeneous plate subjected to uniform tension are solved.