ISRN Mechanical EngineeringVolume 2014 (2014), Article ID 834659, 7 pageshttp://dx.doi.org/10.1155/2014/834659

Research Article

Effect of Density and Pointed Corner Degree of Pore on Local Stress in Welded Structures: Defect in Marine Structures

Farida Bouafia,1,2 Boualem Serier,¹ Nassim Serier,¹ and Sardar Sikandar Hayat³

¹LMPM, Mechanical Engineering Department, University of Sidi Bel Abbes, 22000 Sidi Bel Abbes, Algeria²Institute of Science and Technology, University of Ain Temouchent, BP 284 RP, 46000 Ain Temouchent, Algeria³Department of Physics, Hazara University, Mansehra 21300, Pakistan

Received 30 July 2013; Accepted 23 October 2013; Published 16 January 2014

Academic Editors: N. Anifantis, R. Brighenti, and S. Cartmell

Copyright © 2014 Farida Bouafia et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The process of assembly by welding, in marine structures, led to the creation of microstructural heterogeneities zones. Welded structures generally contain defects such as lack of penetration, slag inclusion, porosity, and misalignment. Generally these zones can be discontinuities geometrical. At the origins of stress concentration, these regions are favourable sites where fatigue cracks can initiate and propagate. In this study, the 3D finite element method is used to analyze the distribution of stress (strain) generated by the presence in the weld nugget of a pore formed during the welding process. This analysis was made in the matrix-pore interface. The effects of density and pointed corner degree of pore on the stress levels were also analyzed.

1. Introduction

The major weld imperfections, in marine structures, include porosity, lack of penetration, lack of fusion, slag inclusion, undercut, misalignment, and so forth [1]. The size and frequency of imperfection depend on the welding process, weld procedure, geometry of weldment including ease of access for welding, and the care exercised in making the weld. These imperfections have different characteristics and, in most cases, are difficult and costly to detect and define nondestructively.

These weld defects can significantly influence the local stress field in the vicinity of welds when the welded component is subjected to cyclic fatigue load. In most cases, weld defects lead to severe stress concentration and thus accelerate fatigue crack growth. In agreement with many authors, fatigue and fracture are the most critical failure mode of welded joints [2–4]. Hence, it is up of importance to study the stress distribution around these zones.

This behavior is attributed to the fact that for a given fatigue life, a much larger embedded imperfection can be tolerated than a surface imperfections [5]. Different welds imperfections as described by Maddox [1] are shown in Figure 1.

Figure 1: Different weld imperfections in a butt-joint [1].

Porosity is defined as cavity type discontinuities (voids) formed by gas entrapment during solidification. The shape of the void is dependent on the relative rates of solidification of the weld metal and the nucleation of the entrapped gas. The resultant stress field surrounding the pore depends upon the pore shape and the loading.

Porosity, though generally spherical in shape, can assume many shapes and configurations. These include elongated pores, rows of single pore or collinear pores, and pore clusters. Determining the effects of various sizes and shapes of pores is an important factor affecting the structural integrity of weldments.

Fatigue cracks generally initiate at a geometrical discontinuity such as a notch or pore. These act as stress concentrations, raising the stress in the region of the notch to levels above the nominal stresses. The material at the notch root may deform plastically while the rest of the component remains essentially elastic. Subjecting the region to cyclic loading resulting in plastic deformation will eventually result in a fatigue crack.

Fatigue crack propagation in marine structures is obviously governed by the mechanics of considerably different four types of multiscale problems. Problems of structural response to environmental loads have a length scale of several hundred meters, whereas the possible detectable size of cracks from initial defects in a weld is of the order of millimeters [6].

Many of the researchers have studied the welding defects [7–9]. Bouafia et al. [9] utilizes a numerical analysis for the computation of stress concentration factor generated by the presence in the weld nugget of a pore formed during the welding process. Welded structure containing porosity is subjected to uniaxial tensile stress. The effects of geometrical parameters of the pore and the interaction pore-defect on the stress concentration factor variation have been analyzed by [9].

