Top Banner
Mechanics and Mechanical Engineering Vol. 22, No. 4 (2018) 1121–1131 c Technical University of Lodz Optimization of Welding Parameters and Microstructure and Fracture Mode Characterization of GMA Welding by Using Taguchi Method on SS304H Austenitic Steel Saadat Ali Rizvi Mechanical Engineering Department, Indian Institute of Technology (BHU), Varanasi, India e-mail: [email protected] Wajahat Ali Mechanical Engineering Department, CCS University, Meerut, UP, India Received (28 September 2017) Revised (12 March 2018) Accepted (20 November 2018) This study is centre on optimizing different welding parameters which affect the weld- ability of SS304H, Taguchi technique was employed to optimize the welding parameters and fracture mode characterization was studied. A number of experiments have been conducted. L9 orthogonal array (3x3) applied for it. Analysis of variance (ANOVA) and signal to noise ratio (SNR), a statistical technique was applied to determine the effect of different welding parameters such as welding current, wire feed speed and gas flow rate on weldability of SS304H.Tensile strength, toughness, micro hardness and mode of fracture was examined to determine weldability of SS304H and it was observed from result that welding voltage have major impact whereas gas flow rate has minour impact on ultimate tensile strength of the welded joints and optimum process parameters were found to be 23 V, 350 IPM travel speed of wire and 15 l/min gas flow rate for tensile strength and mode of fracture was ductile fracture for tensile test specimen. Keywords : GMAW/MIG, mechanical properties, ANOVA, stainless steel, mode of frac- ture, SEM weldability. 1. Introduction Metal inert gas welding or gas metal arc welding (MIG/GMAW) is an advanced version of electric arc welding in which no pressure is applied during the weld- ing process and arc is created between a continuous copper coated wire and work piece [1]. This GMAW commonly used method for joining of steels structural, com- ponents for the automotive industry [2, 3]. SS304 or 304L is the modern evolution of the noval “18-8” austenitic stainless steel. This steel is very inexpensive and versatile anti corrosion stainless steel, suitable for a broad range of general purpose use. SS304H with a higher chromium and lower carbon content. Lower carbon con- https://doi.org/10.2478/mme-2018-0088
11

Optimization of Welding Parameters and Microstructure and ...

Oct 20, 2021

Download

Documents

dariahiddleston
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: Optimization of Welding Parameters and Microstructure and ...

Mechanics and Mechanical EngineeringVol. 22, No. 4 (2018) 1121–1131c© Technical University of Lodz

Optimization of Welding Parameters and Microstructure and FractureMode Characterization of GMA Welding by Using Taguchi Method on

SS304H Austenitic Steel

Saadat Ali Rizvi

Mechanical Engineering Department, Indian Institute of Technology (BHU), Varanasi, Indiae-mail: [email protected]

Wajahat Ali

Mechanical Engineering Department, CCS University, Meerut, UP, India

Received (28 September 2017)Revised (12 March 2018)

Accepted (20 November 2018)

This study is centre on optimizing different welding parameters which affect the weld-ability of SS304H, Taguchi technique was employed to optimize the welding parametersand fracture mode characterization was studied. A number of experiments have beenconducted. L9 orthogonal array (3x3) applied for it. Analysis of variance (ANOVA) andsignal to noise ratio (SNR), a statistical technique was applied to determine the effectof different welding parameters such as welding current, wire feed speed and gas flowrate on weldability of SS304H.Tensile strength, toughness, micro hardness and mode offracture was examined to determine weldability of SS304H and it was observed fromresult that welding voltage have major impact whereas gas flow rate has minour impacton ultimate tensile strength of the welded joints and optimum process parameters werefound to be 23 V, 350 IPM travel speed of wire and 15 l/min gas flow rate for tensilestrength and mode of fracture was ductile fracture for tensile test specimen.

Keywords: GMAW/MIG, mechanical properties, ANOVA, stainless steel, mode of frac-ture, SEM weldability.

1. Introduction

Metal inert gas welding or gas metal arc welding (MIG/GMAW) is an advancedversion of electric arc welding in which no pressure is applied during the weld-ing process and arc is created between a continuous copper coated wire and workpiece [1]. This GMAW commonly used method for joining of steels structural, com-ponents for the automotive industry [2, 3]. SS304 or 304L is the modern evolutionof the noval “18-8” austenitic stainless steel. This steel is very inexpensive andversatile anti corrosion stainless steel, suitable for a broad range of general purposeuse. SS304H with a higher chromium and lower carbon content. Lower carbon con-

