• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Multi criteria decision making through TOPSIS and COPRAS on drilling parameters of magnesium AZ91

    2022-12-26 02:36:06VaratharajuluMuthukannanDuraiselvamBhuvaneshKumarJayaprakashBaskar
    Journal of Magnesium and Alloys 2022年10期

    M.Varatharajulu,Muthukannan Duraiselvam,M.Bhuvanesh Kumar,G.Jayaprakash,N.Baskar

    aDepartment of Production Engineering,National Institute of Technology,Tiruchirappalli 620015,Tamilnadu,India

    b Department of Mechanical Engineering,A.V.C.College of Engineering,Mannampandal,609305,Tamilnadu,India

    cDepartment of Mechanical Engineering,Saranathan College of Engineering,Tiruchirapalli 620012,Tamilnadu,India

    d Department of Mechanical Engineering,Kongu Engineering College,Perundurai 638060,Tamilnadu,India

    Abstract Magnesium(Mg)alloys are extensively used in the automotive and aircraft industries due to their prominent properties.The selection of appropriate process parameters is an important decision to be made because of the cost reduction and quality improvement.This decision entails the selection of suitable process parameters concerning various conflicting factors,so it has to be addressed with the Multiple Criteria Decision Making(MCDM)method.Therefore,this work addresses the MCDM problem through the TOPSIS(Technique for Order Preference by Similarity to Ideal Solution)and COPRAS(COmplex PRoportional ASsessment)methods.The assessment carried out in the material Mg AZ91 with the Solid Carbide(SC)drill bit.The dependent parameters like drilling time,burr height,burr thickness,and roughness are considered with the independent parameters like spindle speed and feed rate.Drilling alternatives are ranked using the above said two methods and the results are evaluated.The optimum combination was found on the basis of TOPSIS and COPRAS for simultaneous minimization of all the responses which is found with a spindle speed of 4540 rpm and a feed rate of 0.076 mm/rev.The identical sequencing order was observed in TOPSIS and COPRAS method.The empirical model was developed through Box-Behnken design for each response.Superior empirical model developed for drilling time which is 3.959 times accurate than the conventional equation.The trends of various dependents based on the heterogeneity of various independents are not identical,these complex mechanisms are identified and reported.The optimized results of the Desirability Function Approach are greater accordance with the TOPSIS and COPRAS top rank.The confirmation results are observed with lesser deviation suggesting the selection of the above independent parameters.

    Keywords:Machining processes;Magnesium;TOPSIS;COPRAS;MCDM;Burr;Roughness.

    1.Introduction

    The applications of various advanced engineering materials have gained significant attention in the fields of automobile,aerospace,electronics,electrical and medical in recent years[1-3].This is due to the fact that they possess superior mechanical properties with less weight[4].Few specific materials and their alloys related to aluminum(Al),magnesium(Mg)are majorly used in structural applications because of high specific strength and low density.They are also difficult to machine due to their brittleness in nature[5-8].With these afore mentioned requirements,aviation industries demand other unique physical-mechanical properties such as rigidity and fatigue strength.The growing requirements necessitate the applications of composite materials[9-14].One of the major challenges while processing advanced materials is machining using conventional tools.The primary concern for researchers while machining advanced materials is the out-put characteristics in terms of Material Removal Rate(MRR),surface finish,tool wear,surface metallurgical characteristics and so on[15-18].Material deformation during machining plays a significant role in deciding the quality of the final component[19].Drilling is one of the fundamental machining operation carried out earlier to reaming,tapping and boring[20].The consistency and dimensional accuracy of a drilled hole is primarily affected by the cutting edge of a drill bit[21].This phenomenon depends on the cutting tool profile and material.The carbide twist drills[20]are performing comparatively better than any other tools in terms of life and surface damage[22,23].For assessing the machinability of Mg alloys using drilling,the tools are further coated and used[24].During drilling,some left out noticed along the tool periphery in the work piece due to the plastic deformation usually called as burr[25].It has many adverse feature-like minor wounds of assembly labors,inappropriate mating and functional problems[26].Hence,it requires deburring process,which additional increases approximately 20–30% of manufacturing cost[27].Controlling the burr height and thickness during machining itself is a necessary study which is highlighted in the present work.

    Multi-criteria decision-making(MCDM)methods are useful resources to offer solutions to decision-making problems,including contrasting and multipurpose goals,both in real life and in professional environments[28,29].These challenges involve a small range of solutions,evidently recognized at the outset of the decision-making process.Each alternative is shown by its success in a variety of parameters.The complexity of a choice-based problem results from the existence of more than one criterion for the determination of alternatives.The approach depends heavily on the desires of the decision-makers(DMs)[30,31].To decide the favorable and unfavorable ideal alternatives from the decision matrix,the Technique for Order Preference by Similarity to Ideal Solution(TOPSIS)method was used[32].The TOPSIS process is applied successfully for the selection of sequential parameters for the drilling process over hybrid polymer matrix composite[33].The combination of CCD-TOPSIS with the surface response method demonstrates its superiority in optimizing process parameters while enhancing the efficiency characteristics of the GFRP marine-grade drilling process[34].The Taguchi is combined with TOPSIS for multi-response optimization of rotary ultrasonic drilling of BK7 concerning taper,chipping width and MRR[35].The process parameters of micro EDM of AISI 304 steel analyzed by the combination of TOPSIS and ANOVA,establishes its efficiency[36].The geometrical characteristics and performance measures are optimized in cryogenically chilled micro EDM drilling process using TOPSIS[37].Solving a multi-response problem in the drilling process,the TOPSIS identified as an effective method and it was applied in cryogenic drilling on Ti-6Al-4V[38].The hole drilling of thermal barrier coated Nickel superalloy C263 using a laser was optimized by TOPSIS[39].The DOE based experimental study used in Fuzzy-TOPSIS methodology,successfully applied to determine the optimized parameters in the laser drilling process[40].

    The multi-objective Taguchi and TOPSIS were applied to optimize the process parameters in Carbon Fiber Reinforced Polymer(CFRP)composite polymer and identified the quadratic model using Response Surface Analysis(RSA)[41].Identification of appropriate cutting fluids beneficial for machining operations and increases the productivity,life of the tool and quality of the components.In order to do so,the modified similarity-based method was applied effectively via TOPSIS[42].The different conventional and nonconventional manufacturing processes successfully applied the TOPSIS algorithm for optimizing their process parameter.In addition,this algorithm is very useful in chemical engineering,maintenance engineering and civil engineering[43].Petkovi′c,Madi′c and Radenkovi′c[44]carried out a comparative analysis of the use of Weighted Aggregated Sum Product ASsessment(WASPAS)and COmplex PRoportional ASsessment(COPRAS)in MCDM methods to assist in the identification of the most effective non-conventional machining processes for ceramic machining.Taguchi based COPRAS was presented to identify the optimal combination of parameters in turning operation over stainless steel 304.The COPRAS method is successfully applied to convert the multi-objective into a single objective to carry out the examination[45].Five different normalization procedures evaluated by Yazdani,Jahan and Zavadskas[46]and suggested the superior one.The COPRAS method is applied for finding the best alternative in the drilling process on Al alloy using a solid carbide drill bit with a high-pressure coolant[47].In order to address the shortcomings of the fuzzy AHP methodology,the integration of a stable fuzzy AHP and fuzzy COPRAS was proposed for the machine tool selection.The specification of the consistency ratio(CR)is prevented while the fuzzy linguistic choice relationship is used to incorporate into the AHP.The suggested method of machine tool selection is intended to be conveniently extended[48].The stochastic type of COPRAS is effectively applied for selecting the suitable service provider in the cargo industry[31].

    Development of relation between independent and dependent parameters are essential in the case of machining with newer material like Mg AZ91[2].Minimal cutting forces observed in the case of AZ91 when it was assessed with AZ31 during turning operation,however the machining characteristics of AZ alloys during drilling was not described so far[5].Assessment by past researches reveals that,not much research has been done in the Mg alloy during the drilling process through TOPSIS and COPRAS while MCDM is a concern.It is worth looking at the machining behavior of Mg AZ91.This study focuses on the identification of significant aspects through TOPSIS and COPRAS and comparing the ranking with Desirability Function Approach(DFA).Though both the methods are adopted widely in the other machining parameters for optimization,ranking using TOPSIS is based on the relative closeness coefficient and COPRAS is based on utility degree[49].Hence there may be differences with their results that may not be a standard one.Different MCDM tools may result in a different sequence of results that may not be acceptable for generic concern.The novel aspect of this presentstudy is the comparative analysis of two MCDM tools namely TOPSIS and COPRAS for the selection of generic drilling parameters to machine Mg alloys.The complex behavior mechanism was analyzed with the appropriate technical reason.

    Fig.1.Methodology.

    2.Materials and methods

    The methodology employed in the current study is shown in Fig.1.The proposed method begins with an experimental work based on the coded and actual factors as shown in Table 1.The X6323 model ALTO make three axes vertical machining center was used to perform the experimental work.Addison make micro grain solid carbide straight shank twistdrill bit of 6 mm diameter(Jobber series)(IS5101–1991/DIN 338)is used for the drilling process.[Flute length 57 mm,Overall length 93 mm,helix angle 300,point angle 1180].The work material used is Mg AZ91 alloy and the elemental contents are presented in Table 2.To confirm the designation of the selected material,energy dispersive spectroscopy has beenanalyzed as shown in Fig.2.Further,the material properties are confirmed by the basic microstructure,in which the second phase precipitate across the grain boundary can be noted,as shown in Fig.3.The dimensions of the prepared work specimen are 150 mm(L)×70 mm(W)×10 mm(T).The range of machining parameters is decided based on the tool manufacturer’s recommendation and the competency of the machine tool.The different set of independent combinations(17 run)derived based on Box-Behnken Design(BBD)and experimented(Table 3).For each experimental run,drilling time is measured using a stop watch.Before experimentation,the drilling time is calculated based on the basic formula.The ratio of length of workpiece(in mm)with cross product of spindle speed(in rpm)and feed rate(mm/rev.)providing drilling time.However,these values are noticed with larger disparities and the average deviation found larger with 16.687%(with experimental drilling time),which inquisitive to find a superior model.

    Table 1Levels and ranges of input parameters.

    Table 2Chemical composition of Mg AZ91.

    Table 3Number of run versus coded and actual factors(Box-Behnken Design).

    Fig.2.Energy Dispersive Spectroscopy report.

    Fig.3.SEM images of work piece with(a)200X Magnification(b)500X Magnification.

    After the drilling process,the burr height at the entry and exit of the hole is measured using a coordinated measuring machine named HELMEL 80–19(Model 216–142).The measurement of burr height is determined based on the average taken from four values measured over different places of both entry and exit portions respectively.The measurement of burr thickness is done using a digital vernier caliper(Mitutoyo CD 6 CS3).The thickness of the burr is considered as the difference between the external and internal diameter of the projected burr.The roughness of the resulting surface inside the hole is measured in different places using a surface roughness testing machine Mitutoyo SJ 210 and the average is taken for analysis.The responses measured during and after the experiments are presented in Table 4.The complete experimental setup and process flow are illustrated in Fig.4.The determination of entry and exit burr height using CMM,gauging surface roughness using surface roughness tester and the burr development in the workpiece material were illustrated in Fig.5.

    Table 4Experimental data.

    The objective of this study is to minimize multiple responses as stated above by identifying appropriate spindle speed and feed rate,which is more likely to be a multiresponse optimization problem.Hence to make a decision on this multi-response problem,TOPSIS and COPRAS techniques are incorporated.

    3.Result and discussion

    3.1.TOPSIS methodology

    TOPSIS is an efficient and advanced MCDM methodology,which was first introduced by Hwang and Yoon[50]to obtain the finest choice based on the compromise solution principle.The acceptable solution can be seen as preferring the answer with the shortest route from the favorable ultimate limit and the longest route from the unfavorable ultimate limit[51].

    The subsequent procedures are used for choosing the right alternatives using the TOPSIS algorithm[52,53].

    Step 1:The identification of essential attributes(dependent and independent parameters)has to be done primary.Generally,the dependent parameters desire maximization considered as most preferable and the dependent parameters wishes minimization reckoned as least preferable.

    Surmising the present study,all the parameters are considered are least preferable attribute(for minimization)while none of the parameter considered as a most preferable attribute.Because the diminishing the responses leads to an increase in the productivity and quality of the component.

    Step 2:A matrix more often termed as decision matrix used to represent all the information,which has i row(m-alternatives)and j column(n-criteria).

    The decision matrix[D17×6]is given in Eq.(1).

    Step 3:The following formula[Eq.(2)]is used to determine the elements in the normalized matrixNij[53].

    The calculation of normalized data[Eq.(3)]for first and last element is as follows

    The calculated normalized decision matrix[N17×6]is presented in Eq.(4).

    Step 4:The cross multiplication of individual weight with the respective column of normalized decision matrix will provide weighted normalized decision matrix.

    whereWjis the weight criteria andNijis the normalized matrix.The weight(Wj)of each criterion is decided based on the expertise,which is 0.3 for drilling time,surface roughness and for rest it was 0.1.The weighted normalized value W17×6is calculated using Eq.(5)is given as Eq.(6).

    The weighted normalized matrix[W17×6]is presented in Eq.(7).

    Step 5:Estimation of the positive ultimate solution(A??)and the negative ultimate solution(A?).These are analyzed using following formula,Eq.(8)and Eq.(9):

    J=1,2,3,…,n-whereJis related with the benefit criteriaJ’=1,2,3,…,n-whereJ’is related with the cost criteria.The positive ideal solution(A??)is considered based on minimization whereas the unfavorable ideal solution(A?)is reckoned based on maximization.Deeming the present study all the six responses have to be minimized for increasing the productivity and quality of the component.The respective favorable and unfavorable ideal solution are given below[Eq.(10)]

    Step 6:The separation measure is obtained using Eq.(11).The same for each preference from the favorable ultimate one is specified by Eq.(12).

    Similarly,the unfavorable ultimate one is calculated using Eq.(13)and it is specified by Eq.(14).

    Step 7:The relative closeness is examined using Eq.(15)for each alternative and it is described in Eq.(16).

    The best alternative is chosen based onCi?value and closeness to ultimate solution.The separation measure of favorable,unfavorable ultimate solution and relative closeness value are presented in Table 5.

    Table 5Separation measure of positive,negative ideal solutions and relative closeness value.

    Table 6Relative weight,utility degree and its rank.

    Fig.4.Experimental setup and process flow.

    Step 8:Ranking based on relative closeness value.

    According to the calculated results,the complete ranking of the alternative is obtained as 17–16–15–12–13–14–5–7–2–10–6–4–3–8–1–9–11.The optimal experimental design is sequenced as follows 15>9>13>12>7>11>8>14>16>10>17>4>5>6>3>2>1.It means that the best-combined input is in run 15,and the worst is in run 1.Run 15 has 91.33% relative closeness and run 1 has 33.7% relative closeness.The superior combination to minimize the responses identified in run 15 succeeding to run 9 and 13.The optimal combination was identified with a spindle speed of 4540 rpm and a feed rate of 0.076 mm/rev.Therefore,it is suggested to select the above sequence in order to diminish the manufacturing cost and to enhance the quality of the component.

    3.2.COPRAS methodology

    Zavadskas,Kaklauskas and Sarka[54]conceived the COPRAS which entails direct and proportional reliance on the importance and usefulness of the alternatives available in the presence of mutually contrasting parameters.COPRAS en-compasses the success of variants in terms of various parameters and related weights by using the step-by-step ranking and determining the system of choice of variants in terms of their relevance and their degree of usefulness.Variety of dynamic decision-making problems COPRAS proposes to take advantage of this approach in the application of engineering[55].The primary benefit of this technique is its ease of use and friendliness.However,when dealing with qualitative metrics and characteristics,it has its drawbacks[56].

    Fig.5.Response measurement(a)burr height measurement using CMM(b)surface roughness measurement using SR tester(c)burr development at the tool exit portion(d)sample roughness measurement for run number 9.

    The subsequent steps are used for choosing the right alternatives using the COPRAS algorithm[57].

    Step 1:Identification of the objective and the significant attributes(same as that of TOPSIS).

    Step 2:Representation of decision matrix[Eq.(1)].

    Step 3:The following formula[Eq.(17)]is used to calculate the elements in the normalized matrixNij[57].In order to bring the unit free responses,the normalization is essential for comparing them[58].

    The calculation of normalized data for first and last element is as follows[Eq.(18)].

    Eq.(19)shows the normalized decision matrix[N17×6].The calculated decision matrix correctness can be crosschecked on the following basis.The summation of an individual attribute should be equal to the ratio of unity to the number of runs.Referring Eq.(19),the average of each attribute is 0.0588 is equal to the ratio of unit to number of run(1/17=0.0588).

    Step 4:The procedure followed in TOPSIS is same in determining the weighted normalized decision matrix[Eq.(20)],even the same weight is considered in order to keep them constant.

    The perceived weighted normalized matrix[W17×6]is given in Eq.(21).

    Step 5:Calculation ofPi,which is identified by the summation of attributes those need to be maximized.That is,the summation over the beneficial criteria[59].

    In the above Eq.(22),k is the number of responses which necessitate to be maximized.Concerning the current study,none of the attribute has to be maximized;therefore,all thePivalue is equal to zero.

    Step 6:Calculation ofRi,which is calculated by the summation of attributes those need to be minimized.That is,the summation over the non-beneficial criteria[59].

    In the above formula Eq.(23),(m-k)is the number of responses which required to be lessened.Conceiving the current study,all of the attributes have to be minimized,due to their adverse features.The sample calulation is presented in Eq.(24).

    Step 7:Perceiving the diminutive value ofRi[Eq.(25)].

    Step 8:Determination of the relative weight of each responsesQi.The relative noteworthy value of a response shows the measure of contentment achieves by that response.The larger is the priority of the alternative is based on the higher value ofQi.The superior selection among the alternatives is based on the largest relative valueQmax[59].

    Applying the values calculated through the Eq.(26),Eq.(28)and Eq.(29)in Eq.(27),the relative weight can be calculated,which is presented in Eq.(30).

    Rather,Eq.(25)some researcher will follow the following Eq.(31),however both the equation will return the same results.

    Step 9:Computation of the optimality criterionQmax

    Step 10:Computation of the precedence of the work.The major weight(relative weight of response)Qi,the longer is the primacy(rank)of the work.SurmisingQmax,the agreement degree is the largest.

    Step 11:Estimation of the utility degree of each response.The degree of utility of the responses,which assist to a complete ranking of responses,is calculated by evaluating the preferences of all responses with the most effective[Eq.(34)]and can be defined as follows[59]:

    From the application of Eq.(31)and Eq.(32),the above one is examined[Eq.(35)].These utility degree values of each alternative vary from 0% to 100%.

    Based on the evidence presented on the requirements for the response of the drilling process,reasonable strategies can be sought to increase efficiency and minimize costs.According to the calculated results,the complete ranking of the alternative is obtained as 17–16–15–12–13–14–6–7–2–11–5–3–4–8–1–9–10.The optimal experimental design is sequenced as follows 15>9>12>13>11>7>8>14>16>17>10>4>5>6>3>2>1.It means that the bestcombined input is in run 15,and the worst is in run 1.Run 15 is the best alternative with 100% utility degree and run 1 is the worst alternative with 44% utility degree(Table 6).The superior combination to minimize the responses identified in run 15 succeeding to run 9 and 12.The optimal combination was identified with a spindle speed of 4540 rpm and a feed rate of 0.076 mm/rev.Therefore,it is suggested to select the above sequence in order to diminish the manufacturing cost and to augment the quality of the component.

    Table 7Confirmation test result.

    The ranking performance of TOPSIS has greater concurrence with the COPRAS methods.There are 64.71% of ranking coincide with the COPRAS,deeming the higher and lower order it was identical.The 1st,2nd,7th,8th,9th,12th to 17th rankings are the same(Fig.6).Distinguishing ranking was observed for 35.29% of data;however,the sequence is coherent.The significant similarities between the two ranking methods are observed.The computational techniques of the two approaches are simple and easy in terms of recognizing and implementing these methods for comparing the alternatives and choosing the drilling parameters.

    Fig.6.Comparison of TOPSIS and COPRAS ranking.

    3.3.Modeling through Response Surface Methodology(RSM)

    The empirical models have been developed for the present study using BBD[2,3].The coefficient of determination(R2)value obtained for various responses are 99.68,97.64,89.2,94.77,99.57 and 94.13 respectively,seem better the prediction ability of the developed model.The developed individual responses are presented in Eqs.(36)to(41).In the following equation,dtdenotes drilling time,bhenrepresents entry burr height,bhexspecifies exit burr height,btenindicates entry burr thickness,btexdesignates exit burr thickness,srspecifies surface roughness,Ndenotes spindle speed andFrepresents feed rate.

    The accordance between the experimental and RSM data is evident from Fig.7.The 17 number of runs are converted into 9 by taking the average in the middle value of the feed rate.Reckoning the BBD,the run 2 and 3 are similar,run 5 and 6 are the same,from 7 to 11 are identical,run 12 and 13 are equal,run 15 and 16 are indistinguishable.Therefore,Fig.7 has a 9 number of runs in the abscissa.Based on the accordance,these empirical models can be applied for the identification of suitable inputs by process planning engineer.Once the superior model is developed,the BBD is further extended to optimization through DFA.The procedure for DFA followed from the past literature[3,60].The result of DFA is greater concurrence with the top rank of TOPSIS and COPRAS.Conceiving Eq.(36),the RSM drilling time is determined,subsequently the average percentage of deviation is calculated with 4.215% which is 3.959 times superior than drilling time(Fig.7a)which is calculated based on basic formula.

    3.4.Characterization

    The behavior study of the spindle speed and feed rate over all the responses are performed through RSM 3D graph and individual value plot(Fig.8).Considering drilling time(Fig.8a,b),an increase in spindle speed decreases the drilling time at all levels of feed rate.This is because of increasedshearing rate over the material during machining[25,61].The recent study found that,enhancing spindle speed propagating incisaling action.Incisaling action is introduced by the recent work which includes material removal along the shearing action[62].During the lower range of feed rate sudden drop-in drilling time was observed and for a medium and larger range of feed rate,reduction in drilling time is gradual.Increasing feed rate declines drilling time till intermittent range thereafter increases while the spindle speed is in minimal and intermittent range.Till intermittent range of feed rate the augmented dribbling action propagating material removal causes drop-off in drilling time,beyond intermittent range due to the same phenomena the material removal found larger which leads to poor extrusion of chips causes ascent in drilling time.While the spindle speed is high,the amplifying feed rate lessens drilling time that too with plodding way.The combination of minor spindle speed and minor feed rate produces larger drilling time and on the other hand shorter drilling time identified at the combination of major spindle speed and major feed rate.

    Fig.7.Validation of experimental results for various responses(a)drilling time(b)entry burr height(c)exit burr height(d)entry burr thickness(e)exit burr thickness(f)surface roughness.

    Fig.8.Spindle speed versus feed rate on various responses(a)Drilling time:3D graph(b)Drilling time:Individual value plot(c)Entry burr height:3D graph(d)Entry burr height:Individual value plot(e)Exit burr height:3D graph(f)Exit burr height:Individual value plot(g)Entry burr thickness:3D graph(h)Entry burr thickness:Individual value plot(i)Exit burr thickness:3D graph(j)Exit burr thickness:Individual value plot(k)Surface roughness:3D graph(l)Surface roughness:Individual value plot.

    Fig.8.Continued

    In view of entry burr height(Fig.8c & d),entry burr thickness(Fig.8g & h),exit burr thickness(Fig.8i & j)and surface roughness(Fig.8k & l),the spindle speed has a sinusoidal relation and larger is the spindle speed shorter is the entry burr height till 1960 rpm thereafter matures till 3680 rpm;beyond this point,the response condenses.This occurrence is found in all ranges of feed rate.Increase in spindle speed till 1960 rpm increases incisaling rate causes minimal responses,thereafter further increase in incisaling rate increases more material removal results in lack of room for movement of chips and consequences larger responses.Once,it reaches 3680 rpm,due to higher centrifugal material ejection the chips are thrown away from the base material causes lesser responses.Type of chip generated during perforation process has a significant role in the response behavior pattern[16,61].Types of chips reported in the past study[63]are ribbon type,tubular type,cork screw type,helical type,spiral type,arc type,elemental type and needle type.All these types are in three different form and they are short,long and snarled.Concerning the present study,chips are identified with elemental type,short and snarled helical type.Feed rate has erratic relation over the responses[other than drilling time and exit burr height].For the lower and middle range of spindle speed;the longer is the feed rate,the shorter is the responses.This could be because of high temperature generation due to faster tool penetration which leads to thermal softening resulting minimal response[64].For the higher range of spindle speed;lengthening the feed rate shrinkages the response till the middle range thereafter extends.When both the independent are high,due to the increased tool vibration and wobbling the chips are left around the periphery even it has larger material ejection and higher thermal softening.It has greater concurrence with the past study[65],due to high heat distribution on surface causes shorter chips which left around the periphery.Least responses are observed during the lower spindle speed(1960 rpm)at different feed rate and greater responses are noticed at higher spindle speed(3680 rpm)at any ranges of feed rate.The minimal entry burr height(Fig.8d)observed in the combination of higher speed and medium feed,larger entry burr height attained in the lower speed and feed.The arrangement of diminutive independents are producing least entry burr thickness(Fig.8h)and it was largest when those independents are maximal.The mixture of intermittent speed and higher feed produced minimal exit burr thickness(Fig.8j).The minimal surface roughness(Fig.8l)observed in the amalgamation of higher speed and lower feed.These two responses(exit burr thickness and surface roughness)are found larger with lower speed and intermittent feed.

    For all ranges of spindle speed;increase in feed rate,decreases exit burr height till medium range,later increases(Fig.8e & f).This could be due to high temperature generation due to faster tool penetration which leads to thermal softening resulting minimal exit burr height till intermittent feed range,subsequently due to the amplified frictional force among workpiece material and drill bit causing deficient strain,due to which breaking of chips decreased and resulted in retentive burrs at exit of holes[64].The poorer exit burr height is produced during all ranges of the spindle speed at medium feed rate and greater exit burr height is produced by all ranges of the spindle speed and lower range of feed rate.Burr generation found minimal at entry side when it was evaluated with exit side.Initially shearing noticed due to plastic deformation,when drill bit penetrates in the material with chisel edge followed by cutting lip and peripheral cutting edge[66].Till cutting tool engages with workpiece,it produces elemental chips and fan chips which ousted at the top surface causes minimal burr height and thickness.Once,it engages with the workpiece surface,sheared material starts to move through flute and some chips are left with the periphery causing noticeable burr height and thickness in the entry side.The thin lamina of material is bent and exiled when drill bit near the bottom due to tearing[67].Due to negligible shearing and ousting,the chips are left around the periphery causing higher burr height and thickness in the exit side.

    3.5.Confirmation test

    To verify the enhancement of the output characteristics,validation test has been carried out of the drilling process with an optimum combination of parameters.Based on TOPSIS,COPRAS and DFA,the optimal drilling independents is found at experiment number 15 with the spindle speed of 4540 rpm and a feed rate of 0.076 mm/rev.With these suggested independents,the confirmation test carried out with three trails and the average value taken for evaluation.The results of the confirmation test are presented in Table 7.Theconfirmation test results are concurrent with the initial values having an error percent of 4.622,confirms the validation of suggested inputs.

    4.Conclusion

    In the current study TOPSIS and COPRAS algorithms are applied for multi-objective optimization to determine the superior amalgamation of drilling parameters such as spindle speed and feed rate for simultaneous minimization of all the responses while drilling Mg AZ91 with carbide tools.

    ·The TOPSIS method is worn to choose the greatest arrangement of drilling parameters.Based on the value of relative closeness,the responses are sorted as follows:17–16–15–12–13–14–5–7–2–10–6–4–3–8–1–9–11.The optimal experimental design is sequenced as follows 15>9>13>12>7>11>8>14>16>10>17>4>5>6>3>2>1.

    ·The same arrangement based on COPRAS is as follows:17–16–15–12–13–14–6–7–2–11–5–3–4–8–1–9–10.The optimal experimental design is sequenced as follows 15>9>12>13>11>7>8>14>16>17>10>4>5>6>3>2>1.

    ·From both the algorithm,the best-combined input is observed in run 15,and the worst is observed in run 1.The sequencing performances of TOPSIS and COPRAS methods are identical.

    ·The optimal arrangement identified based on TOPSIS and COPRAS of drilling independent for simultaneous minimization of all the responses are identified with a spindle speed of 4540 rpm and feed rate of 0.076 mm/rev.The same process parameter identified as optimum input through DFA.

    ·The result of the confirmation test validates the proposed independent parameter with minimal error percentage.

    ·The empirical model is developed with a higher reliability measure through BBD[Eqs.(36)to(41)],which can be applied for finding the independents without experimentation.

    ·For drilling time,a superior mathematical model was established that is 3.959 times more accurate than the standard model drilling time.

    ·The rate of incisaling found larger with an addition of spindle speed causes larger material removal rate,resulting drilling in a faster rate.Poor extrusion of chips making the process time consuming.Centrifugal material ejection found with larger speed is the reason for minimal responses in general.

    ·Dribbling action found with feed rate influencing on material removal,extending the study with lack of room for movement of chips affecting response.Faster tool penetration found with larger feed rate resulting minimal responses.

    ·The experiment can be done at a faster rate with a higher combination of both the response.The lower combination of both the responses produces minimal entry burr thickness.The combination of lower spindle speed with a higher feed rate produces entry burr height in a minimal way.

    ·The higher surface quality components achieved through the mixture of higher spindle speed with a lower feed rate.The intermittent spindle speed associated with a higher feed rate minimizes the exit burr height and exit burr thickness.

    The effect of process parameters of drilling parameters using coated tools is not considered in the present work,might be worth some investigation.Two different process input parameters alone considered in the present investigation.Inclusion of other parameters like material thickness,drill diameter,drill bit point angle,drill type,cutting fluid may enhance the rate of prediction.The responses considered are also limited.Consideration of thrust force,torque,material removal rate,tool wear and study of chip formation would be a good extension.

    国产精品综合久久久久久久免费 | 一二三四社区在线视频社区8| 在线观看一区二区三区激情| 国产在线一区二区三区精| 老熟妇仑乱视频hdxx| 国产精品久久久久久精品古装| 国产精品一区二区精品视频观看| av网站免费在线观看视频| videosex国产| av有码第一页| 亚洲精品美女久久久久99蜜臀| 1024香蕉在线观看| 国产不卡av网站在线观看| 欧美精品高潮呻吟av久久| 女人高潮潮喷娇喘18禁视频| 日韩免费av在线播放| 悠悠久久av| 激情视频va一区二区三区| 水蜜桃什么品种好| 免费在线观看视频国产中文字幕亚洲| 日韩免费高清中文字幕av| 国产亚洲精品一区二区www | 少妇被粗大的猛进出69影院| 美女午夜性视频免费| 波多野结衣av一区二区av| 色在线成人网| 午夜影院日韩av| 人人妻,人人澡人人爽秒播| 久久久久久久国产电影| 国产高清视频在线播放一区| 天天躁夜夜躁狠狠躁躁| 亚洲少妇的诱惑av| 在线观看一区二区三区激情| 国产又色又爽无遮挡免费看| 99re在线观看精品视频| 久热这里只有精品99| 热99久久久久精品小说推荐| 91麻豆精品激情在线观看国产 | 人妻一区二区av| 久热这里只有精品99| 亚洲,欧美精品.| 99在线人妻在线中文字幕 | 一级毛片精品| 久久国产精品人妻蜜桃| 夜夜躁狠狠躁天天躁| 精品一区二区三区av网在线观看| 看免费av毛片| 久久国产亚洲av麻豆专区| 中文字幕高清在线视频| 午夜日韩欧美国产| 精品国产国语对白av| 午夜福利影视在线免费观看| 国产乱人伦免费视频| 亚洲五月色婷婷综合| 丰满饥渴人妻一区二区三| 亚洲国产欧美日韩在线播放| 免费在线观看影片大全网站| 涩涩av久久男人的天堂| 亚洲第一av免费看| 天天影视国产精品| 在线观看一区二区三区激情| av欧美777| 人人妻人人爽人人添夜夜欢视频| 亚洲av成人av| 五月开心婷婷网| 国产男女超爽视频在线观看| 香蕉丝袜av| 丝袜美腿诱惑在线| 国产男女超爽视频在线观看| 亚洲av熟女| 色精品久久人妻99蜜桃| 桃红色精品国产亚洲av| 一进一出抽搐动态| e午夜精品久久久久久久| 婷婷丁香在线五月| 国产精品欧美亚洲77777| 欧美国产精品一级二级三级| 免费在线观看日本一区| 窝窝影院91人妻| 好看av亚洲va欧美ⅴa在| 亚洲av成人一区二区三| 亚洲中文日韩欧美视频| 国产成人精品久久二区二区免费| 亚洲一区二区三区不卡视频| 黄网站色视频无遮挡免费观看| 免费日韩欧美在线观看| 在线十欧美十亚洲十日本专区| 成年人午夜在线观看视频| 欧美 亚洲 国产 日韩一| 一边摸一边做爽爽视频免费| 国产无遮挡羞羞视频在线观看| 大型黄色视频在线免费观看| 亚洲国产毛片av蜜桃av| 亚洲国产看品久久| 宅男免费午夜| 久久久水蜜桃国产精品网| 亚洲精品粉嫩美女一区| 一二三四社区在线视频社区8| 91成年电影在线观看| 狠狠婷婷综合久久久久久88av| 69av精品久久久久久| 如日韩欧美国产精品一区二区三区| 一进一出抽搐动态| 国产精品欧美亚洲77777| 精品高清国产在线一区| av免费在线观看网站| 国产精品亚洲一级av第二区| 村上凉子中文字幕在线| 精品少妇一区二区三区视频日本电影| 麻豆成人av在线观看| 美女视频免费永久观看网站| 12—13女人毛片做爰片一| 十八禁高潮呻吟视频| 免费黄频网站在线观看国产| 国产在线观看jvid| 国产精品久久久人人做人人爽| 精品人妻熟女毛片av久久网站| 19禁男女啪啪无遮挡网站| 久久中文字幕人妻熟女| 亚洲av熟女| 国产欧美日韩一区二区精品| 色婷婷av一区二区三区视频| 夫妻午夜视频| 亚洲在线自拍视频| 9191精品国产免费久久| 高清欧美精品videossex| 免费在线观看影片大全网站| 亚洲美女黄片视频| 人妻 亚洲 视频| 欧美黑人欧美精品刺激| 国产三级黄色录像| 久久亚洲真实| 99久久国产精品久久久| 亚洲国产中文字幕在线视频| 人妻丰满熟妇av一区二区三区 | 国产精品一区二区精品视频观看| 女人被躁到高潮嗷嗷叫费观| 免费少妇av软件| 久久中文看片网| 国产精品自产拍在线观看55亚洲 | 午夜免费鲁丝| 在线观看免费高清a一片| 一级a爱视频在线免费观看| 一进一出抽搐gif免费好疼 | 国产男靠女视频免费网站| 午夜福利乱码中文字幕| 大香蕉久久网| a级毛片在线看网站| av片东京热男人的天堂| 久久热在线av| 欧美 亚洲 国产 日韩一| 欧美亚洲日本最大视频资源| 黑人巨大精品欧美一区二区蜜桃| 国产在线观看jvid| 高清视频免费观看一区二区| 香蕉丝袜av| 18在线观看网站| 久久久久久亚洲精品国产蜜桃av| 两个人看的免费小视频| a级毛片黄视频| 韩国精品一区二区三区| 亚洲精品一二三| 亚洲一码二码三码区别大吗| 免费观看a级毛片全部| 亚洲avbb在线观看| 亚洲精品av麻豆狂野| 一二三四在线观看免费中文在| 极品少妇高潮喷水抽搐| 天堂动漫精品| 亚洲九九香蕉| 久久午夜亚洲精品久久| 成人三级做爰电影| 国产1区2区3区精品| 色婷婷久久久亚洲欧美| 久久久久久人人人人人| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲,欧美精品.| 大陆偷拍与自拍| 免费观看a级毛片全部| 夫妻午夜视频| 好看av亚洲va欧美ⅴa在| 午夜两性在线视频| 国产精品久久久av美女十八| 在线天堂中文资源库| 在线观看免费日韩欧美大片| 久久久久视频综合| 人人澡人人妻人| 激情在线观看视频在线高清 | 久久午夜亚洲精品久久| 一级片'在线观看视频| 一级毛片精品| 老司机靠b影院| videos熟女内射| 亚洲熟妇熟女久久| 欧美精品人与动牲交sv欧美| 亚洲精品乱久久久久久| 精品国内亚洲2022精品成人 | 国产一区有黄有色的免费视频| 日本黄色日本黄色录像| 欧美日韩成人在线一区二区| 99久久综合精品五月天人人| 黄色怎么调成土黄色| 亚洲精品成人av观看孕妇| 亚洲美女黄片视频| 女性生殖器流出的白浆| 久久精品成人免费网站| 中文字幕色久视频| 精品高清国产在线一区| 999久久久精品免费观看国产| 在线观看免费视频网站a站| 性色av乱码一区二区三区2| 精品一区二区三区四区五区乱码| 精品久久久久久久久久免费视频 | 99热网站在线观看| 看片在线看免费视频| 成人特级黄色片久久久久久久| 交换朋友夫妻互换小说| 极品人妻少妇av视频| 一级毛片高清免费大全| 悠悠久久av| 亚洲精品粉嫩美女一区| av福利片在线| 欧美国产精品一级二级三级| 免费日韩欧美在线观看| 精品人妻1区二区| 亚洲,欧美精品.| 在线观看午夜福利视频| 欧美成人免费av一区二区三区 | 国产男靠女视频免费网站| 亚洲精品国产精品久久久不卡| 午夜精品在线福利| av一本久久久久| 女性被躁到高潮视频| 国产日韩欧美亚洲二区| 午夜福利免费观看在线| 麻豆国产av国片精品| 久久久久精品国产欧美久久久| 亚洲自偷自拍图片 自拍| 精品亚洲成a人片在线观看| 欧美黑人欧美精品刺激| 在线观看午夜福利视频| 国产亚洲欧美精品永久| 日韩免费高清中文字幕av| 欧美精品高潮呻吟av久久| 国产成人精品在线电影| 国产欧美日韩一区二区精品| 亚洲精品国产色婷婷电影| 在线观看66精品国产| 国产蜜桃级精品一区二区三区 | 美女福利国产在线| 超碰97精品在线观看| 美女高潮喷水抽搐中文字幕| 一级片免费观看大全| 久久草成人影院| 欧美精品啪啪一区二区三区| 日韩成人在线观看一区二区三区| 精品国产乱码久久久久久男人| 国产区一区二久久| 欧美日韩国产mv在线观看视频| av天堂在线播放| 国产三级黄色录像| 欧美日韩一级在线毛片| 国产人伦9x9x在线观看| 久久国产精品影院| 国产av一区二区精品久久| 老鸭窝网址在线观看| 欧洲精品卡2卡3卡4卡5卡区| 国产精品一区二区精品视频观看| 成年动漫av网址| 亚洲一码二码三码区别大吗| av网站在线播放免费| 国产精品99久久99久久久不卡| avwww免费| a级毛片在线看网站| 大片电影免费在线观看免费| 亚洲午夜精品一区,二区,三区| 黄色女人牲交| 99riav亚洲国产免费| 黑人欧美特级aaaaaa片| 妹子高潮喷水视频| 人人妻人人爽人人添夜夜欢视频| 两性夫妻黄色片| 国产99久久九九免费精品| 伊人久久大香线蕉亚洲五| av有码第一页| 色婷婷av一区二区三区视频| 国产免费男女视频| 亚洲欧美色中文字幕在线| 男女免费视频国产| 成人三级做爰电影| 一二三四在线观看免费中文在| 国产精品免费一区二区三区在线 | 欧美久久黑人一区二区| 热re99久久精品国产66热6| 亚洲男人天堂网一区| 香蕉国产在线看| 中文字幕av电影在线播放| av线在线观看网站| 日韩视频一区二区在线观看| 国产区一区二久久| 下体分泌物呈黄色| av视频免费观看在线观看| 18在线观看网站| 在线播放国产精品三级| 国产成人精品久久二区二区免费| 久久久久久久久免费视频了| 久久久精品国产亚洲av高清涩受| 巨乳人妻的诱惑在线观看| av国产精品久久久久影院| 日韩精品免费视频一区二区三区| 国产又色又爽无遮挡免费看| 亚洲av第一区精品v没综合| 正在播放国产对白刺激| 精品电影一区二区在线| 久久精品人人爽人人爽视色| 久久久精品国产亚洲av高清涩受| 亚洲午夜理论影院| 免费少妇av软件| 在线观看午夜福利视频| 人人澡人人妻人| 日本五十路高清| 俄罗斯特黄特色一大片| 成熟少妇高潮喷水视频| 国产视频一区二区在线看| 黄色视频不卡| 女人被狂操c到高潮| 欧美日韩av久久| 久久中文看片网| 国产欧美日韩综合在线一区二区| 久久国产亚洲av麻豆专区| 日韩制服丝袜自拍偷拍| 99精国产麻豆久久婷婷| 亚洲精品乱久久久久久| 亚洲精品一卡2卡三卡4卡5卡| 热99久久久久精品小说推荐| 成年人黄色毛片网站| 午夜日韩欧美国产| 一级a爱片免费观看的视频| 国产精品偷伦视频观看了| 亚洲在线自拍视频| 亚洲精品在线观看二区| 交换朋友夫妻互换小说| 丰满迷人的少妇在线观看| 婷婷精品国产亚洲av在线 | 9191精品国产免费久久| 免费看a级黄色片| 老鸭窝网址在线观看| 制服人妻中文乱码| 久久久久久久国产电影| 午夜亚洲福利在线播放| 多毛熟女@视频| 亚洲精品中文字幕在线视频| 亚洲av电影在线进入| 久久久国产一区二区| 久久久精品区二区三区| 色94色欧美一区二区| 成年女人毛片免费观看观看9 | 在线免费观看的www视频| 别揉我奶头~嗯~啊~动态视频| 亚洲专区字幕在线| 国产精品国产av在线观看| 麻豆国产av国片精品| 国产成+人综合+亚洲专区| 另类亚洲欧美激情| 99香蕉大伊视频| 看片在线看免费视频| 村上凉子中文字幕在线| 18禁裸乳无遮挡免费网站照片 | av一本久久久久| 成人av一区二区三区在线看| 久久人妻福利社区极品人妻图片| 国产精品 欧美亚洲| 一进一出抽搐gif免费好疼 | 男女下面插进去视频免费观看| 精品国产乱子伦一区二区三区| 亚洲国产精品sss在线观看 | 日本精品一区二区三区蜜桃| 女警被强在线播放| 欧美 亚洲 国产 日韩一| 一本一本久久a久久精品综合妖精| 免费在线观看亚洲国产| a级毛片黄视频| 19禁男女啪啪无遮挡网站| 久久精品国产99精品国产亚洲性色 | 亚洲情色 制服丝袜| videos熟女内射| a在线观看视频网站| av电影中文网址| 亚洲男人天堂网一区| 两人在一起打扑克的视频| 9热在线视频观看99| 国产免费现黄频在线看| 日韩欧美在线二视频 | 男人的好看免费观看在线视频 | 99热网站在线观看| 曰老女人黄片| 日本黄色日本黄色录像| 免费不卡黄色视频| 国产精品亚洲一级av第二区| 亚洲av日韩精品久久久久久密| 少妇的丰满在线观看| 999久久久精品免费观看国产| 男女之事视频高清在线观看| 精品久久久久久久久久免费视频 | 一级黄色大片毛片| 午夜福利在线免费观看网站| 久久午夜综合久久蜜桃| 男女午夜视频在线观看| 亚洲五月婷婷丁香| 亚洲va日本ⅴa欧美va伊人久久| 深夜精品福利| 欧美另类亚洲清纯唯美| 女人被躁到高潮嗷嗷叫费观| av网站免费在线观看视频| 一级毛片精品| 亚洲片人在线观看| 亚洲成a人片在线一区二区| 999久久久精品免费观看国产| 男女之事视频高清在线观看| 黑人欧美特级aaaaaa片| 亚洲国产精品sss在线观看 | 国产一区二区三区综合在线观看| 成人三级做爰电影| 免费在线观看影片大全网站| 久久99一区二区三区| 色播在线永久视频| 国产蜜桃级精品一区二区三区 | 91字幕亚洲| 一边摸一边抽搐一进一小说 | 国产一区二区三区在线臀色熟女 | 精品国产乱子伦一区二区三区| 女人高潮潮喷娇喘18禁视频| 国产av又大| 亚洲色图av天堂| 美女高潮到喷水免费观看| 在线播放国产精品三级| 新久久久久国产一级毛片| 亚洲熟女精品中文字幕| 美国免费a级毛片| 一二三四社区在线视频社区8| 色尼玛亚洲综合影院| 精品人妻1区二区| 桃红色精品国产亚洲av| 亚洲午夜精品一区,二区,三区| 日韩视频一区二区在线观看| 国产精品一区二区在线不卡| 国产在线一区二区三区精| 99riav亚洲国产免费| 女人爽到高潮嗷嗷叫在线视频| 国产欧美日韩精品亚洲av| 国产aⅴ精品一区二区三区波| 男人舔女人的私密视频| 男女免费视频国产| 免费一级毛片在线播放高清视频 | 免费一级毛片在线播放高清视频 | 丁香六月欧美| 国产精品香港三级国产av潘金莲| 美女 人体艺术 gogo| 亚洲精品国产区一区二| 久久亚洲真实| 国产在视频线精品| 视频区欧美日本亚洲| 久久ye,这里只有精品| 三上悠亚av全集在线观看| 亚洲av日韩在线播放| av电影中文网址| 亚洲精品粉嫩美女一区| 久久久久精品国产欧美久久久| tube8黄色片| 99热国产这里只有精品6| 黑人巨大精品欧美一区二区mp4| 性色av乱码一区二区三区2| 国产99久久九九免费精品| 久久精品国产清高在天天线| 精品久久蜜臀av无| 婷婷精品国产亚洲av在线 | 建设人人有责人人尽责人人享有的| 在线看a的网站| 亚洲欧美一区二区三区久久| 在线免费观看的www视频| 一级片'在线观看视频| 十八禁高潮呻吟视频| 久久九九热精品免费| 亚洲精品成人av观看孕妇| 日韩人妻精品一区2区三区| 一本大道久久a久久精品| 久久精品人人爽人人爽视色| 怎么达到女性高潮| 国产91精品成人一区二区三区| 99久久综合精品五月天人人| 亚洲精品在线观看二区| 一级作爱视频免费观看| 久久久精品区二区三区| 亚洲av欧美aⅴ国产| 老司机亚洲免费影院| 精品国产乱码久久久久久男人| 久久 成人 亚洲| 亚洲精品美女久久av网站| 美女扒开内裤让男人捅视频| 91成人精品电影| 三级毛片av免费| 国产欧美日韩一区二区三区在线| 嫁个100分男人电影在线观看| 国产精品1区2区在线观看. | 三级毛片av免费| 三上悠亚av全集在线观看| a在线观看视频网站| 欧美日韩精品网址| 久久99一区二区三区| 中文字幕精品免费在线观看视频| 久久国产精品大桥未久av| 一级毛片高清免费大全| 岛国毛片在线播放| 一级毛片女人18水好多| 日韩欧美三级三区| 久久午夜综合久久蜜桃| 两个人免费观看高清视频| 人妻久久中文字幕网| 啦啦啦免费观看视频1| 成人三级做爰电影| 一级a爱片免费观看的视频| 久久草成人影院| 国产真人三级小视频在线观看| 久久久国产成人精品二区 | 欧美日韩亚洲高清精品| 国产精品成人在线| 亚洲成人免费电影在线观看| 天天操日日干夜夜撸| 成年人午夜在线观看视频| 亚洲aⅴ乱码一区二区在线播放 | 一级作爱视频免费观看| 婷婷成人精品国产| 91国产中文字幕| 精品国产超薄肉色丝袜足j| 男人的好看免费观看在线视频 | 一级片'在线观看视频| 免费少妇av软件| 亚洲综合色网址| 午夜福利欧美成人| 中出人妻视频一区二区| tube8黄色片| 精品欧美一区二区三区在线| 天天影视国产精品| 欧美丝袜亚洲另类 | 美国免费a级毛片| 久久香蕉国产精品| 91精品国产国语对白视频| 国产三级黄色录像| 成人精品一区二区免费| 成在线人永久免费视频| 日韩免费高清中文字幕av| 少妇裸体淫交视频免费看高清 | 叶爱在线成人免费视频播放| 好看av亚洲va欧美ⅴa在| 国内毛片毛片毛片毛片毛片| 国产精品久久久久成人av| 叶爱在线成人免费视频播放| 国产在视频线精品| 欧美在线一区亚洲| 一二三四在线观看免费中文在| av网站免费在线观看视频| 变态另类成人亚洲欧美熟女 | 国产成人系列免费观看| 亚洲第一欧美日韩一区二区三区| 亚洲第一av免费看| 久久久国产一区二区| 国产蜜桃级精品一区二区三区 | 精品人妻在线不人妻| 国内久久婷婷六月综合欲色啪| 宅男免费午夜| 亚洲在线自拍视频| 国产精品.久久久| 精品国产国语对白av| 国内毛片毛片毛片毛片毛片| 欧美精品啪啪一区二区三区| 精品福利永久在线观看| 在线观看66精品国产| 一区二区三区精品91| 12—13女人毛片做爰片一| 美女 人体艺术 gogo| 搡老岳熟女国产| 多毛熟女@视频| www.999成人在线观看| 精品卡一卡二卡四卡免费| 99国产综合亚洲精品| 在线十欧美十亚洲十日本专区| 亚洲av熟女| 麻豆成人av在线观看| 制服诱惑二区| 18禁观看日本| 国产精品秋霞免费鲁丝片| 午夜福利视频在线观看免费| 精品午夜福利视频在线观看一区| 精品人妻在线不人妻| 国产aⅴ精品一区二区三区波| 国产精品二区激情视频| 色综合婷婷激情| а√天堂www在线а√下载 | 啪啪无遮挡十八禁网站| 变态另类成人亚洲欧美熟女 | 国产精品影院久久| 国产单亲对白刺激| 午夜福利一区二区在线看| 国产精品一区二区免费欧美| 女人久久www免费人成看片| 亚洲一区二区三区不卡视频| 啪啪无遮挡十八禁网站| 久久久久久免费高清国产稀缺| 成人黄色视频免费在线看| 伊人久久大香线蕉亚洲五| 亚洲第一欧美日韩一区二区三区| 国产欧美日韩一区二区三| 日日爽夜夜爽网站| 国产亚洲精品久久久久5区| 欧美日韩av久久|