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

    Comparative study on high temperature deformation behavior and processing maps of Mg-4Zn-1RE-0.5Zr alloy with and without in-situ sub-micron sized TiB2 reinforcement

    2023-01-08 10:22:34SahooPanigrahi
    Journal of Magnesium and Alloys 2022年12期

    S.K.Sahoo,S.K.Panigrahi,*

    a Department of Mechanical Engineering,Indian Institute of Technology Madras,Chennai-600036,India

    b Applied Magnesium Research Group,Center of Excellence in Materials and Manufacturing for Futuristic Mobility,Indian Institute of Technology Madras, Chennai-600036,India

    Abstract Mg-4Zn-1RE-0.5Zr(ZE41)Mg alloy is extensively used in the aerospace and automobile industries.In order to improve the applicability and performance,this alloy was engineered with in-situ TiB2 reinforcement to form TiB2/ZE41 composite.The high temperature deformation behavior and manufacturability of the newly developed TiB2/ZE41 composite and the parent ZE41 Mg alloy were studied via establishing constitutive modeling of fl w stress,deformation activation energy and processing map over a temperature range of 250 °C-450 °C and strain rate range of 0.001 s-1-10 s-1.The predicted fl w stress behavior of both materials were found to be well consistent with the experimental values.A significan improvement in activation energy was found in TiB2/ZE41 composite(171.54 kJ/mol)as compared to the ZE41 alloy(148.15 kJ/mol)due to the dispersed strengthening of in-situ TiB2 particles.The processing maps were developed via dynamic material modeling.A wider workability domain and higher peak efficien y(45%)were observed in TiB2/ZE41 composite as compared to ZE41 alloy(41%).The Dynamic recrystallization is found to be the dominating deformation mechanism for both materials;however,particle stimulated nucleation was found to be an additional mode of deformation in TiB2/ZE41 composite.The twinning and stress induced cracks were observed in both the materials at low temperature and high strain rate.A narrow range of instability zone is found in the present TiB2/ZE41 composite among the existing published literature on Mg based composites.The detailed microstructural characterization was carried out in both workability and instability domains to establish the governing deformation mechanisms.? 2022 Chongqing University.Publishing services provided by Elsevier B.V.on behalf of KeAi Communications Co.Ltd.This is an open access article under the CC BY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/)Peer review under responsibility of Chongqing University

    Keywords:ZE41 Mg composite;In-situ TiB2 particles;ZE41 Mg alloy;Constitutive equations;Processing map;Dynamic recrystallization.

    1.Introduction

    Magnesium(Mg)matrix in-situ reinforced composites have attracted enormous interest in automobile and aerospace industries because of their inherent properties like high specifi strength,high specifi stiffness,good thermal stability and high wear resistance[1-4].Among various cast Mg alloys,the rare earth contained Mg alloys have proven as potential materials in several applications[5].ZE41(Mg-4Zn-1RE-0.5Zr;where RE=rare earth)is a commercial rare earth contained Mg alloy used for aircraft gearbox and generator housings,particularly in military helicopters[6].The performance,mechanical properties and tribological properties of ZE41 alloy can be improved significantl by reinforcing in-situ/ex-situ particles into ZE41 matrix.Among the conventional in-situ ceramic particles,Titanium diboride(TiB2)particles exhibit superior physical and mechanical properties for the development of Mg metal matrix composites(MMCs)because of high elastic modulus,high hardness,better thermal stability,and wear resistance[7].Apart from this,TiB2particles act as a grain refine for Mg metal matrix composites.In Mg matrix composites,researchers found that TiB2particles reinforced by the in-situ route exhibited superior mechanical properties than reinforced by the ex-situ process[8-10].Therefore insitu TiB2reinforced ZE41 Mg composite will be a prominent material for automobile and aerospace industries.

    The addition of TiB2particles into soft magnesium alloy matrix often causes metallurgical defects during liquid state processing routes like casting.Generally,as cast particle reinforced Mg matrix composites have coarse grains and hardly homogeneous microstructure can be obtained.A secondary bulk processing is necessary to tackle these problems.The additional secondary bulk plastic forming,such as forging,rolling,and extrusion is required to refin microstructure in order to enhance the mechanical properties.However,limited number of slip systems due to hexagonal close packed(HCP)crystal structure of Mg and the addition of high modulus brittle ceramic reinforcements limits plastic forming ability of Mg composites at room temperatures.Therefore,to improve workability,magnesium alloys and their composites are usually deformed at high temperatures,by which the additional slip systems are activated[11].It is necessary to understand the high temperature deformation behavior of newly developed materials prior to the implementation of any kind of secondary thermo-mechanical/manufacturing processes.Many researchers have used constitutive models to describe the fl w behavior material with strain rate,stress,and deformation temperature[12-16].Processing map is used to evaluate the processing window for any kind of alloys and composites.The dynamic material model(DMM)is the most popular model to develop processing map.Based on the processing map,stable and unstable zones for a newly developed material can be identified

    The approach of processing map has been used by many researchers to fin optimum processing windows for Mg based alloys and composites.The primary hot deformation mechanisms causing the workability and instability in the material are identifie based on the microstructural observation of deformed materials.Many researchers found the workability domains for the processing maps of Mg based composites at deformation temperatures greater than 300 °C.They observed continuous dynamic recrystallization(CDRX)and discontinuous dynamic recrystallization(DDRX)are the hot deformation mechanisms for workability domains based on the microstructural characterization[17,18].During the hot deformation of Mg based composites,the dislocation pile-up generally forms near the particles.This leads to the formation of particle deformation zone,which is favorable for the DRX nucleation.The hot deformation mechanism can be predicted through the activation energy(Q)and the stress exponent(n)through constitutive equations[20].The hot deformation activation energy of the SiCp/AZ91 composite[21]is found to be close to the value for grain boundary diffusion in pure Mg,whereas,with the addition of nano-sized in-situ TiB2particle,a much higher activation energy is found in TiB2/AZ91 composite as compared to that of the unreinforced AZ91 Mg alloy.The deformation mechanism based on the constitutive equation for Mg matrix composites is found to be the climb of dislocation[21,22].

    Though several researchers have studied high temperature deformation behavior of Mg based composites via processing map,till now none of the researchers have established the processing map of ZE41 Mg based composite.In our previous research work,the manufacturing methodology to develop sub-micron sized in-situ TiB2reinforced ZE41 Mg MMCs(TiB2/ZE41)along with the scientifi knowhow was highlighted[23].The development of processing map and constitutive model via hot deformation for a novel in-situ sub-micron sized TiB2/ZE41 composite was carried out in this work for the firs time.Apart from these,none of the researchers have studied the comparison of processing map between any form of in-situ Mg composite vs their parent unreinforced counterpart so far.In the current work a firs attempt was carried out to compare the manufacturability via processing map of an in-situ Mg composite(TiB2/ZE41 composite)vs it’s unreinforced counterpart(ZE41).

    The prime objectives of this present research work are(i)to investigate the high temperature deformation behavior of both ZE41 Mg alloy and in-situ sub-micron sized TiB2/ZE41 composite and study the influenc of in-situ TiB2particles on hot workability,(ii)to fin optimum hot working conditions for manufacturability via processing map,(iii)establishment of constitutive based models to describe the fl w behavior,(iv)to understand the deformation mechanisms via microstructural characterization in both stable and unstable zones of processing maps.

    2.Experimental procedure

    2.1.Fabrication of TiB2/ZE41 composite

    A commercial ZE41 rare earth contained magnesium alloy was taken as a matrix material for the development of insitu composite.Titanium(Ti)and Boron(B)powders were used as the starting alloying elements with a weight fraction of 0.68:0.32 for the development of in-situ TiB2reinforcement.The initial metallic powder system and it’s weight fraction were identifie based on the thermodynamic calculation.The reaction temperature for the formation of TiB2phase was found by the differential thermal analysis(DTA)analysis as 900 °C.The Ti and B powers were mixed in a planetary ball milling and the milling was carried out for an optimized time period of 8 h at a milling speed of 350 rpm to form Ti+B clusters.The preheated Ti+B cluster powders were fille into the custom designed drilled holes of ZE41 Mg ingots slices and were placed inside the crucible of casting unit.The melt temperature was increased to 900 °C and held for 2 h after stirring.The Ti+B powder clusters were fragmented and diffused with each other to form in-situ TiB2in the molten ZE41 matrix.The pouring of the melt was done by a custom designed bottom pouring system in a preheated rectangular mold.The entire fabrication process was carried out under the argon gas environment to avoid any chances of fir and oxidation.The detailed fabrication process and in-situ formation mechanism of the developed TiB2/ZE41 composite were reported in the author’s earlier work[23].Similarly,ZE41 Mg alloy was casted to compare with TiB2/ZE41 composite.The weight fraction of TiB2reinforcement was 10 wt%(i.e.4.23 vol fraction)for the developed composite.The size of in-situ reinforcement was found in the range of 400 nm to 1.5μm.The average TiB2particle size was found to be 765 nm.To eliminate the casting defects like porosity,the developed base alloy and composite were partially warm rolled.The samples were heat-treated at 330 °C for 2 h before hot compression.

    Fig.1.Schematic representation of the hot compression test.

    2.2.Hot compression test

    Hot compression tests of unreinforced ZE41 and in-situ TiB2reinforced ZE41 composite were carried out in a high precision INSTRON 3365,which has an integrated high temperature furnace.The hot compression cylindrical samples were cut from both ZE41 alloy and composite as per the ASTM-E9 standard.To ensure the repeatability of results,three samples in each condition were tested.The hot compression was performed at different deformation temperatures of 250,300,350,400 and 450 °C,and the strain rates of 0.001,0.01,0.1,1 and 10 s-1with a constant true strain of 0.4.The samples were heated to the deformation temperature and kept for 5 min before hot compression to get uniform distribution of temperature throughout the sample.The samples were quenched in cold water immediately after hot compression to preserve the deformed microstructure for further analysis.Fig.1 shows the schematic of hot compression process.

    2.3.Microstructural characterization

    The microstructural analysis of unreinforced ZE41 and TiB2/ZE41 composite were studied by using scanning electron microscopy(SEM),optical microscopy(OM)and transmission electron microscopy(TEM).The samples were cut parallel to the compressive loading axis and mounted before polishing.OM and SEM samples were etched by using picric acid after mechanical polishing with emery sheets of different grades followed by diamond polishing of 1 and 0.5μm.The samples for TEM characterization were mechanical polished till the thickness reach to 80μm with emery sheets of different grades and then ion milled.The TEM was carried out by using TECHNAI electron microscopy at 200 kV.

    3.Constitutive equations and processing map

    3.1.Constitutive equations

    Generally,in hot working processes,constitutive equations based on mathematical models are significan to describe the relationship between processing.There are many constitutive mathematical models used by researchers for the prediction of fl w stress in hot deformation of a material.The Arrhenius model is one of the famous model to describe the thermo mechanical behavior of a material based on the activation energy[24,25].There are three types of power law relationship applied in hot deformation,such as(i)power law for low level stress,(ii)exponential law for high level stress and(iii)hyperbolic-sine law for a wide range of stress[26-29].The equations are as follows:

    where,˙εis the strain rate(s-1);σis peak stress(MPa);Ris the universal gas constant(8.314 Jmol-1K-1);Qis the activation energy(kJ mol-1);Tis the temperature(K);nis the stress exponent.A1,A2,n1,β,αare the material constants.The termαis a stress multiplier and can be expressed asβ/n1[27].The use of above equation in hot deformation depends upon the stress level of the material.Generally,the power law relationship Eq.(1))is suitable when the stress level is low(ασ<0.8),the exponential relationship(Eq.(2))is used in high stress level(ασ>1.2),and the hyperbolic sine relationship(Eq.(3))is applied when there is a wide range of stress[25].Therefore,for metal matrix composites,a hyperbolic sine constitutive equation is used.To calculate the above constants,natural logarithm is taken on both sides of Eqs.(1)to(3)and expressed as follows:

    Fig.2.Microstructures of(a,c)ZE41 alloy;(b,d)TiB2/ZE41 composite;(e)Mg7Zn3RE phase;(g)sub-micron sized in-situ TiB2 reinforcement,and(f,h)EDS results of corresponding selected areas.

    From Eqs.(4)and(5),the valuesn1andβof and can be calculated by the linear regression ln˙εvs lnσand ln˙εvsσof respectively.The value ofαcan be get byβ/n1.The slope of in Eq.(6)is the stress exponentnat a specifi temperature.The expression fornis as follow:

    The stress exponent(n)values represents the governing mechanism for hot deformation,n=2 for grain boundaries sliding,n=3 for the viscous glide of dislocation,n=5 for the climb of dislocation andn=8 for the cross-slip screw dislocation/constant substructure model[30-32].

    The apparent activation energy(Q)required for deformation can be derived from Eq.(6)as follows:wheresrepresents the slope of ln[sinh(ασ)]and 1000/Tat a constant strain rate.

    Besides this,the Zener-Hollomon equation can be used to combine the effect of deformation temperature and strain rate on the fl w behavior of a material[33].Zener-Hollomon parameter(Z)is the Arrhenius function and can be expressed as follows:

    After the natural logarithm of Eq.(9)and combing with Eq.(6);the equation will be as follows:

    From Eq.(10),lnAis the intercept of lnZvs.ln[sinh(ασ)].By substituting the values ofn,α,Q,andA,the constitutive model for a new novel material can be developed.

    3.2.Principle of processing map

    The processing map of a particular material is used to understand the workability at different hot deformation conditions in terms of microstructural mechanisms.The processing map has been used widely by many industries to optimize hot deformation parameters in order to avoid any defects during deformation.The processing map can be developed based on Dynamic material modeling(DMM).It is constructed by the superimpose of instability map over power dissipation efficien y map at various deformation temperatures and strain rates.The work piece dissipate power in terms of heat while material under goes hot deformation.This power dissipation in the material results in microstructural changes like dynamic recrystallization,phase transition,and superplastic fl w,which are reflecte by the power dissipation efficien y.

    According to the DMM theory[34],the dissipated power(P)can be split into two complementary parts,such asGcontent andJco-content and can be expressed as follows:

    The contentGexpresses the power dissipation due to plastic deformation and most of the heat is lost in the form of deformation temperature.In contrast,the co-contentJsignifie power dissipation due to microstructural transformations like dynamic recrystallization,dynamic recovery,phase transformation.

    When deformation temperature is constant,the fl w stress of a material according to power law can be expressed as:

    wherekandmrepresent material constant and strain rate sensitivity respectively.

    The strain rate sensitivity(m)can be calculated in terms ofJandGas follows:

    Eq.(12)can be used in Eq.(11)to obtained the co-contentJas:

    For ideal linear dissipation,the value ofmis equal to one.So,the J co-content reaches to the maximum asJmax=σ˙ε/2.For non-linear dissipation,the power dissipation efficien y can be represented by a dimensionless parameter(η)and is the ratio of power dissipation to the maximum power dissipation via microstructural transformation.The expression for power dissipation efficien y as follows:

    The higher the efficien y is the better the workability of a material.However,in some cases,there is deformation instability such as wedge cracking,cavitation,and fracture during hot working,which needs to be avoided.So,it is necessary to fin out the unstable domain in the processing map via instability criteria.Instability criteria are expressed based on extremum principles of thermodynamic irreversibility as applied to large plastic fl w continuum mechanics.Prasad et al.combine the Ziegler’s fl w theory with DMM to establish the fl w instability during the hot deformation process[35].According to Ziegler,the plastic fl w becomes unstable when:

    whereDrepresents the dissipation function,which is the characteristic of the fl w behavior of the material,Prasad replacedDwith power dissipation co-contentJand the equation becomes as:

    Which implies that

    Fig.3.XRD analysis of unreinforced ZE41 Mg alloy and TiB2/ZE41 composite.

    By substituting the value ofJfrom Eq.(14)into Eq.(18),the fl w instability criterion can be expressed as:

    whereξ(˙ε)is a dimensionless fl w instability parameter developed by Prasad.Negative valuesξ(˙ε)show the f ow instability domain or unsafe domain in the processing map.

    4.Results

    4.1.Initial microstructure

    The microstructures prior to the hot compression test for unreinforced ZE41 and TiB2/ZE41 composite are shown in Fig.2.Both the ZE41 alloy and composite show equiaxed grains with the distribution of ternary Mg7Zn3RE phases shown by the arrow mark.Fig.2a shows the OM image of ZE41 alloy havingα-Mg and ternary Mg7Zn3RE phases.Fig.2c shows the SEM image of the ZE41 alloy.As compared to the ZE41 alloy,the ternary phases in composite material are broken(Fig.2b,14c).The average grain size for ZE41 and TiB2/ZE41 are 41±2.3μm and 26±1.9μm respectively.The presence of TiB2particle cluster in TiB2/ZE41 composite is shown in Fig.2d.The grain size of TiB2/ZE41 composite material is smaller compared to the unreinforced ZE41 alloy.The high magnifie SEM image of Mg7Zn3RE phase is shown in Fig.2e with corresponding EDS result in Fig.2f.The EDS result confirm the presence of cerium(Ce)and neodymium(Nd)as rare earth elements.Fig.2g shows the high magnifica tion SEM image of in-situ TiB2particle and Fig.2h shows the corresponding EDS result.The EDS analysis result confirm the presence of Ti and B in-situ reinforcement as elemental content.The TiB2particle is mostly hexagonal in shape.Fig.3 shows the XRD results of ZE41 Mg alloy and in-situ TiB2/ZE41 composite.The XRD peaks of ZE41 show mainlyα-Mg and the presence of ternary phase i.e.Mg7Zn3RE.The XRD spectra of TiB2/ZE41 composite reveals the presence of TiB2phase includingα-Mg,Mg7Zn3RE and MgO.

    4.2.Flow behavior of ZE41 and TiB2/ZE41

    The hot compression test was carried to understand the fl w behavior of both ZE41 alloy and in-situ TiB2/ZE41 composite at different temperatures and strain rates.Fig.4 and Fig.5 show the true stress-true strain curve of ZE41 alloy and TiB2/ZE41 composite at various temperatures(250-450 °C)and at different strain rates(0.001-10 s-1)respectively.The fl w stress significantl depends on hot deformation parameters(mainly temperature and strain rate).The fl w stresses of both ZE41 alloy and in-situ TiB2/ZE41 composite increased with increasing in strain rate and decrease with increase in deformation temperature.It could be seen from Figs.4 and 5 that,during initial stage of compression i.e.up to a very low true strain,the fl w stress of both alloy and composite increase steadily irrespective of the deformation temperature.This rapid increase in fl w stress is the characteristics of work hardening of the alloy and composite at the initial stage of deformation.The generation and accumulation of dislocation at the earlier stage results in increase of fl w stress in the material.From the fl w behavior of ZE41 and TiB2/ZE41,it can also be observed that,the rate of increase in fl w stress gradually decreases with progress in hot compression i.e.increase in true strain.The ZE41 Mg alloy deformed at 350 °C with strain rate of 0.001 s-1and for TiB2/ZE41 composite deformed at 300 °C with strain rate of 1 s-1show dynamic softening type fl w behavior after work hardening(Fig.4c,Fig.5b).The fl w stress in the composite is found to be higher than the unreinforced ZE41 alloy(Fig.5).This is due to the presence of fin TiB2particles.These fin particles(i)enhance strengthening effect(ii)provide better thermal stability to the material and(iii)also oppose the dislocation movement during deformation.In addition to this,the average grain size of TiB2/ZE41 composite is decreased by 37% than that of unreinforced ZE41 Mg alloy.This is due to the pinning action of TiB2particles,which restricts the grain growth in composite.According to the Hall-Petch relationship,the strength of material increases with decrease in grain size[11].Hence,the increase in fl w stress of TiB2/ZE41 is also due to the decrease in grain size.

    The 3D contour map of true stress vs.temperature vs.strain rate(log scale)is shown in Fig.6 for both ZE41 alloy and TiB2/ZE41 composite.With the increase in strain rate and decrease in deformation temperature,fl w stress increases in both the materials.This shows that,strain rate and deformation temperature have significan impact on the true stress of the material.Similar fl w behavior has also been observed by many researchers for Mg alloys and it’s composites[17,28].

    4.3.Constitutive analysis

    Fig.4.True stress-true strain curves of ZE41 alloy at different temperatures(250 to 450 °C)and various strain rate(0.001 to 10 s-1).

    Fig.5.True stress-true strain curves of TiB2/ZE41 composite at different temperatures(250 to 450 °C)and various strain rate(0.001 to 10 s-1).

    The quantitative relationship between hot deformation conditions and fl w stress can be predicted for both ZE41 Mg alloy and TiB2/ZE41 composite by using the constitutive equations.The compression testing was carried out up to a true strain of 0.4 for both materials and the fl w stresses are selected at a true strain of 0.25 in order to establish the constitutive equations.The constitutive constantsn,β,andαcan be calculated from the above explained constitutive equations.From Eq.(4),n1has been calculated from the linear regression of ln˙εvs.lnσas shown in Fig.7.The value ofn1for ZE41 and TiB2/ZE41 are found to be 8.48 and 11.19 respectively.Similarly,the value ofβis obtained from the slope of ln˙εvs.σas shown in Fig.8 and are found to be 0.065 and 0.074 for ZE41 and TiB2/ZE41 respectively.Then,αfor both alloy and composite are found to be 0.0076 and 0.0067 respectively.From the linear fittin between ln˙εvs.ln[sinh(ασ)](Fig.9)the average stress exponent(n)values for both ZE41 and TiB2/ZE41 are calculated and found to be 5.96 and 8.09 respectively.So,the dominating deformation mechanism for hot compression of ZE41 alloy and TiB2/ZE41 composite are climb of dislocation and cross-slip of screw dislocation/constant substructure model respectively.The average values ofsare calculated from the linear fittin of 1000/Tvs ln[sinh(ασ)](Fig.10).The activation energy of both ZE41 and TiB2/ZE41 are calculated and the quantitative values are represented in Table 1.

    Fig.6.3D contour maps of true stress vs.strain rate(log scale)and temperature for(a)ZE41 alloy and(b)TiB2/ZE41 composite.

    Fig.7.Relationship between ln˙εand lnσfor(a)ZE41 alloy and(b)TiB2/ZE41 composite.

    The activation energy of ZE41 is found to be 148.15 kJ/mol,whereas for TiB2/ZE41 the value is 171.54 kJ/mol.The lattice self-diffusion activation energy of pure magnesium is 135 kJ/mol[36].The ZE41 alloy has higher activation energy of than pure Mg due to the presence of a stable Mg7Zn3RE ternary phase.These ternary phases obstruct the dislocation movement by giving back stress.The activation energy of the composite is higher than the pure Mg and ZE41 alloy.These fin thermally stable in-situ TiB2particles obstruct the dislocation movement during hot deformation,which improves the activation energy of deformation.Apart from this,there is a miss match in the coeffi cient of thermal expansion(CTE)between in-situ TiB2particles and ZE41 matrix,resulting in the generation of dislocations around particles[37]and increase in activation energy.Many researchers also reported the higher activation energy of Mg matrix composites during hot deformation[22,38].

    The lnAvalue was calculated from lnZand ln[sinh(ασ)]and shown in Fig.11.The lnAfor both ZE41 and TiB2/ZE41 are 26.37 and 31.07 respectively.Based on the above constitutive analysis,the constitutive parameters(n,β,andα)are calculated and represented in Table 1.

    Fig.8.Relationship between ln˙εandσfor(a)ZE41 alloy and(b)TiB2/ZE41 composite.

    Fig.9.Relationship between ln˙εand ln[sinh(ασ)]for(a)ZE41 alloy and(b)TiB2/ZE41 composite.

    Table 1Calculated values of constitutive equation constants for both ZE41 and TiB2/ZE41 composite.

    So,the constitutive equation during hot deformation of ZE41 Mg alloy can be expressed as follows:

    Similarly,the constitutive equation of sub-micron sized insitu TiB2reinforced ZE41 composite during hot deformation is expressed as:

    Fig.10.Relationship between 1000/T and ln[sinh(ασ)]for(a)ZE41 alloy and(b)TiB2/ZE41 composite.

    Fig.11.Relationship between lnZ and ln[sinh(ασ)]for both ZE41 alloy and TiB2/ZE41 composite.

    The fl w stress values can be predicted based on the developed Arrhenius based constitutive equations(Eq.(20)and Eq.(21))for both ZE41 Mg alloy and TiB2/ZE41 composite.Fig.12 represents the predicted fl w stress values and the experimental stress values.Most of the stress values for both ZE41 Mg alloy and TiB2/ZE41 composite are near the fittin line of the plot.The R2values for ZE41 and TiB2/ZE41 are 0.967 and 0.959 respectively.These high values of R2show a good correlation of experimental values with the predicted values.This indicates the proposed Arrhenius based constitutive model equations for both materials can well signify the fl w behavior.

    4.4.Processing maps

    The processing map for hot deformation can be developed by superimpose of instability map on power dissipation map.Fig.13 shows the processing maps of unreinforced ZE41 and TiB2/ZE41 composite at a true strain of 0.25.The contour lines represent the iso power dissipation efficien y and the higher the efficien y is,the better the workability of the material will be.It has been reported by many researchers for Mg alloys that,the efficien y greater than 30-35% is suitable for workability of materials[39].The higher values of power dissipation efficien y enhance the nucleation and growth of DRX.So,the workability domains are chosen by considering power dissipation efficien y greater than 33%.The peak effi ciency zones for both ZE41 alloy and TiB2/ZE41 composite are marked with dotted rectangles.The instability or unsafe zones are shown by shaded areas,where the fl w instability parameter is negative(ξ(˙ε)<0)(Fig.13).

    The efficien y values and processing parameters in different zones of processing map for both ZE41 alloy and TiB2/ZE41 composite are shown in Table 2.In the ZE41 alloy,the power dissipation efficien y value decreases with increase in strain rate at low deformation temperature(250-300 °C).Whereas,at moderate and high temperatures with an increase in strain rate,the efficien y value increases.There are three workability domains(Zone-I,II and III)found in processing map of ZE41 alloy.In Zone-I,the efficien y values were found to be 0.33 to 0.41 at a low strain rate(0.001-0.0055 s-1)and a temperature range of 300 °C-375 °C.The efficien y value in Zone-II is found to be 0.33 at a strain rate of 0.1 to 3 s-1(408-434 °C)and in Zone-III the peak efficien y of 0.33-0.39 is observed at a strain rate range of 0.1 to 3 s-1(350-400 °C)(Table 2).The Zone-A represents the instability domain in the processing map of ZE41 alloy(Fig.13a).The instability/unsafe domain is observed at a strain rate ranging from 0.025 to 10 s-1with a temperature range of 250-350 °C(Table.2).Fig.13b shows the processing map for in-situ TiB2/ZE41 composite.Similar to ZE41 alloy,the efficien y value decreases with an increase in strain rate at low temperatures.Two peak efficien y zones(Zone-I and Zone-II)are observed in the case of composite material.In Zone-I,the efficien y value lies in the range of 0.33-0.45 at a stain rate range of 0.001-0.008 s-1(292-384 °C)and in Zone-II,the values range is 0.33-0.39 at strain rate range of 0.0058 to 1 s-1(384-447 °C)(Table.2).These zones in composite material are considered to be the safe workability domain.The instability domains(Zone-A and Zone-B)in the composite are smaller than that of ZE41 alloy.The instability Zone-A is observed at a strain rate of 0.0675-10 s-1(250-300 °C),whereas in Zone B,the instability is found at a high strain rate(6.5-10 s-1)and high temperature(385-450 °C).

    Fig.12.Comparative analysis of the experimental and predicted fl w stress data from the proposed Arrhenius based constitutive model equation for(a)ZE41 alloy and(b)TiB2/ZE41 composite.

    Table 2Different zones in processing maps of both ZE41 and TiB2/ZE41 composite.

    5.Discussion

    The present research work emphasized on the manufacturability of novel ZE41 with and without addition of sub-micron sized in-situ TiB2particles.The hot deformation behavior was studied by developing processing map.In the processing map,each domain is related to different microstructural evolutions and deformation mechanisms.Therefore,for each workability domain and instability domain,microstructural characterization was analyzed in depth for both ZE41 and TiB2/ZE41 composite to understand the hot deformation mechanism discussed in the following section.

    5.1.Workability domain

    The maximum power dissipation efficien y zones in the processing map represent the workability domain for a particular material with certain hot deformation parameters.There are three workability domains(Zone-I,Zone-II and Zone-III)found for ZE41 and two workability domains(Zone-I and Zone-II)in case of TiB2/ZE41 composite(Table 2).Each zone of workability domains is discussed as follows:

    5.1.1.Zone-I

    The Zone-I workability domain possess at medium deformation temperature range(290-385 °C)and at slow strain rate(0.001 s-1)for both the materials(Table.2).However,the Zone-I in composite is wider(almost 1.5 times)than un-reinforced ZE41 alloy.At this workability zone,both base and composite exhibit peak power dissipation efficiencies So,the microstructures at this zone(350 °C and 0.001 s-1)have equiaxed DRX grain structure for both base and composite(Fig.14a,b).Theα-Mg of ZE41 alloy is significantl refine(around 80%)with equiaxed fin grains and average grain size of 8μm(Fig.14c).The ternary Mg7Zn3RE phases are fragmented(Fig.14e).The microstructure of composite also reveals equiaxed fin grains with well-define grain boundaries.The composite microstructure shows the distribution of Mg7Zn3RE phases with some TiB2particle clusters(Fig.14f).The EDS results gives the evidence of presence of phases and particles(Fig.14g,h).These particle clusters are present mainly around the grain boundaries.As compared to the grain size of deformed ZE41 in Zone-I,the average grain size of composite found to be smaller(~6μm)(Fig.14d).The decrease in grain size is mainly due to the in-situ particle pinning effect on grain boundaries which limit the grain growth.The significan refinemen ofα-Mg grains in Zone-I hot deformation condition is primarily due to the fully dynamic recrystallization of material.The DRX is a thermally activated process[41].The extent of DRX depends on temperature and time.In zone-I deformation condition(300-384 °C;0.001-0.008 s-1),the combined action of slow strain rate and moderate temperature helps the material fully recrystallize.For Mg based alloys and composites,the DRX is often predominant due to the activation of dislocation in non-basal slip systems and grain boundary migration at lower strain rate(<0.01 s-1)and relatively higher temperature(>320 °C)[41,42].The fin grain structure with random orientation in zone-I is mainly due to dislocation clime and cross-slip,which provides nuclei for continuous dynamic recrystallization(CDRX).In addition to this,the value of stress component(n)for ZE41 and TiB2/ZE41 are 5.9 and 8.09(Table.1),which indicates the hot deformation is mainly controlled by the climb of dislocation for ZE41 and cross-slip of screw dislocation for TiB2/ZE41.The CDRX is the primary hot deformation mechanism for the zone-I deformation condition.

    Fig.13.Processing maps for(a)ZE41 alloy and(b)TiB2/ZE41 composite.

    The TEM image of ZE41 alloy deformed at Zone-I(350 °C,0.001 s-1)reveals the presence of DRX grains(Fig.15a).The segregation of Ce,Nd rare earth content in term of T-Mg7Zn3RE phases effectively reduce the DRX grains size by solute drag effect during hot deformation[43].Apart from this,Ce,Nd content rare earth phases(Mg7Zn3RE)have better thermal stability and creep properties at this deformation temperature range.Meanwhile,the activation energy during hot deformation has also increased to 148 kJ/mol which is higher than the self-diffusion activation energy of pure Mg(135 kJ/mol).This is due to the presence of rare earth content T phase present in the matrix,which hinder the dislocation motion.The TEM observation shows the dislocation accumulation around the Mg7Zn3RE phase(Fig.16a).It also observed that,there is pile-up of dislocations around this rare earth phase and these phases act as barrier to the dislocation movement(Fig.16a),which helps to promote DRX in ZE41.However,in composite,along with Mg7Zn3RE phases,the presence of sub-micron sized in-situ TiB2particles helps to generate and accumulate more dislocations than ZE41 alloy.Due to the mismatch in the CTE between in-situ TiB2and matrix alloy,a strain fiel is created and hence dislocation clusters are formed around the particles.In addition to this,these in-situ particles act as barrier to dislocation movement,which results in additional pile-up of dislocations in TiB2/ZE41(Fig.16b).The insert image in Fig.16b shows the selected area electron diffraction(SAED)pattern of TiB2particle along zone axis of[0-111].This results in the formation of high dislocation density at the vicinity of particles.These regions are known as particle deformation zone(PDZ)[44]and is significan when the particle size is in the sub-micron or micron level.The DRX and nucleation of new grains are promoted with these high energy stored in PDZ and this mechanism is called as particle stimulated nucleation(PSN)[44].The microstructure reveals the presence of fin DRX grains surrounded by the in-situ TiB2particle(Fig.15b).Under similar deformation conditions(350 °C;0.001s-1)Zhou et al.[45]have observed fin DRX grains for AZ91-5SiC-0.5CNTs composite.They have also found the growth of these fin DRX grains were restricted due to the pinning action of particles.Several researchers have also reported that with the addition of fin reinforcements in Mg metal matrix composite,the nucleation of DRX is also promoted[19,46].

    Fig.14.OM image,grain size distribution plot,SEM image and EDS result of samples deformed at 350 °C;0.001 s-1 for;(a,c,e,g)ZE41 alloy;and(b,d,f,h)TiB2/ZE41 composite.

    Fig.15.TEM microstructures of samples deformed at 350 °C;0.001s-1(a)ZE41 alloy and(b)TiB2/ZE41 composite.

    Fig.16.TEM images of dislocations arrested by(a)Mg7Zn3RE in ZE41 alloy and(b)in-situ TiB2 particles in TiB2/ZE41 composite;(c)TEM images for fragmentation of phase due to particle.

    5.1.2.Zone-II

    Fig.17.OM image,grain size distribution plots,SEM image and EDS results of(a,c,e,g)ZE41 alloy at 425 °C,1 s-1 and(b,d,f,h)TiB2/ZE41 composite at;425 °C,0.1 s-1.

    Another workability domain i.e.Zone-II was found for hot deformation of both ZE41 alloy and TiB2/ZE41 composite at high temperature range(385-447 °C)and medium strain rate(Table.2).As compared to the workability zone-I,the power dissipation efficien y is lower for both materials.The power dissipation efficiencie are 0.33 and 0.33-0.39 for base and composite respectively.But,the zone-II of composite is much wider(almost twice)than base alloy.This is due to the existence of thermally stable in-situ TiB2particles in the material.The microstructure of ZE41 deformed at 425 °C and 1 s-1reveals equiaxed grain and the Mg7Zn3RE phases are mostly present in grain boundaries(Fig.17a,e).However,the average grain size for ZE41 at this zone is 31μm(Fig.17c),which is much higher than that of zone-I,where as in composite it is 20μm(Fig.17d).The grain size is lower in composite because of the presence of in-situ particle.Some particle cluster also present in the composite material(Fig.17b),which is shown in high magnificatio SEM image(Fig.17f).The corresponding EDS results revels the presence of phases and particles(Fig.17g,h).DRX for Mg alloy and its composite is a thermally activated process and the extent of DRX depends on time and temperature.As the deformation progress with higher temperature,the DRX grains grow up faster.The grain size of deformed ZE41 at Zone-II is much larger than that of the sample at Zone-I.At high deformation temperature(>400°C),the DRX grains completely grow with a homogeneous microstructure.The size of grains in Zone-II is increased with an increase in deformation temperature due to the enhancement of mobility of grain boundaries.The main mechanism of DRX is grain boundary migration because the energy supply is high at higher deformation temperature.Similar kind of grain growth characteristics are also observed at deformation temperature range of 400-450 °C and stain rate range of 0.01-1 s-1for different Mg based alloys and composites like AZ31,Mg-5Zn and SiCp/Mg-5Zn composite[46,47].

    5.1.3.Zone-III

    Apart from the two workability zones in both ZE41 and TiB2/ZE41 materials,another peak power dissipation efficien y was found in case of un-reinforced ZE41 alloy(Fig.13a).In comparison to the other two workability domains,this Zone-III workability domain is smaller.The Zone-III workability domain have power dissipation efficien y of 0.33-0.35 and this domain is observed at deformation temperature of 350-400 °C with strain rate range of 5 to 10 s-1(Table 2).Fig.18 shows the typical optical microstructure of ZE41 sample deformed at Zone-III workability domain(375 °C and strain rate of 10 s-1).The ZE41 deformed at this condition was characterized with a mixture of fin dynamic recrystallized and coarse grains.These fin DRX grains are formed as necklace structure(Fig.18).As per the microstructural observation and grain size distribution plot(Fig.18b),a duplex nature of microstructural variation(bimodal grain size)is found in Zone-III i.e.(i)fin necklace type DRX grains(<5μm)with a number fraction of 65% along the grain boundaries and(ii)coarse grain structure.It is also clearly observed that the volume fraction of course grains(>20μm)is much higher than the fin grains.The necklace structure DRX grains are observed as nucleated around the bulging grains.This types of DRX grains are mainly discontinuous dynamic recrystallization(DDRX)grains.Mostly,DDRX are found in the hot deformation of low stacking fault energy materials like Mg[36].The fin DDRX grains were also observed at the vicinity of the rare earth phases(T-Mg7Zn3RE).This is due to more fraction of dislocation accumulation around the T-phases.These rare earth phases are mostly found along the grain boundaries.During hot deformation,the nucleation of new grains starts the areas having high store energy.At high strain rate and high temperature,Mg7Zn3RE had high stored energy than theα-Mg.So,the new grains start nucleating at the vicinity of rare earth phases.With the constant rise in deformation temperature in hot working,more new fin grains are nucleated and form necklace structure.A similar observation was also found for hot deformation of Mg alloys at deformation temperatures greater than 300 °C with a higher strain rate range(>1 s-1)[28,33,48].

    Based on the above study,it is found that,for both ZE41 alloy and TiB2/ZE41 composite,DRX is the dominating mechanism for hot deformation.The optimum workability domain is suggested with high power dissipation efficien y.All the workability domains are at temperature ranges from medium to high domain(300 to 450 °C).In the case of ZE41 alloy,the workability domains are existing at low,medium and high strain rates,whereas in composite,the workability domains are only limited to low and medium strain rates.As compared to the unreinforced ZE41 alloy,the workability domains in the composite are wider.This shows,incorporation of sub-micron sized in-situ TiB2reinforcements expands the workability of ZE41.The DRX phenomenon is a thermally activated process and it happens as long as the activation energy of DRX is satisfied During hot deformation,the material’s internal energy has key role to provide activation energy for DRX.The calculated activation energy of composite is higher than that of ZE41 alloy(Table 1).Apart from the activation energy,the presence TiB2particles promotes the dislocation generation and pile up.This provides additional internal energy required for DRX.Hence TiB2/ZE41 composite possesses better workability domain than unreinforced ZE41 alloy.

    Based on the above microstructural observations,the governing mechanisms of each workability zone for both ZE41 Mg alloy and TiB2/ZE41 Mg composite are identified The continuous dynamic recrystallization(CDRX)is the deformation mechanism for zone-I workability domain.The CDRX is mainly due to the climb of dislocations for ZE41 Mg alloy and cross-slip of screw dislocation for TiB2/ZE41 composite.The PSN effect in composite also helps to promote the DRX in TiB2/ZE41 Mg composite.The hot deformation mechanism in zone-II workability for both materials is mainly due to the grain growth kinetics of DRX grains at high deformation temperature domain(400-450 °C).The zone-III workability of ZE41 Mg alloy forms a typical necklace kind of fin grain structure and this is due to discontinuous dynamic recrystallization(DDRX).

    5.2.Instability domain

    The instability in processing map represents the instantaneous workability of a material.The shaded area in the processing map shows the instability domain of the material.The processing map of un-reinforced ZE41 alloy has one instability domains(Zone-A),however,in case of TiB2/ZE41 composite,there are two instability domains(Zone-A and Zone-B)(Fig.13).The instability of both alloy and composite have occurred mainly at low deformation temperature domain(250 °C-300 °C)and high strain rate(Table.2).The instability area in composite is smaller than that of base alloy.Generally,fl w localization and cracking are the phenomenon for instability in a Mg alloy.Typical optical microstructure in the instability zone for ZE41 sample deformed at 250 °C and 10 s-1is shown in Fig.19a.The microstructure is characterized with the number of mechanical twins,intergranular cracks.The mechanical twins are shown by red arrows while the intergranular cracks with yellow arrows.Fig.19c shows the corresponding SEM image.This shows the presence of micro crack in ZE41 alloy deformed at instability domain.These cracks are generated near to the TMg7Zn3RE phases of the ZE41 alloy.The dislocation movement in Mg alloy can be achieved when the applied stress will be greater than the critical resolved shear stress(CRSS).The required stress for hot deformation is temperature and time-dependent.For Mg alloy,at low temperature,only basal mode of slip systems is activated since the CRSS of basal slip system is lower compare to other slip systems.The additional twins get activated due to lack of slip systems in low-temperature deformation.But,with increase in temperature,the other slip system gets activated,which suppresses the extent of twin formation.So,in low deformation temperature twin mode of deformation is dominated.Apart from this,deformation at a high strain rate increases stress concentration,resulting in the nucleation of voids and cracks.This is because at high strain rate,there is not sufficien time for the material to dissipate heat during hot deformation[49].

    In the case of in-situ TiB2/ZE41 composite,the instability domains are smaller than those of the ZE41 alloy.Fig.19b shows the optical microstructure of TiB2/ZE41 deformed at 250 °C and 1 s-1.The microstructure shows presence of intergranular cracks and mechanical twins in the material.However,intergranular cracks are present like the ZE41 alloy and these cracks are mostly initiated from the particle cluster,which is marked as white ellipse(Fig.19b).Fig.19d shows the SEM image of TiB2/ZE41 sample deformed at instability domain(250 °C,1 s-1).The presence of void and crack around the in-situ TiB2particle clusters are clearly visible.This is due to the dislocation generation and accumulation because of the mismatch in CTE and elastic modulus between the TiB2particle and alloy matrix.During deformation,the generated dislocation may concentrate around the in-situ particles,which further leads the stress concentration at the interface.As the stress concentration exceeds the interface bond strength,intergranular cracking happens.Due to this stress concentration,the interface debonding occurs between the in-situ TiB2particle and the matrix.Fig.19e shows the debonding of a sub-micron size in-situ TiB2particle with the matrix.The corresponding EDS result proofs the evidence of TiB2(Fig.19f).The debondings are the initiation of cracks and voids near the particles.During the deformation of SiCp/Mg composite at deformation temperature of 270 °C and strain rate of 1s-1,similar cracks are observed mainly at the interface between Mg matrix and SiCp[27].The TEM image shows the presence of twins along with high dislocation density present in TiB2/ZE41 composite deformed at instability domain(250 °C and 1s-1)(Fig.20).During deformation at low temperature and high strain rate,twining plays a key role in accommodation of dislocations due to limited number of slip systems.Meanwhile increase in dislocation fraction in composite material leads to an increase in the difficult of plastic deformation and results in fl w instability.So,mechanical twinning is the main cause of fl w instability for the ZE41 Mg alloy.Whereas stress concentration and debonding between TiB2particle and ZE41 matrix interface are causes the fl w instability in composite material.

    Fig.19.Microstructure of sample deformed at 250 °C,10 s-1(a,c)ZE41 alloy;(b,d)TiB2/ZE41 composite and(e,f)debonding between TiB2 and matrix with corresponding EDS result.

    In order to compare our result of instability domains for TiB2/ZE41 composite with the published literature of other Mg matrix composites[18,22,26,27,40,46,50,51],a thorough literature survey was carried out and the results are shown in Fig.21.Most of the researchers conducted hot compression up to a maximum strain rate limit of 1s-1,so the comparison is made within the strain rate range of 0.001s-1to 1s-1.In most of the cases,the instability domains are found in low to medium deformation temperature range(250 to 350 °C).The instability zone for TiB2/ZE41 composite is also found in low temperature and high strain rate as like most of the composites.The fl w instability in most of the composites is due to the micro cracks,voids,wedge cracking and mechanical twinning.They found the cracks and voids at the vicinity of the particles due to the stress concentration.The wider instability domain was found in case of high-volume fraction of reinforcement for the processing map of Mg composites[27,50].Apart from this,micron sized Mg composites exhibit high instability than that of sub-micron/nano size composites.But in few case studies,it was observed that,even though addition of nano reinforcement to Mg alloy matrix,there is increase in fl w instability domain.It was found that,all these composites were developed by ex-situ processing route.So,the size of the reinforcement,volume fraction and processing route play important roles in hot deformation of Mg composites.In this present work,the instability domain of TiB2/ZE41 has become narrow in comparison with other published literature.Because of the novel in-situ processing route and sub-micron sized TiB2particles(average particle size of 765 nm),the developed composite shows wide range of workability domain with narrow instability zone.These in-situ particles provide strong interfacial bonding betweenα-Mg and TiB2.Apart from this,in-situ TiB2particles give better thermal stability than that of ex-situ Mg composite[3].During high temperature deformation of Mg composites,the matrix alloy is getting thermal softening faster specificall at high strain rate than the reinforcements.The matrix alloy undergoes plastic deformation and accumulate stress surrounding the particles.This generates the stress concentration and hence exceeds the critical stress value.If the interface is not strong as in case of in-situ particles,there will be thermal residual stress at the interface.This leads to fl w instability and cracking near the particles.Apart from this,the nano/submicron size particle reinforced Mg composites developed via ex-situ route have more particle clusters than the in-situ Mg composite.This is due to the agglomeration and surface oxidation of starting fin powders,which results in poor wettability with the matrix.During hot deformation of ex-situ nano/sub-micron composites,these particle clusters increase the stress concentration which results cracks and voids.So,the developed sub-micron sized in-situ TiB2/ZE41 composite gives a better choice for manufacturability.

    Fig.20.TEM image of TiB2/ZE41 composite at instability domain.

    Fig.21.Comparison of instability domains of different Mg composites with TiB2/ZE41.

    6.Conclusions

    In this present research work,the hot deformation behavior of both ZE41 and in-situ TiB2/ZE41 composite are studied at different temperatures(250 °C-450 °C)and strain rates(0.001 s-1-10 s-1)based on constitutive analysis and processing maps.The detailed analysis of deformation mechanisms forboth workability and instability are discussed based on the microstructural characterization.The following important conclusions are drawn:

    1)Compared with monolithic ZE41 alloy,the addition of in-situ sub-micron sized TiB2particles effectively enhanced the fl w stress of TiB2/ZE41 composite by hindering the strain-induced dislocations.

    2)The apparent activation energy of TiB2/ZE41 composite is increased from 148.15 kJ/mol to 171.54 kJ/mol due to the effective dispersion strengthening action of in-situ TiB2particles.

    3)The developed constitutive model equations for both ZE41 Mg alloy and TiB2/ZE41 composite are as follows:

    The predicted fl w stress values from the constitutive equations of both materials are correlated well with the experimental results.

    4)According to the developed processing maps,the workability domains at peak efficiencie are found to be at 300-375 °C,0.001-0.0055 s-1for ZE41 Mg alloy and at 292-384 °C,0.001-0.008 s-1for in-situ TiB2/ZE41 composite.The DDRX is the primary governing mechanism in the workability domain of both the materials.The fl w instability for both the materials occurs at a low temperature with a high strain rate,and stressinduced cracking,twinning are the main characteristics of fl w instability.

    5)The peak efficien y of TiB2/ZE41 composite(45%)is higher than that of ZE41 alloy(41%).The processing map of TiB2/ZE41 composite results in wider workability domains and higher peak efficien y(45%)as compared to ZE41 Mg alloy.The presence of these fin in-situ TiB2particles promotes the DRX of the composite due to the additional PSN effect in addition to DDRX.Hence,the deformation of TiB2/ZE41 composite is more stable and homogenous with a wide range of manufacturability.

    6)The presently developed TiB2/ZE41 composite showed a smaller range of instability domain as compared to that of existing Mg based composites published in the literature till date.The higher interfacial bonding between ZE41 matrix and in-situ TiB2reinforcement due to in-situ reaction prevent earlier debonding during high temperature deformation.

    Acknowledgment

    The authors show gratitude to Department of Science and Technology,India[grant number of DST/TDT/AMT/2017/211(G)(MEE/18-19/412/DSTX/SUSH)for the financia support and FIST grant,Department of Science and Technology,India[grant number SR/FST/ET11-059/2012(G)]for funding electron microscope facility.This research work is a part of Center of Excellence(CoE)in Applied Magnesium Research(A Vertical of Center for Materials and Manufacturing for Futuristic Mobility),IIT Madras.The authors acknowledge the Ministry of Human Resource and Development for funding this CoE through grant number-SB20210992MEMHRD008517.

    亚洲最大成人中文| 亚洲内射少妇av| 久久人人爽人人爽人人片va| 97热精品久久久久久| 亚洲成a人片在线一区二区| 级片在线观看| 三级经典国产精品| 亚洲美女黄片视频| 欧美性猛交╳xxx乱大交人| 国产欧美日韩一区二区精品| 日本一本二区三区精品| 1000部很黄的大片| 国产成年人精品一区二区| 18禁在线无遮挡免费观看视频 | 亚洲乱码一区二区免费版| 亚洲内射少妇av| 99久国产av精品国产电影| av黄色大香蕉| av.在线天堂| 色综合亚洲欧美另类图片| av视频在线观看入口| 久久99热6这里只有精品| 精品人妻一区二区三区麻豆 | 国产在线男女| 99久久九九国产精品国产免费| 亚洲乱码一区二区免费版| or卡值多少钱| 午夜激情福利司机影院| 日韩成人伦理影院| 精品人妻偷拍中文字幕| 国产国拍精品亚洲av在线观看| 亚洲欧美清纯卡通| 别揉我奶头 嗯啊视频| 天天一区二区日本电影三级| 国产一区二区在线观看日韩| 观看美女的网站| 国产高清有码在线观看视频| 成年女人看的毛片在线观看| 亚洲成人av在线免费| 2021天堂中文幕一二区在线观| 亚洲欧美精品自产自拍| 悠悠久久av| 三级国产精品欧美在线观看| 欧美+亚洲+日韩+国产| 亚洲激情五月婷婷啪啪| 国产不卡一卡二| 亚洲一级一片aⅴ在线观看| 狂野欧美激情性xxxx在线观看| 精品久久国产蜜桃| 一个人看的www免费观看视频| 日韩在线高清观看一区二区三区| 亚洲丝袜综合中文字幕| 国产精品久久电影中文字幕| av黄色大香蕉| 一级毛片aaaaaa免费看小| 亚洲最大成人中文| 国产一级毛片七仙女欲春2| 日韩欧美一区二区三区在线观看| 亚洲精品久久国产高清桃花| 成人一区二区视频在线观看| 亚洲精品一区av在线观看| 俄罗斯特黄特色一大片| 内地一区二区视频在线| 精品一区二区三区视频在线观看免费| 国国产精品蜜臀av免费| 女人被狂操c到高潮| 国产伦一二天堂av在线观看| 女同久久另类99精品国产91| 中文在线观看免费www的网站| 嫩草影院精品99| 久久久久久大精品| 日日干狠狠操夜夜爽| 男人和女人高潮做爰伦理| 一级黄片播放器| 久久久久久久久久久丰满| 免费av不卡在线播放| 日韩制服骚丝袜av| 精品欧美国产一区二区三| av卡一久久| 女的被弄到高潮叫床怎么办| 国产毛片a区久久久久| avwww免费| av卡一久久| 天堂影院成人在线观看| 99久久久亚洲精品蜜臀av| 国产av不卡久久| 国产人妻一区二区三区在| 亚洲人成网站在线观看播放| 男人舔女人下体高潮全视频| 99久国产av精品国产电影| 联通29元200g的流量卡| 三级经典国产精品| 婷婷精品国产亚洲av在线| 91av网一区二区| 少妇丰满av| 亚洲精品粉嫩美女一区| 嫩草影院新地址| 最近手机中文字幕大全| 精品无人区乱码1区二区| 校园春色视频在线观看| 免费无遮挡裸体视频| 亚洲av.av天堂| 亚洲最大成人手机在线| 少妇的逼水好多| 日产精品乱码卡一卡2卡三| 超碰av人人做人人爽久久| 如何舔出高潮| 小说图片视频综合网站| 一级毛片电影观看 | 久久这里只有精品中国| 国产片特级美女逼逼视频| 欧美潮喷喷水| 欧美日韩精品成人综合77777| 国产精品久久久久久av不卡| 老女人水多毛片| 亚洲国产高清在线一区二区三| 日韩欧美三级三区| 亚洲成人久久爱视频| 色5月婷婷丁香| av国产免费在线观看| 久久久精品欧美日韩精品| 最近的中文字幕免费完整| 成人美女网站在线观看视频| 久久久久国产精品人妻aⅴ院| 日韩欧美 国产精品| 日本色播在线视频| 成人特级av手机在线观看| 国产精品乱码一区二三区的特点| 在线观看66精品国产| 日本与韩国留学比较| 国产高清视频在线观看网站| 99国产精品一区二区蜜桃av| 亚洲电影在线观看av| av天堂中文字幕网| 深夜a级毛片| 成人av一区二区三区在线看| 久久久久免费精品人妻一区二区| 欧美激情国产日韩精品一区| 欧美激情久久久久久爽电影| 1000部很黄的大片| 99久久中文字幕三级久久日本| 舔av片在线| 可以在线观看的亚洲视频| 观看免费一级毛片| 久久久a久久爽久久v久久| 在线播放无遮挡| 插逼视频在线观看| 自拍偷自拍亚洲精品老妇| 在线免费观看的www视频| 日韩成人av中文字幕在线观看 | 黄色日韩在线| 久久久精品94久久精品| 99热这里只有精品一区| 最近手机中文字幕大全| 精品一区二区免费观看| 中文字幕av在线有码专区| 色5月婷婷丁香| 亚洲人成网站在线播放欧美日韩| 2021天堂中文幕一二区在线观| 免费大片18禁| 久久精品国产鲁丝片午夜精品| 亚洲欧美精品自产自拍| 最近中文字幕高清免费大全6| 别揉我奶头 嗯啊视频| 99久久久亚洲精品蜜臀av| 一级毛片久久久久久久久女| 亚洲在线自拍视频| 3wmmmm亚洲av在线观看| 久99久视频精品免费| 国产成人freesex在线 | 日韩成人伦理影院| 精品午夜福利在线看| 久久人人爽人人片av| 色综合色国产| www日本黄色视频网| 菩萨蛮人人尽说江南好唐韦庄 | 村上凉子中文字幕在线| 亚洲精品粉嫩美女一区| 六月丁香七月| 国产精品久久久久久久久免| 欧美日本视频| 亚洲综合色惰| 又爽又黄无遮挡网站| 欧美区成人在线视频| 亚洲无线在线观看| 国产在线精品亚洲第一网站| 精品久久久久久久久亚洲| 男插女下体视频免费在线播放| 欧美潮喷喷水| av免费在线看不卡| 欧美丝袜亚洲另类| 高清午夜精品一区二区三区 | 秋霞在线观看毛片| 亚洲欧美精品自产自拍| 一区福利在线观看| 桃色一区二区三区在线观看| av在线蜜桃| 91在线精品国自产拍蜜月| 成人美女网站在线观看视频| 97人妻精品一区二区三区麻豆| 春色校园在线视频观看| 亚洲精华国产精华液的使用体验 | 亚洲av第一区精品v没综合| 色视频www国产| 天堂网av新在线| 免费无遮挡裸体视频| 欧美一区二区精品小视频在线| 亚洲欧美成人综合另类久久久 | 精品熟女少妇av免费看| 在线免费观看的www视频| 人妻少妇偷人精品九色| 精品不卡国产一区二区三区| 日本爱情动作片www.在线观看 | av在线天堂中文字幕| 天堂影院成人在线观看| 色噜噜av男人的天堂激情| 亚洲国产精品成人久久小说 | 亚洲美女视频黄频| 亚洲成a人片在线一区二区| 全区人妻精品视频| 一区二区三区四区激情视频 | 嫩草影视91久久| 国产精品久久久久久av不卡| 内射极品少妇av片p| 免费不卡的大黄色大毛片视频在线观看 | 成年女人毛片免费观看观看9| 永久网站在线| 欧美成人一区二区免费高清观看| 日本熟妇午夜| 舔av片在线| 男女那种视频在线观看| 男女下面进入的视频免费午夜| 国产精品一区www在线观看| 欧美不卡视频在线免费观看| 免费看美女性在线毛片视频| 又黄又爽又刺激的免费视频.| 美女大奶头视频| 国产欧美日韩一区二区精品| 插逼视频在线观看| av卡一久久| 欧美日本视频| 久久久精品94久久精品| 亚洲三级黄色毛片| 亚洲欧美成人综合另类久久久 | 久久久久国内视频| 午夜久久久久精精品| 国产亚洲精品久久久久久毛片| 五月玫瑰六月丁香| 欧美高清性xxxxhd video| 日韩一本色道免费dvd| 亚洲精品久久国产高清桃花| 成年版毛片免费区| 日本a在线网址| 亚洲丝袜综合中文字幕| 亚洲18禁久久av| 99热这里只有精品一区| 国产 一区精品| 麻豆国产av国片精品| 成年版毛片免费区| 亚洲七黄色美女视频| 国产成人a区在线观看| 久久热精品热| 国产精品精品国产色婷婷| 久久久久性生活片| 中文字幕av在线有码专区| 人人妻,人人澡人人爽秒播| 日韩强制内射视频| aaaaa片日本免费| 欧美+日韩+精品| 欧美bdsm另类| 俺也久久电影网| 久久久久免费精品人妻一区二区| 天堂√8在线中文| 亚洲不卡免费看| 此物有八面人人有两片| 最近在线观看免费完整版| 国语自产精品视频在线第100页| 嫩草影院精品99| 老师上课跳d突然被开到最大视频| 国产高潮美女av| 搡老熟女国产l中国老女人| 亚洲国产高清在线一区二区三| 亚洲精品成人久久久久久| 免费看美女性在线毛片视频| 久久久国产成人免费| 国产aⅴ精品一区二区三区波| 国产精品久久久久久久电影| АⅤ资源中文在线天堂| 干丝袜人妻中文字幕| 亚洲中文字幕日韩| 日本免费a在线| 午夜久久久久精精品| 午夜日韩欧美国产| 精品久久久久久久人妻蜜臀av| 亚洲av中文av极速乱| 淫秽高清视频在线观看| 熟女人妻精品中文字幕| 中文字幕av成人在线电影| 免费电影在线观看免费观看| 国产亚洲精品久久久com| 亚洲丝袜综合中文字幕| 国产成人91sexporn| 亚洲欧美精品综合久久99| 欧美又色又爽又黄视频| 午夜爱爱视频在线播放| 免费在线观看影片大全网站| 亚洲四区av| 免费av不卡在线播放| а√天堂www在线а√下载| 国产精品不卡视频一区二区| 日韩欧美一区二区三区在线观看| 嫩草影院精品99| 久久精品人妻少妇| a级毛片a级免费在线| 成年免费大片在线观看| 大香蕉久久网| 久久中文看片网| 一级黄片播放器| 99久久久亚洲精品蜜臀av| 久久久久国内视频| 99视频精品全部免费 在线| 女人十人毛片免费观看3o分钟| 久久亚洲国产成人精品v| 女生性感内裤真人,穿戴方法视频| 亚洲在线观看片| 香蕉av资源在线| 色尼玛亚洲综合影院| 一a级毛片在线观看| 免费看av在线观看网站| 亚洲国产精品sss在线观看| 少妇熟女欧美另类| 99国产极品粉嫩在线观看| 亚洲人与动物交配视频| 天天躁夜夜躁狠狠久久av| 麻豆国产av国片精品| 国产av在哪里看| 亚洲人与动物交配视频| 久久欧美精品欧美久久欧美| 秋霞在线观看毛片| 自拍偷自拍亚洲精品老妇| 国产片特级美女逼逼视频| 亚洲国产高清在线一区二区三| 中国美女看黄片| 亚洲精品一卡2卡三卡4卡5卡| 日日摸夜夜添夜夜爱| 久久久久久久久久成人| 狠狠狠狠99中文字幕| 国产探花极品一区二区| 日本爱情动作片www.在线观看 | 国产一区二区在线av高清观看| 国产 一区 欧美 日韩| 草草在线视频免费看| 91精品国产九色| 人人妻人人澡欧美一区二区| 国产单亲对白刺激| 少妇被粗大猛烈的视频| 亚洲图色成人| 国产精品一及| 亚洲欧美日韩高清在线视频| 女生性感内裤真人,穿戴方法视频| 欧美潮喷喷水| 日韩强制内射视频| 国产伦在线观看视频一区| 激情 狠狠 欧美| 亚洲丝袜综合中文字幕| 91久久精品国产一区二区三区| 免费一级毛片在线播放高清视频| a级毛片免费高清观看在线播放| 可以在线观看毛片的网站| 久久人人爽人人爽人人片va| 赤兔流量卡办理| a级毛片免费高清观看在线播放| 午夜视频国产福利| 偷拍熟女少妇极品色| 欧美+日韩+精品| 一夜夜www| 一个人看视频在线观看www免费| 免费看美女性在线毛片视频| 又黄又爽又免费观看的视频| 亚洲aⅴ乱码一区二区在线播放| 伦理电影大哥的女人| 男女视频在线观看网站免费| 在线a可以看的网站| 性色avwww在线观看| 高清午夜精品一区二区三区 | 国内精品久久久久精免费| 国产成年人精品一区二区| 精品久久久久久久久久免费视频| 精品久久国产蜜桃| 久久亚洲国产成人精品v| 激情 狠狠 欧美| 欧美日韩国产亚洲二区| 狠狠狠狠99中文字幕| 少妇高潮的动态图| 插阴视频在线观看视频| 国内精品美女久久久久久| 麻豆国产av国片精品| 国内少妇人妻偷人精品xxx网站| 亚洲人与动物交配视频| 色在线成人网| 老司机福利观看| 欧美bdsm另类| 国产aⅴ精品一区二区三区波| 深夜精品福利| 赤兔流量卡办理| 黄色视频,在线免费观看| 欧美日本视频| 日日摸夜夜添夜夜爱| 亚洲经典国产精华液单| 午夜福利视频1000在线观看| 99久久久亚洲精品蜜臀av| 最近2019中文字幕mv第一页| 亚洲成人精品中文字幕电影| 亚洲美女搞黄在线观看 | 波多野结衣巨乳人妻| 噜噜噜噜噜久久久久久91| 中国美白少妇内射xxxbb| 一夜夜www| 看十八女毛片水多多多| 国产亚洲av嫩草精品影院| 寂寞人妻少妇视频99o| 别揉我奶头~嗯~啊~动态视频| 亚洲精品亚洲一区二区| 国产v大片淫在线免费观看| 久久精品夜夜夜夜夜久久蜜豆| 在线播放国产精品三级| 国产私拍福利视频在线观看| 老女人水多毛片| 哪里可以看免费的av片| 日日摸夜夜添夜夜添av毛片| 卡戴珊不雅视频在线播放| av福利片在线观看| 麻豆国产av国片精品| 亚洲精品日韩在线中文字幕 | 国产精品乱码一区二三区的特点| 97热精品久久久久久| 永久网站在线| 麻豆乱淫一区二区| 午夜a级毛片| 久久久a久久爽久久v久久| 国产精品,欧美在线| 亚洲图色成人| 看非洲黑人一级黄片| 日本撒尿小便嘘嘘汇集6| 国产女主播在线喷水免费视频网站 | 色综合亚洲欧美另类图片| 夜夜夜夜夜久久久久| 男女啪啪激烈高潮av片| 亚洲最大成人手机在线| 成年版毛片免费区| 狠狠狠狠99中文字幕| 国产午夜精品论理片| 亚洲aⅴ乱码一区二区在线播放| 九九在线视频观看精品| 简卡轻食公司| 国产一区二区三区av在线 | 免费av不卡在线播放| 亚洲欧美日韩高清专用| av中文乱码字幕在线| 成人亚洲精品av一区二区| 国产精品,欧美在线| 免费一级毛片在线播放高清视频| 色尼玛亚洲综合影院| 一区福利在线观看| 午夜精品一区二区三区免费看| 国国产精品蜜臀av免费| 日本熟妇午夜| 国产精华一区二区三区| 日产精品乱码卡一卡2卡三| 亚洲av不卡在线观看| 国产精品久久久久久精品电影| 亚洲在线观看片| 久久久精品欧美日韩精品| 日韩精品中文字幕看吧| 在现免费观看毛片| 全区人妻精品视频| 亚洲精品久久国产高清桃花| 伦理电影大哥的女人| 少妇裸体淫交视频免费看高清| 一级毛片aaaaaa免费看小| 免费不卡的大黄色大毛片视频在线观看 | 嫩草影视91久久| 久久99热6这里只有精品| 久久精品人妻少妇| 久久久国产成人免费| 亚洲成人av在线免费| 日本免费a在线| 一个人看视频在线观看www免费| 97热精品久久久久久| 精品人妻偷拍中文字幕| 色av中文字幕| 哪里可以看免费的av片| 日本熟妇午夜| 一本一本综合久久| 最好的美女福利视频网| 直男gayav资源| 国产精品三级大全| 国产精品国产高清国产av| 亚洲四区av| 国产精品福利在线免费观看| 久久久久久久久久久丰满| 久久精品久久久久久噜噜老黄 | 久久精品夜夜夜夜夜久久蜜豆| 国产69精品久久久久777片| .国产精品久久| 狠狠狠狠99中文字幕| 亚洲欧美日韩东京热| 亚洲av第一区精品v没综合| 亚洲,欧美,日韩| 久久综合国产亚洲精品| 亚洲,欧美,日韩| 校园人妻丝袜中文字幕| 国产成人一区二区在线| 欧美性猛交╳xxx乱大交人| 波多野结衣巨乳人妻| 一区二区三区四区激情视频 | 十八禁网站免费在线| 久久久久久久久久成人| 禁无遮挡网站| 国产一区二区三区av在线 | 悠悠久久av| 99久国产av精品| 天堂动漫精品| 少妇熟女aⅴ在线视频| 精品欧美国产一区二区三| 美女大奶头视频| 成人美女网站在线观看视频| 精品免费久久久久久久清纯| 国产午夜福利久久久久久| 免费看av在线观看网站| 亚洲自偷自拍三级| 日本熟妇午夜| 亚洲国产色片| 色噜噜av男人的天堂激情| 一级毛片aaaaaa免费看小| 欧美激情在线99| 国产爱豆传媒在线观看| 亚洲精品乱码久久久v下载方式| 国国产精品蜜臀av免费| 免费av不卡在线播放| 99热6这里只有精品| 国产伦精品一区二区三区视频9| 美女被艹到高潮喷水动态| 国产高清有码在线观看视频| 91久久精品电影网| 校园春色视频在线观看| 色综合站精品国产| 波野结衣二区三区在线| 99热全是精品| 天美传媒精品一区二区| 不卡一级毛片| 午夜精品国产一区二区电影 | 卡戴珊不雅视频在线播放| 国内少妇人妻偷人精品xxx网站| 亚洲美女搞黄在线观看 | 国产黄片美女视频| 最近的中文字幕免费完整| www.色视频.com| 尾随美女入室| 国产亚洲av嫩草精品影院| 91av网一区二区| 精品国内亚洲2022精品成人| 一个人观看的视频www高清免费观看| 日韩一本色道免费dvd| 国内久久婷婷六月综合欲色啪| 亚洲最大成人av| 小说图片视频综合网站| 干丝袜人妻中文字幕| 国产欧美日韩精品亚洲av| 国产精品国产高清国产av| 国产麻豆成人av免费视频| 亚洲综合色惰| 国产免费男女视频| 亚洲国产色片| 国产蜜桃级精品一区二区三区| 亚洲国产日韩欧美精品在线观看| 国产亚洲av嫩草精品影院| 3wmmmm亚洲av在线观看| 小蜜桃在线观看免费完整版高清| 日韩欧美精品免费久久| 精品人妻视频免费看| 啦啦啦韩国在线观看视频| 男女视频在线观看网站免费| 国产av在哪里看| 久久久精品大字幕| 免费搜索国产男女视频| 欧美又色又爽又黄视频| 精品一区二区三区av网在线观看| 69av精品久久久久久| 波多野结衣高清无吗| 国产精品一区二区免费欧美| 欧美zozozo另类| 日日啪夜夜撸| 久久久久久久亚洲中文字幕| 欧美在线一区亚洲| 蜜臀久久99精品久久宅男| 国产在线精品亚洲第一网站| 日韩av在线大香蕉| 91午夜精品亚洲一区二区三区| 国产aⅴ精品一区二区三区波| 狠狠狠狠99中文字幕| 日韩精品青青久久久久久| 又黄又爽又刺激的免费视频.| 欧美zozozo另类| 日产精品乱码卡一卡2卡三| 欧美一级a爱片免费观看看| 男女视频在线观看网站免费| 国产精品福利在线免费观看| 欧美最黄视频在线播放免费| 国模一区二区三区四区视频| 国产 一区精品| 少妇的逼水好多| 久久午夜亚洲精品久久| 久久精品国产亚洲av香蕉五月| 特级一级黄色大片|