Stress intensity factors for asymmetric branched cracks in plane extension by using crack-tip displacement discontinuity elements

2005 ◽  
Vol 32 (4) ◽  
pp. 375-384 ◽  
Author(s):  
Xiangqiao Yan
2005 ◽  
Vol 72 (3) ◽  
pp. 330-340 ◽  
Author(s):  
Xiangqiao Yan

Based on the analytical solution of Crouch to the problem of a constant discontinuity in displacement over a finite line segment in an infinite elastic solid, in the present paper, the crack-tip displacement discontinuity elements, which can be classified as the left and the right crack-tip elements, are presented to model the singularity of stress near a crack tip. Furthermore, the crack-tip elements together with the constant displacement discontinuity elements presented by Crouch and Starfied are used to develop a numerical approach for calculating the stress intensity factors (SIFs) of general plane cracks. In the boundary element implementation, the left or the right crack-tip element is placed locally at the corresponding left or right crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. The method is called the hybrid displacement discontinuity method (HDDM). Numerical examples are given and compared with the available solutions. It can be found that the numerical approach is simple, yet very accurate for calculating the SIFs of branched cracks. As a new example, cracks emanating from a rhombus hole in an infinite plate under biaxial loads are taken into consideration. The numerical results indicate the efficiency of the present numerical approach and can reveal the effect of the biaxial load on the SIFs. In addition, the hybrid displacement discontinuity method together with the maximum circumferential stress criterion (Erdogan and Sih) becomes a very effective numerical approach for simulating the fatigue crack propagation process in plane elastic bodies under mixed-mode conditions. In the numerical simulation, for each increment of crack extension, remeshing of existing boundaries is not required because of an intrinsic feature of the HDDM. Crack propagation is simulated by adding new boundary elements on the incremental crack extension to the previous crack boundaries. At the same time, the element characters of some related elements are adjusted according to the manner in which the boundary element method is implemented.


1990 ◽  
Vol 57 (1) ◽  
pp. 117-127 ◽  
Author(s):  
Chien-Ching Ma

The dynamic stress intensity factors of an initially stationary semi-infinite crack in an unbounded linear elastic solid which kinks at some time tf after the arrival of a stress wave is obtained as a function of kinking crack tip velocity v, kinking angle δ, incident stress wave angle α, time t, and the delay time tf. A perturbation method, using the kinking angle δ as the perturbation parameter, is used. The method relies on solving simple problems which can be used with linear superposition to solve the problem of a kinked crack. The solutions can be compared with numerical results and other approximate results for the case of tf = 0 and give excellent agreement for a large range of kinking angles. The elastodynamic stress intensity factors of the kinking crack tip are used to compute the corresponding fluxes of energy into the propagating crack-tip, and these results are discussed in terms of an assumed fracture criterion.


2019 ◽  
Vol 9 (17) ◽  
pp. 3581 ◽  
Author(s):  
Jin-Rae Cho

This paper presents the numerical prediction of stress intensity factors (SIFs) of 2-D inhomogeneous functionally graded materials (FGMs) by an enriched Petrov-Galerkin natural element method (PG-NEM). The overall trial displacement field was approximated in terms of Laplace interpolation functions, and the crack tip one was enhanced by the crack-tip singular displacement field. The overall stress and strain distributions, which were obtained by PG-NEM, were smoothened and improved by the stress recovery. The modified interaction integral M ˜ ( 1 , 2 ) was employed to evaluate the stress intensity factors of FGMs with spatially varying elastic moduli. The proposed method was validated through the representative numerical examples and the effectiveness was justified by comparing the numerical results with the reference solutions.


Author(s):  
A G Philipps ◽  
S Karuppanan ◽  
N Banerjee ◽  
D A Hills

Crack tip stress intensity factors are found for the problem of a short crack adjacent to the apex of a notch, and lying perpendicular to one of the notch faces. Loading is represented by the two Williams eigensolutions, the ratio between which provides a reference length scale and permits a comprehensive display of the solution. The results are applied to the problem of a crack starting from the edge of a notionally adhered complete contact, and conditions for the avoidance of crack development are found.


1987 ◽  
Vol 22 (4) ◽  
pp. 203-207 ◽  
Author(s):  
M H Aliabadi ◽  
D P Rooke ◽  
D J Cartwright

In order to compute stress intensity factors accurately, the standard boundary element method is modified to take explicit account of the singularity in the stresses at a crack-tip. The known expansion terms of the crack tip displacement and stress fields are subtracted to remove the numerical difficulties associated with the representation of a singular stress field at the crack-tip. Hence the accuracy of calculation is much improved, without appreciably increasing the amount of computation involved. Furthermore, the stress intensity factor is directly obtained as a part of a solution and no extrapolations are required. The improved formulation is applied to a configuration, which is representative of a part of the wing in a civil transport aeroplane. This configuration consists of a pair of circular cut-outs (supply ports) near to which smaller holes exist; these small holes are particularly susceptible to cracking.


2011 ◽  
Vol 214 ◽  
pp. 192-196 ◽  
Author(s):  
Al Emran Ismail ◽  
Ahmad Kamal Ariffin ◽  
Shahrum Abdullah ◽  
Mariyam Jameelah Ghazali ◽  
Ruslizam Daud

This study presents a numerical investigation on the stress intensity factors (SIF), K of surface cracks in round bars that were obtained under pure torsion loadings or mode III. ANSYS finite element analysis (FEA) was used to determine the SIFs along the crack front of surface cracks embedded in the solid circular bars. 20-node isoparametric singular elements were used around the crack tip by shifting the mid-side node ¼-position close to a crack tip. Different crack aspect ratio, a/b were used ranging between 0.0 to 1.2 and relative crack depth, a/D were ranged between 0.1 to 0.6. Mode I SIF, KI obtained under bending moment was used to validate the proposed model and it was assumed this proposed model validated for analyzing mode III problems. It was found that, the mode II SIF, FII and mode III SIF, FIII were dependent on the crack geometries and the sites of crack growth were also dependent on a/b and a/D.


Sign in / Sign up

Export Citation Format

Share Document