An extremly accurate and efficient numerical method for solving the problem was presented, which is mainly by means of the crack isolating analysis technique, stress superposition principle, the Legendre polynomial expansion of the pseudo-traction as well as the segmental average collocation technique. The singular equations of the Cauchy type and Fredholm integral equations of the first type were deduced. In the process of dealing with the superposition of infinite number of kinked cracks, the crack boundary conditions are satisfied. Many complex computing examples were given, and for some typical examples, numerical results were compared with the analytic solutions and the numerical solutions obtained by a boundary element method. The numerical results show that the stress intensity factors depend on the crack configuration, and on the geometrical and physical parameters.