https://doi.org/10.2478/mme-2018-0088

Page 2: Optimization of Welding Parameters and Microstructure and ...

1122 Optimization of Welding Parameters and Microstructure and . . .

tents reduces chromium carbide precipitation during process and its susceptibilityto intergranular corrosion. SS304H frequently used in various industries such asChemical and Petrochemical, Processing industries, pressure vessels, tanks, valvesand pumps, heat exchangers, piping systems, flanges, fittings, medical, pharmaceu-tical processing, food, beverage processing and nuclear industries due to its excellenttensile strength, good weldability, and better corrosion resistance properties [4].

Dinesh Mohan arya et al. [5] investigated process parameters for Metal InertGas welding and they reported that welding current is having maximum percentagecontribution in experimental work. Nabendu Ghosh et al. [6] optimized the metalinert gas welding parameters, by Grey-Based Taguchi technique and they reportedin their result that Current having major impact in influencing the tensile strengthof welded joint as compared to gas flow rate.

Chauhan and Jadoun [7] studied the joining of two dissimilar metals SS304 andLow Carbon Steel by metal inert gas Welding (MIG) and they optimized the processparameter by using Taguchi Design Method and finally they informed that the effectof welding parameters on the tensile strength can be arranged in reducing manneras given: voltage >speed > current. Prakash et al. [8] determined the (welding)process parameters which influence the mechanical properties by using the Taguchimethod and they produced a result that Welding Current has the greatest influenceon Tensile and Hardness in the Weldability of welded joint followed by wire feedspeed and arc voltage. Bayazid et al. [9] predicted welding variables like travelspeed, rotational speed and position of plates on mechanical and microstructuralproperties of Friction Stir Welded joint of two dissimilar Aluminum alloys, i.e.,AA6063 and AA7075 with Taguchi technique and they reported that rotationalspeed, travel speed and plates position have 59, 30 and 7% influence on tensilestrength of welded joint respectively.

Saurav Datta et al. [10] developed a multi-response problem to optimize param-eters by combining to yield favourable bead geometry of submerged arc weldingbead on-plate weldment and they coupled the Taguchi optimization method withGrey relation technique to evaluate the optimal parametric combination for deeperpenetration, minimum bead height and depth HAZ of welded part.

Kalita and Barua [11] investigated the effect of the process parameters of MetalInert Gas Welding such as welding current, arc voltage and shielding gas flow rateon tensile strength of welded joints by the Taguchi optimization method and theyconcluded that welding voltage has significant effect, both on mean and variationof the tensile strength of the weld having 87.019% and 85.398% contribution re-spectively, whereas welding current has significant effect on mean only (10.807%contribution) where as shielding gas flow rate has insignificant effect on the tensilestrength of the welded joint.

Therefore in this research article an attempt has to be made to optimize theprocess parameters of metal inert gas (MIG) welding. Chikhale et al. [13] predictedthe mechanical performance of AA 6061-T6 by metal inert gas welding and theyconsider the welding current, arc voltage and wire feed speed as welding parametersand finally they optimized the parameters by reporting that welding current havingprincipal impact on the tensile strength, depth of penetration and toughness of weldjoint. Rizvi et al. [4] optimized different welding process parameters by applicationof Taguchi technique on MIG welding during bonding of IS2062 steel and they

Page 3: Optimization of Welding Parameters and Microstructure and ...

S.A. Rizivi and A. Wajahat 1123

mentioned in their outcome that welding voltage and welding current have majorimpact on tensile strength of welded joint whereas gas flow rate have least significanteffect on tensile strength of the weldment.

Liu et al. [16] investigate the tensile behavior and fracture characteristics of SSclad plate by vacuum rolling and they observed that tensile ductility increases withincreasing rolling temperature. Buddu et al. [17] investigate mechanical properties,microstructure and fracture morphology of SS304L plate welded by laser weldingand they reported in their result that tensile fracture has revealed ductile fracturemode with fine dimples and impact fracture test having lower value for weld zone,HAZ as compared to base metal. Doddamani and Kaleemulla [18] investigate thefracture toughness of Al 6061-graphite particle composite and authors told thatmaximum fracture toughness was found for Al6061- 9%Gr for a/w = 0.45 and thevalue is 16.74 MPa