In this study, a finite element model is developed to analyze numerically the distribution of local stresses generated by the presence of a pore formed during the welding process. The analysis has been extended to the effect of density and pointed corner degree of pore on the stress levels.

Through this study, we show that the formation of porosity due welding process can weaken the weld by notch effect. In the real case air bubbles or moisture trapped during welding may lead to the formation of porosity.

2. Finite Element Model

The distribution of the pores in the matrix was modeled using unit cells in which the pores are placed perfectly uniformly within the matrix and unit cells containing large numbers of randomly spaced pores (Figure 2(a)). For the unit cells with the uniformly spaced pore arrangement, the pores were packed in a square array as shown in Figure 2(b). This packing is the most common arrangement used in many so-called “unit cell” analyses [9]. In addition to the uniformly spaced pore arrangements in Figure 2, numerical simulations were conducted in this work using repeating square cells of the matrix which contained pore. In this investigation, a three-dimensional finite element analysis was developed. For reasons of symmetry in geometry, we studied a quarter of the complete model (Figure 2(c)) in order to reduce the calculation time. Hence, the boundary conditions, must be such that the symmetry is retained, that is, all the nodes at the faces of symmetries are constrained to have a null normal displacement (Figure 3(a)), all other boundaries were free [9]. The structure was subjected to a tensile mechanical load (Figure 3(b)). The choice of tensile stresses in the welded joint subjected is justified by the analysis of the mechanical behavior in open mode (mode 1), the most dangerous mode. During commissioning of the seal, the stresses are highly concentrated in the vicinity of this very close porosity. The initiation of the crack from a pore may be considered at the moment when the seal reaches a level of damage. After this level, this joint begins to undergo a softening behavior and gradually loses its rigidity.

(a)

(b)

(c)

Figure 2: (a) Square packing of pores within the matrix and the square periodic unit cell. (b) Unit cell, (c) quarter of unit cell.

(a)

(b)

(c)

Figure 3: (a) The three-dimensional finite element models, (b) boundary condition and loading condition, and (c) finite element models mesh.

The porosity is located in the center of the weld nugget and simulated as a hole with a diameter φ of 50 μm, characterizing an average size selected from the literature [9]. The matrix is considered as a Steel matrix with a Young modulus  GPa and a Poisson coefficient .

The finite-element computations was performed, using commercial FE software package ABAQUS [10]. The precision of numerical computation is strongly related to the quality of the mesh in the structure. In addition, due to the stress concentrations expected at the matrix-pore interface, the mesh was refined at this zone, A 8-node linear brick (C3D8) finite elements are used for the model (Figure 3(c)). The finite element model is shown in Figure 3(c) with 10000 elements, and it has a fine grid at the matrix-pore interface.

3. Results and Discussion

3.1. Stress Analysis

We studied the distribution of the von Mises stress and normal stress along -axis in the matrix and the pore-matrix interface for four different positions of path numbered Path 1, Path 2, Path 3, and Path 4 (Figure 4). This figure shows, respectively, the variation of equivalent and normal stress according to the peripheral angle θ (peripheral angle θ is the angle around the pore. It is between 0° and 90°), for different paths. The higher stresses are located at the vicinity of the interface between the matrix and pore. Far from the interface the matrix is completely free of this stress. Indeed, these stresses increase gradually the amplitude from θ = 0° to θ = 90°. These stresses are positive in the matrix and the pore-matrix interface. This clearly shows that the matrix is ​ ​in tension. These stresses are of greater intensity at the matrix-pore interface (Path 1), and far from the interface, this stress is relaxed. They are the largest in the immediate vicinity of the interface. This stress can constitute a risk of damage to the welding joint.

(a)

(b)

Figure 4: Variation of equivalent and normal stresses according to the peripheral angle θ, for different paths, for  MPa.

3.2. Strain Analysis

