Microsoft Word - numero_50_art_21_2580 L. Romanin et al., Frattura ed Integrità Strutturale, 50 (2019) 251-263; DOI: 10.3221/IGF-ESIS.50.21 251 Focussed on IGF25 – Fracture and Structural Integrity International Conference 2019 A metallurgical and thermal analysis of Inconel 625 electron-beam welded joints L. Romanin, P. Ferro, A. Fabrizi University of Padova, Italy luca.romanin@phd.unipd.it paolo.ferro@unipd.it, https://orcid.org/0000-0001-8682-3486 alberto.fabrizi@unipd.it F. Berto University of Trondheim (NTNU), Norway filippo.berto@ntnu.no, https://orcid.org/0000-0002-4207-0109 ABSTRACT. Inconel 625 is a nickel-based superalloy that finds application in many industrial sectors thanks to its high strength, excellent fabricability, and outstanding corrosion resistance. It is characterized by a good weldability and often used in the as-welded conditions. In the following, a metallurgical and thermal analysis of Inconel 625 Electron- Beam welded joints is described. The thermal analysis was supported by a numerical model that uses a superimposition of a spherical and a conical shape heat source with Gaussian power density distribution in order to reproduce the nail-shape of the fusion zone. The heat source parameters were calibrated by using experimental data coming from metallurgical observations and temperature measurements. Numerical and experimental results were found in good agreement. KEYWORDS. Welding; Inconel 625; Thermal Analysis; Numerical Modelling; Microstructure. Citation: Romanin, L., Ferro, P., Fabrizi, A., Berto, F., A metallurgical and thermal analysis of Inconel 625 electron-beam welded joints, Frattura ed Integrità Strutturale, 50 (2019) 251-263. Received: 27.07.2019 Accepted: 17.08.2019 Published: 01.10.2019 Copyright: © 2019 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. INTRODUCTION ickel-base superalloys for aerospace and industrial applications are of great interest in the materials research field. Their creep and oxidation resistance properties make them the base material for aero-engine components in hot sections, such as turbine discs, blades and casings. In such applications, the possibility to repair or assemble different parts by welding is a fundamental task for the design and industrialization of these products. Therefore, studies related to the weldability of this kind of material are of primary importance in order to obtain sound welds. Their hot cracking sensitivity, arising from a synergy of phenomena such as wide solidification temperature interval, chemical N http://www.gruppofrattura.it/VA/50/2580.mp4 L. Romanin et al., Frattura ed Integrità Strutturale, 50 (2019) 251-263; DOI: 10.3221/IGF-ESIS.50.21 252 composition and thermal stresses, is the most important problem to solve since it enormously compromises the low cycle fatigue resistance [1]. For this reason, high power density welding technologies [2–4], with a particular attention to electron-beam welding (EBW), are the most favorite. In fact, they ensure a very narrow fusion zone (FZ) as well as heat- affected zone (HAZ) and a reduction in distortion and residual stresses. In particular, the high vacuum, which characterizes the EBW, avoids material contamination and yields deep penetration. Experimental results carried out on a limited number of nickel-base alloys [4–6], showed that it is possible to obtain sound welds by optimizing the welding speed, pre-heating, beam focus and welding current. Thanks to computer development, it is possible to avoid expensive experimental work and use numerical simulation for predicting the stress evolution during welding [7,8]. A typical keyhole weld is usually accompanied by a distinct ‘nail head’ appearance at the top. The first attempts to model this shape were analytical. Klykow et al. [9] combined a depth source, which describes the keyhole effects and the energy related transfer, and a surface source, representing the plume radiation on the workpiece surface. Both sources had a Gaussian power distribution. Here, it was found that, in laser welding, plasma may be regarded as an independent heat source in the surface, and the laser beam as a volume source resulting in a dagger shape form of the penetration zone. Steen et al. [10] combined the point and line source of Rosenthal [11] to model more effectively a keyhole weld and estimate the power actually absorbed by the weld. The line source represents absorption down the keyhole and the point source represents the plasma radiation from the plume. Bonollo et al. [12] studied the effect of the operative pressure on the weld beads of different materials welded using a CO2 laser beam. An analytical model of the process was developed using two heat sources that allowed the quantitative evaluation of the power distribution between the keyhole and the plume. Experimental and analytical results were in good agreement. Instead of using two sources, Binda et al. [13] proposed a semi-empirical model of the temperature field in laser welding, based on a modification of the Rosenthal solution. The hypothesis of a constant heat rate along the thickness was replaced by a generic unknown function, to be determined from experimental data. All the analytical models do not consider the variation of the thermal properties with temperature and phase change. Besides, they are not able to predict the thermal and residual stress fields induced by the welding process. For this reason, numerical thermo–mechanical models were developed. In high power density welding technology, a conical source shape with Gaussian power density distribution is often used [7,14]. Though the experimental and numerical results are in sufficient overall agreement, a single source is insufficient to reproduce the correct ‘nail’ shape of the fusion zone and thus the local stress development at the FZ–HAZ interface. Du et al. [15], developed a mathematical model for flow simulation of full penetration laser beam welding of titanium alloy. The model is based on a plane heat source on the top surface and a cylindrical heat source along the z-direction, which take into account the plasma effects and the keyhole absorption. The model, based on the numerical solution of the fluid mechanics equations, gave a good prediction of the FZ shape but was not used for the stress evaluation arising from the welding process. Nickel-base superalloys suffer, above all, from weld-cracking tendency [5]. The idea to correlate weld cracking with the stress–strain evolution during weld cooling was developed by Feng et al. [16]. They considered only transversal and longitudinal stresses as a function of weld cooling. Cracking will be promoted if a weak microstructure and/or a sufficiently high tensile stress exists. However, the model does not consider the correct FZ shape of the bead and the 3D distribution of the stress field. Dye et al. [17] proposed a numerical method for the prediction of the processing conditions that are liable for producing defects during welding such as constitutional liquation, solidification cracking and a centreline grain formation. In particular, solidification cracking is assumed to arise due to the generation of a positive transverse stress at the point behind the heat source where the liquid fraction is still significant. The model was applied to TIG welding of Inconel 718 using a semi-analytical solution for the temperature field and a thermal elasto-plastic analysis in two-dimension, assuming a state of plane stress. In such a way, the real shape of the FZ through the thickness and thus the stress component in z-direction was neglected. For this reason, the model is not suitable for the prediction of microfissures formation along the thickness. It is noted that a lack of knowledge is still present in the specific literature on EBW process applied to alloy systems such as Nickel-base superalloys. Mayor [18] investigated the characteristics of Inconel 718 and 706 joints, though welded by GTAW (gas tungsten-arc welding), finding an excellent weldability and good tensile strength at both room and elevated temperature. The thermal effects of welding operations have been observed to affect the desirable structure of the heat- treatable alloys by producing heat-affected zones with poor stress-rupture properties. The properties are substantially restored to the required level, however, by post-weld heat-treatment involving, in general, solution and ageing treatments at conventional temperatures. In this case, age hardening treatment following welding was found to be more than adequate, without a need for a re-solution treatment. Ferro et al. [19] investigated the electron beam welding of Inconel 706. They found a good agreement between thermal-mechanical analysis and experimental data. The numerical model was L. Romanin et al., Frattura ed Integrità Strutturale, 50 (2019) 251-263; DOI: 10.3221/IGF-ESIS.50.21 253 helpful in correlating microfissures occurring at the grain boundaries, under the nail head of the bead, with process parameters. Laki et al. [20] studied the EBW of 30HGSA steel using a conical heat source. They also found good consistency with empirical data. In a more recent investigation, Laki et al. [21], from a series of 49 weld experiments, created a Partial Least Square model of the FZ and subsequently defined the heat source utilizing mesh segments and constant power density on each segment. The phenomenological approach has been used by Palmer et al. [22] when they simulated the EBW of 304L stainless steel. They included thermal conductivity and viscosity to account for enhanced heat and mass transfer due to turbulence in the weld pool and found that convective heat transfer was very significant in determining the weld geometry. However, this type of analysis is not convenient for a subsequent mechanical analysis because it involves results mapping from a Finite Volume Method to a Finite Element Method code. It is pointed out that only few researchers have investigated in the past the effects of high-power density welding techniques on Ni-base wrought superalloys. Some works can be found about laser welds in wrought Waspaloy [2] but specific information on the high-energy beam welding of other Ni-based superalloy have not yet fully established. In the recent literature some works focused their attention to the precipitation of Laves phase particles [23] or equilibrium and metastable phase/grain morphology [24] during laser additive manufacturing or rapid solidification of Ni-based superalloys. Such work uses quite complex models solved by the phase-field method or multi-scale model combining the finite element method and stochastic analysis. The JMAK-type model was instead applied by different researcher [25–27] with the aim to predict the microstructure evolution in nickel-based superalloys. In particular the Avrami exponent for ’ precipitation was determined to be in the 1.5-2.3 range, suggesting spherical or irregular growth. A methodology was also developed to take into account the temperature dependence of the rate coefficient k(T) in the non-isothermal JMA equation. Finally, it is interesting to mention a recent review paper about welding characteristics of aerospace materials that highlight the potentiality of high-power density welding processes (LW, EBW) compare to the traditional arc welding technologies [28]. In this work, a numerical model of EBW of Inconel 625 has been developed. The phenomenological approach proposed is based on the calibration of the source parameters on the basis of macrographs and temperature measurements. Experimental and numerical results were found in good agreement and suitable for a future mechanical analysis. The obtained results will be used also for a further development of the thermal-metallurgical model that includes phase precipitation and dissolution [29]. Moreover, the utilized welding process did not produce any grain growth or detrimental phases such as σ, μ and Laves Phases in more than trace amounts as it will be described afterwards. Hardness measurement has been presented, as well. MATERIAL AND EXPERIMENTAL PROCEDURE he chemical composition of the specimens has been measured by Energy Dispersive X-ray spectroscopy (EDS) and is summarized in Tab. 1. Ni Cr Fe Mo Nb+Ta C Mn Si Al Ti Co 58.30 20.39 3.94 8.91 3.76 - - 0.32 - 0.22 - Table 1. Chemical composition calculated by EDS on specimen Test 1 (wt%) Four butt-welded joints were produced, each one obtained by electron beam welding two plates which dimensions are 16 mm x 72 mm x 2.5 mm (Fig. 1). Since no filler material has been used, the joining edges have been prepared to obtain the best contact. In order to validate thermal analysis results, two K-type thermocouples have been inserted at 2.5 mm and 6 mm far from the weld seam, as shown in Fig. 1. The holes, with a diameter of 1 mm, have been realized with Electro Discharge Machining. Picotech TC-08 Termocouple Data Logger has been used for data recording. The plates to be welded were first clamped as shown in Fig. 2 and then tack welded in six points, one every 10 mm, in order to avoid plates detachment due to thermal expansion. Weld tacking increased also the overall specimen temperature because of the small mass of the specimens. Welding tests have been carried out according to the parameters defined in Tab. 2 and using a vacuum chamber with a pressure of 4 10-4 mbar. Full penetration was not reached in the first 10 mm of the welding line. However, starting and ending points are not taken into account in this work. In Tests 2 and 3, the current proved to be too high since a crater formed at about 2/3 of the T L. Romanin et al., Frattura ed Integrità Strutturale, 50 (2019) 251-263; DOI: 10.3221/IGF-ESIS.50.21 254 welding line. Because starting and ending lack of full penetration could be neglected, the current has been therefore slightly lowered to 14.5 mA in the last welding test (Test 4). Optical and Scanning Electron Microscopes (with microprobe EDS) were used for the metallurgical investigation. For optical microscope (OM) analysis, the specimens have been etched with a solution made of nitric acid (HNO3) and chloridic acid (HCl) (1/3 parts in volume). Precipitates were analyzed in terms of shape, dimension and chemical composition. Detected X-ray spectra have been compared with previous ones published in literature with the aim to find matching patterns and recognize phases with a reasonable degree of accuracy. Micro-segregation patterns across the dendritic substructure of the FZ were finally checked. Figure 1: 3D representation (a) and top view (b) of the specimen showing the holes position for the thermocouples. Figure 2: Clamping configuration and EBW in progress. Accelerating voltage (kV) Beam current (mA) Focus current (mA) Spot diam. (mm) Freq. (Hz) Speed (mm/s) Notes Test 1 150 15 347 1.5 300 12 - Test 2 150 15 347 1.5 300 12 - Test 3 150 15 347 1.5 300 12 - Test 4 150 14,5 347 1.5 300 12 Thermocouples Table 2: EBW tests parameters EXPERIMENTAL RESULTS Metallurgical results The FZ shape is shown in Fig. 3a. The grains are dendritic and oriented along the heat flux direction (Fig. 3b). Fig. 4 reveals both the characteristic epitaxial grains grow at the FZ/HAZ interface and a narrow HAZ of about 0.2 mm. Grains, in the parent material, maintains an equiassic morphology. No grain growth was detected in the HAZ, as it could be observed in Fig. 4. This is because EBW is characterized by a rapid heating and cooling, greatly reducing the material L. Romanin et al., Frattura ed Integrità Strutturale, 50 (2019) 251-263; DOI: 10.3221/IGF-ESIS.50.21 255 alteration compared to arc welding. It is also noted that in the heat affected zone, grains boundaries are less evident than those in the parent material (PM). Figure 3: FZ shape (a) and microstructure (b) Figure 4: HAZ details Grain boundaries in the PM were carefully analyzed by SEM (Fig. 5). A dark precipitate (Label 1), of about 3 μm, and a series of lighter and smaller precipitates (Label 2) with dimensions in the order of 0.1 μm were detected (Fig. 5b). The dark precipitate contains mostly Titanium and Niobium (Fig. 5c) and is therefore classified, according to literature [30], as MC (TiC) carbide. The smaller precipitates (Label 2) are distributed along the grain boundaries and have a fine globular aspect. Moreover, they contain mostly molybdenum and chromium like the M23C6 carbides in references [31,32]. Niobium carbides and Leave’s precipitates are also found in the PM as shown in Fig. 6 [30,33]. No other type of precipitate has been found in the parent material. By comparing the spectra coming from EDS measurements in the dendritic and interdendritic zones, Mo was found to segregate in the interdendritic space (label 10) (Fig. 7). For the sake of simplicity, spectra coming from points 8 and 9 of Fig. 7 are not reported in this work. Anyway, because of their content of nickel, molybdenum, niobium and chromium, they were classified as Leave’s phase [33,34]. L. Romanin et al., Frattura ed Integrità Strutturale, 50 (2019) 251-263; DOI: 10.3221/IGF-ESIS.50.21 256 Figure 5: Grain boundaries precipitates at different magnifications (10000x (a), 40000x (b)) and EDS spectra in points 1 (c) and 2 (d) Figure 6: Precipitates in PM: NbC (a), Leave’s phase (b) Microhardness profiles Two Vickers microhardness profiles across the FZ are shown in Fig. 8. The FZ extends from about -1500 μm to +1700 μm. In that region, hardness is more variable reflecting the inhomogeneous dendritic microstructure. However, microhardness results are contained around a mean value of about 243 HV. This confirms both the very narrow HAZ and the good electron beam weldability of the analyzed alloy. As a matter of fact, no microhardness profile depressions were observed that could suggest a sensitive reduction of welded joints mechanical properties compared to those of the parent metal. L. Romanin et al., Frattura ed Integrità Strutturale, 50 (2019) 251-263; DOI: 10.3221/IGF-ESIS.50.21 257 Figure 7: Precipitates and Mo segregation in FZ. Figure 8: Microhardness profiles NUMERICAL MODEL non-linear thermal transient analysis of the electron beam welding process has been performed. The fast motion and the extremely high power concentration of the heat source cause severe temperature gradients and heating/cooling rates which have to be captured by a fine mesh and a small time-step. Taking advantage of geometrical and load symmetry, only one plate has been modelled. A time step of 0.05 s was set during the passage of the heat source over the piece. A L. Romanin et al., Frattura ed Integrità Strutturale, 50 (2019) 251-263; DOI: 10.3221/IGF-ESIS.50.21 258 The thermal model has been validated using data of the specimen named “Test 4” in which two thermocouples have been inserted and a macrograph of the middle cross section has been obtained (Fig. 3a). About 40000 3D linear elements were used. The mesh has been graded both in the transversal and in the longitudinal direction to optimize the number of elements. On the symmetry plane, along the thickness, 8 elements of around 0.3 mm have been used (Fig. 9, cross section). Moving far from the weld bead, gradually, only 4 elements along the thickness have been used. In the longitudinal direction, the mesh is finer at the start and end of the weld seam and coarser in the middle of the plate where thermal condition can be assumed quasi-stationary (Fig. 9, top view). For this reason, in the longitudinal direction, a parabolic bias has been imposed to have a finer mesh in the starting and ending regions because of the higher temperature gradients due to the start and stop of the weld. The elements size ranges from 0.17 mm (welding start/end) to 0.77 mm in the middle of the plate. Figure 9: Mesh of the numerical model Boundary Conditions Heat loss by radiation and conduction trough the clamps has been taken into account. Thermal radiation was simulated by using the Stefan-Boltzman law, while the heat loss by conduction by using an equivalent convection coefficient. As a matter of fact, because of the small dimensions of the specimen, the effect of the clamps is non negligible and the advantage of using the convection coefficient is to avoid the modelling of the contact resistance which would have required a phenomenological calibration and introduced another source of error. The value of the convection coefficient was empirically determined to be 10 W/m2, while the room temperature has been set to 25°C. The initial temperature has been set to 95°C, according to experimental data. This is due to the tack welding operation carried out before EBW. Heat Source During EBW, in the beam impingement zone, a portion of material evaporates creating a cavity containing metallic plasma because of the high temperature reached. L. Romanin et al., Frattura ed Integrità Strutturale, 50 (2019) 251-263; DOI: 10.3221/IGF-ESIS.50.21 259 Various authors have studied the keyhole formation and evolution in EBW or laser welding by using CFD analysis with the aim to predict gas porosity. In this case, a phenomenological approach has been adopted that uses thermocouples data and micrographs of the FZ to calibrate a power density distribution function. In particular, the characteristic nail-shape of the FZ has been simulated by superimposing a spherical and a conical heat source. The spherical heat source is defined by a constant power density distribution (q0) in the range between r = 0 and r = R1 (inner radius) and a linear power density decreasing from R1 to R2 (outer radius) (Eq. 1). 0 1 2 0 1 2 2 1 0q if r R q R r q if R r R R R        (1) The conical heat source, on the other hand, is characterized by a Gaussian power density distribution centered in its axis (Eq. 2) 2 2 0 0 r rq q e   (2) In Eq. (2) r2=(z-vt)2+x2 and r0 =Re-(Re-Ri)(ye-y)/(ye-yi) for yi> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice