Multiple Response Optimization of Dimensional Accuracy of Nimonic Alloy Miniature Gear Machined on Wire Edm Using Entropy Topsis Andpareto Anova

2021-04-29 07:50TinaChaudharyArshadNoorSiddiqueeArindamKumarChandaShafiAhmadIrfanAnjumBadruddinandZahidKhan

Tina Chaudhary,Arshad Noor Siddiquee,Arindam Kumar Chanda,Shafi Ahmad,Irfan Anjum Badruddinand Zahid A.Khan

1Department of Mechanical and Automation Engineering,Indira Gandhi Delhi Technical University for Women,Kashmere Gate,110006,India

2Department of Mechanical Engineering,Jamia Millia Islamia(A Central University),New Delhi,110025,India

3Department of Mechanical and Automation Engineering,G.B.Pant Engineering College,New Delhi,110020,India

4Mechanical Engineering Department,College of Engineering,King Khalid University,Abha-61421,Asir,Saudi Arabia

ABSTRACT The purpose of this research is to obtain the optimum cutting parameters to achieve the dimensional accuracy of Nimonic alloy miniature gear manufactured using Wire EDM.The cutting parameters investigated in this study are current,pulse on time(PON),pulse offtime(POFF),wire tension(WT)and dielectric fluids.Ethylene glycol, nanopowder of alumina and oxygen are mixed to demineralized water to prepare novel dielectric fluids.Deviation in inner diameter, deviation in outer diameter, deviation in land and deviation in tooth width are considered to check the dimensional accuracy.Taguchi L16 is employed for experimental design and multiple response optimization is performed using Entropy TOPSIS and Pareto ANOVA.Results indicate that pulse on time is the most notable parameter for good dimensional accuracy followed by dielectric fluid,current,pulse offtime and wire tension.Ethylene glycol mixed demineralized water is preferred for low dimensional deviation.The optimum WEDM parameters are pulse on time at 20 µs,Ethylene glycol mixed demineralized water dielectric fluid,current at 3 A,pulse offtime at 4µs,and wire tension at 18 N.

KEYWORDS Miniature gear; multiple optimization; dimensional accuracy; entropy TOPSIS; Pareto ANOVA

1 Introduction

Among various micro parts, micro gear is the extensively used components in MEMS and NEMS technology, watches, turbines, pumps, harmonic drives, dental and medical devices, micromotors, precision measuring instruments, and electronic home appliances, etc.[1].Use of the miniature gears is continuously rising.Gears of less than 10 mm diameter is considered as miniature gears and subdivided into meso gear with outer diameter between 1-10 mm and micro gear with outer diameter less than 1 mm [2].Fabrication of miniature gears were initially done by traditional processes like hobbing, gear rolling, stamping, die casting, and powder metallurgy etc.There are several limitations of traditional processes such as dimensional inaccuracy, tool wear, requirement of post finishing operation, and more time for setup etc.[2-4].Non-traditional methods such as electric discharge machining (EDM), Wire EDM, wire electric grinding, etc.are in trend now a days.These nonconventional processes are preferred due to high accuracy,repeatability, good surface integrity, less time for setup and less stresses, etc.[5-7].In WEDM,the material removal rate (MRR) is driven by electro-thermal mechanism.Wire feeds from upper spool using a wire guide; the spark formation takes place between wire and work due to high thermal flux density as the electrode wire crosses near workpiece, which results material removal through erosion and evaporation.Due to its distinct characteristics enable is suitable to use for high end applications in automotive and aerospace industries [8].Wire-EDM is the variant of EDM and is most suited for the miniature gear manufacturing.It has ability to machines electrically conductive material with ease irrespective of its high hardness and toughness which make is difficult to machine via conventional methods.Furthermore, gears fabricated by Wire-EDM process are marked with good surface finishing, high accuracy, debris free surfaces and as well cost effective in comparison to other processes [9].

Mainly, miniature products including mini gears are made from a variety of materials like aluminium, bronze, brass, nickel, steels, and titanium based alloys, etc.[10].Wide varieties of materials have been machined by WEDM and machining performance has been evaluated and reported.Many difficult to cut materials like nickel, in conel alloys, tungsten and titanium can be machined by WEDM machining [11,12].Experiment based investigation on micromachining of titanium alloys have been conducted widely [13-16].Titanium alloys show excellent properties corrosion resistance and high specific strength, etc.[17].Nimonic alloys show excellent high temperature properties but are difficult to machines and the work on micro manufacturing are less reported [18-20].

Dimensional accuracy (DA) plays crucial role in deciding the quality of die manufactured using Wire EDM [21].During machining by Wire EDM, material removal mechanism is based on melting and vaporization from micro-sized local affected volume on the workpiece surface.Thermal energy during WEDM is evolved due to spark generation and thermally affects upper layers of machined surface.These layers have different properties from original material [22].Therefore, DA is a major issue in WEDM which needs considerable research attention.Good DA enhances the fatigue strength, wear resistance, and corrosion resistance [23-26].The wire-EDM parameters, like cutting speed, current, pulse on time (PON) and pulse off time (POFF), wire tension (WT) on DA have been studied [27-32].Another issue during WEDM is overcut due to spark generation from the lateral surface and bottom surface of the wire electrode.This causes concerns of dimensional inaccuracy of WEDM manufactured component.Peak current and PON time are identified as major parameters to affect the dimensional overcut [33-35].

It is important to note that more than one output responses are always involved when defining the quality of surfaces machined by WEDM.At times the effects of a single parameter on different output response are often contracting.This makes the multiple response optimization important in WEDM research.In recent researches, fuzzy logic multiple response decision making technique has become popular to optimize various manufacturing processes.Fuzzy-TOPSIS was applied to get the optimum responses like machining time, electrode wire wear rate, and dielectric fluid consumption during WEDM machining of steel [36].Grey-fuzzy logic oriented technique was performed to optimize MRR, and SR while machining Al-4.5%Cu-TiC MMCs [37].Principal component analysis (PCA) and fuzzy with Taguchi technique was applied to get the optimum values of WEDM parameters [38,39].Combined approach of TOPSIS and fuzzy TOPSIS was conducted to optimize the multiple responses characteristics of the turning operation on the pure titanium grade 2 super alloys [40,41].Taguchi and TOPSIS combined approach was applied to check the impact of EDM process parameters on SR of the machined components [42].The importance of multi-response optimization and selections of new yet important process parameters which are not extensively investigated and affect the responses will prove useful to micromachining of gears.Based on the literature survey, it has been observed that TOPSIS methodology may prove to be effective to perform optimization on multi responses of WEDM process by using important process parameters and by addition dielectric fluid as one of them [43].

Dielectric fluid helps to remove material and debris from surface after machining.It acts as an insulating medium between the wire and workpiece and generates required energy.Mostly traditional dielectric fluid such as Deionized water, kerosene, white spirit has been used as a dielectric but researchers have identified new additive mixed dielectric fluids as well.Mixture of Aluminum and SiC in Kerosene improves the surface finishing by enlarging the space between wire electrode and workpiece which scattered the spark energy.Researchers have demonstrated that the mixture of aluminium powder with kerosene is better than mixture of silicon carbide with kerosene [44].Graphite mixed dielectric reduces the surface tension of itself which results in improving the MRR.The study also revealed that the graphite percentage in dielectric directly affected SR and recast layer thickness.The high graphite percentage in dielectric fluid increased the recast layer thickness [16].Oxygen has also been mixed with dielectric to enhance the sparking by oxidation process and improve the MRR, DA and cutting speed [45].Oxygen gas improves cutting speed due to reduction in the inner electrode gap [46].Nitrogen has also been used as a dielectric medium but its cutting speed is lower comparing with other gaseous medium.Air is a good dielectric medium with high cutting speed because of presence of oxygen [47].The studies on WEDM have been reported in which investigations on specific hybrid dielectric have been performed but the use of dielectric as a process parameter is not used.Further, the constituents in powder, gas and liquid form have distinct effect but a hybridization of each of these is done separately.When the constituent in liquid/gas/powder form are mixed and hybrid dielectric is prepared a synergy among constituent may be advantageous.Further, by making different dielectric fluids as process parameter and analyzing their effect on responses simultaneously will prove greatly useful.

In the maiden present work hybridization of dielectric by mixing additive such as metals and ceramic powder, gas and liquid form in basic dielectric “DM water” have not been investigated.The hybridization by mixing additives in all three physical phases solid, liquid, gas (such as alumina ceramic powder, ethylene glycol and oxygen) in the basic dielectric DM water has been carried out.A hybridization of such as kind is also one of the novelty characteristics of this work.