The stress induced locally in the matrix around the pore lead to the strain of this matrix in the vicinity close to this defect. This strain depends on the amplitude of the stress applied to the structure and thus of the level of the stress generated in the joint of welding close to the pore [9]. Variation of Max principal stress () according to the Max principal strain (), for Path 1 is represented in Figure 5. It shows clearly that the amplitude of Max principal stress grows with the increase in Max principal strain and that the level of this stress and strain reaches its maximum close to the matrix-pore interface. This behavior shows that the pore initially spherical form becomes deformed elastically for tending towards an elliptical form. Under these local pressures the pore tends to become deformed preferentially, according to its axis , towards its center of gravity [9].

Figure 5: Variation of Max principal stress () according to the Max principal strain (), for Path 1, for  MPa.

3.3. Effect of Pore Density

The pore density is the significant parameter in the analysis of the local stresses in the weld joint. We will study the effect of 1, 4, and 7 pores on the equivalent and local stresses.

Figure 6 presents the distribution of equivalent stress in the matrix and pore for two models with 4 and 7 porosities for  MPa. The higher stress is located at the vicinity of the interface between the matrix and pores. Far from the interface, the matrix is completely free of this stress. The equivalent stress leads to traction in the pores and matrix.

Figure 6: Equivalent stress distribution for 4 and 7 porosities, for  MPa.

In Figure 7, the effect of pore density on the level of the von Mises equivalent stress is represented. It shows clearly that the amplitude of this stress grows with the increase in the pore density to the level of the matrix. A more significant pore density induced stresses of strong intensities at the interface matrix-pore. In addition, reaches its highest level at 90°.

Figure 7: Variation of equivalent stress according to the peripheral angle θ, for different pore density, for path 1 and  MPa.

Figure 8 presents the distribution of normal stress in a matrix with 4 and 7 pores for  MPa. The higher stresses are located at the vicinity of the interface between the matrix and pores. The effect of pore density on the level of the normal stress is represented in Figure 9. It shows clearly that the intensity of this stress grows with the increase in pore density. The normal stress takes the highest level at the edge of the pore of θ = 90°.

Figure 8: Normal stress distribution for 4 and 7 porosities, for  MPa.

Figure 9: Variation of normal stress according to the peripheral angle θ, for different pore density, for path 1 and  MPa.

The effect of the interaction of two pores on the level of the internal stresses induced in the ferritic matrix in the center of the distance between two pores was studied by Bouafia et al. [11]. Figure 10 presents the distribution of this stress and their intensity in the ferritic matrix and two pores for  MPa. The higher stress occurs in the vicinity of the matrix-pore interface. The stress is computed along a distance . The analysis of this figure shows that the stress along axis in the matrix in the center of the distance between two pores increase progressively when the distance between the two pores decreases. This stress is strongly localised in the matrix close to the surface of the pores. The decrease of the distance between pores involves an increase of the stress.

Figure 10: Variation of normal stress as functions of the distance from the model’s center: cases of two pores, for  MPa [11].

This result shows that an increase of pores density generates a strong stress concentration in the matrix being able to initiate and propagate cracks. This behavior constitutes a damage risk of the weld joint. A very low stresses level does not present any risk of failure of the weld joint [11].

3.4. Effect of Pointed Corner Degree of Pore on Stress and Strain Distributions

In the present study, finite element method was introduced to analyze the equivalent stress and strain fields in five models for different shaped pores (Figure 11). The pointed corner degree of porosity varies from 180° to 30°. It is supposed that the five pores have equal areas.

Figure 11: A 2D schematic view of FEM model: (a) spherical pore; (b) hexagon pore; (c) square pore; (d) triangle pore; (e) shuttle-shaped pore.

Distribution of equivalent stress in the pore-matrix interface varies with the pore shape (Figure 11). As far as the spherical pore is concerned, the equivalent stress distributes uniformly and remains a low level of 106 MPa. When the pore becomes angular, there is a stress concentration in the pointed pore corner. With decreasing pointed corner degree of the pore, this concentration increases rapidly and becomes more and more intense. In the shuttle-shaped pore, the maximum equivalent stress exists at the pointed corner and is up to 301 MPa (Figure 12(a)).

