19welding technology review vol. 91 1/2019 The simulation of the stress field of welded joint containing discontinuity No 1014 mgr inż. Jerzy Nawrocki – Bilfinger-Tebodin Poland; prof. dr hab. inż. Jacek Słania – Czestochowa University of Technology Corresponding author: jerzy.nawrocki@bilfinger.com Introduction Models and their mathematical methods of description for centuries have had a very limited application due to the huge amount of algebraic operations necessary to obtain a satisfactory result. This state of affairs was changed by the invention of the semiconductor, the construction of the transistor in 1947 and the first computer (FEM) in the early 1960s. Since the mid-twentieth century, thanks to the theo- retical work of many outstanding scientist, the FEM method has been developed continuously. Their effort was realized only thanks to unprecedented progress in the field of com- puter techniques. One of the most important articles in this field from the mid-1960s was the publication of „Comput- er Experiments in Fluid Dynamic” by Harlow and Fromm. Consisting in the creation in very short time the third pillar of the science. From that time, in addition to the theory and experiment, a computer simulation is accessible. This arti- cle presents the FEM method in the field of welding. Description of the welded joint For the static simulation the model of 200 mm length butt-welded joint was used. It consists of two 40 mm thick- ness bars with accordance figure 1. For the purpose of sim- ulation the model of base material non-alloy steel S235JRG2 according to EN 10025-2 [1] and consumable GMAW wire Jerzy Nawrocki, Jacek Słania electrode G3Si1 according to EN ISO 14341 [2] was prepared. Regardless of that cracks is not allowed at any quality level EN ISO 5817 [3], the analysis was performed to figure out the stress fields around this type discontinuity. The pa- rameters of the welded joint model shown in Figures 2 and 3 Fig. 1. The welded joint Keywords: FEM Simulation; stress concentration in the welds; discontinuities; fracture mechanic DOI: http://dx.doi.org/10.26628/wtr.v91i1.999 Abstract This paper contains the results simulation and linear elastic fracture mechanic (LEFM) calculation of the welded joint contains it discontinuity 1014. The model base material C-Mn non-alloy steel S235JRG2 according to EN 10025-2 and consumable GMAW wire electrode G3Si1 according to EN ISO 14341 was prepared. Due to stress singularity effect located at the tip of the crack, the Huber-Mises-Hencky (H-M-H) hy- pothesis cannot be used, so that the stress intensity factor (KI) and fracture tough- ness (KIC) was applied. The LEFM analysis showed the range where the existence of crack can be acceptable. Based on well-known LEFM issue the verification of simula- tion program was performed. The stress intensity factor (KI) values was obtained from FEM simulation process. 20 welding technology review vol. 91 1/2019 Fig. 4. Strength of the joint according to the fracture mechanics S tr es s [M P a] Actual crack depth [cm] Fig. 3. Example of density of the grid in the discontinuities Fig. 2. The draft of the discontinuities No 1014 is as hereunder: – finite element size = 0.3 mm, – number of nodes = 156 948, – number of finite elements = 113 815. A division into tetragonal finite elements was applied, taking into account the principle of density of the grid in pla- ces of special importance. It has always been the area of occurrence of the analysed discontinuity. kr kr The nominal cross-section is equal to the total area of the cross-section reduced by the area of discontinuity: Anom=Sb-l1014 h1014 (1) where: S – width of the joint [mm], b – height of the joint [mm], l1014 – crack length [mm], h1014 – crack depth [mm]. The condition of the necessary balance is: σA=σnom Anom (2) As a consequence, the joint’s strength is: σnom=UTS (1-Anom/Sb) (3) To evaluate the stresses created by the discontinuity no. 1014 the single edge notch and criterion failure KI = KIC is adapted. The value of the fracture toughness KIC = 217 MPa•m0.5 suitable for C-Mn steel was adopted from publication [5]. The capacity of the welded joint with discontinuity can be describe as hereunder: σKI=KIC/√(Πh1014)×1/[1,12-0,23(h1014/b)+10,55(h1014/b)2- -21,72(h1014/b)3+30,39(h1014/b)4] (4) The calculation were made for the depth of the crack in the range from 5 to 75% of the joint’s thickness. The graph (Fig. 4) shows three curves: – red – tensile strength of S235 steel equal to 360 MPa, – green – static critical strength σnom, – blue – critical strength from the point of view of fracture mechanics σKI. In Figure 4, that the welded joint with crack up to 20 mm deep will break due to exceeding of the tensile strength of the base material. From the depth 20 mm the maximum strength Analysis of discontinuities using the fracture mechanics method Due to stress singularity effect located at the tip of the crack, the Huber-Mises-Hencky (H-M-H) hypothesis cannot be used. The fracture mechanics method which operates on another type approach shall be used. It is stress intensi- ty factor (KI) and fracture toughness (KIC). The finite dimen- sions of the welded joint was taken into account in this sim- ulations. The ready-made formulas for specific cracks-load cases configuration can be found in the literature [4]. kr kr 21welding technology review vol. 91 1/2019 Fig. 5. The shape of the butt joint stress field for 50% WWN – top view will decrease significantly. For such case the safe working range is below the blue line. Plastic zone of discontinuities FEM simulation results as well as experimental data indi- cates the plastic zone at the tip of the crack where the stress- es are so high, that the approaches of elasticity theory and linear-elastic fracture mechanic can not be used. However implementation of the concept small plastic area (SPA) locat- ed at the tip of the crack (Fig. 5) may solve such problems. There are, however, two conditions: firstly, the SPA area must be much smaller than the crack depth and width of the joint itself, and secondly, the load may not be too large, so that the tension generated on the tip of the crack does not af- fect the stress field too much [6]. Assuming that the metal of the joint is perfectly elastic and plastic, the actual stress- es cannot exceed the yield strength of the base material. The main problem in the correct interpretation of longitudi- nal crack simulation is therefore to determine the length of the plastic zone. According to the so-called Irwin’s first ap- proximation, the crack and its direct proximity to the stress will be approximately constant at the apex, while outside this area they will decrease as it moves away from it. If the values of the stress field exceed the yield point, the only variable will be the area of the SPA zone. The defini- tion of a conventional yield stress YS0.05 according to EN ISO 6892-1:2016-09 [7], which determines this stress as giving 0.05% strain, was used as a criterion. For the butt-joint ana- lysed, where the crack is in the native material of S235 steel with typical parameters (E = 210 GPa, v = 0.3), the elonga- tion ε for the yield point is 225/210000 = 0.00107 (0.107%). Therefore, it was assumed that the plasticisation of the na- tive material will start at the place where the simulation will obtain elongation values above 0.00107. Simulation of welded joint model For the purpose of the simulation the model of the crack whose width at the tip was 0.02 mm was prepared. It is known from the practice that the depth of the 1014 discontinuity should not exceed 2 mm, so this value has been accepted for calculations. The series of simulation was made for two load combi- nations: axial tension and torsion (Fig. 5). In order to exam- ine the dependence of the shape and surface of the plastic zones of PSA on the applied load, eight load steps from 10 to 80% YS0.2 of the base material determined as the stress utilization coefficient were used (Fig. 6). Stress intensity factor The Anysy 18.1 academic program used in this paper has been widely used since early 1970’s, therefore the verifica- tion of result quality is formality only. Nevertheless it was done to fulfilled the principles of the simulations process. For the purpose of the verification the well-known issue from fracture mechanic was used (Fig. 7). The model was divid- ed into 5160 elements and 4323 nodes using hex-dominant method with more density at the tip of the crack. The series of simulations was made at the close contour of the crack. In the context of FEM simulation, where the issue of error of discretization plays a significant role, 6 con- tours were used, after which the said parameter was calcu- lated. The contour no. 1 is closest to the face of the crack, where the grid of elements and the stress field are the most distorted by them. However, the results of the KI coefficient for contours 2 to 6 will be similar to each other. As a rule, the outermost contour is rejected because it is farthest from the front of the crack. The value of the stress intensity coefficient for the ana- lysed case can be determined in accordance with [8]: KI=σ√Πa [1,122-0,231(a/b)+10,550(a/b)2-21,710(a/b)3+ +30,382(a/b)4] (5) where: σ – tensile stress [MPa], a – crack depth [mm], b – plate width [mm], h – distance from applying the tensile force [mm]. TENSION TORSION 22 welding technology review vol. 91 1/2019 Fig. 9. KI coefficient, contour 2, WWN 80%, crack depth 10 mm DEFORMATION LINES ISOBARS OF REDUCED STRESS Fig. 7. The problem of fracture mechanics used to verify the pro- gram Fig. 8. Verification of simulation. Comparison of the values of KITEOR and KISYM coefficients Fig. 6. The shape of the butt joint stress field for 50% WWN – cross- -section ISOBARS OF TENSION STRESSES ISOBARS OF TORSION STRESSES K I [ M P a• m 0, 5 ] Tensile stress [MPa] 23welding technology review vol. 91 1/2019 Fig. 10. Distribution of the coefficient KI for contours 1 to 6 The values used for the simulation are: a = 10 mm, b = 40 mm and h = 100 mm. For such values, the above pattern is redu- ced to the formula: KI=0,265σ (6) Figure 8 shows a comparison graph for two values of stress intensity coefficients. Values calculated from formu- la 5 are designated as KITEOR, and obtained from simulation – KISYM. The analysis of the results allowed to estimate the maximum simulation error at 7.57%, with a minimum of 2.70%. The stress intensity factors KI, KII,KIII, taking into account the shape of the body and the slit, as well as the combina- tion of the load, is sufficient to fully describe the state of the stress field in front of the crack’s tip. Figure 9 shows graphi- cal simulation results of the setting with the force extending the discontinuity 1014 of the longitudinal crack in the basic material. Figure 10 shows the values of the KI factor for individu- al contours. As described earlier, the results for contour 1, the nearest to the front of the crack, significantly differ from the others. For this reason, it was not taken into account and was given only for informational purposes. The KI factor is the smallest in the middle of the crack, the largest at both ends and shows a certain symmetry with respect to the centre of the crack. This is due to the existence of a PSS condition from the not cracked material adjacent to the crack, which by stiffening both ends of the crack generates a three-axis stress condition at both ends of the crack. This coincides with a lot of experimental data, which indicates that along the entire front of the crack, in its half the tension before K I [ M P a• m 0. 5 ] The length of the gap [mm] Conclusions 1. The shape of the stress field generated by the discontinuity no 11014 significantly depends on the load combination. 2. From the fracture mechanic point of view the welded joint with discontinuity 1014 up to 20 mm depth will break in base material. In other words, the existence of discontinuity does not reduce the load capacity of the joint under static load conditions. The literature states that the stress intensity coefficient KIC applies to the crack depth value equal to ~ 0.7 of thickness of the native material. This means that for cracks more than 28 mm deep, the stress values obtained should not be taken into account. From the point of view of the fracture mechanics, there is a safe working area of the joint with a short discontinuity no. 1014 axially stretched. This is the area below the blue line in figure 4 limited to a depth of discon- tinuity of 20 mm. 3. A series of joint simulations performed in various configurations of crack’s depth-tensile stress shows that from the point of view of static loads there is a safe work area. However, it should be borne in mind that the acceptance of simulation re- sults of stress reduced directly from the tip of discontinuity leads to significantly overstated results. Therefore, according to H-M-H hypothesis, it is proposed to use stress outside the plastic deformation zone. the front of the crack is smaller than at its ends. It should be noted that in the available literature there are no calcu- lated values of the KI factor for a gap which is not crossing throughout the entire thickness of the element. References [1] PN-EN 10025-2:2007 Wyroby walcowane na gorąco ze stali konstrukcyj- nych Część 2: Warunki techniczne dostawy stali konstrukcyjnych niesto- powych. 2007. [2] PN-EN ISO 14341. Materiały dodatkowe do spawania. Druty elektrodowe i stopiwo do spawania łukowego elektrodą metalową w osłonie gazu stali niestopowych i drobnoziarnistych. Klasyfikacja. 2017. [3] PN-EN ISO 5817 Spawanie. Złącza spawane ze stali, niklu i tytanu i ich stopów. Poziomy jakości według niezgodności spawalniczych. 2009. [4] J. German, Podstawy mechaniki pękania, Politechnika Krakowska, Kra- ków, (2011). [5] D. Roylance, Introduction to Fracture Mechanic, MIT, Cambrige (2001). [6] J.W. Hutchinson, Fundamentals of the Phenomenological Theory of Non- linear Fracture Mechanics, Journal of Appiled Mechanics, ASME (1981), Tom 50. [7] PN-EN ISO 6892-1:2016-09. Metale. Próba rozciągania. Metoda badania w temperaturze pokojowej, PKN, Warszawa (2016). [8] H. Tada, P.C. Paris i G.R. Irwin, The stress analysis of cracks handbook, ASME Press, New York (2000). © 2019 by the authors. Submitted for possible open access publication under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).