The purpose of the study is to find out the optimum cutting parameters to achieve improved DA.Deviation in inner diameter (DID), Deviation in outer diameter (DOD), Deviation in land(DLAND), and Deviation in tooth width (DTW) are the parameters to investigate the DA.The scheme for the final experiments was chosen and the final experiments were performed to analyze the influence of process parameters and different dielectric fluids on the DA.Entropy TOPSIS and Pareto ANOVA were applied for multiple optimization of cutting parameters.

2 Experiments

2.1 Materials and Methods

Nimonic alloy is used to manufacture miniature gear of inner diameter of 5 mm, outer diameter of 7 mm and land of 0.5 mm.At present, Nimonic alloy is extensively used in many fields such as aerospace, transportation, and medical, etc.Main reason for extensive use of Nimonic alloys is its high temperature properties like wear resistance, strength and hardness, however it is difficult to cut.Further, little literature is available machining on Nimonic alloy for miniature gear.The chemical constituents of Nimonic alloy is shown in Tab.1.

Table 1:Chemical composition of Nimonic alloy (wt.%)

Wire EDM machine DK7712 CNC of Steer Corporation was used in experimentation.Schematic diagram of Wire EDM is represented in Fig.1.The specifications of the Wire EDM employed in the work are given in Tab.2.Molybdenum wire of 0.18 mm diameter was employed as electrode due to its high tensile strength and good conductivity.Special arrangements to provide measured variation to the wire tension were made by installing a push pull gauge which produced tension in electrode wire and were used as one of the process parameters.

The dielectric fluid is one of the most important factors in WEDM which affects all the response factors.Usually WEDM employs demineralized water (DM) as dielectric.In the present work the dielectric fluid was used as a base-dielectric.Constituents such as, Ethylene Glycol (EG,an organic liquid) Al2O3(alumina nanopowder) and oxygen gas were blended in the base-dielectric and four different blends of dielectric fluids were prepared as per the details given in Tab.3.The constituents and their concentrations were selected based on literature, dielectric strengths and trial experiments.The process parameters “dielectric fluid” was considered as a process parameter and four blends of dielectric fluid were used as its four levels.

All the hybrid dielectric fluids hold different dielectric strength and play significant role in spark gap and discharge energy.The dielectric strengths of DM water, oxygen, alumina nano powder and ethylene glycol are 70 kV/mm, 0.92 kV/mm, 15 kV/mm and 35 kV/mm respectively.Dielectric fluid which has high strength delivers high spark energy and the one with low strength deliver slow spark energy.After machining dimensions of the machined miniature gear were measured using optical projectors (diameter 300 mm) of Banbros, JT series as shown in Fig.2.The parameters and their values were chosen on the basis of literature and exhaustive trial experiments and they are given in Tab.3.

The WEDM process works melting evaporation and erosion by electro-thermal energy which removes the material by distinct sparks between the workpiece and electrode wire.The sparking takes place with a present periodicity controlled by spark on and spark off time.During spark on the dielectric breaks down and releases the energy resulting in material removal.During spark off the dielectric fluid flushes out the debris from the gap and also cools down the workpiece.Dielectric is an insulating medium which avoids electrolysis effect on the wire electrode during machining.For every pulse, discharge takes place at single point where wire electrode material releases the energy which melts and vaporizes the material.Thus, small crater forms on electrode and workpiece surface [48].The eroded particles form the debris which are removed by dielectric flow.

Figure 1:Schematic of wire EDM system

2.2 Entropy TOPSIS Method

There are various methods to access the weight of decision making processes like eigenvectors method, entropy method, weighted least square method, and linear programming method for multiple dimensional analyses of preferences (LINMAP) [49].However, LINMAP and entropy methods are highly recommended when the decision matrix data is known.Entropy methodology provides faster results compared to LINMAP [50].The entropy technique is mostly preferable to analyze the discrepancy between the set of information.In the entropy methodology, if identical values are obtained for some attributes for the dissimilar alternatives, attributes has to be removed.

Table 2:Wire EDM specifications

Table 3:Range and levels of machining parameters

