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

    Improved mechanical properties and strengthening mechanism with the altered precipitate orientation in magnesium alloys

    2022-07-13 03:09:50WnZengDouHuQinZengSunQun
    Journal of Magnesium and Alloys 2022年5期

    Y.J.Wn, Y.Zeng,?, Y.C.Dou, D.C.Hu, X.Y.Qin, Q.Zeng, K.X.Sun, G.F.Qun

    a Key Laboratory of Advanced Technologies of Materials, Ministry of Education, School of Material Science and Engineering, Southwest Jiaotong University,Chengdu 610031, China

    b College of Materials Science and Engineering, Sichuan University of Science and Engineering, Zigong 643000, China

    c Department of Materials, The University of Manchester, Oxford Road, Manchester M13 9PL, United Kingdom

    Abstract

    Keywords: Precipitate orientation; Precipitation strengthening; Grain boundary strengthening; Yield strength model.

    1.Introduction

    Magnesium (Mg) alloys, served as the lightest metallic structural materials, have great potential for improving the fuel efficiency and reducing carbon dioxide (CO2) emission of vehicles [1].Mg-3wt.%Al-1wt.%Zn (AZ31) alloy, a common commercial Mg alloy, is widely utilized in the field of transportation industry and aerospace applications due to the high ductility and good formability [2].However, the extensive applications of AZ31 alloy are restricted due to its low strength.

    In general, precipitation strengthening is considered as a feasible method to improve the strength of Mg alloys [3].However, AZ31 Mg alloy is acknowledged as a non-heattreatable Mg alloy [4].During the process of aging treatment,β-Mg17Al12particles would precipitate in Mg matrix.The most commonly observed orientation relationship (OR)between Mg17Al12precipitates and Mg matrix is the Burgers OR,i.e., (011)β//(0001)α, [11]β//[20]α[5–7].It is noteworthy that the majority of Mg17Al12precipitates in Mg-Al alloys have habit planes (HPs) parallel or approximately parallel to the(0002)basal planes of Mg matrix(hereafter named basal plates) [7–9].Such basal plates do not raise appreciable strengthening due to its limited impediment of the dominant slip system in Mg alloys at room temperature (RT),i.e., basal slip [4].Thus, the strengthening effect of AZ31 Mg alloy is limited by the conventional aging treatment.So, it is urgentto find an effective way to improve the strength of AZ31 Mg alloy.

    In terms of precipitation strengthening, the orientation of precipitate with Mg matrix plays an important role in the hardening effect [9,10].Nie [11]proposed a physical model to quantitatively investigate the effect of precipitate shape and orientation on dispersion strengthening based on Orowan strengthening theory for hexagonal-closed packed (HCP)metals and thus predicated that prismatic plates have a strong blocking effect on basal slip.As a result, the basal plates formed in the as-cast AZ91 alloy generally give poor strengthening compared to the precipitates formed in the prismatic plane of the as-cast Mg-Y-RE alloys (hereafter named prismatic plates) [4].However, such prismatic plates were predicted to have a weak blocking effect on prismatic slip which is weaker than basal plates.In addition, Robson et al.[12]also predicted that rod-shaped precipitates with their long-axes parallel to thec-axes of Mg crystals (hereafter namedc-axis rods) are more effective in hardening against basal slip and prismatic slip compared to sphere and basal plates.Although the Orowan model can well estimate the hardening effects on various slip modes, it is not suitable for twinning deformation because the Orowan model was found to underestimate grately the hardening effect on twinning.Recently, Fan et al.[10,13]developed a novel physical model to study the precipitation hardening effect of different types of precipitates (basal plates, prismatic plates,c-axis rods and spherical precipitates) on twinning deformation.It was found that the prismatic plates have the best hardening effect on both twinning and basal slip.From the aforementioned researches,it can be concluded that different types of precipitates lead to different strengthening responses.In this sense, precipitation is considered as an important tool for controlling the strength of individual slip and twin modes,thereby adjusting the asymmetry and improving the overall strength and ductility of Mg alloys [10].Inspired by the provious researches, the strength of AZ31 Mg alloy could be improved by regulating the orientation between precipitate and Mg matrix (hereafter named precipitate orientation).

    For HCP magnesium and its alloys, deformation twinning plays an important role during plastic deformation process [14].Two common twinning modes were observed in Mg alloys.They were classified as {102} entension twins and {101} contraction twins since they result in an extension or contraction alongc-axis, respectively [14].In general,twinning can change the grain orientation of the Mg alloys[15,16].For extension twins, it reorients the lattice by ~86.3°along 〈11ˉ20〉 crystal directions [15].Moreover, {102} extension twins can be easily activated during compression along rolling direction (RD) or transverse direction (TD) or tension along normal direction(ND)in strong basal texture Mg alloys sheet [17].Recently, the interactions of precipitates and twins were reported in Refs.[10,18,19].These researches indicated that the precipitates in Mg alloys might increase the stress for twin growth or detwinning but does not prohibit them, and twinning generally does not shear precipitates but engulf them and lead only to a rigid body rotation of the precipitates by 3–4° with repect to its original position.That is, during the process of {102} extension twinning, the lattice for matrix would be reoriented by ~86.3°, while the precipitates would close to their original position.

    Therefore, it is possible to tailor the precipitate orientation through controlling twinning and precipitation process.More recently, Liu et al.[8]reported a method for regulating the precipitate orientation in AZ80 alloys by coupling twinning, aging and detwinning processes (TAD processes).The results of microstructural observations indicated that the precipitate orientation was altered from basal precipitation to prismatic precipitation.Compression and tensile tests were performed to confirm that such change in precipitate orientation can remarkably improve the strength of as-rolled AZ80 Mg alloy sheet.Furthermore, a crystallography-based algorithm was proposed to predict the possible crystallographic orientation relationships (ORs) between the prismatic precipitates and the detwinned matrix, and the predictions were subsequently validated through experimental observation[20].The results showed that three new ORs betweenβ-Mg17Al12andα-Mg matrix were determined through experimental observation, which agreed well with the predictions through the algorithm after considering the 3.69° rotation of precipitates at twin boundary.It could be further deduced that the precipitates is not exactly parallel to prismatic plane because the lattice ortation angle is 86.3° and the precipitates still have a small rotation at twin boundary [20].Although this TAD process can regulate precipitate orientation and thus improve the overall strength, it seems to be cumbersome and the underlying strengthening mechanism was not discussed in depth.

    Therefore, in the present work, an aging prior to twinning deformation process was applied to improve the mechanical properties of AZ31 alloy through altering the precipitate orientation.The mechanical properties between the aged sheet and the twinned sheet were compared to assess the strengthening effect of the altered precipitate orientation and twin boundary.Specifically, all samples used for testing mechanical properties were conducted along 45°to ND(in the ND-TD plane) to avoid the influence of texture on properties.Furthermore, the inherent mechanism of the improved strength and ductility in the sample containing the prismatic plates was analyzed in detail.

    2.Materials and methods

    The material used in this study is a commercial hot-rolled AZ31 Mg alloy sheet with a thickness of 10mm and a typical strong basal texture, as shown in Fig.1b.The experimental process is illustrated in Fig.1a, and described as follows.Firstly, to obtain uniformed precipitates, the as-received sheet was solution treated at 410°C for 8h, which still possessed a typical strong basal texture with a slightly decreased texture intensity, as shown in Fig.1c.Then, this sheet was aged at 180°C for 14h (T6 sheet).Cubic block samples with 10mm×10mm×15mm along ND, TD and RD were machined from the aged sheet (hereafter named sample A).Secondly, to transform the precipitate orientation, sampleA was compressed along TD at a constant strain rate of 10-3s-1(hereafter named sample AC).Two strain values (8 and 13%) were carried out to obtain the maximum twinning proportion, which expected most grains were completely twinned.Thus, the grain orientation will change fromc-axis// ND toc-axis // TD.Finally, the compressed samples were annealed at 250°C for 30min (hereafter named sample ACA).This treatment could have a slight influence on the texture, but the main grain orientation remainsc-axis // TD.

    Fig.1.(a) A schematic diagram showing the sample produced in each step.The macro-texture of (b) the as-rolled AZ31 sheet and (c) the as-solutionized AZ31 sheet.(d) A schematic diagram showing the compression test process, and the grain orientation is illustrated in the test sample.In this condition, the c-axis of grain is 45° to the compressive axis.

    To avoid the influence of texture on mechanical properties, all the compression test samples with a dimension of 4.5×4.5×7.5mm (length×width×height) were machined along 45° to ND (in the ND-TD plane).The test samples produced in each step were denoted as sample A-45°, sample AC-45° and sample ACA-45°, respectively (see Fig.1a).Mechanical properties were evaluated by compression test with a strain rate of 1×10-3s-1at RT using a MTS-CMT5105.Note that the compressive axis is along the height direction,as demonstrated in Fig.1d.All compression tests were repeated 3 times for each type of samples to get the representative results.

    Morphology observations of the various samples were investigated by using an optical microscopy (OM, ZEISS Axio Lab A1) and a backscatter scanning electron microscopy(BSEM, TESCAN VEGA II LMU).The grain size (d) was determined by linear intercept method withd=1.74L (Lis the size of line intercept) [21].More than 200 grains were counted for each sample to obtain an average grain size.The evolution of the crystal orientation was analyzed by using an electron backscatter diffraction (EBSD, JEOL JSM-7800F)technique.Sample preparation for EBSD incorporated mechanical grinding, polishing with a diamond and colloidal silica slurry as well as electro-polishing at a voltage of 20V and an electric current of 0.1 A for 60s at a temperature of-40 °C in an AC2 electrolyte.All EBSD data were analyzed by using the channel 5 software.The samples for transmission electron microscopy (TEM) observations were prepared by ion-beam milling using a Gatan PIPS 695 and the precipitates were examined by using a JEOL2100 TEM operating at 200kV.

    3.Results

    3.1.Microstructure and texture evolution

    Fig.2 shows the BSEM images of various samples.The microstructure of as-rolled sheet consists of the equiaxed grain structure ofα-Mg and a few irregular granular secondary phases (see Fig.2a).The microstructure of the solution treated AZ31 sheet is free of precipitates (see Fig.2b).After aging at 180 °C for 14h, some irregular granular precipitates disperse on theα-Mg matrix (see Fig.2c).Note that such a heat treatment has no obvious influence on the grain size of Mg-Al alloys, which is consistent with the results reported in Refs.[22,23].

    Fig.3 shows the OM images of sample A, sample AC with a compression strain of 8 and 13% and sample ACA,respectively.The corresponding (0002) macro-pole figures of sample AC with a compression strain of 8 and 13% are also presented as insets.It is found that the average grain size of sample A is ~21.70μm.After compression along TD, the lamellar twins formed on the matrix.Fig.3b shows the microstructure of sample AC with a compression strain of 8%,the grain size is refined to 9.44μm, which is caused by the twin boundaries existing in the matrix.Similarly, the grain size of sample AC with a compression strain of 13% is refined to 11.85μm, as shown in Fig.3c.

    Fig.2.BSEM images of the samples: as-rolled sheet (a), solid solution treated at 410°C for 8h (b) and aged at 180°C for 14h (c).

    Fig.3.OM images of the samples: (a) sample A, sample AC with a compression strain of (b) 8% and (c) 13%, respectively, and (d) sample ACA.The twin boundaries are indicated by red arrows in Fig.3 (b) and (c).The corresponding (0002) macro-pole figures are presented in (b) and (c) as insets.The maximum texture intensity is given and marked by red word in the corresponding (0002) macro-pole figure (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.).

    It was reported [24]that the compression strain of 8% almost causes a completely twinning deformation (with ~94%volume fraction) for a rolled AZ80 Mg alloy, which significantly decreases the twin boundaries in most grains.However,in the present work, profuse twin boundaries still exist inside the grains with the compression strain of 8%.Besides, as demonstrated in Fig.3b, although most of (0002) poles tilt to TD, a few (0002) poles still remain in the center in the sample AC with a compression strain of 8%.With further strain to 13%, the twin boundaries are remarkably decreased and the (0002) poles completely tilt to TD, resulting in the formation of a twin-texture withc-axis // TD, as shown in Fig.3c.Combined the results of the macro-texture and twin morphology, it proves that the sample AC with a compression strain of 13% is suitable for subsequent research.

    In addition, annealing process was subjected to eliminate the hardening effect of dislocations or residual twin boundaries induced in the previous steps.Fig.3d shows the twin boundaries in sample ACA are significantly disappeared due to the annealing process, accompanied with coarser grains(average grain size of 22.44μm).

    Variations of the texture and crystal orientation of the various samples are compared in Fig.4.Here, Fig.4a shows the inverse pole figure (IPF) map of sample A, and the corresponding (0002) pole figure is presented as inset.This sample exhibits a typically strong basal texture, where the basal(0002) planes are mostly parallel to RD-TD plane.Fig.4b presents the pole figure and IPF map of the sample A-45°.Note that the examined plane is vertical to the compressive axis with a cross-section of 4.5×4.5 mm2.Compared with sample A, the pole axis of (0002) pole figure of sample A-45° has a titled angle of 45° from the center to ND.It reveals that thec-axes of grains in sample A-45° are 45° to the compressive axis (see Fig.1d).Moreover, it is well known that{102} extension twin can be easily induced by compression along TD or RD for as-rolled Mg alloys sheet and {102}extension twin can trigger a drastic change ofc-axis with a crystallographic lattice reorientation of 86.3° [15].Therefore,as shown in Fig.4d, the major texture component in sample AC isc-axis // TD.It means that thec-axes of grains are largely oriented parallel to TD.Besides, according to the pole figure shown in Fig.4e, it can be found that the basal poles in sample AC-45° largely tilt towards ND with an angle of 45°, and thus thec-axes of grains are 45° to the compressive axis.It indicates the grain orientation in sample AC-45° is similar to that in sample A-45°.

    Since the texture has a major effect on mechanical properties, the Schmid factor (SF) distribution for basal slip on< a >{0002}<11ˉ20>as well as the mean SF of sample A-45° and sample AC-45° were calculated from the EBSD raw data and are presented in Fig.4c and f.The mean SF for the basal slip of sample A-45° and sample AC-45° is 0.43 and 0.40, respectively.The distribution of the SF also indicates that the dominated SF of the grains in both sample A-45° and sample AC-45° concentrates between 0.4 and 0.5 for basal slip.Such grain orientations favor the activation of basal slip.

    It was reported that the precipitates in Mg-Al alloys,proved to be Mg17Al12, usually presented as the shape of plates and distributed in the basal plane of matrix [25].In general, these precipitates would display a plate-shaped morphology when viewed along [0001]αdirection while present a lath-shaped morphology when viewed along [100]αdirection.Analogously, the possible morphologies of prismaticplates in AZ31 observed along different directions could also be predicted.The schematic morphologies of prismatic plates are depicted in Fig.5a.Compared to basal plates,prismatic plates would exhibit the lath shape when viewed along[0001]αdirection, while show plate shape when viewed along [100]αdirection.

    Fig.4.EBSD inverse pole figure (IPF) maps of (a) sample A, (b) sample A-45°, (d) sample AC and (e) sample AC-45° The corresponding (0002) micro-pole figures are presented as insets.Distributions of the Schmid factor for basal slip in (c) sample A-45°and (f) sample AC-45°.

    To identify the orientation of the precipitates in sample AC,TEM bright-field images viewed along [100]αand [0001]αdirection were obtained.Fig.5b shows the bright-field TEM images viewed along [100]αdirection of sample AC.It indicates that the plate-shaped Mg17Al12precipitates are the major morphologies with their broad face approximately parallel to the prismatic plane of the matrix.Meanwhile, the morphologies of these precipitates viewed [0001]αdirection are lath-shaped (see Fig.5c).The corresponding electron diffraction patterns are presented as insets.As aforementioned, the lattice rotation is 86.3°, and the precipitates usually have a small rotation (<4°) at twin boundry [10].Therefore, the precipitates are not exactly parallel to the prismatic planes of the twinned region.In the present work, it can be concluded that the plate-shaped precipitates in sample AC have the HPs with their broad face approximately parallel to one of the{100}αprismatic plane of the twinnedα-Mg matrix with a small rotation angle (<4°) during the interaction with twin boundaries.

    3.2.Mechanical properties

    Fig.6 shows the compressive properties of sample A-45°,sample AC-45° and sample ACA-45°.The details of compressive properties are summarized in Table 1.It can be clearly seen that, compared to sample A-45°, the compressive yield strength (CYS) of sample AC-45° and sample ACA-45° increased by~40MPa and ~20MPa, respectively.Moreover, the compression ratio of sample AC-45° prolonged by~22% compared to sample A-45°, while the compression ratio of sample ACA-45° is at the same level as sample A-45°(~23%).

    Table 1The mechanical properties consist of the ultimate compressive strength(UCS), compressive yield strength (CYS) and compression ratio of sample A-45°, sample AC-45° and sample ACA-45°.

    3.3.Fracture properties

    As mentioned above, profuse twin boundaries still existed inside the grains in the sample AC-45°.These twin boundaries would exert an extra influence on fracture performance,such as void formation and growth [26].Therefore, the fracture performance of sample A-45° and sample ACA-45° was proposed to examine.

    Fig.7 shows the SEM images of fracture performance of sample A-45° and sample ACA-45°.Particularly, the morphologies of sample A-45° and sample ACA-45° at a compression strain of 20% are displayed in Fig.7a and d.It can be seen that the micro-cracks initiate at grain boundaries in both samples.Some slip traces were observed on the surface of sample ACA-45° (marked with white arrow in Fig.7d) whereas not observed in sample A-45°, as confirmed in Fig.7a and d.Besides, no deformation twins are found in either sample.

    Fig.7b, c and e, f show the morphologies of the fracture surface of sample A-45° and sample ACA-45°, respectively.The fracture surface of both samples composes by tear ridges,cleavage steps, cleavage surfaces and microvoids.However,compared to sample A-45°, more microvoids are observed in sample ACA-45°.

    Fig.5.(a) A schematic diagram showing the morphologic characteristics of Mg17Al12 precipitates viewed along [100]α and [0001]α.Bright-field TEM images of sample AC (b) viewed along [100]α and (c) viewed along [0001]α.

    Fig.6.Compressive stress-strain curve of sample A-45°, sample AC-45° and sample ACA-45°.

    4.Discussion

    4.1.The mechanism for the improved strength

    As mentioned above, the precipitate orientation was altered from basal plane to prismatic plane by the aging prior to twinning deformation process, accompanied with the CYS improvement.In the present work, the influence of texture on the CYS can be avoided due to the similar grain orientation and close SF.Therefore, the possible CYS strengthening mechanism can be mainly divided into 3 parts, namely,(1) solid solution strengthening contribution,σss, (2) precipitation strengthening contribution,σpptand (3) grain boundary strengthening contribution,σgb.The contributions of these strengthening mechanisms to CYS are further quantitatively discuss as follows:

    Fig.7.SEM images of (a) sample A-45°and (d) sample ACA-45° with a compression strain of 20%.The compression direction is indicated by blue arrow in (a).The fracture surface of (b, c) sample A-45° and (e, f) sample ACA-45° The regions surronded by red lines in (b) and (e) are magnified and shown in(c) and (f), respectively.The microcracks, cleavage steps, cleavage surfaces, microvoids and tear ridge are marked by yellow arrows, while the slip traces are marked by whilte arrows (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.).

    4.1.1.Solid solution strengthening

    Solid solution strengthening by Al [27], Zn [28], Sn [29],Rare-Earths (RE) [30–32]and other solutes is an important mechanism to improve the room and elevated temperature(RT and ET) strength of Mg alloys.To estimate the solid solution strengthening effect on basal slip, Yasi et al.[33]developed a model ofσssas a function of the solute concentration and a potency factor for each solute.Taken with Taylor factor,σssis given as follows:

    wherem=4.5 is the Taylor factor,Mxis the potency factor of solutexandcxis the solute concentration in the matrix.However,the Al solutes would precipitate in the form of Mg17Al12particles during the aging process.Thus, in the present work,theσssof Al and Zn can be neglected due to their few solid solubility present in the matrix.

    4.1.2.Precipitation strengthening

    For precipitation strengthening, the structure, morphology and orientation of precipitates play important roles to determine the YS of Mg alloys [9,34].Specifically, the different precipitate orientation can bring significantly different strengthening effect on Mg matrix.Many researches[11,12,35,36]based on Orowan strengthening theory were carried out to predict the strengthening effect of different precipitate orientation.At the same time, these predictions were also verified by the experimental results [8,12,37].Many experimental researches of Mg alloys deformed at RT and theoretical considerations based on single-crystal data suggest that basal slip will dominate the RT deformation process due to the fact that the critical resolved shear stress (CRSS) of basal slip is much lower than non-basal slip systems [4].Specifically, as shown in Fig.4, the average SF for basal slip of the compression samples are both over 0.4, it comfirms that the deformation of AZ31 Mg alloys at yielding occurs by basal slip [4].This suggests that the CYS improvement mainly imputes to the hardening arising from basal slip.In order to quantify the hardening effect of basal plates and prismatic plates on basal slip, the increment of CRSS (Δτ) originated by the need for dislocations to bypass the two different precipitates was estimated as follows.TheΔτis given as follows[11]:

    whereGis the shear modulus of the Mg matrix,bis the magnitude of the Burgers vector for the gliding dislocations,λis the effective planar inter-particle spacing on the slip plane,vis the Poisson’s ratio of the Mg matrix,dpis the mean planar diameter of the particles on the slip plane, andr0is the core radius of the dislocations, which is approximated to be the magnitude ofb[22].

    Based on Eq.(2), the appropriate version of the Orawan equations predicting the strengthening effect of basal plates and prismatic plates on basal slip were developed by Nie[11].For basal plates, Eq.(2) is rewritten in the form

    For prismatic plates, Eq.(2) is rewritten in the form

    whereD, Tandfis the uniform diameter, thickness and volume fraction of the plate-shaped precipitates, respectively.

    According to these equations,theΔτis usually determined by the size (i.e., the uniform diameterDand thicknessTof plates) and the volume fraction (f) of precipitates.In the present study, according to the TEM images, theD, T, andfof plates were confirmed as 32.5nm, 22.4nm and 0.04, respectively.For Mg alloys, theGandvis 16.5GPa and 0.35,respectively.As a consequence, the calculatedΔτof basal slip produced by basal plates and prismatic plates is 40.87 and 59.67MPa, respectively.Clearly, theΔτof basal slip produced by prismatic plates is ~20MPa higher than that produced by basal plates.This is the reason why the sample ACA-45° containing prismatic plates shows the higher CYS than sample A-45° containing basal plates.It is noteworthy that the CYS of sample AC-45° containing prismatic plates is higher than that of sample ACA-45° also containing prismatic plates.This additional strengthening effect should be attributed to the grain refinement and introduced dislocation caused by twinning deformation [38].This part will be discussed in detail in the next section.After annealing, the twin boundaries of sample AC-45° were significantly disappeared,as shown in Fig.3d.Therefore, the hardening effect of sample ACA-45° should be largely attributed to theΔτof basal slip produced by prismatic plates.

    Although, in the present study, the influence of texture on the CYS can be avoided due to the similar grain orientation and close SF, it has great affect on the YS of Mg alloys indeed.According to the Schmid law [39], theσpptis given as follows:

    wheremsis the SF for basal slip.Clearly, the lower SF for basal slip is, the higherσpptis.According to the discussion above, different precipitate orientations make the distinct strengthening effect on basal slip.Hence, for basal plates precipitated alloys,e.g.Mg-Al alloys etc., substituting Eq.(3) into Eq.(5), the appropriateσpptcan be rewritten as follows:

    Similarly, for prismatic plates precipitated alloys,e.g.Mg-RE alloys etc., the appropriateσpptcan be rewritten as follows:

    As discussed before, it is speculated that theσpptproduced by prismatic plates may always be larger than theσpptproduced by basal plates.To further clarify the different hardening effects on basal slip produced by basal plates and prismatic plates,a schematic diagram is plotted to show the effect of these plates on dislocation gilding on (0002) basal plane of the Mg matrix.As described in Fig.8, the slip plane represents the(0002)basal plane of Mg matrix;the red arrows represent the dislocation motion direction; the black dislocations represent the dislocations before slipping; the red dislocations represent the dislocations after slipping and the basal plates and prismatic plates are inserted in Mg matrix,respectively.It implies that basal plates exert a weak blocking effect on dislocation motion, thereby causing a relatively lowΔτ.Thus,the basal slip is easy to occur and the Mg alloys containing such precipitates usually exhibit a lower YS, such as Mg-Al alloy systems [4].However, according to the description of prismatic plates, it usually exerts a strong blocking effect on dislocation motion, thereby causing a highΔτ[11].Thus,the basal slip is hard to occur and the Mg alloys containing prismatic plates usually exhibit a higher YS, such as novel Mg-Gd-Y alloy systems [40,41]and Mg-Nd alloy systems[9,42].

    4.1.3.Grain boundary strengthening

    It is well-known that grain boundary strengthening is a feasible way to enhance the YS of Mg alloys due to the considerable strengthening coefficient [43].According to Hall-Petch relationship [43], theσgbis given as follows:

    wherekis the Hall-Petch coefficient (158MPa)[43]anddis the average grain size.Clearly, the smaller grain size causes the great contribution to YS[43].Therefore, compared to sample A-45°and sample ACA-45°,the CYS improvement of sample AC-45° originated from grain refinement can be given as follows:

    The average grain sizes of sample A-45°, sample AC-45°and sample ACA-45° were illustrated in Fig.3.Accordingly,the grain refinement contribution to YS of sample AC-45° is~12MPa higher than sample A-45° and sample ACA-45°.To reveal the potential mechanism, a schematic diagram describing the microstructure evolution during compression along TD is depicted in Fig.9.For compression along TD, the precipitate orientation was altered, accompanied with the formation of twin boundaries (red line) and dislocations (green symbol),thereby causing the grain refinement strengthening and dislocation strengthening [44], which make extra contributions to the YS.Thus, the sample AC-45° exhibits a greater strengthening effect.

    4.1.4.A yield strength model

    Combined the above strengthening mechanisms, a total YS(σy) model with the dominated basal slip in Mg alloy was developed as:

    Fig.8.A schematic diagram showing the effect of the orientation of precipitate plates on dislocation gilding on basal plane of the magnesium matrix.The brown plates represent the plate-shaped precipitates.The F denotes the force applied on matrix and λ is the angle between the force axis and the slip direction.The red arrow represents the slip direction.The black “dislocation symbol” represents the state of the dislocation before slipping, while the red “dislocation symbol” represents the state of the dislocation after slipping (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.).

    Fig.9.A schematic diagram describing the microstructure evolution during compression along TD.The brown plates represent the plate-shaped precipitates.The red lines denote the twin boundaries and the green “dislocation symbol” represents the dislocation caused by twinning deformation.The Mg crystal shows the grain orientation evolution during the compression process (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.).

    Note that theσ0can be determined via the SF law[39]giving as follows:

    whereτ0is the initial CRSS of basal slip(5MPa)[45],Δτ/msis considered asσppt, and thusσ0equals toτ0/ms.

    The YS model is proposed to analyze the effect of multiple factors, including solute, precipitation behavior, grain boundary and texture, on the YS of Mg alloys.As for precipitation behavior, it refers to precipitate orientation, size and volume fraction.Hence, for basal plates precipitated alloys with dominated basal slip,e.g.Mg-Al alloysetc., substituting Eqs.(1),(6) and (8) into Eq.(10), the appropriate YS model can be rewritten as follows:

    Similarly,for prismatic plates precipitated alloys with dominated basal slip,e.g.Mg-RE alloysetc., substituting Eqs.(1),(7) and (8) into Eq.(10), the appropriate YS model can be rewritten as follows

    For a micro-alloyed or a peak aged Mg alloy, such as the aged AZ31 investigated here, the solid solution strengthening contribution can be neglected due to the rare solute concentration in the matrix.Conbined with the present experimental results, when we discuss the variation of the YS as a function ofDand SF, thefanddare assumed to be 0.04 and 22μm,respectively.Besides, the ratio ofTtoD, namely aspect ratio,is supposed to be 0.7.When we discuss the variation of the YS as a function ofdand SF, the precipitate size is assumed thatD=32nm andT=22nm.The SF is assumed to be 0.4 when we discuss the variation of the YS as a function ofDandd.Accordingly, the trend of the YS as the function ofD,dand SF is depicted in Fig.10.

    As shown in Fig.10, the theoretical YS increases with the decreasingD, dand SF, whether the alloy containing pris-matic plates or basal plates.When the SF for basal slip is close to 0, the theoretical YS will extremely increase to more than 1000MPa.This is due to other deformation systems will dominate the plastic deformation process [46].For instance,the extruded Mg alloy bar owns a strong<100>fiber texture in which the basal planes largely parallel to the extrusion direction (ED) [37].When tested in tension or compression along ED, the SF for basal slip is 0.However, it demonstrated that the activation of prismatic slip controls the YS in tension and the activation of {102} twinning controls the YS in compression [37].To compare the YS of the alloy containing prismatic plates and basal plates, the YS ratio of the alloy containing prismatic plates to the alloy containing basal plates are depicted in Fig.10g–i.It is found that the YS ratio is always larger than 1 with the same varites of SF,dandD.It can further imply that the YS of the alloy containing prismatic plates is greater than that containing basal plates continuously, confirming that the Mg alloys formed the prismatic plates always exhibit the appreciable YS[9,42].

    Fig.10.The variation of the YS as a function of D, d and SF of (a)–(c) the alloy containing basal plates and (d)–(e) the alloy containing prismatic plates.(g)–(i) The YS ratio of the alloy containing prismatic plates to the alloy containing basal plates as a function of D, d and SF.Specifically, (a), (d) and (g) the YS/YS ratio as a function of D and SF; (b), (e) and (h) the YS/YS ratio as a function of d and SF; (c), (f), (i) the YS/YS ratio as a function of d and D.

    4.2.Effect of precipitate type on ductility

    It was reported that although prismatic plates have a strong blocking effect on basal slip, the blocking effect on prismatic slip is relatively weak [11].Wang et al.[36]developed the appropriate Orowan equations to quantify the strengthening effect on prismatic slip in HCP metals.For basal plates, the appropriate version of the Orowan equation is given as

    For prismatic plates,the appropriate version of the Orowan equation is given as

    Based on the counted results ofD, Tandffrom the TEM images, theΔτof basal and prismatic slip produced by basal plates and prismatic plates are calculated and listed in Table 2.Furthermore, the total CRSS (τ) of basal and prismatic slip can be obtained via the equation [22]:τ=τ0+Δτ, whereτ0is the initial CRSS of basal slip and prismatic slip.Note that theτ0of basal and prismatic slip adopted in the present study equals to 5 and 85MPa [45], respectively.

    Table 2The calculated results of Δτ and τ of basal and prismatic slip according to the developed Orowan equations in HCP structure.

    Fig.11.The histograms represent the CRSS for basal and prismatic slip which hardened by basal and prismatic plates, respectively.The black dashed line shows a descreasing trend of the CRSS ratio of prismatic slip to basal slip with the precipitate orientation changed from basal plates to prismatic plates.

    It is found that the CRSS increment for prismatic slip caused by basal plates and prismatic plates is 42.68 and 41.91MPa, respectively.Subsequently, the total CRSS for prismatic slip strengthened by prismatic plates (126.91MPa)is lower than that strengthened by basal plates (127.68MPa).

    The calculated results of CRSS are also plotted in Fig.11.Obviously, the CRSS ratio of prismatic slip to basal slip(τprism/τbasal) in the sample containing prismatic plates (1.96)decreases compared to the sample containing basal plates(2.78).It was reported [27,47]that the descreased CRSS ratio of non-basal slip to basal slip decreased, provided more opportunity to activate non-basal slips, which can enhance the ductility of Mg alloys at RT.Specifically, there is a high activity of prismatic slip at RT with aτprism/τbasalbetween 1.1 and 2.5 [48].Eventually, prismatic slip might be activated in sample ACA-45°, resulting in maintaining excellent ductility.It could account for the phenomenon that the ductility and strength of sample AC-45° improved simultaneously, and the strength of sample ACA-45° enhanced but the ductility not decreased.

    Generally, the presence of the cleavage surface is usually the characteristic of brittle fracture.On the contrary, the presence of micro void is usually the characteristic of ductile fracture [49].These two fracture characteristics were employed simultaneously in sample A-45°and sample ACA-45°.However, compared to sample A-45°, more microvoids were formed on the fracture surface of sample ACA-45°, which indicates that sample ACA-45° has better plastic deformation capacity than sample A-45°.This phenomenon might also confirm the activation of non-basal slip systems.More slip systems could provide more mobile dislocations and less deformation resistance in the latter period of the deformation,resulting in excellent ductility.In terms of sample AC-45°, it shows better ductility than sample ACA-45°.This additional improvement in ductility can be attributed to the grain refinement caused by twin boundary [38].

    5.Conclusion

    In the present study, the effects of altered precipitate orientation on mechanical properties in AZ31 Mg alloy were investigated, and the corresponding strengthening mechanisms were discussed in detial.The key conclusions are summarized as follows:

    1.The grain orientation was oriented fromc-axis//ND tocaxis//TD due to {101ˉ2} extension twinning process.Thus,the precipitate orientation in AZ31 Mg alloy was altered from basal plates to prismatic plates.

    2.The CYS of sample AC-45°and sample ACA-45°with the prismatic plates improved by ~40MPa and ~20MPa than sample A-45°, respectively.Meanwhile, the compression ratio of sample AC-45° raised by 22% than sample A-45.

    3.A YS model was developed to analyze the strengthening mechanisms.It was revealed that the optimized strength was attritubed to the altered precipitate orientation, refined grain size and dislocations produced by aging prior to twinning deformation process.Particularly, prismatic plates always exert a stronger blocking effect on basal slip than basal plates under the same variates ofD, dand SF.

    4.Theτprism/τbasal of the sample AC-45° strengthened by prismatic plate was lower than that of the sample A-45°strengthened by basal plate,which gives the chance to activate prismatic slip.These slip systems could provide more mobile dislocations and less deformation resistance in thelatter period of the deformation, resulting in excellent ductility.

    Declaration of Competing Interests

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgments

    The authors gratefully acknowledge the Fundamental Research Funds for the Project of Science & Technology Department of Sichuan Province (2018HH0026), National Natural Science Foundation of China (51701132, U1764253) and China Scholarship Council (201907005018).We also thank the Analytical and Testing Center of Southwest Jiaotong University for assistance with SEM and EBSD experiments.

    国产 一区精品| 亚洲,一卡二卡三卡| 99九九线精品视频在线观看视频| 国产高清不卡午夜福利| 久久精品国产亚洲av天美| 99热这里只有是精品在线观看| 国产亚洲精品久久久com| av又黄又爽大尺度在线免费看| 日韩大片免费观看网站| 免费av不卡在线播放| 国产精品不卡视频一区二区| 国产精品国产av在线观看| 人妻制服诱惑在线中文字幕| 午夜视频国产福利| 极品少妇高潮喷水抽搐| 18禁动态无遮挡网站| 三级国产精品片| 欧美xxxx性猛交bbbb| 精品人妻熟女av久视频| 最近手机中文字幕大全| 一边亲一边摸免费视频| 少妇丰满av| 男女免费视频国产| 久久精品熟女亚洲av麻豆精品| 国产无遮挡羞羞视频在线观看| 久久精品国产鲁丝片午夜精品| 国产大屁股一区二区在线视频| 又爽又黄a免费视频| 黄色欧美视频在线观看| 在线观看人妻少妇| 国产精品久久久久久久久免| 熟女人妻精品中文字幕| 秋霞在线观看毛片| 麻豆成人午夜福利视频| 在线免费观看不下载黄p国产| 精品久久国产蜜桃| 亚洲国产精品专区欧美| 99热国产这里只有精品6| 精品一区二区免费观看| 国产精品一及| 亚洲一级一片aⅴ在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美xxxx黑人xx丫x性爽| 亚洲欧洲日产国产| 九色成人免费人妻av| 中文精品一卡2卡3卡4更新| 日韩av免费高清视频| 国产亚洲5aaaaa淫片| 国内少妇人妻偷人精品xxx网站| 永久免费av网站大全| 边亲边吃奶的免费视频| 亚洲av国产av综合av卡| 在线观看免费日韩欧美大片 | 下体分泌物呈黄色| 在线观看美女被高潮喷水网站| 亚洲国产精品专区欧美| 亚洲欧美一区二区三区黑人 | 自拍偷自拍亚洲精品老妇| 秋霞伦理黄片| 成人美女网站在线观看视频| 伊人久久国产一区二区| 免费黄网站久久成人精品| 少妇人妻精品综合一区二区| 女性被躁到高潮视频| 亚洲国产最新在线播放| 毛片一级片免费看久久久久| 深爱激情五月婷婷| 少妇裸体淫交视频免费看高清| 伦理电影免费视频| 免费高清在线观看视频在线观看| 99久久中文字幕三级久久日本| 亚洲国产精品专区欧美| 日本黄色日本黄色录像| 国产一级毛片在线| 国产黄色视频一区二区在线观看| 久久精品国产亚洲网站| 春色校园在线视频观看| 色婷婷久久久亚洲欧美| 国产片特级美女逼逼视频| 毛片女人毛片| 男男h啪啪无遮挡| av卡一久久| 色婷婷久久久亚洲欧美| 国产 精品1| 丝袜脚勾引网站| 国产片特级美女逼逼视频| 久久久色成人| 麻豆乱淫一区二区| 国产精品久久久久久精品古装| 国产人妻一区二区三区在| 亚洲色图综合在线观看| 欧美xxxx黑人xx丫x性爽| 久久97久久精品| 好男人视频免费观看在线| 色婷婷久久久亚洲欧美| 我要看黄色一级片免费的| 日日摸夜夜添夜夜爱| av国产精品久久久久影院| 精品少妇久久久久久888优播| 十分钟在线观看高清视频www | 亚洲人成网站在线播| 国产美女午夜福利| 王馨瑶露胸无遮挡在线观看| 黑人猛操日本美女一级片| 久久久久久伊人网av| 少妇丰满av| 日韩欧美一区视频在线观看 | 亚洲av成人精品一二三区| 永久免费av网站大全| 中文欧美无线码| 亚洲国产高清在线一区二区三| av福利片在线观看| 高清黄色对白视频在线免费看 | 国产欧美亚洲国产| 九草在线视频观看| 国产精品.久久久| av天堂中文字幕网| 久久久久久久久久人人人人人人| 亚洲成人av在线免费| 波野结衣二区三区在线| 亚洲av日韩在线播放| 国产黄频视频在线观看| 深夜a级毛片| 大片电影免费在线观看免费| 亚洲精品久久午夜乱码| 欧美日韩国产mv在线观看视频 | 老熟女久久久| 亚洲丝袜综合中文字幕| 男女无遮挡免费网站观看| 国产又色又爽无遮挡免| 久久亚洲国产成人精品v| 国产在线视频一区二区| www.av在线官网国产| 国产精品.久久久| 久久久久久九九精品二区国产| 这个男人来自地球电影免费观看 | 国产精品蜜桃在线观看| 一级av片app| 亚洲人与动物交配视频| 婷婷色综合www| 精品久久久精品久久久| 亚洲综合精品二区| 久久久久国产精品人妻一区二区| 日本欧美国产在线视频| 中文乱码字字幕精品一区二区三区| 国产精品99久久99久久久不卡 | 国模一区二区三区四区视频| 六月丁香七月| 国产综合精华液| 国产黄频视频在线观看| 亚洲人成网站在线观看播放| 免费黄频网站在线观看国产| 日产精品乱码卡一卡2卡三| 插逼视频在线观看| 国产精品久久久久久精品电影小说 | 一级片'在线观看视频| 91久久精品电影网| 国产在线视频一区二区| 身体一侧抽搐| 国产精品一区二区性色av| 免费观看av网站的网址| 国产精品.久久久| 尾随美女入室| 女性生殖器流出的白浆| 国产一区二区三区综合在线观看 | 国产高清国产精品国产三级 | 欧美3d第一页| av在线播放精品| 王馨瑶露胸无遮挡在线观看| 婷婷色综合www| 久热久热在线精品观看| 婷婷色av中文字幕| 在线 av 中文字幕| 校园人妻丝袜中文字幕| 久久久久人妻精品一区果冻| 国产黄频视频在线观看| 高清午夜精品一区二区三区| 免费人妻精品一区二区三区视频| 亚洲欧美一区二区三区黑人 | 99热这里只有精品一区| 久久久久网色| 熟妇人妻不卡中文字幕| 亚洲图色成人| 国产精品爽爽va在线观看网站| 久久久午夜欧美精品| 国产午夜精品久久久久久一区二区三区| 高清欧美精品videossex| 天美传媒精品一区二区| 国产精品久久久久成人av| 美女内射精品一级片tv| 亚洲欧美清纯卡通| 精品久久久久久久末码| 欧美一区二区亚洲| 国产精品秋霞免费鲁丝片| 日本-黄色视频高清免费观看| 在线精品无人区一区二区三 | 中文字幕免费在线视频6| 国产中年淑女户外野战色| 亚洲欧美日韩另类电影网站 | 成年免费大片在线观看| 春色校园在线视频观看| 好男人视频免费观看在线| 国产乱来视频区| 女的被弄到高潮叫床怎么办| 亚洲av成人精品一区久久| 男人添女人高潮全过程视频| 99re6热这里在线精品视频| 久久久久久久久久久免费av| 国产精品.久久久| 国产亚洲5aaaaa淫片| 少妇人妻一区二区三区视频| 国产av精品麻豆| 国产有黄有色有爽视频| 国产一区二区三区av在线| 免费久久久久久久精品成人欧美视频 | 岛国毛片在线播放| 日日啪夜夜撸| 最新中文字幕久久久久| 日韩成人av中文字幕在线观看| 国产成人a∨麻豆精品| av国产久精品久网站免费入址| 日韩欧美精品免费久久| 亚洲av综合色区一区| 久久99热这里只频精品6学生| 欧美另类一区| 国产成人精品婷婷| 高清日韩中文字幕在线| 国产成人免费无遮挡视频| 国产免费一区二区三区四区乱码| 国产免费一区二区三区四区乱码| 嫩草影院新地址| 男男h啪啪无遮挡| 亚洲av综合色区一区| 99久久人妻综合| 欧美成人精品欧美一级黄| 久久精品国产鲁丝片午夜精品| 色综合色国产| 亚洲最大成人中文| 看十八女毛片水多多多| 在线看a的网站| 国产高清不卡午夜福利| 国产一区二区三区综合在线观看 | 少妇人妻久久综合中文| 精品国产乱码久久久久久小说| 在线观看一区二区三区| 国产精品爽爽va在线观看网站| 色哟哟·www| 久久影院123| 久久99热6这里只有精品| 乱码一卡2卡4卡精品| 日韩欧美一区视频在线观看 | 丰满乱子伦码专区| 91久久精品国产一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 日韩 亚洲 欧美在线| 欧美 日韩 精品 国产| 中文乱码字字幕精品一区二区三区| 高清欧美精品videossex| 热99国产精品久久久久久7| 99久久精品国产国产毛片| 欧美成人精品欧美一级黄| 欧美极品一区二区三区四区| 99久久人妻综合| 精品一区二区三卡| 男人添女人高潮全过程视频| 亚洲电影在线观看av| 日韩av不卡免费在线播放| 精华霜和精华液先用哪个| 在线观看国产h片| 国产亚洲午夜精品一区二区久久| 丰满迷人的少妇在线观看| 嫩草影院入口| 国产高清国产精品国产三级 | 国产爱豆传媒在线观看| 国产一区亚洲一区在线观看| 校园人妻丝袜中文字幕| 久久久久久久久大av| 亚洲精品乱久久久久久| 色婷婷av一区二区三区视频| 亚洲成人手机| 中文字幕精品免费在线观看视频 | 高清毛片免费看| 亚洲不卡免费看| 直男gayav资源| 97超碰精品成人国产| 男人舔奶头视频| 日韩中字成人| 男女下面进入的视频免费午夜| 亚洲自偷自拍三级| 国产精品三级大全| 少妇丰满av| 久久青草综合色| 黑人高潮一二区| 国产黄频视频在线观看| 综合色丁香网| 国产精品熟女久久久久浪| 啦啦啦啦在线视频资源| 色视频www国产| 欧美精品亚洲一区二区| 中文在线观看免费www的网站| 97超碰精品成人国产| 欧美区成人在线视频| 久久99精品国语久久久| 亚洲精品国产色婷婷电影| 久久6这里有精品| 在现免费观看毛片| 久久久久网色| 大香蕉97超碰在线| 99久久人妻综合| 人体艺术视频欧美日本| 国产 一区精品| 国产成人精品婷婷| 高清不卡的av网站| 成人18禁高潮啪啪吃奶动态图 | 综合色丁香网| 久久久久网色| 久久久久久人妻| 国产午夜精品一二区理论片| av在线观看视频网站免费| 亚洲国产精品专区欧美| 日韩欧美精品免费久久| 久热久热在线精品观看| 中国美白少妇内射xxxbb| 国产精品国产三级专区第一集| av播播在线观看一区| 97在线视频观看| 日本av手机在线免费观看| 久久精品国产鲁丝片午夜精品| 久久精品国产亚洲av涩爱| 1000部很黄的大片| 日本wwww免费看| 能在线免费看毛片的网站| 午夜激情福利司机影院| 韩国高清视频一区二区三区| 久久这里有精品视频免费| 亚洲精品自拍成人| 免费不卡的大黄色大毛片视频在线观看| 久久人人爽人人爽人人片va| 精品久久久久久久久亚洲| 一区二区三区乱码不卡18| 一级爰片在线观看| 国产熟女欧美一区二区| 欧美少妇被猛烈插入视频| 麻豆乱淫一区二区| 亚洲第一av免费看| 欧美老熟妇乱子伦牲交| 少妇被粗大猛烈的视频| 有码 亚洲区| 久久97久久精品| 欧美丝袜亚洲另类| 国产亚洲最大av| 精品久久久久久久末码| 三级国产精品片| 国产 一区 欧美 日韩| 国产成人精品福利久久| 日韩一本色道免费dvd| 两个人的视频大全免费| 色视频在线一区二区三区| 日韩 亚洲 欧美在线| 国产精品99久久久久久久久| av在线观看视频网站免费| 九九久久精品国产亚洲av麻豆| 五月天丁香电影| 国产久久久一区二区三区| 日本欧美国产在线视频| 校园人妻丝袜中文字幕| 少妇熟女欧美另类| 纯流量卡能插随身wifi吗| 天堂中文最新版在线下载| 99久久中文字幕三级久久日本| 国产乱来视频区| 久热久热在线精品观看| 亚洲在久久综合| 汤姆久久久久久久影院中文字幕| 亚洲自偷自拍三级| 欧美xxxx性猛交bbbb| 精品亚洲成国产av| 亚洲国产日韩一区二区| 亚洲欧美成人综合另类久久久| 91久久精品国产一区二区三区| 免费观看无遮挡的男女| 免费看日本二区| 亚洲国产高清在线一区二区三| 99热国产这里只有精品6| 高清av免费在线| 国内少妇人妻偷人精品xxx网站| 一区在线观看完整版| 一级毛片我不卡| 涩涩av久久男人的天堂| 男的添女的下面高潮视频| 99久久中文字幕三级久久日本| 中国国产av一级| 中文乱码字字幕精品一区二区三区| 精品一品国产午夜福利视频| 一级毛片我不卡| 欧美高清性xxxxhd video| 欧美激情极品国产一区二区三区 | 国产精品久久久久成人av| 国产一级毛片在线| h日本视频在线播放| 国产精品偷伦视频观看了| 亚洲精品国产成人久久av| 六月丁香七月| 特大巨黑吊av在线直播| 免费观看在线日韩| 黄片wwwwww| 久久青草综合色| 中文在线观看免费www的网站| 国产色婷婷99| 人人妻人人看人人澡| 亚洲欧美一区二区三区黑人 | 国产精品人妻久久久影院| 国产精品久久久久成人av| 尤物成人国产欧美一区二区三区| 欧美变态另类bdsm刘玥| 精品酒店卫生间| 午夜激情久久久久久久| 成人无遮挡网站| 又爽又黄a免费视频| 男女边摸边吃奶| 午夜激情久久久久久久| 久久久久精品性色| 亚洲国产日韩一区二区| 亚洲精品,欧美精品| 丰满迷人的少妇在线观看| 51国产日韩欧美| 十八禁网站网址无遮挡 | 国产精品99久久久久久久久| 亚洲精品乱码久久久久久按摩| 少妇被粗大猛烈的视频| 国产视频首页在线观看| 中文在线观看免费www的网站| 欧美三级亚洲精品| 亚洲精品色激情综合| 色视频www国产| 日韩国内少妇激情av| 久久久久国产网址| 国产亚洲欧美精品永久| 免费观看性生交大片5| 欧美 日韩 精品 国产| 国产精品久久久久久久电影| 一区二区三区精品91| 亚洲国产av新网站| 高清不卡的av网站| 天堂俺去俺来也www色官网| 插逼视频在线观看| av国产免费在线观看| 免费不卡的大黄色大毛片视频在线观看| av国产精品久久久久影院| 少妇人妻精品综合一区二区| 成年女人在线观看亚洲视频| 黄色视频在线播放观看不卡| 精品亚洲成a人片在线观看 | 欧美老熟妇乱子伦牲交| 午夜免费男女啪啪视频观看| 人人妻人人澡人人爽人人夜夜| 97超视频在线观看视频| 亚洲国产成人一精品久久久| 高清午夜精品一区二区三区| 91午夜精品亚洲一区二区三区| 王馨瑶露胸无遮挡在线观看| av线在线观看网站| 精品酒店卫生间| 亚洲激情五月婷婷啪啪| 五月天丁香电影| 欧美97在线视频| 性高湖久久久久久久久免费观看| 精品一区在线观看国产| 成年人午夜在线观看视频| 成年免费大片在线观看| 国产日韩欧美亚洲二区| 哪个播放器可以免费观看大片| 国产一区亚洲一区在线观看| 免费高清在线观看视频在线观看| 日韩av免费高清视频| 国产亚洲一区二区精品| 欧美激情国产日韩精品一区| 成人黄色视频免费在线看| 少妇人妻精品综合一区二区| 亚洲经典国产精华液单| 亚洲国产欧美人成| 五月天丁香电影| 欧美97在线视频| tube8黄色片| 男人舔奶头视频| 亚洲内射少妇av| 中文字幕亚洲精品专区| 2022亚洲国产成人精品| 日韩制服骚丝袜av| 色视频在线一区二区三区| 国语对白做爰xxxⅹ性视频网站| 少妇熟女欧美另类| 国产精品国产av在线观看| 亚洲第一av免费看| 国内揄拍国产精品人妻在线| 国产乱人视频| 色吧在线观看| 一本久久精品| 人妻系列 视频| 中文天堂在线官网| 观看av在线不卡| 久久精品久久精品一区二区三区| 久久国产精品男人的天堂亚洲 | 亚州av有码| 久久青草综合色| 国产深夜福利视频在线观看| 免费看不卡的av| 一区在线观看完整版| 免费大片18禁| 亚洲丝袜综合中文字幕| 日本与韩国留学比较| 欧美性感艳星| 99久久精品一区二区三区| 国产精品一区二区性色av| 久久影院123| 少妇丰满av| 国产乱来视频区| 国产 一区 欧美 日韩| 国产精品国产三级专区第一集| 国产片特级美女逼逼视频| 欧美极品一区二区三区四区| 大又大粗又爽又黄少妇毛片口| 只有这里有精品99| 国产精品嫩草影院av在线观看| 欧美 日韩 精品 国产| 永久网站在线| 天天躁日日操中文字幕| 一级爰片在线观看| 国产精品一区二区性色av| 午夜福利影视在线免费观看| 人人妻人人看人人澡| 2018国产大陆天天弄谢| 偷拍熟女少妇极品色| 极品少妇高潮喷水抽搐| 日韩伦理黄色片| 春色校园在线视频观看| 精品视频人人做人人爽| 欧美精品国产亚洲| 国产精品蜜桃在线观看| 成人毛片60女人毛片免费| 91在线精品国自产拍蜜月| 免费人成在线观看视频色| 国产精品三级大全| 久久97久久精品| 精品视频人人做人人爽| 26uuu在线亚洲综合色| 亚洲av综合色区一区| 黄色配什么色好看| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲一区二区三区欧美精品| 丰满人妻一区二区三区视频av| 一级黄片播放器| 午夜免费鲁丝| 黄色怎么调成土黄色| 成人影院久久| 丰满少妇做爰视频| 一区二区三区乱码不卡18| 五月天丁香电影| 色婷婷久久久亚洲欧美| 久久97久久精品| 观看美女的网站| 下体分泌物呈黄色| 成人免费观看视频高清| 日韩av在线免费看完整版不卡| 国产精品一区www在线观看| 免费看光身美女| 又粗又硬又长又爽又黄的视频| 最近最新中文字幕大全电影3| 91aial.com中文字幕在线观看| 久久久久性生活片| av在线老鸭窝| 日本与韩国留学比较| 亚洲,一卡二卡三卡| 亚洲最大成人中文| 亚洲国产精品专区欧美| 女性被躁到高潮视频| 亚洲成人中文字幕在线播放| 日韩亚洲欧美综合| 黄色配什么色好看| 欧美成人a在线观看| 一级爰片在线观看| 亚洲av.av天堂| 如何舔出高潮| 亚洲人成网站在线观看播放| 精品国产三级普通话版| 欧美成人一区二区免费高清观看| 一个人看的www免费观看视频| av一本久久久久| 久久精品国产自在天天线| 下体分泌物呈黄色| 成人综合一区亚洲| 五月玫瑰六月丁香| 最近手机中文字幕大全| 九草在线视频观看| 亚洲精品国产成人久久av| 欧美老熟妇乱子伦牲交| 色婷婷av一区二区三区视频| 久久久久精品性色| 久久99精品国语久久久| 91精品伊人久久大香线蕉| 人妻 亚洲 视频| 黑人高潮一二区| 亚洲人成网站高清观看| 精品一区二区免费观看| 国产成人a区在线观看| 亚洲欧美精品专区久久| 又大又黄又爽视频免费| 在线观看国产h片| 免费看光身美女| 久久久久视频综合| 美女视频免费永久观看网站| 久久精品国产自在天天线| 亚洲精品第二区| 国产一级毛片在线| 亚洲欧美一区二区三区黑人 |