Abstract: Fracture of Kirchhoff plates is analyzed by the theory of complex variables and boundary collocation method. The deflections, moments and shearing forces of the plates are assumed to be the functions of complex variables. The functions can satisfy a series of basic equations and governing conditions, such as the equilibrium equations in the domain, the boundary conditions on the crack surfaces and stress singularity at the crack tips. Thus, it is only necessary to consider the boundary conditions on the external boundaries of the plate, which can be approximately satisfied by the collocation method and least square technique. Different boundary conditions and loading cases of the cracked plates are analyzed and calculated. Compared to other methods, the numerical examples show that the present method has many advantages such as good accuracy and less computer time. This is an effective semi-analytical and semi-numerical method.