The basic approach of TOPSIS methodology is the perfect alternative having minimum distance from positive ideal solutions (PIS) and maximum distance from negative ideal solutions(NIS) [51].The benefit of TOPSIS technique is that it indicates the rational of human choice including both good and bad attributes of the alternatives concurrently, and the clarity on presentation and computation [52].Number of steps is not affected by number of attributes thus it provides a quick solution [53].In present scenario, TOPSIS method is used successfully as a decision-making tool into different areas such as water management [54-56], human resource [52], transportation planning [50], mechanical engineering [57], policies development [56],and manufacturing engineering [51].

Figure 2:Optical profile projector

Entropy TOPSIS method is used in this work to obtain the optimum cutting parameters for improved DA.Since, DA is defined in terms of DID, DOD, DLAND and DTW.Entropy TOPSIS method transforms multiple DA measures viz.DID, DOD, DLAND and DTW into a single parameter Pi.Entropy method is used to determine the weights of the responses.Whereas, TOPSIS method is employed to compute the single replicated parameter Pi for DA.The computational steps involved in Entropy TOPSIS method are as follows:

Step 1:Formulate a decision matrix by arranging the alternatives in rows and attributes in columns as shown in Eq.(1).For better understanding of the computations steps,number of alternatives and attributes are denoted as “n” and “J” respectively.

Step 2:Normalize the decision matrix using Eq.(2):

Step 3:Compute entropy of the attributes using Eq.(3).

Step 4:Compute weight of the attributes using Eq.(4).

where,dj=1-Ej

Step 5:Calculate weighted normalized decision matrix as per Eq.(5).

where,wjis weight ofjth attribute and

Step 6:Determine Positive and negative ideal solutions by using Eqs.(6) and (7):

Step 7:Calculate separation measures by using Eqs.(8) and (9).Ideal solution for separation of alternative is given as:

Step 8:Calculate relative closeness from the ideal solutions as per Eq.(10).

Step 9:Ranking the order on the basis of descending order ofis further defined as an index which indicates the combined performance values of the alternatives.It is to be noted that for better performance of the attributes, higher values ofis required.

2.3 Pareto ANOVA Optimization Technique

Pareto ANOVA is mostly used to examine the data for process optimization [58].Pareto method is depended on Pareto principles which explain that 80% output is obtained by 20%input and it is comparatively easier than the statistical or classical ANOVA developed by Fisher.Therefore, ANOVA table and F-ratio are not required in this technique.Hence, this method is recommended when there is lack of sufficient knowledge of statistical/classical ANOVA approach.The advantage of this technique is that it provides simplified technique to analyze the results and this is the reason it is highly useful for the industrial practitioners.

3 Results and Discussions

The experiments were accomplished using L16orthogonal array (OA) of the Taguchi’s design of experiment.Three replicates of each experiment were performed and average values of DID,DOD, DLAND, and DTW were recorded.The experimental data for deviation in inner diameter,outer diameter, and tooth width are shown in Tab.4.One of the miniature gears which were machined is shown in Fig.3.

Table 4:Experimental data for deviation in inner diameter, outer diameter, land and tooth width

Deviations in the inner diameter, land, outer diameter and tooth width are the factors that were taken into account to define the dimensional accuracy of the miniature gear.A single multiresponse comprising of all key elements of miniature-gear geometry are investigated in this maiden work.Hence, Entropy-TOPSIS multiple response characteristics was applied to assess the deviation of gear and Pareto ANOVA analysis was conducted to obtain the optimum combination of parameters shown from Tab.5 to Tab.11.

Normalized decision matrix andrijnormalized values were determined.Decision matrix was converted into dimensionless matrix by usingrijmatrix.rijmatrix for multiple responses is shown in Tab.5.

Figure 3:Nimonic alloy miniature gear

Table 5: rij matrix for DA

Weighted dimensionless matrix was obtained by the entropy methodology.Response variable weight for multiple responses is shown in Tab.6.

Then positive and negative ideal solutions were identified.Values of the distances from PIS and NIS have been obtained as shown in Tab.7.Closeness coefficient (Pi) was calculated after calculating PIS and NIS as shown in Tab.8.The closeness coefficient represents a single index for simultaneous optimization of the multiple responses of DA considered in the present study.

Table 6:Weights of the response variable for DA