√m. Gryguc et al. [19] studied tensile and fatigue behaviour of

as-forged AZ31B extrusion and it was observed that fracture was ductile fracturein forged and as extruded samples. More dimples and plastic deformation wasidentified in fractures surface in forged samples.

Patil et al. [20] conduct experiment to identify the hardness of Friction Stir weld-ing and Tungsten Inert Gas welded joints of two different Aluminium alloys, i.e.,AA7075 and AA6061 and they recorded that for two different Aluminium alloyswelded by FSW, on increasing the rotation speed and transverse speed hardnessvalue reduces and hardness is too much affected by precipitate distribution. Thevoids presence in the TIG welds contributes to reduced hardness. Authors con-cluded welded joint prepared by FSW having more hardness value as compared towelded joint prepared via TIG welding. Singh et al. [21] applied Taguchi techniqueto determine the effect of FSW on tensile properties of AA6063 under differentwelding condition and authors observed from their result that tensile strength ofAA6063 welded joint increases with increasing the rotational speed and reduceswith increasing transverse feed. Rizvi and Tewari [22] coupled Taguchi techniquewith grey relational analysis to optimize the process parameter during the weldingof SS304 by MIG welding. It was observed that wire feed speed had most significanteffect followed by voltage and gas flow rate.

2. Design and Experimental Work

2.1. Work Piece Material

In this research article, i.e., stainless steel 304H was used as raw material. SS304H plate of dimension 150 × 60 × 5 mm were bonded by using gas metal arcwelding (GMAW) machine, with polarity direct current electrode negative [DCEN].A schematic diagram of GMAW is shown in Fig. 1.

After completed weld, welded plates were machined on horizontal milling ma-chine for generating V-groove. Chemical composition of parent metal plate andfiller rod used in process for welding purpose is given in Tab. 1 respectively.

Table 1 Composition of parent metal and filler wireMaterial % C % Cr % Ni % Mn % Si % P % S Cu FeBase plate 0.06 18.68 8.54 1.9 0.41 0.031 0.005 – restFiller wire 0.08 19.5-22.0 9.0-11.0 1.0-2.5 0.3-0.65 0.030 0.03 0.75 rest

Page 4: Optimization of Welding Parameters and Microstructure and ...

1124 Optimization of Welding Parameters and Microstructure and . . .

Figure 1 Schematic diagram of GMA welding setup [15]

Taguchi design of experiment is a simple and foolproof approach which is used toreduce the analysis up to an optimal level. In this research article, three factors wereselected with their levels as shown in Tab. 2. Selection of orthogonal array was basedon DOF. Degree of freedom for all three factors is 6 in this design and welded byGMAW process with different welding parameters, nine tensile test specimens werecut from welded piece longitudinally; vertical milling machine is used to producean arc of R 12.5mm and to produced “V” notch in charpy impact test specimens.Tensile test specimens are prepared as per ASTM standard and a standard tensiletest sample is shown in Fig. 2.

Figure 2 Tensile test specimen as per ASTM [12]

All tensile test specimens are tested on UTM-40 T at room temperature. Tensiletest specimens after fracture are shown in the Fig. 3. Three different weldingparameters are used to perform the welding. In entire, this research work pure Arwas used as shielding gas to avoid any contamination of weld pool, as pure argonhaving special characteristics.

2.2. Welding Parameters

In this research article, parameters which are used for welding purpose are given inTab. 2.

Page 5: Optimization of Welding Parameters and Microstructure and ...

S.A. Rizivi and A. Wajahat 1125

Figure 3 Tensile test specimen after fracture

Table 2 Input (welding) parameters with their levelVariable Level I Level II Level III

Voltage(V) 21 22 23Gas flow rate (l/min) 10 15 20Wire feed rate(IPM) 300 350 400

In present research, an L9 OA with 3 columns and 3 rows was used. This arraycan handle three level process parameters. Nine experiments conducted to study thewelding parameters using the L9 OA. OA and the corresponding values of weldingparameters are listed in Tab. 3.

3. Results

3.1. Evaluating the Signal-to-Noise (S/N) Ratios

