The general approach of numerical treatment of integro-differential equation of the flat crack problem is considered. It consists in presenting the crack surface loading as the set of the polynomial functions of two Cartesian coordinates while the corresponding crack surface displacements are chosen as the similar polynomials multiplied by the function of form (FoF) which reflects the required singularity of their behavior. To find the relations matrixes between these two sets a new effective numerical procedure for the integration over the area of arbitrary shape crack is developed. In based on the classical hyper-singular method, i.e. Laplace operator is initially analytically applied to the integral part of equation and the resulting hyper singular equation is subsequently considered.
The presented approach can be implemented with any variant of FoF, but Oore-Burns FoF, which was earlier suggested in their famous 3D weight function method, is supposed to be the most accurate and universal. It takes into account all points of crack contour, which provides perfect physical conditionality of the solution, but such FoF is relatively heavy in implementation and of low computational speed. The special procedure is developed for the approximation of the crack contour of arbitrary shape by the circular and straight segments. It allows to easily obtain analytical expression for Oore-Burns FoF, which greatly increases the calculation speed and accuracy. The accuracy of the considered method is confirmed by the examples of the circular, elliptic, semicircular and square cracks at different polynomial laws of loading.
The developed methods are used in the implemented procedure for crack growth simulation. It allows to model growth of crack of arbitrary shape at arbitrary polynomial loading, at that all contour points are taken into account and can expand with their own speeds each. Procedure has high accuracy and don’t need complex and high-cost re-meshing process between the iterations unlike FEM or other numerical methods. At that usage of Oore-Burns FoF provides high flexibility of the presented approach: unlike similar theoretical methods, where FoF calculation procedure is rigidly connected with the crack shape, which complicates the adequate crack growth modeling, the used FoF automatically takes into account all points of crack contour, even if its shape became complex during the growth. Presented crack growth procedure can be effectively used to test accuracy and correctness of correspondent numerical methods, including the newest XFEM approach.