Table 7:Positive and negative ideal solutions for DA

Table 8:Closeness coefficient

It can be seen from Tab.8 that based on the closeness coefficient (Pi) value, experiment number 10 is the top ranked experiment (rank 1) which implies that the combination of machining parameters of this experiment, i.e., A3B2C4D3E1is the best one for multiple responses but it cannot be inferred that whether it is optimum or not.Consequently, optimum combination of the machining parameters for multiple responses was obtained through further analysis of closeness coefficient (Pi) including computation of S/N ratios, analysis of mean (ANOM) in terms of response table, and Pareto ANOVA as suggested in the literature [59].S/N ratio was calculated by considering higher-the-better characteristic for Pi as the objective is always to achieve its maximum value.The results of S/N ratio are shown in Tab.9.In ANOM, mean value of the S/N ratio at each level of each machining parameter was computed and the results are presented in the response table, i.e., Tab.10.The results of Pareto ANOVA are given in Tab.11 and its plot is shown in Fig.4 which reveals that the most important machining parameter is PON time (B)with the contribution ratio of 30.229%, second most significant parameter is dielectric fluid (E)with contribution ratio of 26.239%, third important parameter is current (A) with contribution ration of 22.045%.The least significant parameters are POFF time (C) and wire tension (D) with contribution ratio of 12.297% and 9.190% respectively.The levels of these parameters that yield optimum Pi value are selected from Tab.10 which shows significance of the response table.

Table 9:S/N ratios for Pi

Table 10:Response Table for S/N ratio for Pi

It may be noted that relationship formula of S/N ratio calculates its values which makes it independent of the characteristic type, i.e., whether it is lower-the-better or higher-the-better or nominal-the-better, maximum value of S/N ratio becomes the representative of desired performance.Consequently, it is evident from the results shown in Tabs.10 and 11 that the optimum combination of machining parameters for DA is B3E1A3C4D4.Thus, current at third level (3 A),PON time at third level (20µs), POFF time at fourth level (4µs), WT at fourth level (18 N) and dielectric fluid represented at first level [DM water(60%)+Ethylene Glycol(40%)] is the optimum combination for DA.

Table 11:Parto ANOVA analysis for Pi

Fig.4 shows that the most significant factor for dimensional deviation is PON time.The values of PON is the indicator for the duration of time for which the sparking will be on.This factor is therefore the sole factor responsible for the pulse energy dissipation in the machining zone and consequently is the most important factor for WEDM.The pulse on time at level 3 produced best condition for lowest DA and optimum accuracy.While during PON time period sustained machining takes place, but it is not paused, flushing-out of the debris will suffer and it will have detrimental effect on the machining accuracy and erosion of electrode wire.Thus, the PON time at level 3 optimum situation leading to better DA.

The dielectric fluid is the battery of electric power which is discharged as pulse energy during spark.In case of hybrid the nature of constituents such as oxygen and Al2O3nanoparticles in the present case will alter the manner, this energy is dissipated.The presence of oxygen involves surface oxidation and adjustment of thermal expansions of constituents.The presence of oxygen, thus, creates the violent dissipation of energy and may increase the dimensional deviations.Similarly, the nanoparticles being solid insulating particles alter the average spark gap during travel of spark plasma between the electric poles.In the present case the pure DM water could be able to create a moderate spark and limit the deviation and enhanced the DA.

The current creates an effect of third largest contribution.The effect that a parameter creates on a response mainly depends on:(a) The range of parameters values and (b) The types of selected response.Further, the current is a factor which is responsible for the intensity of energy which is discharged in the spark gap.Under the range of parameters selected in this study the current at level 3, i.e., 3 A created optimum deviation which gave the best over accuracy of the multi-performance characteristics.

The fourth most significant parameter for dimensional deviation turns out to be POFF time.The DA is poor at low POFF time.As the POFF time increases, the dimensional deviation is lower because the pause of spark dictated by higher POFF the available time for the flushing-out of the debris is sufficient and the spark gap is properly cleared.The trapped debris being bad conductor also makes the spark erratic and leading to reduced accuracy.Thus, POFF time at level 4 gave best dimensional accuracy.

Figure 4:Pareto ANOVA plot represents contribution ratio of parameters for Pi