(a)

(b)

Figure 12: Variation of equivalent stress and Max principal strain () according to the pointed corner degree of pore, for  MPa.

In Figure 12(b), it is obvious that the strain near the interface varies with pore shape. When pointed corner degree of pore decreases there is an increasing strain concentration in the matrix around the pointed corner. The maximum strain in the shuttle-shaped pore is much higher than that in the spherical pore where the strain distributes uniformly around the pore (Figure 12(b)).

4. Conclusions

The porosity in the steel matrix resulting from a process of welding in marine structures generates normal stress whose level reaches its maximum value with the interface matrix porosity. Far from this interface the matrix is completely released of this stress. These stresses can constitute a risk of damage to the welding joint.

The level and the distribution of the equivalent and normal stress depend on the density of porosity. These stresses are all the stronger as the density of the pore grows. The pore density has a great effect on stress distributions of the matrix-pore model.

The pore shape has a significant influence on the level of the developed stress and strain. The pointed pore corner can result in intense stress concentration in the pointed pore corner and in situ strain concentration in the matrix around the pointed pore corner.

Distribution of strain in pore-matrix interface varies with the pore shape. There is a strain concentration around the pointed corner degree of pore and this concentration increases with decreasing pointed corner degree of pore.

The distribution of stress in the pore also varies with the pore shape. There is a serious stress concentration in the pointed corner of pore and this concentration increases rapidly with decreasing pointed corner degree of pore.

List of Symbols

: | Diameter of porosity (μm)
: | Young modulus (GPa)
: | Poisson coefficient.
: | Tensile mechanical load (MPa)
: | von Mises stress (MPa)
: | Normal stress along -axis (MPa)
: | Peripheral angle (°)
: | Max principal stress (MPa)
: | Max principal strain (MPa).

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

S. J. Maddox, “Applying fitness-for purpose concepts to the fatigue assessment of welded joints,” in Proceesings of the International Conference on Fatigue, vol. 1, pp. 72–81, Toronto, Canada, 1994.

W. Fricke, “Fatigue analysis of welded joints: state of development,” Marine Structures, vol. 16, no. 3, pp. 185–200, 2003.

Y. Zhang and D. Taylor, “Sheet thickness effect of spot welds based on crack propagation,” Engineering Fracture Mechanics, vol. 67, no. 1, pp. 55–63, 2000.

H. Kang, M. E. Barkey, and Y. Lee, “Evaluation of multiaxial spot weld fatigue parameters for proportional loading,” International Journal of Fatigue, vol. 22, no. 8, pp. 691–702, 2000.

J. M. Barsom, “Fracture mechanics analysis of fatigue crack initiation and growth,” in Proceedings of the International Conference on Fatigue, pp. 88–98, Toronto, Canada, 1994.

Y. Sumi and T. Inoue, “Multi-scale modeling of fatigue crack propagation applied to random sequence of clustered loading,” Marine Structures, vol. 24, no. 2, pp. 117–131, 2011.

E. E. Gdoutos, Fracture Mechanics: An Introduction, vol. 14 of Solid Mechanics and Its Applications, Springer, 1993.

S. Liu, A. M. Pope, and R. Daemen, “Welding consumables and weldability,” in Proceedings of the International Workshop on Underwater Welding of Marine Structures, pp. 321–350, New Orleans, La, USA, December 1994.

F. Bouafia, S. Boualem, M. M. E. Amin, and B. Benali, “3-D finite element analysis of stress concentration factor in spot-welded joints of steel: the effect of process-induced porosity,” Computational Materials Science, vol. 50, no. 4, pp. 1455–1464, 2011.

ABAQUS, User's Manual, 6. 5, Hibbit, Karlsson & Sorensen.

F. Bouafia, B. Serier, B. Boutabout, and B. A. B. Bouiadjra, “Numerical analysis of the effect of pore on local stresses in spot welds,” Computational Materials Science, vol. 50, no. 8, pp. 2517–2529, 2011.