Noise factors are those uncontrollable factors which affect the process result (Out-put), where as derived response is known as signal. The variation of index is knownas S/N ratio. Variations are usually three types, i.e., “lower is better”, “higher isbetter” and “normal is better”. In the present experimental work the UTS, microhardness, toughness (impact strength) were output (weld quality). For good qualityof weld hardness, impact strength and UTS were considered as “higher is better”.In order to evaluate the influence of each selected factor on the responses, S/N

Table 3 Result for tensile strength hardness and impact toughnessExp. Voltage Gas wire UTS SNRA1 Micro SNRA2 Impact SNRA3No. flow feed hardness strength

rate speed1 21 10 300 540 54.6479 171 44.6599 170 44.60902 21 15 350 610 55.7066 183 45.2490 234 47.38433 21 20 400 446 52.9867 200 46.0206 210 46.44444 22 10 350 574 55.1782 233 47.3471 222 46.92715 22 15 400 571 55.1327 210 46.4444 142 43.04586 22 20 300 491 53.8216 187 45.4368 184 45.29647 23 12 400 596 55.5049 205 46.2351 216 46.68918 23 15 300 620 55.8478 176 44.9103 160 44.08249 23 20 350 569 55.1022 216 46.6891 196 45.8451

Page 6: Optimization of Welding Parameters and Microstructure and ...

1126 Optimization of Welding Parameters and Microstructure and . . .

ratios for each control factor was calculated.In this present research work tensile strength, micro hardness and impact strength

of welded pieces were acknowledged as the responses, hence, “higher the better” con-sider for ultimate tensile strength and “nominal the best” for hardness propertiesfor the analysis purpose.

S

N= −10 log

n∑i−0

1

y2i

higher is better

S

N= −10 log

1

n

n∑i−0

y2i lower is better

S

N= −10 log

1

n

n∑i−0

(yi −m)2

normal is better

For tensile strength. Response table or signal to noise is shown in Tab. 4.

Table 4 Response table for signal to noise ratiosLevel voltage gas flow rate wire feed speed

1 54.45 55.11 54.772 54.71 55.56 55.333 55.49 53.97 54.54

Delta 1.04 1.59 0.79Rank 2 1 3

Figure 4 Main effects lot for SN ratios (UTS)

3.1.1. Tensile Strength

UTS were calculated experimentally, Taguchi technique was applied for analysiswith support of ANOVA and mode of fracture was studied. On the basis of dataanalyzed, plots for S/N ratio are shown in Fig 4. It is much cleared from Fig. 4

Page 7: Optimization of Welding Parameters and Microstructure and ...

S.A. Rizivi and A. Wajahat 1127

that third level of voltage (23 V), second level of gas flow rate (15 l/min) and thirdlevel of wire feed speed (350 IPM) gives higher tensile strength.

3.1.2. Hardness

VHN of test samples also experimentally calculated, Taguchi method was appliedto find out the optimal welding process parameters with support of ANOVA. Onthe basis of data summarized, plots for the S/N ratio are expressed in Fig. 5. FromFig. 5 it is observed that and voltage (22 V), second level of gas flow rate (10 l/min)and third level of wire feed speed (350 IPM) gives normal values of hardness.

Figure 5 Main effects plot for SN ratios (VHN)

3.2. ANOVA

Analysis of variance is a statistic tool or technique, which is applied to evaluate thedifferences between the mean and their associated procedure. ANOVA result forUTS is given in Tab. 5, shows that gas flow rate has the major principal effect with59% contribution followed by arc voltage 25% contribution, while wire feed speedhaving least effect.

Table 5 Analysis of variance for SNRA1 (UTS), using adjusted SS for testsSource DF Adj. SS Adj. MS F P % contrib.

Arc voltage 2 6613.6 3306.8 11.58 0.079 25.5Gas flow rate 2 15213.6 7606.8 26.65 0.036 58.7

Wire feed speed 2 3494.2 1747.1 6.12 0.140 13.5Error 1 570.9 285.4 2.2Total 8 25892.2

ANOVA for S/N ration for hardness is summarized in Tab. 6 and it is verycleared from it that gas flow rate has the major significant impact with 58% contri-bution followed by arc voltage 25% contribution, while wire feed speed having leasteffect.

Page 8: Optimization of Welding Parameters and Microstructure and ...