Among all the parameters chosen for the investigation in this work, WT is the only electrically neutral parameter.But is produces an indirect effect on the quality of machined surface.WT should be at optimum value to achieve better accuracy.Forces due to sparking act in the transverse direction on wire electrode, and this affects the dimensional inaccuracy in the manufacturing parts.The wire experiences transverse force every times the spark occurs.Under transverse forces,the wire vibrates and the frequency of vibration is proportional to the square root of wire tension.Here lies the main reason, that the wire vibrates following a wave-form, and at low tension the amplitude of vibrations is large.Large amplitude results in greater change in spark gap.Thus, every wire tension value results in a path-change of the spark front.This issue becomes compounded by all factors which influence the transverse force on the wire (e.g., peak current,type of dielectric, PON and POFF time, etc.).The wire tension at low level results in amplitude of vibration of wire to increase and significantly reduces the accuracy.As WT increases, deflection in wire reduce and frequency increases which consequently makes the path of spark front stabilized and results in improvement the DA.

3.1 Validation of Results

The correctness of the obtained results is ensured by comparing the experimental value of the Pi with that of its corresponding predicted value at optimum combination.In order to estimate the predicted value of the Pi, Eq.(11) was used [60].

where, (PA3,PB3,PC4,PD4,PE1) represents the optimum level average values of Pi.Pavgstate the average of all the Pi values obtained from the experimental study.As a result of the calculations,it was estimated thatPopt=0.8779.

To calculate the experimental value of Pi, it was required to perform the experiment for these settings.However, Pi value cannot be determined by performing only one experiment.Therefore,to determine the experimental values of Pi at optimum combination, regression analysis was performed [61].Regression equation thus obtained for Pi is given in Eq.(12).The value of the coefficient of determination (R2) for Eq.(12) was found to be 0.9803 i.e., 98.03 % which shows a very strong relation between Pi and the various deviations (DID, DLAND, DOD and DTW).

Further, a confirmatory experiment was performed by setting input parameters at the optimum combination and the values of the DA were obtained which were as follows:DID=0.075 mm, DLAND=0.0021 mm, DOD=0.118 mm and DTW=0.0048 mm.Subsequently, Eq.(11)was employed to obtain the experimental value of Pi which was found to be 0.8652.Since, the predicted (0.8779) and experimental (0.8652) values of Pi are in close agreement and thus, it validates the optimal result obtained in this study.

4 Conclusion

Multiple response optimization for Wire EDM Nimonic alloy miniature gear was performed in this research.Taguchi L16orthogonal array was applied for design of experiments and Entropy TOPSIS was performed for multiple optimization of parameters.Pareto plot was employed to indicate the cumulative percentage contribution of the parameters to improve dimensional accuracy.The cutting parameters included in this research were current, pulse on time, pulse off time,wire tension and dielectric fluids while deviation in inner diameter, deviation in outer diameter,deviation in land and deviation in tooth width were main responses.Results indicated that pulse on time was the highly significant parameter for dimensional accuracy followed by dielectric fluid,current, pulse off time and wire tension.The PON time supplies time for the steady discharge of pulse energy and is the most significant parameter.The PON time at level 3 provided sufficient pulse energy discharge time to develop optimum accuracy.The dielectric fluid was next higher parameter to affect dimensional accuracy significantly.The pure DM water created conditions favorable for the optimum accuracy with the selected range all the parameters.The current which is responsible for the intensity of energy discharge in the gap was the third most significant parameter and the current at level 3 (3 A) created condition for minimum deviation and highest accuracy.The pulse OFF time provides essential pause in the discharge and makes way for the cooling and flushing-out of debris and increase the accuracy.The POFF time at level 4 produced conditions for the best dimensional accuracy.The only non-electric parameters, i.e., wire tension affected the machining performance mainly because of the deflections caused by the transverse forces generated due to spark.At lower tensions the magnitude of deflections are large and cause greater variation in the dimensions of machined part.The wire tension at level 4 created minimum deviation and best accuracy.Based on the study, the optimum machining parameters are pulse on time at 20µs, current at 3 A, pulse off time at 4µs, wire tension at 18 N, and dielectric fluid at first level, i.e., DM water(60%)+Ethylene Glycol(40%).

Funding Statement:The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University, for funding this work through research groups program under Grant No.(R.G.P.1/197/41).

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.