Stress intensity factors computations using the singularity subtraction technique incorporated with the Tau Method
Abstract
In this paper we discuss the use of the singularity subtraction technique incorporated with the Tau Method for the numerical solution of singular partial differential equations which are relevant to the linear elastic fracture mechanics. To treat the singularity, we apply the singularity subtraction technique to the singular boundary value problems. The problems arising in this application are not in the standard form required by the Tau software. By introducing the pseudo-differential equations, Ak~ = 0, k = l(l)m, to determine the stress intensity and higher order factors Ak results in the standard boundary value problems. We consider two model crack problems including Motz' anti-plane crack problem and a plane strain problem defined by the biharmonic equation. We obtain results of considerable accuracy which compare favorably with those published in the recent literature.
Keywords
References
[1] U. Ascher, R.D. Russell. Reformulation of boundary value problems into standard form. SIAM Review, 23: 238- 254, 1981.[2] M.J.M. Bernal, J .R. Whiteman. Numerical treatment of biharmonic boundary value problems with re-entrant boundaries. Comp. J., 13: 87-91, 1970.
[3] S.K. Chan, I.S. Tuba, W.K. Wilson. On the finite element method in linear fracture mechanics. Engng. Fracture Mech., 2: 1-17, 1970.
[4] B. Gross, J.E. Srawley, W.F. Brown Jr. Stress-intensity factors for a single-edge-notch tension specimen by" boundary collocation of a stress function. Technical Note D-2395, N.A.S.A., Washington, D.C., 1964.
[5] C. Lánczos. Trigonometric interpolation of empirical and analytical functions. J. Math. Phys., 17: 123-199, 1938.
This work is licensed under a Creative Commons Attribution 4.0 International License.