1128 Optimization of Welding Parameters and Microstructure and . . .

Table 6 Analysis of variance for SNRA2 (hardness), using adjusted SS for testsSource/parameters DF Adj. SS Adj. MS F P % contrib.

Arc voltage 2 968.2 484.11 4.88 0.170 25Gas flow rate 3 310.2 155.11 1.57 0.390 58

Wire feed speed 2 1828.2 914.11 9.22 0. 0985 14Error 1 198.2 99.11 3Total 8 3304.9

4. Mode of fracture

4.1. Fracture mode of tensile test specimen

Figure 6 Fractograph morphology (SEM) of tensile test specimens

Fracture mode of tensile test specimens were studied by SEM apparatus at roomtemperature. Fracture surface of SS304H welded joint obtained with MIG weldingare shown in Fig. 6. Mode of fracture of tensile test specimens are try to understandby given figure. Fig. 6(b) shows the fractographs of tensile fractured surface. It isvery clear from figure that mode of fracture was ductile with numerous dimples. Itis also observed from result that format of fracture is not uniform; cleavage fractureis a tear type fracture. Cleavage and secondary cleavage are clearly visible in SEMimage of fracture surface.

5. Fracture mode of toughness test specimen

To determine the toughness “V” notch cahrpy samples were prepared as per ASTM.Specimens after fracture are shown in Fig. 7. Impact fractured specimens examinedto determine the surface morphology.

Page 9: Optimization of Welding Parameters and Microstructure and ...

S.A. Rizivi and A. Wajahat 1129

Figure 7 Impact test samples after fracture

Figure 8(a) and (b) shows the images of fractograph of impact chirpy “V” notchtest piece. Sample 1 in Fig. 8 shows the ductile fracture with coarse dimples.Enough finer dimples are observed in sample 3 in Fig. 8. Combined ductile andbrittle fracture mode formation is responsible for poor absorbed energy. Sample 2shows shallow dimples.

Figure 8 Macro images and SEM images of impact fracture samples for morphology

Page 10: Optimization of Welding Parameters and Microstructure and ...

1130 Optimization of Welding Parameters and Microstructure and . . .

6. Conclusions

In the present research work it is tried to investigate the effect of gas metal arcwelding (GMAW) processes on various welding variables and fracture mode chara-cterization was studied. The welding process parameters in this experiment werearc voltage, wire feed speed and gas flow rate. In this conclusion it found that:

• in this research article, the choice of the process parameters for MIG weldingof stainless steel (SS) with the optimal weldability has been marked;

• the modified Taguchi technique is applied to determine the optimal weldabilitywith three higher-the better quality characteristics;

• analysis shows that gas flow rate having significant parameter that effect theUTS and VHN followed by arc voltage and wire feed speed;

• ductile fracture mode observed with fine dimples for tensile test samples;

• SEM analysis of fractured test pieces, fracture morphology shows rough dim-ples with combination of ductile and brittle fracture in impact toughness weld-ment samples;

• all chirpy “V” notch test samples break from the centre line.

References

[1] Rizvi, S.A., Tewari, S.P., Wajahat, A.: A text book of Advanced Welding Tech-nology, Kataria & Sons (P) Ltd., New Delhi 2009.

[2] Ghazvinlo, H.R., Honarbakhsh, A., Shadfar, N.: Effect of Arc Voltage, WeldingCurrent and Welding Speed on Fatigue Life, Impact Energy and Bead Penetrationof AA6061 joints Produced by robotic MIG Welding, Indian Journal of Science andTechnology, 3(2), 2010.

[3] Ramazani, A., Mukherjee, K., Abdurakhmanov,A., Prahl, U., Schleser,M., Reisgen, U., Bleck, W.: Micromacro-characterization and modeling of me-chanical properties of gas metal arc welded (GMAW) DP600 steel, Material Science&Engineering A, 589, 1–14, 2014.

[4] http://www.aperam.com/uploads/stainlesseurope/TechnicalDataSheet/FT%20aust%C3%A9nitiques/Anglais/304%20UK.pdf

[5] Dinesh, M.A., Vedansh, C., Jyoti, V.: Parametric Optimization of MIG pro-cess parameters using Taguchi and Grey Taguchi analysis, International Journal ofResearch in Engineering and Applied Sciences, 3(6), 1–17, 2015.

