Stress intensity factor solutions for two-dimensional elastostatic problems by the hypersingular boundary integral equation

2009 ◽  
Vol 44 (4) ◽  
pp. 235-247 ◽  
Author(s):  
A Sahli ◽  
O Rahmani

The boundary element method is a numerical method that can be advantageously used for a wide range of engineering problems, including the stress concentration problems encountered in fracture mechanics. In linear elastic fracture mechanics (LEFM), the stress intensity factor is an important parameter. Cracks, if present in the region experiencing the modes of deformation, increase the stress amplitude significantly and this high stress may lead to premature failure of the engineering components. If the value of the SIF is known, it is possible to predict whether the crack will propagate or not. As the conventional boundary integral equation (CBIE) degenerates when a mathematical crack is modelled, a previously developed dual boundary integral equation approach has been adopted in the current work. It utilizes the hypersingular boundary integral equation (HBIE) along with the CBIE. A weakly singular form of HBIE is utilized in the current work to eliminate the hypersingularity analytically. Stress intensity factors are evaluated using the crack opening displacement (COD), displacement extrapolation (DE), and the J-integral approaches. A stand-alone code has been developed for calculating the stress intensity factors of general two-dimensional domains with cracks. The code has been validated by evaluating the stress intensity factors for the standard components, for which the stress intensity factor values are available in the literature. Accurate and well-converged results are obtained proving the robustness of the code. A linear combination of the CBIE and HBIE was applied at the crack and a significant (87–97 per cent) reduction in the condition numbers for the system of equations was observed for the examples studied. Again, the results obtained are accurate and well converged.

1988 ◽  
Vol 110 (3) ◽  
pp. 408-413 ◽  
Author(s):  
L. J. Ghosn

Crack propagation in a rotating inner raceway of a high-speed roller bearing is analyzed using the boundary integral method. The model consists of an edge plate under plane strain condition upon which varying Hertzian stress fields are superimposed. A multidomain boundary integral equation using quadratic elements was written to determine the stress intensity factors KI and KII at the crack tip for various roller positions. The multidomain formulation allows the two faces of the crack to be modeled in two different subregions making it possible to analyze crack closure when the roller is positioned on or close to the crack line. KI and KII stress intensity factors along any direction were computed. These calculations permit determination of crack growth direction along which the average KI times the alternating KI is maximum.


1991 ◽  
Vol 113 (1) ◽  
pp. 22-27 ◽  
Author(s):  
J. L. Desjardins ◽  
D. J. Burns ◽  
R. Bell ◽  
J. C. Thompson

Finite elements and two-dimensional photoelasticity have been used to analyze thick-walled cylinders which contain arrays of straight-fronted, longitudinal-radial cracks of unequal depth. The stress intensity factor K1 has been computed for the dominant crack and for some of the surrounding cracks. Cylinders with 2, 4, 6, 8, 16, 36 and 40 cracks have been considered. Good agreement has been obtained between the experimental and the numerical results and, for cylinders with 2 or 4 cracks, with previously published predictions. The results for all of the foregoing cases are used to develop simple, approximate techniques for estimating K1 for the dominant crack, when the total number of cracks is different from those that have been considered herein. Estimates of K1 obtained by these techniques agree well with corresponding finite element results.


Author(s):  
D. J. Shim ◽  
S. Tang ◽  
T. J. Kim ◽  
N. S. Huh

Stress intensity factor solutions are readily available for flaws found in pipe to pipe welds or shell to shell welds (i.e., circumferential/axial crack in cylinder). In some situations, flaws can be detected in locations where an appropriate crack model is not readily available. For instance, there are no practical stress intensity factor solutions for circular-arc cracks which can form in circular welds (e.g., nozzle to vessel shell welds and storage cask closure welds). In this paper, stress intensity factors for circular-arc cracks in finite plates were calculated using finite element analysis. As a first step, stress intensity factors for circular-arc through-wall crack under uniform tension and crack face pressure were calculated. These results were compared with the analytical solutions which showed reasonable agreement. Then, stress intensity factors were calculated for circular-arc semi-elliptical surface cracks under the lateral and crack face pressure loading conditions. Lastly, to investigate the applicability of straight crack solutions for circular-arc cracks, stress intensity factors for circular-arc and straight cracks (both through-wall and surface cracks) were compared.


Author(s):  
Mayumi Ochi ◽  
Kiminobu Hojo ◽  
Itaru Muroya ◽  
Kazuo Ogawa

Alloy 600 weld joints have potential for primary water stress corrosion cracks (PWSCC). At the present time it has been understood that PWSCC generates and propagates in the Alloy 600 base metal and the Alloy 600 weld metal and there has been no observation of cracking the stainless and the low alloy steel. For the life time evaluation of the pipes or components the crack extension analysis is required. To perform the axial crack extension analysis the stress intensity database or estimation equation corresponding to the extension crack shape is needed. From the PWSCC extension nature mentioned above, stress intensity factors of the conventional handbooks are not suitable because most of them assume a semi-elliptical crack and the maximum aspect ratio crack depth/crack half length is one (The evaluation in this paper had been performed before API 579-1/ASME FFS was published). Normally, with the advance of crack extension in the thickness direction at the weld joint, the crack aspect ratio exceeds one and the K-value of the conventional handbook can not be applied. Even if those equations are applied, the result would be overestimated. In this paper, considering characteristics of PWSCC’s extension behavior in the welding material, the axial crack was modeled in the FE model as a rectangular shape and the stress intensity factors at the deepest point were calculated with change of crack depth. From the database of the stress intensity factors, the simplified equation of stress intensity factor with parameter of radius/thickness and thickness/weld width was proposed.


Sign in / Sign up

Export Citation Format

Share Document