[6] Ghosh, N., Kumar Pal, P., Nandi, G.: Parametric Optimization of MIG Weldingon 316L Austenitic Stainless Steel by Grey-Based Taguchi Method, Procedia Tech-nology, 25, 1038–1048, 2016.

[7] Chauhan, V., Jadoun, R.S.: Parametric Optimization of Mig Welding For StainlessSteel (SS-304) and Low Carbon Steel Using Taguchi Design Method, InternationalJournal of Recent Scientific Research, 6(2), 2662–2666, 2015.

[8] Prakash, A., Kumar, B.R., Ohdar, P., Sankar, S.R.: Parametric Optimizationof Metal Inert Gas Welding By Using Taguchi Approach, IJRET, 5(2), 176–181, 2016.

[9] Bayazida, M., Farhangia, H., Ghahramani, A.: Investigation of friction stirwelding parameters of 6063-7075 Aluminum alloys by Taguchi method, Procedia Ma-terials Science, 11, 6–11, 2015.

Page 11: Optimization of Welding Parameters and Microstructure and ...

S.A. Rizivi and A. Wajahat 1131

[10] Datta, S., Bandyopadhyay, A., Kumar, P.: Grey-based Taguchi method foroptimization of bead geometry in submerged arc bead-on-plate welding, InternationalJournal of Advanced Manufacturing Technology, 39, 1136–1143, 2008.

[11] Kalita, D., Barua, P.B.: Taguchi Optimization of MIG Welding Parameters Af-fecting Tensile Strength of C20 Welds, IJETT, 26, 43–49, 2015.

[12] ASTM E8/E8M-11. Standard Test Methods for Tension Testing of Metallic Materials.ASM International; 2013.

[13] Chikhale, S.R., Kolhe, K.P. and Kumar, P.: Prediction of Mechanical propertiesof Al Alloy 6061-T6 by using GMAW, International Journal of Current Engineeringand Technology, 5, 300–306, 2016.

[14] Rizvi, A.S., Tewari, S.P., Wajahat, A.: Application of Taguchi technique tooptimize the Process Parameters of MIG Wedging on IS2062 Steel, InternationalJournal on Emerging Trends in Mechanical & Production Engineering, 2(2), 1–11,2016.

[15] Ganjigatti, P., Pratihar, D.K. and Roy Choudhury, A.: Modeling of the MIGwelding process using statistical approaches, International Journal of Advanced Manu-facturing Technology, 35, 1166–1190, 2008.

[16] Liu, B.X., Yin, F.X., Dai, X.L., He, J.N., Fang, W., Chena, C.X., Dong,Y.C.: The tensile behaviors and fracture characteristics of stainless steel clad plateswith different interfacial status, Materials Science & Engineering A, 679, 172–182,2017.

[17] Kumar, R.B., Chauhan, N., Raole, P.M., Natu, H.: Studies on mechanicalproperties, microstructure and fracture morphology details of laser beam welded thickSS304L plates for fusion reactor applications, Fusion Engineering and Design, 95,34–43, 2015.

[18] Doddamani, S. and Kaleemulla, M.: Fracture toughness investigations of Al6061-Graphite particulate composite using compact specimens, Frattura ed Integrita Strut-turale, 41, 484–490, 2017.

[19] Gryguc, A., Shaha, S.K., Jahed, H., Wells, M., Williams, B., McKinley,J.: Tensile and fatigue behaviour of as-forged AZ31B extrusion, Frattura ed IntegritaStrutturale, 38, 251–258, 2016.

[20] Patil, C., Patil, H., Patil, H.: Experimental investigation of hardness of FSWand TIG joints of Aluminium alloys of AA7075 and AA6061, Frattura ed IntegritaStrutturale, 37, 325–332, 2016.

[21] Singh, R., Rizvi, S.A., Tewari, S.P.: Effect of Friction Stir Welding on the TensileProperties of AA6063 under Different Conditions, IJE Transactions A: Basics, 30(4),597–603, April 2017.

[22] Rizvi, S.A., Tewari, S.P.: Optimization of Gas Metal Arc welding parameters ofSS304 Austenatic Steel by Taguchi–Grey relational analysis, Journal of Computa-tional and Applied Research in Mechanical Engineering, accepted for publishing in 8,2018.